Endogenous pararetroviral sequences (EPRVs) are a recently discovered class of repetitive sequences that is broadly distributed in the plant kingdom. The potential contribution of EPRVs to plant pathogenicity or, conversely, to virus resistance is just beginning to be explored. Some members of the family Solanaceae are particularly rich in EPRVs. In previous work, EPRVs have been characterized molecularly in various species of Nicotiana including N.tabacum (tobacco) and Solanum tuberosum (potato). Here we describe a family of EPRVs in cultivated tomato (Solanum lycopersicum L.) and a wild relative (S.habrochaites).
Molecular cloning and DNA sequence analysis revealed that tomato EPRVs (named LycEPRVs) are most closely related to those in tobacco. The sequence similarity of LycEPRVs in S.lycopersicum and S.habrochaites indicates they are potentially derived from the same pararetrovirus. DNA blot analysis revealed a similar genomic organization in the two species, but also some independent excision or insertion events after species separation, or flanking sequence divergence. LycEPRVs share with the tobacco elements a disrupted genomic structure and frequent association with retrotransposons. Fluorescence in situ hybridization revealed that copies of LycEPRV are dispersed on all chromosomes in predominantly heterochromatic regions. Methylation of LycEPRVs was detected in CHG and asymmetric CHH nucleotide groups. Although normally quiescent EPRVs can be reactivated and produce symptoms of infection in some Nicotiana interspecific hybrids, a similar pathogenicity of LycEPRVs could not be demonstrated in Solanum L. section Lycopersicon [Mill.] hybrids. Even in healthy plants, however, transcripts derived from multiple LycEPRV loci and short RNAs complementary to LycEPRVs were detected and were elevated upon infection with heterologous viruses encoding suppressors of PTGS.
The analysis of LycEPRVs provides further evidence for the extensive invasion of pararetroviral sequences into the genomes of solanaceous plants. The detection of asymmetric CHH methylation and short RNAs, which are hallmarks of RNAi in plants, suggests that LycEPRVs are controlled by an RNA-mediated silencing mechanism.
Non-retroviral RNA virus sequences (NRVSs) have been found in the chromosomes of vertebrates and fungi, but not plants. Here we report similarly endogenized NRVSs derived from plus-, negative-, and double-stranded RNA viruses in plant chromosomes. These sequences were found by searching public genomic sequence databases, and, importantly, most NRVSs were subsequently detected by direct molecular analyses of plant DNAs. The most widespread NRVSs were related to the coat protein (CP) genes of the family Partitiviridae which have bisegmented dsRNA genomes, and included plant- and fungus-infecting members. The CP of a novel fungal virus (Rosellinia necatrix partitivirus 2, RnPV2) had the greatest sequence similarity to Arabidopsis thaliana ILR2, which is thought to regulate the activities of the phytohormone auxin, indole-3-acetic acid (IAA). Furthermore, partitivirus CP-like sequences much more closely related to plant partitiviruses than to RnPV2 were identified in a wide range of plant species. In addition, the nucleocapsid protein genes of cytorhabdoviruses and varicosaviruses were found in species of over 9 plant families, including Brassicaceae and Solanaceae. A replicase-like sequence of a betaflexivirus was identified in the cucumber genome. The pattern of occurrence of NRVSs and the phylogenetic analyses of NRVSs and related viruses indicate that multiple independent integrations into many plant lineages may have occurred. For example, one of the NRVSs was retained in Ar. thaliana but not in Ar. lyrata or other related Camelina species, whereas another NRVS displayed the reverse pattern. Our study has shown that single- and double-stranded RNA viral sequences are widespread in plant genomes, and shows the potential of genome integrated NRVSs to contribute to resolve unclear phylogenetic relationships of plant species.
Eukaryotic genomes contain sequences that have originated from DNA viruses and reverse-transcribing viruses, i.e., retroviruses, pararetroviruses (DNA viruses), and transposons. However, the sequences of non-retroviral RNA viruses, which are unable to convert their genomes to DNA, were until recently considered not to be integrated into eukaryotic nuclear genomes. We present evidence for multiple independent events of horizontal gene transfer from a wide range of RNA viruses, including plus-sense, minus-sense, and double-stranded RNA viruses, into the genomes of distantly related plant lineages. Some non-retroviral integrated RNA viral sequences are conserved across genera within a plant family, whereas others are retained only in a limited number of species in a genus. Integration profiles of non-retroviral integrated RNA viral sequences demonstrate the potential of these sequences to serve as powerful molecular tools for deciphering phylogenetic relationships among related plants. Moreover, this study highlights plants co-opting non-retroviral RNA virus sequences, and provides insights into plant genome evolution and interplay between non-reverse-transcribing RNA viruses and their hosts.
Cauliflower mosaic virus (CaMV) and Rice tungro bacilliform virus (RTBV) belong to distinct genera of pararetroviruses infecting dicot and monocot plants, respectively. In both viruses, polycistronic translation of pregenomic (pg) RNA is initiated by shunting ribosomes that bypass a large region of the pgRNA leader with several short (s)ORFs and a stable stem-loop structure. The shunt requires translation of a 5′-proximal sORF terminating near the stem. In CaMV, mutations knocking out this sORF nearly abolish shunting and virus viability.
Here we show that two distant regions of the CaMV leader that form a minimal shunt configuration comprising the sORF, a bottom part of the stem, and a shunt landing sequence can be replaced by heterologous sequences that form a structurally similar configuration in RTBV without any dramatic effect on shunt-mediated translation and CaMV infectivity. The CaMV-RTBV chimeric leader sequence was largely stable over five viral passages in turnip plants: a few alterations that did eventually occur in the virus progenies are indicative of fine tuning of the chimeric sequence during adaptation to a new host.
Our findings demonstrate cross-species functionality of pararetroviral cis-elements driving ribosome shunting and evolutionary conservation of the shunt mechanism.
We are grateful to Matthias Müller and Sandra Pauli for technical assistance. This work was initiated at Friedrich Miescher Institute (Basel, Switzerland). We thank Prof. Thomas Boller for hosting the group at the Institute of Botany.
RNA silencing is a natural defence mechanism against viruses in plants, and transgenes expressing viral RNA-derived sequences were previously shown to confer silencing-based enhanced resistance against the cognate virus in several species. However, RNA silencing was shown to dysfunction at low temperatures in several species, questioning the relevance of this strategy in perennial plants such as grapevines, which are often exposed to low temperatures during the winter season. Here, we show that inverted-repeat (IR) constructs trigger a highly efficient silencing reaction in all somatic tissues in grapevines. Similarly to other plant species, IR-derived siRNAs trigger production of secondary transitive siRNAs. However, and in sharp contrast to other species tested to date where RNA silencing is hindered at low temperature, this process remained active in grapevine cultivated at 4°C. Consistently, siRNA levels remained steady in grapevines cultivated between 26°C and 4°C, whereas they are severely decreased in Arabidopsis grown at 15°C and almost undetectable at 4°C. Altogether, these results demonstrate that RNA silencing operates in grapevine in a conserved manner but is resistant to far lower temperatures than ever described in other species.
A graft-transmissible disease displaying red veins, red blotches and total reddening of leaves in red-berried wine grape (Vitis vinifera L.) cultivars was observed in commercial vineyards. Next-generation sequencing technology was used to identify etiological agent(s) associated with this emerging disease, designated as grapevine redleaf disease (GRD). High quality RNA extracted from leaves of grape cultivars Merlot and Cabernet Franc with and without GRD symptoms was used to prepare cDNA libraries. Assembly of highly informative sequence reads generated from Illumina sequencing of cDNA libraries, followed by bioinformatic analyses of sequence contigs resulted in specific identification of taxonomically disparate viruses and viroids in samples with and without GRD symptoms. A single-stranded DNA virus, tentatively named Grapevine redleaf-associated virus (GRLaV), and Grapevine fanleaf virus were detected only in grapevines showing GRD symptoms. In contrast, Grapevine rupestris stem pitting-associated virus, Hop stunt viroid, Grapevine yellow speckle viroid 1, Citrus exocortis viroid and Citrus exocortis Yucatan viroid were present in both symptomatic and non-symptomatic grapevines. GRLaV was transmitted by the Virginia creeper leafhopper (Erythroneura ziczac Walsh) from grapevine-to-grapevine under greenhouse conditions. Molecular and phylogenetic analyses indicated that GRLaV, almost identical to recently reported Grapevine Cabernet Franc-associated virus from New York and Grapevine red blotch-associated virus from California, represents an evolutionarily distinct lineage in the family Geminiviridae with genome characteristics distinct from other leafhopper-transmitted geminiviruses. GRD significantly reduced fruit yield and affected berry quality parameters demonstrating negative impacts of the disease. Higher quantities of carbohydrates were present in symptomatic leaves suggesting their possible role in the expression of redleaf symptoms.
Virus-infected plants accumulate abundant, 21–24 nucleotide viral siRNAs which are generated by the evolutionary conserved RNA interference (RNAi) machinery that regulates gene expression and defends against invasive nucleic acids. Here we show that, similar to RNA viruses, the entire genome sequences of DNA viruses are densely covered with siRNAs in both sense and antisense orientations. This implies pervasive transcription of both coding and non-coding viral DNA in the nucleus, which generates double-stranded RNA precursors of viral siRNAs. Consistent with our finding and hypothesis, we demonstrate that the complete genomes of DNA viruses from Caulimoviridae and Geminiviridae families can be reconstructed by deep sequencing and de novo assembly of viral siRNAs using bioinformatics tools. Furthermore, we prove that this ‘siRNA omics’ approach can be used for reliable identification of the consensus master genome and its microvariants in viral quasispecies. Finally, we utilized this approach to reconstruct an emerging DNA virus and two viroids associated with economically-important red blotch disease of grapevine, and to rapidly generate a biologically-active clone representing the wild type master genome of Oilseed rape mosaic virus. Our findings show that deep siRNA sequencing allows for de novo reconstruction of any DNA or RNA virus genome and its microvariants, making it suitable for universal characterization of evolving viral quasispecies as well as for studying the mechanisms of siRNA biogenesis and RNAi-based antiviral defense.
Rice tungro disease is caused by synergistic interaction of an RNA picorna-like virus Rice tungro spherical virus (RTSV) and a DNA pararetrovirus Rice tungro bacilliform virus (RTBV). It is spread by insects owing to an RTSV-encoded transmission factor. RTBV has evolved a ribosome shunt mechanism to initiate translation of its pregenomic RNA having a long and highly structured leader. We found that a long leader of RTSV genomic RNA remarkably resembles the RTBV leader: both contain several short ORFs (sORFs) and potentially fold into a large stem-loop structure with the first sORF terminating in front of the stem basal helix. Using translation assays in rice protoplasts and wheat germ extracts, we show that, like in RTBV, both initiation and proper termination of the first sORF translation in front of the stem are required for shunt-mediated translation of a reporter ORF placed downstream of the RTSV leader. The base pairing that forms the basal helix is required for shunting, but its sequence can be varied. Shunt efficiency in RTSV is lower than in RTBV. But in addition to shunting the RTSV leader sequence allows relatively efficient linear ribosome migration, which also contributes to translation initiation downstream of the leader. We conclude that RTSV and RTBV have developed a similar, sORF-dependent shunt mechanism possibly to adapt to the host translation system and/or coordinate their life cycles. Given that sORF-dependent shunting also operates in a pararetrovirus Cauliflower mosaic virus and likely in other pararetroviruses that possess a conserved shunt configuration in their leaders it is tempting to propose that RTSV may have acquired shunt cis-elements from RTBV during their co-existence.
Ribosome shunting, first discovered in plant pararetroviruses, is a translation initiation mechanism that combines 5′ end-dependent scanning and internal initiation and allows a bypass of highly-structured leaders of certain viral and cellular mRNAs. Here we demonstrate that a similar shunt mechanism has been developed by the RNA picorna-like virus RTSV and the DNA pararetrovirus RTBV that form a disease complex in rice. Leader sequences of the RTSV genomic RNA and the RTBV pregenomic RNA possess a conserved shunt configuration with a 5′-proximal short ORF (sORF1) terminating in front of a large stem-loop structure. Like in RTBV and a related pararetrovirus Cauliflower mosaic virus, shunt-mediated translation downstream of the RTSV leader depends on initiation and proper termination of sORF1 translation and on formation of the basal helix of the downstream secondary structure. Given that RTBV-like shunt elements with identical sequence motifs are present in all RTSV isolates but absent in related picorna-like viruses, it is likely that RTSV could have acquired these elements after its encounter with RTBV. Alternatively, the RTSV shunt elements could have evolved independently to adapt to the rice translation machinery. Our study highlights on-going genetic exchange and co-adaptation to the host in emerging viral disease complexes.
Endogenous plant pararetroviruses (EPRVs) are viral sequences of the family Caulimoviridae integrated into the nuclear genome of numerous plant species. The ability of some endogenous sequences of Banana streak viruses (eBSVs) in the genome of banana (Musa sp.) to induce infections just like the virus itself was recently demonstrated (P. Gayral et al., J. Virol. 83:6697-6710, 2008). Although eBSVs probably arose from accidental events, infectious eBSVs constitute an extreme case of parasitism, as well as a newly described strategy for vertical virus transmission in plants. We investigated the early evolutionary stages of infectious eBSV for two distinct BSV species—GF (BSGFV) and Imové (BSImV)—through the study of their distribution, insertion polymorphism, and structure evolution among selected banana genotypes representative of the diversity of 60 wild Musa species and genotypes. To do so, the historical frame of host evolution was analyzed by inferring banana phylogeny from two chloroplast regions—matK and trnL-trnF—as well as from the nuclear genome, using 19 microsatellite loci. We demonstrated that both BSV species integrated recently in banana evolution, circa 640,000 years ago. The two infectious eBSVs were subjected to different selective pressures and showed distinct levels of rearrangement within their final structure. In addition, the molecular phylogenies of integrated and nonintegrated BSVs enabled us to establish the phylogenetic origins of eBSGFV and eBSImV.
Grapevine leafroll disease (GLD) is one of the most important grapevine viral diseases affecting grapevines worldwide. The impact on vine health, crop yield, and quality is difficult to assess due to a high number of variables, but significant economic losses are consistently reported over the lifespan of a vineyard if intervention strategies are not implemented. Several viruses from the family Closteroviridae are associated with GLD. However, Grapevine leafroll-associated virus 3 (GLRaV-3), the type species for the genus Ampelovirus, is regarded as the most important causative agent. Here we provide a general overview on various aspects of GLRaV-3, with an emphasis on the latest advances in the characterization of the genome. The full genome of several isolates have recently been sequenced and annotated, revealing the existence of several genetic variants. The classification of these variants, based on their genome sequence, will be discussed and a guideline is presented to facilitate future comparative studies. The characterization of sgRNAs produced during the infection cycle of GLRaV-3 has given some insight into the replication strategy and the putative functionality of the ORFs. The latest nucleotide sequence based molecular diagnostic techniques were shown to be more sensitive than conventional serological assays and although ELISA is not as sensitive it remains valuable for high-throughput screening and complementary to molecular diagnostics. The application of next-generation sequencing is proving to be a valuable tool to study the complexity of viral infection as well as plant pathogen interaction. Next-generation sequencing data can provide information regarding disease complexes, variants of viral species, and abundance of particular viruses. This information can be used to develop more accurate diagnostic assays. Reliable virus screening in support of robust grapevine certification programs remains the cornerstone of GLD management.
grapevine leafroll disease; GLRaV-3; ampelovirus; Closteroviridae; genetic variants
We have characterized the virome in single grapevines by 454 high-throughput sequencing of double-stranded RNA recovered from the vine stem. The analysis revealed a substantial set of sequences similar to those of fungal viruses. Twenty-six putative fungal virus groups were identified from a single plant source. These represented half of all known mycoviral families including the Chrysoviridae, Hypoviridae, Narnaviridae, Partitiviridae, and Totiviridae. Three of the mycoviruses were associated with Botrytis cinerea, a common fungal pathogen of grapes. Most of the rest appeared to be undescribed. The presence of viral sequences identified by BLAST analysis was confirmed by sequencing PCR products generated from the starting material using primers designed from the genomic sequences of putative mycoviruses. To further characterize these sequences as fungal viruses, fungi from the grapevine tissue were cultured and screened with the same PCR probes. Five of the mycoviruses identified in the total grapevine extract were identified again in extracts of the fungal cultures.
Electronic supplementary material
The online version of this article (doi:10.1007/s00705-010-0869-8) contains supplementary material, which is available to authorized users.
We have used multiplexed high-throughput sequencing to characterize changes in small RNA populations that occur during viral infection in animal cells. Small RNA-based mechanisms such as RNA interference (RNAi) have been shown in plant and invertebrate systems to play a key role in host responses to viral infection. Although homologs of the key RNAi effector pathways are present in mammalian cells, and can launch an RNAi-mediated degradation of experimentally targeted mRNAs, any role for such responses in mammalian host-virus interactions remains to be characterized. Six different viruses were examined in 41 experimentally susceptible and resistant host systems. We identified virus-derived small RNAs (vsRNAs) from all six viruses, with total abundance varying from “vanishingly rare” (less than 0.1% of cellular small RNA) to highly abundant (comparable to abundant micro-RNAs “miRNAs”). In addition to the appearance of vsRNAs during infection, we saw a number of specific changes in host miRNA profiles. For several infection models investigated in more detail, the RNAi and Interferon pathways modulated the abundance of vsRNAs. We also found evidence for populations of vsRNAs that exist as duplexed siRNAs with zero to three nucleotide 3′ overhangs. Using populations of cells carrying a Hepatitis C replicon, we observed strand-selective loading of siRNAs onto Argonaute complexes. These experiments define vsRNAs as one possible component of the interplay between animal viruses and their hosts.
Short RNAs derived from invading viruses with RNA genomes are important components of antiviral immunity in plants, worms and flies. The regulated generation of these short RNAs, and their engagement by the immune apparatus, is essential for inhibiting viral growth in these organisms. Mammals have the necessary protein components to generate these viral-derived short RNAs (“vsRNAs”), raising the question of whether vsRNAs in mammals are a general feature of infections with RNA viruses. Our work with Hepatitis C, Polio, Dengue, Vesicular Stomatitis, and West Nile viruses in a broad host repertoire demonstrates the generality of RNA virus-derived vsRNA production, and the ability of the cellular short RNA apparatus to engage these vsRNAs in mammalian cells. Detailed analyses of vsRNA and host-derived short RNA populations demonstrate both common and virus-specific features of the interplay between viral infection and short RNA populations. The vsRNA populations described in this work represent a novel dimension in both viral pathogenesis and host response.
Grapevine fanleaf virus (GFLV) and Arabis mosaic virus (ArMV) from the genus Nepovirus, family Secoviridae, cause a severe degeneration of grapevines. GFLV and ArMV have a bipartite RNA genome and are transmitted specifically by the ectoparasitic nematodes Xiphinema index and Xiphinema diversicaudatum, respectively. The transmission specificity of both viruses maps to their respective RNA2-encoded coat protein (CP). To further delineate the GFLV CP determinants of transmission specificity, three-dimensional (3D) homology structure models of virions and CP subunits were constructed based on the crystal structure of Tobacco ringspot virus, the type member of the genus Nepovirus. The 3D models were examined to predict amino acids that are exposed at the external virion surface, highly conserved among GFLV isolates but divergent between GFLV and ArMV. Five short amino acid stretches that matched these topographical and sequence conservation criteria were selected and substituted in single and multiple combinations by their ArMV counterparts in a GFLV RNA2 cDNA clone. Among the 21 chimeric RNA2 molecules engineered, transcripts of only three of them induced systemic plant infection in the presence of GFLV RNA1. Nematode transmission assays of the three viable recombinant viruses showed that swapping a stretch of (i) 11 residues in the βB-βC loop near the icosahedral 3-fold axis abolished transmission by X. index but was insufficient to restore transmission by X. diversicaudatum and (ii) 7 residues in the βE-αB loop did not interfere with transmission by the two Xiphinema species. This study provides new insights into GFLV CP determinants of nematode transmission.
Diaporthe perjuncta is a pathogen of grapevines worldwide. A positive-strand RNA virus, Diaporthe RNA virus (DaRV), occurs in hypovirulent isolates of this fungus. A virus-free isolate from a South African grapevine was transfected with in vitro-transcribed positive strands of DaRV. Based on reverse transcription-PCR and partial sequence analysis, the transfected virus was identified as DaRV. The in vitro-transcribed RNA transcripts used to transfect fungal spheroplasts contained parts of the vector at their distal ends. These vector sequences were separated from the DaRV genome during replication in the new host. The transfected isolate had morphological features that differed from those of the isogenic virus-free strain, including production of a yellow pigment, a decreased growth rate, and lack of sporulation. An apple-based pathogenicity test did not reveal any differences in virulence between the virus-free and DaRV-transfected isolates. This study showed that virus-free fungal hosts can be successfully transfected with viruses other than the Cryphonectria parasitica hypovirus.
Viroids are circular, highly structured, non-protein-coding RNAs that, usurping cellular enzymes and escaping host defense mechanisms, are able to replicate and move through infected plants. Similarly to viruses, viroid infections are associated with the accumulation of viroid-derived 21–24 nt small RNAs (vd-sRNAs) with the typical features of the small interfering RNAs characteristic of RNA silencing, a sequence-specific mechanism involved in defense against invading nucleic acids and in regulation of gene expression in most eukaryotic organisms.
To gain further insights on the genesis and possible role of vd-sRNAs in plant-viroid interaction, sRNAs isolated from Vitis vinifera infected by Hop stunt viroid (HSVd) and Grapevine yellow speckle viroid 1 (GYSVd1) were sequenced by the high-throughput platform Solexa-Illumina, and the vd-sRNAs were analyzed. The large majority of HSVd- and GYSVd1-sRNAs derived from a few specific regions (hotspots) of the genomic (+) and (−) viroid RNAs, with a prevalence of those from the (−) strands of both viroids. When grouped according to their sizes, vd-sRNAs always assumed a distribution with prominent 21-, 22- and 24-nt peaks, which, interestingly, mapped at the same hotspots.
These findings show that different Dicer-like enzymes (DCLs) target viroid RNAs, preferentially accessing to the same viroid domains. Interestingly, our results also suggest that viroid RNAs may interact with host enzymes involved in the RNA-directed DNA methylation pathway, indicating more complex scenarios than previously thought for both vd-sRNAs genesis and possible interference with host gene expression.
Alignment analysis of the Vv-miRNAs identified from various grapevine cultivars indicates that over 30% orthologous Vv-miRNAs exhibit a 1–3 nucleotide discrepancy only at their ends, suggesting that this sequence discrepancy is not a random event, but might mainly derive from divergence of cultivars. With advantages of miR-RACE technology in determining precise sequences of potential miRNAs from bioinformatics prediction, the precise sequences of vv-miRNAs predicted computationally can be verified with miR-RACE in a different grapevine cultivar. This presents itself as a new approach for large scale discovery of precise miRNAs in different grapevine varieties.
Among 88 unique sequences of Vv-miRNAs from bioinformatics prediction, 83 (96.3%) were successfully validated with MiR-RACE in grapevine cv. ‘Summer Black’. All the validated sequences were identical to their corresponding ones obtained from deep sequencing of the small RNA library of ‘Summer Black’. Quantitative RT-PCR analysis of the expressions levels of 10 Vv-miRNA/target gene pairs in grapevine tissues showed some negative correlation trends. Finally, comparison of Vv-miRNA sequences with their orthologs in Arabidopsis and study on the influence of divergent bases of the orthologous miRNAs on their targeting patterns in grapevine were also done.
The validation of precise sequences of potential Vv-miRNAs from computational prediction in a different grapevine cultivar can be a new way to identify the orthologous Vv-miRNAs. Nucleotide discrepancy of orthologous Vv-miRNAs from different grapevine cultivars normally does not change their target genes. However, sequence variations of some orthologous miRNAs in grapevine and Arabidopsis can change their targeting patterns. These precise Vv-miRNAs sequences validated in our study could benefit some further study on grapevine functional genomics.
Replicating circular RNAs are independent plant pathogens known as viroids, or act to modulate the pathogenesis of plant and animal viruses as their satellite RNAs. The rate of discovery of these subviral pathogens was low over the past 40 years because the classical approaches are technical demanding and time-consuming. We previously described an approach for homology-independent discovery of replicating circular RNAs by analysing the total small RNA populations from samples of diseased tissues with a computational program known as progressive filtering of overlapping small RNAs (PFOR). However, PFOR written in PERL language is extremely slow and is unable to discover those subviral pathogens that do not trigger in vivo accumulation of extensively overlapping small RNAs. Moreover, PFOR is yet to identify a new viroid capable of initiating independent infection. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was 3 to 8 times faster than PFOR. A new computational program was further developed and incorporated into PFOR2 to allow the identification of circular RNAs by deep sequencing of long RNAs instead of small RNAs. PFOR2 analysis of the small RNA libraries from grapevine and apple plants led to the discovery of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA (AHVd-like RNA), respectively. GLVd was proposed as a new species in the genus Apscaviroid, because it contained the typical structural elements found in this group of viroids and initiated independent infection in grapevine seedlings. AHVd-like RNA encoded a biologically active hammerhead ribozyme in both polarities, and was not specifically associated with any of the viruses found in apple plants. We propose that these computational algorithms have the potential to discover novel circular RNAs in plants, invertebrates and vertebrates regardless of whether they replicate and/or induce the in vivo accumulation of small RNAs.
Viroids are a unique class of subviral pathogens found in plants, and they are difficult to identify since they are free circular non-coding RNAs and often replicate to low levels in host cells. We previously described the computational algorithm PFOR that discovers viroids by analyzing total small RNAs of the infected plants obtained by next-generation sequencing platforms. However, the algorithm written in PERL language is very slow, and viroid identification depends on the in vivo accumulation of extensively overlapping sets of small RNAs to target viroids. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was significantly faster than PFOR. We also describe a simple computational program that after incorporation into PFOR2 is capable of identifying viroids from deep sequencing of long RNAs instead of small RNAs. Moreover, we report the identification of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA by the computational approach. Since our new algorithms do not depend on the analysis of viroid-derived small RNAs produced in vivo, it is possible to discover viroids in a wide range of host species including plants, invertebrates and vertebrates.
Plants have followed a reticulate type of evolution and taxa have frequently merged via allopolyploidization. A polyploid structure of sequenced genomes has often been proposed, but the chromosomes belonging to putative component genomes are difficult to identify. The 19 grapevine chromosomes are evolutionary stable structures: their homologous triplets have strongly conserved gene order, interrupted by rare translocations. The aim of this study is to examine how the grapevine nucleotide-binding site (NBS)-encoding resistance (NBS-R) genes have evolved in the genomic context and to understand mechanisms for the genome evolution. We show that, in grapevine, i) helitrons have significantly contributed to transposition of NBS-R genes, and ii) NBS-R gene cluster similarity indicates the existence of two groups of chromosomes (named as Va and Vc) that may have evolved independently. Chromosome triplets consist of two Va and one Vc chromosomes, as expected from the tetraploid and diploid conditions of the two component genomes. The hexaploid state could have been derived from either allopolyploidy or the separation of the Va and Vc component genomes in the same nucleus before fusion, as known for Rosaceae species. Time estimation indicates that grapevine component genomes may have fused about 60 mya, having had at least 40–60 mya to evolve independently. Chromosome number variation in the Vitaceae and related families, and the gap between the time of eudicot radiation and the age of Vitaceae fossils, are accounted for by our hypothesis.
Sequencing of plant nuclear genomes reveals the widespread presence of integrated viral sequences known as endogenous pararetroviruses (EPRVs). Banana is one of the three plant species known to harbor infectious EPRVs. Musa balbisiana carries integrated copies of Banana streak virus (BSV), which are infectious by releasing virions in interspecific hybrids. Here, we analyze the organization of the EPRV of BSV Goldfinger (BSGfV) present in the wild diploid M. balbisiana cv. Pisang Klutuk Wulung (PKW) revealed by the study of Musa bacterial artificial chromosome resources and interspecific genetic cross. cv. PKW contains two similar EPRVs of BSGfV. Genotyping of these integrants and studies of their segregation pattern show an allelic insertion. Despite the fact that integrated BSGfV has undergone extensive rearrangement, both EPRVs contain the full-length viral genome. The high degree of sequence conservation between the integrated and episomal form of the virus indicates a recent integration event; however, only one allele is infectious. Analysis of BSGfV EPRV segregation among an F1 population from an interspecific genetic cross revealed that these EPRV sequences correspond to two alleles originating from a single integration event. We describe here for the first time the full genomic and genetic organization of the two EPRVs of BSGfV present in cv. PKW in response to the challenge facing both scientists and breeders to identify and generate genetic resources free from BSV. We discuss the consequences of this unique host-pathogen interaction in terms of genetic and genomic plant defenses versus strategies of infectious BSGfV EPRVs.
Hepadnaviruses are DNA viruses but, as pararetroviruses, their morphogenesis initiates with the encapsidation of an RNA pregenome, and these viruses have therefore evolved mechanisms to exclude nucleocapsids that contain incompletely matured genomes from participating in budding and secretion. We provide here evidence that binding of hepadnavirus core particles from the cytosol to their target membranes is a distinct step in morphogenesis, discriminating among different populations of intracellular capsids. Using the duck hepatitis B virus (DHBV) and a flotation assay, we found about half of the intracellular capsids to be membrane associated due to an intrinsic membrane-binding affinity. In contrast to free cytosolic capsids, this subpopulation contained largely mature, double-stranded DNA genomes and lacked core protein hyperphosphorylation, both features characteristic for secreted virions. Against expectation, however, the selective membrane attachment observed did not require the presence of the large DHBV envelope protein, which has been considered to be crucial for nucleocapsid-membrane interaction. Furthermore, removal of surface-exposed phosphate residues from nonfloating capsids by itself did not suffice to confer membrane affinity and, finally, hyperphosphorylation was absent from nonenveloped nucleocapsids that were released from DHBV-transfected cells. Collectively, these observations argue for a model in which nucleocapsid maturation, involving the viral genome, capsid structure, and capsid dephosphorylation, leads to the exposure of a membrane-binding signal as a step crucial for selecting the matured nucleocapsid to be incorporated into the capsid-independent budding of virus particles.
A large amount of short interfering RNA (vsiRNA) is generated from plant viruses during infection, but the function, structure and biogenesis of these is not understood. We profiled vsiRNAs using two different high-throughput sequencing platforms and also developed a hybridisation based array approach. The profiles obtained through the Solexa platform and by hybridisation were very similar to each other but different from the 454 profile. Both deep sequencing techniques revealed a strong bias in vsiRNAs for the positive strand of the virus and identified regions on the viral genome that produced vsiRNA in much higher abundance than other regions. The hybridisation approach also showed that the position of highly abundant vsiRNAs was the same in different plant species and in the absence of RDR6. We used the Terminator 5′-Phosphate-Dependent Exonuclease to study the 5′ end of vsiRNAs and showed that a perfect control duplex was not digested by the enzyme without denaturation and that the efficiency of the Terminator was strongly affected by the concentration of the substrate. We found that most vsiRNAs have 5′ monophosphates, which was also confirmed by profiling short RNA libraries following either direct ligation of adapters to the 5′ end of short RNAs or after replacing any potential 5′ ends with monophosphates. The Terminator experiments also showed that vsiRNAs were not perfect duplexes. Using a sensor construct we also found that regions from the viral genome that were complementary to non-abundant vsiRNAs were targeted in planta just as efficiently as regions recognised by abundant vsiRNAs. Different high-throughput sequencing techniques have different reproducible sequence bias and generate different profiles of short RNAs. The Terminator exonuclease does not process double stranded RNA, and because short RNAs can quickly re-anneal at high concentration, this assay can be misleading if the substrate is not denatured and not analysed in a dilution series. The sequence profiles and Terminator digests suggest that CymRSV siRNAs are produced from the structured positive strand rather than from perfect double stranded RNA or by RNA dependent RNA polymerase.
Viral RNA is processed into short RNAs in plants, which guide a complex to the viral RNA and cause cleavage of the viral RNA. We profiled Cymbidium ringspot virus (CymRSV) derived short RNAs using three different methods. Profiling of viral short interfering RNAs revealed a different sequence bias for the 454 and Solexa high-throughput sequencing platforms. We also found that viral short RNAs are primarily produced from the positive strand of the virus and produced with very different frequency along the viral genome. The hybridisation approach showed that the profile of viral short RNAs is determined by the virus itself because the profiles were the same in different species and it also showed that the process was RDR6 independent. We used the Terminator exonuclease to study the 5′ end of viral short RNAs and discovered that this enzyme cannot digest double stranded RNA. A control perfect duplex was only partially processed even after denaturation. Since double stranded short RNAs can quickly re-anneal, this assay must be carried out using different concentrations of the substrate. We found that most of the CymRSV short RNAs had 5′ monophosphate and were not perfect duplexes. Taken together, these results suggest that CymRSV short RNAs are produced from the structured positive strand rather than from perfect double stranded RNA or by RNA dependent RNA polymerase. We also found that regions from the viral genome that are not complementary to highly abundant viral short RNAs were targeted in the plant just as efficiently as regions recognised by abundant short RNAs.
The improvement of the agricultural and wine-making qualities of the grapevine (Vitis vinifera) is hampered by adherence to traditional varieties, the recalcitrance of this plant to genetic modifications, and public resistance to genetically modified organism (GMO) technologies. To address these challenges, we developed an RNA virus-based vector for the introduction of desired traits into grapevine without heritable modifications to the genome. This vector expresses recombinant proteins in the phloem tissue that is involved in sugar transport throughout the plant, from leaves to roots to berries. Furthermore, the vector provides a powerful RNA interference (RNAi) capability of regulating the expression of endogenous genes via virus-induced gene-silencing (VIGS) technology. Additional advantages of this vector include superb genetic capacity and stability, as well as the swiftness of technology implementation. The most significant applications of the viral vector include functional genomics of the grapevine and disease control via RNAi-enabled vaccination against pathogens or invertebrate pests.
Botrytis cinerea CCg378 is a wild-type strain infected with two types of double-stranded RNA (dsRNA) mycoviruses and which presents hypovirulence-associated traits. The objectives of the present study were to characterize the mycoviruses and investigate their relationship with the low virulence degree of the fungal host.
B. cinerea CCg378 contains five dsRNA molecules that are associated with two different types of isometric viral particles of 32 and 23 nm in diameter, formed by structural polypeptides of 70-kDa and 48-kDa, respectively.
The transfection of spheroplasts of a virus-free strain, B. cinerea CKg54, with viral particles purified from the CCg378 strain revealed that the 2.2-kbp dsRNAs have no dependency on the smaller molecules for its stable maintenance in the fungal cytoplasm, because a fungal clone that only contains the 2.2-kbp dsRNAs associated with the 32-nm particles was obtained, which we named B. cinerea CKg54vi378. One of the 2.2 kbpdsRNA segments (2219 bp) was sequenced and corresponds to the gene encoding the capsid protein of B. cinerea CCg378 virus 1 (Bc378V1), a putative new member of the Partitiviridae family.
Furthermore, physiological parameters related to the degree of virulence of the fungus, such as the sporulation rate and laccase activity, were lower in B. cinerea CCg378 and B. cinerea CKg54vi378 than in B. cinerea CKg54. Additionally, bioassays performed on grapevine leaves showed that the CCg378 and CKg54vi378 strains presented a lower degree of invasiveness on the plant tissue than the CKg54 strain.
The results show that B. cinerea CCg378 is coinfected by two mycoviruses and that the 2.2-kbp dsRNAs correspond to the 32-nm mycovirus genome, which would be a new member of the Partitiviridae family as it has the typical pattern of partitiviruses. On the other hand, the results suggest that the hypovirulence of B. cinerea CCg378 could be conferred by both mycoviruses, since the fungal clone B. cinerea CKg54vi378 presents an intermediate virulence between the CKg54 and CCg378 strains. Therefore, the putative partitivirus would be partially contributing to the hypovirulence phenotype of the CCg378 strain.
Mycovirus; Partitivirus; Double-stranded RNA; Botrytis cinerea; Transfection; Hypovirulence
It is known that crustaceans and insects can persistently carry one or more viral pathogens at low levels, without signs of disease. They may transmit them to their offspring or to naïve individuals, often with lethal consequences. The underlying molecular mechanisms have not been elucidated, but the process has been called viral accommodation. Since tolerance to one virus does not confer tolerance to another, tolerance is pathogen-specific, so the requirement for a specific pathogen response mechanism (memory) was included in the original viral accommodation concept. Later, it was hypothesized that specific responses were based on the presence of viruses in persistent infections. However, recent developments suggest that specific responses may be based on viral sequences inserted into the host genome.
Presentation of the hypothesis
Non-retroviral fragments of both RNA and DNA viruses have been found in insect and crustacean genomes. In addition, reverse-transcriptase (RT) and integrase (IN) sequences are also common in their genomes. It is hypothesized that shrimp and other arthropods use these RT to recognize "foreign" mRNA of both RNA and DNA viruses and use the integrases (IN) to randomly insert short cDNA sequences into their genomes. By chance, some of these sequences result in production of immunospecific RNA (imRNA) capable of stimulating RNAi that suppresses viral propagation. Individuals with protective inserts would pass these on to the next generation, together with similar protective inserts for other viruses that could be amalgamated rapidly in individual offspring by random assortment of chromosomes. The most successful individuals would be environmentally selected from billions of offspring.
This hypothesis for immunity based on an imRNA generation mechanism fits with the general principle of invertebrate immunity based on a non-host, "pattern recognition" process. If proven correct, understanding the process would allow directed preparation of vaccines for selection of crustacean and insect lines applicable in commercial production species (e.g., shrimp and bees) or in control of insect-borne diseases. Arising from a natural host mechanism, the resulting animals would not be artificially, genetically modified (GMO).
This article was reviewed by Akria Shibuya, Eugene V. Koonin and L. Aravind.
The version of this article published in BMC Genomics 2009, 10:558, contains data in Table 1 which are now known to be unreliable, and an illustration, in Figure 1, of unusual miRNA processing events predicted by these unreliable data. In this full-length correction, new data replace those found to be unreliable, leading to a more straightforward interpretation without altering the principle conclusions of the study. Table 1 and associated methods have been corrected, Figure 1 deleted, supplementary file 1 added, and modifications made to the sections "Deep sequencing of small RNAs from grapevine leaf tissue" and "Microarray analysis of miRNA expression". The editors and authors regret the inconvenience caused to readers by premature publication of the original paper.
MicroRNAs are short (~21 base) single stranded RNAs that, in plants, are generally coded by specific genes and cleaved specifically from hairpin precursors. MicroRNAs are critical for the regulation of multiple developmental, stress related and other physiological processes in plants. The recent annotation of the genome of the grapevine (Vitis vinifera L.) allowed the identification of many putative conserved microRNA precursors, grouped into multiple gene families.
Here we use oligonucleotide arrays to provide the first indication that many of these microRNAs show differential expression patterns between tissues and during the maturation of fruit in the grapevine. Furthermore we demonstrate that whole transcriptome sequencing and deep-sequencing of small RNA fractions can be used both to identify which microRNA precursors are expressed in different tissues and to estimate genomic coordinates and patterns of splicing and alternative splicing for many primary miRNA transcripts.
Our results show that many microRNAs are differentially expressed in different tissues and during fruit maturation in the grapevine. Furthermore, the demonstration that whole transcriptome sequencing can be used to identify candidate splicing events and approximate primary microRNA transcript coordinates represents a significant step towards the large-scale elucidation of mechanisms regulating the expression of microRNAs at the transcriptional and post-transcriptional levels.
RNA silencing in plants is a natural defense system against foreign genetic elements including viruses. This natural antiviral mechanism has been adopted to develop virus-resistant plants through expression of virus-derived double-stranded RNAs or hairpin RNAs, which in turn are processed into small interfering RNAs (siRNAs) by the host's RNA silencing machinery. While these virus-specific siRNAs were shown to be a hallmark of the acquired virus resistance, the functionality of another set of the RNA silencing-related small RNAs, microRNAs (miRNAs), in engineering plant virus resistance has not been extensively explored. Here we show that expression of an artificial miRNA, targeting sequences encoding the silencing suppressor 2b of Cucumber mosaic virus (CMV), can efficiently inhibit 2b gene expression and protein suppressor function in transient expression assays and confer on transgenic tobacco plants effective resistance to CMV infection. Moreover, the resistance level conferred by the transgenic miRNA is well correlated to the miRNA expression level. Comparison of the anti-CMV effect of the artificial miRNA to that of a short hairpin RNA-derived small RNA targeting the same site revealed that the miRNA approach is superior to the approach using short hairpin RNA both in transient assays and in transgenic plants. Together, our data demonstrate that expression of virus-specific artificial miRNAs is an effective and predictable new approach to engineering resistance to CMV and, possibly, to other plant viruses as well.