Tracheophyte genomes keep track of the deep evolution of the Caulimoviridae

Characterization of Caulimovirid-like Sequences from Upland Cotton (Gossypium hirsutum L.) Exhibiting Terminal Abortion in Georgia, USA

In this study, we investigated the potential involvement of endogenous viral elements (EVEs) in the development of apical tissue necrosis, resulting in the terminal abortion of upland cotton (Gossypium hirsutum L.) in Georgia. The high-throughput sequence analysis of symptomatic and asymptomatic plant tissue samples revealed near-complete EVE-Georgia (EVE-GA) sequences closely related to caulimoviruses. The analysis of EVE-GA's putative open reading frames (ORFs) compared to cotton virus A and endogenous cotton pararetroviral elements (eCPRVE) revealed their similarity in putative ORFs 1-4. However, in the ORF 5 and ORF 6 encoding putative coat protein and reverse transcriptase, respectively, the sequences from EVE-GA have stop codons similar to eCPRVE sequences from Mississippi. In silico mining of the cotton genome database using EVE-GA as a query uncovered near-complete viral sequence insertions in the genomes of G. hirsutum species (~7 kb) but partial in G. tomentosum (~5.3 kb) and G. mustelinum (~5.1 kb) species. Furthermore, cotton EVEs' episomal forms and messenger RNA (mRNA) transcripts were detected in both symptomatic and asymptomatic plants collected from cotton fields. No significant yield difference was observed between symptomatic and asymptomatic plants of the two varieties evaluated in the experimental plot. Additionally, EVEs were also detected in cotton seeds and seedlings. This study emphasizes the need for future research on EVE sequences, their coding capacity, and any potential role in host immunity or pathogenicity.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record

BACKGROUND

Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes.

RESULTS

Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses.

CONCLUSIONS

We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

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.

Endogenous Caulimovirids: Fossils, Zombies, and Living in Plant Genomes

The Caulimoviridae is a family of double-stranded DNA viruses that infect plants. The genomes of most vascular plants contain endogenous caulimovirids (ECVs), a class of repetitive DNA elements that is abundant in some plant genomes, resulting from the integration of viral DNA in the chromosomes of germline cells during episodes of infection that have sometimes occurred millions of years ago. In this review, we reflect on 25 years of research on ECVs that has shown that members of the Caulimoviridae have occupied an unprecedented range of ecological niches over time and shed light on their diversity and macroevolution. We highlight gaps in knowledge and prospects of future research fueled by increased access to plant genome sequence data and new tools for genome annotation for addressing the extent, impact, and role of ECVs on plant biology and the origin and evolutionary trajectories of the Caulimoviridae.

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.

CAULIFINDER: a pipeline for the automated detection and annotation of caulimovirid endogenous viral elements in plant genomes

Plant, animal and protist genomes often contain endogenous viral elements (EVEs), which correspond to partial and sometimes entire viral genomes that have been captured in the genome of their host organism through a variety of integration mechanisms. While the number of sequenced eukaryotic genomes is rapidly increasing, the annotation and characterization of EVEs remains largely overlooked. EVEs that derive from members of the family Caulimoviridae are widespread across tracheophyte plants, and sometimes they occur in very high copy numbers. However, existing programs for annotating repetitive DNA elements in plant genomes are poor at identifying and then classifying these EVEs. Other than accurately annotating plant genomes, there is intrinsic value in a tool that could identify caulimovirid EVEs as they testify to recent or ancient host-virus interactions and provide valuable insights into virus evolution. In response to this research need, we have developed CAULIFINDER, an automated and sensitive annotation software package. CAULIFINDER consists of two complementary workflows, one to reconstruct, annotate and group caulimovirid EVEs in a given plant genome and the second to classify these genetic elements into officially recognized or tentative genera in the Caulimoviridae. We have benchmarked the CAULIFINDER package using the Vitis vinifera reference genome, which contains a rich assortment of caulimovirid EVEs that have previously been characterized using manual methods. The CAULIFINDER package is distributed in the form of a Docker image.

Rearranged Endogenized Plant Pararetroviruses as Evidence of Heritable RNA-based Immunity

Eukaryotic genomics frequently revealed historical spontaneous endogenization events of external invading nucleic acids, such as viral elements. In plants, an extensive occurrence of endogenous plant pararetroviruses (EPRVs) is usually believed to endow hosts with an additional layer of internal suppressive weaponry. However, an actual demonstration of this activity remains speculative. We analyzed the EPRV component and accompanying silencing effectors of Solanum lycopersicum, documenting that intronic/intergenic pararetroviral integrations bearing inverted-repeats fuel the plant's RNA-based immune system with suitable transcripts capable of evoking a silencing response. A surprisingly small set of rearrangements explained a substantial fraction of pararetroviral-derived endogenous small-interfering (si)RNAs, enriched in 22-nt forms typically associated with anti-viral post-transcriptional gene silencing. We provide preliminary evidence that such genetic and immunological signals may be found in other species outside the genus Solanum. Based on molecular dating, bioinformatics, and empirical explorations, we propose that homology-dependent silencing emerging from particular immuno-competent rearranged chromosomal areas that constitute an adaptive heritable trans-acting record of past infections, with potential impact against the unlocking of plant latent EPRVs and cognate-free pararetroviruses.

A glimpse into the DNA virome of the unique "living fossil" Welwitschia mirabilis

Here, we report the identification and characterization of four novel DNA viruses from Welwitschia mirabilis transcriptomic and genomic datasets. Complete circular virus-like sequences with affinity to members of the Caulimoviridae and Geminiviridae families were detected and characterized from Welwitschia mirabilis genomic data. The two newly members of the Caulimoviridae family have been tentatively named as Welwitschia mirabilis virus 1 and 2 (WMV1-WMV2); whereas the two identified geminiviruses were named as Welwitschia mirabilis associated geminivirus A and B (WMaGVA-WMaGVB). Phylogenetic analysis suggests that WMV1-2 belong to a proposed genus of Caulimoviridae-infecting gymnosperms. WMaGVA-B are phylogenetically related with both mastreviruses and capulaviruses and likely represent a distinct evolutionary lineage within geminiviruses. Additionally, we detected several endogenous virus-like elements (EVE) linked to the discovered viruses in the recently reported W. mirabilis genome, suggesting a shared ancient evolutionary history of these viruses and the Welwithschia.

Expanding the environmental virome: Infection profile in a native rainforest tree species

Agroforestry systems (AFS) for cocoa production combine traditional land-use practices with local biodiversity conservation, resulting in both ecological and agricultural benefits. The cacao-cabruca AFS model is widely implemented in regions of the Brazilian Atlantic Forest. Carpotroche brasiliensis (Raddi) A. Gray (Achariaceae) is a tree found in cabruca landscapes that is often used for reforestation and biotechnological applications. Despite its importance, we still lack information about viruses circulating in C. brasiliensis, particularly considering the possibility of spillover that could affect cocoa production. In our study, we analyzed the Carpotroche brasiliensis virome from Atlantic Forest and cacao-cabruca AFS regions using metatranscriptomics from several vegetative and reproductive organs. Our results revealed a diverse virome detecting near-complete or partial coding sequences of single- and double-stranded DNA and RNA viruses classified into at least six families (Botourmiaviridae, Bromoviridae, Caulimoviridae, Genomoviridae, Mitoviridae, and Rhabdoviridae) plus unclassified elements. We described with high confidence the near-complete and the partial genomes of two tentative novel viruses: Carpotroche-associated ilarvirus and Carpotroche-associated genomovirus, respectively. Interestingly, we also described sequences likely derived from a rhabdovirus, which could represent a novel member of the genus Gammanucleorhabdovirus. We observed higher viral diversity in cacao-cabruca AFS and reproductive organs of C. brasiliensis with preferential tropism to fruits, which could directly affect production. Altogether, our results provide data to better understand the virome in this unexplored agroecological interface, such as cacao-cabruca AFS and forest ecosystem, providing information on the aspects of virus–plant interactions.

Risk Assessment of Infectious Endogenous Banana Streak Viruses in Guadeloupe

Infectious alleles of endogenous banana streak viruses (eBSVs) are present in the genome of all banana interspecific cultivars, including plantains and cooking types. Activation of these infectious eBSV alleles by biotic and abiotic stresses leads to spontaneous infections by cognate viruses and raises concerns about their ability to promote outbreaks of banana streak viruses under field cultivation conditions. We undertook a comprehensive risk assessment study of infectious eBSV alleles of species BSOLV, BSGFV and BSIMV in banana interspecific cultivars in Guadeloupe, a tropical island of the Caribbean where bananas are grown for export and local markets. We carried out a prevalence survey of BSOLV, BSGFV and BSIMV species in a range of cultivars grown in Guadeloupe. Our results suggest that BSOLV and BSGFV infections arise from the activation of infectious eBSVs rather than vector-borne transmission and point to a correlation between altitude and infection rates in interspecific hybrids with AAB genotypes. We studied the dynamics of activation of infectious eBSOLV and eBSGFV alleles by tissue culture and field cultivation in a range of cultivars. We showed that tissue culture and field cultivation trigger distinct activation pathways, resulting in distinct activation patterns. We also showed that activation decreased over time during cell culture and field cultivation and that BSV infections arising from the activation of infectious eBSV alleles cause symptomless infections in the most cultivated plantain in Guadeloupe, French Clair. Overall, our study shows that the risk of BSV outbreaks resulting from the activation of infectious eBSVs in plantain originating from vegetative multiplication is negligible in Guadeloupe.

The diversity of endogenous viral elements in insects

We provide an overview of the currently known diversity of viral sequences integrated into insect genomes. Such endogenous viral elements (EVE) have so far been annotated in at least eight insect orders and can be assigned to at least three families of large double-stranded (ds) DNA viruses, at least 22 families of RNA viruses, and three families of single-stranded DNA viruses. The study of these EVE has already produced important insights into insect-virus interactions, including the discovery of a new form of adaptive antiviral immunity. Insect EVE diversity will continue to increase as new insect genomes and exogenous viruses are sequenced, which will continue to make paleovirology a vibrant research field in this group of animals in the years to come.

Viruses Infecting Trees and Herbs That Produce Edible Fleshy Fruits with a Prominent Value in the Global Market: An Evolutionary Perspective

Trees and herbs that produce fruits represent the most valuable agricultural food commodities in the world. However, the yield of these crops is not fully achieved due to biotic factors such as bacteria, fungi, and viruses. Viruses are capable of causing alterations in plant growth and development, thereby impacting the yield of their hosts significantly. In this work, we first compiled the world's most comprehensive list of known edible fruits that fits our definition. Then, plant viruses infecting those trees and herbs that produce fruits with commercial importance in the global market were identified. The identified plant viruses belong to 30 families, most of them containing single-stranded RNA genomes. Importantly, we show the overall picture of the host range for some virus families following an evolutionary approach. Further, the current knowledge about plant-virus interactions, focusing on the main disorders they cause, as well as yield losses, is summarized. Additionally, since accurate diagnosis methods are of pivotal importance for viral diseases control, the current and emerging technologies for the detection of these plant pathogens are described. Finally, the most promising strategies employed to control viral diseases in the field are presented, focusing on solutions that are long-lasting.

Phylogenetic Distribution of WhiB- and Lsr2-Type Regulators in Actinobacteriophage Genomes

Viruses that infect different actinobacterial host species are known as actinobacteriophages. They are composed of highly divergent and mosaic genomes due to frequent gene exchange between their bacterial hosts and related viral species. This is also reflected by the adaptive incorporation of host transcription factors (TFs) into phage regulatory networks. Previous studies discovered Lsr2-type and WhiB-type regulators encoded by actinobacteriophage genomes. However, limited information is available about their distribution, evolution, and impact on host species. In this study, we computationally screened the distribution of known bacterial and phage TFs inside 2951 complete actinobacteriophage genomes and identified 13 different TF domains. Among those, WhiB, Lsr2, MerR, and Cro/CI-like proteins were widespread and found in more than 10% of the analyzed actinobacteriophage genomes. Neighboring genomic context analysis of the whiB and lsr2 loci showed group-specific conservation of gene synteny and potential involvement of these genes in diverse regulatory functions. Both genes were significantly enriched in temperate phages, and the Lsr2-encoding genomes featured an overall lower GC content. Phylogenetic analysis of WhiB and Lsr2 proteins showed the grouping of phage sequences within bacterial clades, suggesting gene acquisition by phages from their bacterial host species or by multiple, independent acquisition events. Overall, our study reports the global distribution of actinobacteriophage regulatory proteins and sheds light on their origin and evolution.

IMPORTANCE Actinobacteriophages are viruses that infect bacterial species of the diverse phylum of Actinobacteria. Phages engage in a close relationship with their bacterial host. This is also reflected by the adoption of genetic material from their host and its incorporation into phage regulatory circuits. In this study, we systematically searched the genomes of actinobacteriophages for the presence of transcription factor domains. We show that proteins belonging to the regulator families of WhiB and Lsr2 belong to the most abundant regulatory proteins encoded by actinobacteriophages. Further phylogenetic analysis shed light on their origin and evolution. Altogether, this study provides an important basis for further experimental investigation of their role in the coordination of the phage life cycle and their interaction with the host regulatory network in this important bacterial phylum.

Characterisation of novel endogenous geminiviral elements in macadamia

Background

The presence of geminivirus sequences in a preliminary analysis of sRNA sequences from the leaves of macadamia trees with abnormal vertical growth (AVG) syndrome was investigated.

Results

A locus of endogenous geminiviral elements (EGE) in the macadamia genome was analysed, and the sequences revealed a high level of deletions and/or partial integrations, thus rendering the EGE transcriptionally inactive. The replication defective EGE in the macadamia genome indicates its inability to be the source of new viral infections and thus cause AVG or any other disease in macadamia. The EGE sequences were detected in two edible Macadamia species that constitute commercial cultivars and the wild germplasm of edible and inedible species of Macadamia. This strongly suggests that the integration preceded speciation of the genus Macadamia. A draft genome of a locus of EGE in Macadamia was developed. The findings of this study provide evidence to suggest the endogenization of the geminiviral sequences in the macadamia genome and the ancestral relationship of EGE with Macadamia in the Proteaceae family. Random mutations accumulating in the EGE inform that the sequence is evolving.

Conclusions

The EGE in Macadamia is inactive and thus not a direct cause of any diseases or syndromes including AVG in macadamia. The insertion of the EGE in the macadamia genome preceded speciation of the genus Macadamia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08174-0.

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.

Genome-wide identification of endogenous viral sequences in alfalfa (Medicago sativa L.)

Endogenous viral elements (EVEs) have been for the most part described in animals and to a less extent in plants. The endogenization was proposed to contribute toward evolution of living organisms via horizontal gene transfer of novel genetic material and resultant genetic diversity. During the last two decades, several full-length and fragmented EVEs of pararetroviral and non-retroviral nature have been identified in different plant genomes, both monocots and eudicots. Prior to this work, no EVEs have been reported in alfalfa (Medicago sativa L.), the most cultivated forage legume in the world. In this study, taking advantage of the most recent developments in the field of alfalfa research, we have assessed alfalfa genome on the presence of viral-related sequences. Our analysis revealed segmented EVEs resembling two dsDNA reverse-transcribing virus species: Soybean chlorotic mottle virus (family Caulimoviridae, genus Soymovirus) and Figwort mosaic virus (family Caulimoviridae, genus Caulimovirus). The EVEs appear to be stable constituents of the host genome and in that capacity could potentially acquire functional roles in alfalfa’s development and response to environmental stresses.

Paleovirology of the DNA viruses of eukaryotes

Paleovirology is the study of ancient viruses and how they have coevolved with their hosts. An increasingly detailed understanding of the diversity, origins, and evolution of the DNA viruses of eukaryotes has been obtained through the lens of paleovirology in recent years. Members of multiple viral families have been found integrated in the genomes of eukaryotes, providing a rich fossil record to study. These elements have extended our knowledge of exogenous viral diversity, host ranges, and the timing of viral evolution, and are revealing the existence of entire new families of eukaryotic integrating dsDNA viruses and transposons. Future work in paleovirology will continue to provide insights into antiviral immunity, viral diversity, and potential applications, and reveal other secrets of the viral world.

Frontiers in the Standardization of the Plant Platform for High Scale Production of Vaccines

The recent COVID-19 pandemic has highlighted the value of technologies that allow a fast setup and production of biopharmaceuticals in emergency situations. The plant factory system can provide a fast response to epidemics/pandemics. Thanks to their scalability and genome plasticity, plants represent advantageous platforms to produce vaccines. Plant systems imply less complicated production processes and quality controls with respect to mammalian and bacterial cells. The expression of vaccines in plants is based on transient or stable transformation systems and the recent progresses in genome editing techniques, based on the CRISPR/Cas method, allow the manipulation of DNA in an efficient, fast, and easy way by introducing specific modifications in specific sites of a genome. Nonetheless, CRISPR/Cas is far away from being fully exploited for vaccine expression in plants. In this review, an overview of the potential conjugation of the renewed vaccine technologies (i.e., virus-like particles-VLPs, and industrialization of the production process) with genome editing to produce vaccines in plants is reported, illustrating the potential advantages in the standardization of the plant platforms, with the overtaking of constancy of large-scale production challenges, facilitating regulatory requirements and expediting the release and commercialization of the vaccine products of genome edited plants.

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.

A Primer on the Analysis of High-Throughput Sequencing Data for Detection of Plant Viruses

High-throughput sequencing (HTS) technologies have become indispensable tools assisting plant virus diagnostics and research thanks to their ability to detect any plant virus in a sample without prior knowledge. As HTS technologies are heavily relying on bioinformatics analysis of the huge amount of generated sequences, it is of utmost importance that researchers can rely on efficient and reliable bioinformatic tools and can understand the principles, advantages, and disadvantages of the tools used. Here, we present a critical overview of the steps involved in HTS as employed for plant virus detection and virome characterization. We start from sample preparation and nucleic acid extraction as appropriate to the chosen HTS strategy, which is followed by basic data analysis requirements, an extensive overview of the in-depth data processing options, and taxonomic classification of viral sequences detected. By presenting the bioinformatic tools and a detailed overview of the consecutive steps that can be used to implement a well-structured HTS data analysis in an easy and accessible way, this paper is targeted at both beginners and expert scientists engaging in HTS plant virome projects.

Insertion of Badnaviral DNA in the Late Blight Resistance Gene (R1a) of Brinjal Eggplant (Solanum melongena)

Endogenous viral elements (EVEs) are widespread in plant genomes. They result from the random integration of viral sequences into host plant genomes by horizontal DNA transfer and have the potential to alter host gene expression. We performed a large-scale search for co-transcripts including caulimovirid and plant sequences in 1,678 plant and 230 algal species and characterized 50 co-transcripts in 45 distinct plant species belonging to lycophytes, ferns, gymnosperms and angiosperms. We found that insertion of badnavirus EVEs along with Ty-1 copia mobile elements occurred into a late blight resistance gene (R1) of brinjal eggplant (Solanum melongena) and wild relatives in genus Solanum and disrupted R1 orthologs. EVEs of two previously unreported badnaviruses were identified in the genome of S. melongena, whereas EVEs from an additional novel badnavirus were identified in the genome of S. aethiopicum, the cultivated scarlet eggplant. Insertion of these viruses in the ancestral lineages of the direct wild relatives of the eggplant would have occurred during the last 3 Myr, further supporting the distinctiveness of the group of the eggplant within the giant genus Solanum.

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.

ICTV Virus Taxonomy Profile: Caulimoviridae

Caulimoviridae is a family of non-enveloped reverse-transcribing plant viruses with non-covalently closed circular dsDNA genomes of 7.1-9.8 kbp in the order Ortervirales. They infect a wide range of monocots and dicots. Some viruses cause economically important diseases of tropical and subtropical crops. Transmission occurs through insect vectors (aphids, mealybugs, leafhoppers, lace bugs) and grafting. Activation of infectious endogenous viral elements occurs in Musa balbisiana, Petunia hybrida and Nicotiana edwardsonii. However, most endogenous caulimovirids are not infectious. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Caulimoviridae, which is available at ictv.global/report/caulimoviridae.

Large-scale survey reveals pervasiveness and potential function of endogenous geminiviral sequences in plants

The family Geminiviridae contains viruses with single-stranded DNA genomes that have been found infecting a wide variety of angiosperm species. The discovery within the last 25 years of endogenous geminivirus-like (EGV) elements within the nuclear genomes of several angiosperms has raised questions relating to the pervasiveness of EGVs and their impacts on host biology. Only a few EGVs have currently been characterized and it remains unclear whether any of these have influenced, or are currently influencing, the evolutionary fitness of their hosts. We therefore undertook a large-scale search for evidence of EGVs within 134 genome and 797 transcriptome sequences of green plant species. We detected homologues of geminivirus replication-associated protein (Rep) genes in forty-two angiosperm species, including two monocots, thirty-nine dicots, and one ANITA-grade basal angiosperm species (Amborella trichopoda). While EGVs were present in the members of many different plant orders, they were particularly common within the large and diverse order, Ericales, with the highest copy numbers of EGVs being found in two varieties of tea plant (Camellia sinensis). Phylogenetic and clustering analyses revealed multiple highly divergent previously unknown geminivirus Rep lineages, two of which occur in C.sinensis alone. We find that some of the Camellia EGVs are likely transcriptionally active, sometimes co-transcribed with the same host genes across several Camellia species. Overall, our analyses expand the known breadths of both geminivirus diversity and geminivirus host ranges, and strengthens support for the hypothesis that EGVs impact the biology of their hosts.

Genomic characterisation of a newly identified badnavirus infecting ivy (Hedera helix)

High-throughput sequencing (HTS) was used to investigate ringspots on ivy (Hedera helix) leaves. De novo assembly of HTS data generated from a total RNA extract from these leaves yielded a contig with sequence similarity to viruses of the genus Badnavirus, family Caulimoviridae. The complete genome sequence of this virus consists of 8,885 nucleotides and has three open reading frames (ORFs). Genome organisation and phylogenetic analysis identifies this newly identified virus as a new member of the genus Badnavirus for which we propose the name "ivy ringspot-associated virus" (IRSaV).

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.

Genetic differences between Korean and American isolates of Petunia vein clearing virus

Petunia plants are used for urban landscaping in many parts of the world, including South Korea. In this study, we aimed to investigate the occurrence of petunia vein clearing virus (PVCV) infection in petunia plants in Seoul, South Korea. PVCV was detected from 23 of 79 petunia samples collected from Seoul. We obtained the complete genome sequences of the Korean isolates in this study (called PVCV-Kr, Kr2, and Kr3), which were compared with the genome sequence of the USA isolate of the virus (PVCV-USA). The genomic DNA of the three PVCV isolates was found to comprise 7210-7267 nucleotides (nts), which is 4-15 nts longer than the PVCV-USA genome. The genomes of the Kr and Kr2 isolates encode a large polyprotein of 252 kDa (2180 amino acids (aa)). The genome of the Kr3 isolate encodes a large polyprotein of 255 kDa (2203 aa). The polyprotein has six protein domains: a movement protein (MP; 72 aa), a coiled-coil domain (CC; 33 aa), an RNA-binding domain (RB; 18 aa), a protease (PR; 21 aa), a reverse transcriptase (RT; 196 aa), and an RNase H (RH; 121 aa). The large polyprotein and six domains of the three isolates showed 93.9-100.0% sequence homology with those of PVCV-USA. Furthermore, the polymerase polyprotein gene (PR, RT, and RH) of the four PVCV isolates containing the USA isolate grouped with those of Rice tungro bacilliform virus and Soybean chlorotic mottle virus, which belong to the same family (Caulimoviridae). Our findings suggested that the Korean isolates represent a new isolate of PVCV. To our knowledge, this is the first report of PVCV detection in South Korea.

Revisiting the Concept of Host Range of Plant Pathogens

Strategies to manage plant disease-from use of resistant varieties to crop rotation, elimination of reservoirs, landscape planning, surveillance, quarantine, risk modeling, and anticipation of disease emergences-all rely on knowledge of pathogen host range. However, awareness of the multitude of factors that influence the outcome of plant-microorganism interactions, the spatial and temporal dynamics of these factors, and the diversity of any given pathogen makes it increasingly challenging to define simple, all-purpose rules to circumscribe the host range of a pathogen. For bacteria, fungi, oomycetes, and viruses, we illustrate that host range is often an overlapping continuum-more so than the separation of discrete pathotypes-and that host jumps are common. By setting the mechanisms of plant-pathogen interactions into the scales of contemporary land use and Earth history, we propose a framework to assess the frontiers of host range for practical applications and research on pathogen evolution.

A New Putative Caulimoviridae Genus Discovered through Air Metagenomics

Members of the Caulimoviridae family are important plant pathogens. These circular double-stranded DNA viruses may integrate into the host genome, although this integration is not required for the viral replication cycle. Here, we describe three complete genomes belonging to a new putative Caulimoviridae genus discovered through air metagenomics.

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.

Ancient Endogenous Pararetroviruses in Oryza Genomes Provide Insights into the Heterogeneity of Viral Gene Macroevolution

Endogenous viral sequences in eukaryotic genomes, such as those derived from plant pararetroviruses (PRVs), can serve as genomic fossils to study viral macroevolution. Many aspects of viral evolutionary rates are heterogeneous, including substitution rate differences between genes. However, the evolutionary dynamics of this viral gene rate heterogeneity (GRH) have been rarely examined. Characterizing such GRH may help to elucidate viral adaptive evolution. In this study, based on robust phylogenetic analysis, we determined an ancient endogenous PRV group in Oryza genomes in the range of being 2.41-15.00 Myr old. We subsequently used this ancient endogenous PRV group and three younger groups to estimate the GRH of PRVs. Long-term substitution rates for the most conserved gene and a divergent gene were 2.69 × 10-8 to 8.07 × 10-8 and 4.72 × 10-8 to 1.42 × 10-7 substitutions/site/year, respectively. On the basis of a direct comparison, a long-term GRH of 1.83-fold was identified between these two genes, which is unexpectedly low and lower than the short-term GRH (>3.40-fold) of PRVs calculated using published data. The lower long-term GRH of PRVs was due to the slightly faster rate decay of divergent genes than of conserved genes during evolution. To the best of our knowledge, we quantified for the first time the long-term GRH of viral genes using paleovirological analyses, and proposed that the GRH of PRVs might be heterogeneous on time scales (time-dependent GRH). Our findings provide special insights into viral gene macroevolution and should encourage a more detailed examination of the viral GRH.

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.

Pest categorisation of 'Blight and blight-like' diseases of citrus

The EFSA Panel on Plant Health performed a pest categorisation of 'Blight and blight-like' for the EU territory. Blight is a major disease of citrus. Similar 'blight-like' diseases are also known (e.g. declinio, declinamiento) and are addressed simultaneously with Blight in the present categorisation. The causal agent(s) remain(s) unknown and the potential role of a recently identified citrus endogenous pararetrovirus (Citrus Blight-associated pararetrovirus, CBaPRV) remains to be established. Transmissibility and ability to produce consistent (although poorly specific) symptoms have been demonstrated and a combination of indirect approaches is used, with limits, for diagnosis. There are large uncertainties on the biology of the causal agent(s) and on the epidemiology of the disease, including the transmission mechanism(s) responsible for the observed field spread. Blight has been reported from North, Central and South America, Africa and Oceania but is not known to occur in the EU. It is listed in Annex IIA of Directive 2000/29EC. It has the potential to enter, establish and spread in the EU territory. The main entry pathway (citrus plants for planting) is closed by existing legislation and entry is only possible on minor pathways (such as illegal import). Blight is a severe disease and a negative impact is expected should it be introduced in the EU, but the magnitude of this negative impact is very difficult to estimate. 'Blight and blight like' satisfies all criteria evaluated by EFSA to qualify as a Union quarantine pest. It does not meet the criterion of being present in the EU to qualify as a Union regulated non-quarantine pest (RNQP). Since the identity of the causal agent(s) of the Blight and blight-like disease(s) and the existence and efficiency of natural spread mechanism(s) remain unknown, large uncertainties affect all aspects of the present pest categorisation.

Exploring the role of cellular homologous of the 30K-superfamily of plant virus movement proteins

Genes orthologous to the 30K-superfamily of movement proteins (MP) from plant viruses have been recently discovered by bioinformatics analyses as integrated elements in the genome of most vascular plants. However, their functional relevance for plants is still unclear. Here, we undertake some preliminary steps into the functional characterization of one of these putative MP genes found in Arabidopsis thaliana. We found that the AtMP gene is expressed at different stages of the plant development, with accumulation being highest in flowers but lowest in mature siliques. We also found down-regulation of the gene may result in a small delay in plant development and in an exacerbation of the negative effect of salinity in germination efficiency. We have also explored whether changes in expression of the endogenous AtMP have any effect on susceptibility to infection with several viruses, and found that the infectivity of tobacco rattle tobravirus was strongly dependent on the expression of the endogenous AtMP. Finally, we have cloned the endogenous MP from four different plant species into an expression vector that allows for specifically assessing their activity as cell-to-cell movement proteins and have shown that though some may still retain the ancestral activity, they do so in a quite inefficient manner, thus suggesting they have acquired a novel function during adaptation to the host genome.