[en] Cotton leaf curl disease (CLCuD), caused by cotton leaf curl viruses (CLCuVs), is among the most devastating diseases in cotton. While the widely cultivated cotton species Gossypium hirsutum is generally susceptible, the diploid species G. arboreum is a natural source for resistance against CLCuD. However, the influence of CLCuD on the G. arboreum transcriptome and the interaction of CLCuD with G. arboreum remains to be elucidated. Here we have used an RNA-Seq based study to analyze differential gene expression in G. arboreum under CLCuD infestation. G. arboreum plants were infested by graft inoculation using a CLCuD infected scion of G. hirsutum. CLCuD infested asymptomatic and symptomatic plants were analyzed with RNA-seq using an Illumina HiSeq. 2500. Data analysis revealed 1062 differentially expressed genes (DEGs) in G. arboreum. We selected 17 genes for qPCR to validate RNA-Seq data. We identified several genes involved in disease resistance and pathogen defense. Furthermore, a weighted gene co-expression network was constructed from the RNA-Seq dataset that indicated 50 hub genes, most of which are involved in transport processes and might have a role in the defense response of G. arboreum against CLCuD. This fundamental study will improve the understanding of virus-host interaction and identification of important genes involved in G. arboreum tolerance against CLCuD.
Lee, J. J., Woodward, A. W. & Chen, Z. J. Gene expression changes and early events in cotton fibre development. Annals of botany 100, 1391-1401 (2007).
Chen, Z. J. et al. Toward sequencing cotton (Gossypium) genomes. Plant physiology 145, 1303-1310 (2007).
Paterson, A. H. et al. Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature 492, 423-427 (2012).
Maqbool, A. et al. Gossypium arboreum GHSP26 enhances drought tolerance in Gossypium hirsutum. Biotechnology progress 26, 21-25 (2010).
Robinson, A. F. Reniform in US Cotton: When, Where, Why, and Some Remedies 1. Annual review of phytopathology 45, 263-288 (2007).
Li, R. & Erpelding, J. E. Genetic diversity analysis of Gossypium arboreum germplasm accessions using genotyping-by-sequencing. Genetica 144, 535-545 (2016).
Gong, Q. et al. Salicylic acid-related cotton (Gossypium arboreum) ribosomal protein GaRPL18 contributes to resistance to Verticillium dahliae. BMC plant biology 17, 59 (2017).
Razaq, M., Abbas, G., Farooq, M., Aslam, M. & Athar, H.-U.-R. Effect of Insecticidal Application on Aphid Population, Photosynthetic Parameters and Yield Components of Late Sown Varieties of Canola, Brassica napus L. Pakistan journal of zoology 46, 661-668 (2014).
Zubair, M., Zaidi, S. S., Shakir, S., Amin, I. & Mansoor, S. An insight into Cotton leaf curl Multan betasatellite, the most important component of cotton leaf curl disease complex. Viruses 9, 280 (2017).
Zaidi, S. S., Briddon, R. W. & Mansoor, S. Engineering dual Begomovirus-Bemisia tabaci resistance in plants. Trends in Plant Science 22, 6-8 (2017).
Zaidi, S. S., Martin, D. P., Amin, I., Farooq, M. & Mansoor, S. Tomato leaf curl New Delhi virus: a widespread bipartite begomovirus in the territory of monopartite begomoviruses. Molecular Plant Pathology 18, 901-911 (2017).
Zaidi, S. S.-e.-A., Tashkandi, M., Mansoor, S. & Mahfouz, M. M. Engineering Plant Immunity: Using CRISPR/Cas9 to Generate Virus Resistance. Frontiers in plant science 7 (2016).
Zaidi, S. S. et al. Frequent occurrence of Tomato leaf curl New Delhi virus in cotton leaf curl disease affected cotton in Pakistan. Plos one 11, e0155520 (2016).
Zubair, M. et al. Multiple begomoviruses found associated with cotton leaf curl disease in Pakistan in early 1990 are back in cultivated cotton. Scientific reports 7, 680 (2017).
Gatehouse, J. A. Plant resistance towards insect herbivores: a dynamic interaction. New phytologist 156, 145-169 (2002).
Park, W., Scheffler, B. E., Bauer, P. J. & Campbell, B. T. Genome-wide identification of differentially expressed genes under water deficit stress in upland cotton (Gossypium hirsutum L.). BMC plant biology 12, 90 (2012).
Wang, L. et al. RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biology 14, 169 (2014).
Heidel-Fischer, H. M., Musser, R. O. & Vogel, H. Plant transcriptomic responses to herbivory. Annual plant reviews 47, 155-196 (2014).
Thompson, G. A. & Goggin, F. L. Transcriptomics and functional genomics of plant defence induction by phloem-feeding insects. Journal of experimental botany 57, 755-766 (2006).
Artico, S. et al. Transcriptome analysis of Gossypium hirsutum flower buds infested by cotton boll weevil (Anthonomus grandis) larvae. BMC genomics 15, 854 (2014).
Dubey, N. K. et al. Comparative transcriptome analysis of Gossypium hirsutum L. in response to sap sucking insects: aphid and whitefly. BMC genomics 14, 241 (2013).
Li, J. et al. Transcriptome analysis reveals a comprehensive insect resistance response mechanism in cotton to infestation by the phloem feeding insect Bemisia tabaci (whitefly). Plant biotechnology journal (2016).
Babu, M., Gagarinova, A. G., Brandle, J. E. & Wang, A. Association of the transcriptional response of soybean plants with soybean mosaic virus systemic infection. Journal of general virology 89, 1069-1080 (2008).
Zhang, F. et al. Genetic regulation of salt stress tolerance revealed by RNA-Seq in cotton diploid wild species, Gossypium davidsonii. Scientific reports 6, 20582, https://doi.org/10.1038/srep20582 (2016).
Du, P. et al. Viral infection induces expression of novel phased microRNAs from conserved cellular microRNA precursors. PLoS pathogens 7, e1002176 (2011).
Lu, J. et al. Transcriptome analysis of Nicotiana tabacum infected by Cucumber mosaic virus during systemic symptom development. PLoS one 7, e43447 (2012).
Yang, C. et al. Spatial analysis of Arabidopsis thaliana gene expression in response to Turnip mosaic virus infection. Molecular plant-microbe interactions 20, 358-370 (2007).
Yang, J., Zhang, F., Li, J., Chen, J.-P. & Zhang, H.-M. Integrative analysis of the microRNAome and transcriptome illuminates the response of susceptible rice plants to rice stripe virus. PLoS one 11, e0146946 (2016).
Zhang, Z. et al. De novo transcriptome sequencing in Frankliniella occidentalis to identify genes involved in plant virus transmission and insecticide resistance. Genomics 101, 296-305 (2013).
Akhtar, K. P. et al. First symptomatic evidence of infection of Gossypium arboreum with Cotton leaf curl Burewala virus through grafting. Internationa Journal of Agricultural Biology 15, 157-160 (2013).
Akhtar, K. et al. Evaluation of Gossypium species for resistance to cotton leaf curl Burewala virus. Annals of applied biology 157, 135-147 (2010).
Ullah, R. et al. An analysis of the resistance of Gossypium arboreum to cotton leaf curl disease by grafting. European journal of plant pathology 139, 837-847 (2014).
Azhar, M. T., Amin, I., Bashir, A. & Mansoor, S. Characterization of resistance gene analogs from Gossypium arboreum and their evolutionary relationships with homologs from tetraploid cottons. Euphytica 178, 351-362 (2011).
Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature protocols 7, 562-578 (2012).
Artico, S., Nardeli, S. M., Brilhante, O., Grossi-de-Sa, M. F. & Alves-Ferreira, M. Identification and evaluation of new reference genes in Gossypium hirsutum for accurate normalization of real-time quantitative RT-PCR data. BMC plant biology 10, 49 (2010).
Henry, E., Yadeta, K. A. & Coaker, G. Recognition of bacterial plant pathogens: local, systemic and transgenerational immunity. New phytologist 199, 908-915 (2013).
Zhang, C. et al. Overexpression of a novel peanut NBS-LRR gene AhRRS5 enhances disease resistance to Ralstonia solanacearum in tobacco. Plant biotechnology journal (2016).
Lim, H.-S. et al. The coat protein of Alternanthera mosaic virus is the elicitor of a temperature-sensitive systemic necrosis in Nicotiana benthamiana, and interacts with a host boron transporter protein. Virology 452, 264-278 (2014).
Chandran, D. Co-option of developmentally regulated plant SWEET transporters for pathogen nutrition and abiotic stress tolerance. IUBMB life 67, 461-471 (2015).
Antony, G. et al. Rice xa13 recessive resistance to bacterial blight is defeated by induction of the disease susceptibility gene Os-11N3. The plant cell 22, 3864-3876 (2010).
Wang, Y. & Wu, W.-H. Plant sensing and signaling in response to K +-deficiency. Molecular plant 3, 280-287 (2010).
Zhou, L. et al. Overexpression of GmAKT2 potassium channel enhances resistance to soybean mosaic virus. BMC plant biology 14, 154 (2014).
Zou, J. et al. Expression profiling soybean response to Pseudomonas syringae reveals new defense-related genes and rapid HR-specific downregulation of photosynthesis. Molecular plant-microbe interactions 18, 1161-1174 (2005).
Kim, M. J., Kim, H. R. & Paek, K.-H. Arabidopsis tonoplast proteins TIP1 and TIP2 interact with the cucumber mosaic virus 1a replication protein. Journal of general virology 87, 3425-3431 (2006).
Katiyar, A. et al. Genome-wide classification and expression analysis of MYB transcription factor families in rice and Arabidopsis. BMC genomics 13, 544 (2012).
Wang, J. et al. Genome-wide analysis of bHLH transcription factor and involvement in the infection by yellow leaf curl virus in tomato (Solanum lycopersicum). BMC genomics 16, 39 (2015).
Huang, P.-Y., Catinot, J. & Zimmerli, L. Ethylene response factors in Arabidopsis immunity. Journal of experimental botany, erv518 (2015).
Thomma, B. P., Nürnberger, T. & Joosten, M. H. Of PAMPs and effectors: the blurred PTI-ETI dichotomy. The plant cell 23, 4-15 (2011).
Chen, J.-Y. et al. Genome-wide analysis of the gene families of resistance gene analogues in cotton and their response to Verticillium wilt. BMC plant biology 15, 148 (2015).
Rowland, O. et al. Functional analysis of Avr9/Cf-9 rapidly elicited genes identifies a protein kinase, ACIK1, that is essential for full Cf-9-dependent disease resistance in tomato. The Plant cell 17, 295-310 (2005).
Song, W.-Y., Wang, G.-L., Chen, L.-L. & Kim, H.-S. A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science 270, 1804 (1995).
Kørner, C. J., Du, X., Vollmer, M. E. & Pajerowska-Mukhtar, K. M. Endoplasmic reticulum stress signaling in plant immunity-at the crossroad of life and death. International journal of molecular sciences 16, 26582-26598 (2015).
Chen, M.-H., Tian, G.-W., Gafni, Y. & Citovsky, V. Effects of calreticulin on viral cell-to-cell movement. Plant physiology 138, 1866-1876 (2005).
Stange, C. Plant-virus interactions during the infective process. Cien. Inv. Agr. (in English) 33 (1): 1-18. Ciencia e investigacion agraria 33, 1-18 (2006).
Liu, J. Z. & Whitham, S. A. Overexpression of a soybean nuclear localized type-III DnaJ domain-containing HSP40 reveals its roles in cell death and disease resistance. The plant journal 74, 110-121 (2013).
Taguwa, S. et al. Defining Hsp70 subnetworks in dengue virus replication reveals key vulnerability in flavivirus infection. Cell 163, 1108-1123 (2015).
Gorovits, R. & Czosnek, H. The involvement of heat shock proteins in the establishment of Tomato yellow leaf curl virus infection. Frontiers in plant science 8 (2017).
Citovsky, V., Zaltsman, A., Kozlovsky, S. V., Gafni, Y. & Krichevsky, A. In Seminars in cell & developmental biology. 1048-1054 (Elsevier).
Mandadi, K. K. & Scholthof, K.-B. G. Plant immune responses against viruses: how does a virus cause disease? The plant cell 25, 1489-1505 (2013).
Hua, Z. & Vierstra, R. D. The cullin-RING ubiquitin-protein ligases. Annual review of plant biology 62, 299-334 (2011).
Sánchez-Durán, M. A. et al. Interaction between geminivirus replication protein and the SUMO-conjugating enzyme is required for viral infection. Journal of virology 85, 9789-9800 (2011).
Castillo, A. et al. Identification of Plant Genes Involved in TYLCV Replication. Tomato Yellow Leaf Curl Virus Disease, 207-221 (2007).
Lozano-Durán, R. et al. Geminiviruses subvert ubiquitination by altering CSN-mediated derubylation of SCF E3 ligase complexes and inhibit jasmonate signaling in Arabidopsis thaliana. The Plant Cell Online 23, 1014-1032 (2011).
Eini, O. et al. Interaction with a host ubiquitin-conjugating enzyme is required for the pathogenicity of a geminiviral DNA β satellite. Molecular plant-microbe interactions 22, 737-746 (2009).
Garbutt, C. C., Bangalore, P. V., Kannar, P. & Mukhtar, M. Getting to the edge: protein dynamical networks as a new frontier in plant-microbe interactions. Frontiers in plant science 5 (2014).
Tully, J. P. et al. Expression-based network biology identifies immune-related functional modules involved in plant defense. BMC genomics 15, 421 (2014).
Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC bioinformatics 9, 559 (2008).
Zhang, B. & Horvath, S. A general framework for weighted gene co-expression network analysis. Statistical applications in genetics and molecular biology 4 (2005).
Langfelder, P. & Horvath, S. Eigengene networks for studying the relationships between co-expression modules. BMC systems biology 1, 54 (2007).
Wang, J., Chen, G., Li, M. & Pan, Y. Integration of breast cancer gene signatures based on graph centrality. BMC systems biology 5, S10 (2011).
Nour-Eldin, H. H. et al. NRT/PTR transporters are essential for translocation of glucosinolate defence compounds to seeds. Nature 488, 531-534 (2012).
Sade, D. et al. Involvement of the hexose transporter gene LeHT1 and of sugars in resistance of tomato to tomato yellow leaf curl virus. Molecular plant 6, 1707-1710 (2013).
Delprato, M. L., Krapp, A. R. & Carrillo, N. Green Light to Plant Responses to Pathogens: The Role of Chloroplast Light-Dependent Signaling in Biotic Stress. Photochemistry and photobiology 91, 1004-1011 (2015).
Liu, X. et al. Reactive oxygen species are involved in plant defense against a gall midge. Plant physiology 152, 985-999 (2010).
Bolton, M. D. Primary metabolism and plant defense-fuel for the fire. Molecular plant-microbe Interactions 22, 487-497 (2009).
Rojas, C. M., Senthil-Kumar, M., Tzin, V. & Mysore, K. Regulation of primary plant metabolism during plant-pathogen interactions and its contribution to plant defense. Frontiers in plant science 5, 17 (2014).
Trouvelot, S. et al. Carbohydrates in plant immunity and plant protection: roles and potential application as foliar sprays. Frontiers in plant science 5 (2014).
Wu, J. & Baldwin, I. T. New insights into plant responses to the attack from insect herbivores. Annual review of genetics 44, 1-24 (2010).
Howe, G. A. & Jander, G. Plant immunity to insect herbivores. Annu. Rev. Plant Biol. 59, 41-66 (2008).
Zaidi, S. S., Tashkandi, M. & Mahfouz, M. M. Engineering molecular immunity against plant viruses. Progress in Molecular Biology and Translation Science 149, 167-186 (2017).
Zaidi, S. S., Mansoor, S., Ali, Z., Tashkandi, M. & Mahfouz, M. M. Engineering plants for geminivirus resistance with CRISPR/Cas9 system. Trends in Plant Science 21, 279-281 (2016).
Ali, Z., Ali, S., Tashkandi, M., Zaidi, S. S. & Mahfouz, M. M. CRISPR/Cas9-mediated immunity to geminiviruses: Differential interference and evasion. Scientific reports 6, 26912 (2016).
Zhong, S. et al. High-throughput illumina strand-specific RNA sequencing library preparation. Cold spring harbor protocols 2011, pdb. prot5652 (2011).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, btu170 (2014).
Andrews, S. FastQC: a quality control tool for high throughput sequence data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (2010).
Li, F. et al. Genome sequence of the cultivated cotton Gossypium arboreum. Nature genetics 46, 567-572 (2014).
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nature protocols 11, 1650-1667 (2016).
Zhang, X. et al. mRNA-seq analysis of the Gossypium arboreum transcriptome reveals tissue selective signaling in response to water stress during seedling stage. PLoS one 8, e54762 (2013).
Hagberg, A., Swart, P. J. & Schult, D. A. Exploring network structure, dynamics, and function using NetworkX. (Los Alamos National Laboratory (LANL), 2008).
Saito, R. et al. A travel guide to Cytoscape plugins. Nature methods 9, 1069-1076 (2012).
Demchak, B. et al. Cytoscape: the network visualization tool for GenomeSpace workflows. F1000Research 3 (2014).
Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y. & Morishima, K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic acids research. 45, D353-D361 (2017).
Marino, G. & Funk, C. Matrix metalloproteinases in plants: a brief overview. Physiologia plantarum 145, 196-202 (2012).
Lacerda, A., Vasconcelos, É. A. R., Pelegrini, P. B. & Grossi-de-Sa, M. F. Antifungal defensins and their role in plant defense. Frontiers in microbiology 5, 116 (2014).
Thomma, B. P., Cammue, B. P. & Thevissen, K. Plant defensins. Planta 216, 193-202 (2002).
Spoel, S. H. & Dong, X. How do plants achieve immunity? Defence without specialized immune cells. Nature reviews immunology 12, 89-100 (2012).
Igarashi, D., Tsuda, K. & Katagiri, F. The peptide growth factor, phytosulfokine, attenuates pattern-triggered immunity. The plant journal 71, 194-204 (2012).
Loivamäki, M. et al. A role for PSK signaling in wounding and microbial interactions in Arabidopsis. Physiologia plantarum 139, 348-357 (2010).
Wang, J. et al. Allosteric receptor activation by the plant peptide hormone phytosulfokine. Nature 525, 265-268 (2015).
Park, C.-J. & Seo, Y.-S. Heat shock proteins: a review of the molecular chaperones for plant immunity. The plant pathology journal 31, 323 (2015).
Verchot, J. Cellular chaperones and folding enzymes are vital contributors to membrane bound replication and movement complexes during plant RNA virus infection. Frontiers in plant science 3, 275 (2012).
Hofius, D. et al. Capsid protein-mediated recruitment of host DnaJ-like proteins is required for Potato virus Y infection in tobacco plants. Journal of virology 81, 11870-11880 (2007).
Rawat, R. et al. REVEILLE1, a Myb-like transcription factor, integrates the circadian clock and auxin pathways. Proceedings of the national academy of sciences 106, 16883-16888 (2009).
Zhong, Y. et al. Digital gene expression analysis of Ponkan Mandarin (Citrus reticulata Blanco) in response to Asia citrus psyllid-vectored Huanglongbing infection. International journal of molecular sciences 17, 1063 (2016).
Maimbo, M., Ohnishi, K., Hikichi, Y., Yoshioka, H. & Kiba, A. Induction of a small heat shock protein and its functional roles in Nicotiana plants in the defense response against Ralstonia solanacearum. Plant physiology 145, 1588-1599 (2007).
Stevenson, C. E. et al. Crystal structure of the MYB domain of the RAD transcription factor from Antirrhinum majus. Proteins: structure, function, and bioinformatics 65, 1041-1045 (2006).
Van Der Vossen, E. et al. An ancient R gene from the wild potato species Solanum bulbocastanum confers broad-spectrum resistance to Phytophthora infestans in cultivated potato and tomato. The plant journal 36, 867-882 (2003).
Song, J. et al. Gene RB cloned from Solanum bulbocastanum confers broad spectrum resistance to potato late blight. Proceedings of the national academy of sciences 100, 9128-9133 (2003).
Bourdenx, B. et al. Overexpression of Arabidopsis ECERIFERUM1 promotes wax very-long-chain alkane biosynthesis and influences plant response to biotic and abiotic stresses. Plant physiology 156, 29-45 (2011).
Meng, D., Walsh, M. & Fricke, W. Rapid changes in root hydraulic conductivity and aquaporin expression in rice (Oryza sativa L.) in response to shoot removal-xylem tension as a possible signal. Annals of botany 118, 809-819 (2016).
Matsumoto, T. et al. Role of the aquaporin PIP1 subfamily in the chilling tolerance of rice. Plant and cell physiology 50, 216-229 (2009).
Tian, S. et al. Plant aquaporin AtPIP1; 4 links apoplastic H2O2 induction to disease immunity pathways. Plant physiology, pp. 01237.02015 (2016).
Nagaraj, S., Senthil-Kumar, M., Ramu, V. S., Wang, K. & Mysore, K. S. Plant ribosomal proteins, RPL12 and RPL19, play a role in nonhost disease resistance against bacterial pathogens. Frontiers in plant science 6, 1192 (2016).
Yang, C., Zhang, C., Dittman, J. D. & Whitham, S. A. Differential requirement of ribosomal protein S6 by plant RNA viruses with different translation initiation strategies. Virology 390, 163-173 (2009).
Figueiredo, A., Monteiro, F. & Sebastiana, M. Subtilisin-like proteases in plant-pathogen recognition and immune priming: a perspective. Frontiers in plant science 5, 739 (2014).
Xia, Y. Proteases in pathogenesis and plant defence. Cellular microbiology 6, 905-913 (2004).
Fagard, M. et al. Nitrogen metabolism meets phytopathology. Journal of experimental botany 65, 5643-5656 (2014).
Yamamoto-Katou, A., Katou, S., Yoshioka, H., Doke, N. & Kawakita, K. Nitrate reductase is responsible for elicitin-induced nitric oxide production in Nicotiana benthamiana. Plant and cell physiology 47, 726-735 (2006).
Modolo, L. V. et al. Decreased arginine and nitrite levels in nitrate reductase-deficient Arabidopsis thaliana plants impair nitric oxide synthesis and the hypersensitive response to Pseudomonas syringae. Plant science 171, 34-40 (2006).
