An Intact Retroviral Gene Conserved in Spiny-Rayed Fishes for over 100 My

The diverse evolutionary histories of domesticated metaviral capsid genes in mammals

Selfish genetic elements and their remnants comprise at least half of the human genome. Active transposons duplicate by inserting copies at new sites in a host genome. Following insertion, transposons can acquire mutations that render them inactive; the accrual of additional mutations can render them unrecognizable over time. However, in rare instances, segments of transposons become useful for the host, in a process called gene domestication. Using the first complete human genome assembly and 25 additional vertebrate genomes, we analyzed the evolutionary trajectories and functional potential of genes domesticated from the capsid genes of Metaviridae, a retroviral-like retrotransposon family. Our analysis reveals four families of domesticated capsid genes in placental mammals with varied evolutionary outcomes, ranging from universal retention to lineage-specific duplications or losses and from purifying selection to lineage-specific rapid evolution. The four families of domesticated capsid genes have divergent amino-terminal domains, inherited from four distinct ancestral metaviruses. Structural predictions reveal that many domesticated genes encode a previously unrecognized RNA-binding domain retained in multiple paralogs in mammalian genomes both adjacent to and independent from the capsid domain. Collectively, our study reveals diverse outcomes of domestication of diverse metaviruses, which led to structurally and evolutionarily diverse genes that encode important, but still largely-unknown functions in placental mammals. (207).

Evolution and Genetic Diversity of the Retroviral Envelope in Anamniotes

Retroviruses are widely distributed in all vertebrates, as are their endogenous forms, endogenous retroviruses (ERV), which serve as "fossil" evidence to trace the ancient origins and history of virus-host interactions over millions of years. The retroviral envelope (Env) plays a significant role in host range determination, but major information on their genetic diversification and evolution in anamniotes is lacking. Here, by incorporating multiple-round in silico similarity search and phylogenomic analysis, more than 30,000 copies of ERV lineages with gamma-type Env (GTE), covalently associated Env, were discovered by searching against all fish and amphibian genomes and transcriptomic assemblies, but no beta-type Env (BTE), noncovalently associated Env, was found. Furthermore, a nine-type classification system of anamniote GTE was proposed by combining phylogenetic and domain/motif analyses. The elastic genomic organization and overall phylogenetic incongruence between anamniotic Env and its neighboring polymerase (Pol) implied that early retroviral diversification in anamniotic vertebrates was facilitated by frequent recombination. At last, host cellular opioid growth factor receptor (OGFr) gene capturing by anamniotic ERVs with GTE was reported for the first time. Overall, our findings overturn traditional Pol genotyping and reveal a complex evolutionary history of anamniotic retroviruses inferred by Env evolution. IMPORTANCE Although the retroviral envelope (Env) protein in amniotes has been well studied, its evolutionary history in anamniotic vertebrates is ambiguous. By analyzing more than 30,000 copies of ERV lineages with gamma-type Env (GTE) in anamniotes, several important evolutionary features were identified. First, GTE was found to be widely distributed among different amphibians and fish. Second, nine types of GTE were discovered and defined, revealing their great genetic diversity. Third, the incongruence between the Env and Pol phylogenies suggested that frequent recombination shaped the early evolution of anamniote retroviruses. Fourth, an ancient horizontal gene transfer event was discovered from anamniotes to ERVs with GTE. These findings reveal a complex evolution pattern for retroviral Env in anamniotes.

Deep-Time Structural Evolution of Retroviral and Filoviral Surface Envelope Proteins

The retroviral surface envelope protein subunit (SU) mediates receptor binding and triggers membrane fusion by the transmembrane (TM) subunit. SU evolves rapidly under strong selective conditions, resulting in seemingly unrelated SU structures in highly divergent retroviruses. Structural modeling of the SUs of several retroviruses and related endogenous retroviral elements with AlphaFold 2 identifies a TM-proximal SU β-sandwich structure that has been conserved in the orthoretroviruses for at least 110 million years. The SU of orthoretroviruses diversified by the differential expansion of the β-sandwich core to form domains involved in virus-host interactions. The β-sandwich domain is also conserved in the SU equivalent GP₁ of Ebola virus although with a significantly different orientation in the trimeric envelope protein structure relative to the β-sandwich of human immunodeficiency virus type 1 gp120, with significant evidence for divergent rather than convergent evolution. The unified structural view of orthoretroviral SU and filoviral GP₁ identifies an ancient, structurally conserved, and evolvable domain underlying the structural diversity of orthoretroviral SU and filoviral GP₁. IMPORTANCE The structural relationships of SUs of retroviral groups are obscured by the high rate of sequence change of SU and the deep-time divergence of retroviral lineages. Previous data showed no structural or functional relationships between the SUs of type C gammaretroviruses and lentiviruses. A deeper understanding of structural relationships between the SUs of different retroviral lineages would allow the generalization of critical processes mediated by these proteins in host cell infection. Modeling of SUs with AlphaFold 2 reveals a conserved core domain underlying the structural diversity of orthoretroviral SUs. Definition of the conserved SU structural core allowed the identification of a homologue structure in the SU equivalent GP₁ of filoviruses that most likely shares an origin, unifying the SU of orthoretroviruses and GP₁ of filoviruses into a single protein family. These findings will allow an understanding of the structural basis for receptor-mediated membrane fusion mechanisms in a broad range of biomedically important retroviruses.

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.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

A novel class III endogenous retrovirus with a class I envelope gene in African frogs with an intact genome and developmentally regulated transcripts in Xenopus tropicalis

BACKGROUND

Retroviruses exist as exogenous infectious agents and as endogenous retroviruses (ERVs) integrated into host chromosomes. Such endogenous retroviruses (ERVs) are grouped into three classes roughly corresponding to the seven genera of infectious retroviruses: class I (gamma-, epsilonretroviruses), class II (alpha-, beta-, delta-, lentiretroviruses) and class III (spumaretroviruses). Some ERVs have counterparts among the known infectious retroviruses, while others represent paleovirological relics of extinct or undiscovered retroviruses.

RESULTS

Here we identify an intact ERV in the Anuran amphibian, Xenopus tropicalis. XtERV-S has open reading frames (ORFs) for gag, pol (polymerase) and env (envelope) genes, with a small additional ORF in pol and a serine tRNA primer binding site. It has unusual features and domain relationships to known retroviruses. Analyses based on phylogeny and functional motifs establish that XtERV-S gag and pol genes are related to the ancient env-less class III ERV-L family but the surface subunit of env is unrelated to known retroviruses while its transmembrane subunit is class I-like. LTR constructs show transcriptional activity, and XtERV-S transcripts are detected in embryos after the maternal to zygotic mid-blastula transition and before the late tailbud stage. Tagged Gag protein shows typical subcellular localization. The presence of ORFs in all three protein-coding regions along with identical 5' and 3' LTRs (long terminal repeats) indicate this is a very recent germline acquisition. There are older, full-length, nonorthologous, defective copies in Xenopus laevis and the distantly related African bullfrog, Pyxicephalus adspersus. Additional older, internally deleted copies in X. tropicalis carry a 300 bp LTR substitution.

CONCLUSIONS

XtERV-S represents a genera-spanning member of the largely env-less class III ERV that has ancient and modern copies in Anurans. This provirus has an env ORF with a surface subunit unrelated to known retroviruses and a transmembrane subunit related to class I gammaretroviruses in sequence and organization, and is expressed in early embryogenesis. Additional XtERV-S-related but defective copies are present in X. tropicalis and other African frog taxa. XtERV-S is an unusual class III ERV variant, and it may represent an important transitional retroviral form that has been spreading in African frogs for tens of millions of years.

The evolution of endogenous retroviral envelope genes in bats and their potential contribution to host biology

Bats are the primary reservoirs and carriers of a wide range of viruses of unknown infectivity and pathogenic potential. Some of those if transmitted to other species can cause enormous economic losses in agriculture, and mortality in humans. Bats can be persistently infected with viruses while not showing any symptoms of disease, despite having high virus titre levels in their tissues and shedding virions for months or years after primary infection. It has been suggested that the lack of symptoms of viral infections and low mortality rate in bats might be due to immune adaptations that result from their long-term co-evolution with viruses. In this study, we screened all publicly available bat genomes from six bat families within which we have identified several envelope sequences of retroviral origin (gammaretroviruses). We analysed the identified sequences with Bayesian methods and maximum-likelihood inference to generate a phylogenetic tree with additional reference sequences of known endogenous and exogenous viral envelope genes. We also identified groups of orthologous viral envelopes and analysed them to determine if any of them might be an EVE (endogenous virus element) with an EDI (EVE- derived immunity) function or a candidate for a bat syncytin gene, which is an endogenized viral envelope, mostly known from its function in placentation in animals. Our study shows that bat genomes contain a substantial number of large, intact envelopes with open reading frames, which were found clustering closely on a phylogenetic tree reconstruction with syncytin sequences of other species. That might indicate that such sequences are good candidates for further bat-syncytin/EDI search.

Convergent Co-option of the Retroviral gag Gene during the Early Evolution of Mammals

Endogenous retroviruses, records of past retroviral infections, are ubiquitous in vertebrate genomes. On occasion, vertebrate hosts have co-opted retroviral genes for their own biological functions. Here, we perform a phylogenomic survey of retroviral gag gene homologs within vertebrate genomes and identify two ancient co-opted retroviral gag genes, designated wucaishi1 (wcs1) and wucaishi2 (wcs2), in mammals. Conserved synteny and evolutionary analyses suggest that the wcs1 and wcs2 co-options occurred before the origin of modern placental mammals (∼100 million years ago) and before the origin of modern marsupials (∼80 million years ago), respectively. We found that the wcs genes were lost or pseudogenized multiple times during the evolutionary course of mammals. While the wcs1 gene is mainly subject to negative selection in placental mammals (except in Perissodactyla), the wcs2 gene underwent positive selection in marsupials. Moreover, analyses of transcriptome-sequencing (RNA-seq) data suggest that the wcs1 and the wcs2 genes are expressed in a wide range of tissues. The convergent wcs co-option in mammals implies the retroviral gag gene might have been repurposed more frequently than previously thought.IMPORTANCE Retroviruses occasionally can infect host germ lines, forming endogenous retroviruses. Vertebrates, in turn, recruited retroviral genes for their own biological functions, a process formally known as co-option or exaptation. To date, co-opted retroviral gag genes have rarely been reported. In this study, we identified two co-opted retroviral gag genes, designated wucaishi1 (wcs1) and wucaishi2 (wcs2), in mammals. The co-option of wcs1 and wcs2 occurred before the origin of modern placentals and before the origin of modern marsupials, respectively. Our study indicates that retroviral gag gene co-option might have occurred more frequently than previously thought during the evolutionary course of vertebrates.

The placenta goes viral: Retroviruses control gene expression in pregnancy

The co-option of endogenous retroviruses (ERVs) is increasingly recognized as a recurrent theme in placental biology, which has far-reaching implications for our understanding of mammalian evolution and reproductive health. Most research in this area has focused on ERV-derived proteins, which have been repeatedly co-opted to promote cell–cell fusion and immune modulation in the placenta. ERVs also harbor regulatory sequences that can potentially control placental gene expression, but there has been limited evidence to support this role. In a recent study, Dunn-Fletcher and colleagues discover a striking example of an ERV-derived enhancer element that has been co-opted to regulate a gene important for human pregnancy. Using genomic and experimental approaches, they firmly establish that a primate-specific ERV functions as a placenta-specific enhancer for corticotropin-releasing hormone (CRH), a hormone linked to the control of birth timing in humans. Their findings implicate an extensive yet understudied role for retroviruses in shaping the evolution of placental gene regulatory networks.

An endogenous retroviral envelope syncytin and its cognate receptor identified in the viviparous placental Mabuya lizard

Syncytins are envelope genes from endogenous retroviruses that have been captured during evolution for a function in placentation. They have been found in all placental mammals in which they have been searched, including marsupials. Placental structures are not restricted to mammals but also emerged in some other vertebrates, most frequently in lizards, such as the viviparous Mabuya Scincidae. Here, we performed high-throughput RNA sequencing of a Mabuya placenta transcriptome and screened for the presence of retroviral env genes with a full-length ORF. We identified one such gene, which we named "syncytin-Mab1," that has all the characteristics expected for a syncytin gene. It encodes a membrane-bound envelope protein with fusogenic activity ex vivo, is expressed at the placental level as revealed by in situ hybridization and immunohistochemistry, and is conserved in all Mabuya species tested, spanning over 25 My of evolution. Its cognate receptor, required for its fusogenic activity, was searched for by a screening assay using the GeneBridge4 human/Chinese hamster radiation hybrid panel and found to be the MPZL1 gene, previously identified in mammals as a signal-transducing transmembrane protein involved in cell migration. Together, these results show that syncytin capture is not restricted to placental mammals, but can also take place in the rare nonmammalian vertebrates in which a viviparous placentotrophic mode of reproduction emerged. It suggests that similar molecular tools have been used for the convergent evolution of placentation in independently evolved and highly distant vertebrates.

Co-option of endogenous viral sequences for host cell function

Eukaryotic genomes are littered with sequences of diverse viral origins, termed endogenous viral elements (EVEs). Here we used examples primarily drawn from mammalian endogenous retroviruses to document how the influx of EVEs has provided a source of prefabricated coding and regulatory sequences that were formerly utilized for viral infection and replication, but have been occasionally repurposed for cellular function. While EVE co-option has benefited a variety of host biological functions, there appears to be a disproportionate contribution to immunity and antiviral defense. The mammalian embryo and placenta offer opportunistic routes of viral transmission to the next host generation and as such they represent hotbeds for EVE cooption. Based on these observations, we propose that EVE cooption is initially driven as a mean to mitigate conflicts between host and viruses, which in turn acts as a stepping-stone toward the evolution of cellular innovations serving host physiology and development.