Could the Human Endogenous Retrovirus-Derived Syncytialization Inhibitor, Suppressyn, Limit Heterotypic Cell Fusion Events in the Decidua?

Proper placental development relies on tightly regulated trophoblast differentiation and interaction with maternal cells. Human endogenous retroviruses (HERVs) play an integral role in modulating cell fusion events in the trophoblast cells of the developing placenta. Syncytin-1 (ERVW-1) and its receptor, solute-linked carrier family A member 5 (SLC1A5/ASCT2), promote fusion of cytotrophoblast (CTB) cells to generate the multi-nucleated syncytiotrophoblast (STB) layer which is in direct contact with maternal blood. Another HERV-derived protein known as Suppressyn (ERVH48-1/SUPYN) is implicated in anti-fusogenic events as it shares the common receptor with ERVW-1. Here, we explore primary tissue and publicly available datasets to determine the distribution of ERVW-1, ERVH48-1 and SLC1A5 expression at the maternal-fetal interface. While SLC1A5 is broadly expressed in placental and decidual cell types, ERVW-1 and ERVH48-1 are confined to trophoblast cell types. ERVH48-1 displays higher expression levels in CTB and extravillous trophoblast, than in STB, while ERVW-1 is generally highest in STB. We have demonstrated through gene targeting studies that suppressyn has the ability to prevent ERVW-1-induced fusion events in co-culture models of trophoblast cell/maternal endometrial cell interactions. These findings suggest that differential HERV expression is vital to control fusion and anti-fusogenic events in the placenta and consequently, any imbalance or dysregulation in HERV expression may contribute to adverse pregnancy outcomes.

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

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.

HERVs establish a distinct molecular subtype in stage II/III colorectal cancer with poor outcome

Colorectal cancer (CRC) is one of the most lethal malignancies. The extreme heterogeneity in survival rate is driving the need for new prognostic biomarkers. Human endogenous retroviruses (hERVs) have been suggested to influence tumor progression, oncogenesis and elicit an immune response. We examined multiple next-generation sequencing (NGS)-derived biomarkers in 114 CRC patients with paired whole-exome and whole-transcriptome sequencing (WES and WTS, respectively). First, we demonstrate that the median expression of hERVs can serve as a potential biomarker for prognosis, relapse, and resistance to chemotherapy in stage II and III CRC. We show that hERV expression and CD8+ tumor-infiltrating T-lymphocytes (TILs) synergistically stratify overall and relapse-free survival (OS and RFS): the median OS of the CD8-/hERV+ subgroup was 29.8 months compared with 37.5 months for other subgroups (HR = 4.4, log-rank P < 0.001). Combing NGS-based biomarkers (hERV/CD8 status) with clinicopathological factors provided a better prediction of patient survival compared to clinicopathological factors alone. Moreover, we explored the association between genomic and transcriptomic features of tumors with high hERV expression and establish this subtype as distinct from previously described consensus molecular subtypes of CRC. Overall, our results underscore a previously unknown role for hERVs in leading to a more aggressive subtype of CRC.

Unbiased examination of genome-wide human endogenous retrovirus transcripts in MS brain lesions

Background: Human endogenous retrovirus (HERV) expression in multiple sclerosis (MS) brain lesions may contribute to chronic inflammation, but expression of genome-wide HERVs in different MS lesions is unknown. Objective: We examined the HERV expression landscape in different MS lesions compared to control brains. Methods: Transcripts from 71 MS brain samples and 25 control WM were obtained by next-generation RNA sequencing and mapped against HERV transcripts across the human genome. Differential expression of mapped HERV-W and HERV-H reads between MS lesion types and controls was analysed. Results: Out of 6.38 billion high-quality paired end reads, 174 million reads (2.73%) mapped to HERV transcripts. There was no difference in HERVs expression level between MS and control brains, but HERV-W transcripts were significantly reduced in chronic active lesions. Of the four HERV-W transcripts exclusively present in MS, ERV3633503 located on chromosome 7q21.13 close to the MS genetic risk locus had the highest number of reads. In the HERV-H family, 75% of transcripts located to nearby 7q21-22 were overrepresented in MS, and ERV3643914 was expressed more than 16 times in MS compared to control brains. Conclusion: Novel HERV-W and HERV-H transcripts located at chromosome 7 regions were uniquely expressed in MS lesions, indicating their potential role in brain lesion evolution.

Human Endogenous Retrovirus Expression Is Associated with Head and Neck Cancer and Differential Survival

Human endogenous retroviruses (HERVs) have been implicated in a variety of human diseases including cancers. However, technical challenges in analyzing HERV sequence data have limited locus-specific characterization of HERV expression. Here, we use the software Telescope (developed to identify expressed transposable elements from metatranscriptomic data) on 43 paired tumor and adjacent normal tissue samples from The Cancer Genome Atlas Program to produce the first locus-specific retrotranscriptome of head and neck cancer. Telescope identified over 3000 expressed HERVs in tumor and adjacent normal tissue, and 1078 HERVs were differentially expressed between the two tissue types. The majority of differentially expressed HERVs were expressed at a higher level in tumor tissue. Differentially expressed HERVs were enriched in members of the HERVH family. Hierarchical clustering based on HERV expression in tumor-adjacent normal tissue resulted in two distinct clusters with significantly different survival probability. Together, these results highlight the importance of future work on the role of HERVs across a range of cancers.

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.

Activation of Transposable Elements in Immune Cells of Fibromyalgia Patients

Advancements in nucleic acid sequencing technology combined with an unprecedented availability of metadata have revealed that 45% of the human genome constituted by transposable elements (TEs) is not only transcriptionally active but also physiologically necessary. Dysregulation of TEs, including human retroviral endogenous sequences (HERVs) has been shown to associate with several neurologic and autoimmune diseases, including Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS). However, no study has yet addressed whether abnormal expression of these sequences correlates with fibromyalgia (FM), a disease frequently comorbid with ME/CFS. The work presented here shows, for the first time, that, in fact, HERVs of the H, K and W types are overexpressed in immune cells of FM patients with or without comorbid ME/CFS. Patients with increased HERV expression (N = 14) presented increased levels of interferon (INF-β and INF-γ) but unchanged levels of TNF-α. The findings reported in this study could explain the flu-like symptoms FM patients present with in clinical practice, in the absence of concomitant infections. Future work aimed at identifying specific genomic loci differentially affected in FM and/or ME/CFS is warranted.

Activation of Transposable Elements in Immune Cells of Fibromyalgia Patients

The development of nucleic acid sequencing technology and the unprecedented availability of metadata has evidenced that 45% of human genome constituted by transposable elements (TEs) is not only transcriptionally active but also physiologically needed. Aberrant regulation of TEs, and of human retroviral endogenous sequences (HERVs) in particular, associates with several neurologic and autoimmune diseases, including the Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) frequently comorbid with fibromyalgia (FM). However, no study has yet addressed whether abnormal expression of these sequences correlates with FM. The work presented here shows, for the first time, that in fact HERVs of the H, K and W types are overexpressed in the cells of the immune system of FM patients with or without comorbid ME/CFS. The patients with increased HERV expression (N=14) presented increased levels of interferon (INF-&beta; and INF-&gamma;) but unchanged levels of TNF-&alpha;. In support of our proposal that TE activation is a contributor to FM, we find that the tRNA pools are decreased in comparison to matched healthy participants (N=14). The findings reported here could explain the flu-like symptoms FM patients present with in the absence of concomitant infections. Future work towards identifying specific genomic loci differentially affected in FM and ME/CFS is granted.

Telescope: Characterization of the retrotranscriptome by accurate estimation of transposable element expression

Characterization of Human Endogenous Retrovirus (HERV) expression within the transcriptomic landscape using RNA-seq is complicated by uncertainty in fragment assignment because of sequence similarity. We present Telescope, a computational software tool that provides accurate estimation of transposable element expression (retrotranscriptome) resolved to specific genomic locations. Telescope directly addresses uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We demonstrate the utility of our approach through single locus analysis of HERV expression in 13 ENCODE cell types. When examined at this resolution, we find that the magnitude and breadth of the retrotranscriptome can be vastly different among cell types. Furthermore, our approach is robust to differences in sequencing technology, and demonstrates that the retrotranscriptome has potential to be used for cell type identification. Telescope performs highly accurate quantification of the retrotranscriptomic landscape in RNA-seq experiments, revealing a differential complexity in the transposable element biology of complex systems not previously observed. Telescope is available at github.com/mlbendall/telescope.Author SummaryAlmost half of the human genome is composed of Transposable elements (TEs), but their contribution to the transcriptome, their cell-type specific expression patterns, and their role in disease remains poorly understood. Recent studies have found many elements to be actively expressed and involved in key cellular processes. For example, human endogenous retroviruses (HERVs) are reported to be involved in human embryonic stem cell differentiation. Discovering which exact HERVs are differentially expressed in RNA-seq data would be a major advance in understanding such processes. However, because HERVs have a high level of sequence similarity it is hard to identify which exact HERV is differentially expressed. To solve this problem, we developed a computer program which addressed uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We call this program, “Telescope”. We then used Telescope to identify HERV expression in 13 well-studied cell types from the ENCODE consortium and found that different cell types could be characterized by enrichment for different HERV families, and for locus specific expression. We also showed that Telescope performed better than other methods currently used to determine TE expression. The use of this computational tool to examine new and existing RNA-seq data sets may lead to new understanding of the roles of TEs in health and disease.