Endogenous florendoviruses are major components of plant genomes and hallmarks of virus evolution

Construction of a comprehensive library of repeated sequences for the annotation of Citrus genomes

BACKGROUND

The comprehensive annotation of repeated sequences in genomes is an essential prerequisite for studying the dynamics of these sequences over time and their involvement in gene regulation. Currently, the diversity of repeated sequences in Citrus genomes is only partially characterized because the annotations have been performed using heterogeneous bioinformatics tools, each with its specificity and dedicated only to the annotation of transposable elements.

RESULTS

We combined complementary repeat-finding programs including REPET, CAULIFINDER, and TAREAN, to enable the identification of all types of repetitive sequences found in plant genomes, including transposable elements, endogenous caulimovirids, and satellite DNAs. A fine-grained annotation method was first developed to create a consensus sequence library of repeated sequences identified in the genome assemblies of C. medica, C. micrantha, C. reticulata, and C. maxima, the four ancestral parental species involved in the formation of economically valuable cultivated Citrus varieties. A second, faster annotation method was developed to enrich the dataset by adding new repeated sequences retrieved from genome assemblies of other Citrus species and closely related species belonging to the Aurantioideae subfamily. The final reference library contains 3,091 consensus sequences, of which 94.5% are transposable elements. The diversity of endogenous caulimovirids was characterized for the first time within the genus Citrus, contributing 160 consensus sequences to the final reference library. Finally, 10 satellite DNAs were also identified.

CONCLUSION

Combining multiple repeat detection methods enables the comprehensive annotation of all repeated sequences in Citrus genomes. Using the final reference library reported in this work will improve our understanding of the dynamics of repeated sequences during Citrus speciation, particularly following the genome duplication and hybridization events that led to modern cultivars. The exploration of repeat position insertions along chromosomes using the developed web interface, RepeatLoc Citrus, will also make it possible to further investigate the role of transposable elements and endogenous caulimovirids in genome structure and gene regulation in Citrus species.

A multidisciplinary and integrative review of the structural genome and epigenome of Capsicum L. species

MAIN CONCLUSION

We revised and integrated the genomic and epigenomic data into a comparative Capsicum ideogram, evidencing the advances and future perspectives. Capsicum L. (Solanaceae) genome has been characterized concerning karyotype, nuclear and chromosomal genome size, genome sequencing and physical mapping. In addition, the epigenome has been investigated, showing chromosomal distribution of epimarks in histone amino acids. Genetic and epigenetic discoveries have given light to understanding the structure and organization of the Capsicum "omics". In addition, interspecific and intraspecific similarities and diversities have been identified, characterized and compared in taxonomic and evolutive scenarios. The journey through Capsicum studies allows us to know the 2n = 2x = 24 and 2n = 2x = 26 chromosome numbers, as well as the relatively homomorphic karyotype, and the 1C chromosomal DNA content. In addition, Capsicum "omics" diversity has mainly been evidenced from the nuclear 1C value, as well as from repeatome composition and mapping. Like this, Capsicum provides several opportunities for "omics", ecological, agronomic and conservation approaches, as well as subjects that can be used at different levels of education. In this context, we revisit and integrate Capsicum data about the genome size, karyotype, sequencing and cytogenomics, pointing out the progress and impact of this knowledge in taxonomic, evolutive and agronomic contexts. We also noticed gaps, which can be a focus of further studies. From this multidisciplinary and integrative review, we intend to show the beauty and intrigue of the Capsicum genome and epigenome, as well as the outcomes of these similarities and differences.

Synoptic Variation Drives Genetic Diversity and Transmission Mode of Airborne DNA Viruses in Urban Space

Airborne viruses are ubiquitous and play critical roles in maintaining ecosystem balance, however, they remain unexplored. Here, it is aimed to demonstrate that highly diverse airborne viromes carry out specific metabolic functions and use different transmission modes under different air quality conditions. A total of 263.5-Gb data are collected from 13 air samples for viral metagenomic analysis. After assembly and curation, a total of 12 484 viral contigs (1.5-184.2 kb) are assigned to 221 genus-level clades belonging to 47 families, 19 orders, and 15 classes. The composition of viral communities is influenced by weather conditions, with the main biomarker being Caudoviricetes. The most dominant viruses in these air samples belong to the dsDNA Caudoviricetes (54.0%) and ssDNA Repensiviricetes (31.2%) classes. Twelve novel candidate viruses are identified at the order/family/genus levels by alignment of complete genomes and core genes. Notably, Caudoviricetes are highly prevalent in cloudy and smoggy air, whereas Repensiviricetes are highly dominant in sunny and rainy air. Diverse auxiliary metabolic genes of airborne viruses are mainly involved in deoxynucleotide synthesis, implying their unique roles in atmosphere ecosystem. These findings deepen the understanding of the meteorological impacts on viral composition, transmission mode, and ecological roles in the air that we breathe.

Three endogenous pararetrovirus genomes in Medicago sativa

The genome sequences of three related endogenous pararetroviruses were obtained by high-throughput genomic sequencing of Medicago sativa. The genomes were found to be integrated within plant genes. The phylogeny revealed that Caulimovirus-MSa3 was closely related to caulimoviruses of petunia, whereas Caulimovirus-MSa1 and Caulimovirus-MSa2 were distinct from constructed clades.

Novel genetic variants of banana streak MY virus and banana streak IM virus naturally infecting banana in Northeast India

Divergent banana streak viruses (BSV) were characterized from banana plants exhibiting diverse symptoms in the Northeast region (NER) of India. Using rolling circle amplification (RCA), the complete genome sequences of seven episomal banana streak MY virus (BSMYV) isolates, including two novel variants, and two new banana streak IM virus (BSIMV) isolates were characterized. The novel BSMYV genetic variants were associated with conspicuous necrosis on newly emerged leaves, peduncle distortion, pseudostem internal necrosis, in addition to common streak symptoms. For complete genome nucleotide sequences, BSMYV-IN4 and IN5 shared 77-79% identity with other BSMYVs, while BSMYV-IN7 and IN8 exhibited identities of 77-97%. This study reports for the first time, the complete genomes of two banana streak IM virus (BSIMV-IN1 and -IN2) infecting triploid banana hybrids exhibiting leaf distortion, stunted rosette-like growth, and necrosis, sharing 87% sequence identity with reference BSIMV genome (GenBank accession no. HQ593112). Phylogenetic inference based on complete genomes revealed the distinct and congruent placement of BSMYV-IN4 and IN5 within the BSMYV cluster. Pairwise sequence comparisons of the conserved RT/RNase H nucleotide (nt) sequences revealed that the BSMYV-IN7 and IN4 isolates showed 85% and 97% identity to BSMYV (AY805074), respectively, which shared highest nt identity with BSMYV-IN6, IN9, and IN10, at 100%. The RT/RNase H nt sequences of BSIMV-IN1 and IN2 had 98% identity with the BSIMV (HQ593112), but were characterized as novel variants of BSIMV based on complete genomes. An analysis of relative synonymous codon usage (RSCU) pattern in the ORFIII polyprotein of BSMYV and BSIMV isolates revealed AGA and AGG (arginine) as the most frequently overrepresented codons (>1.5), evolutionary conserved in the genome of both species. A total of 14 recombination events were detected among the 36 BSV genomes, with recombination breakpoints mainly located in the ORFI, III, and IGR genomic regions. A novel phylogenetic cluster, comprised of BSMYV-IN4 and IN5 within the clade I was probably derived from heterologous recombination between parents resembling banana streak VN virus (BSVNV; AY750155) and banana streak GF virus (BSGFV; KJ013507) isolates. The present study conclusively reports the infection of genetically and symptomatically distinct variants of BSMYV and BSIMV infecting banana hybrids in NER India.

High-Throughput Sequencing Identified Multiple Fig Viruses and Viroids Associated with Fig Mosaic Disease in Iraq

Mosaic is the most common viral disease affecting fig plants. Although the Fig mosaic virus is the leading cause of mosaic disease, other viruses are also involved. High-throughput sequencing was used to assess viral infections in fig plants with mosaic. The genomic DNA and total RNAseq of mosaic-symptomatic fig leaves were sequenced using the Illumina platform. The analysis revealed the presence of fig badnavirus 1 (FBV-1), grapevine badnavirus 1 (GBV-1), citrus exocortis viroid (CEVd), and apple dimple fruit viroid (ADFVd). The FBV-1 and GBV-1 sequences were 7,140 bp and 7,239 bp long, respectively. The two genomes encode one open reading frame containing five major protein domains. The viroids, CEVd and ADFVd, were 397 bp and 305 bp long. Phylogenetic analyses revealed a close relationship between FBV-1 and Iranian isolates of the same species, while GBV-1 was closely related to Russian grapevine badnavirus isolates (Tem64, Blu17, KDH48, and Pal9). CEVd was closely related to other Iraqi isolates, while ADFVd was strongly related to a Spanish isolate. A registered endogenous pararetrovirus, caulimovirus-Fca1, with a size of 7,556 bp, was found in the RNA transcripts with a low expression level. This integrant was also detected in the genomes of the two lines 'Horaishi' (a female line) and 'Caprifig 6085' (a male line). Phylogenetic analyses revealed that caulimovirus-Fca1 was distinct from two other clades of different endogenous virus genera.

Study of Endogenous Viruses in the Strawberry Plants

Endogenous viral elements (EVEs) have been reported to exist widely in the genomes of eukaryotic organisms, and they are closely associated with the growth, development, genetics, adaptation, and evolution of their hosts. In this study, two methods-homologous sequence search and genome alignment-were used to explore the endogenous viral sequences in the genomes of Fragaria species. Results revealed abundant endogenous pararetroviruses (EPRVs) in the genomes of Fragaria species, including 786 sequences belonging to five known taxa such as Caulimovirus and other unclassified taxa. Differences were observed in the detected EPRVs between the two methods, with the homologous sequence search having a greater number of EPRVs. On the contrary, genome alignment identified various types and sources of virus-like sequences. Furthermore, through genome alignment, a 267-bp sequence with 95% similarity to the gene encoding the aphid-transmitted protein of Strawberry vein banding virus (Caulimovirus venafragariae) was discovered in the F. chiloensis genome, which was likely a recent insertion. In addition, the statistical analysis of the genome alignment results indicated a remarkably higher abundance of virus-like sequences in the genomes of polyploid strawberries compared with diploid ones. Moreover, the differences in virus-like sequences were observed between the genomes of Fragaria species and those of their close relatives. This study enriched the diversity of viruses that infect strawberries, and laid a theoretical foundation for further research on the origin of endogenous viruses in the strawberry genome, host-virus interactions, adaptation, evolution, and their functions.

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.

Distribution and Location of BEVs in Different Genotypes of Bananas Reveal the Coevolution of BSVs and Bananas

Members of the family Caulimoviridae contain abundant endogenous pararetroviral sequences (EPRVs) integrated into the host genome. Banana streak virus (BSV), a member of the genus Badnavirus in this family, has two distinct badnaviral integrated sequences, endogenous BSV (eBSV) and banana endogenous badnavirus sequences (BEVs). BEVs are distributed widely across the genomes of different genotypes of bananas. To clarify the distribution and location of BEVs in different genotypes of bananas and their coevolutionary relationship with bananas and BSVs, BEVs and BSVs were identified in 102 collected banana samples, and a total of 327 BEVs were obtained and categorized into 26 BEVs species with different detection rates. However, the majority of BEVs were found in Clade II, and a few were clustered in Clade I. Additionally, BEVs and BSVs shared five common conserved motifs. However, BEVs had two unique amino acids, methionine and lysine, which differed from BSVs. BEVs were distributed unequally on most of chromosomes and formed hotspots. Interestingly, a colinear relationship of BEVs was found between AA and BB, as well as AA and SS genotypes of bananas. Notably, the chromosome integration time of different BEVs varied. Based on our findings, we propose that the coevolution of bananas and BSVs is driven by BSV Driving Force (BDF), a complex interaction between BSVs, eBSVs, and BEVs. This study provides the first clarification of the relationship between BEVs and the coevolution of BSVs and bananas in China.

Introduction of Plant Transposon Annotation for Beginners

Transposons are mobile DNA sequences that contribute large fractions of many plant genomes. They provide exclusive resources for tracking gene and genome evolution and for developing molecular tools for basic and applied research. Despite extensive efforts, it is still challenging to accurately annotate transposons, especially for beginners, as transposon prediction requires necessary expertise in both transposon biology and bioinformatics. Moreover, the complexity of plant genomes and the dynamic evolution of transposons also bring difficulties for genome-wide transposon discovery. This review summarizes the three major strategies for transposon detection including repeat-based, structure-based, and homology-based annotation, and introduces the transposon superfamilies identified in plants thus far, and some related bioinformatics resources for detecting plant transposons. Furthermore, it describes transposon classification and explains why the terms 'autonomous' and 'non-autonomous' cannot be used to classify the superfamilies of transposons. Lastly, this review also discusses how to identify misannotated transposons and improve the quality of the transposon database. This review provides helpful information about plant transposons and a beginner's guide on annotating these repetitive sequences.

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.

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.

Systems biology of the genomes' microsatellite signature of Orthopoxvirus including the Monkeypox virus

This study is an attempt to extract and analyse the microsatellites or simple sequence repeats (SSRs) from the genomes of eight species of the genus Orthopoxvirus. The average size of genomes included in the study was 205 kb while the GC% was 33% for all but one. A total of 10,584 SSRs and 854 cSSRs were observed. POX2 with the largest genome of 224.499 kb had maximum of 1493 SSRs and 121 cSSRs (compound SSR) while POX7 with the smallest genome of 185.578 kb had minimum incident SSRs and cSSRs at 1181 and 96, respectively. There was significant correlation between genome size and SSR incidence. Di-nucleotide repeats were the most prevalent (57.47%) followed by mono- at 33% and tri- at 8.6%. Mono-nucleotide SSRs were predominantly T (51%) and A (48.4%). A majority of 80.32% SSRs were in the coding region. The three most similar genomes as per heat map POX1, POX7 and POX5 (93% similarity) are adjacent to one another in the phylogenetic tree. Ankyrin/Ankyrin like protein and Kelch protein which are associated with host determination and divergence have the highest SSR density in almost all studied viruses. Thus, SSRs are involved in genome evolution and host determination of viruses.

Viromic approach reveals differences in the composition, diversity and relative abundance of pumpkin viruses across main growing regions of China

China is the leading country for pumpkin production in the world. As other cucurbits, diseases caused by viruses are among the serious threats to pumpkin production, but our knowledge on the virus species infecting pumpkin plants is fragmentary. To understand the occurrence of viral diseases on pumpkin, we determined the geographical distribution characteristics, relative abundance and evolutionary relationship of pumpkin infected viruses by meta-transcriptome sequencing (RNA-seq) and viromic analysis of 159 samples exhibited typical viral symptoms collected across China in this study. Totally, 11 known and 3 new viruses were identified. Interestingly, 3 new viruses identified in this study should be positive-sense single-stranded RNA virus whose hosts are prokaryotes. The viruses identified in different sampling locations exhibit significant variations in term of virus species and relative abundance. Here, the results provide valuable information for understanding the virus species and their diversity in cultivated pumpkin across major growing regions of China.

Cacao Swollen Shoot Viruses in Ghana

Cacao swollen shoot virus causes cacao swollen shoot disease of Theobroma cacao (cacao) plants. At least six cacao-infecting Badnavirus species-Cacao swollen shoot Togo A virus, Cacao swollen shoot Togo B virus (previously known as Cacao swollen shoot virus), Cacao swollen shoot CE virus, Cacao swollen shoot Ghana M virus, Cacao swollen shoot Ghana N virus, and Cacao swollen shoot Ghana Q virus-are responsible for the swollen shoot disease of cacao in Ghana. Each of these species consists of a multiplicity of strains. The New Juaben strain, the most virulent cacao swollen shoot virus strain in Ghana, belongs to the Cacao swollen shoot Togo B virus species, and is a commonly used strain in laboratory transmission assays. Infection of cacao trees with multiple strains of the virus is common and new evidence suggests that these coinfections may have resulted in the emergence of recombinant strains of the virus. The impact of these emerging recombinant strains on disease severity is uncertain. This review focuses largely on the discovery of cacao swollen shoot virus in Ghana, diversity of the virus strains, molecular characterization, propagation of virus infection in cacao plants, emergence of recombinant virus strains, vector-mediated transmission of the virus, and the management of the cacao swollen shoot disease in Ghana. It also contains sections on the botany and origin of the cacao tree, its introduction to Ghana, the role of cacao swollen shoot disease in facilitating Ghana's independence from Britain, and a brief history of chocolate.

Teratorn and its relatives - a cross-point of distinct mobile elements, transposons and viruses

Mobile genetic elements (e.g., transposable elements and plasmids) and viruses display significant diversity with various life cycles, but how this diversity emerges remains obscure. We previously reported a novel and giant (180 kb long) mobile element, Teratorn, originally identified in the genome of medaka, Oryzias latipes. Teratorn is a composite DNA transposon created by a fusion of a piggyBac-like DNA transposon (piggyBac) and a novel herpesvirus of the Alloherpesviridae family. Genomic survey revealed that Teratorn-like herpesviruses are widely distributed among teleost genomes, the majority of which are also fused with piggyBac, suggesting that fusion with piggyBac is a trigger for the life-cycle shift of authentic herpesviruses to an intragenomic parasite. Thus, Teratorn-like herpesvirus provides a clear example of how novel mobile elements emerge, that is to say, the creation of diversity. In this review, we discuss the unique sequence and life-cycle characteristics of Teratorn, followed by the evolutionary process of piggyBac-herpesvirus fusion based on the distribution of Teratorn-like herpesviruses (relatives) among teleosts. Finally, we provide other examples of evolutionary associations between different classes of elements and propose that recombination could be a driving force generating novel mobile elements.

Teratorn and Its Related Elements – a Novel Group of Herpesviruses Widespread in Teleost Genomes

Herpesviruses are a large family of DNA viruses infecting vertebrates and invertebrates, and are important pathogens in the field of aquaculture. In general, herpesviruses do not have the ability to integrate into the host genomes since they do not have a chromosomal integration step in their life cycles. Recently, we identified a novel group of herpesviruses, "Teratorn" and its related elements, in the genomes of various teleost fish species. At least some of the Teratorn-like herpesviruses are fused with a piggyBac-like DNA transposon, suggesting that they have acquired the transposon-like intragenomic lifestyle by hijacking the transposon system. In this review, we describe the sequence characteristics of Teratorn-like herpesviruses and phylogenetic relationships with other herpesviruses. Then we discuss the process of transposon-herpesvirus fusion, their life cycle, and the generality of transposon-virus fusion. Teratorn-like herpesviruses provide a piece of concrete evidence that even non-retroviral elements can become intragenomic parasites retaining replication capacity, by acquiring transposition machinery from other sources.

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.

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.

An atlas of small RNAs from potato

Small RNAs, including microRNAs (miRNAs), phased secondary small interfering RNAs (phasiRNA), and heterochromatic small interfering RNAs (hc-siRNA) are an essential component of gene regulation. To establish a broad potato small RNA atlas, we constructed an expression atlas of leaves, flowers, roots, and tubers of Desiree and Eva, which are commercially important potato (Solanum tuberosum) cultivars. All small RNAs identified were observed to be conserved between both cultivars, supporting the hypothesis that small RNAs have a low evolutionary rate and are mostly conserved between lineages. However, we also found that a few miRNAs showed differential accumulation between the two potato cultivars, and that hc-siRNAs have a tissue specific expression. We further identified dozens of reproductive and non-reproductive phasiRNAs originating from coding and noncoding regions that appeared to exhibit tissue-specific expression. Together, this study provides an extensive small RNA profiling of different potato tissues that might be used as a resource for future investigations.

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.

Exploring the Emergence and Evolution of Plant Pathogenic Microbes Using Historical and Paleontological Sources

Biotechnological advances now permit broad exploration of past microbial communities preserved in diverse substrates. Despite biomolecular degradation, high-throughput sequencing of preserved materials can yield invaluable genomic and metagenomic data from the past. This line of research has expanded from its initial human- and animal-centric foci to include plant-associated microbes (viruses, archaea, bacteria, fungi, and oomycetes), for which historical, archaeological, and paleontological data illuminate past epidemics and evolutionary history. Genetic mechanisms underlying the acquisition of microbial pathogenicity, including hybridization, polyploidization, and horizontal gene transfer, can now be reconstructed, as can gene-for-gene coevolution with plant hosts. Epidemiological parameters, such as geographic origin and range expansion, can also be assessed. Building on published case studies with individual phytomicrobial taxa, the stage is now set for broader, community-wide studies of preserved plant microbiomes to strengthen mechanistic understanding of microbial interactions and plant disease emergence.

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.

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.

Diversity of the virome associated with alfalfa (Medicago sativa L.) in the U.S. Pacific Northwest

Alfalfa (Medicago sativa L.) is one of the most extensively cultivated forage legumes in the world. It is currently the third most valuable field crop in the United States with an estimated value of over $9.3 billion. Alfalfa productivity is limited by various infectious diseases that can reduce forage yield and quality and shorten stand life. The crop can frequently be infected with a diverse array of pathogens and other organisms that have distinct life cycles, biology, and mode of action. Among them are many coinfecting viruses, that greatly contribute to the heterogeneity of within-host pathogenic communities, representing a ubiquitous and abundant background for all other host–pathogen interactions. Regrettably, the impact of viral diseases, their role in alfalfa health and involvement in the severity of multi-pathogen infections are often underestimated and not well understood. As high-throughput sequencing approaches have been developed, opportunities to delve into these complex interactions can be realized. In this work, we have characterized a diversity of viral populations in several commercial alfalfa production fields located in the U.S. Pacific Northwest. At least 45 distinct viruses have been identified in all alfalfa samples. Among them some were known to infect the crop prior to this study, and others were designated as emerging, novel and viruses integrated into the alfalfa genome. Known viruses included alfalfa mosaic virus, pea streak virus and bean leafroll virus, while among emerging and novel agents were alfalfa virus S, cherry virus Trakiya, several rhabdoviruses and others. Additional biological and impact studies will be needed to determine if newly identified viruses, especially those that have not been reported from alfalfa before, should be considered pathogens of this crop.

Homing in on Endogenous Badnaviral Elements: Development of Multiplex PCR-DGGE for Detection and Rapid Identification of Badnavirus Sequences in Yam Germplasm

Viruses of the genus Badnavirus (family Caulimoviridae) are double-stranded DNA-reverse transcribing (dsDNA-RT) plant viruses and have emerged as serious pathogens of tropical and temperate crops globally. Endogenous badnaviral sequences are found integrated in the genomes of several economically important plant species. Infection due to activation of replication-competent integrated copies of the genera Badnavirus, Petuvirus and Cavemovirus has been described. Such endogenous badnaviral elements pose challenges to the development of nucleic acid-based diagnostic methods for episomal virus infections and decisions on health certification for international movement of germplasm and seed. One major food security crop affected is yam (Dioscorea spp.). A diverse range of Dioscorea bacilliform viruses (DBVs), and endogenous DBV (eDBV) sequences have been found to be widespread in yams cultivated in West Africa and other parts of the world. This study outlines the development of multiplex PCR-dependent denaturing gradient gel electrophoresis (PCR-DGGE) to assist in the detection and analysis of eDBVs, through the example of analysing yam germplasm from Nigeria and Ghana. Primers targeting the three most prevalent DBV monophyletic species groups in West Africa were designed to improve DGGE resolution of complex eDBV sequence fingerprints. Multiplex PCR-DGGE with the addition of a tailor-made DGGE sequence marker enables rapid comparison of endogenous badnaviral sequence diversity across germplasm, as illustrated in this study for eDBV diversity in yam.

A haplotype resolved chromosomal level avocado genome allows analysis of novel avocado genes

Avocado (Persea americana) is a member of the magnoliids, an early branching lineage of angiosperms that has high value globally with the fruit being highly nutritious. Here, we report a chromosome-level genome assembly for the commercial avocado cultivar Hass, which represents 80% of the world's avocado consumption. The DNA contigs produced from Pacific Biosciences HiFi reads were further assembled using a previously published version of the genome supported by a genetic map. The total assembly was 913 Mb with a contig N50 of 84 Mb. Contigs assigned to the 12 chromosomes represented 874 Mb and covered 98.8% of benchmarked single-copy genes from embryophytes. Annotation of protein coding sequences identified 48 915 avocado genes of which 39 207 could be ascribed functions. The genome contained 62.6% repeat elements. Specific biosynthetic pathways of interest in the genome were investigated. The analysis suggested that the predominant pathway of heptose biosynthesis in avocado may be through sedoheptulose 1,7 bisphosphate rather than via alternative routes. Endoglucanase genes were high in number, consistent with avocado using cellulase for fruit ripening. The avocado genome appeared to have a limited number of translocations between homeologous chromosomes, despite having undergone multiple genome duplication events. Proteome clustering with related species permitted identification of genes unique to avocado and other members of the Lauraceae family, as well as genes unique to species diverged near or prior to the divergence of monocots and eudicots. This genome provides a tool to support future advances in the development of elite avocado varieties with higher yields and fruit quality.

Piper DNA virus 1 and 2 are endogenous pararetroviruses integrated into chromosomes of black pepper (Piper nigrum L)

A previous study named 7178 and 892 bp contigs obtained through high-throughput sequencing (HTS) of black pepper as Piper DNA virus 1 (PDV-1) and PDV-2 respectively. In the present study, HTS results were confirmed through polymerase chain reaction and Sanger sequencing. The sequenced region of both PDV-1 and PDV-2 contained partial genomes with motifs characteristic of pararetroviruses. BLAST analysis of PDV-1 and PDV-2 against the whole genome sequence of the black pepper showed integration of the PDV-1 at 22 loci in chromosome number 14, and PDV-2 at two loci in chromosome number 12 of black pepper. The integration was confirmed through amplification and sequencing of the junction regions. The present study suggests that both PDV-1 and PDV-2 occur as endogenous viruses in black pepper. Further studies are needed to determine whether these endogenous viruses occur in episomal forms, their complete genome sequence and whether they are activable under abiotic stress conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-021-00752-w.

A vertically transmitted amalgavirus is present in certain accessions of the bryophyte Physcomitrium patens

In the last few years, next-generation sequencing techniques have started to be used to identify new viruses infecting plants. This has allowed to rapidly increase our knowledge on viruses other than those causing symptoms in economically important crops. Here we used this approach to identify a virus infecting Physcomitrium patens that has the typical structure of the double-stranded RNA endogenous viruses of the Amalgaviridae family, which we named Physcomitrium patens amalgavirus 1, or PHPAV1. PHPAV1 is present only in certain accessions of P. patens, where its RNA can be detected throughout the cell cycle of the plant. Our analysis demonstrates that PHPAV1 can be vertically transmitted through both paternal and maternal germlines, in crosses between accessions that contain the virus with accessions that do not contain it. This work suggests that PHPAV1 can replicate in genomic backgrounds different from those that actually contain the virus and opens the door for future studies on virus-host coevolution.

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.

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.

Horizontal Gene Transfers in Plants

In plants, as in all eukaryotes, the vertical transmission of genetic information through reproduction ensures the maintenance of the integrity of species. However, many reports over the past few years have clearly shown that horizontal gene transfers, referred to as HGTs (the interspecific transmission of genetic information across reproductive barriers) are very common in nature and concern all living organisms including plants. The advent of next-generation sequencing technologies (NGS) has opened new perspectives for the study of HGTs through comparative genomic approaches. In this review, we provide an up-to-date view of our current knowledge of HGTs in plants.

Diverse Trajectories Drive the Expression of a Giant Virus in the Oomycete Plant Pathogen Phytophthora parasitica

Giant viruses of amoebas, recently classified in the class Megaviricetes, are a group of viruses that can infect major eukaryotic lineages. We previously identified a set of giant virus sequences in the genome of Phytophthora parasitica, an oomycete and a devastating major plant pathogen. How viral insertions shape the structure and evolution of the invaded genomes is unclear, but it is known that the unprecedented functional potential of giant viruses is the result of an intense genetic interplay with their hosts. We previously identified a set of giant virus sequences in the genome of P. parasitica, an oomycete and a devastating major plant pathogen. Here, we show that viral pieces are found in a 550-kb locus and are organized in three main clusters. Viral sequences, namely RNA polymerases I and II and a major capsid protein, were identified, along with orphan sequences, as a hallmark of giant viruses insertions. Mining of public databases and phylogenetic reconstructions suggest an ancient association of oomycetes and giant viruses of amoeba, including faustoviruses, African swine fever virus (ASFV) and pandoraviruses, and that a single viral insertion occurred early in the evolutionary history of oomycetes prior to the Phytophthora–Pythium radiation, estimated at ∼80 million years ago. Functional annotation reveals that the viral insertions are located in a gene sparse region of the Phytophthora genome, characterized by a plethora of transposable elements (TEs), effectors and other genes potentially involved in virulence. Transcription of viral genes was investigated through analysis of RNA-Seq data and qPCR experiments. We show that most viral genes are not expressed, and that a variety of mechanisms, including deletions, TEs insertions and RNA interference may contribute to transcriptional repression. However, a gene coding a truncated copy of RNA polymerase II along a set of neighboring sequences have been shown to be expressed in a wide range of physiological conditions, including responses to stress. These results, which describe for the first time the endogenization of a giant virus in an oomycete, contribute to challenge our view of Phytophthora evolution.

Sixty Years from the First Disease Description, a Novel Badnavirus Associated with Chestnut Mosaic Disease

Although chestnut mosaic disease (ChMD) was described several decades ago, its etiology is still not clear. Using classical approaches and high-throughput sequencing (HTS) techniques, we identified a novel Badnavirus that is a strong etiological candidate for ChMD. Two disease sources from Italy and France were submitted to HTS-based viral indexing. Total RNAs were extracted, ribodepleted, and sequenced on an Illumina NextSeq500 (2 × 150 nt or 2 × 75 nt). In each source, we identified a single contig of ≈7.2 kb that corresponds to a complete circular viral genome and shares homologies with various badnaviruses. The genomes of the two isolates have an average nucleotide identity of 90.5%, with a typical badnaviral genome organization comprising three open reading frames. Phylogenetic analyses and sequence comparisons showed that this virus is a novel species; we propose the name Chestnut mosaic virus (ChMV). Using a newly developed molecular detection test, we systematically detected the virus in symptomatic graft-inoculated indicator plants (chestnut and American oak) as well in chestnut trees presenting typical ChMD symptoms in the field (100 and 87% in France and Italy surveys, respectively). Datamining of publicly available chestnut sequence read archive transcriptomic data allowed the reconstruction of two additional complete ChMV genomes from two Castanea mollissima sources from the United States as well as ChMV detection in C. dentata from the United States. Preliminary epidemiological studies performed in France and central eastern Italy showed that ChMV has a high incidence in some commercial orchards and low within-orchard genetic diversity.

Identification and distribution of novel badnaviral sequences integrated in the genome of cacao (Theobroma cacao)

Theobroma cacao is one of the most economically important tropical trees, being the source of chocolate. As part of an ongoing study to understand the diversity of the badnavirus complex, responsible for the cacao swollen shoot virus disease in West Africa, evidence was found recently of virus-like sequences in asymptomatic cacao plants. The present study exploited the wealth of genomic resources in this crop, and combined bioinformatic, molecular, and genetic approaches to report for the first time the presence of integrated badnaviral sequences in most of the cacao genetic groups. These sequences, which we propose to name eTcBV for endogenous T. cacao bacilliform virus, varied in type with each predominating in a specific genetic group. A diagnostic multiplex PCR method was developed to identify the homozygous or hemizygous condition of one specific insert, which was inherited as a single Mendelian trait. These data suggest that these integration events occurred before or during the species diversification in Central and South America, and prior to its cultivation in other regions. Such evidence of integrated sequences is relevant to the management of cacao quarantine facilities and may also aid novel methods to reduce the impact of such viruses in this crop.

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.

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.