Endogenous Retroviruses in Fish Genomes: From Relics of Past Infections to Evolutionary Innovations?

A reference genome for the Harpy Eagle reveals steady demographic decline and chromosomal rearrangements in the origin of Accipitriformes

The Harpy Eagle (Harpia harpyja) is an iconic species that inhabits forested landscapes in Neotropical regions, with decreasing population trends mainly due to habitat loss, and currently classified as vulnerable. Here, we report on a chromosome-scale genome assembly for a female individual combining long reads, optical mapping, and chromatin conformation capture reads. The final assembly spans 1.35 Gb, with N50scaffold equal to 58.1 Mb and BUSCO completeness of 99.7%. We built the first extensive transposable element (TE) library for the Accipitridae to date and identified 7,228 intact TEs. We found a burst of an unknown TE ~ 13-22 million years ago (MYA), coincident with the split of the Harpy Eagle from other Harpiinae eagles. We also report a burst of solo-LTRs and CR1 retrotransposons ~ 31-33 MYA, overlapping with the split of the ancestor to all Harpiinae from other Accipitridae subfamilies. Comparative genomics with other Accipitridae, the closely related Cathartidae and Galloanserae revealed major chromosome-level rearrangements at the basal Accipitriformes genome, in contrast to a conserved ancient genome architecture for the latter two groups. A historical demography reconstruction showed a rapid decline in effective population size over the last 20,000 years. This reference genome serves as a crucial resource for future conservation efforts towards the Harpy Eagle.

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

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

Exploring the role of endogenous retroviruses in seasonal reproductive cycles: a case study of the ERV-V envelope gene in mink

Introduction

Endogenous retroviruses (ERVs), which originated from exogenous retroviral infections of germline cells millions of years ago and were inherited by subsequent generations as per Mendelian inheritance patterns, predominantly comprise non-protein-coding sequences due to the accumulation of mutations, insertions, deletions, and truncations. Nevertheless, recent studies have revealed that ERVs play a crucial role in diverse biological processes by encoding various proteins.

Methods

In this study, we successfully identified an ERV envelope (env) gene in a mink species. A phylogenetic tree of mink ERV-V env and reference sequences was constructed using Bayesian methods and maximum-likelihood inference.

Results

Phylogenetic analyses indicated a significant degree of sequence conservation and positive selection within the env-surface open reading frame. Additionally, qRT-PCR revealed diverse patterns of mink ERV-V env expression in various tissues. The expression of mink ERV-V env gene in testicular tissue strongly correlated with the seasonal reproductive cycles of minks.

Discussion

Our study suggests that the ERV-V env gene in mink may have been repurposed for host functions.

Comparison of a Genotype 1 and a Genotype 2 Macaque Foamy Virus env Gene Indicates Distinct Infectivity and Cell-Cell Fusion but Similar Tropism and Restriction of Cell Entry by Interferon-Induced Transmembrane Proteins

Foamy viruses (FVs) are naturally found in many different animals and also in primates with the notable exception of humans, but zoonotic infections are common. In several species, two different envelope (env) gene sequence clades or genotypes exist. We constructed a simian FV (SFV) clone containing a reporter gene cassette. In this background, we compared the env genes of the SFVmmu-DPZ9524 (genotype 1) and of the SFVmmu_R289hybAGM (genotype 2) isolates. SFVmmu_R289hybAGM env-driven infection was largely resistant to neutralization by SFVmmu-DPZ9524-neutralizing sera. While SFVmmu_R289hybAGM env consistently effected higher infectivity and cell-cell fusion, we found no differences in the cell tropism conferred by either env across a range of different cells. Infection by both viruses was weakly and non-significantly enhanced by simultaneous knockout of interferon-induced transmembrane proteins (IFITMs) 1, 2, and 3 in A549 cells, irrespective of prior interferon stimulation. Infection was modestly reduced by recombinant overexpression of IFITM3, suggesting that the SFV entry step might be weakly restricted by IFITM3 under some conditions. Overall, our results suggest that the different env gene clades in macaque foamy viruses induce genotype-specific neutralizing antibodies without exhibiting overt differences in cell tropism, but individual env genes may differ significantly with regard to fitness.

Evolutionary analysis of endogenous intronic retroviruses in primates reveals an enrichment in transcription binding sites associated with key regulatory processes

Background

Endogenous retroviruses (ERVs) are the result of the integration of retroviruses into host DNA following germline infection. Endogenous retroviruses are made up of three main genes: gag, pol, and env, each of which encodes viral proteins that can be conserved or not. ERVs have been observed in a wide range of vertebrate genomes and their functions are associated with viral silencing and gene regulation.

Results

In this work, we studied the evolutionary history of endogenous retroviruses associated with five human genes (INPP5B, DET1, PSMA1, USH2A, and MACROD2), which are located within intron sections. To verify the retroviral origin of the candidates, several approaches were used to detect and locate ERV elements. Both orthologous and paralogous genes were identified by Ensembl and then analyzed for ERV presence using RetroTector. A phylogenetic tree was reconstructed to identify the minimum time point of ERV acquisition. From that search, we detected ERVs throughout the primate lineage and in some other groups. Also, we identified the minimum origin of the ERVs from the parvorder Catarrhini to the Homininae subfamily.

Conclusions

With the data collected, and by observing the transcription factors annotated inside ERVs, we propose that these elements play a relevant role in gene expression regulation and they probably possess important features for tumorigenesis control.

Cross-species transmission of an ancient endogenous retrovirus and convergent co-option of its envelope gene in two mammalian orders

Endogenous retroviruses (ERVs) found in vertebrate genomes are remnants of retroviral invasions of their ancestral species. ERVs thus represent molecular fossil records of ancient retroviruses and provide a unique opportunity to study viral-host interactions, including cross-species transmissions, in deep time. While most ERVs contain the mutated remains of the original retrovirus, on rare occasions evolutionary selection pressures lead to the co-option/exaptation of ERV genes for a host function. Here, we report the identification of two ancient related non-orthologous ERV env genes, ARTenvV and CARenvV, that are preserved with large open reading frames (ORFs) in the mammalian orders Artiodactyla and Carnivora, respectively, but are not found in other mammals. These Env proteins lack a transmembrane motif, but phylogenetic analyses show strong sequence preservation and positive selection of the env surface ORF in their respective orders, and transcriptomic analyses show a broad tissue expression pattern for both ARTenvV and CARenvV, suggesting that these genes may be exapted for a host function. Multiple lines of evidence indicate that ARTenvV and CARenvV were derived from an ancient ancestral exogenous gamma-like retrovirus that was independently endogenized in two mammalian orders more than 60 million years ago, which roughly coincides with the K-Pg mass extinction event and subsequent mammalian diversification. Thus, these findings identify the oldest known retroviral cross-ordinal transmission of a gamma-like retrovirus with no known extant infectious counterpart in mammals, and the first discovery of the convergent co-option of an ERV gene derived from the same ancestral retrovirus in two different mammalian orders.

Global network mapping research landscape and trends of the endogenous retroviruses: a look through bibliometric analysis

Endogenous retrovirus (ERV) research amalgamates host-retroviral coevolutionary, phylogenomic, infection, immunity, and cellular studies in various hosts ranging from fish to humans. Henceforth, a bibliometric analysis of these publications may aid in the identification of trends in ERV research. It was the foremost bibliographic study, with the key aim to conduct the bibliometric network analysis (e.g. co-authorship, co-occurrence, citation, bibliographic coupling, and co-citation analysis) to find the most prolific authors, organizations, and countries in ERV research, based on the mapping of bibliographic data. Second, the mapping based on text data comprised to chalk out the research trend over the time. The global literature about endogenous retroviruses published between 1985 and Sep 2021 was searched in the Web of Science (Core Collection) database using the “ENDOGENOUS RETROVIRUS” keyword. The bibliometric analysis of this dataset was carried out using VOSviewer version 1.6.17. According to findings, English was the de facto language of these publications, and 2157 were original articles. Among 2939 published documents, “endogenous retrovirus” was the most frequent keyword. Moreover, it revealed the United States as a core contributor to studies on the ERV. The Journal of Virology published a substantial amount of manuscripts in ERV. Robert Koch Institute and Harvard University were leading organizations for research in this field. The application of ERV research from China could be the research hotspot to follow in the coming years. Current bibliometric analysis provides a comprehensive picture of ERV research progress and has highlighted the contribution of different stakeholders.

Mapping the m1A, m5C, m6A and m7G methylation atlas in zebrafish brain under hypoxic conditions by MeRIP-seq

BACKGROUND

The epigenetic modifications play important regulatory roles in tissue development, maintenance of physiological functions and pathological process. RNA methylations, including newly identified m1A, m5C, m6A and m7G, are important epigenetic modifications. However, how these modifications are distributed in the transcriptome of vertebrate brains and whether their abundance is altered under pathological conditions are still poorly understood. In this study, we chose the model animal of zebrafish to conduct a systematic study to investigate the mRNA methylation atlas in the brain.

RESULTS

By performing unbiased analyses of the m1A, m5C, m6A and m7G methylation of mRNA, we found that within the whole brain transcriptome, with the increase of the gene expression levels, the overall level of each of these four modifications on the related genes was also progressively increased. Further bioinformatics analysis indicated that the zebrafish brain has an abundance of m1A modifications. In the hypoxia-treated zebrafish brains, the proportion of m1A is decreased, affecting the RNA splicing and zebrafish endogenous retroviruses.

CONCLUSIONS

Our study presents the first comprehensive atlas of m1A, m5C, m6A and m7G in the epitranscriptome of the zebrafish brain and reveals the distribution of these modifications in mRNA under hypoxic conditions. These data provide an invaluable resource for further research on the involvement of m1A, m5C, m6A and m7G in the regulation of miRNA and repeat elements in vertebrates, and provide new thoughts to study the brain hypoxic injury on the aspect of epitranscriptome.

Genome-Wide Characterization of Zebrafish Endogenous Retroviruses Reveals Unexpected Diversity in Genetic Organizations and Functional Potentials

Endogenous retroviruses (ERVs) occupy a substantial fraction of mammalian genomes. However, whether ERVs extensively exist in ancient vertebrates remains unexplored. Here, we performed a genome-wide characterization of ERVs in a zebrafish (Danio rerio) model. Approximately 3,315 ERV-like elements (DrERVs) were identified as Gypsy, Copia, Bel, and class I−III groups. DrERVs accounted for approximately 2.3% of zebrafish genome and were distributed in all 25 chromosomes, with a remarkable bias on chromosome 4. Gypsy and class I are the two most abundant groups with earlier insertion times. The vast majority of the DrERVs have varied structural defects. A total of 509 gag and 71 env genes with coding potentials were detected. The env-coding elements were well-characterized and classified into four subgroups. A ERV-E4.8.43-DanRer element shows high similarity with HERV9NC-int in humans and analogous sequences were detected in species spanning from fish to mammals. RNA-seq data showed that hundreds of DrERVs were expressed in embryos and tissues under physiological conditions, and most of them exhibited stage and tissue specificity. Additionally, 421 DrERVs showed strong responsiveness to virus infection. A unique group of DrERVs with immune-relevant genes, such as fga, ddx41, ftr35, igl1c3, and tbk1, instead of intrinsic viral genes were identified. These DrERVs are regulated by transcriptional factors binding at the long terminal repeats. This study provided a survey of the composition, phylogeny, and potential functions of ERVs in a fish model, which benefits the understanding of the evolutionary history of ERVs from fish to mammals.

IMPORTANCE Endogenous retroviruses (ERVs) are relics of past infection that constitute up to 8% of the human genome. Understanding the genetic evolution of the ERV family and the interplay of ERVs and encoded RNAs and proteins with host function has become a new frontier in biology. Fish, as the most primitive vertebrate host for retroviruses, is an indispensable integral part for such investigations. In the present study, we report the genome-wide characterization of ERVs in zebrafish, an attractive model organism of ancient vertebrates from multiple perspectives, including composition, genomic organization, chromosome distribution, classification, phylogeny, insertion time, characterization of gag and env genes, and expression profiles in embryos and tissues. The result helps uncover the evolutionarily conserved and fish-specific ERVs, as well as the immune-relevant ERVs in response to virus infection. This study demonstrates the previously unrecognized abundance, diversification, and extensive activity of ERVs at the early stage of ERV evolution.

Social Networking of Quasi-Species Consortia drive Virolution via Persistence

The emergence of cooperative quasi-species consortia (QS-C) thinking from the more accepted quasispecies equations of Manfred Eigen, provides a conceptual foundation from which concerted action of RNA agents can now be understood. As group membership becomes a basic criteria for the emergence of living systems, we also start to understand why the history and context of social RNA networks become crucial for survival and function. History and context of social RNA networks also lead to the emergence of a natural genetic code. Indeed, this QS-C thinking can also provide us with a transition point between the chemical world of RNA replicators and the living world of RNA agents that actively differentiate self from non-self and generate group identity with membership roles. Importantly the social force of a consortia to solve complex, multilevel problems also depend on using opposing and minority functions. The consortial action of social networks of RNA stem-loops subsequently lead to the evolution of cellular organisms representing a tree of life.

Unexpected Discovery and Expression of Amphibian Class II Endogenous Retroviruses

Endogenous retroviruses (ERVs) are the remnants of past retroviral infections. Fossil records of class II retroviruses have been discovered in a range of vertebrates, with the exception of amphibians, which are known only to possess class I and class III-like ERVs. Through genomic mining of all available amphibian genomes, we identified, for the first time, class II ERVs in amphibians. The class II ERVs were found only in Gymnophiona (caecilians) and not in the genomes of the other amphibian orders, Anura (frogs and toads) and Caudata (salamanders and newts), which are phylogenetically closely related. Therefore, the ERV endogenization occurred after the split of Gymnophiona, Anura, and Caudata (323 million years ago). Investigation of phylogenetic relationship and genomic structure revealed that the ERVs may originate from alpha- or betaretroviruses. We offer evidence that class II ERVs infiltrated amphibian genomes recently and may still have infectious members. Remarkably, certain amphibian class II ERVs can be expressed in diverse tissues. This discovery closes the major gap in the retroviral fossil record of class II ERVs and provides important insights into the evolution of class II ERVs in vertebrates.IMPORTANCE Class II retroviruses, largely distributed among mammals and birds, are of particular importance for medicine and economics. Class II ERVs have been discovered in a range of vertebrates, with the exception of amphibians, which are known only to possess class I and class III-like ERVs. Here, for the first time, we discovered class II ERVs in amphibians. We also revealed that the ERVs may originate from alpha- or betaretroviruses. We revealed that class II ERVs were integrated into amphibian genomes recently and certain amphibian class II ERVs can be expressed in diverse tissues. Our discovery closes the major gap in the retroviral fossil record of class II ERVs, and also indicates that amphibians may be still infected by class II retroviruses.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

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.

GGERV20, a recently integrated, segregating endogenous retrovirus in Gallus gallus

Endogenous retroviruses (ERVs) are widespread in vertebrate genomes. The recent availability of whole eukaryotic genomes has enabled their characterization in many organisms, including Gallus gallus (red jungle fowl), the progenitor of the domesticated chicken. Our bioinformatics analysis of a G. gallus ERV previously designated GGERV20 identified 35 proviruses with complete long terminal repeats (LTRs) and gag-pol open reading frames (ORFs) in the Genome Reference Consortium Chicken Build 6a, of which 8 showed potential for translation of functional retroviral polyproteins, including the integrase and reverse transcriptase enzymes. No elements were discovered with an env gene. Fifteen loci had LTR sequences with 100 % identity, indicative of recent integration. Chicken embryo fibroblast RNA-seq datasets showed reads representing the entire length of the GGERV20 provirus, supporting their potential for expressing viral proteins. To investigate the possibility that GGERV20 elements may not be fixed in the genome, we assessed the integration status of five loci in a meat-type chicken. PCRs targeting a GGERV20 locus on G. gallus chromosome one (GGERV20_1-1) reproducibly amplified both LTRs and the preintegration state, indicating that the bird from which the DNA was sampled was hemizygous at this locus. The four other loci examined only produced the preintegration state amplicons. These results reveal that GGERV20 is not fixed in the G. gallus population, and taken together with the lack of mutations seen in several provirus LTRs and their transcriptional activity, suggest that GGERV20 retroviruses have recently been and continue to be active in the chicken genome.

Viral Symbiosis in the Origins and Evolution of Life with a Particular Focus on the Placental Mammals

Advances in understanding over the last decade or so highlight the need for a reappraisal of the role of viruses in relation to the origins and evolution of cellular life, as well as in the homeostasis of the biosphere on which all of life depends. The relevant advances have, in particular, revealed an important contribution of viruses to the evolution of the placental mammals, while also contributing key roles to mammalian embryogenesis, genomic evolution, and physiology. Part of this reappraisal will include the origins of viruses, a redefinition of their quintessential nature, and a suggestion as to how we might view viruses in relation to the tree of life.

Endogenous retroviruses of non-avian/mammalian vertebrates illuminate diversity and deep history of retroviruses

The deep history and early diversification of retroviruses remains elusive, largely because few retroviruses have been characterized in vertebrates other than mammals and birds. Endogenous retroviruses (ERVs) documented past retroviral infections and thus provide ‘molecular fossils’ for studying the deep history of retroviruses. Here we perform a comprehensive phylogenomic analysis of ERVs within the genomes of 92 non-avian/mammalian vertebrates, including 72 fishes, 4 amphibians, and 16 reptiles. We find that ERVs are present in all the genomes of jawed vertebrates, revealing the ubiquitous presence of ERVs in jawed vertebrates. We identify a total of >8,000 ERVs and reconstruct ~450 complete or partial ERV genomes, which dramatically expands the phylogenetic diversity of retroviruses and suggests that the diversity of exogenous retroviruses might be much underestimated in non-avian/mammalian vertebrates. Phylogenetic analyses show that retroviruses cluster into five major groups with different host distributions, providing important insights into the classification and diversification of retroviruses. Moreover, we find retroviruses mainly underwent frequent host switches in non-avian/mammalian vertebrates, with exception of spumavirus-related viruses that codiverged with their ray-finned fish hosts. Interestingly, ray-finned fishes and turtles appear to serve as unappreciated hubs for the transmission of retroviruses. Finally, we find retroviruses underwent many independent water-land transmissions, indicating the water-land interface is not a strict barrier for retrovirus transmission. Our analyses provide unprecedented insights into and valuable resources for studying the diversification, key evolutionary transitions, and macroevolution of retroviruses.

Retroviruses infect a wide range of vertebrates and cause many diseases, such as AIDS and cancers. To date, retroviruses have been rarely characterized in vertebrates other than mammals and birds, impeding our understanding of the diversity and early evolution of retroviruses. Retroviruses can occasionally integrate into host genomes and become endogenous retroviruses (ERVs), which provide molecular fossils for studying the long-term evolution of retroviruses. Here we performed comparative genomic and evolutionary analyses of ERVs within 92 non-avian/mammalian vertebrates (fishes, amphibians, and reptiles) and uncovered extraordinary diversity of retroviruses in non-avian/mammalian vertebrates. Our analyses reveal an ancient aquatic origin of retroviruses and retroviruses underwent frequent host-switching. Our findings have important implications in understanding the deep history and evolutionary mode of retroviruses.