Discovering host genes involved in the infection by the Tomato Yellow Leaf Curl Virus complex and in the establishment of resistance to the virus using Tobacco Rattle Virus-based post transcriptional gene silencing

Insights into the early transcriptomic response against watermelon mosaic virus in melon

BACKGROUND

Watermelon mosaic virus (WMV) is one of the most prevalent viruses affecting melon worldwide. Recessive resistance to WMV in melon has previously been reported in the African accession TGR-1551. Moreover, the genomic regions associated to the resistance have also been described. Nevertheless, the transcriptomic response that might infer the resistance to this potyvirus has not been explored.

RESULTS

We have performed a comparative transcriptomic analysis using mock and WMV-inoculated plants of the susceptible cultivar "Bola de oro" (BO) and a resistant RIL (Recombinant inbred line) derived from the initial cross between "TGR-1551" and BO. In total, 616 genes were identified as differentially expressed and the weighted gene co-expression network analysis (WGCNA) detected 19 gene clusters (GCs), of which 7 were differentially expressed for the genotype x treatment interaction term. SNPs with a predicted high impact on the protein function were detected within the coding regions of most of the detected DEGs. Moreover, 3 and 16 DEGs were detected within the QTL regions previously described in chromosomes 11 and 5, respectively. In addition to these two specific genomic regions, we also observde large transcriptomic changes from genes spread across the genome in the resistant plants in response to the virus infection. This early response against WMV implied genes involved in plant-pathogen interaction, plant hormone signal transduction, the MAPK signaling pathway or ubiquitin mediated proteolysis, in detriment to the photosynthetic and basal metabolites pathways. Moreover, the gene MELO3C021395, which coded a mediator of RNA polymerase II transcription subunit 33A (MED33A), has been proposed as the candidate gene located on chromosome 11 conferring resistance to WMV.

CONCLUSIONS

The comparative transcriptomic analysis presented here showed that, even though the resistance to WMV in TGR-1551 has a recessive nature, it triggers an active defense response at a transcriptomic level, which involves broad-spectrum resistance mechanisms. Thus, this study represents a step forward on our understanding of the mechanisms underlaying WMV resistance in melon. In addition, it sheds light into a broader topic on the mechanisms of recessive resistances.

Viruses and Cajal Bodies: A Critical Cellular Target in Virus Infection?

Nuclear bodies (NBs) are dynamic structures present in eukaryotic cell nuclei. They are not bounded by membranes and are often considered biomolecular condensates, defined structurally and functionally by the localisation of core components. Nuclear architecture can be reorganised during normal cellular processes such as the cell cycle as well as in response to cellular stress. Many plant and animal viruses target their proteins to NBs, in some cases triggering their structural disruption and redistribution. Although not all such interactions have been well characterised, subversion of NBs and their functions may form a key part of the life cycle of eukaryotic viruses that require the nucleus for their replication. This review will focus on Cajal bodies (CBs) and the viruses that target them. Since CBs are dynamic structures, other NBs (principally nucleoli and promyelocytic leukaemia, PML and bodies), whose components interact with CBs, will also be considered. As well as providing important insights into key virus-host cell interactions, studies on Cajal and associated NBs may identify novel cellular targets for development of antiviral compounds.

N-terminal acetylation of the βC1 protein encoded by the betasatellite of tomato yellow leaf curl China virus is critical for its viral pathogenicity

N-terminal acetylation (N-acetylation) is one of the most common protein modifications and plays crucial roles in viability and stress responses in animals and plants. However, very little is known about N-acetylation of viral proteins. Here, we identified the Thr residue at position 2 (Thr-2) in the βC1 protein encoded by the betasatellite of tomato yellow leaf curl China virus (TYLCCNB-βC1) as a novel N-acetylation site. Furthermore, the effects of TYLCCNB-βC1 N-acetylation on its function as a pathogenicity factor were determined via N-acetylation mutants in Nicotiana benthamiana plants. We found that N-acetylation of TYLCCNB-βC1 is critical for its self-interaction in the nucleus and viral pathogenesis, and that removal of N-acetylation of TYLCCNB-βC1 attenuated tomato yellow leaf curl China virus-induced symptoms and led to accelerated degradation of TYLCCNB-βC1 through the ubiquitin-proteasome system. Our data reveal a protective effect of N-acetylation of TYLCCNB-βC1 on its pathogenesis and demonstrate an antagonistic crosstalk between N-acetylation and ubiquitination in this geminiviral protein.

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.

RNA-seq analysis reveals an early defense response to tomato leaf curl New Delhi virus in potato cultivar Kufri Bahar

Potatoes in India are very susceptible to apical leaf curl disease, which causes severe symptoms and greater yield losses. Because the majority of potato cultivars are susceptible to the virus, it is crucial to discover sources of resistance and investigate the mechanism of resistance/susceptibility in potato cultivars. In this study, the gene expression profile of two potato cultivars, Kufri Bahar (resistant) and Kufri Pukhraj (susceptible), varying in their level of resistance to ToLCNDV, was analyzed using RNA-Seq. The Ion ProtonTM system was used to sequence eight RiboMinus RNA libraries from inoculated and uninoculated potato plants at 15 and 20 days after inoculation (DAI). The findings indicated that the majority of differentially expressed genes (DEGs) were cultivar-or time-specific. These DEGs included genes for proteins that interact with viruses, genes linked with the cell cycle, genes for proteins involved in defense, transcription and translation initiation factors, and plant hormone signaling pathway genes. Interestingly, defense responses were generated early in Kufri Bahar, at 15 DAI, which may have impeded the replication and spread of ToLCNDV. This research provides a genome-wide transcriptional analysis of two potato cultivars with variable levels of ToLCNDV resistance. At an early stage, we observed suppression of genes that interact with viral proteins, induction of genes associated with restriction of cell division, genes encoding defense proteins, AP2/ERF transcription factors, and altered expression of zinc finger protein genes, HSPs, JA, and SA pathway-related genes. Our findings add to a greater comprehension of the molecular basis of potato resistance to ToLCNDV and may aid in the development of more effective disease management techniques.

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.

Fungal F8-Culture Filtrate Induces Tomato Resistance against Tomato Yellow Leaf Curl Thailand Virus

Tomato (Solanum lycopersicum) is an important economic crop worldwide. However, tomato production is jeopardized by the devastating tomato yellow leaf curl disease caused by whitefly-transmitted begomoviruses (WTBs). In this study, we evaluated the efficacy of our previously developed plant antiviral immunity inducer, fungal F8-culture filtrate, on tomato to combat tomato yellow leaf curl Thailand virus (TYLCTHV), the predominant WTB in Taiwan. Our results indicated that F8-culture filtrate treatment induced strong resistance, did not reduce the growth of tomato, and induced prominent resistance against TYLCTHV both in the greenhouse and in the field. Among TYLCTHV-inoculated Yu-Nu tomato grown in the greenhouse, a greater percentage of plants treated with F8-culture filtrate (43-100%) were healthy-looking compared to the H₂O control (0-14%). We found that TYLCTHV cannot move systemically only on the F8-culture filtrate pretreated healthy-looking plants. Tracking the expression of phytohormone-mediated immune maker genes revealed that F8-culture filtrate mainly induced salicylic acid-mediated plant immunity. Furthermore, callose depositions and the expression of the pathogen-induced callose synthase gene, POWDERY MILDEW RESISTANT 4 were only strongly induced by TYLCTHV on tomato pretreated with F8-culture filtrate. This study provides an effective way to induce tomato resistance against TYLCTHV.

Plant responses to geminivirus infection: guardians of the plant immunity

BACKGROUND

Geminiviruses are circular, single-stranded viruses responsible for enormous crop loss worldwide. Rapid expansion of geminivirus diversity outweighs the continuous effort to control its spread. Geminiviruses channelize the host cell machinery in their favour by manipulating the gene expression, cell signalling, protein turnover, and metabolic reprogramming of plants. As a response to viral infection, plants have evolved to deploy various strategies to subvert the virus invasion and reinstate cellular homeostasis.

MAIN BODY

Numerous reports exploring various aspects of plant-geminivirus interaction portray the subtlety and flexibility of the host-pathogen dynamics. To leverage this pool of knowledge towards raising antiviral resistance in host plants, a comprehensive account of plant's defence response against geminiviruses is required. This review discusses the current knowledge of plant's antiviral responses exerted to geminivirus in the light of resistance mechanisms and the innate genetic factors contributing to the defence. We have revisited the defence pathways involving transcriptional and post-transcriptional gene silencing, ubiquitin-proteasomal degradation pathway, protein kinase signalling cascades, autophagy, and hypersensitive responses. In addition, geminivirus-induced phytohormonal fluctuations, the subsequent alterations in primary and secondary metabolites, and their impact on pathogenesis along with the recent advancements of CRISPR-Cas9 technique in generating the geminivirus resistance in plants have been discussed.

CONCLUSIONS

Considering the rapid development in the field of plant-virus interaction, this review provides a timely and comprehensive account of molecular nuances that define the course of geminivirus infection and can be exploited in generating virus-resistant plants to control global agricultural damage.

The Ubiquitin Proteasome System as a Double Agent in Plant-Virus Interactions

The ubiquitin proteasome is a rapid, adaptive mechanism for selective protein degradation, crucial for proper plant growth and development. The ubiquitin proteasome system (UPS) has also been shown to be an integral part of plant responses to stresses, including plant defence against pathogens. Recently, significant progress has been made in the understanding of the involvement of the UPS in the signalling and regulation of the interaction between plants and viruses. This review aims to discuss the current knowledge about the response of plant viral infection by the UPS and how the viruses counteract this system, or even use it for their own benefit.

The Global Dimension of Tomato Yellow Leaf Curl Disease: Current Status and Breeding Perspectives

Tomato yellow leaf curl disease (TYLCD) caused by tomato yellow leaf curl virus (TYLCV) and a group of related begomoviruses is an important disease which in recent years has caused serious economic problems in tomato (Solanum lycopersicum) production worldwide. Spreading of the vectors, whiteflies of the Bemisia tabaci complex, has been responsible for many TYLCD outbreaks. In this review, we summarize the current knowledge of TYLCV and TYLV-like begomoviruses and the driving forces of the increasing global significance through rapid evolution of begomovirus variants, mixed infection in the field, association with betasatellites and host range expansion. Breeding for host plant resistance is considered as one of the most promising and sustainable methods in controlling TYLCD. Resistance to TYLCD was found in several wild relatives of tomato from which six TYLCV resistance genes (Ty-1 to Ty-6) have been identified. Currently, Ty-1 and Ty-3 are the primary resistance genes widely used in tomato breeding programs. Ty-2 is also exploited commercially either alone or in combination with other Ty-genes (i.e., Ty-1, Ty-3 or ty-5). Additionally, screening of a large collection of wild tomato species has resulted in the identification of novel TYLCD resistance sources. In this review, we focus on genetic resources used to date in breeding for TYLCVD resistance. For future breeding strategies, we discuss several leads in order to make full use of the naturally occurring and engineered resistance to mount a broad-spectrum and sustainable begomovirus resistance.

RNA-Seq Transcriptome Analysis Provides Candidate Genes for Resistance to Tomato Leaf Curl New Delhi Virus in Melon

Tomato leaf curl New Delhi virus (ToLCNDV) emerged in the Mediterranean Basin in 2012 as the first DNA bipartite begomovirus (Geminiviridae family), causing severe yield and economic losses in cucurbit crops. A major resistance locus was identified in the wild melon accession WM-7 (Cucumis melo kachri group), but the mechanisms involved in the resistant response remained unknown. In this work, we used RNA-sequencing to identify disease-associated genes that are differentially expressed in the course of ToLCNDV infection and could contribute to resistance. Transcriptomes of the resistant WM-7 genotype and the susceptible cultivar Piñonet Piel de Sapo (PS) (C. melo ibericus group) in ToLCNDV and mock inoculated plants were compared at four time points during infection (0, 3, 6, and 12 days post inoculation). Different gene expression patterns were observed over time in the resistant and susceptible genotypes in comparison to their respective controls. Differentially expressed genes (DEGs) in ToLCNDV-infected plants were classified using gene ontology (GO) terms, and genes of the categories transcription, DNA replication, and helicase activity were downregulated in WM-7 but upregulated in PS, suggesting that reduced activity of these functions reduces ToLCNDV replication and intercellular spread and thereby contributes to resistance. DEGs involved in the jasmonic acid signaling pathway, photosynthesis, RNA silencing, transmembrane, and sugar transporters entail adverse consequences for systemic infection in the resistant genotype, and lead to susceptibility in PS. The expression levels of selected candidate genes were validated by qRT-PCR to corroborate their differential expression upon ToLCNDV infection in resistant and susceptible melon. Furthermore, single nucleotide polymorphism (SNPs) with an effect on structural functionality of DEGs linked to the main QTLs for ToLCNDV resistance have been identified. The obtained results pinpoint cellular functions and candidate genes that are differentially expressed in a resistant and susceptible melon line in response to ToLCNDV, an information of great relevance for breeding ToLCNDV-resistant melon cultivars.

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.

Tomato yellow leaf curl virus (TYLCV)-resistant tomatoes share molecular mechanisms sustaining resistance with their wild progenitor Solanum habrochaites but not with TYLCV-susceptible tomatoes

The wild tomato species Solanum habrochaites (Sh) has been used as a source for tomato yellow leaf curl virus (TYLCV) resistance in a breeding program to generate a TYLCV-resistant tomato line. Susceptible (S) and resistant (R) lines have been developed through this program. We compared the behavior of R, S and Sh tomato plants upon infection to find out whether the resistant phenotype of R plants originated from Sh. Results showed that mechanisms involving sugar-signaling (i.e., LIN6/HT1), water channels (i.e., TIP1;1), hormone homeostasis (i.e., ABA and SA) and urea accumulation were shared by S. habrochaites and R plants, but not by S. habrochaites and S tomatoes. This finding supports the hypothesis that these mechanisms were introgressed in the R genotype from the wild tomato progenitor during breeding for TYLCV resistance. Hence, identification of genes contributing to resistance to biotic stress from wild tomato species and their introgression into domestic plants ensures tomato supply and food security.

Superinfection by PHYVV Alters the Recovery Process in PepGMV-Infected Pepper Plants

Geminiviruses are important plant pathogens that affect crops around the world. In some geminivirus-host interactions, infected plants show recovery, a phenomenon characterized by symptom disappearance in newly emerging leaves. In pepper-Pepper golden mosaic virus (PepGMV) interaction, the host recovery process involves a silencing mechanism that includes both post-transcriptional (PTGS) and transcriptional (TGS) gene silencing pathways. Under field conditions, PepGMV is frequently found in mixed infections with Pepper huasteco yellow vein virus (PHYVV), another bipartite begomovirus. Mixed infected plants generally show a synergetic phenomenon and do not present recovery. Little is known about the molecular mechanism of this interaction. In the present study, we explored the effect of superinfection by PHYVV on a PepGMV-infected pepper plant showing recovery. Superinfection with PHYVV led to (a) the appearance of severe symptoms, (b) an increase of the levels of PepGMV DNA accumulation, (c) a decrease of the relative methylation levels of PepGMV DNA, and (d) an increase of chromatin activation marks present in viral minichromosomes. Finally, using heterologous expression and silencing suppression reporter systems, we found that PHYVV REn presents TGS silencing suppressor activity, whereas similar experiments suggest that Rep might be involved in suppressing PTGS.

Changes in maize transcriptome in response to maize Iranian mosaic virus infection

Background

Maize Iranian mosaic virus (MIMV, genus Nucleorhabdovirus, family Rhabdoviridae) causes an economically important disease in maize and other gramineous crops in Iran. MIMV negative-sense RNA genome sequence of 12,426 nucleotides has recently been completed. Maize Genetics and Genomics database shows that 39,498 coding genes and 4,976 non-coding genes of maize have been determined, but still some transcripts could not be annotated. The molecular host cell responses of maize to MIMV infection including differential gene expression have so far not been elucidated.

Methodology/Principal findings

Complementary DNA libraries were prepared from total RNA of MIMV-infected and mock-inoculated maize leaves and sequenced using Illumina HiSeq 2500. Cleaned raw transcript reads from MIMV-infected maize were mapped to reads from uninfected maize and to a maize reference genome. Differentially expressed transcripts were characterized by gene ontology and biochemical pathway analyses. Transcriptome data for selected genes were validated by real-time quantitative PCR.

Conclusion/Significance

Approximately 42 million clean reads for each treatment were obtained. In MIMV-infected maize compared to uninfected plants, 1689 transcripts were up-regulated and 213 transcripts were down-regulated. In response to MIMV infection, several pathways were activated in maize including immune receptor signaling, metabolic pathways, RNA silencing, hormone-mediated pathways, protein degradation, protein kinase and ATP binding activity, and fatty acid metabolism. Also, several transcripts including those encoding hydrophobic protein RCI2B, adenosylmethionine decarboxylase NAC transcription factor and nucleic acid binding, leucine-rich repeat, heat shock protein, 26S proteasome, oxidoreductases and endonuclease activity protein were up-regulated. These data will contribute to the identification of genes and pathways involved in plant-virus interactions that may serve as future targets for improved disease control.

Development of transcriptome based web genomic resources of yellow mosaic disease in Vigna mungo

Vigna mungo (Urdbean) is cultivated in the tropical and sub-tropical continental region of Asia. It is not only important source of dietary protein and nutritional elements, but also of immense value to human health due to medicinal properties. Yellow mosaic disease caused by Mungbean Yellow Mosaic India Virus is known to incur huge loss to crop, adversely affecting crop yield. Contrasting genotypes are ideal source for knowledge discovery of plant defence mechanism and associated candidate genes for varietal improvement. Whole genome sequence of this crop is yet to be completed. Moreover, genomic resources are also not freely accessible, thus available transcriptome data can be of immense use. V. mungo Transcriptome database, accessible at http://webtom.cabgrid.res.in/vmtdb/ has been developed using available data of two contrasting varieties viz., cv. VM84 (resistant) and cv. T9 (susceptible). De novo assembly was carried out using Trinity and CAP3. Out of total 240,945 unigenes, 165,894 (68.8%) showed similarity with known genes against NR database, and remaining 31.2% were found to be novel. We found 22,101 differentially expressed genes in all datasets, 44,335 putative genic SSR markers, 4105 SNPs and Indels, 64,964 transcriptional factor, 546 mature miRNA target prediction in 703 differentially expressed unigenes and 137 pathways. MAPK, salicylic acid-binding protein 2-like, pathogenesis-related protein and NBS-LRR domain were found which may play an important role in defence against pathogens. This is the first web genomic resource of V. mungo for future genome annotation as well as ready to use markers for future variety improvement program.

Geminiviruses and Plant Hosts: A Closer Examination of the Molecular Arms Race

Geminiviruses are plant-infecting viruses characterized by a single-stranded DNA (ssDNA) genome. Geminivirus-derived proteins are multifunctional and effective regulators in modulating the host cellular processes resulting in successful infection. Virus-host interactions result in changes in host gene expression patterns, reprogram plant signaling controls, disrupt central cellular metabolic pathways, impair plant's defense system, and effectively evade RNA silencing response leading to host susceptibility. This review summarizes what is known about the cellular processes in the continuing tug of war between geminiviruses and their plant hosts at the molecular level. In addition, implications for engineered resistance to geminivirus infection in the context of a greater understanding of the molecular processes are also discussed. Finally, the prospect of employing geminivirus-based vectors in plant genome engineering and the emergence of powerful genome editing tools to confer geminivirus resistance are highlighted to complete the perspective on geminivirus-plant molecular 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.

Characterization of virus-derived small interfering RNAs in Apple stem grooving virus-infected in vitro-cultured Pyrus pyrifolia shoot tips in response to high temperature treatment

Background

Heat treatment (known as thermotherapy) together with in vitro culture of shoot meristem tips is a commonly used technology to obtain virus-free germplasm for the effective control of virus diseases in fruit trees. RNA silencing as an antiviral defense mechanism has been implicated in this process. To understand if high temperature-mediated acceleration of the host antiviral gene silencing system in the meristem tip facilitates virus-derived small interfering RNAs (vsiRNA) accumulation to reduce the viral RNA titer in the fruit tree meristem tip cells, we used the Apple stem grooving virus (ASGV)–Pyrus pyrifolia pathosystem to explore the possible roles of vsiRNA in thermotherapy.

Results

At first we determined the full-length genome sequence of the ASGV-Js2 isolate and then profiled vsiRNAs in the meristem tip of in vitro-grown pear (cv. ‘Jinshui no. 2’) shoots infected by ASGV-Js2 and cultured at 24 and 37 °C. A total of 7,495 and 7,949 small RNA reads were obtained from the tips of pear shoots cultured at 24 and 37 °C, respectively. Mapping of the vsiRNAs to the ASGV-Js2 genome revealed that they were unevenly distributed along the ASGV-Js2 genome, and that 21- and 22-nt vsiRNAs preferentially accumulated at both temperatures. The 5′-terminal nucleotides of ASGV-specific siRNAs in the tips cultured under different temperatures had a similar distribution pattern, and the nucleotide U was the most frequent. RT-qPCR analyses suggested that viral genome accumulation was drastically compromised at 37 °C compared to 24 °C, which was accompanied with the elevated levels of vsiRNAs at 37 °C. As plant Dicer-like proteins (DCLs), Argonaute proteins (AGOs), and RNA-dependent RNA polymerases (RDRs) are implicated in vsiRNA biogenesis, we also cloned the partial sequences of PpDCL2,4, PpAGO1,2,4 and PpRDR1 genes, and found their expression levels were up-regulated in the ASGV-infected pear shoots at 37 °C.

Conclusions

Collectively, these results showed that upon high temperature treatment, the ASGV-infected meristem shoot tips up-regulated the expression of key genes in the RNA silencing pathway, induced the biogenesis of vsiRNAs and inhibited viral RNA accumulation. This study represents the first report on the characterization of the vsiRNA population in pear plants infected by ASGV-Js2, in response to high temperature treatment.

Electronic supplementary material

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

Next-Generation Sequencing and Genome Editing in Plant Virology

Next-generation sequencing (NGS) has been applied to plant virology since 2009. NGS provides highly efficient, rapid, low cost DNA, or RNA high-throughput sequencing of the genomes of plant viruses and viroids and of the specific small RNAs generated during the infection process. These small RNAs, which cover frequently the whole genome of the infectious agent, are 21-24 nt long and are known as vsRNAs for viruses and vd-sRNAs for viroids. NGS has been used in a number of studies in plant virology including, but not limited to, discovery of novel viruses and viroids as well as detection and identification of those pathogens already known, analysis of genome diversity and evolution, and study of pathogen epidemiology. The genome engineering editing method, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been successfully used recently to engineer resistance to DNA geminiviruses (family, Geminiviridae) by targeting different viral genome sequences in infected Nicotiana benthamiana or Arabidopsis plants. The DNA viruses targeted include tomato yellow leaf curl virus and merremia mosaic virus (begomovirus); beet curly top virus and beet severe curly top virus (curtovirus); and bean yellow dwarf virus (mastrevirus). The technique has also been used against the RNA viruses zucchini yellow mosaic virus, papaya ringspot virus and turnip mosaic virus (potyvirus) and cucumber vein yellowing virus (ipomovirus, family, Potyviridae) by targeting the translation initiation genes eIF4E in cucumber or Arabidopsis plants. From these recent advances of major importance, it is expected that NGS and CRISPR-Cas technologies will play a significant role in the very near future in advancing the field of plant virology and connecting it with other related fields of biology.

CLCuMuB βC1 Subverts Ubiquitination by Interacting with NbSKP1s to Enhance Geminivirus Infection in Nicotiana benthamiana

Viruses interfere with and usurp host machinery and circumvent defense responses to create a suitable cellular environment for successful infection. This is usually achieved through interactions between viral proteins and host factors. Geminiviruses are a group of plant-infecting DNA viruses, of which some contain a betasatellite, known as DNAβ. Here, we report that Cotton leaf curl Multan virus (CLCuMuV) uses its sole satellite-encoded protein βC1 to regulate the plant ubiquitination pathway for effective infection. We found that CLCuMu betasatellite (CLCuMuB) βC1 interacts with NbSKP1, and interrupts the interaction of NbSKP1s with NbCUL1. Silencing of either NbSKP1s or NbCUL1 enhances the accumulation of CLCuMuV genomic DNA and results in severe disease symptoms in plants. βC1 impairs the integrity of SCFCOI1 and the stabilization of GAI, a substrate of the SCFSYL1 to hinder responses to jasmonates (JA) and gibberellins (GA). Moreover, JA treatment reduces viral accumulation and symptoms. These results suggest that CLCuMuB βC1 inhibits the ubiquitination function of SCF E3 ligases through interacting with NbSKP1s to enhance CLCuMuV infection and symptom induction in plants.

The receptor-like kinase SlSOBIR1 is differentially modulated by virus infection but its overexpression in tobacco has no significant impact on virus accumulation

KEY MESSAGE

The role of the tomato receptor-like kinase SlSOBIR1 in antiviral defense was investigated. SlSOBIR1 was transcriptionally modulated by unrelated viruses but its ectopic expression had no effect on virus accumulation. Leucine-rich repeat receptor-like kinases (LRR-RLK) constitute a diverse group of proteins allowing the cell to recognize and respond to the extracellular environment. In the present study we focused on a gene encoding a tomato LRR-RLK (named SlSOBIR1) involved in the host defense against fungal pathogens. Curiously, SlSOBIR1 has been previously reported to be down-regulated by Pepper yellow mosaic virus (PepYMV) infection. Here, we show that SlSOBIR1 is responsive to wounding and differentially modulated by unrelated virus infection, i.e., down-regulated by PepYMV and up-regulated by Tomato chlorotic spot virus (TCSV). Despite these divergent expression profiles, SlSOBIR1 overexpression in transgenic tobacco plants had no evident effect on TCSV and PepYMV accumulation. On the other hand, overexpression of SlSOBIR1 significantly increased the expression of selected defense genes (PR-1a and PR-6) and exacerbated superoxide production in wounded leaves. Our data indicate that the observed modulation of SlSOBIR1 expression is probably triggered by secondary effects of the virus infection process and suggest that SlSOBIR1 is not directly involved in antiviral signaling response.

Quantitative proteomics identifies 38 proteins that are differentially expressed in cucumber in response to cucumber green mottle mosaic virus infection

BACKGROUND

Since it was first reported in 1935, Cucumber green mottle mosaic virus (CGMMV) has become a serious pathogen in a range of cucurbit crops. The virus is generally transmitted by propagation materials, and to date no effective chemical or cultural methods of control have been developed to combat its spread. The current study presents a preliminary analysis of the pathogenic mechanisms from the perspective of protein expression levels in an infected cucumber host, with the objective of elucidating the infection process and potential strategies to reduce both the economic and yield losses associated with CGMMV.

METHODS

Isobaric tags for relative and absolute quantitation (iTRAQ) technology coupled with liquid chromatography-tandem mass spectrometric (LC-MS/MS) were used to identify the differentially expressed proteins in cucumber plants infected with CGMMV compared with mock-inoculated plants. The functions of the proteins were deduced by functional annotation and their involvement in metabolic processes explored by KEGG pathway analysis to identify their interactions during CGMMV infection, while their in vivo expression was further verified by qPCR.

RESULTS

Infection by CGMMV altered both the expression level and absolute quantity of 38 proteins (fold change >0.6) in cucumber hosts. Of these, 23 were found to be up-regulated, while 15 were down-regulated. Gene ontology (GO) analysis revealed that 22 of the proteins had a combined function and were associated with molecular function (MF), biological process (BP) and cellular component (CC). Several other proteins had a dual function with 1, 7, and 2 proteins being associated with BP/CC, BP/MF, CC/MF, respectively. The remaining 3 proteins were only involved in MF. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 18 proteins that were involved in 13 separate metabolic pathways. These pathways were subsequently merged to generate three network diagrams illustrating the interactions between the different pathways, while qPCR was used to track the changes in expression levels of the proteins identified at 3 time points during CGMMV infection. Taken together these results greatly expand our understanding of the relationships between CGMMV and cucumber hosts.

CONCLUSIONS

The results of the study indicate that CGMMV infection significantly changes the physiology of cucumbers, affecting the expression levels of individual proteins as well as entire metabolic pathways. The bioinformatic analysis also identified several pathogenesis-related (PR) proteins that could be useful in the development of disease-resistant plants.

Mechanism of regulation of tomato TRN1 gene expression in late infection with tomato leaf curl New Delhi virus (ToLCNDV)

Tomato leaf curl disease caused by geminiviruses is manifested by curling and puckering of leaves and thickening of veins, resembling developmental defects. This is probably due to the long-term altered regulation of expression of development related gene(s). Our results show that in the infected leaves the transcript level of TORNADO1 (SlTRN1), a gene important for cell expansion and vein formation, increased significantly. SlTRN1 is transcribed from two start sites. The preferential usage of one start site governs its expression in viral-stressed plants. To investigate the role of specific promoter elements in mediating differential expression of SlTRN1, we performed SlTRN1 promoter analysis. The promoter-regulatory sequences harbor multiple W-boxes. The SlWRKY16 transcription factor actively interacts with one of the W-boxes. WRKY proteins are commonly induced by salicylic acid (SA), and consequently SA treatment increased transcript level of SlWRKY16 and SlTRN1. Further mutational analyses confirmed the role of W-boxes in mediating SlTRN1 induction during ToLCNDV infection or SA treatment. We postulate that the activation of SA pathway during stress-response in tomato induces WRKY16, which in turn modulates transcription of SlTRN1 gene. This study unravels the mechanism of regulation of a developmental gene during stress-response, which may affect the severity of symptoms.

Reference Gene Selection for qPCR Analysis in Tomato-Bipartite Begomovirus Interaction and Validation in Additional Tomato-Virus Pathosystems

Quantitative Polymerase Chain Reaction (qPCR) is currently the most sensitive technique used for absolute and relative quantification of a target gene transcript, requiring the use of appropriated reference genes for data normalization. To accurately estimate the relative expression of target tomato (Solanum lycopersicum L.) genes responsive to several virus species in reverse transcription qPCR analysis, the identification of reliable reference genes is mandatory. In the present study, ten reference genes were analyzed across a set of eight samples: two tomato contrasting genotypes ('Santa Clara', susceptible, and its near-isogenic line 'LAM 157', resistant); subjected to two treatments (inoculation with Tomato chlorotic mottle virus (ToCMoV) and its mock-inoculated control) and in two distinct times after inoculation (early and late). Reference genes stability was estimated by three statistical programs (geNorm, NormFinder and BestKeeper). To validate the results over broader experimental conditions, a set of ten samples, corresponding to additional three tomato-virus pathosystems that included tospovirus, crinivirus and tymovirus + tobamovirus, was analyzed together with the tomato-ToCMoV pathosystem dataset, using the same algorithms. Taking into account the combined analyses of the ranking order outputs from the three algorithms, TIP41 and EF1 were identified as the most stable genes for tomato-ToCMoV pathosystem, and TIP41 and EXP for the four pathosystems together, and selected to be used as reference in the forthcoming expression qPCR analysis of target genes in experimental conditions involving the aforementioned tomato-virus pathosystems.

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.

Mapping in the era of sequencing: high density genotyping and its application for mapping TYLCV resistance in Solanum pimpinellifolium

BACKGROUND

A RIL population between Solanum lycopersicum cv. Moneymaker and S. pimpinellifolium G1.1554 was genotyped with a custom made SNP array. Additionally, a subset of the lines was genotyped by sequencing (GBS).

RESULTS

A total of 1974 polymorphic SNPs were selected to develop a linkage map of 715 unique genetic loci. We generated plots for visualizing the recombination patterns of the population relating physical and genetic positions along the genome.This linkage map was used to identify two QTLs for TYLCV resistance which contained favourable alleles derived from S. pimpinellifolium. Further GBS was used to saturate regions of interest, and the mapping resolution of the two QTLs was improved. The analysis showed highest significance on Chromosome 11 close to the region of 51.3 Mb (qTy-p11) and another on Chromosome 3 near 46.5 Mb (qTy-p3). Furthermore, we explored the population using untargeted metabolic profiling, and the most significant differences between susceptible and resistant plants were mainly associated with sucrose and flavonoid glycosides.

CONCLUSIONS

The SNP information obtained from an array allowed a first QTL screening of our RIL population. With additional SNP data of a RILs subset, obtained through GBS, we were able to perform an in silico mapping improvement to further confirm regions associated with our trait of interest. With the combination of different ~ omics platforms we provide valuable insight into the genetics of S. pimpinellifolium-derived TYLCV resistance.

Comparative transcriptome profiling of a resistant vs. susceptible tomato (Solanum lycopersicum) cultivar in response to infection by tomato yellow leaf curl virus

Tomato yellow leaf curl virus (TYLCV) threatens tomato production worldwide by causing leaf yellowing, leaf curling, plant stunting and flower abscission. The current understanding of the host plant defense response to this virus is very limited. Using whole transcriptome sequencing, we analyzed the differential gene expression in response to TYLCV infection in the TYLCV-resistant tomato breeding line CLN2777A (R) and TYLCV-susceptible tomato breeding line TMXA48-4-0 (S). The mixed inoculated samples from 3, 5 and 7 day post inoculation (dpi) were compared to non-inoculated samples at 0 dpi. Of the total of 34831 mapped transcripts, 209 and 809 genes were differentially expressed in the R and S tomato line, respectively. The proportion of up-regulated differentially expressed genes (DEGs) in the R tomato line (58.37%) was higher than that in the S line (9.17%). Gene ontology (GO) analyses revealed that similar GO terms existed in both DEGs of R and S lines; however, some sets of defense related genes and their expression levels were not similar between the two tomato lines. Genes encoding for WRKY transcriptional factors, R genes, protein kinases and receptor (-like) kinases which were identified as down-regulated DEGs in the S line were up-regulated or not differentially expressed in the R line. The up-regulated DEGs in the R tomato line revealed the defense response of tomato to TYLCV infection was characterized by the induction and regulation of a series of genes involved in cell wall reorganization, transcriptional regulation, defense response, ubiquitination, metabolite synthesis and so on. The present study provides insights into various reactions underlining the successful establishment of resistance to TYLCV in the R tomato line, and helps in the identification of important defense-related genes in tomato for TYLCV disease management.