Comprehensive Characterization of the Human Endogenous Retrovirus HERV-K(HML-6) Group: Overview of Structure, Phylogeny, and Contribution to the Human Genome

Comprehensive characterization of ERV-K (HML-8) in the chimpanzee genome revealed less genomic activity than humans

Endogenous retroviruses (ERVs) originate from ancestral germline infections caused by exogenous retroviruses. Throughout evolution, they have become fixed within the genome of the animals into which they were integrated. As ERV elements coevolve with the host, they are normally epigenetically silenced and can become upregulated in a series of physiological and pathological processes. Generally, a detailed ERV profile in the host genome is critical for understanding the evolutionary history and functional performance of the host genome. We previously characterized and cataloged all the ERV-K subtype HML-8 loci in the human genome; however, this has not been done for the chimpanzee, the nearest living relative of humans. In this study, we aimed to catalog and characterize the integration of HML-8 in the chimpanzee genome and compare it with the integration of HML-8 in the human genome. We analyzed the integration of HML-8 and found that HML-8 pervasively invaded the chimpanzee genome. A total of 76 proviral elements were characterized on 23/24 chromosomes, including detailed elements distribution, structure, phylogeny, integration time, and their potential to regulate adjacent genes. The incomplete structure of HML-8 proviral LTRs will undoubtedly affect their activity. Moreover, the results indicated that HML-8 integration occurred before the divergence between humans and chimpanzees. Furthermore, chimpanzees include more HML-8 proviral elements (76 vs. 40) and fewer solo long terminal repeats (LTR) (0 vs. 5) than humans. These results suggested that chimpanzee genome activity is less than the human genome and that humans may have a better ability to shape and screen integrated proviral elements. Our work is informative in both an evolutionary and a functional context for ERVs.

Wildlife endogenous retroviruses: colonization, consequences, and cooption

Endogenous retroviruses (ERVs) are inherited genomic remains of past germline retroviral infections. Research on human ERVs has focused on medical implications of their dysregulation on various diseases. However, recent studies incorporating wildlife are yielding remarkable perspectives on long-term retrovirus-host interactions. These initial forays into broader taxonomic analysis, including sequencing of multiple individuals per species, show the incredible plasticity and variation of ERVs within and among wildlife species. This demonstrates that stochastic processes govern much of the vertebrate genome. In this review, we elaborate on discoveries pertaining to wildlife ERV origins and evolution, genome colonization, and consequences for host biology.

In Vitro Expression Analysis Reveals HML6-c14 to Be an Attractive Research Target

HML6-c14, a long terminal repeat (LTR)-type retrotransposon identified by expressed sequence tag (EST) database screening, was found to undergo RNA processing resembling that of placental tissue by in vitro expression analysis. Previous in situ hybridization studies using normal placental tissue showed that the transcript remained in the nucleus. However, among the transcripts forcedly expressed in cultured cells, the transcript that retained the 3.3 kb intron was observed in the nucleus, and a part of the spliced transcript was observed outside the nucleus. To verify whether this cytoplasmic transcript could be translated, we examined the coding potential of the open reading frame (ORF), consisting of 109 codons on the spliced transcript, along with two other putative ORFs detected in the intronic region. As a result, none of the ORF-derived products could be detected by Western blotting as fusion proteins tagged with the FLAG epitope, suggesting that HML6-c14 belongs to a group of long non-coding RNA (lncRNA) genes. Promoter analysis of the upstream 6.4 kb genomic region also suggested that the 5'-LTR itself potentially retains high promoter activity. Despite losing the ability to produce functional proteins, HML6-c14 continues to retain its transcriptional ability while converting to an lncRNA gene, which is an interesting subject for future research.

Anti-HERV-K Drugs and Vaccines, Possible Therapies against Tumors

The footprint of human endogenous retroviruses (HERV), specifically HERV-K, has been found in malignancies, such as melanoma, teratocarcinoma, osteosarcoma, breast cancer, lymphoma, and ovary and prostate cancers. HERV-K is characterized as the most biologically active HERV due to possession of open reading frames (ORF) for all Gag, Pol, and Env genes, which enables it to be more infective and obstructive towards specific cell lines and other exogenous viruses, respectively. Some factors might contribute to carcinogenicity and at least one of them has been recognized in various tumors, including overexpression/methylation of long interspersed nuclear element 1 (LINE-1), HERV-K Gag, and Env genes themselves plus their transcripts and protein products, and HERV-K reverse transcriptase (RT). Therapies effective for HERV-K-associated tumors mostly target invasive autoimmune responses or growth of tumors through suppression of HERV-K Gag or Env protein and RT. To design new therapeutic options, more studies are needed to better understand whether HERV-K and its products (Gag/Env transcripts and HERV-K proteins/RT) are the initiators of tumor formation or just the disorder’s developers. Accordingly, this review aims to present evidence that highlights the association between HERV-K and tumorigenicity and introduces some of the available or potential therapies against HERV-K-induced tumors.

Identification and Characterization of the HERV-K (HML-8) Group of Human Endogenous Retroviruses in the Genome

Human endogenous retroviruses (HERVs) can be vertically transmitted in a Mendelian fashion, are stably maintained in the human genome, and are estimated to constitute ∼8% of the genome. HERVs affect human physiology and pathology through their provirus-encoded protein or long terminal repeat (LTR) element effect. Characterization of the genomic distribution is an essential step to understanding the relationships between endogenous retrovirus expression and diseases. However, the poor characterization of human MMTV-like (HML)-8 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, the definition of an accurate and updated HERV-K HML-8 genomic map is urgently needed. In this study, we report the results of a comprehensive analysis of HERV-K HML-8 sequence presence and distribution within the human genome and hominoids, with a detailed description of the different structural and phylogenetic aspects characterizing the group. A total of 40 proviruses and 5 solo LTR elements for human were characterized, which included a detailed description of provirus structure, integration time, potentially regulated genes, transcription factor-binding sites, and primer-binding site features. Besides, 9 chimpanzee sequences, 8 gorilla sequences, and 10 orangutan sequences belonging to the HML-8 subgroup were identified. The integration time results showed that the HML-8 elements were integrated into the primate lineage around 35 and 42 million years ago (mya), during primates evolutionary speciation. Overall, the results clarified the composition of the HML-8 groups, providing an exhaustive background for subsequent functional studies.

Human Endogenous Retrovirus (HERV) Transcriptome Is Dynamically Modulated during SARS-CoV-2 Infection and Allows Discrimination of COVID-19 Clinical Stages

SARS-CoV-2 infection is known to trigger an important inflammatory response, which has a major role in COVID-19 pathogenesis. In infectious and inflammatory contexts, the modulation of human endogenous retroviruses (HERV) has been broadly reported, being able to further sustain innate immune responses due to the expression of immunogenic viral transcripts, including double-stranded DNA (dsRNA), and eventually, immunogenic proteins. To gain insights on this poorly characterized interplay, we performed a high-throughput expression analysis of ~3,300 specific HERV loci in the peripheral blood mononuclear cells (PBMCs) of 10 healthy controls and 16 individuals being either convalescent after the infection (6) or retesting positive after convalescence (10) (Gene Expression Onmibus [GEO] data set GSE166253). Results showed that the exposure to SARS-CoV-2 infection modulates HERV expression according to the disease stage and reflecting COVID-19 immune signatures. The differential expression analysis between healthy control (HC) and COVID-19 patients allowed us to identify a total of 282 differentially expressed HERV loci (deHERV) in the individuals exposed to SARS-CoV-2 infection, independently from the clinical form. In addition, 278 and 60 deHERV loci that were specifically modulated in individuals convalescent after COVID19 infection (C) and patients that retested positive to SARS-CoV-2 after convalescence (RTP) as individually compared to HC, respectively, as well as 164 deHERV loci between C and RTP patients were identified. The identified HERV loci belonged to 36 different HERV groups, including members of all three classes. The present study provides an exhaustive picture of the HERV transcriptome in PBMCs and its dynamic variation in the presence of COVID-19, revealing specific modulation patterns according to the infection stage that can be relevant to the disease clinical manifestation and outcome. IMPORTANCE We report here the first high-throughput analysis of HERV loci expression along SARS-CoV-2 infection, as performed with peripheral blood mononuclear cells (PBMCs). Such cells are not directly infected by the virus but have a crucial role in the plethora of inflammatory and immune events that constitute a major hallmark of COVID-19 pathogenesis. Results provide a novel and exhaustive picture of HERV expression in PBMCs, revealing specific modulation patterns according to the disease condition and the concomitant immune activation. To our knowledge, this is the first set of deHERVs whose expression is dynamically modulated across COVID-19 stages, confirming a tight interplay between HERV and cellular immunity and revealing specific transcriptional signatures that can have an impact on the disease clinical manifestation and outcome.

Identification, comprehensive characterization, and comparative genomics of the HERV-K(HML8) integrations in the human genome

Around 8% of the human genome is composed by Human Endogenous Retroviruses (HERVs), ancient viral sequences inherited from the primate germ line after their infection by now extinct retroviruses. Given the still underexplored physiological and pathological roles of HERVs, it is fundamental to increase our information about the genomic composition of the different groups, to lay reliable foundation for functional studies. Among HERVs, the most characterized elements belong to the beta-like superfamily HERV-K, comprising 10 groups (HML1-10) with HML2 being the most recent and studied one. Among HMLs, the HML8 group is the only one still lacking a comprehensive genomic description. In the present work, we investigated HML8 sequences' distribution in the human genome (GRCh38/hg38), identifying 23 novel proviruses and characterizing the overall 78 HML8 proviruses in terms of genome structure, phylogeny, and integration pattern. HML8 elements were significantly enriched in human chromosomes 8 and X (p<0.005) while chromosomes 17 and 20 showed fewer integrations than expected (p<0.025 and p<0.005, respectively). Phylogenetic analyses classified HML8 members into 3 clusters, corresponding to the three LTR types MER11A, MER11B and MER11C. Besides different LTR types, common signatures in the internal structure suggested the potential existence of three different ancestral HML8 variants. Accordingly, time of integration estimation coupled with comparative genomics revealed that these three clusters have a different time of integration in the primates' genome, with MER11C elements being significantly younger than MER11A- and MER11B associated proviruses (p<0.005 and p<0.05, respectively). Approximately 30% of the HML8 elements were found co-localized within human genes, sometimes in exonic portions and with the same orientation, deserving further studies for their possible effects on gene expression. Overall, we provide the first detailed picture of the HML8 group distribution and variety among the genome, creating the backbone for the specific analysis of their transcriptional activity in healthy and diseased conditions.

Prolonged activity of HERV-K(HML2) in Old World Monkeys accounts for recent integrations and novel recombinant variants

Around 8% of the human genome comprises Human Endogenous Retroviruses (HERVs) acquired over primate evolution. Some are specific to primates such as HERV-K, consisting of 10 HML subtypes and including the most recently acquired elements. Particularly, HML2 is the youngest clade, having some human-specific integrations, and while it has been widely described in humans its presence and distribution in non-human primates remain poorly characterized. To investigate HML2 distribution in non-human primates, the present study focused on the characterization of HML2 integrations in Macaca fascicularis and Macaca mulatta which are the most evolutionarily distant species related to humans in the Catarrhini parvorder. We identified overall 208 HML2 proviruses for M. fascicularis (77) and M. mulatta (131). Among them, 46 proviruses are shared by the two species while the others are species specific. Only 12 proviruses were shared with humans, confirming that the major wave of HML2 diffusion in humans occurred after macaques' divergence. Phylogenetic analysis confirmed structural variations between HML2 macaques' species-specific proviruses, and the ones shared between macaques and humans. The HML2 loci were characterized in terms of structure, focusing on potential residual open reading frames (ORFs) for gag, pol, and env genes for the latter being reported to be expressed in human pathological conditions. The analysis identified highly conserved gag and pol genes, while the env genes had a very divergent nature. Of the 208 HML2 proviral sequences present in Macaca species, 81 sequences form a cluster having a MER11A, a characteristic HML8 LTR sequence, insertion in the env region indicating a recombination event that occurred between the HML2 env gene and the HML8 LTR. This recombination event, which was shown to be present only in a subset of macaques' shared sequences and species-specific sequences, highlights a recent viral activity leading to the emergence of an env variant specific to the Old World Monkeys (OWMs). We performed an exhaustive analysis of HML2 in two species of OWMs, in terms of its evolutionary history, structural features, and potential residual coding capacity highlighting recent activity of HML2 in macaques that occurred after its split from the Catarrhini parvorder, leading to the emergence of viral variants, hence providing a better understanding of the endogenization and diffusion of HML2 along primate evolution.

Comprehensive identification and characterization of the HERV-K (HML-9) group in the human genome

Background

Human endogenous retroviruses (HERVs) result from ancestral infections caused by exogenous retroviruses that became incorporated into the germline DNA and evolutionarily fixed in the human genome. HERVs can be transmitted vertically in a Mendelian fashion and be stably maintained in the human genome, of which they are estimated to comprise approximately 8%. HERV-K (HML1-10) transcription has been confirmed to be associated with a variety of diseases, such as breast cancer, lung cancer, prostate cancer, melanoma, rheumatoid arthritis, and amyotrophic lateral sclerosis. However, the poor characterization of HML-9 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, a precise and updated HERV-K HML-9 genomic map is urgently needed to better evaluate the role of these elements in human health.

Results

We report a comprehensive analysis of the presence and distribution of HERV-K HML-9 elements within the human genome, with a detailed characterization of the structural and phylogenetic properties of the group. A total of 23 proviruses and 47 solo LTR elements were characterized, with a detailed description of the provirus structure, integration time, potential regulated genes, transcription factor binding sites (TFBS), and primer binding site (PBS) features. The integration time results showed that the HML-9 elements found in the human genome integrated into the primate lineage between 17.5 and 48.5 million years ago (mya).

Conclusion

The results provide a clear characterization of HML-9 and a comprehensive background for subsequent functional studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12977-022-00596-2.

Differential expression of an endogenous retroviral element [HERV-K(HML-6)] is associated with reduced survival in glioblastoma patients

Comprising approximately 8% of our genome, Human Endogenous RetroViruses (HERVs) represent a class of germline retroviral infections that are regulated through epigenetic modifications. In cancer cells, which often have epigenetic dysregulation, HERVs have been implicated as potential oncogenic drivers. However, their role in gliomas is not known. Given the link between HERV expression in cancer cell lines and the distinct epigenetic dysregulation in gliomas, we utilized a tailored bioinformatic pipeline to characterize and validate the glioma retrotranscriptome and correlate HERV expression with locus-specific epigenetic modifications. We identified robust overexpression of multiple HERVs in our cell lines, including a retroviral transcript, HML-6, at 19q13.43b in glioblastoma cells. HERV expression inversely correlated with loci-specific DNA methylation. HML-6 contains an intact open reading frame encoding a small envelope protein, ERVK3-1. Increased expression of ERVK3-1 in GBM patients is associated with a poor prognosis independent of IDH-mutational status. Our results suggest that not only is HML-6 uniquely overexpressed in highly invasive cell lines and tissue samples, but also its gene product, ERVK3-1, may be associated with reduced survival in GBM patients. These results may have implications for both the tumor biology of GBM and the role of ERVK3-1 as a potential therapeutic target.

An assessment of bioinformatics tools for the detection of human endogenous retroviral insertions in short-read genome sequencing data

There is a growing interest in the study of human endogenous retroviruses (HERVs) given the substantial body of evidence that implicates them in many human diseases. Although their genomic characterization presents numerous technical challenges, next-generation sequencing (NGS) has shown potential to detect HERV insertions and their polymorphisms in humans. Currently, a number of computational tools to detect them in short-read NGS data exist. In order to design optimal analysis pipelines, an independent evaluation of the available tools is required. We evaluated the performance of a set of such tools using a variety of experimental designs and datasets. These included 50 human short-read whole-genome sequencing samples, matching long and short-read sequencing data, and simulated short-read NGS data. Our results highlight a great performance variability of the tools across the datasets and suggest that different tools might be suitable for different study designs. However, specialized tools designed to detect exclusively human endogenous retroviruses consistently outperformed generalist tools that detect a wider range of transposable elements. We suggest that, if sufficient computing resources are available, using multiple HERV detection tools to obtain a consensus set of insertion loci may be ideal. Furthermore, given that the false positive discovery rate of the tools varied between 8% and 55% across tools and datasets, we recommend the wet lab validation of predicted insertions if DNA samples are available.

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.

Human endogenous retroviruses in development and disease

Human endogenous retroviruses (HERVs) represent ∼8% of human genome, deriving from exogenous retroviral infections of germ line cells occurred millions of years ago and being inherited by the offspring in a Mendelian fashion. Most of HERVs are nonprotein-coding because of the accumulation of mutations, insertions, deletions, and/or truncations. It has been long thought that HERVs were “junk DNA”. However, it is now known that HERVs are involved in various biological processes through encoding proteins, acting as promoters/enhancers, or lncRNAs to affect human health and disease. In this review, we summarized recent findings about HERVs, with implications in embryonic development, pluripotency, cancer, aging, and neurodegenerative diseases.

Infection by Diverse HIV-1 Subtypes Leads to Different Elevations in HERV-K Transcriptional Levels in Human T Cell Lines

Human endogenous retroviruses (HERVs) make up ~8% of the human genome, and for millions of years, they have been subject to strict biological regulation. Many HERVs do not participate in normal physiological activities in the body. However, in some pathological conditions, they can be abnormally activated. For example, HIV infection can cause abnormal activation of HERVs, and under different infection conditions, HERV expression may be different. We observed significant differences in HERV-K transcription levels among HIV-1 subtype-infected individuals. The transcriptional levels in the HERV-K gag region were significantly increased in HIV-1 B subtype-infected patients, while the transcriptional levels in the HERV-K pol region were significantly increased in CRF01_AE and CRF07_BC subtype-infected patients. In vitro, the transcriptional levels of HEVR-K were increased 5-fold and 15-fold in MT2 cells transfected with two different HIV-1 strains (B and CRF01_AE, respectively). However, there was no significant difference in transcriptional levels among regions of HERV-K. When MT2 cells were infected with different subtypes of HIV-1 Tat proteins (B, CRF01_AE), which is constructed by lentiviruses, and the transcription levels of HERV-K were increased 4-fold and 2-fold, respectively. Thus, different subtypes of HIV-1 have different effects on HERV-K transcription levels, which may be caused by many factors, not only Tat protein.

HERV-K(HML7) Integrations in the Human Genome: Comprehensive Characterization and Comparative Analysis in Non-Human Primates

Simple Summary

The human genome is not human at all, but it includes a multitude of sequences inherited from ancient viral infections that affected primates’ germ line. These elements can be seen as the fossils of now-extinct retroviruses, and are called Human Endogenous Retroviruses (HERVs). View as “junk DNA” for a long time, HERVs constitute 4 times the amount of DNA needed to produce all cellular proteins, and growing evidence indicates their crucial role in primate brain evolution, placenta development, and innate immunity shaping. HERVs are also intensively studied for a pathological role, even if the incomplete knowledge about their exact number and genomic position has thus far prevented any causal association. Among possible relevant HERVs, the HERV-K supergroup is of particular interest, including some of the oldest (HML5) as well as youngest (HML2) integrations. Among HERV-Ks, the HML7 group still lack a detailed description, and the present work thus aimed to identify and characterize all HML7 elements in the human genome. Results showed that the HML7 group includes 23 elements and an additional 160 “scars” of past infection that invaded in primates mostly between 20 and 30 million years ago, providing an exhaustive background to study their impact on human pathophysiology.

Abstract

Endogenous Retroviruses (ERVs) are ancient relics of infections that affected the primate germ line and constitute about 8% of our genome. Growing evidence indicates that ERVs had a major role in vertebrate evolution, being occasionally domesticated by the host physiology. In addition, human ERV (HERV) expression is highly investigated for a possible pathological role, even if no clear associations have been reported yet. In fact, on the one side, the study of HERV expression in high-throughput data is a powerful and promising tool to assess their actual dysregulation in diseased conditions; but, on the other side, the poor knowledge about the various HERV group genomic diversity and individual members somehow prevented the association between specific HERV loci and a given molecular mechanism of pathogenesis. The present study is focused on the HERV-K(HML7) group that—differently from the other HERV-K members—still remains poorly characterized. Starting from an initial identification performed with the software RetroTector, we collected 23 HML7 proviral insertions and about 160 HML7 solitary LTRs that were analyzed in terms of genomic distribution, revealing a significant enrichment in chromosome X and the frequent localization within human gene introns as well as in pericentromeric and centromeric regions. Phylogenetic analyses showed that HML7 members form a monophyletic group, which based on age estimation and comparative localization in non-human primates had its major diffusion between 20 and 30 million years ago. Structural characterization revealed that besides 3 complete HML7 proviruses, the other group members shared a highly defective structure that, however, still presents recognizable functional domains, making it worth further investigation in the human population to assess the presence of residual coding potential.

Human Endogenous Retroviruses (HERVs) and Mammalian Apparent LTRs Retrotransposons (MaLRs) Are Dynamically Modulated in Different Stages of Immunity

Human Endogenous retroviruses (HERVs) and Mammalian Apparent LTRs Retrotransposons (MaLRs) are remnants of ancient retroviral infections that represent a large fraction of our genome. The HERV and MaLR transcriptional activity is regulated in developmental stages, adult tissues, and pathological conditions. In this work, we used a bioinformatics approach based on RNA-sequencing (RNA-seq) to study the expression and modulation of HERVs and MaLR in a scenario of activation of the immune response. We analyzed transcriptome data from subjects before and after the administration of an inactivated vaccine against the Hantaan orthohantavirus, the causative agent of Korean hemorrhagic fever, to investigate the HERV and MaLR expression and differential expression in response to the administration of the vaccine. Specifically, we described the HERV transcriptome in PBMCs and identified HERV and MaLR loci differentially expressed after the 2nd, 3rd, and 4th inactivated vaccine administrations. We found that the expression of 545 HERV and MaLR elements increased in response to the vaccine and that the activation of several individual HERV and MaLR loci is specific for each vaccine administration and correlated to different genes and immune-related pathways.

Retrotransposons as Drivers of Mammalian Brain Evolution

Retrotransposons, a large and diverse class of transposable elements that are still active in humans, represent a remarkable force of genomic innovation underlying mammalian evolution. Among the features distinguishing mammals from all other vertebrates, the presence of a neocortex with a peculiar neuronal organization, composition and connectivity is perhaps the one that, by affecting the cognitive abilities of mammals, contributed mostly to their evolutionary success. Among mammals, hominids and especially humans display an extraordinarily expanded cortical volume, an enrichment of the repertoire of neural cell types and more elaborate patterns of neuronal connectivity. Retrotransposon-derived sequences have recently been implicated in multiple layers of gene regulation in the brain, from transcriptional and post-transcriptional control to both local and large-scale three-dimensional chromatin organization. Accordingly, an increasing variety of neurodevelopmental and neurodegenerative conditions are being recognized to be associated with retrotransposon dysregulation. We review here a large body of recent studies lending support to the idea that retrotransposon-dependent evolutionary novelties were crucial for the emergence of mammalian, primate and human peculiarities of brain morphology and function.

Enhanced antitumoral activity of TLR7 agonists via activation of human endogenous retroviruses by HDAC inhibitors

In this work, we are reporting that "Shock and Kill", a therapeutic approach designed to eliminate latent HIV from cell reservoirs, is extrapolatable to cancer therapy. This is based on the observation that malignant cells express a spectrum of human endogenous retroviral elements (HERVs) which can be transcriptionally boosted by HDAC inhibitors. The endoretroviral gene HERV-V2 codes for an envelope protein, which resembles syncytins. It is significantly overexpressed upon exposure to HDAC inhibitors and can be effectively targeted by simultaneous application of TLR7/8 agonists, triggering intrinsic apoptosis. We demonstrated that this synergistic cytotoxic effect was accompanied by the functional disruption of the TLR7/8-NFκB, Akt/PKB, and Ras-MEK-ERK signalling pathways. CRISPR/Cas9 ablation of TLR7 and HERV-V1/V2 curtailed apoptosis significantly, proving the pivotal role of these elements in driving cell death. The effectiveness of this new approach was confirmed in ovarian tumour xenograft studies, revealing a promising avenue for future cancer therapies.

RNA-Seq Transcriptome Analysis Reveals Long Terminal Repeat Retrotransposon Modulation in Human Peripheral Blood Mononuclear Cells after In Vivo Lipopolysaccharide Injection

Human endogenous retroviruses (HERVs) and mammalian apparent long terminal repeat (LTR) retrotransposons (MaLRs) are retroviral sequences that integrated into germ line cells millions of years ago. Transcripts of these LTR retrotransposons are present in several tissues, and their expression is modulated in pathological conditions, although their function remains often far from being understood. Here, we focused on the HERV/MaLR expression and modulation in a scenario of immune system activation. We used a public data set of human peripheral blood mononuclear cells (PBMCs) RNA-Seq from 15 healthy participants to a clinical trial before and after exposure to lipopolysaccharide (LPS), for which we established an RNA-Seq workflow for the identification of expressed and modulated cellular genes and LTR retrotransposon elements.IMPORTANCE We described the HERV and MaLR transcriptome in PBMCs, finding that about 8.4% of the LTR retrotransposon loci were expressed and identifying the betaretrovirus-like HERVs as those with the highest percentage of expressed loci. We found 4,607 HERV and MaLR loci that were modulated as a result of in vivo stimulation with LPS. The HERV-H group showed the highest number of differentially expressed most intact proviruses. We characterized the HERV and MaLR loci as differentially expressed, checking their genomic context of insertion and observing a general colocalization with genes that are involved and modulated in the immune response, as a consequence of LPS stimulation. The analyses of HERV and MaLR expression and modulation show that these LTR retrotransposons are expressed in PBMCs and regulated in inflammatory settings. The similar regulation of HERVs/MaLRs and genes after LPS stimulation suggests possible interactions of LTR retrotransposons and the immune host response.

Contribution of Human Retroviruses to Disease Development-A Focus on the HIV- and HERV-Cancer Relationships and Treatment Strategies

Animal retroviruses are known for their transforming potential, and this is also true for the ones hosted by humans, which have gathered expanding attention as one of the potent causative agents in various disease, including specific cancer types. For instance, Human T Lymphotropic virus (HTLV) is a well-studied class of oncoviruses causing T cell leukemia, while human immunodeficiency virus (HIV) leads to acquired immunodeficiency syndrome (AIDS), which is linked to a series of defining cancers including Kaposi sarcoma, certain types of non-Hodgkin lymphoma, and cervical cancer. Of note, in addition to these "modern" exogenous retroviruses, our genome harbors a staggering number of human endogenous retroviruses (HERVs). HERVs are the genetic remnants of ancient retroviral germline infection of human ancestors and are typically silenced in normal tissues due to inactivating mutations and sequence loss. While some HERV elements have been appropriated and contribute to human physiological functions, others can be reactivated through epigenetic dysregulations to express retroviral elements and promote carcinogenesis. Conversely, HERV replication intermediates or protein products can also serve as intrinsic pathogen-associated molecular patterns that cause the immune system to interpret it as an exogenous infection, thereby stimulating immune responses against tumors. As such, HERVs have also been targeted as a potential internal strategy to sensitize tumor cells for promising immunotherapies. In this review, we discuss the dynamic role of human retroviruses in cancer development, focusing on HIV and HERVs contribution. We also describe potential treatment strategies, including immunotherapeutic targeting of HERVs, inhibiting DNA methylation to expose HERV signatures, and the use of antiretroviral drugs against HIV and HERVs, which can be employed as prospective anti-cancer modalities.

A human MMTV-like betaretrovirus linked to breast cancer has been present in humans at least since the copper age

The betaretrovirus Mouse Mammary Tumor Virus (MMTV) is the well characterized etiological agent of mammary tumors in mice. In contrast, the etiology of sporadic human breast cancer (BC) is unknown, but accumulating data indicate a possible viral origin also for these malignancies. The presence of MMTVenv-like sequences (MMTVels) in the human salivary glands and saliva supports the latter as possible route of inter-human dissemination. In the absence of the demonstration of a mouse-man transmission of MMTV, we considered the possibility that a cross-species transmission could have occurred in ancient times. Therefore, we investigated MMTVels in the ancient dental calculus, which originates from saliva and is an excellent material for paleovirology. The calculus was collected from 36 ancient human skulls, excluding any possible mouse contamination. MMTV-like sequences were identified in the calculus of 6 individuals dated from the Copper Age to the 17^th century. The MMTV-like sequences were compared with known human endogenous betaretroviruses and with animal exogenous betaretroviruses, confirming their exogenous origin and relation to MMTV. These data reveal that a human exogenous betaretrovirus similar to MMTV has existed at least since 4,500 years ago and indirectly support the hypothesis that it could play a role in human breast cancer.

Human Endogenous Retrovirus K in Cancer: A Potential Biomarker and Immunotherapeutic Target

In diseases where epigenetic mechanisms are changed, such as cancer, many genes show altered gene expression and inhibited genes become activated. Human endogenous retrovirus type K (HERV-K) expression is usually inhibited in normal cells from healthy adults. In tumor cells, however, HERV-K mRNA expression has been frequently documented to increase. Importantly, HERV-K-derived proteins can act as tumor-specific antigens, a class of neoantigens, and induce immune responses in different types of cancer. In this review, we describe the function of the HERV-K HML-2 subtype in carcinogenesis as biomarkers, and their potential as targets for cancer immunotherapy.

High-Throughput Sequencing is a Crucial Tool to Investigate the Contribution of Human Endogenous Retroviruses (HERVs) to Human Biology and Development

Human Endogenous retroviruses (HERVs) are remnants of ancient retroviral infections that represent a large fraction of our genome. Their transcriptional activity is finely regulated in early developmental stages and their expression is modulated in different cell types and tissues. Such activity has an impact on human physiology and pathology that is only partially understood up to date. Novel high-throughput sequencing tools have recently allowed for a great advancement in elucidating the various HERV expression patterns in different tissues as well as the mechanisms controlling their transcription, and overall, have helped in gaining better insights in an all-inclusive understanding of the impact of HERVs in biology of the host.

Comprehensive Analysis of HERV Transcriptome in HIV+ Cells: Absence of HML2 Activation and General Downregulation of Individual HERV Loci

Human endogenous retrovirus (HERV) expression is currently studied for its possible activation by HIV infection. In this context, the HERV-K(HML2) group is the most investigated: it has been proposed that HIV-1 infection can prompt HML2 transcription, and that HML2 proteins can affect HIV-1 replication, either complementing HIV or possibly influencing antiretroviral therapy. However, little information is available on the expression of other HERV groups in HIV infection. In the present study, we used a bioinformatics pipeline to investigate the transcriptional modulation of approximately 3250 well-characterized HERV loci, comparing their expression in a public RNA-seq profile, including a HIV-1-infected and a control T cell culture. In our pilot study, we found approximately 200 HERV loci belonging to 35 HERV groups that were expressed in one or both conditions, with transcripts per million (TPM) values from 1 to >500. Intriguingly, HML2 elements constituted only the 3% of expressed HERV loci, and in most cases (160) HERV expression was downregulated in the HIV-infected culture, showing from a 1- to 14-fold decrease as compared to uninfected cells. HERV transcriptome has been inferred de novo and employed to predict a total of about 950 HERV open reading frames (ORFs). These have been validated according to the coding potential and estimated abundance of the corresponding transcripts, leading to a set of 57 putative proteins potentially encoded by 23 HERV loci. Analysis showed that some individual loci have a coding potential that deserves further investigation. Among them, a HML6 provirus at locus 19q13.43 was predicted to produce a transcript showing the highest TPM among HERV-derived transcripts, being upregulated in HIV+ cells and inferred to produce Gag and Env puteins with possible biological activity.

Identification and characterization of ERV-W-like sequences in Platyrrhini species provides new insights into the evolutionary history of ERV-W in primates

Background

Endogenous Retroviruses (ERVs) constitute approximately 8% of every human genome and are relics of ancestral infections that affected the germ line cells. The ERV-W group contributed to primate physiology by providing an envelope protein (Syncytin-1) that has been adopted for placenta development in hominoids. Expression of Human ERV-W (HERV-W) sequences is investigated for a pathological role in various human diseases.

Results

We previously characterized ERV-W group genomic sequences in human and non-human Catarrhini species. We now investigated ERV-W-like sequences in the parvorder Platyrrhini, especially regarding two species with complete genome assemblies, namely marmoset (Callithrix jacchus) and squirrel monkey (Saimiri boliviensis). We identified in both species proviral sequences, annotated as ERV1–1 in respective genome assemblies, sharing high sequence similarities with Catarrhini ERV-W. A total of 130 relatively intact proviruses from the genomes of marmoset and squirrel monkey were characterized regarding their structural and evolutionarily relationships with Catarrhini ERV-W elements. Platyrrhini ERV-W sequences share several structural features with Catarrhini ERV-W elements and are closely related phylogenetically with the latter as well as with other ERV-W-related gammaretrovirus-like ERVs. The ERV-W group colonized Platyrrhini primates of both Callitrichidae and Atelidae lineages, with provirus formations having occurred mostly between 25 and 15 mya. Two LTR subgroups were associated with monophyletic proviral bodies. A pre-gag region appears to be a sequence feature common to the ERV-W group: it harbors a putative intron sequence that is missing in some ERV-W loci, holding a putative ORF as well. The presence of a long pre-gag portion was confirmed among all gammaretroviral ERV analyzed, suggesting a role in the latter biology. It is noteworthy that, contrary to Catarrhini ERV-W, there was no evidence of L1-mediated mobilization for Platyrrhini ERV-W sequences.

Conclusions

Our data establish that ERV-W is not exclusive to Catarrhini primates but colonized both parvorders of Simiiformes, providing further insight into the evolution of ERV-W and the colonization of primate genomes.