Transcriptome analysis of Nicotiana benthamiana infected by Tobacco curly shoot virus

Cold suppresses virus accumulation and alters the host transcriptomic response in the turnip mosaic virus-Arabidopsis halleri system

Since plant viruses cause lifelong infections, virus-plant interactions are exposed to large temperature fluctuations in evergreen perennials. In such circumstances, virus-plant interactions are expected to change significantly between the warm and cold seasons. However, few studies have investigated the effects of cold temperatures on virus-plant interactions. Here, we show that in a persistent infection system of the turnip mosaic virus (TuMV)-Arabidopsis halleri, cold temperatures lead to slow viral replication/spreading within the host, attenuated host symptoms, and cold-specific transcriptomic responses. Many differentially expressed genes (DEGs) were detected between virus-inoculated and mock-inoculated plants under warm and cold conditions; however, the sets of DEGs and response timings were temperature-dependent. At cold temperatures, the expression of photosynthesis-related genes decreased in the early stages of infection. However, it recovered to the same level as that in uninfected plants in the later stages. In contrast, the transcriptomic changes under warm conditions suggest that viral infections cause auxin signaling disruption. These responses coincided with the inhibition of host growth. We identified 6 cold- and 38 warm-specific DEGs, which changed their expression in response to TuMV infection under more than half of the conditions for either cold or warm temperatures. Further validation of the putative relationships between transcriptomic and phenotypic responses of the host is required. Our findings on temperature-dependent host responses at both symptomatic and transcriptomic levels help us understand how warm and cold temperatures affect virus-plant interactions in seasonal environments.

Mutating the cysteine residues within the WCCH motif of tobacco curly shoot virus (TbCSV) V2 protein attenuates viral infection and confers cross-protection against infection by severe strain

BACKGROUND

Tobacco curly shoot virus (TbCSV, Begomovirus, Geminiviridae) is a viral agent that poses a serious threat to economically important crops like tomato and tobacco in Southwest China. The geminiviral V2 protein has been demonstrated to play multifunctional roles essential for viral infection and pathogenicity. This study comprehensively investigated the biological importance of five cysteines within the conserved 'WCCH' motif of the TbCSV V2 protein, and evaluated the effects of cysteine mutations on viral pathogenicity and the efficacy for cross-protection against wild-type severe TbCSV train.

RESULTS

Site-directed mutagenesis revealed that among the five cysteines in the V2 protein, C84, C86 and C89 residues are critical for its RNA silencing suppression activity, with C84 being the only residue whose mutation impaired its ability to induce tissue necrosis when expressed by potato virus X (PVX). Mutations at C84, C86 or C89 attenuated viral pathogenicity, and the mild symptoms were accompanied by reduced viral accumulations in Nicotiana benthamiana plants. Notably, pre-inoculation with the attenuated TbCSV mutant C84A in N. benthamiana and tomato plants conferred effective cross-protection against subsequent challenge infection by wild-type TbCSV.

CONCLUSION

The cysteine residues within the 'WCCH' motif of TbCSV V2 protein play distinct roles in its functions and viral pathogenicity. Importantly, the attenuated C84A mutant generated in this study not only offers insights into the pathogenic mechanisms of TbCSV but also holds promise for application in crop protection through cross-protection strategy in the field. © 2025 Society of Chemical Industry.

Interplay between tobacco curly shoot virus vsiRNA24 and triosephosphate isomerase: implications for Nicotiana benthamiana viral defense

Virus-derived small interfering RNAs (vsiRNAs) play an important role in viral infection by regulating the expression of host genes. At present, research on the regulation of plant primary metabolic pathways by vsiRNAs is very limited. TvsiRNA24 derived from tobacco curly shoot virus (TbCSV) was amplified by reverse transcription polymerase chain reaction, and its target gene NbTPI (triosephosphate isomerase) was verified using reverse transcription quantitative polymerase chain reaction and GFP fluorescence observation. The effect of the interaction between TvsiRNA24 and NbTPI on TbCSV infection was analyzed by virus mediated, genetic transformation, western blotting, and quantitative detection. The expression of TvsiRNA24 retards the growth of Nicotiana benthamiana and enhances TbCSV accumulation within N. benthamiana. The overexpression of NbTPI attenuates the accumulation of TbCSV, and the silencing of NbTPI leads to the growth retardation of N. benthamiana and intensifies symptoms post-TbCSV infection. Moreover, the expression of some genes related to photosynthesis, primary metabolism and immune response is regulated by NbTPI. Our results unveil the specific role of TvsiRNA24-NbTPI in the pathogenicity of TbCSV, resulting in hindrance to plant growth and facilitation of viral infection. The identification of this regulatory pathway provides valuable insights that can be utilized to devise novel antiviral approaches targeting the reduction of viral pathogenicity.

Ethylene response transcription factor 5 (ERF5) enhances defense against tobacco curly shoot virus and associated betasatellite (TbCSV/TbCSB) in Nicotiana benthamiana

Begomovirus/betasatellite disease complex significantly threatens global crop production. Identifying potential plant antiviral genes is crucial for disease control. Nicotiana benthamiana is susceptible to viruses and contains 266 ethylene response transcription factors (ERFs). This study identified 29 NbERFs that were differentially upregulated in tobacco curly shoot virus and its associated betasatellite (TbCSV/TbCSB) infection, with ERF5 being the most common. Nine NbERF5s cluster phylogenetically and Niben101Scf00163g22002 (NbERF5) responds significantly to exogenous ACC treatment. Further analysis confirms the nuclear localization and transcriptional activation activity of NbERF5. Protein interaction assays demonstrate that NbERF5 has no self-interaction and does not interact with the βC1 protein of TbCSB. Silencing NbERF5 enhances TbCSV/TbCSB infection, and overexpression of NbERF5 inhibits TbCSV/TbCSB infection. Importantly, NbERF5 positively regulates the expression of the pathogenesis-related (PR) genes, NbPR1a and NbNPR1. Our findings suggest that NbERF5 enhances TbCSV/TbCSB resistance by activating the PR genes, indicating that NbERF5 is a potential antiviral gene.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato

BACKGROUND

Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV).

RESULTS

We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes.

CONCLUSION

We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

Size variations of mesoporous silica nanoparticle control uptake efficiency and delivery of AC2-derived dsRNA for protection against tomato leaf curl New Delhi virus

KEY MESSAGE

We report the size dependent uptake of dsRNA loaded MSNPs into the leaves and roots of Nicotiana benthamiana plants and accessed for their relative reduction in Tomato leaf curl New Delhi viral load. A non-GMO method of RNA interference (RNAi) has been recently in practice through direct delivery of double stranded RNA into the plant cells. Tomato leaf curl New Delhi virus (ToLCNDV), a bipartitie begomovirus, is a significant viral pathogen of many crops in the Indian subcontinent. Conventional RNAi cargo delivery strategies for instance uses viral vectors and Agrobacterium-facilitated delivery, exhibiting specific host responses from the plant system. In the present study, we synthesized three different sizes of amine-functionalized mesoporous silica nanoparticles (amino-MSNPs) to mediate the delivery of dsRNA derived from the AC2 (dsAC2) gene of ToLCNDV and showed that these dsRNA loaded nanoparticles enabled effective reduction in viral load. Furthermore, we demonstrate that amino-MSNPs protected the dsRNA molecules from nuclease degradation, while the complex was efficiently taken up by the leaves and roots of Nicotiana benthamiana. The real time gene expression evaluation showed that plants treated with nanoparticles of different sizes ~ 10 nm (MSNPDEA), ~ 32 nm (MSNPTEA) and ~ 66 nm (MSNPNH3) showed five-, eleven- and threefold reduction of ToLCNDV in N. benthamiana, respectively compared to the plants treated with naked dsRNA. This work clearly demonstrates the size dependent internalization of amino-MSNPs and relative efficacy in transporting dsRNA into the plant system, which will be useful in convenient topical treatment to protect plants against their pathogens including viruses. Mesoporous silica nanoparticles loaded with FITC, checked for its uptake into Nicotiana benthamiana.

Transcriptome Sequencing and WGCNA Reveal Key Genes in Response to Leaf Blight in Poplar

Leaf blight is a fungal disease that mainly affects the growth and development of leaves in plants. To investigate the molecular mechanisms of leaf blight defense in poplar, we performed RNA-Seq and enzyme activity assays on the Populus simonii × Populus nigra leaves inoculated with Alternaria alternate fungus. Through weighted gene co-expression network analysis (WGCNA), we obtained co-expression gene modules significantly associated with SOD and POD activities, containing 183 and 275 genes, respectively. We then constructed a co-expression network of poplar genes related to leaf blight resistance based on weight values. Additionally, we identified hub transcription factors (TFs) and structural genes in the network. The network was dominated by 15 TFs, and four out of them, including ATWRKY75, ANAC062, ATMYB23 and ATEBP, had high connectivity in the network, which might play important functions in leaf blight defense. In addition, GO enrichment analysis revealed a total of 44 structural genes involved in biotic stress, resistance, cell wall and immune-related biological processes in the network. Among them, there were 16 highly linked structural genes in the central part, which may be directly involved in poplar resistance to leaf blight. The study explores key genes associated with leaf blight defense in poplar, which further gains an understanding of the molecular mechanisms of biotic stress response in plants.

Molecular insights into the responses of barley to yellow mosaic disease through transcriptome analysis

BACKGROUND

Barley (Hordeum vulgare L.) represents the fourth most essential cereal crop in the world, vulnerable to barley yellow mosaic virus (BaYMV) and/or barley mild mosaic virus (BaMMV), leading to the significant yield reduction. To gain a better understanding of the mechanisms regarding barley crop tolerance to virus infection, we employed a transcriptome sequencing approach and investigated global gene expression among three barley varieties under both infected and control conditions.

RESULTS

High-throughput sequencing outputs revealed massive genetic responses, reflected by the barley transcriptome after BaYMV and/or BaMMV infection. Significant enrichments in peptidase complex and protein processing in endoplasmic reticulum were clustered through Gene ontology and KEGG analysis. Many genes were identified as transcription factors, antioxidants, disease resistance genes and plant hormones and differentially expressed between infected and uninfected barley varieties. Importantly, general response genes, variety-specific and infection-specific genes were also discovered. Our results provide useful information for future barley breeding to resist BaYMV and BaMMV.

CONCLUSIONS

Our study elucidates transcriptomic adaptations in barley response to BaYMV/BaMMV infection through high-throughput sequencing technique. The analysis outcome from GO and KEGG pathways suggests that BaYMV disease induced regulations in multiple molecular-biology processes and signalling pathways. Moreover, critical DEGs involved in defence and stress tolerance mechanisms were displayed. Further functional investigations focusing on these DEGs contributes to understanding the molecular mechanisms of plant response to BaYMV disease infection, thereby offering precious genetic resources for breeding barley varieties resistant to BaYMV disease.

vsiRNA18 derived from tobacco curly shoot virus can regulate virus infection in Nicotiana benthamiana

Virus-derived small interfering RNAs (vsiRNAs) play important roles in regulating host endogenous gene expression to promote virus infection and induce RNA silencing to suppress virus infection. However, to date, how vsiRNAs affect geminivirus infection in host plants has been less studied. In this study, we found that tobacco curly shoot virus (TbCSV)-derived vsiRNA18 (TvsiRNA18) can regulate TbCSV infection in Nicotiana benthamiana plants. The virus-mediated small RNA expression system and stable transformation technique were used to clarify the molecular role of TvsiRNA18 in TbCSV infection. The results indicate that TvsiRNA18 can aggravate disease symptoms in these plants and enhance viral DNA accumulation. ATP-dependent RNA helicase (ATP-dRH) was proven to be a target of TvsiRNA18, and down-regulation of ATP-dRH in plants was shown to induce virus-like leaf curling symptoms and increase TbCSV infection. These results suggest that TvsiRNA18 is an important regulator of TbCSV infection by suppressing ATP-dRH expression. This is the first report to demonstrate that TbCSV-derived vsiRNA can target host endogenous genes to affect symptom development, which helps to reveal the molecular mechanism of symptom occurrence after the virus infects the host.

Plant Defense and Viral Counter-Defense during Plant-Geminivirus Interactions

Geminiviruses are the largest family of plant viruses that cause severe diseases and devastating yield losses of economically important crops worldwide. In response to geminivirus infection, plants have evolved ingenious defense mechanisms to diminish or eliminate invading viral pathogens. However, increasing evidence shows that geminiviruses can interfere with plant defense response and create a suitable cell environment by hijacking host plant machinery to achieve successful infections. In this review, we discuss recent findings about plant defense and viral counter-defense during plant-geminivirus interactions.

Transcriptome mapping related genes encoding PR1 protein involved in necrotic symptoms to soybean mosaic virus infection

Necrosis caused by soybean mosaic virus (SMV) has not been specifically distinguished from susceptible symptoms. The molecular mechanism for the occurrence of necrosis is largely overlooked in soybean genetic research. Field evaluation reveals that SMV disease seriously influences soybean production as indicated by decreasing 22.4% ~ 77.0% and 8.8% ~ 17.0% of yield and quality production, respectively. To expand molecular mechanism behind necrotic reactions, transcriptomic data obtained from the asymptomatic, mosaic, and necrotic pools were assessed. Compared between asymptomatic and mosaic plants, 1689 and 1752 up- and down-regulated differentially expressed genes (DEGs) were specifically found in necrotic plants. Interestingly, the top five enriched pathways with up-regulated DEGs were highly related to the process of the stress response, whereas the top three enriched pathways with down-regulated DEGs were highly related to the process of photosynthesis, demonstrating that defense systems are extensively activated, while the photosynthesis systems were severely destroyed. Further, results of the phylogenetic tree based on gene expression pattern and an amino acid sequence and validation experiments discovered three PR1 genes, Glyma.15G062400, Glyma.15G062500, and Glyma.15G062700, which were especially expressed in necrotic leaves. Meanwhile, exogenous salicylic acid (SA) but not methyl jasmonate (MeJA) could induce the three PR1 gene expressions on healthy leaves. Contrastingly, exogenous SA obviously decreased the expression level of Glyma.15G062400, Glyma.15G062500, and concentration of SMV, but increased Glyma.15G062700 expression in necrotic leaves. These results showed that GmPR1 is associated with the development of SMV-induced necrotic symptoms in soybean. Glyma.15G062400, Glyma.15G062500, and Glyma.15G062700 is up-regulated in necrotic leaves at the transcriptional levels, which will greatly facilitate a better understanding of the mechanism behind necrosis caused by SMV disease.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-022-01351-3.

Transcriptomic alterations in the sweet orange vasculature correlate with growth repression induced by a variant of citrus tristeza virus

Citrus tristeza virus (CTV, family Closteroviridae) is an economically important pathogen of citrus. CTV resides in the phloem of the infected plants and induces a range of disease phenotypes, including stem pitting and quick decline as well as a number of other deleterious syndromes. To uncover the biological processes underlying the poorly understood damaging symptoms of CTV, we profiled the transcriptome of sweet orange (Citrus sinensis) phloem-rich bark tissues of non-infected, mock-inoculated trees and trees singly infected with two distinct variants of CTV, T36 or T68-1. The T36 and T68-1 variants accumulated in the infected plants at similar titers. With that, young trees infected with T68-1 were markedly repressed in growth, while the growth rate of the trees infected with T36 was comparable to the mock-inoculated trees. Only a small number of differentially expressed genes (DEGs) were identified in the nearly asymptomatic T36-infected trees, whereas almost fourfold the number of DEGs were identified with the growth-restricting T68-1 infection. DEGs were validated using quantitative reverse transcription-PCR. While T36 did not induce many noteworthy changes, T68-1 altered the expression of numerous host mRNAs encoding proteins within significant biological pathways, including immunity and stress response proteins, papain-like cysteine proteases (PLCPs), cell-wall modifying enzymes, vascular development proteins and others. The transcriptomic alterations in the T68-1-infected trees, in particular, the strong and persistent increase in the expression levels of PLCPs, appear to contribute to the observed stem growth repression. On the other hand, analysis of the viral small interfering RNAs revealed that the host RNA silencing-based response to the infection by T36 and that by T68-1 was comparable, and thus, the induction of this antiviral mechanism may not contribute to the difference in the observed symptoms. The DEGs identified in this study promote our understanding of the underlying mechanisms of the yet unexplained growth repression induced by severe CTV isolates in sweet orange trees.

Integrated next-generation sequencing and comparative transcriptomic analysis of leaves provides novel insights into the ethylene pathway of Chrysanthemum morifolium in response to a Chinese isolate of chrysanthemum virus B

Background

Chrysanthemum virus B (CVB), a key member of the genus Carlavirus, family Betaflexiviridae, causes severe viral diseases in chrysanthemum (Chrysanthemum morifolium) plants worldwide. However, information on the mechanisms underlying the response of chrysanthemum plants to CVB is scant.

Methods

Here, an integrated next-generation sequencing and comparative transcriptomic analysis of chrysanthemum leaves was conducted to explore the molecular response mechanisms of plants to a Chinese isolate of CVB (CVB-CN) at the molecular level.

Results

In total, 4934 significant differentially expressed genes (SDEGs) were identified to respond to CVB-CN, of which 4097 were upregulated and 837 were downregulated. Gene ontology and functional classification showed that the majority of upregulated SDEGs were categorized into gene cohorts involved in plant hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, and ribosome metabolism. Enrichment analysis demonstrated that ethylene pathway-related genes were significantly upregulated following CVB-CN infection, indicating a strong promotion of ethylene biosynthesis and signaling. Furthermore, disruption of the ethylene pathway in Nicotiana benthamiana, a model plant, using virus-induced gene silencing technology rendered them more susceptible to cysteine-rich protein of CVB-CN induced hypersensitive response, suggesting a crucial role of this pathway in response to CVB-CN infection.

Conclusion

This study provides evidence that ethylene pathway has an essential role of plant in response to CVB and offers valuable insights into the defense mechanisms of chrysanthemum against Carlavirus.

Transcriptome analyses unveiled differential regulation of AGO and DCL genes by pepino mosaic virus strains

Pepino mosaic virus (PepMV) is a single-stranded (ss), positive-sense (+) RNA potexvirus that affects tomato crops worldwide. We have described an in planta antagonistic interaction between PepMV isolates of two strains in which the EU isolate represses the accumulation of the CH2 isolate during mixed infections. Reports describing transcriptomic responses to mixed infections are scant. We carried out transcriptomic analyses of tomato plants singly and mixed-infected with two PepMV isolates of both strains. Comparison of the transcriptomes of singly infected plants showed that deeper transcriptomic alterations occurred at early infection times, and also that each of the viral strains modulated the host transcriptome differentially. Mixed infections caused transcriptomic alterations similar to those for the sum of single infections at early infection times, but clearly differing at later times postinfection. We next tested the hypothesis that PepMV-EU, in either single or mixed infections, deregulates host gene expression differentially so that virus accumulation of both strains gets repressed. That seemed to be the case for the genes AGO1a, DCL2d, AGO2a, and DCL2b, which are involved in the antiviral silencing pathway and were upregulated by PepMV-EU but not by PepMV-CH2 at early times postinfection. The pattern of AGO2a expression was validated by reverse transcription-quantitative PCR in tomato and Nicotiana benthamiana plants. Using an N. benthamiana ago2 mutant line, we showed that AGO2 indeed plays an important role in the antiviral defence against PepMV, but it is not the primary determinant of the outcome of the antagonistic interaction between the two PepMV strains.

Integrated transcriptomic and transgenic analyses reveal potential mechanisms of poplar resistance to Alternaria alternata infection

BACKGROUND

Populus davidiana × P. bollena is a species of poplar from northeastern China that is characterized by cold resistance and fast growth but now suffers from pathogen infections. Leaf blight caused by Alternaria alternata has become a common poplar disease that causes serious economic impacts, but the molecular mechanisms of resistance to A. alternata in P. davidiana × P. bollena are still unclear.

RESULTS

In this study, the transcriptomic response of P. davidiana × P. bollena to A. alternata infection was determined via RNA-Seq. Twelve cDNA libraries were generated from RNA isolated from three biological replicates at four time points (0, 2, 3, and 4 d post inoculation), and a total of 5,930 differentially expressed genes (DEGs) were detected (| log2 fold change |≥ 1 and FDR values < 0.05). Functional analysis revealed that the DEGs were mainly enriched for the "plant hormone signal transduction" pathway, followed by the "phenylpropanoid biosynthesis" pathway. In addition, DEGs that encode defense-related proteins and are related to ROS metabolism were also identified. Numerous transcription factors, such as the bHLH, WRKY and MYB families, were also induced by A. alternata infection. Among these DEGs, those related to JA biosynthesis and JA signal transduction were consistently activated. Therefore, the lipoxygenase gene PdbLOX2, which is involved in JA biosynthesis, was selected for functional characterization. Overexpression of PdbLOX2 enhanced the resistance of P. davidiana × P. bollena to A. alternata, whereas silencing this gene enhanced susceptibility to A. alternata infection.

CONCLUSIONS

These results provide new insight into the molecular mechanisms of poplar resistance to A. alternata infection and provide candidate genes for breeding resistant cultivars using genetic engineering.

Tomato zonate spot virus induced hypersensitive resistance via an auxin-related pathway in pepper

Tomato zonate spotvirus (TZSV) often incurs significant losses in many food and ornamental crops in Yunnan province, China, and the surrounding areas. The pepper (Capsicum chinensePI152225)can develop hypersensitive resistance following infection with TZSV, through an as yet unknown mechanism. The transcriptome dataset showed a total of 45.81 GB of clean data were obtained from six libraries, and the average percentage of the reads mapped to the pepper genome was over 90.00 %. A total of 1403 differentially expressed genes (DEGs) were obtained after TZSV infection, including 825significantly up-regulated genes and 578 down-regulated genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that most up-regulated DEGs were involved in basal defenses. RT-qPCR, and virus induced gene silencing (VIGS) were used preliminarily to identifyBBC_22506 and BBC_18917, among total of 71 differentially expressed genes (DEGs), that play a key role in mediating the auxin-induced signaling pathway that might take part in hypersensitive response (HR) conferred resistance to viral infection in pepper (PI152225) byTZSV. This is the first study on the mechanism of auxin resistance, involved in defense responses of pepper against viral diseases, which lay the foundation for further study on the pathogenic mechanism of TZSV, as well as the mechanism of resistance to TZSV, in peppers.

The interplay of plant hormonal pathways and geminiviral proteins: partners in disease development

Viruses belonging to the family Geminiviridae infect plants and are responsible for a number of diseases of crops in the tropical and sub-tropical regions of the World. The innate immune response of the plant assists in its defense against such viral pathogens by the recognition of pathogen/microbe-associated molecular patterns through pattern-recognition receptors. Phytohormone signalling pathways play a vital role in plant defense responses against these devastating viruses. Geminiviruses, however, have developed counter-defense strategies that prevail over the above defense pathways. The proteins encoded by geminiviruses act as suppressors of plant immunity by interacting with the signalling components of several hormones. In this review we focus on the molecular interplay of phytohormone pathways and geminiviral infection and try to find interesting parallels with similar mechanisms known in other plant-infecting viruses and strengthen the argument that this interplay is necessary for disease development.

Effects of the noncoding subgenomic RNA of red clover necrotic mosaic virus in virus infection

In recent years, a new class of viral noncoding subgenomic RNA (ncsgRNA) has been identified. This RNA is generated as a stable degradation product via an exoribonuclease-resistant RNA (xrRNA) structure, which blocks the progression of 5′→3′ exoribonuclease on viral RNAs in infected cells. Here, we assess the effects of the ncsgRNA of red clover necrotic mosaic virus (RCNMV), called SR1f, in infected plants. We demonstrate the following: (i) the absence of SR1f reduces symptoms and decreases viral RNA accumulation in Nicotiana benthamiana and Arabidopsis thaliana plants; (ii) SR1f has an essential function other than suppression of RNA silencing; and (iii) the cytoplasmic exoribonuclease involved in mRNA turnover, XRN4, is not required for SR1f production or virus infection. A comparative transcriptomic analysis in N. benthamiana infected with wild-type RCNMV or an SR1f-deficient mutant RCNMV revealed that wild-type RCNMV infection, which produces SR1f and much higher levels of virus, has a greater and more significant impact on cellular gene expression than the SR1f-deficient mutant. Upregulated pathways include plant hormone signaling, plant-pathogen interaction, MAPK signaling, and several metabolic pathways, while photosynthesis-related genes were downregulated. We compare this to host genes known to participate in infection by other tombusvirids. Viral reads revealed a 10- to 100-fold ratio of positive to negative strand, and the abundance of reads of both strands mapping to the 3′ region of RCNMV RNA1 support the premature transcription termination mechanism of synthesis for the coding sgRNA. These results provide a framework for future studies of the interactions and functions of noncoding RNAs of plant viruses.

IMPORTANCE Knowledge of how RNA viruses manipulate host and viral gene expression is crucial to our understanding of infection and disease. Unlike viral protein-host interactions, little is known about the control of gene expression by viral RNA. Here, we begin to address this question by investigating the noncoding subgenomic RNA (ncsgRNA) of red clover necrotic mosaic virus (RCNMV), called SR1f. Similar exoribonuclease-resistant RNAs of flaviviruses are well studied, but the roles of plant viral ncsgRNAs, and how they arise, are poorly understood. Surprisingly, we find the likely exonuclease candidate, XRN4, is not required to generate SR1f, and we assess the effects of SR1f on virus accumulation and symptom development. Finally, we compare the effects of infection by wild-type RCNMV versus an SR1f-deficient mutant on host gene expression in Nicotiana benthamiana, which reveals that ncsgRNAs such as SR1f are key players in virus-host interactions to facilitate productive infection.

RNA sequencing reveals transcriptomic changes in tobacco (Nicotiana tabacum) following NtCPS2 knockdown

BACKGROUND

Amber-like compounds form in tobacco (Nicotiana tabacum) during leaf curing and impact aromatic quality. In particular, cis-abienol, a polycyclic labdane-related diterpenoid, is of research interest as a precursor of these compounds. Glandular trichome cells specifically express copalyl diphosphate synthase (NtCPS2) at high levels in tobacco, which, together with NtABS, are major regulators of cis-abienol biosynthesis in tobacco.

RESULTS

To identify the genes involved in the biosynthesis of cis-abienol in tobacco, we constructed transgenic tobacco lines based on an NtCPS2 gene-knockdown model using CRISPR/Cas9 genome-editing technology to inhibit NtCPS2 function in vitro. In mutant plants, cis-abienol and labdene diol contents decreased, whereas the gibberellin and abscisic acid (ABA) contents increased compared with those in wild-type tobacco plants. RNA sequencing analysis revealed the presence of 9514 differentially expressed genes (DEGs; 4279 upregulated, 5235 downregulated) when the leaves of wild-type and NtCPS2-knockdown tobacco plants were screened. Among these DEGs, the genes encoding cis-abienol synthase, ent-kaurene oxidase, auxin/ABA-related proteins, and transcription factors were found to be involved in various biological and physiochemical processes, including diterpenoid biosynthesis, plant hormone signal transduction, and plant-pathogen interactions.

CONCLUSIONS

The present study provides insight into the unique transcriptome profile of NtCPS2 knockdown tobacco, allowing for a better understanding of the biosynthesis of cis-abienol in tobacco.

Comparative RNA-Seq analysis unfolds a complex regulatory network imparting yellow mosaic disease resistance in mungbean [Vigna radiata (L.) R. Wilczek]

Yellow Mosaic Disease (YMD) in mungbean [Vigna radiata (L.) R. Wilczek] is one of the most damaging diseases in Asia. In the northern part of India, the YMD is caused by Mungbean Yellow Mosaic India Virus (MYMIV), while in southern India this is caused by Mungbean Yellow Mosaic Virus (MYMV). The molecular mechanism of YMD resistance in mungbean remains largely unknown. In this study, RNA-seq analysis was conducted between a resistant (PMR-1) and a susceptible (Pusa Vishal) mungbean genotype under infected and control conditions to understand the regulatory network operating between mungbean-YMV. Overall, 76.8 million raw reads could be generated in different treatment combinations, while mapping rate per library to the reference genome varied from 86.78% to 93.35%. The resistance to MYMIV showed a very complicated gene network, which begins with the production of general PAMPs (pathogen-associated molecular patterns), then activation of various signaling cascades like kinases, jasmonic acid (JA) and brassinosteroid (BR), and finally the expression of specific genes (like PR-proteins, virus resistance and R-gene proteins) leading to resistance response. The function of WRKY, NAC and MYB transcription factors in imparting the resistance against MYMIV could be established. The string analysis also revealed the role of proteins involved in kinase, viral movement and phytoene synthase activity in imparting YMD resistance. A set of novel stress-related EST-SSRs are also identified from the RNA-Seq data which may be used to find the linked genes/QTLs with the YMD resistance. Also, 11 defence-related transcripts could be validated through quantitative real-time PCR analysis. The identified gene networks have led to an insight about the defence mechanism operating against MYMIV infection in mungbean which will be of immense use to manage the YMD resistance in mungbean.

Transcriptomic and functional analyses reveal an antiviral role of autophagy during pepper mild mottle virus infection

BACKGROUND

Pepper mild mottle virus (PMMoV) is a member in the genus Tobamovirus and infects mainly solanaceous plants. However, the mechanism of virus-host interactions remains unclear. To explore the responses of pepper plants to PMMoV infection, we analyzed the transcriptomic changes in pepper plants after PMMoV infection using a high-throughput RNA sequencing approach and explored the roles of host autophagy in regulating PMMoV infection.

RESULTS

A total of 197 differentially expressed genes (DEGs) were obtained after PMMoV infection, including 172 significantly up-regulated genes and 25 down-regulated genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that most up-regulated DEGs were involved in plant abiotic and biotic stresses. Further analyses showed the expressions of multiple autophagy-related genes (ATGs) were increased after PMMoV infection in pepper and Nicotiana benthamiana plants. Through confocal microscopy and transmission electron microscopy, we have found that PMMoV infection in plant can induce autophagy, evidenced by the increased number of GFP-ATG8a fluorescent punctate and the appearance of double membrane autophagic structures in cells of N. benthamiana. Additionally, inhibition of autophagy significantly increased PMMoV RNA accumulation and aggravated systemic PMMoV symptoms through autophagy inhibitor (3-MA and E64d) treatment and silencing of NbATG expressions by a Tobacco rattle virus-induced gene silencing assays. These results indicated that autophagy played a positive role in plant resistance to PMMoV infection.

CONCLUSIONS

Taken together, our results provide a transcriptomic insight into pepper responding to PMMoV infection and reveal that autophagy induced by PMMoV infection has an antiviral role in regulating PMMoV infection. These results also help us to better understand the mechanism controlling PMMoV infection in plants and to develop better strategies for breeding projects for virus-resistant crops.

Tobacco curly shoot virus C3 protein enhances viral replication and gene expression in Nicotiana benthamiana plants

Geminiviruses are single-stranded DNA viruses that cause devastating diseases in many crops worldwide. The replication enhancer proteins (REn), encoded by the C3 (AC3, and AL3) ORFs of geminiviruses, have critical roles in viral DNA accumulation and symptom development in infected plants. In the current study, we have constructed an infectious clone of the Tobacco curly shoot virus (TbCSV) C3 mutant, TbCSVΔC3, that contains two start codon mutations that abrogated C3 ORF expression, but did not alter the amino acid sequence of the C2 ORF. As predicted, the absence of the C3 protein reduced TbCSV DNA accumulation, and over-expression of the C3 protein enhanced TbCSV DNA accumulation in infected leaves of Nicotiana benthamiana. The C3 mutation reduced the expression levels of both virion- and complementary-sense TbCSV genes whereas over-expression of the C3 protein increased TbCSV gene expression. Furthermore, the expression of the wild-type and site-directed mutants of C3 proteins using the potato virus X (PVX) system showed that Y93A mutation reduced the replication enhancement activity of the C3 protein in N. benthamiana. All the available evidence demonstrates that the C3 protein is tightly coupled with TbCSV DNA accumulation. However, the TbCSVΔC3 mutant was nearly as infectious in N. benthamiana as TbCSVWT and only had slightly delayed and attenuated symptom expression. Our findings demonstrate that TbCSV C3 protein enhances viral replication and gene expression, but has only moderate effects on symptom development in N. benthamiana.

Diversity, Function and Regulation of Cell Surface and Intracellular Immune Receptors in Solanaceae

The first layer of the plant immune system comprises plasma membrane-localized receptor proteins and intracellular receptors of the nucleotide-binding leucine-rich repeat protein superfamily. Together, these immune receptors act as a network of surveillance machines in recognizing extracellular and intracellular pathogen invasion-derived molecules, ranging from conserved structural epitopes to virulence-promoting effectors. Successful pathogen recognition leads to physiological and molecular changes in the host plants, which are critical for counteracting and defending against biotic attack. A breadth of significant insights and conceptual advances have been derived from decades of research in various model plant species regarding the structural complexity, functional diversity, and regulatory mechanisms of these plant immune receptors. In this article, we review the current state-of-the-art of how these host surveillance proteins function and how they are regulated. We will focus on the latest progress made in plant species belonging to the Solanaceae family, because of their tremendous importance as model organisms and agriculturally valuable crops.

Suppression of nbe-miR1919c-5p Expression in Nicotiana benthamiana Enhances Tobacco Curly Shoot Virus and Its Betasatellite Co-Infection

MicroRNAs (miRNAs) are non-coding but functional RNA molecules of 21-25 nucleotides in length. MiRNAs play significant regulatory roles in diverse plant biological processes. In order to decipher the relationship between nbe-miR1919c-5p and the accumulations of tobacco curly shoot virus (TbCSV) and its betasatellite (TbCSB) DNAs, as well as viral symptom development, we investigated the function of nbe-miR1919c-5p during TbCSV and TbCSB co-infection in plants using a PVX-and a TRV-based short tandem target mimic (STTM) technology. Suppression of nbe-miR1919c-5p expression using these two technologies enhanced TbCSV and TbCSB co-infection-induced leaf curling symptoms in Nicotiana benthamiana plants. Furthermore, suppression of nbe-miR1919c-5p expression enhanced TbCSV and TbCSB DNA accumulations in the infected plants. Our results can advance our knowledge on the nbe-miR1919c-5p function during TbCSV and TbCSB co-infection.

Transcriptome reprogramming in the shoot apical meristem of CymRSV-infected Nicotiana benthamiana plants associates with viral exclusion and the lack of recovery

In some plant-virus interactions plants show a sign of healing from virus infection, a phenomenon called symptom recovery. It is assumed that the meristem exclusion of the virus is essential to this process. The discovery of RNA silencing provided a possible mechanism to explain meristem exclusion and recovery. Here we show evidence that silencing is not the reason for meristem exclusion in Nicotiana benthamiana plants infected with Cymbidium ringspot virus (CymRSV). Transcriptome analysis followed by in situ hybridization shed light on the changes in gene expression in the shoot apical meristem (SAM) on virus infection. We observed the down-regulation of meristem-specific genes, including WUSCHEL (WUS). However, WUS was not down-regulated in the SAM of plants infected with meristem-invading viruses such as turnip vein-clearing virus (TVCV) and cucumber mosaic virus (CMV). Moreover, there is no connection between loss of meristem function and fast shoot necrosis since TVCV necrotized the shoot while CMV did not. Our findings suggest that the observed transcriptional changes on virus infection in the shoot are key factors in tip necrosis and symptom recovery. We observed a lack of GLYCERALDEHYDE 3-PHOSPHATE DEHYDROGENASE (GAPDH) expression in tissues around the meristem, which likely stops virus replication and spread into the meristem.

Identification of microRNAs regulated by tobacco curly shoot virus co-infection with its betasatellite in Nicotiana benthamiana

BACKGROUND

MicroRNAs (miRNAs) are a class of 21-24 nucleotide endogenous non-coding small RNAs that play important roles in plant development and defense responses to biotic and abiotic stresses. Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus, cause leaf curling and plant stunting symptoms in many Solanaceae plants. The betasatellite of TbCSV (TbCSB) induces more severe symptoms and enhances virus accumulation when co-infect the plants with TbCSV.

METHODS

In this study, miRNAs regulated by TbCSV and TbCSB co-infection in Nicotiana benthamiana were characterized using high-throughput sequencing technology.

RESULTS

Small RNA sequencing analysis revealed that a total of 13 known miRNAs and 42 novel miRNAs were differentially expressed in TbCSV and TbCSB co-infected N. benthamiana plants. Several potential miRNA-targeted genes were identified through data mining and were involved in both catalytic and metabolic processes, in addition to plant defense mechanisms against virus infections according to Gene Ontology (GO) analyses. In addition, the expressions of several differentially expressed miRNAs and their miRNA-targeted gene were validated through quantitative real time polymerase chain reaction (qRT-PCR) approach.

CONCLUSIONS

A large number of miRNAs are identified, and their target genes, functional annotations also have been explored. Our results provide the information on N. benthamiana miRNAs and would be useful to further understand miRNA regulatory mechanisms after TbCSV and TbCSB co-infection.

Phenotyping viral infection in sweetpotato using a high-throughput chlorophyll fluorescence and thermal imaging platform

BACKGROUND

Virus diseases caused by co-infection with Sweet potato feathery mottle virus (SPFMV) and Sweetpotato chlorotic stunt virus (SPCSV) are a severe problem in the production of sweetpotato (Ipomoea batatas L.). Traditional molecular virus detection methods include nucleic acid-based and serological tests. In this study, we aimed to validate the use of a non-destructive imaging-based plant phenotype platform to study plant-virus synergism in sweetpotato by comparing four virus treatments with two healthy controls.

RESULTS

By monitoring physiological and morphological effects of viral infection in sweetpotato over 29 days, we quantified photosynthetic performance from chlorophyll fluorescence (ChlF) imaging and leaf thermography from thermal infrared (TIR) imaging among sweetpotatoes. Moreover, the differences among different treatments observed from ChlF and TIR imaging were related to virus accumulation and distribution in sweetpotato. These findings were further validated at the molecular level by related gene expression in both photosynthesis and carbon fixation pathways.

CONCLUSION

Our study validated for the first time the use of ChlF- and TIR-based imaging systems to distinguish the severity of virus diseases related to SPFMV and SPCSV in sweetpotato. In addition, we demonstrated that the operating efficiency of PSII and photochemical quenching were the most sensitive parameters for the quantification of virus effects compared with maximum quantum efficiency, non-photochemical quenching, and leaf temperature.

Molecular Regulation of Host Defense Responses Mediated by Biological Anti-TMV Agent Ningnanmycin

Ningnanmycin (NNM) belongs to microbial pesticides that display comprehensive antiviral activity against plant viruses. NNM treatment has been shown to efficiently delay or suppress the disease symptoms caused by tobacco mosaic virus (TMV) infection in local-inoculated or systemic-uninoculated tobacco leaves, respectively. However, the underlying molecular mechanism of NNM-mediated antiviral activity remains to be further elucidated. In this study, 414 differentially expressed genes (DEGs), including 383 which were up-regulated and 31 down-regulated, caused by NNM treatment in TMV-infected BY-2 protoplasts, were discovered by RNA-seq. In addition, KEGG analysis indicated significant enrichment of DEGs in the plant-pathogen interaction and MAPK signaling pathway. The up-regulated expression of crucial DEGs, including defense-responsive genes, such as the receptor-like kinase FLS2, RLK1, and the mitogen-activated protein kinase kinase kinase MAPKKK, calcium signaling genes, such as the calcium-binding protein CML19, as well as phytohormone responsive genes, such as the WRKY transcription factors WRKY40 and WRKY70, were confirmed by RT-qPCR. These findings provided valuable insights into the antiviral mechanisms of NNM, which indicated that the agent induces tobacco systemic resistance against TMV via activating multiple plant defense signaling pathways.

Comparison of Transcriptome Differences in Soybean Response to Soybean Mosaic Virus under Normal Light and in the Shade

Shading in the intercropping system is a major abiotic factor which influences soybean growth and development, while soybean mosaic virus (SMV) is a biotic factor that limits the yield and quality of soybean. However, little is known about the defense response of soybean to SMV in the shade. Thus, in the current study, both intensity and quality (red:far-red, R:FR) of the light were changed to simulate the shaded environment and comparative transcriptome analysis was performed. Morphologically, plant growth was inhibited by SMV, which decreased 35.93% of plant height and 8.97% of stem diameter in the shade. A total of 3548 and 4319 differentially expressed genes (DEGs) were identified in soybean plants infected with SMV under normal light and in the shade. Enrichment analysis showed that the plant defense-related genes were upregulated under normal light but downregulated in the shade. Pathways that were repressed include plant-pathogen interaction, secondary metabolism, sugar metabolism, and vitamin metabolism. In addition, genes associated with signaling pathways such as salicylic acid (SA), jasmonic acid (JA), and ethylene (ETH) were also downregulated in the shade. A qRT-PCR assay of 15 DEGs was performed to confirm transcriptome results. According to our knowledge, this is the first report on soybean response to dual stress factors. These results provide insights into the molecular mechanisms in which soybean plants were infected with SMV in the shade.

The significance of methionine cycle enzymes in plant virus infections

Both biotic and abiotic stresses cause changes in the activities of plant methionine cycle (MTC) enzymes. These changes contribute to the ability of the plant to manage stress. On the other hand, viruses utilize MTC enzymes to promote infection. Here, we review the growing but still limited knowledge of the interactions between plant viral proteins and MTC enzymes. Virus-induced changes in S-adenosyl methionine synthetase and S-adenosyl homocysteine hydrolase activities debilitate transcriptional and post-transcriptional RNA silencing and affect antiviral defense reactions connected to ethylene and polyamine biosynthesis pathways. Viral perturbations of host methionine homeostasis couple trans-sulfuration and gluthathione biosynthesis pathways to MTC functions. Large multiprotein complexes, which contain viral proteins and MTC enzymes, may represent metabolons assembled for specific viral functions or host defense responses. Proper understanding of the MTC-associated metabolic and regulatory interactions will reveal those with potential to create virus resistance in plants.

Transcriptional reprogramming caused by the geminivirus Tomato yellow leaf curl virus in local or systemic infections in Nicotiana benthamiana

BACKGROUND

Viruses have evolved to create a cellular environment permissive for viral replication in susceptible hosts. Possibly both enabling and resulting from these virus-triggered changes, infected hosts undergo a dramatic transcriptional reprogramming, the analysis of which can shed light on the molecular processes underlying the outcome of virus-host interactions. The study of the transcriptional changes triggered by the plant DNA viruses geminiviruses is potentially hampered by the low representation of infected cells in the total population, a situation that becomes extreme in those cases, like that of Tomato yellow leaf curl virus (TYLCV), in which the virus is restricted to phloem companion cells.

RESULTS

In order to gain insight into how different the transcriptional landscapes of TYLCV-infected cells or whole tissues of TYLCV-infected plants might be, here we compare the transcriptional changes in leaf patches infected with TYLCV by agroinfiltration or in systemic leaves of TYLCV-infected plants in Nicotiana benthamiana. Our results show that, in agreement with previous works, infection by TYLCV induces a dramatic transcriptional reprogramming; the detected changes, however, are not equivalent in local and systemic infections, with a much larger number of genes differentially expressed locally, and some genes responding in an opposite manner. Interestingly, a transcriptional repression of the auxin signalling pathway and a transcriptional activation of the ethylene signalling pathway were detected in both local and systemically infected samples. A transcriptional activation of defence was also detectable in both cases. Comparison with the transcriptional changes induced by systemic infection by the geminivirus Tobacco curly shoot virus (TbSV) shows common subsets of up- and down-regulated genes similarly affected by both viral species, unveiling a common transcriptional repression of terpenoid biosynthesis, a process also suppressed by the geminivirus Tomato yellow leaf curl China virus.

CONCLUSIONS

Taken together, the results presented here not only offer insight into the transcriptional changes derived from the infection by TYLCV in N. benthamiana, but also demonstrate that the resolution provided by local and systemic infection approaches largely differs, highlighting the urge to come up with a better system to gain an accurate view of the molecular and physiological changes caused by the viral invasion.

Identification and Analysis of Potential Genes Regulated by an Alphasatellite (TYLCCNA) that Contribute to Host Resistance against Tomato Yellow Leaf Curl China Virus and Its Betasatellite (TYLCCNV/TYLCCNB) Infection in Nicotiana benthamiana

Recently, begomovirus/betasatellite disease complexes were found to be associated with alphasatellites, and their presence modulated disease symptoms and/or viral DNA accumulation in infected plants. However, the biological functions of alphasatellites during begomovirus/betasatellite infections remain unclear. Tomato yellow leaf curl China virus (TYLCCNV) associated with a betasatellite (TYLCCNB) is a widespread monopartite begomovirus in China. In the Yunnan province of China, the TYLCCNV/TYLCCNB disease complex is found in association with an alphasatellite (TYLCCNA). In this study, in order to explain the mechanisms underlying TYLCCNV/TYLCCNB infection and reductions in viral DNA accumulation caused by TYLCCNA, we analyzed the transcriptome profiles of Nicotiana benthamiana seedlings challenged by TYLCCNV/TYLCCNB or TYLCCNV/TYLCCNB/TYLCCNA using RNA sequencing. In total, 2272 and 1207 differentially expressed genes (DEGs) were identified to respond to TYLCCNV/TYLCCNB and TYLCCNV/TYLCCNB/TYLCCNA infections, respectively. Compared with the DEGs in the TYLCCNV/TYLCCNB-infected N. benthamiana seedlings, the number of DEGs in plants co-infected with TYLCCNV/TYLCCNB + TYLCCNA was significantly reduced. Additionally, 36 DEGs were identified to be regulated by TYLCCNA, six of which were further analyzed using the virus-induced gene silencing (VIGS) approach. Silencing of these six TYLCCNA responsive DEGs caused more severe disease symptoms and higher viral DNA accumulation levels, suggesting that TYLCCNA responsive DEGs may attenuate TYLCCNV/TYLCCNB infection.

Transcriptomic analysis of cultivated cotton Gossypium hirsutum provides insights into host responses upon whitefly-mediated transmission of cotton leaf curl disease

Cotton is a commercial and economically important crop that generates billions of dollars in annual revenue worldwide. However, cotton yield is affected by a sap-sucking insect Bemisia tabaci (whitefly), and whitefly-borne cotton leaf curl disease (CLCuD). The causative agent of devastating CLCuD is led by the viruses belonging to the genus Begomovirus (family Geminiviridae), collectively called cotton leaf curl viruses. Unfortunately, the extensively cultivated cotton (Gossypium hirsutum) species are highly susceptible and vulnerable to CLCuD. Yet, the concomitant influence of whitefly and CLCuD on the susceptible G. hirsutum transcriptome has not been interpreted. In the present study we have employed an RNA Sequencing (RNA-Seq) transcriptomics approach to explore the differential gene expression in susceptible G. hirsutum variety upon infection with viruliferous whiteflies. Comparative RNA-Seq of control and CLCuD infected plants was done using Illumina HiSeq 2500. This study yielded 468 differentially expressed genes (DEGs). Among them, we identified 220 up and 248 downregulated DEGs involved in disease responses and pathogen defense. We selected ten genes for downstream RT-qPCR analyses on two cultivars, Karishma and MNH 786 that are susceptible to CLCuD. We observed a similar expression pattern of these genes in both susceptible cultivars that was also consistent with our transcriptome data further implying a wider application of our global transcription study on host susceptibility to CLCuD. We next performed weighted gene co-expression network analysis that revealed six modules. This analysis also identified highly co-expressed genes as well as 55 hub genes that co-express with ≥ 50 genes. Intriguingly, most of these hub genes are shown to be downregulated and enriched in cellular processes. Under-expression of such highly co-expressed genes suggests their roles in favoring the virus and enhancing plant susceptibility to CLCuD. We also discuss the potential mechanisms governing the establishment of disease susceptibility. Overall, our study provides a comprehensive differential gene expression analysis of G. hirsutum under whitefly-mediated CLCuD infection. This vital study will advance the understanding of simultaneous effect of whitefly and virus on their host and aid in identifying important G. hirsutum genes which intricate in its susceptibility to CLCuD.