Endogenous pararetroviruses—a reservoir of virus infection in plants

Phylogenetic analysis of endogenous viral elements in the rice genome reveals local chromosomal evolution in Oryza AA-genome species

Introduction

Rice genomes contain endogenous viral elements homologous to rice tungro bacilliform virus (RTBV) from the pararetrovirus family Caulimoviridae. These viral elements, known as endogenous RTBV-like sequences (eRTBVLs), comprise five subfamilies, eRTBVL-A, -B, -C, -D, and -X. Four subfamilies (A, B, C, and X) are present to a limited degree in the genomes of the Asian cultivated rice Oryza sativa (spp. japonica and indica) and the closely related wild species Oryza rufipogon.

Methods

The eRTBVL-D sequences are widely distributed within these and other Oryza AA-genome species. Fifteen eRTBVL-D segments identified in the japonica (Nipponbare) genome occur mostly at orthologous chromosomal positions in other AA-genome species. The eRTBVL-D sequences were inserted into the genomes just before speciation of the AA-genome species.

Results and discussion

Ten eRTBVL-D segments are located at six loci, which were used for our evolutionary analyses during the speciation of the AA-genome species. The degree of genetic differentiation varied among the eRTBVL-D segments. Of the six loci, three showed phylogenetic trees consistent with the standard speciation pattern (SSP) of the AA-genome species (Type A), and the other three represented phylogenies different from the SSP (Type B). The atypical phylogenetic trees for the Type B loci revealed chromosome region-specific evolution among the AA-genome species that is associated with phylogenetic incongruences: complex genome rearrangements between eRTBVL-D segments, an introgression between the distant species, and low genetic diversity of a shared eRTBVL-D segment. Using eRTBVL-D as an indicator, this study revealed the phylogenetic incongruence of local chromosomal regions with different topologies that developed during speciation.

Discovery and Analyses of Caulimovirid-like Sequences in Upland Cotton (Gossypium hirsutum)

Analyses of Illumina-based high-throughput sequencing data generated during characterization of the cotton leafroll dwarf virus population in Mississippi (2020-2022) consistently yielded contigs varying in size (most frequently from 4 to 7 kb) with identical nucleotide content and sharing similarities with reverse transcriptases (RTases) encoded by extant plant pararetroviruses (family Caulimoviridiae). Initial data prompted an in-depth study involving molecular and bioinformatic approaches to characterize the nature and origins of these caulimovirid-like sequences. As a result, here, we report on endogenous viral elements (EVEs) related to extant members of the family Caulimoviridae, integrated into a genome of upland cotton (Gossypium hirsutum), for which we propose the provisional name "endogenous cotton pararetroviral elements" (eCPRVE). Our investigations pinpointed a ~15 kbp-long locus on the A04 chromosome consisting of head-to-head orientated tandem copies located on positive- and negative-sense DNA strands (eCPRVE+ and eCPRVE-). Sequences of the eCPRVE+ comprised nearly complete and slightly decayed genome information, including ORFs coding for the viral movement protein (MP), coat protein (CP), RTase, and transactivator/viroplasm protein (TA). Phylogenetic analyses of major viral proteins suggest that the eCPRVE+ may have been initially derived from a genome of a cognate virus belonging to a putative new genus within the family. Unexpectedly, an identical 15 kb-long locus composed of two eCPRVE copies was also detected in a newly recognized species G. ekmanianum, shedding some light on the relatively recent evolution within the cotton family.

Managing the deluge of newly discovered plant viruses and viroids: an optimized scientific and regulatory framework for their characterization and risk analysis

The advances in high-throughput sequencing (HTS) technologies and bioinformatic tools have provided new opportunities for virus and viroid discovery and diagnostics. Hence, new sequences of viral origin are being discovered and published at a previously unseen rate. Therefore, a collective effort was undertaken to write and propose a framework for prioritizing the biological characterization steps needed after discovering a new plant virus to evaluate its impact at different levels. Even though the proposed approach was widely used, a revision of these guidelines was prepared to consider virus discovery and characterization trends and integrate novel approaches and tools recently published or under development. This updated framework is more adapted to the current rate of virus discovery and provides an improved prioritization for filling knowledge and data gaps. It consists of four distinct steps adapted to include a multi-stakeholder feedback loop. Key improvements include better prioritization and organization of the various steps, earlier data sharing among researchers and involved stakeholders, public database screening, and exploitation of genomic information to predict biological properties.

Emerging and Re-Emerging Diseases Caused by Badnaviruses

New and emerging plant diseases are caused by different pathogens including viruses that often cause significant crop losses. Badnaviruses are pararetroviruses that contain a single molecule of ds DNA genome of 7 to 9 kb in size and infect a large number of economically important crops such as banana and plantains, black pepper, cacao, citrus, grapevine, pineapple, sugarcane, sweet potato, taro, and yam, causing significant yield losses. Many of the species in the genus have a restricted host range and several of them are known to infect a single crop. Combined infections of different virus species and strains offer conditions that favor the development of new strains via recombination, especially in vegetatively propagated crops. The primary spread of badnaviruses is through vegetative propagating materials while for the secondary spread, they depend on insects such as mealybugs and aphids. Disease emerges as a consequence of the interactions between host and pathogens under favorable environmental conditions. The viral genome of the pararetroviruses is known to be integrated into the chromosome of the host and a few plants with integrants when subjected to different kinds of abiotic stress will give rise to episomal forms of the virus and cause disease. Attempts have been made to develop management strategies for badnaviruses both conventionally and using precision breeding techniques such as genome editing. Until 2016 only 32 badnavirus species infecting different crops were known, but in a span of six years, this number has gone up to 68. The current review highlights the emerging disease problems and management options for badnaviruses infecting economically important crops.

The essential role of the quasi-long terminal repeat sequence for replication and gene expression of an endogenous pararetrovirus, petunia vein clearing virus

Petunia vein clearing virus (PVCV) is a type member of the genus Petuvirus within the Caulimoviridae family and is defined as one viral unit consisting of a single open reading frame (ORF) encoding a viral polyprotein and one quasi-long terminal repeat (QTR) sequence. Since some full-length PVCV sequences are found in the petunia genome and a vector for horizontal transmission of PVCV has not been identified yet, PVCV is referred to as an endogenous pararetrovirus. Molecular mechanisms of replication, gene expression and horizontal transmission of endogenous pararetroviruses in plants are elusive. In this study, agroinfiltration experiments using various PVCV infectious clones indicated that the replication (episomal DNA synthesis) and gene expression of PVCV were efficient when the QTR sequences are present on both sides of the ORF. Whereas replacement of the QTR with another promoter and/or terminator is possible for gene expression, it is essential for QTR sequences to be on both sides for viral replication. Although horizontal transmission of PVCV by grafting and biolistic inoculation was previously reported, agroinfiltration is a useful and convenient method for studying its replication and gene expression.

Genome-wide identification of Reverse Transcriptase domains of recently inserted endogenous plant pararetrovirus (Caulimoviridae)

Endogenous viral elements (EVEs) are viral sequences that have been integrated into the nuclear chromosomes. Endogenous pararetrovirus (EPRV) are a class of EVEs derived from DNA viruses of the family Caulimoviridae. Previous works based on a limited number of genome assemblies demonstrated that EPRVs are abundant in plants and are present in several species. The availability of genome sequences has been immensely increased in the recent years and we took advantage of these resources to have a more extensive view of the presence of EPRVs in plant genomes. We analyzed 278 genome assemblies corresponding to 267 species (254 from Viridiplantae) using tBLASTn against a collection of conserved domains of the Reverse Transcriptases (RT) of Caulimoviridae. We concentrated our search on complete and well-conserved RT domains with an uninterrupted ORF comprising the genetic information for at least 300 amino acids. We obtained 11.527 sequences from the genomes of 202 species spanning the whole Tracheophyta clade. These elements were grouped in 57 clusters and classified in 13 genera, including a newly proposed genus we called Wendovirus. Wendoviruses are characterized by the presence of four open reading frames and two of them encode for aspartic proteinases. Comparing plant genomes, we observed important differences between the plant families and genera in the number and type of EPRVs found. In general, florendoviruses are the most abundant and widely distributed EPRVs. The presence of multiple identical RT domain sequences in some of the genomes suggests their recent amplification.

Interspecific hybridization in tomato influences endogenous viral sRNAs and alters gene expression

Background

Hybridization is associated with the activation of transposable elements and changes in the patterns of gene expression leading to phenotypic changes. However, the underlying mechanisms are not well understood.

Results

Here, we describe the changes to the gene expression in interspecific Solanum hybrids that are associated with small RNAs derived from endogenous (para)retroviruses (EPRV). There were prominent changes to sRNA profiles in these hybrids involving 22-nt species produced in the DCL2 biogenesis pathway, and the hybridization-induced changes to the gene expression were similar to those in a dcl2 mutant.

Conclusions

These findings indicate that hybridization leads to activation of EPRV, perturbation of small RNA profiles, and, consequently, changes in the gene expression. Such hybridization-induced variation in the gene expression could increase the natural phenotypic variation in natural evolution or in breeding for agriculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02685-z.

CRISPR/Cas-Mediated Resistance against Viruses in Plants

CRISPR/Cas9 provides a robust and widely adaptable system with enormous potential for genome editing directed towards generating useful products. It has been used extensively to generate resistance against viruses infecting plants with more effective and prolonged efficiency as compared with previous antiviral approaches, thus holding promise to alleviate crop losses. In this review, we have discussed the reports of CRISPR/Cas-based virus resistance strategies against plant viruses. These strategies include approaches targeting single or multiple genes (or non-coding region) in the viral genome and targeting host factors essential for virus propagation. In addition, the utilization of base editing has been discussed to generate transgene-free plants resistant to viruses. This review also compares the efficiencies of these approaches. Finally, we discuss combinatorial approaches, including multiplexing, to increase editing efficiency and bypass the generation of escape mutants.

Extrachromosomal viral DNA produced by transcriptionally active endogenous viral elements in non-infected banana hybrids impedes quantitative PCR diagnostics of banana streak virus infections in banana hybrids

The main edible and cultivated banana varieties are intra- and interspecific hybrids of the two main Musa species, Musa acuminata and Musa balbisiana, having diploid genomes denoted A and B, respectively. The B genome naturally hosts sequences of banana streak virus (BSV) named endogenous BSV (eBSV). Upon stress, eBSVs are identified as the origin of BSV infection for at least three BSV species, causing banana streak disease. For each of the three species, BSV and eBSV share >99.9 % sequence identity, complicating PCR-based diagnosis of viral infection in the B genome-containing bananas. Here, we designed a quantitative PCR-based method to only quantify episomal BSV particles produced, overcoming the limitation of eBSV also being detected by qPCR by using it as a 'calibrator'. However, our results revealed unexpected variation of eBSV amplification in calibrator plants composed of a clonal population of 53 replicating virus-free banana hybrids with the same AAB genotype. Our in-depth molecular analyses suggest that this calibrator variation is due to the variable abundance of non-encapsidated extrachromosomal viral DNA, likely produced via the transcription of eBSVs, followed by occasional reverse transcription. We also present evidence that accumulation of viral transcripts in AAB plants is downregulated both at post-transcriptional and transcriptional levels by an RNA interference mechanism that keeps the plants free of virus infection. Finally, we recommend that such eBSV amplification variation be taken into account to establish a quantitative viral diagnostic for banana plants with the B genome.

Under siege: virus control in plant meristems and progeny

In the arms race between plants and viruses, two frontiers have been utilized for decades to combat viral infections in agriculture. First, many pathogenic viruses are excluded from plant meristems, which allows the regeneration of virus-free plant material by tissue culture. Second, vertical transmission of viruses to the host progeny is often inefficient, thereby reducing the danger of viral transmission through seeds. Numerous reports point to the existence of tightly linked meristematic and transgenerational antiviral barriers that remain poorly understood. In this review, we summarize the current understanding of the molecular mechanisms that exclude viruses from plant stem cells and progeny. We also discuss the evidence connecting viral invasion of meristematic cells and the ability of plants to recover from acute infections. Research spanning decades performed on a variety of virus/host combinations has made clear that, beside morphological barriers, RNA interference (RNAi) plays a crucial role in preventing-or allowing-meristem invasion and vertical transmission. How a virus interacts with plant RNAi pathways in the meristem has profound effects on its symptomatology, persistence, replication rates, and, ultimately, entry into the host progeny.

Broken, silent, and in hiding: tamed endogenous pararetroviruses escape elimination from the genome of sugar beet (Beta vulgaris)

BACKGROUND AND AIMS

Endogenous pararetroviruses (EPRVs) are widespread components of plant genomes that originated from episomal DNA viruses of the Caulimoviridae family. Due to fragmentation and rearrangements, most EPRVs have lost their ability to replicate through reverse transcription and to initiate viral infection. Similar to the closely related retrotransposons, extant EPRVs were retained and often amplified in plant genomes for several million years. Here, we characterize the complete genomic EPRV fraction of the crop sugar beet (Beta vulgaris, Amaranthaceae) to understand how they shaped the beet genome and to suggest explanations for their absent virulence.

METHODS

Using next- and third-generation sequencing data and genome assembly, we reconstructed full-length in silico representatives for the three host-specific EPRVs (beetEPRVs) in the B. vulgaris genome. Focusing on the endogenous caulimovirid beetEPRV3, we investigated its chromosomal localization, abundance and distribution by fluorescent in situ and Southern hybridization.

KEY RESULTS

Full-length beetEPRVs range between 7.5 and 10.7 kb in size, are heterogeneous in structure and sequence, and occupy about 0.3 % of the beet genome. Although all three beetEPRVs were assigned to the florendoviruses, they showed variably arranged protein-coding domains, different fragmentation, and preferences for diverse sequence contexts. We observed small RNAs that specifically target the individual beetEPRVs, indicating stringent epigenetic suppression. BeetEPRV3 sequences occur along all sugar beet chromosomes, preferentially in the vicinity of each other and are associated with heterochromatic, centromeric and intercalary satellite DNAs. BeetEPRV3 members also exist in genomes of related wild species, indicating an initial beetEPRV3 integration 13.4-7.2 million years ago.

CONCLUSIONS

Our study in beet illustrates the variability of EPRV structure and sequence in a single host genome. Evidence of sequence fragmentation and epigenetic silencing implies possible plant strategies to cope with long-term persistence of EPRVs, including amplification, fixation in the heterochromatin, and containment of EPRV virulence.

Badnaviruses and banana genomes: a long association sheds light on Musa phylogeny and origin

Badnaviruses are double-stranded DNA pararetroviruses of the family Caulimoviridae. Badnaviral sequences found in banana are distributed over three main clades of the genus Badnavirus and exhibit wide genetic diversity. Interestingly, the nuclear genome of many plants, including banana, is invaded by numerous badnaviral sequences although badnaviruses do not require an integration step to replicate, unlike animal retroviruses. Here, we confirm that banana streak viruses (BSVs) are restricted to clades 1 and 3. We also show that only BSVs from clade 3 encompassing East African viral species are not integrated into Musa genomes, unlike BSVs from clade 1. Finally, we demonstrate that sequences from clade 2 are definitively integrated into Musa genomes with no evidence of episomal counterparts; all are phylogenetically distant from BSVs known to date. Using different molecular approaches, we dissected the coevolution between badnaviral sequences of clade 2 and banana by comparing badnavirus integration patterns across a banana sampling representing major Musa speciation events. Our data suggest that primary viral integrations occurred millions of years ago in banana genomes under different possible scenarios. Endogenous badnaviral sequences can be used as powerful markers to better characterize the Musa phylogeny, narrowing down the likely geographical origin of the Musa ancestor.

Participation of Multifunctional RNA in Replication, Recombination and Regulation of Endogenous Plant Pararetroviruses (EPRVs)

Pararetroviruses, taxon Caulimoviridae, are typical of retroelements with reverse transcriptase and share a common origin with retroviruses and LTR retrotransposons, presumably dating back 1.6 billion years and illustrating the transition from an RNA to a DNA world. After transcription of the viral genome in the host nucleus, viral DNA synthesis occurs in the cytoplasm on the generated terminally redundant RNA including inter- and intra-molecule recombination steps rather than relying on nuclear DNA replication. RNA recombination events between an ancestral genomic retroelement with exogenous RNA viruses were seminal in pararetrovirus evolution resulting in horizontal transmission and episomal replication. Instead of active integration, pararetroviruses use the host DNA repair machinery to prevail in genomes of angiosperms, gymnosperms and ferns. Pararetrovirus integration - leading to Endogenous ParaRetroViruses, EPRVs - by illegitimate recombination can happen if their sequences instead of homologous host genomic sequences on the sister chromatid (during mitosis) or homologous chromosome (during meiosis) are used as template. Multiple layers of RNA interference exist regulating episomal and chromosomal forms of the pararetrovirus. Pararetroviruses have evolved suppressors against this plant defense in the arms race during co-evolution which can result in deregulation of plant genes. Small RNAs serve as signaling molecules for Transcriptional and Post-Transcriptional Gene Silencing (TGS, PTGS) pathways. Different populations of small RNAs comprising 21-24 nt and 18-30 nt in length have been reported for Citrus, Fritillaria, Musa, Petunia, Solanum and Beta. Recombination and RNA interference are driving forces for evolution and regulation of EPRVs.

Deep Roots and Splendid Boughs of the Global Plant Virome

Land plants host a vast and diverse virome that is dominated by RNA viruses, with major additional contributions from reverse-transcribing and single-stranded (ss) DNA viruses. Here, we introduce the recently adopted comprehensive taxonomy of viruses based on phylogenomic analyses, as applied to the plant virome. We further trace the evolutionary ancestry of distinct plant virus lineages to primordial genetic mobile elements. We discuss the growing evidence of the pivotal role of horizontal virus transfer from invertebrates to plants during the terrestrialization of these organisms, which was enabled by the evolution of close ecological associations between these diverse organisms. It is our hope that the emerging big picture of the formation and global architecture of the plant virome will be of broad interest to plant biologists and virologists alike and will stimulate ever deeper inquiry into the fascinating field of virus-plant coevolution.

Disturbance of floral colour pattern by activation of an endogenous pararetrovirus, petunia vein clearing virus, in aged petunia plants

White areas of star-type bicolour petals of petunia (Petunia hybrida) are caused by post-transcriptional gene silencing (PTGS) of the key enzyme of anthocyanin biosynthesis. We observed blotched flowers and a vein-clearing symptom in aged petunia plants. To determine the cause of blotched flowers, we focused on an endogenous pararetrovirus, petunia vein clearing virus (PVCV), because this virus may have a suppressor of PTGS (VSR). Transcripts and episomal DNAs derived from proviral PVCVs accumulated in aged plants, indicating that PVCV was activated as the host plant aged. Furthermore, DNA methylation of CG and CHG sites in the promoter region of proviral PVCV decreased in aged plants, suggesting that poor maintenance of DNA methylation activates PVCV. In parallel, de novo DNA methylation of CHH sites in its promoter region was also detected. Therefore, both activation and inactivation of PVCV occurred in aged plants. The accumulation of PVCV transcripts and episomal DNAs in blotched regions and the detection of VSR activity support a mechanism in which suppression of PTGS by PVCV causes blotched flowers.

Badnaviruses of Sweet Potato: Symptomless Coinhabitants on a Global Scale

Sweet potato is among the most important root crops worldwide, particularly in developing countries, and its production is affected severely by a variety of virus diseases. During the last decade, a number of new viruses have been discovered in sweet potatoes through next-generation sequencing studies. Among them are viruses belonging to the genus Badnavirus and collectively assigned to the species sweet potato pakakuy virus (SPPV). We determined the complete genome sequence of two SPPV isolates and show the ubiquitous presence of similar viruses in germplasm and field material from around the globe. We show that SPPV is not integrated into the sweet potato genome, occurs only at extremely low titers, and is efficiently transmitted through seeds and cuttings. They are unaffected by virus elimination therapy and do not induce discernible symptoms in sweet potatoes or indicator host plants. They show considerable variation in their nucleotide sequences and correspond to several genetic lineages. Studies of their interaction with the two most important sweet potato viruses showed only limited synergistic increase in the titers of one of two SPPV isolates. We contend that these viruses may pose little threat to sweet potato production and more likely represent a new type of persistent virus in sweet potato.

Recombinase polymerase amplification assay for the detection of piper yellow mottle virus infecting black pepper

Recombinase polymerase amplification (RPA) is a quick, specific, sensitive molecular tool carried out at a constant temperature for pathogen detection. In the present study, RPA and reverse transcription (RT) RPA assays were optimized for the detection of piper yellow mottle virus (PYMoV) infecting black pepper. Out of the eight primer pairs targeted to amplify open reading frames (ORFs) 2 and 3 of the virus, the primer pair targeted to ORF2 gave specific amplification only with DNA isolated from infected plant but not with healthy plant. A magnesium acetate concentration of 18 mM, 40 min of incubation time and a temperature of 37-42 °C was found optimum for detection of the virus in RPA assay. Comparison of sensitivity of detection revealed that RPA could detect the virus up to 10^-5 dilution of the total DNA while PCR could detect the virus up to 10^-4 dilution indicating that RPA is 10 times more sensitive than PCR. RPA was further simplified using crude extract as template which could detect the virus up to 10^-3 dilution. RT-RPA was optimized for the detection of PYMoV using total RNA isolated from infected plants as the template. Both RT-RPA and RPA assays were validated using field samples of black pepper representing different varieties and geographical regions by using CTAB isolated DNA, crude DNA extract and cDNA. Our study showed that RPA and RT-RPA can be successfully adopted as a substitute to PCR for detection of PYMoV infecting black pepper.

Virus Latency and the Impact on Plants

Plant viruses are thought to be essentially harmful to the lives of their cultivated crop hosts. In most cases studied, the interaction between viruses and cultivated crop plants negatively affects host morphology and physiology, thereby resulting in disease. Native wild/non-cultivated plants are often latently infected with viruses without any clear symptoms. Although seemingly non-harmful, these viruses pose a threat to cultivated crops because they can be transmitted by vectors and cause disease. Reports are accumulating on infections with latent plant viruses that do not cause disease but rather seem to be beneficial to the lives of wild host plants. In a few cases, viral latency involves the integration of full-length genome copies into the host genome that, in response to environmental stress or during certain developmental stages of host plants, can become activated to generate and replicate episomal copies, a transition from latency to reactivation and causation of disease development. The interaction between viruses and host plants may also lead to the integration of partial-length segments of viral DNA genomes or copy DNA of viral RNA genome sequences into the host genome. Transcripts derived from such integrated viral elements (EVEs) may be beneficial to host plants, for example, by conferring levels of virus resistance and/or causing persistence/latency of viral infections. Studies on viral latency in wild host plants might help us to understand and elucidate the underlying mechanisms of latency and provide insights into the raison d’être for viruses in the lives of plants.

The Previously Unidentified, Divergent Badnavirus Species Cacao red vein-banding virus is Associated with Cacao Swollen Shoot Disease in Nigeria

Cacao swollen shoot disease (CSSD) of Theobroma cacao was reported in Nigeria in 1944; however, no badnaviral genome sequences have been found associated with the symptomatic trees. In 2017, leaf samples (n = 18) were collected from cacao trees from Osun and Oyo, Nigeria showing foliar symptoms that included red vein-banding and shoot swelling, and variable secondary mosaic, mottling, and fern-like pattern symptoms. Abutting primers designed around previously determined 500-bp intergenic region sequences were used for polymerase chain reaction (PCR) amplification. Of the 18 samples, 9 yielded an approximately 7,000-bp, apparently genome-size product. The nine genomes were sequenced and found to encode four open reading frames, and to share 86 to 99% nucleotide identity. Pairwise analysis of the Nigerian genomes with 21 previously reported CSSD badnaviruses, at the complete genome and reverse-transcription ribonuclease H (1,230 bp) sequence levels, indicated 71 to 75 and 72 to 76% nucleotide identity, respectively. Phylogenetic analysis of the nine complete genomes indicated that the closest relatives of the divergent Nigerian isolates were previously described West African CSSD badnaviruses. Based on pairwise comparisons and phylogenetic analyses, the Nigerian CSSD isolates constitute a previously unrecognized Badnavirus sp., herein named Cacao red vein-banding virus (CRVBV). Primers designed based on the CRVBV genome sequences amplified a 1,068-bp fragment from 16 of 18 field samples tested by PCR, suggesting the possible existence of additional CRVBV variants.

Identification and characterization of water chestnut Soymovirus-1 (WCSV-1), a novel Soymovirus in water chestnuts (Eleocharis dulcis)

BACKGROUND

A disease of unknown etiology in water chestnut plants (Eleocharis dulcis) was reported in China between 2012 and 2014. High throughput sequencing of small RNA (sRNA) combined with bioinformatics, and molecular identification based on PCR detection with virus-specific primers and DNA sequencing is a desirable approach to identify an unknown infectious agent. In this study, we employed this approach to identify viral sequences in water chestnut plants and to explore the molecular interaction of the identified viral pathogen and its natural plant host.

RESULTS

Based on high throughput sequencing of virus-derived small RNAs (vsRNA), we identified the sequence a new-to-science double-strand DNA virus isolated from water chestnut cv. 'Tuanfeng' samples, a widely grown cultivar in Hubei province, China, and analyzed its genomic organization. The complete genomic sequence is 7535 base-pairs in length, and shares 42-52% nucleotide sequence identity with viruses in the Caulimoviridae family. The virus contains nine predicated open reading frames (ORFs) encoding nine hypothetical proteins, with conserved domains characteristic of caulimoviruses. Phylogenetic analyses at the nucleotide and amino acid levels indicated that the virus belongs to the genus Soymovirus. The virus is tentatively named Water chestnut soymovirus-1 (WCSV-1). Phylogenetic analysis of the putative viral polymerase protein suggested that WCSV-1 is distinct to other well established species in the Soymovirus genus. This conclusion was supported by phylogenetic analyses of the amino acid sequences encoded by ORFs I, IV, VI, or VII. The sRNA bioinformatics showed that the majority of the vsRNAs are 22-nt in length with a preference for U at the 5'-terminal nucleotide. The vsRNAs are unevenly distributed over both strands of the entire WCSV-1 circular genome, and are clustered into small defined regions. In addition, we detected WCSV-1 in asymptomatic and symptomatic water chestnut samples collected from different regions of China by using PCR. RNA-seq assays further confirmed the presence of WCSV-1-derived viral RNA in infected plants.

CONCLUSIONS

This is the first discovery of a dsDNA virus in the genus Soymovirus infecting water chestnuts. Data presented also add new information towards a better understanding of the co-evolutionary mechanisms between the virus and its natural plant host.

Small RNA-Omics for Plant Virus Identification, Virome Reconstruction, and Antiviral Defense Characterization

RNA interference (RNAi)-based antiviral defense generates small interfering RNAs that represent the entire genome sequences of both RNA and DNA viruses as well as viroids and viral satellites. Therefore, deep sequencing and bioinformatics analysis of small RNA population (small RNA-ome) allows not only for universal virus detection and genome reconstruction but also for complete virome reconstruction in mixed infections. Viral infections (like other stress factors) can also perturb the RNAi and gene silencing pathways regulating endogenous gene expression and repressing transposons and host genome-integrated endogenous viral elements which can potentially be released from the genome and contribute to disease. This review describes the application of small RNA-omics for virus detection, virome reconstruction and antiviral defense characterization in cultivated and non-cultivated plants. Reviewing available evidence from a large and ever growing number of studies of naturally or experimentally infected hosts revealed that all families of land plant viruses, their satellites and viroids spawn characteristic small RNAs which can be assembled into contigs of sufficient length for virus, satellite or viroid identification and for exhaustive reconstruction of complex viromes. Moreover, the small RNA size, polarity and hotspot profiles reflect virome interactions with the plant RNAi machinery and allow to distinguish between silent endogenous viral elements and their replicating episomal counterparts. Models for the biogenesis and functions of small interfering RNAs derived from all types of RNA and DNA viruses, satellites and viroids as well as endogenous viral elements are presented and discussed.

Identification of viral domains integrated into Arabidopsis proteome

Horizontal gene transfer (HGT) contributes to the genome evolution of living organisms. In particular, several recent studies provide convincing data on the integration of viral sequences into diverse organisms. Here, we identified 101 viral domains integrated into the model plant Arabidopsis proteome. Functional analysis based on gene ontology (GO) terms indicates that viral domains in the Arabidopsis proteome were involved in various stress responses with binding functions. Protein interaction networks support the strong protein interactions of viral domains with other Arabidopsis proteins. A proteome-wide analysis gave a comprehensive evolutionary view of viral domains integrated into 41 plant proteomes, revealing the specific and conserved integration of viral domains into plant proteomes. Phylogenetic analyses revealed the possible HGT between viral domains and plant proteomes. Our results provide an overview of the integration of viral domains into plant proteomes and their possible functional roles associated with plant defense mechanisms.

Endogenous pararetrovirus sequences are widely present in Citrinae genomes

Endogenous pararetroviruses (EPRVs) are characterized in several plant genomes and their biological effects have been reported. In this study, hundreds of EPRV segments were identified in six Citrinae genomes. A total of 1034 EPRV segments were identified in the genomes of sweet orange, 2036 in pummelo, 598 in clementine mandarin, 752 in Ichang papeda, 2060 in citron and 245 in atalantia. Genomic analysis indicated that EPRV segments tend to cluster as hot spots in the genomes, particularly on chromosome 2 and 5. Large numbers of simple repeats and transposable elements were identified in the 2-kb flanking regions of the EPRV segments. Comparative genomic analysis and PCR experiments showed that there are highly conserved EPRV segments and species-specific EPRV segments between the Citrinae genomes. Phylogenetic analysis suggested that the integration events of EPRVs could initiate in a common progenitor of Citrinae species and repeatedly occur during the Citrinae divergence.

Ribosome Shunting, Polycistronic Translation, and Evasion of Antiviral Defenses in Plant Pararetroviruses and Beyond

Viruses have compact genomes and usually translate more than one protein from polycistronic RNAs using leaky scanning, frameshifting, stop codon suppression or reinitiation mechanisms. Viral (pre-)genomic RNAs often contain long 5′-leader sequences with short upstream open reading frames (uORFs) and secondary structure elements, which control both translation initiation and replication. In plants, viral RNA and DNA are targeted by RNA interference (RNAi) generating small RNAs that silence viral gene expression, while viral proteins are recognized by innate immunity and autophagy that restrict viral infection. In this review we focus on plant pararetroviruses of the family Caulimoviridae and describe the mechanisms of uORF- and secondary structure-driven ribosome shunting, leaky scanning and reinitiation after translation of short and long uORFs. We discuss conservation of these mechanisms in different genera of Caulimoviridae, including host genome-integrated endogenous viral elements, as well as in other viral families, and highlight a multipurpose use of the highly-structured leader sequence of plant pararetroviruses in regulation of translation, splicing, packaging, and reverse transcription of pregenomic RNA (pgRNA), and in evasion of RNAi. Furthermore, we illustrate how targeting of several host factors by a pararetroviral effector protein can lead to transactivation of viral polycistronic translation and concomitant suppression of antiviral defenses. Thus, activation of the plant protein kinase target of rapamycin (TOR) by the Cauliflower mosaic virus transactivator/viroplasmin (TAV) promotes reinitiation of translation after long ORFs on viral pgRNA and blocks antiviral autophagy and innate immunity responses, while interaction of TAV with the plant RNAi machinery interferes with antiviral silencing.

Utilization of engineered resistance to viruses in crops of the developing world, with emphasis on sub-Saharan Africa

Viral diseases in crop plants constitute a major obstacle to food security in the developing world. Subsistence crops, including cassava, sweetpotato, potato, banana, papaya, common bean, rice and maize are often infected with RNA and/or DNA viruses that cannot be controlled with pesticides. Hence, healthy planting materials and virus-resistant cultivars are essential for high yields of good quality. However, resistance genes are not available for all viral diseases of crop plants. Therefore, virus resistance engineered in plants using modern biotechnology methods is an important addition to the crop production toolbox.

Symbiosis: Viruses as Intimate Partners

Viruses must establish an intimate relationship with their hosts and vectors in order to infect, replicate, and disseminate; hence, viruses can be considered as symbionts with their hosts. Symbiotic relationships encompass different lifestyles, including antagonistic (or pathogenic, the most well-studied lifestyle for viruses), commensal (probably the most common lifestyle), and mutualistic (important beneficial partners). Symbiotic relationships can shape the evolution of the partners in a holobiont, and placing viruses in this context provides an important framework for understanding virus-host relationships and virus ecology. Although antagonistic relationships are thought to lead to coevolution, this is not always clear in virus-host interactions, and impacts on evolution may be complex. Commensalism implies a hitchhiking role for viruses-selfish elements just along for the ride. Mutualistic relationships have been described in detail in the past decade, and they reveal how important viruses are in considering host ecology. Ultimately, symbiosis can lead to symbiogenesis, or speciation through fusion, and the presence of large amounts of viral sequence in the genomes of everything from bacteria to humans, including some important functional genes, illustrates the significance of viral symbiogenesis in the evolution of all life on Earth.

PCR-DGGE Analysis: Unravelling Complex Mixtures of Badnavirus Sequences Present in Yam Germplasm

Badnaviruses (family Caulimoviridae, genus Badnavirus) have emerged as serious pathogens especially affecting the cultivation of tropical crops. Badnavirus sequences can be integrated in host genomes, complicating the detection of episomal infections and the assessment of viral genetic diversity in samples containing a complex mixture of sequences. Yam (Dioscorea spp.) plants are hosts to a diverse range of badnavirus species, and recent findings have suggested that mixed infections occur frequently in West African yam germplasm. Historically, the determination of the diversity of badnaviruses present in yam breeding lines has been achieved by cloning and sequencing of polymerase chain reaction (PCR) products. In this study, the molecular diversity of partial reverse transcriptase (RT)-ribonuclease H (RNaseH) sequences from yam badnaviruses was analysed using PCR-dependent denaturing gradient gel electrophoresis (PCR-DGGE). This resulted in the identification of complex ‘fingerprints’ composed of multiple sequences of Dioscorea bacilliform viruses (DBVs). Many of these sequences show high nucleotide identities to endogenous DBV (eDBV) sequences deposited in GenBank, and fall into six monophyletic species groups. Our findings highlight PCR-DGGE as a powerful tool in badnavirus diversity studies enabling a rapid indication of sequence diversity as well as potential candidate integrated sequences revealed by their conserved nature across germplasm.

Genomic fossils reveal adaptation of non-autonomous pararetroviruses driven by concerted evolution of noncoding regulatory sequences

The interplay of different virus species in a host cell after infection can affect the adaptation of each virus. Endogenous viral elements, such as endogenous pararetroviruses (PRVs), have arisen from vertical inheritance of viral sequences integrated into host germline genomes. As viral genomic fossils, these sequences can thus serve as valuable paleogenomic data to study the long-term evolutionary dynamics of virus-virus interactions, but they have rarely been applied for this purpose. All extant PRVs have been considered autonomous species in their parasitic life cycle in host cells. Here, we provide evidence for multiple non-autonomous PRV species with structural defects in viral activity that have frequently infected ancient grass hosts and adapted through interplay between viruses. Our paleogenomic analyses using endogenous PRVs in grass genomes revealed that these non-autonomous PRV species have participated in interplay with autonomous PRVs in a possible commensal partnership, or, alternatively, with one another in a possible mutualistic partnership. These partnerships, which have been established by the sharing of noncoding regulatory sequences (NRSs) in intergenic regions between two partner viruses, have been further maintained and altered by the sequence homogenization of NRSs between partners. Strikingly, we found that frequent region-specific recombination, rather than mutation selection, is the main causative mechanism of NRS homogenization. Our results, obtained from ancient DNA records of viruses, suggest that adaptation of PRVs has occurred by concerted evolution of NRSs between different virus species in the same host. Our findings further imply that evaluation of within-host NRS interactions within and between populations of viral pathogens may be important.

The Search for Resistance to Cassava Mosaic Geminiviruses: How Much We Have Accomplished, and What Lies Ahead

The cassava mosaic disease (CMD), which occurs in all cassava growing regions of Africa and the Indian subcontinent, is caused by cassava mosaic geminiviruses (CMGs). CMGs are considered to be the most damaging vector-borne plant pathogens. So far, the most successful approach used to control these viruses has been the transfer of a polygenic recessive resistance locus, designated CMD1, from wild cassava to cassava cultivars. Further progress in harnessing natural resistance to contain CMGs has come from the discovery of the dominant monogenic resistance locus, CMD2, in some West African cassava cultivars. CMD2 has been combined with CMD1 through genetic crosses. Because of the limitations of the cassava breeding approach, especially with regard to time required to produce a variety and the loss of preferred agronomic attributes, efforts have been directed toward the deployment of genetic engineering approaches. Most of these approaches have been centered on RNA silencing strategies, developed mainly in the model plant Nicotiana benthamiana. Early RNA silencing platforms assessed for CMG resistance have been use of viral genes for co-suppression, antisense suppression or for hairpin RNAs-mediated gene silencing. Here, progress and challenges in the deployment of these approaches in the control of CMGs are discussed. Novel functional genomics approaches with potential to overcome some of the drawbacks of the current strategies are also discussed.

Genome-wide analysis of EgEVE_1, a transcriptionally active endogenous viral element associated to small RNAs in Eucalyptus genomes

Endogenous viral elements (EVEs) are the result of heritable horizontal gene transfer from viruses to hosts. In the last years, several EVE integration events were reported in plants by the exponential availability of sequenced genomes. Eucalyptus grandis is a forest tree species with a sequenced genome that is poorly studied in terms of evolution and mobile genetic elements composition. Here we report the characterization of E. grandis endogenous viral element 1 (EgEVE_1), a transcriptionally active EVE with a size of 5,664 bp. Phylogenetic analysis and genomic distribution demonstrated that EgEVE_1 is a newly described member of the Caulimoviridae family, distinct from the recently characterized plant Florendoviruses. Genomic distribution of EgEVE_1 and Florendovirus is also distinct. EgEVE_1 qPCR quantification in Eucalyptus urophylla suggests that this genome has more EgEVE_1 copies than E. grandis. EgEVE_1 transcriptional activity was demonstrated by RT-qPCR in five Eucalyptus species and one intrageneric hybrid. We also identified that Eucalyptus EVEs can generate small RNAs (sRNAs),that might be involved in de novo DNA methylation and virus resistance. Our data suggest that EVE families in Eucalyptus have distinct properties, and we provide the first comparative analysis of EVEs in Eucalyptus genomes.

An internet-based bioinformatics toolkit for plant biosecurity diagnosis and surveillance of viruses and viroids

Background

Detection and preventing entry of exotic viruses and viroids at the border is critical for protecting plant industries trade worldwide. Existing post entry quarantine screening protocols rely on time-consuming biological indicators and/or molecular assays that require knowledge of infecting viral pathogens. Plants have developed the ability to recognise and respond to viral infections through Dicer-like enzymes that cleave viral sequences into specific small RNA products. Many studies reported the use of a broad range of small RNAs encompassing the product sizes of several Dicer enzymes involved in distinct biological pathways. Here we optimise the assembly of viral sequences by using specific small RNA subsets.

Results

We sequenced the small RNA fractions of 21 plants held at quarantine glasshouse facilities in Australia and New Zealand. Benchmarking of several de novo assembler tools yielded SPAdes using a kmer of 19 to produce the best assembly outcomes. We also found that de novo assembly using 21–25 nt small RNAs can result in chimeric assemblies of viral sequences and plant host sequences. Such non-specific assemblies can be resolved by using 21–22 nt or 24 nt small RNAs subsets. Among the 21 selected samples, we identified contigs with sequence similarity to 18 viruses and 3 viroids in 13 samples. Most of the viruses were assembled using only 21–22 nt long virus-derived siRNAs (viRNAs), except for one Citrus endogenous pararetrovirus that was more efficiently assembled using 24 nt long viRNAs. All three viroids found in this study were fully assembled using either 21–22 nt or 24 nt viRNAs. Optimised analysis workflows were customised within the Yabi web-based analytical environment. We present a fully automated viral surveillance and diagnosis web-based bioinformatics toolkit that provides a flexible, user-friendly, robust and scalable interface for the discovery and diagnosis of viral pathogens.

Conclusions

We have implemented an automated viral surveillance and diagnosis (VSD) bioinformatics toolkit that produces improved viruses and viroid sequence assemblies. The VSD toolkit provides several optimised and reusable workflows applicable to distinct viral pathogens. We envisage that this resource will facilitate the surveillance and diagnosis viral pathogens in plants, insects and invertebrates.

Electronic supplementary material

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

Viral reverse transcriptases

Reverse transcriptases (RTs) play a major role in the replication of Retroviridae, Metaviridae, Pseudoviridae, Hepadnaviridae and Caulimoviridae. RTs are enzymes that are able to synthesize DNA using RNA or DNA as templates (DNA polymerase activity), and degrade RNA when forming RNA/DNA hybrids (ribonuclease H activity). In retroviruses and LTR retrotransposons (Metaviridae and Pseudoviridae), the coordinated action of both enzymatic activities converts single-stranded RNA into a double-stranded DNA that is flanked by identical sequences known as long terminal repeats (LTRs). RTs of retroviruses and LTR retrotransposons are active as monomers (e.g. murine leukemia virus RT), homodimers (e.g. Ty3 RT) or heterodimers (e.g. human immunodeficiency virus type 1 (HIV-1) RT). RTs lack proofreading activity and display high intrinsic error rates. Besides, high recombination rates observed in retroviruses are promoted by poor processivity that causes template switching, a hallmark of reverse transcription. HIV-1 RT inhibitors acting on its polymerase activity constitute the backbone of current antiretroviral therapies, although novel drugs, including ribonuclease H inhibitors, are still necessary to fight HIV infections. In Hepadnaviridae and Caulimoviridae, reverse transcription leads to the formation of nicked circular DNAs that will be converted into episomal DNA in the host cell nucleus. Structural and biochemical information on their polymerases is limited, although several drugs inhibiting HIV-1 RT are known to be effective against the human hepatitis B virus polymerase. In this review, we summarize current knowledge on reverse transcription in the five virus families and discuss available biochemical and structural information on RTs, including their biosynthesis, enzymatic activities, and potential inhibition.

Changes in the composition of the RNA virome mark evolutionary transitions in green plants

BACKGROUND

The known plant viruses mostly infect angiosperm hosts and have RNA or small DNA genomes. The only other lineage of green plants with a relatively well-studied virome, unicellular chlorophyte algae, is mostly infected by viruses with large DNA genomes. Thus RNA viruses and small DNA viruses seem to completely displace large DNA virus genomes in late branching angiosperms. To understand better the expansion of RNA viruses in the taxonomic span between algae and angiosperms, we analyzed the transcriptomes of 66 non-angiosperm plants characterized by the 1000 Plants Genomes Project.

RESULTS

We found homologs of virus RNA-dependent RNA polymerases in 28 non-angiosperm plant species, including algae, mosses, liverworts (Marchantiophyta), hornworts (Anthocerotophyta), lycophytes, a horsetail Equisetum, and gymnosperms. Polymerase genes in algae were most closely related to homologs from double-stranded RNA viruses leading latent or persistent lifestyles. Land plants, in addition, contained polymerases close to the homologs from single-stranded RNA viruses of angiosperms, capable of productive infection and systemic spread. For several polymerases, a cognate capsid protein was found in the same library. Another virus hallmark gene family, encoding the 30 K movement proteins, was found in lycophytes and monilophytes but not in mosses or algae.

CONCLUSIONS

The broadened repertoire of RNA viruses suggests that colonization of land and growth in anatomical complexity in land plants coincided with the acquisition of novel sets of viruses with different strategies of infection and reproduction.

Badnaviruses: The Current Global Scenario

Badnaviruses (Family: Caulimoviridae; Genus: Badnavirus) are non-enveloped bacilliform DNA viruses with a monopartite genome containing about 7.2 to 9.2 kb of dsDNA with three to seven open reading frames. They are transmitted by mealybugs and a few species by aphids in a semi-persistent manner. They are one of the most important plant virus groups and have emerged as serious pathogens affecting the cultivation of several horticultural crops in the tropics, especially banana, black pepper, cocoa, citrus, sugarcane, taro, and yam. Some badnaviruses are also known as endogenous viruses integrated into their host genomes and a few such endogenous viruses can be awakened, e.g., through abiotic stress, giving rise to infective episomal forms. The presence of endogenous badnaviruses poses a new challenge for the fool-proof diagnosis, taxonomy, and management of the diseases. The present review aims to highlight emerging disease problems, virus characteristics, transmission, and diagnosis of badnaviruses.

How endogenous plant pararetroviruses shed light on Musa evolution

BACKGROUND AND AIMS

Banana genomes harbour numerous copies of viral sequences derived from banana streak viruses (BSVs) - dsDNA viruses belonging to the family Caulimoviridae.These viral integrants (eBSVs) are mostly defective, probably as a result of 'pseudogenization' driven by host genome evolution. However, some can give rise to infection by releasing a functional viral genome following abiotic stresses. These distinct infective eBSVs correspond to the three main widespread BSV species (BSOLV, BSGFV and BSIMV), fully described within the Musa balbisiana B genomes of the seedy diploid 'Pisang Klutuk Wulung' (PKW).

METHODS

We characterize eBSV distribution among a Musa sampling including seedy BB diploids and interspecific hybrids with Musa acuminate exhibiting different levels of ploidy for the B genome (ABB, AAB, AB). We used representative samples of the two areas of sympatry between M. acuminate and M. balbisiana species representing the native area of the most widely cultivated AAB cultivars (in India and in East Asia, ranging from the Philippines to New Guinea). Seventy-seven accessions were characterized using eBSV-related PCR markers and Southern hybridization approaches. We coded both sets of results to create a common dissimilarity matrix with which to interpret eBSV distribution.

KEY RESULTS

We propose a Musa phylogeny driven by the M. balbisiana genome based on a dendrogram resulting from a joint neighbour-joining analysis of the three BSV species, showing for the first time lineages between BB and ABB/AAB hybrids. eBSVs appear to be relevant phylogenetic markers that can illustrate theM. balbisiana phylogeography story.

CONCLUSION

The theoretical implications of this study for further elucidation of the historical and geographical process of Musa domestication are numerous. Discovery of banana plants with B genome non-infective for eBSV opens the way to the introduction of new genitors in programmes of genetic banana improvement.

Deep sequencing reveals the complete genome and evidence for transcriptional activity of the first virus-like sequences identified in Aristotelia chilensis (Maqui Berry)

Here, we report the genome sequence and evidence for transcriptional activity of a virus-like element in the native Chilean berry tree Aristotelia chilensis. We propose to name the endogenous sequence as Aristotelia chilensis Virus 1 (AcV1). High-throughput sequencing of the genome of this tree uncovered an endogenous viral element, with a size of 7122 bp, corresponding to the complete genome of AcV1. Its sequence contains three open reading frames (ORFs): ORFs 1 and 2 shares 66%–73% amino acid similarity with members of the Caulimoviridae virus family, especially the Petunia vein clearing virus (PVCV), Petuvirus genus. ORF1 encodes a movement protein (MP); ORF2 a Reverse Transcriptase (RT) and a Ribonuclease H (RNase H) domain; and ORF3 showed no amino acid sequence similarity with any other known virus proteins. Analogous to other known endogenous pararetrovirus sequences (EPRVs), AcV1 is integrated in the genome of Maqui Berry and showed low viral transcriptional activity, which was detected by deep sequencing technology (DNA and RNA-seq). Phylogenetic analysis of AcV1 and other pararetroviruses revealed a closer resemblance with Petuvirus. Overall, our data suggests that AcV1 could be a new member of Caulimoviridae family, genus Petuvirus, and the first evidence of this kind of virus in a fruit plant.

Biology, etiology, and control of virus diseases of banana and plantain

Banana and plantain (Musa spp.), produced in 10.3 million ha in the tropics, are among the world's top 10 food crops. They are vegetatively propagated using suckers or tissue culture plants and grown almost as perennial plantations. These are prone to the accumulation of pests and pathogens, especially viruses which contribute to yield reduction and are also barriers to the international exchange of germplasm. The most economically important viruses of banana and plantain are Banana bunchy top virus (BBTV), a complex of banana streak viruses (BSVs) and Banana bract mosaic virus (BBrMV). BBTV is known to cause the most serious economic losses in the "Old World," contributing to a yield reduction of up to 100% and responsible for a dramatic reduction in cropping area. The BSVs exist as episomal and endogenous forms are known to be worldwide in distribution. In India and the Philippines, BBrMV is known to be economically important but recently the virus was discovered in Colombia and Costa Rica, thus signaling its spread into the "New World." Banana and plantain are also known to be susceptible to five other viruses of minor significance, such as Abaca mosaic virus, Abaca bunchy top virus, Banana mild mosaic virus, Banana virus X, and Cucumber mosaic virus. Studies over the past 100 years have contributed to important knowledge on disease biology, distribution, and spread. Research during the last 25 years have led to a better understanding of the virus-vector-host interactions, virus diversity, disease etiology, and epidemiology. In addition, new diagnostic tools were developed which were used for surveillance and the certification of planting material. Due to a lack of durable host resistance in the Musa spp., phytosanitary measures and the use of virus-free planting material are the major methods of virus control. The state of knowledge on BBTV, BBrMV, and BSVs, and other minor viruses, disease spread, and control are summarized in this review.

Endogenous pararetrovirus sequences associated with 24 nt small RNAs at the centromeres of Fritillaria imperialis L. (Liliaceae), a species with a giant genome

Endogenous pararetroviral sequences are the most commonly found virus sequences integrated into angiosperm genomes. We describe an endogenous pararetrovirus (EPRV) repeat in Fritillaria imperialis, a species that is under study as a result of its exceptionally large genome (1C = 42 096 Mbp, approximately 240 times bigger than Arabidopsis thaliana). The repeat (FriEPRV) was identified from Illumina reads using the RepeatExplorer pipeline, and exists in a complex genomic organization at the centromere of most, or all, chromosomes. The repeat was reconstructed into three consensus sequences that formed three interconnected loops, one of which carries sequence motifs expected of an EPRV (including the gag and pol domains). FriEPRV shows sequence similarity to members of the Caulimoviridae pararetrovirus family, with phylogenetic analysis indicating a close relationship to Petuvirus. It is possible that no complete EPRV sequence exists, although our data suggest an abundance that exceeds the genome size of Arabidopsis. Analysis of single nucleotide polymorphisms revealed elevated levels of C→T and G→A transitions, consistent with deamination of methylated cytosine. Bisulphite sequencing revealed high levels of methylation at CG and CHG motifs (up to 100%), and 15-20% methylation, on average, at CHH motifs. FriEPRV's centromeric location may suggest targeted insertion, perhaps associated with meiotic drive. We observed an abundance of 24 nt small RNAs that specifically target FriEPRV, potentially providing a signature of RNA-dependent DNA methylation. Such signatures of epigenetic regulation suggest that the huge genome of F. imperialis has not arisen as a consequence of a catastrophic breakdown in the regulation of repeat amplification.

The genome of African yam (Dioscorea cayenensis-rotundata complex) hosts endogenous sequences from four distinct Badnavirus species

Several endogenous viral elements (EVEs) have been identified in plant genomes, including endogenous pararetroviruses (EPRVs). Here, we report the first characterization of EPRV sequences in the genome of African yam of the Dioscorea cayenensis-rotundata complex. We propose that these sequences should be termed 'endogenous Dioscorea bacilliform viruses' (eDBVs). Molecular characterization of eDBVs shows that they constitute sequences originating from various parts of badnavirus genomes, resulting in a mosaic structure that is typical of most EPRVs characterized to date. Using complementary molecular approaches, we show that eDBVs belong to at least four distinct Badnavirus species, indicating multiple, independent, endogenization events. Phylogenetic analyses of eDBVs support and enrich the current taxonomy of yam badnaviruses and lead to the characterization of a new Badnavirus species in yam. The impact of eDBVs on diagnosis, yam germplasm conservation and movement, and breeding is discussed.

Evasion of short interfering RNA-directed antiviral silencing in Musa acuminata persistently infected with six distinct banana streak pararetroviruses

Vegetatively propagated crop plants often suffer from infections with persistent RNA and DNA viruses. Such viruses appear to evade the plant defenses that normally restrict viral replication and spread. The major antiviral defense mechanism is based on RNA silencing generating viral short interfering RNAs (siRNAs) that can potentially repress viral genes posttranscriptionally through RNA cleavage and transcriptionally through DNA cytosine methylation. Here we examined the RNA silencing machinery of banana plants persistently infected with six pararetroviruses after many years of vegetative propagation. Using deep sequencing, we reconstructed consensus master genomes of the viruses and characterized virus-derived and endogenous small RNAs. Consistent with the presence of endogenous siRNAs that can potentially establish and maintain DNA methylation, the banana genomic DNA was extensively methylated in both healthy and virus-infected plants. A novel class of abundant 20-nucleotide (nt) endogenous small RNAs with 5'-terminal guanosine was identified. In all virus-infected plants, 21- to 24-nt viral siRNAs accumulated at relatively high levels (up to 22% of the total small RNA population) and covered the entire circular viral DNA genomes in both orientations. The hotspots of 21-nt and 22-nt siRNAs occurred within open reading frame (ORF) I and II and the 5' portion of ORF III, while 24-nt siRNAs were more evenly distributed along the viral genome. Despite the presence of abundant viral siRNAs of different size classes, the viral DNA was largely free of cytosine methylation. Thus, the virus is able to evade siRNA-directed DNA methylation and thereby avoid transcriptional silencing. This evasion of silencing likely contributes to the persistence of pararetroviruses in banana plants.

IMPORTANCE

We report that DNA pararetroviruses in Musa acuminata banana plants are able to evade DNA cytosine methylation and transcriptional gene silencing, despite being targeted by the host silencing machinery generating abundant 21- to 24-nucleotide short interfering RNAs. At the same time, the banana genomic DNA is extensively methylated in both healthy and virus-infected plants. Our findings shed light on the siRNA-generating gene silencing machinery of banana and provide a possible explanation why episomal pararetroviruses can persist in plants whereas true retroviruses with an obligatory genome-integration step in their replication cycle do not exist in plants.

Distinct and concurrent pathways of Pol II- and Pol IV-dependent siRNA biogenesis at a repetitive trans-silencer locus in Arabidopsis thaliana

Short interfering RNAs (siRNAs) homologous to transcriptional regulatory regions can induce RNA-directed DNA methylation (RdDM) and transcriptional gene silencing (TGS) of target genes. In our system, siRNAs are produced by transcribing an inverted DNA repeat (IR) of enhancer sequences, yielding a hairpin RNA that is processed by several Dicer activities into siRNAs of 21-24 nt. Primarily 24-nt siRNAs trigger RdDM of the target enhancer in trans and TGS of a downstream GFP reporter gene. We analyzed siRNA accumulation from two different structural forms of a trans-silencer locus in which tandem repeats are embedded in the enhancer IR and distinguished distinct RNA polymerase II (Pol II)- and Pol IV-dependent pathways of siRNA biogenesis. At the original silencer locus, Pol-II transcription of the IR from a 35S promoter produces a hairpin RNA that is diced into abundant siRNAs of 21-24 nt. A silencer variant lacking the 35S promoter revealed a normally masked Pol IV-dependent pathway that produces low levels of 24-nt siRNAs from the tandem repeats. Both pathways operate concurrently at the original silencer locus. siRNAs accrue only from specific regions of the enhancer and embedded tandem repeat. Analysis of these sequences and endogenous tandem repeats producing siRNAs revealed the preferential accumulation of siRNAs at GC-rich regions containing methylated CG dinucleotides. In addition to supporting a correlation between base composition, DNA methylation and siRNA accumulation, our results highlight the complexity of siRNA biogenesis at repetitive loci and show that Pol II and Pol IV use different promoters to transcribe the same template.

The prevalence of badnaviruses in West African yams (Dioscorea cayenensis-rotundata) and evidence of endogenous pararetrovirus sequences in their genomes

Yam (Dioscorea spp.) is an important vegetatively-propagated staple crop in West Africa. Viruses are pervasive in yam worldwide, decreasing growth and yield, as well as hindering the international movement of germplasm. Badnaviruses have been reported to be the most prevalent in yam, and genomes of some other badnaviruses are known to be integrated in their host plant species. However, it was not clear if a similar scenario occurs in Dioscorea yam. This study was conducted to verify the prevalence of badnaviruses, and determine if badnavirus genomes are integrated in the yam genome. Leaf samples (n=58) representing eight species of yam from global yam collections kept at CIRAD, France, and 127 samples of D. rotundata breeding lines (n=112) and landraces (n=15) at IITA, Nigeria, were screened using generic badnavirus PCR primers. Positive amplification of an expected ca. 579bp fragment, corresponding to a partial RT-RNaseH region, was detected in 47 (81%) of 58 samples analysed from CIRAD collections, and 100% of the 127 IITA D. rotundata samples. All the D. cayenensis and D. rotundata samples from the CIRAD and IITA collections tested PCR-positive, and sequencing of a selection of the PCR products confirmed they were typical of the genus Badnavirus. A comparison of serological and nucleic acid techniques was used to investigate whether the PCR-positives were sequences amplified from badnavirus particles or putative endogenous badnavirus sequences in the yam genome. Protein A sandwich-enzyme-linked immunosorbent assay (PAS-ELISA) with badnavirus polyclonal antisera detected cross-reacting viral particles in only 60% (92 of 153) of the CIRAD collection samples analysed, in contrast to the aforementioned 81% by PCR. Immunosorbent electron microscopy (ISEM) of virus preparations of a select set of 16 samples, representing different combinations of positive and negative PCR and PAS-ELISA results, identified bacilliform particles in 11 of these samples. Three PCR-positive yam samples from Burkina Faso (cv. Pilimpikou) were identified in which no viral particles were detected by either PAS-ELISA or ISEM. Southern hybridisation results using a yam badnavirus RT-RNaseH sequence (Gn155Dr) as probe, supported a lack of badnavirus particles in the cv. Pilimpikou and identified their equivalent sequences to be of plant genome origin. Probe Gn155Dr, however, hybridised to viral particles and plant genomic DNA in three D. rotundata samples from Guinea. These results represent the first data demonstrating the presence of integrated sequences of badnaviruses in yam. The implications of this for virus-indexing, breeding and multiplication of seed yams are discussed.

A possible scenario for the evolution of Banana streak virus in banana

Outbreaks of Banana streak virus (BSV) have been recorded worldwide where Musa spp. is grown during the last 20 years with no convincing evidence of epidemics. Epidemics were previously reported in Uganda where BSV is currently endemic. BSV is a plant pararetrovirus of the family Caulimoviridae, genus Badnavirus it causes chlorosis leaf streak disease. The information currently available on banana streak disease makes it possible to identify a complex of distinct BSV species each causing the same disease. BSV exists in two states: one as an episomal form, infecting plant cells; the other as viral DNA integrated within the B genome of banana (endogenous BSV-eBSV) forming a viral genome for de novo viral particles. Both forms can be infectious in banana plants. The BSV phylogeny is polyphyletic with BSV distributed in two clades. Clade 1 clusters BSV species that occur worldwide and may have an eBSV counterpart, whereas Clade 3 only comprises BSV species from Uganda. Clearly, two distinct origins explain such BSV diversity. However, the epidemiology/outbreaks of BSV remains unclear and the role of eBSV needs to be clarified. In this review, the biodiversity of BSV is explained and discussed in the light of field and molecular epidemiology data. A scheme is proposed for the co-evolution of BSV and banana based on old or recent infection hypotheses related to African domestication sites and banana dissemination to explain the disease context.