Assessing Global Transcriptome Changes in Response to South African Cassava Mosaic Virus [ZA-99] Infection in Susceptible Arabidopsis thaliana

Comparative transcriptomics analysis reveals defense mechanisms of Manihot esculenta Crantz against Sri Lanka Cassava MosaicVirus

BACKGROUND

Cassava mosaic disease (CMD), caused by Sri Lankan cassava mosaic virus (SLCMV) infection, has been identified as a major pernicious disease in Manihot esculenta Crantz (cassava) plantations. It is widespread in Southeast Asia, especially in Thailand, which is one of the main cassava supplier countries. With the aim of restricting the spread of SLCMV, we explored the gene expression of a tolerant cassava cultivar vs. a susceptible cassava cultivar from the perspective of transcriptional regulation and the mechanisms underlying plant immunity and adaptation.

RESULTS

Transcriptomic analysis of SLCMV-infected tolerant (Kasetsart 50 [KU 50]) and susceptible (Rayong 11 [R 11]) cultivars at three infection stages-that is, at 21 days post-inoculation (dpi) (early/asymptomatic), 32 dpi (middle/recovery), and 67 dpi (late infection/late recovery)-identified 55,699 expressed genes. Differentially expressed genes (DEGs) between SLCMV-infected KU 50 and R 11 cultivars at (i) 21 dpi to 32 dpi (the early to middle stage), and (ii) 32 dpi to 67 dpi (the middle stage to late stage) were then identified and validated by real-time quantitative PCR (RT-qPCR). DEGs among different infection stages represent genes that respond to and regulate the viral infection during specific stages. The transcriptomic comparison between the tolerant and susceptible cultivars highlighted the role of gene expression regulation in tolerant and susceptible phenotypes.

CONCLUSIONS

This study identified genes involved in epigenetic modification, transcription and transcription factor activities, plant defense and oxidative stress response, gene expression, hormone- and metabolite-related pathways, and translation and translational initiation activities, particularly in KU 50 which represented the tolerant cultivar in this study.

The selection and validation of reference genes for quantitative real-time PCR studies in near-isogenic susceptible and resistant tomato lines, infected with the geminivirus tomato curly stunt virus

Quantitative real-time PCR (qPCR) is a sensitive and commonly used technique for gene expression profiling and provides insight into biological systems. Successful qPCR requires the use of appropriate reference genes for the normalization of data. In the present study, we aimed to identify and assess the best-suited reference genes in near-isogenic resistant (R) and susceptible (S) tomato lines infected with begomovirus Tomato curly stunt virus (ToCSV). Ten candidate reference genes namely, Actin7 (ACT), β-6 Tubulin (TUB), Ubiquitin 3 (UBI), Clathrin adaptor complexes medium subunit (CAC), Phytoene desaturase (PDS), Expressed protein (EXP), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Adenine phosphoribosyl transferase-like protein (APT1), TAP42-interacting protein (TIP41) and Elongation factor 1-alpha (EF1α) were selected and evaluated for their expression stability in resistant and susceptible tomato leaves using the analytical tools geNorm, NormFinder, BestKeeper, and RefFinder. After ranking the reference genes from most to least stable, the results suggested that a combination of ACT, EXP, and EF1α in the S lines and a combination of TIP41, APT1, and ACT in the R line is appropriate for qPCR normalization. Furthermore, to validate the identified reference genes, iron superoxide dismutase (SOD), heat shock protein 70 (HSP70) and Glutathione-S-transferase (GST) were selected as targets for normalization. The relative expression of the target genes varied when normalized against the most stable reference genes in comparison to the least stable genes. These results highlight the importance of careful selection of reference genes for accurate normalization in qPCR studies.

In-silico evolutionary analysis of plant-OBERON proteins during compatible MYMV infection in respect of improving host resistance

Yellow mosaic disease (YMD) of pulses caused by mungbean yellow mosaic virus is a major threat to crop production. An infection that is compatible with regulating and interacting host proteins and the virus causes YMD. Oberon families of proteins OBE1-4 and VIN1-4 are imperative for plants, functions in meristem and vascular development, and were also regulated during compatible disease infection. Furthermore, in-silico expression results suggested the involvement of OBE1 and OBE2 proteins during virus infection of Vigna, Arabidopsis and soybean. Moreover, a common ancestor for the meristem and virus movement related Oberons was inferred through phylogenetic analysis. Protein interaction studies showed three amino acids (Aspartate, glutamate and lysine) in the plant homeodomain (PHD), involved in interaction with the N-terminal region of the virus movement protein and were also conserved in both monocot and dicots. Additionally, major differences in the nuclear localization signal (NLS) showing clade specific conservation and significant variation between dicots and monocots were ascertained in meristem and virus movement related Oberons. Consequently, a combination of PHD, CCD and their interactions with the VPg viral domain increases the susceptibility to YMD. Further, modification in the NLS regions of the viral movement clade Oberons, to knock out allele generation in the OBE1 and OBE2 homologs through genome-editing approaches could be established as alternate strategies for the improvement of host resistance and control yellow mosaic disease in plants, especially in pulse crops.

Encounter With a Selfish Virus Sabotages Its Vector to Orient Toward Requisite Host Plant: A Case Study With Chili Leaf Curl Virus-Whitefly

Interactions of a virus with its vector and host plant have challenged entomologists, pathologists and biologists alike. Phytophagous insects depend on specific host volatile cues to locate suitable host plants for feeding and oviposition. Several studies have revealed that plant viruses modify their insect vector’s orientation toward specific host plants to facilitate their spread and survival. The ecological and molecular basis of this vector behavior modification remains largely unknown and was therefore explored in this study. Interestingly, host volatile preference for non-viruliferous female whiteflies [Bemisia tabaci (Genn.)] was found to be preferentially oriented toward infected chili plant [with chili leaf curl (ChLCV)] volatiles, while viruliferous whiteflies preferred healthy chili plant (Capsicum annum L.) volatiles in olfactometer. The electrophysiological studies involving electroantennogram (EAG) assays exhibited similar trend in EAG response amplitudes. Gas Chromatography linked electroantennodetection (GC EAD) revealed specific plant volatile cues responsible for altered host orientation behavior of the vector. Transcriptome profiling of the viruliferous and non-viruliferous whiteflies and Realtime qPCR validation showed differential expression of certain odorant binding proteins (OBPs) in viruliferous whiteflies. Our results suggest that there is a plant virus mediated altered chemoecological behavior in the vector with respect to orientation toward its host plant. Based on the findings we speculate that the virus mediates such change in the vector for a continued transmission success to the host.

Manipulation of auxin signalling by plant viruses

Compatible plant-virus interactions result in dramatic changes of the plant transcriptome and morphogenesis, and are often associated with rapid alterations in plant hormone homeostasis and signalling. Auxin controls many aspects of plant organogenesis, development, and growth; therefore, plants can rapidly perceive and respond to changes in the cellular auxin levels. Auxin signalling is a tightly controlled process and, hence, is highly vulnerable to changes in the mRNA and protein levels of its components. There are several core nuclear components of auxin signalling. In the nucleus, the interaction of auxin response factors (ARFs) and auxin/indole acetic acid (Aux/IAA) proteins is essential for the control of auxin-regulated pathways. Aux/IAA proteins are negative regulators, whereas ARFs are positive regulators of the auxin response. The interplay between both is essential for the transcriptional regulation of auxin-responsive genes, which primarily regulate developmental processes but also modulate the plant immune system. Recent studies suggest that plant viruses belonging to different families have developed various strategies to disrupt auxin signalling, namely by (a) changing the subcellular localization of Aux/IAAs, (b) preventing degradation of Aux/IAAs by stabilization, or (c) inhibiting the transcriptional activity of ARFs. These interactions perturb auxin signalling and experimental evidence from various studies highlights their importance for virus replication, systemic movement, interaction with vectors for efficient transmission, and symptom development. In this microreview, we summarize and discuss the current knowledge on the interaction of plant viruses with auxin signalling components of their hosts.

A New Type of Satellite Associated with Cassava Mosaic Begomoviruses

Cassava mosaic disease (CMD), which is caused by single-stranded DNA begomoviruses, severely limits cassava production across Africa. A previous study showed that CMD symptom severity and viral DNA accumulation increase in cassava in the presence of a DNA sequence designated SEGS-2 (sequence enhancing geminivirus symptoms). We report here that when SEGS-2 is coinoculated with African cassava mosaic virus (ACMV) onto Arabidopsis thaliana, viral symptoms increase. Transgenic Arabidopsis with an integrated copy of SEGS-2 inoculated with ACMV also display increased symptom severity and viral DNA levels. Moreover, SEGS-2 enables Cabbage leaf curl virus (CaLCuV) to infect a geminivirus-resistant Arabidopsis thaliana accession. Although SEGS-2 is related to cassava genomic sequences, an earlier study showed that it occurs as episomes and is packaged into virions in CMD-infected cassava and viruliferous whiteflies. We identified SEGS-2 episomes in SEGS-2 transgenic Arabidopsis. The episomes occur as both double-stranded and single-stranded DNA, with the single-stranded form packaged into virions. In addition, SEGS-2 episomes replicate in tobacco protoplasts in the presence, but not the absence, of ACMV DNA-A. SEGS-2 episomes contain a SEGS-2 derived promoter and an open reading frame with the potential to encode a 75-amino acid protein. An ATG mutation at the beginning of the SEGS-2 coding region does not enhance ACMV infection in A. thaliana. Together, the results established that SEGS-2 is a new type of begomovirus satellite that enhances viral disease through the action of an SEGS-2-encoded protein that may also be encoded by the cassava genome. IMPORTANCE Cassava is an important root crop in the developing world and a food and income crop for more than 300 million African farmers. Cassava is rising in global importance and trade as the demands for biofuels and commercial starch increase. More than half of the world's cassava is produced in Africa, where it is primarily grown by smallholder farmers, many of whom are from the poorest villages. Although cassava can grow under high temperature, drought, and poor soil conditions, its production is severely limited by viral diseases. Cassava mosaic disease (CMD) is one of the most important viral diseases of cassava and can cause up to 100% yield losses. We provide evidence that SEGS-2, which was originally isolated from cassava crops displaying severe and atypical CMD symptoms in Tanzanian fields, is a novel begomovirus satellite that can compromise the development of durable CMD resistance.

Study of Triticum aestivum Resistome in Response to Wheat dwarf India Virus Infection

Susceptible and resistant germplasm respond differently to pathogenic attack, including virus infections. We compared the transcriptome changes between a resistant wheat cultivar, Sonalika, and a susceptible cultivar, WL711, to understand this process in wheat against wheat dwarf India virus (WDIV) infection. A total of 2760 and 1853 genes were differentially expressed in virus-infected and mock-inoculated Sonalika, respectively, compared to WL711. The overrepresentation of genes involved in signaling, hormone metabolism, enzymes, secondary metabolites, proteolysis, and transcription factors was documented, including the overexpression of multiple PR proteins. We hypothesize that the virus resistance in Sonalika is likely due to strong intracellular surveillance via the action of multiple PR proteins (PR1, RAR1, and RPM1) and ChiB. Other genes such as PIP1, LIP1, DnaJ, defensins, oxalate oxidase, ankyrin repeat protein, serine-threonine kinase, SR proteins, beta-1,3-glucanases, and O-methyltransferases had a significant differential expression and play roles in stress tolerance, may also be contributing towards the virus resistance in Sonalika. In addition, we identified putative genes with unknown functions, which are only expressed in response to WDIV infection in Sonalika. The role of these genes could be further validated and utilized in engineering resistance in wheat and other crops.

WRKY Transcription Factors in Cassava Contribute to Regulation of Tolerance and Susceptibility to Cassava Mosaic Disease through Stress Responses

Among the numerous biological constraints that hinder cassava (Manihot esculenta Crantz) production, foremost is cassava mosaic disease (CMD) caused by virus members of the family Geminiviridae, genus Begomovirus. The mechanisms of CMD tolerance and susceptibility are not fully understood; however, CMD susceptible T200 and tolerant TME3 cassava landraces have been shown to exhibit different large-scale transcriptional reprogramming in response to South African cassava mosaic virus (SACMV). Recent identification of 85 MeWRKY transcription factors in cassava demonstrated high orthology with those in Arabidopsis, however, little is known about their roles in virus responses in this non-model crop. Significant differences in MeWRKY expression and regulatory networks between the T200 and TME3 landraces were demonstrated. Overall, WRKY expression and associated hormone and enriched biological processes in both landraces reflect oxidative and other biotic stress responses to SACMV. Notably, MeWRKY11 and MeWRKY81 were uniquely up and downregulated at 12 and 67 days post infection (dpi) respectively in TME3, implicating a role in tolerance and symptom recovery. AtWRKY28 and AtWRKY40 homologs of MeWRKY81 and MeWRKY11, respectively, have been shown to be involved in regulation of jasmonic and salicylic acid signaling in Arabidopsis. AtWRKY28 is an interactor in the RPW8-NBS resistance (R) protein network and downregulation of its homolog MeWRKY81 at 67 dpi in TME3 suggests a negative role for this WRKY in SACMV tolerance. In contrast, in T200, nine MeWRKYs were differentially expressed from early (12 dpi), middle (32 dpi) to late (67 dpi) infection. MeWRKY27 (homolog AtWRKY33) and MeWRKY55 (homolog AtWRKY53) were uniquely up-regulated at 12, 32 and 67 dpi in T200. AtWRKY33 and AtWRKY53 are positive regulators of leaf senescence and oxidative stress in Arabidopsis, suggesting MeWRKY55 and 27 contribute to susceptibility in T200.

An ERF121 transcription factor from Brassica oleracea is a target for the conserved TAL-effectors from different Xanthomonas campestris pv. campestris strains

Transcription activator-like effectors (TALEs), which induce the expression of specific plant genes to promote infection, are the main pathogenic determinants of various Xanthomonas bacteria. However, investigation of TALEs from Xanthomonas campestris pv. campestris, which causes black rot disease of crucifers, received little attention. In this study, we used PCR-based amplification followed by SMRT amplicon sequencing to identify TALE genes in several X. campestris pv. campestris strains. Computational prediction in conjunction with quantitative reverse transcription PCR analysis was used to find their targets in the Brassica oleracea genome. Transcription factor ERF121, from the AP2/ERF family, was identified as target gene for the conserved TALEs from multiple X. campestris pv. campestris strains. Several members of this family from diverse plants were previously identified as targets of TALEs from different Xanthomonas species. We propose that TALE-dependent activation of AP2/ERF transcription factors promotes susceptibility to Xanthomonas through the misregulation of plant defence pathways.

Differential miRNA profiles in South African cassava mosaic virus-infected cassava landraces reveal clues to susceptibility and tolerance to cassava mosaic disease

Specific miRNA families are involved in susceptibility or antiviral immunity in plants. Manihot esculenta Crantz (cassava) is a perennial plant that is an important food security crop in sub-Saharan Africa. Cassava is susceptible to several begomoviruses that cause cassava mosaic disease (CMD). In this study, we investigated the leaf miRNAome response in a tolerant (TME3) and susceptible (T200) cassava landrace challenged with South African cassava mosaic virus. RNAseq was performed on leaf samples at 12, 32 and 67 days post infection (dpi), representing early, symptomatic and late persistent stages of CMD infection. Significantly, distinct profiles of conserved miRNA family expression between the T200 and TME3 landraces at the three infection stages were observed. Notably at 12 days post SACMV infection, TME3 exhibited significant downregulation (log2fold<2.0) of 42 %, compared to 9% in T200, of the conserved miRNA families. This demonstrates an overall early response to SACMV in TME3 prior to symptom appearance not observed in T200, and expression of a large cohort of miRNA-regulated genes. Notably, at early infection, downregulation of mes-miR162 and 168 that target antiviral posttransriptional gene silencing (PTGS) regulators DCL1 and AGO1, respectively, was observed in TME3, and AGO1 and DCL1 expression was higher compared to T200 post infection. Early rapid responses prior to symptom development, including RNA silencing, may be key to establishing the tolerance/recovery phenotype exhibited by TME3 landrace later on at 67 dpi. At recovery, TME3 was hallmarked by a highly significant down-regulation of mes-miR167. MiR167 targets an auxin responsive factor which plays a role in auxin signaling and adaptive responses to stress, suggesting the importance of the auxin signaling in recovery of SACMV-induced symptoms. The gene targets of these miRNAs and their associated networks may provide clues to the molecular basis of CMD tolerance in perennial hosts such as cassava.

A cassava protoplast system for screening genes associated with the response to South African cassava mosaic virus

Background

The study of transient gene expression in cassava plants during virus infection using existing protocols is laborious and may take approximately fifteen weeks due to cassava’s recalcitrance to transformation. The combination of a protoplast system with CRISPR-mediated gene editing promises to shorten the turnaround time from plant tissue culture to high-throughput gene expression screening for candidate genes. Here, we detail a protocol for screening genes associated with the response to South African cassava mosaic virus (SACMV) in cassava protoplasts, with reference to the ubiquitin E3 ligase gene, MeE3L.

Methods

Cassava protoplasts of model, and SACMV-susceptible and -tolerant genotypes, were transformed with SACMV infectious clones and/or a CRISPR-editing construct targeting the MeE3L using PEG4000-mediated transfection. DNA and RNA were extracted from transformed protoplasts at 24 h post-transfection. Relative SACMV DNA accumulation was determined via qPCR using DpnI-digested total DNA, MeE3L relative expression was determined via reverse transcriptase qPCR, and results were analysed using one-way ANOVA, Tukey’s HSD test and the 2^−ΔΔCTstatistical method. The MeE3L exonic region was sequenced on the ABI 3500XL Genetic Analyzer platform; and sequences were analysed for mutations using MAFTT and MEGA-X software. Construction of a phylogenetic tree was done using the Maximum Likelihood method and Jones-Taylor-Thornton (JTT) matrix-based model.

Results

The differential expression of unedited and mutant MeE3L during SACMV infection of model, susceptible and tolerant cassava protoplasts was determined within 7 weeks after commencement of tissue culture. The study also revealed that SACMV DNA accumulation in cassava protoplasts is genotype-dependent and induces multiple mutations in the tolerant landrace MeE3L homolog. Notably, the susceptible cassava landrace encodes a RINGless MeE3Lwhich is silenced by SACMV-induced mutations. SACMV also induces mutations which silence the MeE3L RING domain in protoplasts from and tolerant cassava landraces.

Conclusions

This protocol presented here halves the turnaround time for high-throughput screening of genes associated with the host response to SACMV. It provides evidence that a cassava E3 ligase is associated with the response to SACMV and forms a basis for validation of these findings by in planta functional and interaction studies.

Characterization of Local and Systemic Impact of Whitefly (Bemisia tabaci) Feeding and Whitefly-Transmitted Tomato Mottle Virus Infection on Tomato Leaves by Comprehensive Proteomics

Tomato mottle virus (ToMoV) is a single-stranded DNA (ssDNA) begomovirus transmitted to solanaceous crops by the whitefly species complex (Bemisia tabaci), causing stunted growth, leaf mottling, and reduced yield. Using a genetic repertoire of seven genes, ToMoV pathogenesis includes the manipulation of multiple plant biological processes to circumvent antiviral defenses. To further understand the effects of whitefly feeding and whitefly-transmitted ToMoV infection on tomato plants (Solanum lycopersicum ‘Florida Lanai’), we generated comprehensive protein profiles of leaves subjected to feeding by either viruliferous whiteflies harboring ToMoV, or non-viruliferous whiteflies, or a no-feeding control. The effects of whitefly feeding and ToMoV infection were measured both locally and systemically by sampling either a mature leaf directly from the site of clip-cage confined whitefly feeding, or from a newly formed leaf 10 days post feeding (dpf). At 3 dpf, tomato’s response to ToMoV included proteins associated with translation initiation and elongation as well as plasmodesmata dynamics. In contrast, systemic impacts of ToMoV on younger leaves 10 dpf were more pronounced and included a virus-specific change in plant proteins associated with mRNA maturation and export, RNA-dependent DNA methylation, and other antiviral plant processes. Our analysis supports previous findings and provides novel insight into tomato’s local and systemic response to whitefly feeding and ToMoV infection.

A plant DNA virus replicates in the salivary glands of its insect vector via recruitment of host DNA synthesis machinery

Viruses pose a great threat to animal and plant health worldwide. Whereas most plant viruses only replicate in plant hosts, some also replicate in their animal (insect) vector. A detailed knowledge of host expansion will give a better understanding of virus evolution, and identification of virus and host components involved in this process can lead to new strategies to combat virus spread. Here, we reveal that a plant DNA virus has evolved to induce and recruit insect DNA synthesis machinery to support its replication in vector salivary glands. Our study sheds light on the understanding of TYLCV–whitefly interactions and provides insights into how a plant virus may evolve to infect and replicate in an insect vector.

Whereas most of the arthropod-borne animal viruses replicate in their vectors, this is less common for plant viruses. So far, only some plant RNA viruses have been demonstrated to replicate in insect vectors and plant hosts. How plant viruses evolved to replicate in the animal kingdom remains largely unknown. Geminiviruses comprise a large family of plant-infecting, single-stranded DNA viruses that cause serious crop losses worldwide. Here, we report evidence and insight into the replication of the geminivirus tomato yellow leaf curl virus (TYLCV) in the whitefly (Bemisia tabaci) vector and that replication is mainly in the salivary glands. We found that TYLCV induces DNA synthesis machinery, proliferating cell nuclear antigen (PCNA) and DNA polymerase δ (Polδ), to establish a replication-competent environment in whiteflies. TYLCV replication-associated protein (Rep) interacts with whitefly PCNA, which recruits DNA Polδ for virus replication. In contrast, another geminivirus, papaya leaf curl China virus (PaLCuCNV), does not replicate in the whitefly vector. PaLCuCNV does not induce DNA-synthesis machinery, and the Rep does not interact with whitefly PCNA. Our findings reveal important mechanisms by which a plant DNA virus replicates across the kingdom barrier in an insect and may help to explain the global spread of this devastating pathogen.

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.

Diverse plant viruses: a toolbox for dissection of cellular pathways

Research in virology has usually focused on one selected host-virus pathosystem to examine the mechanisms underlying a particular disease. However, as exemplified by the mechanistically versatile suppression of antiviral RNA silencing by plant viruses, there may be functionally convergent evolution. Assuming this is a widespread feature, we propose that effector proteins from diverse plant viruses can be a powerful resource for discovering new regulatory mechanisms of distinct cellular pathways. The efficiency of this approach will depend on how deeply and widely the studied pathway is integrated into viral infections. Beyond this, comparative studies using broad virus diversity should increase our global understanding of plant-virus interactions.

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.

Transcriptome analysis of Nicotiana benthamiana infected by Tobacco curly shoot virus

BACKGROUND

Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus associated with betasatellite (Tobacco curly shoot betasatellite, TbCSB), which causes serious leaf curl disease on tomato and tobacco in China. It is interesting that TbCSV induced severe upward leaf curling in Nicotiana benthamiana, but in the presence of TbCSB, symptoms changed to be downward leaf curling. However, the mechanism of interactions between viral pathogenicity, host defense, viral-betasatellite interactions and virus-host interactions remains unclear.

METHODS

In this study, RNA-seq was used to analyze differentially expressed genes (DEGs) in N. benthamiana plants infected by TbCSV (Y35A) and TbCSV together with TbCSB (Y35AB) respectively.

RESULTS

Through mapping to N. benthamiana reference genome, 59,814 unigenes were identified. Transcriptome analysis revealed that a total of 4081 and 3196 DEGs were identified in Y35AB vs CK (control check) and Y35A vs CK, respectively. Both GO and KEGG analyses were conducted to classify the DEGs. Ten of the top 15 GO terms were enriched in both DEGs of Y35AB vs CK and Y35A vs CK, and these enriched GO terms mainly classified into three categories including biological process, cellular component and molecular function. KEGG pathway analysis indicated that 118 and 111 pathways were identified in Y35AB vs CK and Y35A vs CK, respectively, of which nine and six pathways were significantly enriched. Three major pathways in Y35AB vs CK involved in metabolic pathways, carbon metabolism and photosynthesis, while those in Y35A vs CK were related to Ribosome, Glyoxylate and dicarboxylate metabolism and DNA replication. We observed that 8 PR genes were significantly up-regulated and 44 LRR-RLK genes were significantly differentially expressed in Y35A treatment or in Y35AB treatment. In addition, 7 and 13 genes were identified to be significantly changed in biosynthesis and signal transduction pathway of brassinosteroid (BR) and jasmonic acid (JA) respectively.

CONCLUSIONS

These results presented here would be particularly useful to further elucidate the response of the host plant against virus infection.

Plum pox virus capsid protein suppresses plant pathogen-associated molecular pattern (PAMP)-triggered immunity

The perception of pathogen-associated molecular patterns (PAMPs) by immune receptors launches defence mechanisms referred to as PAMP-triggered immunity (PTI). Successful pathogens must suppress PTI pathways via the action of effectors to efficiently colonize their hosts. So far, plant PTI has been reported to be active against most classes of pathogens, except viruses, although this defence layer has been hypothesized recently as an active part of antiviral immunity which needs to be suppressed by viruses for infection success. Here, we report that Arabidopsis PTI genes are regulated upon infection by viruses and contribute to plant resistance to Plum pox virus (PPV). Our experiments further show that PPV suppresses two early PTI responses, the oxidative burst and marker gene expression, during Arabidopsis infection. In planta expression of PPV capsid protein (CP) was found to strongly impair these responses in Nicotiana benthamiana and Arabidopsis, revealing its PTI suppressor activity. In summary, we provide the first clear evidence that plant viruses acquired the ability to suppress PTI mechanisms via the action of effectors, highlighting a novel strategy employed by viruses to escape plant defences.

Plant immunity against viruses: antiviral immune receptors in focus

BACKGROUND

Among the environmental limitations that affect plant growth, viruses cause major crop losses worldwide and represent serious threats to food security. Significant advances in the field of plant-virus interactions have led to an expansion of potential strategies for genetically engineered resistance in crops during recent years. Nevertheless, the evolution of viral virulence represents a constant challenge in agriculture that has led to a continuing interest in the molecular mechanisms of plant-virus interactions that affect disease or resistance.

SCOPE AND CONCLUSION

This review summarizes the molecular mechanisms of the antiviral immune system in plants and the latest breakthroughs reported in plant defence against viruses. Particular attention is given to the immune receptors and transduction pathways in antiviral innate immunity. Plants counteract viral infection with a sophisticated innate immune system that resembles the non-viral pathogenic system, which is broadly divided into pathogen-associated molecular pattern (PAMP)-triggered immunity and effector-triggered immunity. An additional recently uncovered virus-specific defence mechanism relies on host translation suppression mediated by a transmembrane immune receptor. In all cases, the recognition of the virus by the plant during infection is central for the activation of these innate defences, and, conversely, the detection of host plants enables the virus to activate virulence strategies. Plants also circumvent viral infection through RNA interference mechanisms by utilizing small RNAs, which are often suppressed by co-evolving virus suppressors. Additionally, plants defend themselves against viruses through hormone-mediated defences and activation of the ubiquitin-26S proteasome system (UPS), which alternatively impairs and facilitates viral infection. Therefore, plant defence and virulence strategies co-evolve and co-exist; hence, disease development is largely dependent on the extent and rate at which these opposing signals emerge in host and non-host interactions. A deeper understanding of plant antiviral immunity may facilitate innovative biotechnological, genetic and breeding approaches for crop protection and improvement.

iTRAQ analysis of the tobacco leaf proteome reveals that RNA-directed DNA methylation (RdDM) has important roles in defense against geminivirus-betasatellite infection

Geminiviruses have caused serious losses in crop production. To investigate the mechanisms underlying host defenses against geminiviruses, an isobaric tags for relative and absolute quantification (iTRAQ)-based quantitative proteomic approach was used to explore the expression profiles of proteins in Nicotiana benthamiana (N. benthamiana) leaves in response to tomato yellow leaf curl China virus (TYLCCNV) with its betasatellite (TYLCCNB) at an early phase. In total, 4155 proteins were identified and 272 proteins were changed differentially in response to TYLCCNV/TYLCCNB infection. Bioinformatics analysis indicated that S-adenosyl-l-methionine cycle II was the most significantly up-regulated biochemical process during TYLCCNV/TYLCCNB infection. The mRNA levels of three proteins in S-adenosyl-l-methionine cycle II were further analyzed by qPCR, each was found significantly up-regulated in TYLCCNV/TYLCCNB-infected N. benthamiana. This result suggested a strong promotion of the biosynthesis of available methyl groups during geminivirus infections. We further tested the potential role of RdDM in N. benthamiana by virus-induced gene silencing (VIGS) and found that a disruption in RdDM resulted in more severe infectious symptoms and higher accumulation of viral DNA after TYLCCNV/TYLCCNB infection. Although the precise functions of these proteins still need to be determined, our proteomic results enhance the understanding of plant antiviral mechanisms.

BIOLOGICAL SIGNIFICANCE

One of the major limitations to crop growth in the worldwide is the prevalence of geminiviruses. They are able to infect food and cash crops and cause serious crop failures and economic losses worldwide, especially in Africa and Asia. Tomato yellow leaf curl China virus (TYLCCNV), which causes severe viral diseases in China, is a monopartite geminivirus associated with the betasatellite (TYLCCNB). However, the mechanisms underlying the TYLCCNV/TYLCCNB defense in plants are still not fully understood at the molecular level. In this study, the combined proteomic, bioinformatic and VIGS analyses revealed that TYLCCNV/TYLCCNB invasion caused complex proteomic alterations in the leaves of N. benthamiana involving the processes of stress and defense, energy production, photosynthesis, protein homeostasis, metabolism, cell structure, signal transduction, transcription, transportation, and cell growth/division. Promotion of available methyl groups via the S-adenosyl-l-methionine cycle II pathway in N. benthamiana appeared crucial for antiviral responses. These findings enhance our understanding in the proteomic aspects of host antiviral defenses against geminiviruses, and also demonstrate that the combination of proteomics with bioinformatics and VIGS analysis is an effective approach to investigate systemic plant responses to geminiviruses and to shed light on plant-virus interactions.

An insight into differentially regulated genes in resistant and susceptible genotypes of potato in response to tomato leaf curl New Delhi virus-[potato] infection

Apical leaf curl disease, caused by tomato leaf curl New Delhi virus-[potato] (ToLCNDV-[potato]) is one of the most important viral diseases of potato in India. Genetic resistance source for ToLCNDV in potato is not identified so far. However, the cultivar Kufri Bahar is known to show lowest seed degeneration even under high vector levels. Hence, microarray analysis was performed to identify differentially regulated genes during ToLCNDV-[potato] infection in a resistant (Kufri Bahar) and a susceptible cultivar (Kufri Pukhraj). Under artificial inoculation conditions, in Kufri Pukhraj, symptom expressions started at 15days after inoculation (DAI) and then progressed to severe symptoms, whereas no or only very mild symptoms were observed in Kufri Bahar up to 35 DAI. Correspondingly, qPCR assay indicated a high viral load in Kufri Pukhraj and a very low viral load in Kufri Bahar. Microarray analysis showed that a total of 1111 genes and 2588 genes were differentially regulated (|log2 (Fold Change)|>2) in Kufri Bahar and Kufri Pukhraj, respectively, following ToLCNDV-[potato] infection. Gene ontology and mapman analyses revealed that these altered transcripts were involved in various biological & metabolic processes. Several genes with unknown functions were 5 to 100 fold expressed after virus infection and further experiments are necessary to ascertain their role in disease resistance or susceptibility. This study gives an insight into differentially regulated genes in response to ToLCNDV-[potato] infection in resistant and susceptible cultivars and could serve as the basis for the development of new strategies for disease management.

RNA-seq-based digital gene expression analysis reveals modification of host defense responses by rice stripe virus during disease symptom development in Arabidopsis

Background

Virus infection induces and suppresses host gene expression on a global level. Rice stripe virus (RSV) is the type species of the genus Tenuivirus and infects rice and Arabidopsis plants. Microarray-based and next generation sequencing-based transcriptomic approaches have been used to study rice-RSV interactions. However, our knowledge of the response of Arabidopsis plants to RSV infection is limited, and it requires further investigation to determine the similarities (or differences) in virus-host interactions between monocot and dicot hosts infected with RSV.

Methods

We characterized transcriptome changes in Arabidopsis thaliana infected with rice stripe virus (RSV) with RNA-seq based digital gene expression (DGE) analysis. The transcriptomes of RSV-infected samples were compared to those of mock-treated samples at 14 and 21 days post-infection (dpi) during different stages of symptom development.

Results

We identified 624 differentially expressed genes (DEGs) in Arabidopsis influenced by RSV at 14 dpi and 21 dpi, among which at 14 dpi, 255 transcripts were induced, and 38 were repressed; at 21 dpi, 146 were induced, and 237 were repressed. Functional annotation indicated that these DEGs were related to multiple biological functions, including defense response, secondary metabolism, protein amino acid phosphorylation and response to abiotic stress.

Conclusions

Importantly, the transcription of genes related to host defense systems was activated by RSV infection at an early stage of symptom development (14 dpi), whereas over the infection period (21 dpi), the host defense response systems were suppressed. A total of 52 genes were continuously differentially expressed between the two time points, indicating that the majority of DEGs were transient and unique to a particular time point during symptom development. The DEGs, particularly the defense response genes, identified in this study are candidates suitable for further functional analysis during the RSV-Arabidopsis interaction.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0663-7) contains supplementary material, which is available to authorized users.

Transcriptomics of tomato plants infected with TYLCSV or expressing the central TYLCSV Rep protein domain uncover changes impacting pathogen response and senescence

To establish a successful infection viruses need to overcome plant innate immune responses and redirect host gene expression for their multiplication and diffusion. Tomato yellow leaf curl Sardinia virus (TYLCSV) is a geminivirus, which causes significant economic losses in tomato. The multifunctional replication associated geminivirus protein (Rep) has an important role during viral infection. In particular, the Rep central domain spanning from aa 120 to 180 is known to interact with viral and host factors. In this study, we used long serial analysis of gene expression to analyse the transcriptional profiles of transgenic tomato plants expressing the first 210 amino acids of TYLCSV Rep (Rep210) and TYLCSV-infected wild-type tomato plants (Wt-Ty). Also, we compared these profiles with those of transgenic Rep130 tomatoes. Comparison of Wt-Ty and Rep210 libraries with the wild-type one identified 118 and 203 differentially expressed genes (DEGs), respectively. Importantly, 55% of Wt-Ty DEGs were in common with Rep210, and no ones showed opposite expression. Conversely, a negligible overlap was found between Rep130 DEGs and Wt-Ty and Rep210 ones. TYLCSV- and Rep210-repressed genes, but not induced ones, overlapped with the leaf senescence process. Interestingly, TYLCSV upregulates expression of genes involved in the negative regulation of programmed cell death (PCD), several of which were also regulated by the abscisic acid. Rep210 upregulated genes related to defence response, immune system processes and negative regulation of PCD. Collectively, our results support a model in which the Rep central domain has a pivotal role in redirecting host plant gene expression.

Two Novel DNAs That Enhance Symptoms and Overcome CMD2 Resistance to Cassava Mosaic Disease

Cassava mosaic begomoviruses (CMBs) cause cassava mosaic disease (CMD) across Africa and the Indian subcontinent. Like all members of the geminivirus family, CMBs have small, circular single-stranded DNA genomes. We report here the discovery of two novel DNA sequences, designated SEGS-1 and SEGS-2 (forsequencesenhancinggeminivirussymptoms), that enhance symptoms and break resistance to CMD. The SEGS are characterized by GC-rich regions and the absence of long open reading frames. Both SEGS enhanced CMD symptoms in cassava (Manihot esculentaCrantz) when coinoculated withAfrican cassava mosaic virus(ACMV),East African cassava mosaic Cameroon virus(EACMCV), orEast African cassava mosaic virus-Uganda(EACMV-UG). SEGS-1 also overcame resistance of a cassava landrace carrying the CMD2 resistance locus when coinoculated with EACMV-UG. Episomal forms of both SEGS were detected in CMB-infected cassava but not in healthy cassava. SEGS-2 episomes were also found in virions and whiteflies. SEGS-1 has no homology to geminiviruses or their associated satellites, but the cassava genome contains a sequence that is 99% identical to full-length SEGS-1. The cassava genome also includes three sequences with 84 to 89% identity to SEGS-2 that together encompass all of SEGS-2 except for a 52-bp region, which includes the episomal junction and a 26-bp sequence related to alphasatellite replication origins. These results suggest that SEGS-1 is derived from the cassava genome and facilitates CMB infection as an integrated copy and/or an episome, while SEGS-2 was originally from the cassava genome but now is encapsidated into virions and transmitted as an episome by whiteflies.

IMPORTANCE

Cassava is a major crop in the developing world, with its production in Africa being second only to maize. CMD is one of the most important diseases of cassava and a serious constraint to production across Africa. CMD2 is a major CMD resistance locus that has been deployed in many cassava cultivars through large-scale breeding programs. In recent years, severe, atypical CMD symptoms have been observed occasionally on resistant cultivars, some of which carry the CMD2 locus, in African fields. In this report, we identified and characterized two DNA sequences, SEGS-1 and SEGS-2, which produce similar symptoms when coinoculated with cassava mosaic begomoviruses onto a susceptible cultivar or a CMD2-resistant landrace. The ability of SEGS-1 to overcome CMD2 resistance and the transmission of SEGS-2 by whiteflies has major implications for the long-term durability of CMD2 resistance and underscore the need for alternative sources of resistance in cassava.

The Pokeweed Leaf mRNA Transcriptome and Its Regulation by Jasmonic Acid

The American pokeweed plant, Phytolacca americana, is recognized for synthesizing pokeweed antiviral protein (PAP), a ribosome inactivating protein (RIP) that inhibits the replication of several plant and animal viruses. The plant is also a heavy metal accumulator with applications in soil remediation. However, little is known about pokeweed stress responses, as large-scale sequencing projects have not been performed for this species. Here, we sequenced the mRNA transcriptome of pokeweed in the presence and absence of jasmonic acid (JA), a hormone mediating plant defense. Trinity-based de novo assembly of mRNA from leaf tissue and BLASTx homology searches against public sequence databases resulted in the annotation of 59 096 transcripts. Differential expression analysis identified JA-responsive genes that may be involved in defense against pathogen infection and herbivory. We confirmed the existence of several PAP isoforms and cloned a potentially novel isoform of PAP. Expression analysis indicated that PAP isoforms are differentially responsive to JA, perhaps indicating specialized roles within the plant. Finally, we identified 52 305 natural antisense transcript pairs, four of which comprised PAP isoforms, suggesting a novel form of RIP gene regulation. This transcriptome-wide study of a Phytolaccaceae family member provides a source of new genes that may be involved in stress tolerance in this plant. The sequences generated in our study have been deposited in the SRA database under project # SRP069141.

Involvement of host regulatory pathways during geminivirus infection: a novel platform for generating durable resistance

Geminiviruses are widely distributed throughout the world and cause devastating yield losses in almost all the economically important crops. In this review, the newly identified roles of various novel plant factors and pathways participating in plant–virus interaction are summarized with a particular focus on the exploitation of various pathways involving ubiquitin/26S proteasome pathway, small RNA pathways, cell division cycle components, and the epigenetic mechanism as defense responses during plant–pathogen interactions. Capturing the information on these pathways for the development of strategies against geminivirus infection is argued to provide the basis for new genetic approaches to resistance.

Mutualism-parasitism paradigm synthesized from results of root-endophyte models

Plant tissues host a variety of fungi. One important group is the dark septate endophytes (DSEs) that colonize plant roots and form characteristic intracellular structures - melanized hyphae and microsclerotia. The DSE associations are common and frequently observed in various biomes and plant taxa. Reviews suggest that the proportion of plant species colonized by DSE equal that colonized by AM and microscopic studies show that the proportion of the root system colonized by fungi DSE can equal, or even exceed, the colonization by AM fungi. Despite the high frequency and suspected ecological importance, the effects of DSE colonization on plant growth and performance have remained unclear. Here, we draw from over a decade of experimentation with the obscure DSE symbiosis and synthesize across large bodies of published and unpublished data from Arabidopsis thaliana and Allium porrum model systems as well as from experiments that use native plants to better resolve the host responses to DSE colonization. The data indicate similar distribution of host responses in model and native plant studies, validating the use of model plants for tractable dissection of DSE symbioses. The available data also permit empirical testing of the environmental modulation of host responses to DSE colonization and refining the "mutualism-parasitism-continuum" paradigm for DSE symbioses. These data highlight the context dependency of the DSE symbioses: not only plant species but also ecotypes vary in their responses to populations of conspecific DSE fungi - environmental conditions further shift the host responses similar to those predicted based on the mutualism-parasitism-continuum paradigm. The model systems provide several established avenues of inquiry that permit more detailed molecular and functional dissection of fungal endophyte symbioses, identifying thus likely mechanisms that may underlie the observed host responses to endophyte colonization.

Transcriptional analysis of South African cassava mosaic virus-infected susceptible and tolerant landraces of cassava highlights differences in resistance, basal defense and cell wall associated genes during infection

BACKGROUND

Cassava mosaic disease is caused by several distinct geminivirus species, including South African cassava mosaic virus-[South Africa:99] (SACMV). To date, there is limited gene regulation information on viral stress responses in cassava, and global transcriptome profiling in SACMV-infected cassava represents an important step towards understanding natural host responses to plant geminiviruses.

RESULTS

A RNA-seq time course (12, 32 and 67 dpi) study, monitoring gene expression in SACMV-challenged susceptible (T200) and tolerant (TME3) cassava landraces, was performed using the Applied Biosystems (ABI) SOLiD next-generation sequencing platform. The multiplexed paired end sequencing run produced a total of 523 MB and 693 MB of paired-end reads for SACMV-infected susceptible and tolerant cDNA libraries, respectively. Of these, approximately 50.7% of the T200 reads and 55.06% of TME3 reads mapped to the cassava reference genome available in phytozome. Using a log2 fold cut-off (p<0.05), comparative analysis between the six normalized cDNA libraries showed that 4181 and 1008 transcripts in total were differentially expressed in T200 and TME3, respectively, across 12, 32 and 67 days post infection, compared to mock-inoculated. The number of responsive transcripts increased dramatically from 12 to 32 dpi in both cultivars, but in contrast, in T200 the levels did not change significantly at 67 dpi, while in TME3 they declined. GOslim functional groups illustrated that differentially expressed genes in T200 and TME3 were overrepresented in the cellular component category for stress-related genes, plasma membrane and nucleus. Alterations in the expression of other interesting genes such as transcription factors, resistance (R) genes, and histone/DNA methylation-associated genes, were observed. KEGG pathway analysis uncovered important altered metabolic pathways, including phenylpropanoid biosynthesis, sucrose and starch metabolism, and plant hormone signalling.

CONCLUSIONS

Molecular mechanisms for TME3 tolerance are proposed, and differences in patterns and levels of transcriptome profiling between T200 and TME3 with susceptible and tolerant phenotypes, respectively, support the hypothesis that viruses rearrange their molecular interactions in adapting to hosts with different genetic backgrounds.

SnRK1 phosphorylation of AL2 delays Cabbage leaf curl virus infection in Arabidopsis

Geminivirus AL2/C2 proteins play key roles in establishing infection and causing disease in their plant hosts. They are involved in viral gene expression, counter host defenses by suppressing transcriptional gene silencing, and interfere with the host signaling involved in pathogen resistance. We report here that begomovirus and curtovirus AL2/C2 proteins interact strongly with host geminivirus Rep-interacting kinases (GRIKs), which are upstream activating kinases of the protein kinase SnRK1, a global regulator of energy and nutrient levels in plants. We used an in vitro kinase system to show that GRIK-activated SnRK1 phosphorylates recombinant AL2/C2 proteins from several begomoviruses and to map the SnRK1 phosphorylation site to serine-109 in the AL2 proteins of two New World begomoviruses: Cabbage Leaf Curl Virus (CaLCuV) and Tomato mottle virus. A CaLCuV AL2 S109D phosphomimic mutation did not alter viral DNA levels in protoplast replication assays. In contrast, the phosphomimic mutant was delayed for symptom development and viral DNA accumulation during infection of Arabidopsis thaliana, demonstrating that SnRK1 contributes to host defenses against CaLCuV. Our observation that serine-109 is not conserved in all AL2/C2 proteins that are SnRK1 substrates in vitro suggested that phosphorylation of viral proteins by plant kinases contributes to the evolution of geminivirus-host interactions.

IMPORTANCE

Geminiviruses are single-stranded DNA viruses that cause serious diseases in many crops. Dicot-infecting geminiviruses carry genes that encode multifunctional AL2/C2 proteins that are essential for infection. However, it is not clear how AL2/C2 proteins are regulated. Here, we show that the host protein kinase SnRK1, a central regulator of energy balance and nutrient metabolism in plants, phosphorylates serine-109 in AL2 proteins of three subgroups of New World begomoviruses, resulting in a delay in viral DNA accumulation and symptom appearance. Our results support SnRK1's antiviral role and reveal a novel mechanism underlying this function. Phylogenetic analysis suggested that AL2 S109 evolved as begomoviruses migrated from the Old World to the New World and may have provided a selective advantage as begomoviruses adapted to a different environment and different plant hosts. This study provides new insights into the interaction of viral pathogens with their plant hosts at the level of viral protein modification by the host.

Differential soybean gene expression during early phase of infection with Mungbean yellow mosaic India virus

Mungbean yellow mosaic India virus (MYMIV), a bipartite begomovirus, causes yellow mosaic disease to soybean. Studies related to host gene expression in response to begomovirus infection have mostly been performed with systemically infected tissues at a later period of infection. In this study, soybean gene expression analysis has been performed to understand local responses against MYMIV at an early stage of infection before appearance of detectable limit of late viral transcripts. 444 soybean transcripts belonging to eleven functional categories showed significant changes in expression level at two days after infection. MYMIV infection resulted in enhanced expression of genes associated with hypersensitive response, programmed cell death and resistance response pathways and reduced expression of genes for photosynthesis and sugar transport. Comparative expression analysis of selected transcripts in the susceptible and a resistant variety displayed differential expression of host genes involved in intercellular virus movement and long distance signaling of systemic acquired resistance.

SnRK1 phosphorylation of AL2 delays Cabbage leaf curl virus infection in Arabidopsis

Geminivirus AL2/C2 proteins play key roles in establishing infection and causing disease in their plant hosts. They are involved in viral gene expression, counter host defenses by suppressing transcriptional gene silencing, and interfere with the host signaling involved in pathogen resistance. We report here that begomovirus and curtovirus AL2/C2 proteins interact strongly with host geminivirus Rep-interacting kinases (GRIKs), which are upstream activating kinases of the protein kinase SnRK1, a global regulator of energy and nutrient levels in plants. We used an in vitro kinase system to show that GRIK-activated SnRK1 phosphorylates recombinant AL2/C2 proteins from several begomoviruses and to map the SnRK1 phosphorylation site to serine-109 in the AL2 proteins of two New World begomoviruses: Cabbage Leaf Curl Virus (CaLCuV) and Tomato mottle virus. A CaLCuV AL2 S109D phosphomimic mutation did not alter viral DNA levels in protoplast replication assays. In contrast, the phosphomimic mutant was delayed for symptom development and viral DNA accumulation during infection of Arabidopsis thaliana, demonstrating that SnRK1 contributes to host defenses against CaLCuV. Our observation that serine-109 is not conserved in all AL2/C2 proteins that are SnRK1 substrates in vitro suggested that phosphorylation of viral proteins by plant kinases contributes to the evolution of geminivirus-host interactions.

IMPORTANCE

Geminiviruses are single-stranded DNA viruses that cause serious diseases in many crops. Dicot-infecting geminiviruses carry genes that encode multifunctional AL2/C2 proteins that are essential for infection. However, it is not clear how AL2/C2 proteins are regulated. Here, we show that the host protein kinase SnRK1, a central regulator of energy balance and nutrient metabolism in plants, phosphorylates serine-109 in AL2 proteins of three subgroups of New World begomoviruses, resulting in a delay in viral DNA accumulation and symptom appearance. Our results support SnRK1's antiviral role and reveal a novel mechanism underlying this function. Phylogenetic analysis suggested that AL2 S109 evolved as begomoviruses migrated from the Old World to the New World and may have provided a selective advantage as begomoviruses adapted to a different environment and different plant hosts. This study provides new insights into the interaction of viral pathogens with their plant hosts at the level of viral protein modification by the host.

Transcriptomics of the interaction between the monopartite phloem-limited geminivirus tomato yellow leaf curl Sardinia virus and Solanum lycopersicum highlights a role for plant hormones, autophagy and plant immune system fine tuning during infection

Tomato yellow leaf curl Sardinia virus (TYLCSV), a DNA virus belonging to the genus Begomovirus, causes severe losses in tomato crops. It infects only a limited number of cells in the vascular tissues, making difficult to detect changes in host gene expression linked to its presence. Here we present the first microarray study of transcriptional changes induced by the phloem-limited geminivirus TYLCSV infecting tomato, its natural host. The analysis was performed on the midrib of mature leaves, a material naturally enriched in vascular tissues. A total of 2206 genes were up-regulated and 1398 were down-regulated in infected plants, with an overrepresentation of genes involved in hormone metabolism and responses, nucleic acid metabolism, regulation of transcription, ubiquitin-proteasome pathway and autophagy among those up-regulated, and in primary and secondary metabolism, phosphorylation, transcription and methylation-dependent chromatin silencing among those down-regulated. Our analysis showed a series of responses, such as the induction of GA- and ABA-responsive genes, the activation of the autophagic process and the fine tuning of the plant immune system, observed only in TYLCSV-tomato compatible interaction so far. On the other hand, comparisons with transcriptional changes observed in other geminivirus-plant interactions highlighted common host responses consisting in the deregulation of biotic stress responsive genes, key enzymes in the ethylene biosynthesis and methylation cycle, components of the ubiquitin proteasome system and DNA polymerases II. The involvement of conserved miRNAs and of solanaceous- and tomato-specific miRNAs in geminivirus infection, investigated by integrating differential gene expression data with miRNA targeting data, is discussed.

Arabidopsis thaliana NAC083 protein interacts with Mungbean yellow mosaic India virus (MYMIV) Rep protein

Geminiviral replication initiator protein (Rep) is a key player in geminiviral rolling circle mode of replication. However, the virus exploits various host cellular machineries for its replication. Study of these host factors is important to understand the geminiviral DNA replication in greater details. With this view, we screened for the peptides interacting with the Rep protein of a representative of geminivirus, namely, Mungbean yellow mosaic India virus (MYMIV), employing phage display technique. Through this screen, we have identified a host transcription factor, NAC083, as a potential MYMIV-Rep-binding partner. In silico docking studies also suggested possible binding of NAC083 peptide to MYMIV-Rep. We validated the interaction between MYMIV-Rep and Arabidopsis thaliana full-length NAC083 protein using in vitro pull-down assay and yeast two-hybrid analysis. NAC proteins are well-known transcription factors belonging to the largest gene families in plants. This study demonstrates for the first time the interaction of NAC083, a member of NAC transcription factor family, with MYMIV-Rep protein thereby indicating its possible role in MYMIV DNA replication.

Relationship between gene responses and symptoms induced by Rice grassy stunt virus

Rice grassy stunt virus (RGSV) is a serious threat to rice production in Southeast Asia. RGSV is a member of the genus Tenuivirus, and it induces leaf yellowing, stunting, and excess tillering on rice plants. Here we examined gene responses of rice to RGSV infection to gain insight into the gene responses which might be associated with the disease symptoms. The results indicated that (1) many genes related to cell wall synthesis and chlorophyll synthesis were predominantly suppressed by RGSV infection; (2) RGSV infection induced genes associated with tillering process; (3) RGSV activated genes involved in inactivation of gibberellic acid and indole-3-acetic acid; and (4) the genes for strigolactone signaling were suppressed by RGSV. These results suggest that these gene responses to RGSV infection account for the excess tillering specific to RGSV infection as well as other symptoms by RGSV, such as stunting and leaf chlorosis.