On the estimation of the insertion time of LTR retrotransposable elements

A telomere-to-telomere gap-free assembly integrating multi-omics uncovers the genetic mechanism of fruit quality and important agronomic trait associations in pomegranate

Pomegranate is an important perennial fruit tree distributed worldwide. Reference genomes with gaps and limit gene identification controlling important agronomic traits hinder its functional genomics and genetic improvements. Here, we reported a telomere-to-telomere (T2T) gap-free genome assembly of the distinctive cultivar 'Moshiliu'. The Moshiliu reference genome was assembled into eight chromosomes without gaps, totalling ~366.71 Mb, with 32 158 predicted protein-coding genes. All 16 telomeres and eight centromeres were characterized; combined with FISH analysis, we revealed the atypical telomere units in pomegranate as TTTTAGGG. Furthermore, a total of 16 loci associated with 15 important agronomic traits were identified based on GWAS of 146 accessions. Gene editing and biochemical experiments demonstrated that a 37.2-Kb unique chromosome translocation disrupting the coding domain sequence of PgANS was responsible for anthocyanin-less, knockout of PgANS in pomegranate exhibited a defect in anthocyanin production; a unique repeat expansion in the promoter of PgANR may affected its expression, resulting in black peel; notably, the G → A transversion located at the 166-bp coding domain of PgNST3, which caused a E56K mutation in the PgNST3 protein, closely linked with soft-seed trait. Overexpression of PgNST3A in tomato presented smaller and softer seed coats. The E56K mutation in PgNST3 protein, eliminated the binding ability of PgNST3 to the PgMYB46 promoter, which subsequently affected the thickness of the inner seed coat of soft-seeded pomegranates. Collectively, the validated gap-free genome, the identified genes controlling important traits and the CRISPR-Cas9-mediated gene knockout system all provided invaluable resources for pomegranate precise breeding.

High-quality genome of allotetraploid Avena barbata provides insights into the origin and evolution of B subgenome in Avena

Avena barbata, a wild oat species within the genus Avena, is a widely used model for studying plant ecological adaptation due to its strong environmental adaptability and disease resistance, serving as a valuable genetic resource for oat improvement. Here, we phased the high-quality chromosome-level genome assembly of A. barbata (6.88 Gb, contig N50 = 53.74 Mb) into A (3.57 Gb with 47,687 genes) and B (3.31 Gb with 46,029 genes) subgenomes. Comparative genomics and phylogenomic analyses clarified the evolutionary relationships and trajectories of A, B, C and D subgenomes in Avena. We inferred that the A subgenome donor of A. barbata was Avena hirtula, while the B subgenome donor was probably an extinct diploid species closely related to Avena wiestii. Genome evolution analysis revealed the dynamic transposable element (TE) content and subgenome divergence, as well as extensive structure variations across A, B, C, and D subgenomes in Avena. Population genetic analysis of 211 A. barbata accessions from distinct ecotypes identified several candidate genes related to environmental adaptability and drought resistance. Our study provides a comprehensive genetic resource for exploring the genetic basis underlying the strong environmental adaptability of A. barbata and the molecular identification of important agronomic traits for oat breeding.

Crossing Wallace's line: an evolutionarily young gibbon ape leukemia virus like endogenous retrovirus identified from the Philippine flying lemur (Cynocephalus volans)

Wallace's line is a biogeographical barrier to faunal movements between Southeast Asia and the Australo-Papuan region. There are exceptions among rodents and bats, few of which have crossed Wallace's line. The gibbon ape leukemia viruses (GALV) and koala retrovirus (KoRV) have only been identified in wildlife on the Australo-Papuan side of Wallaces's Line with the potential exception of partial sequences identified in two microbat fecal samples from China and a recently described GALV relative in a rodent from Africa. Here we describe a group of GALV-like endogenous retroviral sequences from the Southeast Asian flying lemur (Cynocephalus volans) representing the first known description of a primate relative which has been infected, and the germline colonized, by GALVs on the Southeast Asian side of Wallace's Line.

Mutation Rate Variation and Other Challenges in 2-LTR Dating of Primate Endogenous Retrovirus Integrations

The time of integration of germline-targeting Long Terminal Repeat (LTR) retroposons, such as endogenous retroviruses (ERVs), can be estimated by assessing the nucleotide divergence between the LTR sequences flanking the viral genes. Due to the viral replication mechanism, both LTRs are identical at the moment of integration, when the provirus becomes part of the host genome. After that time, proviral sequences evolve within the host DNA. When the mutation rate is known, nucleotide divergence between the LTRs would then be a measure of time elapsed since integration. Though frequently used, the approach has been complicated by the choice of host mutation rate and, to a lesser extent, by the method selected to estimate nucleotide divergence. As a result, outcomes can be incompatible with, for instance, speciation events identified from the fossil record. The review will give an overview of research reporting LTR-retroposon dating, and a summary of important factors to consider, including the quality, assembly, and alignment of sequences, the mutation rate of foreign DNA in host genomes, and the choice of a distance estimation method. Primates will here be the focus of the analysis because their genomes, ERVs, and fossil record have been extensively studied. However, most of the factors discussed have a wide applicability in the vertebrate field.

Evolutionary Nonindependence Between Human piRNAs and Their Potential Target Sites in Protein-Coding Genes

PIWI-interacting RNAs (piRNAs) are the most diverse small RNAs in animals. These small RNAs have been known to play an important role in the suppression of transposable elements (TEs). Protein-coding genes (PCGs) are the most well-recognized functional genes in genomes. In the present study, we designed and performed a set of statistics-based evolutionary analyses to reveal nonrandom phenomena in the evolution of human piRNA-PCG targeting relationships. Through analyzing the occurrence of single nucleotide variants (SNVs) in potential piRNA target sites in human PCGs, we provide evidence that there exists a mutational force biased to strengthen piRNA-PCG targeting relationships. Through analyzing the allele frequencies of SNVs in potential piRNA target sites in human PCGs, we provide evidence that there exists a piRNA-dependent selective force acting on potential piRNA target sites in human PCGs. Because of these nonrandom evolutionary forces, human piRNAs and their potential target sites in PCGs are not independent in evolution. Additionally, we found evidence that potential piRNA target sites in human PCGs are particularly likely to be present in regions derived from Alu elements. This finding suggests that the aforementioned evolutionary forces acting on piRNA-PCG targeting relationships could be particularly prone to affect Alu-derived regions in human PCGs. Collectively, our findings provide new insights into the evolutionary interplay between piRNAs, PCGs, and Alu elements in the evolution of the human genome.

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

Recombinant origin and interspecies transmission of a HERV-K(HML-2)-related primate retrovirus with a novel RNA transport element

HERV-K(HML-2), the youngest clade of human endogenous retroviruses (HERVs), includes many intact or nearly intact proviruses, but no replication competent HML-2 proviruses have been identified in humans. HML-2-related proviruses are present in other primates, including rhesus macaques, but the extent and timing of HML-2 activity in macaques remains unclear. We have identified 145 HML-2-like proviruses in rhesus macaques, including a clade of young, rhesus-specific insertions. Age estimates, intact open reading frames, and insertional polymorphism of these insertions are consistent with recent or ongoing infectious activity in macaques. 106 of the proviruses form a clade characterized by an ~750 bp sequence between env and the 3' long terminal repeat (LTR), derived from an ancient recombination with a HERV-K(HML-8)-related virus. This clade is found in Old World monkeys (OWM), but not great apes, suggesting it originated after the ape/OWM split. We identified similar proviruses in white-cheeked gibbons; the gibbon insertions cluster within the OWM recombinant clade, suggesting interspecies transmission from OWM to gibbons. The LTRs of the youngest proviruses have deletions in U3, which disrupt the Rec Response Element (RcRE), required for nuclear export of unspliced viral RNA. We show that the HML-8-derived region functions as a Rec-independent constitutive transport element (CTE), indicating the ancestral Rec-RcRE export system was replaced by a CTE mechanism.

Non-B-form DNA is associated with centromere stability in newly-formed polyploid wheat

Non-B-form DNA differs from the classic B-DNA double helix structure and plays a crucial regulatory role in replication and transcription. However, the role of non-B-form DNA in centromeres, especially in polyploid wheat, remains elusive. Here, we systematically analyzed seven non-B-form DNA motif profiles (A-phased DNA repeat, direct repeat, G-quadruplex, inverted repeat, mirror repeat, short tandem repeat, and Z-DNA) in hexaploid wheat. We found that three of these non-B-form DNA motifs were enriched at centromeric regions, especially at the CENH3-binding sites, suggesting that non-B-form DNA may create a favorable loading environment for the CENH3 nucleosome. To investigate the dynamics of centromeric non-B form DNA during the alloploidization process, we analyzed DNA secondary structure using CENH3 ChIP-seq data from newly formed allotetraploid wheat and its two diploid ancestors. We found that newly formed allotetraploid wheat formed more non-B-form DNA in centromeric regions compared with their parents, suggesting that non-B-form DNA is related to the localization of the centromeric regions in newly formed wheat. Furthermore, non-B-form DNA enriched in the centromeric regions was found to preferentially form on young LTR retrotransposons, explaining CENH3's tendency to bind to younger LTR. Collectively, our study describes the landscape of non-B-form DNA in the wheat genome, and sheds light on its potential role in the evolution of polyploid centromeres.

The genome of the glasshouse plant noble rhubarb (Rheum nobile) provides a window into alpine adaptation

Glasshouse plants are species that trap warmth via specialized morphology and physiology, mimicking a human glasshouse. In the Himalayan alpine region, the highly specialized glasshouse morphology has independently evolved in distinct lineages to adapt to intensive UV radiation and low temperature. Here we demonstrate that the glasshouse structure - specialized cauline leaves - is highly effective in absorbing UV light but transmitting visible and infrared light, creating an optimal microclimate for the development of reproductive organs. We reveal that this glasshouse syndrome has evolved at least three times independently in the rhubarb genus Rheum. We report the genome sequence of the flagship glasshouse plant Rheum nobile and identify key genetic network modules in association with the morphological transition to specialized glasshouse leaves, including active secondary cell wall biogenesis, upregulated cuticular cutin biosynthesis, and suppression of photosynthesis and terpenoid biosynthesis. The distinct cell wall organization and cuticle development might be important for the specialized optical property of glasshouse leaves. We also find that the expansion of LTRs has likely played an important role in noble rhubarb adaptation to high elevation environments. Our study will enable additional comparative analyses to identify the genetic basis underlying the convergent occurrence of glasshouse syndrome.

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.

Evolutionarily Young African Rhinoceros Gammaretroviruses

High-throughput sequences were generated from DNA and cDNA from four Southern white rhinoceros (Ceratotherium simum simum) located in the Taronga Western Plain Zoo in Australia. Virome analysis identified reads that were similar to Mus caroli endogenous gammaretrovirus (McERV). Previous analysis of perissodactyl genomes did not recover gammaretroviruses. Our analysis, including the screening of the updated white rhinoceros (Ceratotherium simum) and black rhinoceros (Diceros bicornis) draft genomes identified high-copy orthologous gammaretroviral ERVs. Screening of Asian rhinoceros, extinct rhinoceros, domestic horse, and tapir genomes did not identify related gammaretroviral sequences in these species. The newly identified proviral sequences were designated SimumERV and DicerosERV for the white and black rhinoceros retroviruses, respectively. Two long terminal repeat (LTR) variants (LTR-A and LTR-B) were identified in the black rhinoceros, with different copy numbers associated with each (n = 101 and 373, respectively). Only the LTR-A lineage (n = 467) was found in the white rhinoceros. The African and Asian rhinoceros lineages diverged approximately 16 million years ago. Divergence age estimation of the identified proviruses suggests that the exogenous retroviral ancestor of the African rhinoceros ERVs colonized their genomes within the last 8 million years, a result consistent with the absence of these gammaretroviruses from Asian rhinoceros and other perissodactyls. The black rhinoceros germ line was colonized by two lineages of closely related retroviruses and white rhinoceros by one. Phylogenetic analysis indicates a close evolutionary relationship with ERVs of rodents including sympatric African rats, suggesting a possible African origin of the identified rhinoceros gammaretroviruses. IMPORTANCE Rhinoceros genomes were thought to be devoid of gammaretroviruses, as has been determined for other perissodactyls (horses, tapirs, and rhinoceros). While this may be true of most rhinoceros, the African white and black rhinoceros genomes have been colonized by evolutionarily young gammaretroviruses (SimumERV and DicerosERV for the white and black rhinoceros, respectively). These high-copy endogenous retroviruses (ERVs) may have expanded in multiple waves. The closest relative of SimumERV and DicerosERV is found in rodents, including African endemic species. Restriction of the ERVs to African rhinoceros suggests an African origin for the rhinoceros gammaretroviruses.

PlantLTRdb: An interactive database for 195 plant species LTR-retrotransposons

LTR-retrotransposons (LTR-RTs) are a large group of transposable elements that replicate through an RNA intermediate and alter genome structure. The activities of LTR-RTs in plant genomes provide helpful information about genome evolution and gene function. LTR-RTs near or within genes can directly alter gene function. This work introduces PlantLTRdb, an intact LTR-RT database for 195 plant species. Using homology- and de novo structure-based methods, a total of 150.18 Gbp representing 3,079,469 pseudomolecules/scaffolds were analyzed to identify, characterize, annotate LTR-RTs, estimate insertion ages, detect LTR-RT-gene chimeras, and determine nearby genes. Accordingly, 520,194 intact LTR-RTs were discovered, including 29,462 autonomous and 490,732 nonautonomous LTR-RTs. The autonomous LTR-RTs included 10,286 Gypsy and 19,176 Copia, while the nonautonomous were divided into 224,906 Gypsy, 218,414 Copia, 1,768 BARE-2, 3,147 TR-GAG and 4,2497 unknown. Analysis of the identified LTR-RTs located within genes showed that a total of 36,236 LTR-RTs were LTR-RT-gene chimeras and 11,619 LTR-RTs were within pseudo-genes. In addition, 50,026 genes are within 1 kbp of LTR-RTs, and 250,587 had a distance of 1 to 10 kbp from LTR-RTs. PlantLTRdb allows researchers to search, visualize, BLAST and analyze plant LTR-RTs. PlantLTRdb can contribute to the understanding of structural variations, genome organization, functional genomics, and the development of LTR-RT target markers for molecular plant breeding. PlantLTRdb is available at https://bioinformatics.um6p.ma/PlantLTRdb.

The discovery of multiple active mys-related LTR-retroelements within the Neotominae subfamily of cricetid rodents

Retrotransposon families in the rodent family Cricetidae have been understudied in contrast to Muridae, both taxa classified within the superfamily Muroidea. Therefore, we carried out a study to advance our knowledge of the unique mys LTR-retroelement identified in Peromyscus leucopus, by incorporating intra-ORF PCR, quantitative dot blots, DNA and protein library screens, the generation of molecular phylogenies, and analyses of orthologous LTR-retroelement loci. These analyses led to the discovery of three additional related families of LTR-retroelements, which include a 2900 bp full-length element of mys-related sequences (mysRS), an 8000 bp element containing the mys ORF1 sequence (mORF1) with ERV-related sequences downstream in the reverse orientation, as well as an 1800 bp element primarily consisting of mys ORF2 (mORF2) related sequences flanked by LTRs. Our data revealed only a few full-length mys elements among genera of the Neotominae subfamily of cricetid rodents, most existing as partial copies. The mysRS and mORF1 elements are also limited to the genomes of the Neotominae subfamily, whereas mORF2 appears to be restricted to the Peromyscus genus. Molecular phylogenies demonstrating concerted evolution along with an assessment of orthologous loci in Peromyscus for the presence or absence of elements are consistent with activity of these novel LTR-retroelement families within this genus. Together with known activity of various families of non-LTR retroelements in Peromyscus species, we propose that retrotransposons have been continually contributing to the dynamics of Peromyscus genomes promoting genomic diversity and may be correlated with the evolution of more than 50 identified Peromyscus species.

Identification of Cartilaginous Fish Endogenous Foamy Virus Rooting to Vertebrate Counterparts

Foamy viruses (FVs) are ideal models for studying the long-term evolutionary history between viruses and their hosts. Currently, FVs have been documented in nearly all major taxa of vertebrates, but evidence is lacking for true FV infiltration in cartilaginous fish, the most basal living vertebrates with jaws. Here, we screened 11 available genomes and 10 transcriptome sequence assemblies of cartilaginous fish and revealed a novel endogenous foamy virus, termed cartilaginous fish endogenous foamy virus (CFEFV), in the genomes of sharks and rays. Genomic analysis of CFEFVs revealed feature motifs that were retained among canonical FVs. Phylogenetic analysis using polymerase sequences revealed the rooting nature of CFEFVs to vertebrate FVs, indicating their deep origin. Interestingly, three viral lineages were found in a shark (Scyliorhinus torazame), one of which was clustered with ray-finned fish foamy-like viruses, indicating that multiple episodes of viral infiltrations had occurred in this species. These findings fill a major gap in the Spumaretrovirinae taxon and reveal the aquatic origin of FVs found in terrestrial vertebrates. IMPORTANCE Although foamy viruses (FVs) have been found in major branches of vertebrates, the presence of these viruses in cartilaginous fish, the most basal living vertebrates with jaws, remains to be explored. This study revealed a collection of cartilaginous endogenous FVs in sharks and rays through in silico genomic mining. These viruses were rooted in the polymerase (POL) phylogeny, indicating the ancient aquatic origin of FVs. However, their envelope (ENV) protein grouped with those of amphibian FVs, suggesting different evolutionary histories of different FV genes. Overall, we provide the last missing gap for the taxonomic investigation of Spumaretrovirinae and provide concrete support for the aquatic origin of FVs.

How Well Does Evolution Explain Endogenous Retroviruses?-A Lakatosian Assessment

One of the most sophisticated philosophies of science is the methodology of scientific research programmes (MSRP), developed by Imre Lakatos. According to MSRP, scientists are working within so-called research programmes, consisting of a hard core of fixed convictions and a flexible protective belt of auxiliary hypotheses. Anomalies are accommodated by changes to the protective belt that do not affect the hard core. Under MSRP, research programmes are appraised as 'progressive' if they successfully predict novel facts but are judged as 'degenerative' if they merely offer ad hoc solutions to anomalies. This paper applies these criteria to the evolutionary research programme as it has performed during half a century of ERV research. It describes the early history of the field and the emergence of the endogenization-amplification theory on the origins of retroviral-like sequences. It then discusses various predictions and postdictions that were generated by the programme, regarding orthologous ERVs in different species, the presence of target site duplications and the divergence of long terminal repeats, and appraises how the programme has dealt with data that did not conform to initial expectations. It is concluded that the evolutionary research programme has been progressive with regard to the issues here examined.

Recombination marks the evolutionary dynamics of a recently endogenized retrovirus

All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. Although endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term coevolution of ERV and host because germline infections are rare and expansion is slow, and because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. We previously established the locus-specific distribution of cervid ERV (CrERV) in populations of mule deer. In this study, we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short- and long-insert mate pair reads and CrERV sequence generated at each locus. We report that CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting interlineage recombinants then endogenize and subsequently expand. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population.

Multiple Infiltration and Cross-Species Transmission of Foamy Viruses across the Paleozoic to the Cenozoic Era

Foamy viruses (FVs) are complex retroviruses that can infect humans and other animals. In this study, by integrating transcriptomic and genomic data, we discovered 412 FVs from 6 lineages in amphibians, which significantly increased the known set of FVs in amphibians. Among these lineages, salamander FVs maintained a coevolutionary pattern with their hosts that could be dated back to the Paleozoic era, while in contrast, frog FVs were much more likely acquired from cross-species (class-level) transmission in the Cenozoic era. In addition, we found that three distinct FV lineages had integrated into the genome of a salamander. Unexpectedly, we identified a lineage of endogenous FVs in caecilians that expressed all complete major genes, demonstrating the potential existence of an exogenous form of FV outside of mammals. Our discovery of rare phenomena in amphibian FVs has significantly increased our understanding of the macroevolution of the complex retrovirus. IMPORTANCE Foamy viruses (FVs) represent, more so than other viruses, the best model of coevolution between a virus and a host. This study represents the largest investigation so far of amphibian FVs and reveals 412 FVs of 6 distinct lineages from three major orders of amphibians. Besides a coevolutionary pattern, cross-species and repeated infections were also observed during the evolution of amphibian FVs. Remarkably, expressed FVs including a potential exogenous form were discovered, suggesting that active FVs might be underestimated in nature. These findings revealed that the multiple origins and complex evolution of amphibian FVs started from the Paleozoic era.

Gene Conversion amongst Alu SINE Elements

The process of non-allelic gene conversion acts on homologous sequences during recombination, replacing parts of one with the other to make them uniform. Such concerted evolution is best described as paralogous ribosomal RNA gene unification that serves to preserve the essential house-keeping functions of the converted genes. Transposed elements (TE), especially Alu short interspersed elements (SINE) that have more than a million copies in primate genomes, are a significant source of homologous units and a verified target of gene conversion. The consequences of such a recombination-based process are diverse, including multiplications of functional TE internal binding domains and, for evolutionists, confusing divergent annotations of orthologous transposable elements in related species. We systematically extracted and compared 68,097 Alu insertions in various primates looking for potential events of TE gene conversion and discovered 98 clear cases of Alu-Alu gene conversion, including 64 cases for which the direction of conversion was identified (e.g., AluS conversion to AluY). Gene conversion also does not necessarily affect the entire homologous sequence, and we detected 69 cases of partial gene conversion that resulted in virtual hybrids of two elements. Phylogenetic screening of gene-converted Alus revealed three clear hotspots of the process in the ancestors of Catarrhini, Hominoidea, and gibbons. In general, our systematic screening of orthologous primate loci for gene-converted TEs provides a new strategy and view of a post-integrative process that changes the identities of such elements.

Full-length LTR retroelements in Capsicum annuum revealed a few species-specific family bursts with insertional preferences

Capsicum annuum is a species that has undergone an expansion of the size of its genome caused mainly by the amplification of repetitive DNA sequences, including mobile genetic elements. Based on information obtained from sequencing the genome of pepper, the estimated fraction of retroelements is approximately 81%, and previous results revealed an important contribution of lineages derived from Gypsy superfamily. However, the dynamics of the retroelements in the C. annuum genome is poorly understood. In this way, the present work seeks to investigate the phylogenetic diversity and genomic abundance of the families of autonomous (complete and intact) LTR retroelements from C. annuum and inspect their distribution along its chromosomes. In total, we identified 1151 structurally full-length retroelements (340 Copia; 811 Gypsy) grouped in 124 phylogenetic families in the base of their retrotranscriptase. All the evolutive lineages of LTR retroelements identified in plants were present in pepper; however, three of them comprise 83% of the entire LTR retroelements population, the lineages Athila, Del/Tekay, and Ale/Retrofit. From them, only three families represent 70.8% of the total number of the identified retroelements. A massive family-specific wave of amplification of two of them occurred in the last 0.5 Mya (GypsyCa_16; CopiaCa_01), whereas the third is more ancient and occurred 3.0 Mya (GypsyCa_13). Fluorescent in situ hybridization performed with family and lineage-specific probes revealed contrasting patterns of chromosomal affinity. Our results provide a database of the populations LTR retroelements specific to C. annuum genome. The most abundant families were analyzed according to chromosome insertional preferences, suppling useful tools to the design of retroelement-based markers specific to the species.

Genomic analyses reveal evolutionary and geologic context for the plateau fungus Ophiocordyceps sinensis

Background

Ophiocordyceps sinensis, which is only naturally found in the high-elevation extreme environment of the Tibetan Plateau, has been used in traditional Chinese medicine. Information concerning the evolutionary and geologic context of O. sinensis remains limited, however.

Methods

We constructed the high-quality genome of O. sinensis and provided insight into the evolution and ecology of O. sinensis using comparative genomics.

Results

We mapped the whole genome of the anamorph/asexual form Hirsutella of O. sinensis using Illumina and PacBio sequencing technologies and obtained a well assembled genome of 119.2 Mbp size. Long-read Single Molecule Real Time (SMRT) sequencing technology generated an assembly with more accurate representation of repeat sequence abundances and placement. Evolutionary analyses indicated that O. sinensis diverged from other fungi 65.9 Mya in the Upper Cretaceous, during the uplift of the Tibetan Plateau. Gene family expansions and contractions in addition to genome inflation via long terminal repeat (LTR) retrotransposon insertions were implicated as an important driver of O. sinensis divergence. The insertion rate of LTR sequences into the O. sinensis genome peaked ~ 30-40 Mya, when the Tibetan Plateau rose rapidly. Gene Ontology (GO) enrichment analysis suggested that O. sinensis contained more genes related to ice binding compared to other closely related fungi, which may aid in their adaptability to the cold Tibetan Plateau. Further, heavy metal resistance genes were in low abundance in the O. sinensis genome, which may help to explain previous observations that O. sinensis tissues contain high levels of heavy metals.

Conclusions

Our results reveal the evolutionary, geological, and ecological context for the evolution of the O. sinensis genome and the factors that have contributed to the environmental adaptability of this valuable fungus. These findings suggest that genome inflation via LTR retrotransposon insertions in O. sinensis coincided with the uplift of the Tibetan Plateau. LTRs and the specific genetic mechanisms of O. sinensis contributed to its adaptation to the environment on the plateau.

Discovery and population genomics of structural variation in a songbird genus

Structural variation (SV) constitutes an important type of genetic mutations providing the raw material for evolution. Here, we uncover the genome-wide spectrum of intra- and interspecific SV segregating in natural populations of seven songbird species in the genus Corvus. Combining short-read (N = 127) and long-read re-sequencing (N = 31), as well as optical mapping (N = 16), we apply both assembly- and read mapping approaches to detect SV and characterize a total of 220,452 insertions, deletions and inversions. We exploit sampling across wide phylogenetic timescales to validate SV genotypes and assess the contribution of SV to evolutionary processes in an avian model of incipient speciation. We reveal an evolutionary young (~530,000 years) cis-acting 2.25-kb LTR retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth and evolutionary significance of SV segregating in natural populations and highlight the need for reliable SV genotyping.

Origins and evolutionary consequences of ancient endogenous retroviruses

Retroviruses infect a broad range of vertebrate hosts that includes amphibians, reptiles, fish, birds and mammals. In addition, a typical vertebrate genome contains thousands of loci composed of ancient retroviral sequences known as endogenous retroviruses (ERVs). ERVs are molecular remnants of ancient retroviruses and proof that the ongoing relationship between retroviruses and their vertebrate hosts began hundreds of millions of years ago. The long-term impact of retroviruses on vertebrate evolution is twofold: first, as with other viruses, retroviruses act as agents of selection, driving the evolution of host genes that block viral infection or that mitigate pathogenesis, and second, through the phenomenon of endogenization, retroviruses contribute an abundance of genetic novelty to host genomes, including unique protein-coding genes and cis-acting regulatory elements. This Review describes ERV origins, their diversity and their relationships to retroviruses and discusses the potential for ERVs to reveal virus-host interactions on evolutionary timescales. It also describes some of the many examples of cellular functions, including protein-coding genes and regulatory elements, that have evolved from ERVs.

The Role of Gene Conversion between Transposable Elements in Rewiring Regulatory Networks

Nature has found many ways to utilize transposable elements (TEs) throughout evolution. Many molecular and cellular processes depend on DNA-binding proteins recognizing hundreds or thousands of similar DNA motifs dispersed throughout the genome that are often provided by TEs. It has been suggested that TEs play an important role in the evolution of such systems, in particular, the rewiring of gene regulatory networks. One mechanism that can further enhance the rewiring of regulatory networks is nonallelic gene conversion between copies of TEs. Here, we will first review evidence for nonallelic gene conversion in TEs. Then, we will illustrate the benefits nonallelic gene conversion provides in rewiring regulatory networks. For instance, nonallelic gene conversion between TE copies offers an alternative mechanism to spread beneficial mutations that improve the network, it allows multiple mutations to be combined and transferred together, and it allows natural selection to work efficiently in spreading beneficial mutations and removing disadvantageous mutations. Future studies examining the role of nonallelic gene conversion in the evolution of TEs should help us to better understand how TEs have contributed to evolution.

What Can Long Terminal Repeats Tell Us About the Age of LTR Retrotransposons, Gene Conversion and Ectopic Recombination?

LTR retrotransposons constitute a significant part of plant genomes and their evolutionary dynamics play an important role in genome size changes. Current methods of LTR retrotransposon age estimation are based only on LTR (long terminal repeat) divergence. This has prompted us to analyze sequence similarity of LTRs in 25,144 LTR retrotransposons from fifteen plant species as well as formation of solo LTRs. We found that approximately one fourth of nested retrotransposons showed a higher LTR divergence than the pre-existing retrotransposons into which they had been inserted. Moreover, LTR similarity was correlated with LTR length. We propose that gene conversion can contribute to this phenomenon. Gene conversion prediction in LTRs showed potential converted regions in 25% of LTR pairs. Gene conversion was higher in species with smaller genomes while the proportion of solo LTRs did not change with genome size in analyzed species. The negative correlation between the extent of gene conversion and the abundance of solo LTRs suggests interference between gene conversion and ectopic recombination. Since such phenomena limit the traditional methods of LTR retrotransposon age estimation, we recommend an improved approach based on the exclusion of regions affected by gene conversion.

Identification and evolution of avian endogenous foamy viruses

A history of long-term co-divergence means that foamy viruses (family Retroviridae) provide an ideal framework to understanding virus-host evolution over extended time periods. Endogenous foamy viruses (EndFVs) are rare, and to date have only been described in a limited number of mammals, amphibians, reptiles and fish genomes. By screening 414 avian genomes we identified EndFVs in two bird species: the Maguari Stork (Ciconia maguari) and the Oriental Stork (Ciconia boyciana). Analyses of phylogenetic relationships, genome structures and flanking sequences revealed a single origin of EndFVs in Ciconia species. In addition, the marked incongruence between the virus and host phylogenies suggested that this integration event occurred independently in birds. In sum, by providing evidence that birds can be infected with foamy viruses, we fill the last major gap in the taxonomic distribution of foamy viruses and their animal hosts.

Transcriptionally active HERV-H retrotransposons demarcate topologically associating domains in human pluripotent stem cells

Chromatin architecture has been implicated in cell type-specific gene regulatory programs, yet how chromatin remodels during development remains to be fully elucidated. Here, by interrogating chromatin reorganization during human pluripotent stem cell (hPSC) differentiation, we discover a role for the primate-specific endogenous retrotransposon human endogenous retrovirus subfamily H (HERV-H) in creating topologically associating domains (TADs) in hPSCs. Deleting these HERV-H elements eliminates their corresponding TAD boundaries and reduces the transcription of upstream genes, while de novo insertion of HERV-H elements can introduce new TAD boundaries. The ability of HERV-H to create TAD boundaries depends on high transcription, as transcriptional repression of HERV-H elements prevents the formation of boundaries. This ability is not limited to hPSCs, as these actively transcribed HERV-H elements and their corresponding TAD boundaries also appear in pluripotent stem cells from other hominids but not in more distantly related species lacking HERV-H elements. Overall, our results provide direct evidence for retrotransposons in actively shaping cell type- and species-specific chromatin architecture.

A reptilian endogenous foamy virus sheds light on the early evolution of retroviruses

Endogenous retroviruses (ERVs) represent host genomic 'fossils' of ancient viruses. Foamy viruses, including those that form endogenous copies, provide strong evidence for virus-host co-divergence across the vertebrate phylogeny. Endogenous foamy viruses (EFVs) have previously been discovered in mammals, amphibians, and fish. Here we report a novel endogenous foamy virus, termed ERV-Spuma-Spu, in genome of the tuatara (Sphenodon punctatus), an endangered reptile species endemic to New Zealand. Phylogenetic analyses revealed that foamy viruses have likely co-diverged with their hosts over many millions of years. The discovery of ERV-Spuma-Spu fills a major gap in the fossil record of foamy viruses and provides important insights into the early evolution of retroviruses.

Gorillas have been infected with the HERV-K (HML-2) endogenous retrovirus much more recently than humans and chimpanzees

Human endogenous retrovirus-K (HERV-K) human mouse mammary tumor virus-like 2 (HML-2) is the most recently active endogenous retrovirus group in humans, and the only group with human-specific proviruses. HML-2 expression is associated with cancer and other diseases, but extensive searches have failed to reveal any replication-competent proviruses in humans. However, HML-2 proviruses are found throughout the catarrhine primates, and it is possible that they continue to infect some species today. To investigate this possibility, we searched for gorilla-specific HML-2 elements using both in silico data mining and targeted deep-sequencing approaches. We identified 150 gorilla-specific integrations, including 31 2-LTR proviruses. Many of these proviruses have identical LTRs, and are insertionally polymorphic, consistent with very recent integration. One identified provirus has full-length ORFs for all genes, and thus could potentially be replication-competent. We suggest that gorillas may still harbor infectious HML-2 virus and could serve as a model for understanding retrovirus evolution and pathogenesis in humans.

Variation in proviral content among human genomes mediated by LTR recombination

Background

Human endogenous retroviruses (HERVs) occupy a substantial fraction of the genome and impact cellular function with both beneficial and deleterious consequences. The vast majority of HERV sequences descend from ancient retroviral families no longer capable of infection or genomic propagation. In fact, most are no longer represented by full-length proviruses but by solitary long terminal repeats (solo LTRs) that arose via non-allelic recombination events between the two LTRs of a proviral insertion. Because LTR-LTR recombination events may occur long after proviral insertion but are challenging to detect in resequencing data, we hypothesize that this mechanism is a source of genomic variation in the human population that remains vastly underestimated.

Results

We developed a computational pipeline specifically designed to capture dimorphic proviral/solo HERV allelic variants from short-read genome sequencing data. When applied to 279 individuals sequenced as part of the Simons Genome Diversity Project, the pipeline retrieves most of the dimorphic loci previously reported for the HERV-K(HML2) subfamily as well as dozens of additional candidates, including members of the HERV-H and HERV-W families previously involved in human development and disease. We experimentally validate several of these newly discovered dimorphisms, including the first reported instance of an unfixed HERV-W provirus and an HERV-H locus driving a transcript (ESRG) implicated in the maintenance of embryonic stem cell pluripotency.

Conclusions

Our findings indicate that human proviral content exhibit more extensive interindividual variation than previously recognized, which has important bearings for deciphering the contribution of HERVs to human physiology and disease. Because LTR retroelements and LTR recombination are ubiquitous in eukaryotes, our computational pipeline should facilitate the mapping of this type of genomic variation for a wide range of organisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0142-3) contains supplementary material, which is available to authorized users.

Characterization of a heat-activated retrotransposon in Vigna angularis

In plants, several transposable elements are conserved across species. We found a homolog of ONSEN, which is a heat-activated retrotransposon originally isolated from Arabidopsis thaliana, in Vigna. The ONSEN-like elements (VaONS) were detected in all the analyzed Japanese accessions of Vigna angularis (adzuki bean) by Southern blot analysis. However, VaONS sequences were observed to be polymorphic in the different accessions. Interestingly, extrachromosomal DNA (ecDNA) was detected in some accessions of adzuki bean, indicating the conserved heat-activation of VaONS. Furthermore, we successfully induced retrotransposition of VaONS in adzuki plant regenerated through callus. Findings of our study should provide a new tool for molecular breeding of adzuki bean.

Origin of the retroviruses: when, where, and how?

•

Retroviruses are a virus family of considerable medical and veterinary importance.

•

Until recently, very little was known about deep retroviral origins.

•

New research supports a marine origin of retroviruses, ∼460–550 million years ago.

•

The evolutionary events leading to the origin of retroviruses remain obscure.

•

Improved understanding of Metaviridae diversity and evolution are required for this.

Retroviruses are a virus family of considerable medical and veterinary importance. Additionally, it is now clear that endogenous retroviruses (ERVs) comprise significant portions of vertebrate genomes. Until recently, very little was known about the deep evolutionary origins of retroviruses. However, advances in genomics and bioinformatics have opened the way for great strides in understanding. Recent research employing a wide variety of bioinformatic approaches has demonstrated that retroviruses evolved during the early Palaeozoic Era, between 460 and 550 million years ago, providing the oldest inferred date estimate for any virus group. This finding presents an important framework to investigate the evolutionary transitions that led to the emergence of the retroviruses, offering potential insights into the infectious origins of a major group of pathogenic viruses.

Discovery of an endogenous Deltaretrovirus in the genome of long-fingered bats (Chiroptera: Miniopteridae)

Retroviruses can create endogenous forms on infiltration into the germline cells of their hosts. These forms are then vertically transmitted and can be considered as genetic fossils of ancient viruses. All retrovirus genera, with the exception of deltaretroviruses, have had their representation identified in the host genome as a virus fossil record. Here we describe an endogenous Deltaretrovirus, identified in the germline of long-fingered bats (Miniopteridae). A single, heavily deleted copy of this retrovirus has been found in the genome of miniopterid species, but not in the genomes of the phylogenetically closest bat families, Vespertilionidae and Cistugonidae. Therefore, the endogenization occurred in a time interval between 20 and 45 million years ago. This discovery closes the last major gap in the retroviral fossil record and provides important insights into the history of deltaretroviruses in mammals.

Marine origin of retroviruses in the early Palaeozoic Era

Very little is known about the ancient origin of retroviruses, but owing to the discovery of their ancient endogenous viral counterparts, their early history is beginning to unfold. Here we report 36 lineages of basal amphibian and fish foamy-like endogenous retroviruses (FLERVs). Phylogenetic analyses reveal that ray-finned fish FLERVs exhibit an overall co-speciation pattern with their hosts, while amphibian FLERVs might not. We also observe several possible ancient viral cross-class transmissions, involving lobe-finned fish, shark and frog FLERVs. Sequence examination and analyses reveal two major lineages of ray-finned fish FLERVs, one of which had gained two novel accessory genes within their extraordinarily large genomes. Our phylogenetic analyses suggest that this major retroviral lineage, and therefore retroviruses as a whole, have an ancient marine origin and originated together with, if not before, their jawed vertebrate hosts >450 million years ago in the Ordovician period, early Palaeozoic Era.

Endogenous retroviruses are viruses that have become integrated into the genomes of their hosts. Here, the authors investigate the evolution of foamy-like endogenous retroviruses, and, by taking into account the temporal dynamics of the rate of viral evolution, suggest that retroviruses arose at least 450 million years ago in marine vertebrates.

Impact of repetitive DNA on sex chromosome evolution in plants

Structurally and functionally diverged sex chromosomes have evolved in many animals as well as in some plants. Sex chromosomes represent a specific genomic region(s) with locally suppressed recombination. As a consequence, repetitive sequences involving transposable elements, tandem repeats (satellites and microsatellites), and organellar DNA accumulate on the Y (W) chromosomes. In this paper, we review the main types of repetitive elements, their gathering on the Y chromosome, and discuss new findings showing that not only accumulation of various repeats in non-recombining regions but also opposite processes form Y chromosome. The aim of this review is also to discuss the mechanisms of repetitive DNA spread involving (retro) transposition, DNA polymerase slippage or unequal crossing-over, as well as modes of repeat removal by ectopic recombination. The intensity of these processes differs in non-recombining region(s) of sex chromosomes when compared to the recombining parts of genome. We also speculate about the relationship between heterochromatinization and the formation of heteromorphic sex chromosomes.

Evidence of extensive non-allelic gene conversion among LTR elements in the human genome

Long Terminal Repeats (LTRs) are nearly identical DNA sequences found at either end of Human Endogenous Retroviruses (HERVs). The high sequence similarity that exists among different LTRs suggests they could be substrate of ectopic gene conversion events. To understand the extent to which gene conversion occurs and to gain new insights into the evolutionary history of these elements in humans, we performed an intra-species phylogenetic study of 52 LTRs on different unrelated Y chromosomes. From this analysis, we obtained direct evidence that demonstrates the occurrence of ectopic gene conversion in several LTRs, with donor sequences located on both sex chromosomes and autosomes. We also found that some of these elements are characterized by an extremely high density of polymorphisms, showing one of the highest nucleotide diversities in the human genome, as well as a complex patchwork of sequences derived from different LTRs. Finally, we highlighted the limits of current short-read NGS studies in the analysis of genetic diversity of the LTRs in the human genome. In conclusion, our comparative re-sequencing analysis revealed that ectopic gene conversion is a common event in the evolution of LTR elements, suggesting complex genetic links among LTRs from different chromosomes.

Evolutionary Relationships among Extinct and Extant Sloths: The Evidence of Mitogenomes and Retroviruses

Macroevolutionary trends exhibited by retroviruses are complex and not entirely understood. The sloth endogenized foamy-like retrovirus (SloEFV), which demonstrates incongruence in virus–host evolution among extant sloths (Order Folivora), has not been investigated heretofore in any extinct sloth lineages and its premodern history within folivorans is therefore unknown. Determining retroviral coevolutionary trends requires a robust phylogeny of the viral host, but the highly reduced modern sloth fauna (6 species in 2 genera) does not adequately represent what was once a highly diversified clade (∼100 genera) of placental mammals. At present, the amount of molecular data available for extinct sloth taxa is limited, and analytical results based on these data tend to conflict with phylogenetic inferences made on the basis of morphological studies. To augment the molecular data set, we applied hybridization capture and next-generation Illumina sequencing to two extinct and three extant sloth species to retrieve full mitochondrial genomes (mitogenomes) from the hosts and the polymerase gene of SloEFV. The results produced a fully resolved and well-supported phylogeny that supports dividing crown families into two major clades: 1) The three-toed sloth, Bradypus, and Nothrotheriidae and 2) Megalonychidae, including the two-toed sloth, Choloepus, and Mylodontidae. Our calibrated time tree indicates that the Miocene epoch (23.5 Ma), particularly its earlier part, was an important interval for folivoran diversification. Both extant and extinct sloths demonstrate multiple complex invasions of SloEFV into the ancestral sloth germline followed by subsequent introgressions across different sloth lineages. Thus, sloth mitogenome and SloEFV evolution occurred separately and in parallel among sloths.

Cross-Species Transmission and Differential Fate of an Endogenous Retrovirus in Three Mammal Lineages

Endogenous retroviruses (ERVs) arise from retroviruses chromosomally integrated in the host germline. ERVs are common in vertebrate genomes and provide a valuable fossil record of past retroviral infections to investigate the biology and evolution of retroviruses over a deep time scale, including cross-species transmission events. Here we took advantage of a catalog of ERVs we recently produced for the bat Myotis lucifugus to seek evidence for infiltration of these retroviruses in other mammalian species (>100) currently represented in the genome sequence database. We provide multiple lines of evidence for the cross-ordinal transmission of a gammaretrovirus endogenized independently in the lineages of vespertilionid bats, felid cats and pangolin ~13-25 million years ago. Following its initial introduction, the ERV amplified extensively in parallel in both bat and cat lineages, generating hundreds of species-specific insertions throughout evolution. However, despite being derived from the same viral species, phylogenetic and selection analyses suggest that the ERV experienced different amplification dynamics in the two mammalian lineages. In the cat lineage, the ERV appears to have expanded primarily by retrotransposition of a single proviral progenitor that lost infectious capacity shortly after endogenization. In the bat lineage, the ERV followed a more complex path of germline invasion characterized by both retrotransposition and multiple infection events. The results also suggest that some of the bat ERVs have maintained infectious capacity for extended period of time and may be still infectious today. This study provides one of the most rigorously documented cases of cross-ordinal transmission of a mammalian retrovirus. It also illustrates how the same retrovirus species has transitioned multiple times from an infectious pathogen to a genomic parasite (i.e. retrotransposon), yet experiencing different invasion dynamics in different mammalian hosts.

A Novel Terminal-Repeat Retrotransposon in Miniature (TRIM) Is Massively Expressed in Echinococcus multilocularis Stem Cells

Taeniid cestodes (including the human parasites Echinococcus spp. and Taenia solium) have very few mobile genetic elements (MGEs) in their genome, despite lacking a canonical PIWI pathway. The MGEs of these parasites are virtually unexplored, and nothing is known about their expression and silencing. In this work, we report the discovery of a novel family of small nonautonomous long terminal repeat retrotransposons (also known as terminal-repeat retrotransposons in miniature, TRIMs) which we have named ta-TRIM (taeniid TRIM). ta-TRIMs are only the second family of TRIM elements discovered in animals, and are likely the result of convergent reductive evolution in different taxonomic groups. These elements originated at the base of the taeniid tree and have expanded during taeniid diversification, including after the divergence of closely related species such as Echinococcus multilocularis and Echinococcus granulosus. They are massively expressed in larval stages, from a small proportion of full-length copies and from isolated terminal repeats that show transcriptional read-through into downstream regions, generating novel noncoding RNAs and transcriptional fusions to coding genes. In E. multilocularis, ta-TRIMs are specifically expressed in the germinative cells (the somatic stem cells) during asexual reproduction of metacestode larvae. This would provide a developmental mechanism for insertion of ta-TRIMs into cells that will eventually generate the adult germ line. Future studies of active and inactive ta-TRIM elements could give the first clues on MGE silencing mechanisms in cestodes.

Identification and characterization of expressed retrotransposons in the genome of the Paracoccidioides species complex

BACKGROUND

Species from the Paracoccidioides complex are thermally dimorphic fungi and the causative agents of paracoccidioidomycosis, a deep fungal infection that is the most prevalent systemic mycosis in Latin America and represents the most important cause of death in immunocompetent individuals with systemic mycosis in Brazil. We previously described the identification of eight new families of DNA transposons in Paracoccidioides genomes. In this work, we aimed to identify potentially active retrotransposons in Paracoccidioides genomes.

RESULTS

We identified five different retrotransposon families (four LTR-like and one LINE-like element) in the genomes of three Paracoccidioides isolates. Retrotransposons were present in all of the genomes analyzed. P. brasiliensis and P. lutzii species harbored the same retrotransposon lineages but differed in their copy numbers. In the Pb01, Pb03 and Pb18 genomes, the number of LTR retrotransposons was higher than the number of LINE-like elements, and the LINE-like element RtPc5 was transcribed in Paracoccidioides lutzii (Pb01) but could not be detected in P. brasiliensis (Pb03 and Pb18) by semi-quantitative RT-PCR.

CONCLUSION

Five new potentially active retrotransposons have been identified in the genomic assemblies of the Paracoccidioides species complex using a combined computational and experimental approach. The distribution across the two known species, P. brasiliensis and P. lutzii, and phylogenetics analysis indicate that these elements could have been acquired before speciation occurred. The presence of active retrotransposons in the genome may have implications regarding the evolution and genetic diversification of the Paracoccidioides genus.

A novel endogenous betaretrovirus in the common vampire bat (Desmodus rotundus) suggests multiple independent infection and cross-species transmission events

The Desmodus rotundus endogenous betaretrovirus (DrERV) is fixed in the vampire bat D. rotundus population and in other phyllostomid bats but is not present in all species from this family. DrERV is not phylogenetically related to Old World bat betaretroviruses but to betaretroviruses from rodents and New World primates, suggesting recent cross-species transmission. A recent integration age estimation of the provirus in some taxa indicates that an exogenous counterpart might have been in recent circulation.

The decline of human endogenous retroviruses: extinction and survival

Background

Endogenous Retroviruses (ERVs) are retroviruses that over the course of evolution have integrated into germline cells and eventually become part of the host genome. They proliferate within the germline of their host, making up ~5% of the human and mouse genome sequences. Several lines of evidence have suggested a decline in the rate of ERV integration into the human genome in recent evolutionary history but this has not been investigated quantitatively or possible causes explored.

Results

By dating the integration of ERV loci in 40 mammal species, we show that the human genome and that of other hominoids (great apes and gibbons) have experienced an approximately four-fold decline in the ERV integration rate over the last 10 million years. A major cause is the recent extinction of one very large ERV lineage (HERV-H), which is responsible for most of the integrations over the last 30 million years. The decline however affects most other ERV lineages. Only about 10% of the decline might be attributed to an accompanying increase in body mass (a trait we have shown recently to be negatively correlated with ERV integration rate). Humans are unusual compared to related species – Old World monkeys, great apes and gibbons – in (a) having not acquired any new ERV lineages during the last 30 million years and (b) the possession of an old ERV lineage that has continued to replicate up until at least the last few hundred thousand years – the potentially medically significant HERVK(HML2).

Conclusions

The human genome shares with the genome of other great apes and gibbons a recent decline in ERV integration that is not typical of other primates and mammals. The human genome differs from that of related species both in maintaining up until at least recently a replicating old ERV lineage and in not having acquired any new lineages. We speculate that the decline in ERV integration in the human genome has been exacerbated by a relatively low burden of horizontally-transmitted retroviruses and subsequent reduced risk of endogenization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0136-x) contains supplementary material, which is available to authorized users.

Widespread genome reorganization of an obligate virus mutualist

The family Polydnaviridae is of interest because it provides the best example of viruses that have evolved a mutualistic association with their animal hosts. Polydnaviruses in the genus Bracovirus are strictly associated with parasitoid wasps in the family Braconidae, and evolved ∼100 million years ago from a nudivirus. Each wasp species relies on its associated bracovirus to parasitize hosts, while each bracovirus relies on its wasp for vertical transmission. Prior studies establish that bracovirus genomes consist of proviral segments and nudivirus-like replication genes, but how these components are organized in the genomes of wasps is unknown. Here, we sequenced the genome of the wasp Microplitis demolitor to characterize the proviral genome of M. demolitor bracovirus (MdBV). Unlike nudiviruses, bracoviruses produce virions that package multiple circular, double-stranded DNAs. DNA segments packaged into MdBV virions resided in eight dispersed loci in the M. demolitor genome. Each proviral segment was bounded by homologous motifs that guide processing to form mature viral DNAs. Rapid evolution of proviral segments obscured homology between other bracovirus-carrying wasps and MdBV. However, some domains flanking MdBV proviral loci were shared with other species. All MdBV genes previously identified to encode proteins required for replication were identified. Some of these genes resided in a multigene cluster but others, including subunits of the RNA polymerase that transcribes structural genes and integrases that process proviral segments, were widely dispersed in the M. demolitor genome. Overall, our results indicate that genome dispersal is a key feature in the evolution of bracoviruses into mutualists.

Inter- and intraspecies phylogenetic analyses reveal extensive X-Y gene conversion in the evolution of gametologous sequences of human sex chromosomes

It has long been believed that the male-specific region of the human Y chromosome (MSY) is genetically independent from the X chromosome. This idea has been recently dismissed due to the discovery that X–Y gametologous gene conversion may occur. However, the pervasiveness of this molecular process in the evolution of sex chromosomes has yet to be exhaustively analyzed. In this study, we explored how pervasive X–Y gene conversion has been during the evolution of the youngest stratum of the human sex chromosomes. By comparing about 0.5 Mb of human–chimpanzee gametologous sequences, we identified 19 regions in which extensive gene conversion has occurred. From our analysis, two major features of these emerged: 1) Several of them are evolutionarily conserved between the two species and 2) almost all of the 19 hotspots overlap with regions where X–Y crossing-over has been previously reported to be involved in sex reversal. Furthermore, in order to explore the dynamics of X–Y gametologous conversion in recent human evolution, we resequenced these 19 hotspots in 68 widely divergent Y haplogroups and used publicly available single nucleotide polymorphism data for the X chromosome. We found that at least ten hotspots are still active in humans. Hence, the results of the interspecific analysis are consistent with the hypothesis of widespread reticulate evolution within gametologous sequences in the differentiation of hominini sex chromosomes. In turn, intraspecific analysis demonstrates that X–Y gene conversion may modulate human sex-chromosome-sequence evolution to a greater extent than previously thought.

Evolutionary history of bovine endogenous retroviruses in the Bovidae family

Background

Endogenous retroviruses (ERVs) are genomic elements of retroviral origin that are present in the genomes of almost all vertebrates. In cattle, more than 13,000 elements related to ERVs have been detected, and based on the pol gene, 24 families or groups of bovine ERVs have been described. However, information about ERVs in other bovids and the presence of families of related bovine ERVs in different species of the Bovidae family is scarce.

Results

The 24 families of bovine ERVs previously detected in cattle (Bos taurus) were also detected in zebus (Bos indicus) and yaks (Bos grunniens). In addition, six new families, named BoERV25 to BoERV30, were detected in the three Bos species. Five more ruminant species were screened for related ERVs: 26 families were detected in these species, but four families (BoERV24, BoERV26, BoERV28 and BoERV29) were specific to cattle, zebus, yaks and buffalo. An analysis of the homology of the ERVs of cattle, zebus and yaks revealed that the level of LTR divergence was similar between ERVs from cattle and zebus but was less similar between with ERVs from cattle and yaks. In addition, purifying selection was detected in the genes and retroviral regions of clusters of ERVs of cattle, zebus and yaks.

Conclusions

In this work, the 24 ERV families previously identified in cattle were also found in two other species in the Bos genus. In addition, six new bovine ERV families were detected. Based on LTR divergence, the most recently inserted families are from Class II. The divergence of the LTR, used as an indirect estimate of the ERV insertion time, seemed to be influenced by the differences in genome evolution since the divergence of the species. In addition, purifying selection could be acting on clusters of ERVs from different species.

Three distinct mutational mechanisms acting on a single gene underpin the origin of yellow flesh in peach

Peach flesh color (white or yellow) is among the most popular commercial criteria for peach classification, and has implications for consumer acceptance and fruit nutritional quality. Despite the increasing interest in improving cultivars of both flesh types, little is known about the genetic basis for the carotenoid content diversity in peach. Here we describe the association between genotypes at a locus encoding the carotenoid cleavage dioxygenase 4 (PpCCD4), localized in pseudomolecule 1 of the Prunus persica reference genome sequence, and the flesh color for 37 peach varieties, including two somatic revertants, and three ancestral relatives of peach, providing definitive evidence that this locus is responsible for flesh color phenotype. We show that yellow peach alleles have arisen from various ancestral haplotypes by at least three independent mutational events involving nucleotide substitutions, small insertions and transposable element insertions, and that these mutations, despite being located within the transcribed portion of the gene, also result in marked differences in transcript levels, presumably as a consequence of differential transcript stability involving nonsense-mediated mRNA decay. The PpCCD4 gene provides a unique example of a gene for which humans, in their quest to diversify phenotypic appearance and qualitative characteristics of a fruit, have been able to select and exploit multiple mutations resulting from a variety of mechanisms.

Identification of an ancient endogenous retrovirus, predating the divergence of the placental mammals

The evolutionary arms race between mammals and retroviruses has long been recognized as one of the oldest host-parasite interactions. Rapid evolution rates in exogenous retroviruses have often made accurate viral age estimations highly problematic. Endogenous retroviruses (ERVs), however, integrate into the germline of their hosts, and are subjected to their evolutionary rates. This study describes, for the first time, a retroviral orthologue predating the divergence of placental mammals, giving it a minimum age of 104-110 Myr. Simultaneously, other orthologous selfish genetic elements (SGEs), inserted into the ERV sequence, provide evidence for the oldest individual mammalian-wide interspersed repeat and medium-reiteration frequency interspersed repeat mammalian repeats, with the same minimum age. The combined use of shared SGEs and reconstruction of viral orthologies defines new limits and increases maximum 'lookback' times, with subsequent implications for the field of paleovirology.

Genome-wide characterization of endogenous retroviruses in the bat Myotis lucifugus reveals recent and diverse infections

Bats are increasingly recognized as reservoir species for a variety of zoonotic viruses that pose severe threats to human health. While many RNA viruses have been identified in bats, little is known about bat retroviruses. Endogenous retroviruses (ERVs) represent genomic fossils of past retroviral infections and, thus, can inform us on the diversity and history of retroviruses that have infected a species lineage. Here, we took advantage of the availability of a high-quality genome assembly for the little brown bat, Myotis lucifugus, to systematically identify and analyze ERVs in this species. We mined an initial set of 362 potentially complete proviruses from the three main classes of ERVs, which were further resolved into 13 major families and 86 subfamilies by phylogenetic analysis. Consensus or representative sequences for each of the 86 subfamilies were then merged to the Repbase collection of known ERV/long terminal repeat (LTR) elements to annotate the retroviral complement of the bat genome. The results show that nearly 5% of the genome assembly is occupied by ERV-derived sequences, a quantity comparable to findings for other eutherian mammals. About one-fourth of these sequences belong to subfamilies newly identified in this study. Using two independent methods, intraelement LTR divergence and analysis of orthologous loci in two other bat species, we found that the vast majority of the potentially complete proviruses identified in M. lucifugus were integrated in the last ~25 million years. All three major ERV classes include recently integrated proviruses, suggesting that a wide diversity of retroviruses is still circulating in Myotis bats.

Identification of diverse full-length endogenous betaretroviruses in megabats and microbats

BACKGROUND

Betaretroviruses infect a wide range of species including primates, rodents, ruminants, and marsupials. They exist in both endogenous and exogenous forms and are implicated in animal diseases such as lung cancer in sheep, and in human disease, with members of the human endogenous retrovirus-K (HERV-K) group of endogenous betaretroviruses (βERVs) associated with human cancers and autoimmune diseases. To improve our understanding of betaretroviruses in an evolutionarily distinct host species, we characterized βERVs present in the genomes and transcriptomes of mega- and microbats, which are an important reservoir of emerging viruses.

RESULTS

A diverse range of full-length βERVs were discovered in mega- and microbat genomes and transcriptomes including the first identified intact endogenous retrovirus in a bat. Our analysis revealed that the genus Betaretrovirus can be divided into eight distinct sub-groups with evidence of cross-species transmission. Betaretroviruses are revealed to be a complex retrovirus group, within which one sub-group has evolved from complex to simple genomic organization through the acquisition of an env gene from the genus Gammaretrovirus. Molecular dating suggests that bats have contended with betaretroviral infections for over 30 million years.

CONCLUSIONS

Our study reveals that a diverse range of betaretroviruses have circulated in bats for most of their evolutionary history, and cluster with extant betaretroviruses of divergent mammalian lineages suggesting that their distribution may be largely unrestricted by host species barriers. The presence of βERVs with the ability to transcribe active viral elements in a major animal reservoir for viral pathogens has potential implications for public health.

New insights into nested long terminal repeat retrotransposons in Brassica species

Long terminal repeat (LTR) retrotransposons, one of the foremost types of transposons, continually change or modify gene function and reorganize the genome through bursts of dramatic proliferation. Many LTR-TEs preferentially insert within other LTR-TEs, but the cause and evolutionary significance of these nested LTR-TEs are not well understood. In this study, a total of 1.52Gb of Brassica sequence containing 2020 bacterial artificial chromosomes (BACs) was scanned, and six bacterial artificial chromosome (BAC) clones with extremely nested LTR-TEs (LTR-TEs density: 7.24/kb) were selected for further analysis. The majority of the LTR-TEs in four of the six BACs were found to be derived from the rapid proliferation of retrotransposons originating within the BAC regions, with only a few LTR-TEs originating from the proliferation and insertion of retrotransposons from outside the BAC regions approximately 5-23Mya. LTR-TEs also preferably inserted into TA-rich repeat regions. Gene prediction by Genescan identified 207 genes in the 0.84Mb of total BAC sequences. Only a few genes (3/207) could be matched to the Brassica expressed sequence tag (EST) database, indicating that most genes were inactive after retrotransposon insertion. Five of the six BACs were putatively centromeric. Hence, nested LTR-TEs in centromere regions are rapidly duplicated, repeatedly inserted, and act to suppress activity of genes and to reshuffle the structure of the centromeric sequences. Our results suggest that LTR-TEs burst and proliferate on a local scale to create nested LTR-TE regions, and that these nested LTR-TEs play a role in the formation of centromeres.