Comprehensive genomic analysis reveals dynamic evolution of endogenous retroviruses that code for retroviral-like protein domains

Aberrant expression of human endogenous retrovirus K9-derived elements is associated with better clinical outcome of acute myelocytic leukemia

BACKGROUND

Acute myelocytic leukemia (AML) is a common hematological malignancy in adults. Although several risk stratifications based on cytogenetic and molecular abnormalities are available to guide the indications for allogeneic hematopoietic cell transplantation (allo-HCT), determining optimal treatment strategies for AML remains challenging. In this study, using transcriptome datasets, we investigated the association between event-free survival (EFS) in intensively treated AML patients and the aberrant expression of endogenous viral element (EVE)-derived open reading frames (ORFs), which have been reported to be associated with the pathophysiology of various malignancies and have the potential to serve as neoantigens in specific cancers.

RESULTS

The expression levels of human endogenous retrovirus family K9 (HERVK9) ORFs were associated with EFS, independent of conventional risk stratification. Furthermore, AML cells with higher levels of HERVK9 expression exhibited enhanced antigen processing and presentation, along with increased expression of genes associated with adaptive immune responses and apoptosis, indicating that aberrant HERVK9 expression may initiate an anti-neoplastic immune response via increased antigen presentation.

CONCLUSIONS

HERVK9 expression may have serve as a crucial prognostic indicator that could aid in determining the indications for upfront allo-HCT in AML patients.

Association of a 7.9 kb Endogenous Retrovirus Insertion in Intron 1 of CD36 with Obesity and Fat Measurements in Sheep

BACKGROUND

Endogenous retroviruses (ERVs) enhance genetic diversity in vertebrates, including sheep. This study investigates the role of Ov-ERV-R13-CD36 within CD36 gene and its association with phenotypic traits in sheep. Analyzing 58 sheep genomes revealed that ERVs constitute approximately 6.02% to 10.05% of the genomic content. We identified 31 retroviral insertion polymorphisms (RIPs) from 28 ERV groups. Among these, Ov-ERV-R13-CD36, which is specifically classified as a beta retrovirus, was selected for further analysis due to its location in CD36 gene, known for its role in fat metabolism, obesity (OB), body weight (BW), and body condition score (BCS). We assessed the association of Ov-ERV-R13-CD36 with OB and BCS across six sheep breeds, utilizing data from 1,355 individuals.

RESULTS

Genomic analyses confirmed that Ov-ERV-R13-CD36 is located within CD36 gene on Chromosome 4, with polymorphisms across various sheep genomes. In a subset of 43 genomes, 22 contained the Ov-ERV-R13-CD36 insertion, while 21 exhibited wild-type variants. The studied animals showed variability in BCS and fat content associated with the Ov-ERV-R13-CD36 variant. Notably, Rahmani sheep exhibited a significantly higher BCS (4.62), categorized as obese, while Barki sheep displayed the lowest BCS (2.73), classified as thin to average. The association analysis indicated that sheep with the RIP-/- genotype correlated with higher OB and BCS, particularly in Rahmani and Romanov x Rahmani breeds.

CONCLUSIONS

Findings suggest that Ov-ERV-R13-CD36 within CD36 gene correlates with beneficial economic traits associated with OB and BCS, particularly in Rahmani and Romanov x Rahmani breeds. This indicates that Ov-ERV-R13-CD36 could be a valuable genetic marker for breeding programs aimed at enhancing traits like fat deposition and body condition in sheep.

Endogenous retroviral ERVH48-1 promotes human urine cell reprogramming

Endogenous retroviruses (ERVs), once thought to be mere remnants of ancient viral integrations in the mammalian genome, are now recognized for their critical roles in various physiological processes, including embryonic development, innate immunity, and tumorigenesis. Their impact on host organisms is significant driver of evolutionary changes, offering insight into evolutionary mechanisms. In our study, we explored the functionality of ERVs by examining single-cell transcriptomic profiles from human embryonic stem cells and urine cells. This led to the discovery of a unique ERVH48-1 expression pattern between these cell types. Additionally, somatic cell reprogramming efficacy was enhanced when ERVH48-1 was overexpressed in a urine cell-reprogramming system. Induced pluripotent stem cells (iPSCs) generated with ERVH48-1 overexpression recapitulated the traits of those produced by traditional reprogramming approaches, and the resulting iPSCs demonstrated the capability to differentiate into all three germ layers in vitro. Our research elucidated the role of ERVs in somatic cell reprogramming.

Evolution and subfamilies of HERVL human endogenous retrovirus

Background

Endogenous retroviruses (ERVs), which blur the boundary between virus and transposable element, are genetic material derived from retroviruses and have important implications for evolution. This study examines the diversity and evolution of human endogenous retroviruses (HERVs) of the HERVL family, which has long terminal repeats (LTRs) named MLT2.

Results

By probability-based sequence comparison, we uncover systematic annotation errors that conceal the true complexity and diversity of transposable elements (TEs) in the human genome. Our analysis identifies new subfamilies within the MLT2 group, proposes a refined classification scheme, and constructs new consensus sequences. We present an evolutionary analysis including phylogenetic trees that elucidate the relationships between these subfamilies and their contributions to human evolution. The results underscore the significance of accurate TE annotation in understanding genome evolution, highlighting the potential for misclassified TEs to impact interpretations of genomic studies.

Availability and implementation

Not applicable.

Retroviral Insertion Polymorphism (RIP) of Porcine Endogenous Retroviruses (PERVs) in Pig Genomes

Endogenous retroviruses (ERVs) are one of the superfamilies of long terminal repeat retrotransposons (LTRs) in mice and humans. Approximately 8% of the pig genome is composed of sequences derived from LTRs. While the majority of ERVs in pigs have decayed, a small number of full-length copies can still mobilize within the genome. This study investigated the unexplored retroviral insertion polymorphisms (RIPs) generated by the mobilization of full-length ERVs (Fl-ERVs), and evaluated their impact on phenotypic variation to gain insights into the biological role of Fl-ERVs in pigs. Overall, 39 RIPs (insertions or deletions relative to the pig reference genome) generated by Fl-ERVs were predicted by comparative genomic analysis, and 18 of them were confirmed by PCR detection. Four RIP sites (D5, D14, D15, and D18) were further evaluated by population analysis, and all of them displayed polymorphisms in multiple breeds. The RIP site of ERV-D14, which is a Fl-ERV inserted in the STAB2-like gene, was further confirmed by sequencing. Population analysis of the polymorphic site of ERV-D14 reveals that it presents moderate polymorphism information in the Large White pig breed, and the association analysis reveals that the RIP of ERV-D14 is associated with age variations at 30 kg body weight (p < 0.05) and 100 kg body weight (p < 0.01) in the population of Large White pigs (N = 480). Furthermore, the ERV-D14 RIP is associated with changes in the expression of the target gene STAB2-like in the liver, backfat, and leaf fat in Sushan pigs. These data suggest that some Fl-ERVs are still mobilizing in the pig's genome, and contribute to genomic and phenotypic variations.

Retrotransposon-derived transcripts and their functions in immunity and disease

Retrotransposons, which account for approximately 42% of the human genome, have been increasingly recognized as "non-self" pathogen-associated molecular patterns (PAMPs) due to their virus-like sequences. In abnormal conditions such as cancer and viral infections, retrotransposons that are aberrantly expressed due to impaired epigenetic suppression display PAMPs, leading to their recognition by pattern recognition receptors (PRRs) of the innate immune system and triggering inflammation. This viral mimicry mechanism has been observed in various human diseases, including aging and autoimmune disorders. However, recent evidence suggests that retrotransposons possess highly regulated immune reactivity and play important roles in the development and function of the immune system. In this review, I discuss a wide range of retrotransposon-derived transcripts, their role as targets in immune recognition, and the diseases associated with retrotransposon activity. Furthermore, I explore the implications of chimeric transcripts formed between retrotransposons and known gene mRNAs, which have been previously underestimated, for the increase of immune-related gene isoforms and their influence on immune function. Retrotransposon-derived transcripts have profound and multifaceted effects on immune system function. The aim of this comprehensive review is to provide a better understanding of the complex relationship between retrotransposon transcripts and immune defense.

Ribosomal profiling of human endogenous retroviruses in healthy tissues

Human endogenous retroviruses (HERVs) are the germline embedded proviral fragments of ancient retroviral infections that make up roughly 8% of the human genome. Our understanding of HERVs in physiology primarily surrounds their non-coding functions, while their protein coding capacity remains virtually uncharacterized. Therefore, we applied the bioinformatic pipeline "hervQuant" to high-resolution ribosomal profiling of healthy tissues to provide a comprehensive overview of translationally active HERVs. We find that HERVs account for 0.1-0.4% of all translation in distinct tissue-specific profiles. Collectively, our study further supports claims that HERVs are actively translated throughout healthy tissues to provide sequences of retroviral origin to the human proteome.

Using long-read CAGE sequencing to profile cryptic-promoter-derived transcripts and their contribution to the immunopeptidome

Recent studies have shown that the noncoding genome can produce unannotated proteins as antigens that induce immune response. One major source of this activity is the aberrant epigenetic reactivation of transposable elements (TEs). In tumors, TEs often provide cryptic or alternate promoters, which can generate transcripts that encode tumor-specific unannotated proteins. Thus, TE-derived transcripts (TE transcripts) have the potential to produce tumor-specific, but recurrent, antigens shared among many tumors. Identification of TE-derived tumor antigens holds the promise to improve cancer immunotherapy approaches; however, current genomics and computational tools are not optimized for their detection. Here we combined CAGE technology with full-length long-read transcriptome sequencing (long-read CAGE, or LRCAGE) and developed a suite of computational tools to significantly improve immunopeptidome detection by incorporating TE and other tumor transcripts into the proteome database. By applying our methods to human lung cancer cell line H1299 data, we show that long-read technology significantly improves mapping of promoters with low mappability scores and that LRCAGE guarantees accurate construction of uncharacterized 5' transcript structure. Augmenting a reference proteome database with newly characterized transcripts enabled us to detect noncanonical antigens from HLA-pulldown LC-MS/MS data. Lastly, we show that epigenetic treatment increased the number of noncanonical antigens, particularly those encoded by TE transcripts, which might expand the pool of targetable antigens for cancers with low mutational burden.

Human LINE-1 retrotransposons: impacts on the genome and regulation by host factors

Genome sequencing revealed that nearly half of the human genome is comprised of transposable elements. Although most of these elements have been rendered inactive due to mutations, full-length intact long interspersed element-1 (LINE-1 or L1) copies retain the ability to mobilize through RNA intermediates by a so-called "copy-and-paste" mechanism, termed retrotransposition. L1 is the only known autonomous mobile genetic element in the genome, and its retrotransposition contributes to inter- or intra-individual genetic variation within the human population. However, L1 retrotransposition also poses a threat to genome integrity due to gene disruption and chromosomal instability. Moreover, recent studies suggest that aberrant L1 expression can impact human health by causing diseases such as cancer and chronic inflammation that might lead to autoimmune disorders. To counteract these adverse effects, the host cells have evolved multiple layers of defense mechanisms at the epigenetic, RNA and protein levels. Intriguingly, several host factors have also been reported to facilitate L1 retrotransposition, suggesting that there is competition between negative and positive regulation of L1 by host factors. Here, we summarize the known host proteins that regulate L1 activity at different stages of the replication cycle and discuss how these factors modulate disease-associated phenotypes caused by L1.

TDP-43 Epigenetic Facets and Their Neurodegenerative Implications

Since its initial involvement in numerous neurodegenerative pathologies in 2006, either as a principal actor or as a cofactor, new pathologies implicating transactive response (TAR) DNA-binding protein 43 (TDP-43) are regularly emerging also beyond the neuronal system. This reflects the fact that TDP-43 functions are particularly complex and broad in a great variety of human cells. In neurodegenerative diseases, this protein is often pathologically delocalized to the cytoplasm, where it irreversibly aggregates and is subjected to various post-translational modifications such as phosphorylation, polyubiquitination, and cleavage. Until a few years ago, the research emphasis has been focused particularly on the impacts of this aggregation and/or on its widely described role in complex RNA splicing, whether related to loss- or gain-of-function mechanisms. Interestingly, recent studies have strengthened the knowledge of TDP-43 activity at the chromatin level and its implication in the regulation of DNA transcription and stability. These discoveries have highlighted new features regarding its own transcriptional regulation and suggested additional mechanistic and disease models for the effects of TPD-43. In this review, we aim to give a comprehensive view of the potential epigenetic (de)regulations driven by (and driving) this multitask DNA/RNA-binding protein.

Transposons contribute to the functional diversification of the head, gut, and ovary transcriptomes across Drosophila natural strains

Transcriptomes are dynamic, with cells, tissues, and body parts expressing particular sets of transcripts. Transposable elements (TEs) are a known source of transcriptome diversity; however, studies often focus on a particular type of chimeric transcript, analyze single body parts or cell types, or are based on incomplete TE annotations from a single reference genome. In this work, we have implemented a method based on de novo transcriptome assembly that minimizes the potential sources of errors while identifying a comprehensive set of gene-TE chimeras. We applied this method to the head, gut, and ovary dissected from five Drosophila melanogaster natural strains, with individual reference genomes available. We found that ∼19% of body part-specific transcripts are gene-TE chimeras. Overall, chimeric transcripts contribute a mean of 43% to the total gene expression, and they provide protein domains for DNA binding, catalytic activity, and DNA polymerase activity. Our comprehensive data set is a rich resource for follow-up analysis. Moreover, because TEs are present in virtually all species sequenced to date, their role in spatially restricted transcript expression is likely not exclusive to the species analyzed in this work.

Unique Structure and Distinctive Properties of the Ancient and Ubiquitous Gamma-Type Envelope Glycoprotein

After the onset of the AIDS pandemic, HIV-1 (genus Lentivirus) became the predominant model for studying retrovirus Env glycoproteins and their role in entry. However, HIV Env is an inadequate model for understanding entry of viruses in the Alpharetrovirus, Gammaretrovirus and Deltaretrovirus genera. For example, oncogenic model system viruses such as Rous sarcoma virus (RSV, Alpharetrovirus), murine leukemia virus (MLV, Gammaretrovirus) and human T-cell leukemia viruses (HTLV-I and HTLV-II, Deltaretrovirus) encode Envs that are structurally and functionally distinct from HIV Env. We refer to these as Gamma-type Envs. Gamma-type Envs are probably the most widespread retroviral Envs in nature. They are found in exogenous and endogenous retroviruses representing a broad spectrum of vertebrate hosts including amphibians, birds, reptiles, mammals and fish. In endogenous form, gamma-type Envs have been evolutionarily coopted numerous times, most notably as placental syncytins (e.g., human SYNC1 and SYNC2). Remarkably, gamma-type Envs are also found outside of the Retroviridae. Gp2 proteins of filoviruses (e.g., Ebolavirus) and snake arenaviruses in the genus Reptarenavirus are gamma-type Env homologs, products of ancient recombination events involving viruses of different Baltimore classes. Distinctive hallmarks of gamma-type Envs include a labile disulfide bond linking the surface and transmembrane subunits, a multi-stage attachment and fusion mechanism, a highly conserved (but poorly understood) "immunosuppressive domain", and activation by the viral protease during virion maturation. Here, we synthesize work from diverse retrovirus model systems to illustrate these distinctive properties and to highlight avenues for further exploration of gamma-type Env structure and function.

Monotreme-Specific Conserved Putative Proteins Derived from Retroviral Reverse Transcriptase

Endogenous retroviruses (ERVs) have played an essential role in the evolution of mammals. ERV-derived genes are reported in the therians, many of which are involved in placental development; however, the contribution of the ERV-derived genes in monotremes, which are oviparous mammals, remains to be uncovered. Here, we conducted a comprehensive search for possible ERV-derived genes in platypus and echidna genomes and identified three reverse transcriptase-like genes named RTOM1, RTOM2, and RTOM3 clustered in the GRIP2 intron. Comparative genomic analyses revealed that RTOM1, RTOM2, and RTOM3 are strongly conserved and are under purifying selection between these species. These could be generated by tandem duplications before the divergence of platypus and echidna. All RTOM transcripts were specifically expressed in the testis, possibly suggesting their physiological importance. This is the first study reporting monotreme-specific de novo gene candidates derived from ERVs, which provides new insights into the unique evolution of monotremes.

Mammalian genome innovation through transposon domestication

Since the discovery of transposons, their sheer abundance in host genomes has puzzled many. While historically viewed as largely harmless 'parasitic' DNAs during evolution, transposons are not a mere record of ancient genome invasion. Instead, nearly every element of transposon biology has been integrated into host biology. Here we review how host genome sequences introduced by transposon activities provide raw material for genome innovation and document the distinct evolutionary path of each species.

Viruses Defined by the Position of the Virosphere within the Replicator Space

Originally, viruses were defined as miniscule infectious agents that passed through filters that retain even the smallest cells. Subsequently, viruses were considered obligate intracellular parasites whose reproduction depends on their cellular hosts for energy supply and molecular building blocks. However, these features are insufficient to unambiguously define viruses as they are broadly understood today. We outline possible approaches to define viruses and explore the boundaries of the virosphere within the virtual space of replicators and the relationships between viruses and other types of replicators. Regardless of how, exactly, viruses are defined, viruses clearly have evolved on many occasions from nonviral replicators, such as plasmids, by recruiting host proteins to become virion components. Conversely, other types of replicators have repeatedly evolved from viruses. Thus, the virosphere is a dynamic entity with extensive evolutionary traffic across its boundaries. We argue that the virosphere proper, here termed orthovirosphere, consists of a distinct variety of replicators that encode structural proteins encasing the replicators' genomes, thereby providing protection and facilitating transmission among hosts. Numerous and diverse replicators, such as virus-derived but capsidless RNA and DNA elements, or defective viruses occupy the zone surrounding the orthovirosphere in the virtual replicator space. We define this zone as the perivirosphere. Although intense debates on the nature of certain replicators that adorn the internal and external boundaries of the virosphere will likely continue, we present an operational definition of virus that recently has been accepted by the International Committee on Taxonomy of Viruses.

The Oldest Co-opted gag Gene of a Human Endogenous Retrovirus Shows Placenta-Specific Expression and Is Upregulated in Diffuse Large B-Cell Lymphomas

Vertebrate genomes contain endogenous retroviruses (ERVs) that represent remnants of past germline infections by ancient retroviruses. Despite comprising 8% of the human genome, the human ERVs (HERVs) do not encode a replication competent retrovirus. However, some HERV genes have been co-opted to serve host functions, most notably the viral envelope-derived syncytins involved in placentation. Here, we identify the oldest HERV intact gag gene with an open reading frame, gagV1. Its provirus contains an intact env, envV1, and the first open reading frame found in an HERV gag leader, pre-gagV1, which encodes a novel protein. This HERV is linked to a related gag gene, gagV3, and these three genes all show patterns of evolutionary conservation in primates. gagV1 and pre-gagV1 orthologs are present in all simian primate lineages indicating that this HERV entered the germline of the common simian primate ancestor at least 43 Ma, whereas gagV3 is found in Old and New World monkeys. gagV1 and gagV3 have undergone recurrent gene conversion events and positive selection. Expression of gagV1, gagV3, and pre-gagV1 is restricted to the placenta in humans and macaques suggesting co-option for placenta-specific host functions. Transcriptomic analysis of human tumors also found upregulated levels of gagV1 transcripts in diffuse large B-cell lymphomas. These findings suggest that these HERV-V genes may be useful markers for the most common type of non-Hodgkin's lymphoma and that they may have contributed to the successive domestications of env and gag genes in eutherians involved in the ongoing ERV-driven evolution of the placenta.

Structure of a Ty1 restriction factor reveals the molecular basis of transposition copy number control

Excessive replication of Saccharomyces cerevisiae Ty1 retrotransposons is regulated by Copy Number Control, a process requiring the p22/p18 protein produced from a sub-genomic transcript initiated within Ty1 GAG. In retrotransposition, Gag performs the capsid functions required for replication and re-integration. To minimize genomic damage, p22/p18 interrupts virus-like particle function by interaction with Gag. Here, we present structural, biophysical and genetic analyses of p18m, a minimal fragment of Gag that restricts transposition. The 2.8 Å crystal structure of p18m reveals an all α-helical protein related to mammalian and insect ARC proteins. p18m retains the capacity to dimerise in solution and the crystal structures reveal two exclusive dimer interfaces. We probe our findings through biophysical analysis of interface mutants as well as Ty1 transposition and p18m restriction in vivo. Our data provide insight into Ty1 Gag structure and suggest how p22/p18 might function in restriction through a blocking-of-assembly mechanism.

In Saccharomyces cerevisiae, unchecked proliferation of Ty1 retrotransposons is controlled by the process of copy number control (CNC), which requires the p22/p18 protein, translated from an internal transcript within the Ty1 GAG gene. Here, the authors present the 2.8 Å crystal structure of a minimal p18 from Ty1-Gag that is able to restrict Ty1 transposition and identify two dimer interfaces in p18, whose roles were probed by mutagenesis both in vitro and in vivo. As p22/p18 contains only one of two conserved domains required for retroelement Gag assembly, they propose that p22/p18-Gag interactions block the Ty1 virus-like particle assembly pathway, resulting in defective particles incapable of supporting retrotransposition.

HERV-Derived Ervpb1 Is Conserved in Simiiformes, Exhibiting Expression in Hematopoietic Cell Lineages Including Macrophages

(1) Background: The ERVPb1 gene in humans is derived from an envelope (Env) gene of a human endogenous retrovirus group, HERV-P(b). The ERVPb1 gene reportedly has a conserved open reading frame (ORF) in Old World monkeys. Although its forced expression led to cell-fusion in an ex vivo cell culture system, like other Env-derived genes such as syncytin-1 and -2, its mRNA expression is not placenta-specific, but almost ubiquitous, albeit being quite low in human tissues and organs, implying a distinct role for ERVPb1. (2) Methods: To elucidate the cell lineage(s) in which the ERVPb1 protein is translated in human development, we developed a novel, highly sensitive system for detecting HERV-derived proteins/peptides expressed in the tissue differentiation process of human induced pluripotent stem cells (iPSCs). (3) Results: We first determined that ERVPb1 is also conserved in New World monkeys. Then, we showed that the ERVPb1 protein is translated from a uniquely spliced ERVPb1 transcript in hematopoietic cell lineages, including a subset of macrophages, and further showed that its mRNA expression is upregulated by lipopolysaccharide (LPS) stimulation in primary human monocytes. (4) Conclusions: ERVPb1 is unique to Simiiformes and actually translated in hematopoietic cell lineages, including a subset of macrophages.

Expression of ERV3-1 in leukocytes of acute myelogenous leukemia patients

Acute myelogenous leukemia (AML) is one of the major hematological malignancies. In the human genome, several have been found to originate from retroviruses, and some of which are involved in the progression of various cancers. Hence, to investigate whether retroviral-like genes are associated with AML development, we conducted a transcriptome sequencing analysis of 12 retroviral-like genes of 150 AML patients and 32 healthy donor samples, of which RNA sequencing data were obtained from public databases. We found high expression of ERV3-1, an envelope gene of endogenous retrovirus group 3 member 1, in all AML patients examined in this study. In particular, blood and bone marrow cells of the myeloid lineage in AML patients, exhibited higher expression of ERV3-1 than those of the monocytic AML lineage. We also examined the protein expression of ERV3-1 by immunohistochemical analysis and found expression of the ERV3-1 protein in all 12 myeloid-phenotype patients and 7 out of 12 monocytic-phenotype patients, with a particular concentration observed at the membrane of some leukemic cells. Transcriptome analysis further suggested that upregulated ERV3-1 expression may be associated with chromosome 8 trisomy as anomaly was found to be more common among the high expression group than the low expression group. However, this finding was not corroborated by the immunohistochemical data. This discrepancy may have been caused, in part, by the small number of samples analyzed in this study. Although the precise associated molecular mechanisms remain unclear, our results suggest that ERV3-1 may be involved in AML development.