- Short report
- Open Access
Genome-wide identification of endogenous viral sequences in alfalfa (Medicago sativa L.)
Virology Journal volume 18, Article number: 185 (2021)
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.
Endogenous viral elements (EVEs) are partial or entire viral genes or genomes integrated into host chromosomes and inherited as alleles [1, 2]. EVEs may be eventually removed from the host genomes or reach fixation and retained for millions of years [1, 3, 4]. They are instrumental in a gene flow between viruses and eukaryotes and may have a substantial role in the evolution of their hosts [1, 5, 6]. While the majority of retroviral EVEs originate from a mandatory genome integration stage in the life cycle of reverse-transcribing viruses , insertion mechanisms of plant viruses is debated because integration of a DNA copy into the host genome is not required for their replication .
In plant genomes, two categories of EVEs are recognized: endogenous pararetroviruses (EPRVs), derived from the reverse-transcribing double stranded (ds)DNA viruses of the family Caulimoviridae and endogenous non-retroviral elements (ENREs), originating from dsRNA, ssRNA, and ssDNA viruses [4, 5, 7,8,9]. Several mechanisms for the integration of EPRVs into the plant genome were proposed, among them homologous and illegitimate recombination, integration through reverse-transcribed replication intermediates and via short AT-rich motifs flanking viral segments [8, 10, 11].
Insertion mechanisms of the ENREs remain largely uncharacterized and may involve reverse transcription by reverse transcriptase encoded by retrotransposons, recombination events, and transposon-mediated integration . Localization and distribution of EVEs in plant genomes varies; it was reported that they might be located closer to transposable elements, overlap with genes , reside between the genes or within introns . Numerous EVEs of both classes have been identified in the genomes of different plant species, including a close relative of alfalfa, Medicago truncatula [5, 7, 8, 12, 13].
To our knowledge, prior to this work, none were reported in the genome of alfalfa (M. sativa), widely cultivated perennial legume and important agricultural crop. In this study, taking advantage of the most recent developments in the field of alfalfa genomics, we performed a comprehensive, genome-wide screening aimed to identify virus-related sequences integrated into the alfalfa genome.
Two M. sativa genomes, tetraploid  and diploid (http://www.medicagohapmap.org/downloads/cadl), were used to create a standalone BLAST  database. At the time of this project’s initiation, the tetraploid genome was assembled to the chromosome level, while the diploid genomic sequences were available as scaffolds. To identify possible endogenous viral elements in the M. sativa genome we first used plant viral protein sequences from the Uniprot database (https://www.uniprot.org). These sequences were downloaded and piped through the M. sativa tetraploid and diploid genomes using the tBLASTn program. The standalone BLAST search was also used against reference viral database (https://www.ncbi.nlm.nih.gov/genome/viruses/).
The BLAST analyses revealed sequences resembling two plant virus species of the family Caulimoviridae and integrated into both tetraploid and diploid genomes of M. sativa: Soybean chlorotic mottle virus (SbCMV) of the genus Soymovirus and Figwort mosaic virus (FMV) of the genus Caulimovirus. The resultant parameters of the BLASTn output and the coordinates of the detected homology regions in the M. sativa genomes are shown in Table 1 and Additional file 2: Table S1. The two viral genomes (for Genbank accessions see Table 1) were further used in BLASTn search to confirm the location of the integrated viral sequences. All EVEs found in M. sativa tetraploid and diploid genomes are shown in the Additional file 1: File S1.
The SbCMV-like EPRVs were homologous to the genes encoding reverse transcriptase and movement protein (MP) of SbCMV. Since reverse transcriptase sequences were most probably derived from the multiple integrated transposons and could intervene with the analysis and identification of EPRVs, they were excluded from consideration. The SbCMV-like EPRVs found in both tetraploid and diploid genomes and corresponding to the MP overlapped (percent identity 90%). The top BLAST hit for the longest insert from the diploid genome (563 bp) was SbCMV isolate from Japan (X15828.2; 70%; query cover 100%; E-value = 2e−63). Translation of the SbCMV-like sequence resulted in 163 amino acid-long fragment that had 59% identity with SbCMV movement protein (NP_044299.1; query cover 100%; E-value = 4e−66). Low identity levels (more than 20% of both nucleotide and amino acid differences) indicated that this sequence could belong to a new pararetrovirus, which we provisionally named alfalfa endogenous pararetrovirus (AePV).
To confirm an identity of the discovered AePV EPRVs, we performed a phylogenetic analysis with the representative viruses. Phylogenetic trees were deduced from ClustalW alignments of the AePV sequence and the complete nucleotide sequences of the representative species from the family Caulimoviridae and built using MEGA 7 software  with Maximum Likelihood method based on the Tamura-Nei model and bootstrap analysis of 1000 replicates. Phylogenetic analysis clustered AePV together with SbCMV, thus confirming their relationship (Fig. 1). The group branched to another species in the genus Soymovirus, Blueberry red ringspot virus. All other members of the family grouped correspondingly to their respective genera.
Screening the AePV fragment against the genome of Medicago truncatula, a close relative of alfalfa, using the Plant Genomic Resource Phytozome 12 (https://phytozome.jgi.doe.gov/pz/portal.html), resulted in the identification of orthologous sequences located in all eight chromosomes of the species as well as in several scaffolds (Table 2). Representative orthologous sequence from the chromosome 3 of M. truncatula clustered together with the AePV, supporting their close relationship (Fig. 1).
When the translated portion of the AePV or the orthologous sequence of M. truncatula were searched against the genome of M. truncatula (taxid 3880), a first BLASTp hit corresponded to the “viral movement protein” (AES67244.1) that was 97% identical (100% query cover; E-value = 3e-111) to the AePV. This indicates a presence of the AePV EPRVs in both species.
Short FMV-like EPRVs were found in most chromosomes of the tetraploid genome and scaffold assemblies of the diploid genome (Table 1). The sequences overlapped with each other, forming a 93 nt consensus sequence with a minimum/maximum percent match 85%/100% (Additional file 3: Figure S1). In BLASTn search optimized for high similarity, the consensus sequence had only one hit—with the FMV capsid protein (query cover 97%; E-value = 2e−09; identity 80.43%; X06166.1). When BLASTn was optimized to search for the somewhat similar sequences, the best hits also included portions of the capsid proteins of other species in the family Caulimoviridae.
Phylogenetic trees, deduced from the alignment of the FMV-like consensus sequence and complete genomes of the representative members of the family Caulimoviridae, placed FMV-like consensus in the same cluster with the reference sequence of figwort mosaic virus (NC_003554.1; Fig. 2). Thus, while the FMV-like EPRVs certainly belong to Caulimoviridae and are in all probability resembling FMV, the sequences are too short to make any definite conclusions on their exact identity. The search of the FMV-like consensus against genome of M. truncatula resulted in multiple hits (Table 2), implying that the EPRVs are present in the genome of this species as well.
To ensure that identified EPRVs are not artifacts or contaminants in sequence databases, they were randomly examined by PCR following by Sanger sequencing of the PCR products. DNA was extracted using DNeasy Plant kit (Qiagen) from one-week-old seedlings of alfalfa (M. sativa L) cv. Regency SY and common pea (Pisum sativum) cv. Lincoln, germinated in Petri dishes on moistened sterile filter paper. Prior to germination, alfalfa seeds were scarified with H2SO4, surface-sterilized with 70% ethanol for 1 min and rinsed with distilled water. PCR was performed using AmpliTaq Gold 360 DNA Polymerase under conditions recommended by the manufacturer (Thermo Fisher Scientific, USA), in two–three technical replicates. PCR primers used in the assays were specific to the alfalfa sequences located outside of the selected endogenous elements (Additional file 2: Table S1). DNA from common pea (P. sativum cv. Lincoln) served as a negative control. PCR products were purified with Qiagen PCR Purification kit and sequenced at Macrogen USA facility (Macrogen, MD USA). In all cases, PCR led to the amplification of the products of the expected size (Fig. 3). Sequencing of the purified reactions confirmed that they contained targeted regions incorporating EPRVs. These experiments verified computational findings.
In conclusion, we report the identification of two types of EPRVs in publicly available genomes of tetraploid and diploid alfalfa: SbCMV-like endogenous segments, presumably representing a novel virus in the family Caulimoviridae, tentatively named AePV; and FMV-like elements, which could not be classified assuredly due to the small size, although their homology to FMV and other members of the family Caulimoviridae is evident.
While great many EVEs have been found in the genomes of various plant species, as far as we know, none have been reported in alfalfa prior to this work [5, 7, 8, 12, 13]. This is likely due to the absence of the alfalfa genome, which only recently became available .
Pararetroviruses of the family Caulimoviridae identified in this study represent the most common group of EVEs in plants. They have been proposed to be major components of various plant genomes, possibly contributing to evolution of their hosts as sources of novel genetic material . Despite that two viruses, SbCMV and FMV, which are similar to the EPRVs found in this study, have not currently been reported to infect alfalfa, these infections might occur in nature but so far had gone unnoticed. As a matter of fact, our unpublished results indicate that SbCMV and FMV are present among other viral genomes detected by high-throughput sequencing in field alfalfa samples (Nemchinov et al. unpublished).
Based on the BLASTn results and phylogenetic analysis, it can be assumed that the integration events recorded here are unlikely of recent nature and most probably occurred at some point during the species evolution. The diploid genome used in this work originated from CADL alfalfa (Cultivated Alfalfa at the Diploid Level), developed from cultivated tetraploids using haploidy, with breeding and selection . It is, therefore, plausible that EPRVs identified in the CADL genome could be traced back to the tetraploid germplasm. This implies that they have reached fixation in the alfalfa genome, were inherited as alleles, and may presumably carry a functional load [12, 18].
Of two classes of EPRVs, only FMV-like sequences were found in most chromosomal sets of the tetraploid genome (Table 1). This may suggest, although speculatively due to the small size of the FMV inserts, that the integration events occurred before the whole genome duplication and both diploid progenitors of the tetraploid genome  contained EPRVs.
Hypothetically, uncovering the orthologous viral sequences in the genome of M. truncatula may further postpone the timetable of the integration events to at least ~ 5.3 million years ago, before alfalfa divergence from M. truncatula . Therefore, it appears that EPRVs are stable constituents of the host genome. As a result, they could potentially acquire functional roles in alfalfa’s normal growth, organ development, metabolism, and response to environmental stresses . It cannot be ruled out that EPRVs in alfalfa may also represent a source of infection.
Availability of data and materials
Nucleotide sequences supporting reported results can be found in the Additional file 1: File S1.
Endogenous viral elements
Endogenous non-retroviral elements
Soybean chlorotic mottle virus
Figwort mosaic virus
Alfalfa endogenous pararetrovirus
Cultivated Alfalfa at the Diploid Level
Katzourakis A, Gifford RJ, Endogenous viral elements in animal genomes. PLoS Genet. 2010; 6: e1001191.
Holmes EC. The evolution of endogenous viral elements. Cell Host Microbe. 2011;10:368–77.
Bejarano ER, Khashoggi A, Witty M, Lichtenstein C. Integration of multiple repeats of geminiviral DNA into the nuclear genome of tobacco during evolution. Proc Natl Acad Sci USA. 1996;93:7.
Diop SI, Geering ADW, Alfama-Depauw F, Loaec M, Teycheney PY, Maumus F. Tracheophyte genomes keep track of the deep evolution of the Caulimoviridae. Sci Rep. 2018;8:572.
Chiba, S, Kondo H, Tani A, Saisho D, Sakamoto W. et al. Widespread endogenization of genome sequences of non-retroviral RNA viruses into plant genomes. PLoS Pathog. 2011; 7: e1002146.
Skirmuntt EC, Escalera-Zamudio M, Teeling EC, Smith A, Katzourakis A. The Potential role of endogenous viral elements in the evolution of bats as reservoirs for zoonotic viruses. Annu Rev Virol. 2020;7:103–19.
Takahashi H, Fukuhara T, Kitazawa H, Kormelink R. Virus latency and the impact on plants. Front Microbiol. 2019;10:2764.
Staginnus C, Richert-Pöggeler KR. Endogenous pararetroviruses: two-faced travelers in the plant genome. Trends Plant Sci. 2006;11:485–91.
Chu H, Jo Y, Cho WK. Evolution of endogenous non-retroviral genes integrated into plant genomes. Curr Plant Biol. 2014;1:55–9.
Harper G, Hull R, Lockhart L, Olszewski N. Viral sequences integrated into plant genomes. Ann Rev Phytopathol. 2002;40:119–36.
Kunii M, Kanda M, Nagano H, Uyeda I, Kishima Y, Sano Y. Reconstruction of putative DNA virus from endogenous rice tungro bacilliform virus-like sequences in the rice genome: implications for integration and evolution. BMC Genomics. 2004;5:80.
Geering ADW, Maumus F, Copetti D, Choisne N, Zwickl DJ, et al. Endogenous florendoviruses are major components of plant genomes and hallmarks of virus evolution. Nat Commun. 2014;5:5269.
Liu H, Fu Y, Jiang D, Li G, Xie J, et al. Widespread horizontal gene transfer from double-stranded RNA viruses to eukaryotic nuclear genomes. J Virol. 2010;84:11876–87.
Chen H, Zeng Y, Yang Y, Huang L, Tang B, et al. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat Comm. 2020;11:2494.
Altschul SF, Gish W, Miller W, Myers EW. Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990; 215:403–410.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for Bigger Datasets. Mol Biol Evol. 2016; 33:1870–4.
Bingham ET, McCoy TJ. Cultivated alfalfa at the diploid level: origin, reproductive stability, and yield of seed and forage. Crop Sci. 1979;19:97–100.
Mette MF, Kanno T, Aufsatz W, Jakowitsch J, Winden J, Der V, Matzke MA, et al. Endogenous viral sequences and their potential contribution to heritable virus resistance in plants. EMBO J. 2002;21:461–9.
Quiros CF, Bauchan GR. The genus Medicago and the origin of the Medicago sativa Complex. In: Hanson AA, Barnes DK, Hill RR, editors. Alfalfa and Alfalfa improvement. Madison: American Society of Agronomy, Crop Science Society of America, Soil Science Society of America; 1988; p. 93–124.
This project was supported by USDA-ARS, Beltsville Agricultural Research Center, Beltsville, Maryland, USA, CRIS project # 8042-21000-300-00D and by the Institute of Basic Biological Problems, Russian Academy of Sciences, Pushchino, Moscow Region, Russia.
Ethics approval and consent to participate
Consent for publication
All authors consent to the publication of the manuscript.
The authors declare that they have no competing interests.
This paper is dedicated to the memory of A.M. Boutanaev (1956–2021).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. File S1: Endogenous viral sequences identified in tetraploid and diploid genomes of M. sativa
. Table S1, page 1: Coordinates of the M. sativa genomic regions homologous to the SbCMV and FMV sequences. Table S1, page 2: PCR primers used for confirmation of randomly selected endogenous viral sequences
. Figure S1. A, a consensus sequence of the FMV-like EPRVs assembled using SeqMan tool of the DNASTAR software (DNASTAR, Inc. Madison, Wisconsin USA). B, a fragment of the DNASTAR alignment that generated a consensus sequence.
About this article
Cite this article
Boutanaev, A.M., Nemchinov, L.G. Genome-wide identification of endogenous viral sequences in alfalfa (Medicago sativa L.). Virol J 18, 185 (2021). https://doi.org/10.1186/s12985-021-01650-9
- Medicago sativa L.
- Endogenous viral elements
- Soybean chlorotic mottle virus
- Figwort mosaic virus