Human Endogenous Retroviruses Loci Research Articles

Limitations of current high-throughput sequencing technologies lead to biased expression estimates of endogenous retroviral elements.

Human endogenous retroviruses (HERVs), the remnants of ancient germline retroviral integrations, comprise almost 8% of the human genome. The elucidation of their biological roles is hampered by our inability to link HERV mRNA and protein production with specific HERV loci. To solve the riddle of the integration-specific RNA expression of HERVs, several bioinformatics approaches have been proposed; however, no single process seems to yield optimal results due to the repetitiveness of HERV integrations. The performance of existing data-bioinformatics pipelines has been evaluated against real world datasets whose true expression profile is unknown, thus the accuracy of widely-used approaches remains unclear. Here, we simulated mRNA production from specific HERV integrations to evaluate second and third generation sequencing technologies along with widely used bioinformatic approaches to estimate the accuracy in describing integration-specific expression. We demonstrate that, while a HERV-family approach offers accurate results, per-integration analyses of HERV expression suffer from substantial expression bias, which is only partially mitigated by algorithms developed for calculating the per-integration HERV expression, and is more pronounced in recent integrations. Hence, this bias could erroneously result into biologically meaningful inferences. Finally, we demonstrate the merits of accurate long-read high-throughput sequencing technologies in the resolution of per-locus HERV expression.

A transposable element atlas of aggressive B‐cell non‐Hodgkin lymphomas defines novel classifications of Burkitt lymphoma independent of EBV status

Introduction: Clinical and pathological heterogeneity in aggressive B-cell non-Hodgkin lymphomas may be driven by several mechanisms of oncogenesis, such as translocations involving the immunoglobulin locus and proto-oncogenes, the ‘cell-of-origin’ from the germinal center (GC), and viral associations. Previous classifications of B-cell lymphomas have been defined by one or more of these characteristics and transforming events. However, transposable elements (TEs), which make up roughly half of the human genome and have key regulatory functions, have been previously overlooked in both the healthy GC and in B-cell lymphomas. We hypothesized that unique TE signatures, which may in part be driven by human endogenous retroviruses (HERVs), would allow a novel categorization of B-cell lymphomas. Here we present a comprehensive, locus-specific atlas of TE expression in healthy GC B cells, diffuse large B-cell lymphoma (DLBCL), Epstein-Barr Virus (EBV)-positive and negative Burkitt lymphoma (BL), and follicular lymphoma (FL). Methods: We obtained RNA-seq data from 529 DLBCLs belonging to the TCGA and NCICCR cohorts, 113 sporadic and endemic BL, and 12 FL from the CGCI cohort for gene and locus-specific TE quantification. RNA-seq data from sorted cells from the healthy GC were obtained from two recent studies (Holmes et al., 2020; Agirre et al., 2019). TEs were quantified with Telescope, which addresses the computational challenges of repetitive and interspersed TE reads, and low read counts using an expectation-maximization algorithm. We used STAR for indexing and alignment, HTSeq for quantifying EBV viral reads in endemic BL samples, DESeq2 to identify differentially expressed TEs, and LRT, Boruta, and Lasso for feature selection. Results: TE-driven clustering of healthy B cells showed that HERVs could independently distinguish stages of B cell differentiation and specific GC subsets, with plasmablasts and bone marrow plasma cells having the highest number of differentially expressed HERVs. We found HERV loci of interest upstream of POU5F1B and MYC that can differentiate plasmablasts from other GC cells. Furthermore, we describe a map of locus-specific TE expression in DLBCL, BL, and FL. Strikingly, we find that BL can be subdivided into three HERV-driven clusters, which are not obtained with gene-only clustering, and identify 4 HERV loci sufficient to distinguish between clusters. Clusters are independent of EBV status, clinical variant, location of MYC translocation, and display a high number of differentially expressed LINE elements, lncRNAs, and snoRNAs. Conclusions: Here, we report an atlas of TE expression in BL, DLBCL, and FL, and provide evidence for disease-specific changes in TE expression in B-cell lymphomas. TEs that are selectively expressed in lymphoma subtypes provide opportunities for novel subclassification that may have biological and pathological relevance. The research was funded by: The work was supported in part by US National Institutes of Health (NIH) grant CA260691 (DFN), and UM1AI164559 (DFN). MLB is supported in part by the Department of Medicine, Fund for the Future program at Weill Cornell Medicine sponsored by the Elsa Miller Foundation. J.L.M. was supported in part by a Medical Scientist Training Program grant to the Weill Cornell–Rockefeller–Sloan Kettering Tri-Institutional MD-PhD Program (T32GM007739). Keywords: Aggressive B-cell non-Hodgkin lymphoma, Bioinformatics, Computational and Systems Biology, Pathology and Classification of Lymphomas No conflicts of interests pertinent to the abstract.

Epstein-Barr Virus-Induced Genes and Endogenous Retroviruses in Immortalized B Cells from Patients with Multiple Sclerosis.

The immune pathogenesis of multiple sclerosis (MS) is thought to be triggered by environmental factors in individuals with an unfavorable genetic predisposition. Epstein-Barr virus (EBV) infection is a major risk factor for subsequent development of MS. Human endogenous retroviruses (HERVs) can be activated by EBV, and might be a missing link between an initial EBV infection and the later onset of MS. In this study, we investigated differential gene expression patterns in EBV-immortalized lymphoblastoid B cell lines (LCL) from MS-affected individuals (MSLCL) and controls by using RNAseq and qRT-PCR. RNAseq data from LCL mapped to the human genome and a virtual virus metagenome were used to identify possible biomarkers for MS or disease-relevant risk factors, e.g., the relapse rate. We observed that lytic EBNA-1 transcripts seemed to be negatively correlated with age leading to an increased expression in LCL from younger PBMC donors. Further, HERV-K (HML-2) GAG was increased upon EBV-triggered immortalization. Besides the well-known transactivation of HERV-K18, our results suggest that another six HERV loci are up-regulated upon stimulation with EBV. We identified differentially expressed genes in MSLCL, e.g., several HERV-K loci, ERVMER61-1 and ERV3-1, as well as genes associated with relapses. In summary, EBV induces genes and HERV in LCL that might be suitable as biomarkers for MS or the relapse risk.

EPCO-33. DIFFERENTIAL EXPRESSION OF AN ENDOGENOUS RETROVIRAL ELEMENT [HERV-K(HML-6)] IS ASSOCIATED WITH REDUCED SURVIVAL IN GLIOBLASTOMA PATIENTS

Abstract INTRODUCTION 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.Method:A custom workflow was used to quantify HERV expression in our cell lines of interest. Cell-line methylation was quantified using a custom script. We generated primers specific for the Human endogenous Mouse mammary tumor (MMTV)-Like virus 6 (HML-6). Visualization of RNA transcripts was performed using RNA-scope. Clinical data was obtained using the R package, TCGABiolinks. RESULTS The A172 cell line had significantly higher mean overall HERV expression relative to the M059J and H4 cell lines (p< 0.0001 for both). A172 cells had significantly lower mean number of CpG islands relative to M059J cells and H4 cells (p< 0.0001 for both). There was a significant inverse correlation between mean beta value and FC HERV expression (R=-0.57, p=0.01). qPCR confirmed robust expression of the HML-6 locus in cell culture and neurospheres. Elevated ERVK3-1 expression was associated reduced survival among IDHwt GBM patients (18.3 vs. 15.1 months, p=0.039). This was preserved among IDH mutant (IDHm) GBM as well (17.9 months vs. 14.0 months, p=0.0088). CONCLUSION In gliomas, HERV expression correlates with loss of DNA methylation at HERV loci. HML-6 is overexpressed in highly invasive glioblastoma cell lines and patient-derived neurospheres. We have demonstrated a potential survival detriment associated with elevated expression of the HML-6 product, ERVK3-1.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV’s roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species’ genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates’ germlines after the segregation of Gibbons’ species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold–fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells’ populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein–protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer’s, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells’ differentiation programs.

Specific human endogenous retroviruses predict metastatic potential in uveal melanoma.

Uveal melanoma (UM) is a unique disease in that patients with primary UM are well stratified based on their risk of developing metastasis, yet there are limited effective treatments once metastases occur. There is an urgent need to better understand the distinct molecular pathogenesis of UM and the characteristics of patients at high risk for metastasis to identify neoantigenic targets that can be used in immunotherapy and to develop novel therapeutic strategies that may effectively target this lethal transition. An important and overlooked area of molecular pathogenesis and neoantigenic targets in UM comes from human endogenous retroviruses (HERVs). We investigated the HERV expression landscape in primary UM and found that tumors were stratified into 4 HERV-based subsets that provide clear delineation of risk outcome and support subtypes identified by other molecular indicators. Specific HERV loci are associated with the risk of uveal melanoma metastasis and may offer mechanistic insights into this process, including dysregulation of HERVs on chromosomes 3 and 8. A HERV signature composed of 17 loci was sufficient to classify tumors according to subtype with greater than 95% accuracy, including at least 1 intergenic HERV with coding potential (HERVE_Xp11.23) that could represent a potential HERV E target for immunotherapy.

Transcriptome Analyses Implicate Endogenous Retroviruses Involved in the Host Antiviral Immune System through the Interferon Pathway

Human endogenous retroviruses (HERVs) are the remains of ancient retroviruses that invaded our ancestors’ germline cell and were integrated into the genome. The expression of HERVs has always been a cause for concern because of its association with various cancers and diseases. However, few previous studies have focused on specific activation of HERVs by viral infections. Our previous study has shown that dengue virus type 2 (DENV-2) infection induces the transcription of a large number of abnormal HERVs loci; therefore, the purpose of this study was to explore the relationship between exogenous viral infection and HERV activation further. In this study, we retrieved and reanalyzed published data on 21 transcriptomes of human cells infected with various viruses. We found that infection with different viruses could induce transcriptional activation of HERV loci. Through the comparative analysis of all viral datasets, we identified 43 key HERV loci that were up-regulated by DENV-2, influenza A virus, influenza B virus, Zika virus, measles virus, and West Nile virus infections. Furthermore, the neighboring genes of these HERVs were simultaneously up-regulated, and almost all such neighboring genes were interferon-stimulated genes (ISGs), which are enriched in the host’s antiviral immune response pathways. Our data supported the hypothesis that activation of HERVs, probably via an interferon-mediated mechanism, plays an important role in innate immunity against viral infections.

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

Simple SummaryThe 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.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.

Simple SummaryHuman Endogenous retroviruses (HERVs) and Mammalian Apparent LTRs Retrotransposons (MaLRs) are remnants of ancient retroviral infections that make up about 8% of the human genome. HERV and MaLR expression is regulated in immunity, and, in particular, is known for their up-regulation after innate immune activation. In this work, we analyzed the differential expression of HERVs and MaLRs in different stages of immunity, triggered by the administration of an inactivated vaccine.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.

Human endogenous retroviruses form a reservoir of T cell targets in hematological cancers

Human endogenous retroviruses (HERV) form a substantial part of the human genome, but mostly remain transcriptionally silent under strict epigenetic regulation, yet can potentially be reactivated by malignant transformation or epigenetic therapies. Here, we evaluate the potential for T cell recognition of HERV elements in myeloid malignancies by mapping transcribed HERV genes and generating a library of 1169 potential antigenic HERV-derived peptides predicted for presentation by 4 HLA class I molecules. Using DNA barcode-labeled MHC-I multimers, we find CD8+ T cell populations recognizing 29 HERV-derived peptides representing 18 different HERV loci, of which HERVH-5, HERVW-1, and HERVE-3 have more profound responses; such HERV-specific T cells are present in 17 of the 34 patients, but less frequently in healthy donors. Transcriptomic analyses reveal enhanced transcription of the HERVs in patients; meanwhile DNA-demethylating therapy causes a small and heterogeneous enhancement in HERV transcription without altering T cell recognition. Our study thus uncovers T cell recognition of HERVs in myeloid malignancies, thereby implicating HERVs as potential targets for immunotherapeutic therapies.

Identification of human endogenous retrovirus transcripts in Hodgkin Lymphoma cells.

During the last decades, the prognosis for patients with Hodgkin Lymphoma (HL) has been steadily improved. Nevertheless, new and less toxic therapy strategies have to be developed especially for patients with advanced disease. The activation of human endogenous retroviruses (HERV) is suspected to occur in HL and therefore, HERV might represent interesting target structures. In order to identify transcribed HERV of the HERV-H and HERV-K families in HL we used a reverse transcription-polymerase chain reaction based cloning approach. In addition to unspliced HERV-H and HERV-K transcripts, we detected spliced HERV-K transcripts that matched genomic sequences with the expected splicing-donor and splicing-acceptor sites. Of particular interest was the expression of HERV-K18 related transcripts at the CD48 locus. Our data indicate transcriptional activity of several HERV loci in HL cells.

Genetic Evidence for Involvement of Human Endogenous Retrovirus Herv-Fc1 in the Pathogenesis of MPNs

Introduction: Human endogenous retroviruses (HERVs) are genomic sequences that result from ancestral germ-line infections by exogenous RVs and therefore are transmitted in a Mendelian fashion. Several lines of genetic evidence suggest the involvement of HERVs in the etiology of autoimmune diseases and cancer. Here, by means of genetic epidemiology, we describe a case-control association study of SNPs located in the proximity of HERVs elements in patients with myeloproliferative neoplasms (MPNs).Methods: Genomic DNA from two cohorts recruited from 4 university hospitals, in total 376 patients with MPN and 527 control subjects, was isolated from whole blood. 148 SNPs related to 49 HERV loci were tested using the Sequenom platform. To quantify genomic DNA copy number variation, real-time PCR was performed using the LightCycler system. Data was analyzed using the χ2-test. For synergy, we recoded controls and cases as 0 and 1, respectively, and used ANOVA with type III sum of squares. Adjustment for multiple hypothesis testing was done using the Bonferroni correction (pB <0.05).Results: In cohort 1, analysis of 148 SNPs in 188 MPN patients and 453 controls demonstrated a significant association in MPNs of the two SNPs rs219077 and rs5993426 belonging to HERV-K and HERV-9, respectively. The two SNPs were validated in an additional cohort of 188 MPNs and were again highly significant, rs219077 (OR 1.97 (1.43 - 2.7), combined pB = 0.005) and rs5993426 (OR 3.94 (2.72 - 5.71), combined pB = 9x10-11). Additionally, we report that in patients carrying the SNP risk allele (G), rs219077, qPCR analysis results in significantly higher copy numbers compared to patients not carrying the risk allele.Discussion and conclusion: About 30 years ago HERV-K particles were reported in megakaryocytes cultured from patients with ET. Surprisingly, since then no studies have focused upon a potential link between virus infection and the development of MPNs. In the present study, we revive the HERV-story in MPNs. Thus, we show that two SNPs, rs5993426 and rs219077, close to HERVs correlate with risk of MPNs supporting a role of these genetics elements in MPN pathogenesis. In our study, rs219077, located on the chromosome 2, associates with the disease and maps in the proximity of the HERV-9 sequence. Additionally, the SNP risk allele (G) in rs219077 results in significantly higher copy numbers. There are two possible explanations for the amplification of rs219077 in MPNs. First, it might be an amplification of a part of chromosome 2, as genomic features of MPNs often include clonal, non-random chromosomal amplifications. Second, it might be due to the formation of a transducing virus from the right LTR of HERV-9, possibly now (after a few rounds of replication) with HERV-K LTRs. The observed amplification of risk allele (G) in MPN may indicate that the recombinant virus pathogenic copy can retrotranspose in PBMCs with an active viral reverse transcriptase and, eventually, integrate through recombination/integration events. rs219077 is also located in genetic proximity to MIR3 SINE elements. Thus, one fascinating possibility is that the putative transducing ERV-9 virus could encode or regulate microRNA. The second marker showing significant association with MPN risk was rs5993426. This marker in the proximity of HERV-K encoding Gag, Pol has been reported to be associated with rheumatoid arthritis. Our genetic data are currently exclusively based on genetic association and cannot as such predict any functional relationship beyond presence or absence of SNPs. Nevertheless, statistical associations are most easily explained by functional relations and it is therefore pertinent to pursue the potential molecular mechanisms mediated by HERVs in MPNs. As an example, the induction and progression of rodent erythroleukemias by Friend spleen-forming virus are mainly due to the ability of proviruses to activate cellular oncogenes or inactivate tumor suppressor genes. In this regard, we have most recently shown downregulation of TP53 in MPNs to be significantly upregulated during treatment with interferon-alpha2. In conclusion, we have for the first time demonstrated a genetic association of MPNs with HERV elements. The association supports a pathogenetic role of these loci in MPNs. Our observations stimulate to further studies on HERVs in the pathogenesis of MPNs and our results need to be confirmed in larger cohorts of patients. DisclosuresHasselbalch:Novartis: Research Funding.

Human endogenous retroviruses in the aetiology of MS.

Several lines of investigation have provided strong indications for an association between the immune-mediated, neurologic disease multiple sclerosis (MS) and human endogenous retroviruses (HERVs). Whether the relationship is causal is yet to be established. Endogenous retroviruses are pathogenic-in other species than the human. Several aspects of the activation and involvement of specific HERV families (HERV-H/F and HERV-W/MSRV) have been documented, both for cells in the periphery and in the central nervous system. Specific HERV-encoded genes and certain gene products (envelope proteins, Envs) appear strongly associated with the disease and have pathogenic potential. Most HERV sequences are non-functional, whereas some HERV loci have coding potential but remain quiescent in non-pathological conditions, so the importance of regulatory pathways and epigenetics involved in regulating HERV activation, de-repression, and also involvement of retroviral restriction factors, is emerging. Disease intervention by means of antiretrovirals has potential as a novel therapeutic strategy in MS treatment; this is compounded by the apparently reduced risk of MS in HIV infection as a consequence of therapy. Extensive studies of HERVs, their role in neurologic diseases, and their potential as therapeutic targets are needed.

IFITM1 suppresses expression of human endogenous retroviruses in human embryonic stem cells.

Interferon‐induced transmembrane protein 1 (IFITM1), a member of the IFITM protein family, is a component of a multimeric complex involved in the transduction of antiproliferation and cell adhesion signals. IFITM1 is thought to play a role in antiproliferation and immune surveillance, and has been shown to restrict infection by numerous viruses. It is highly expressed in human embryonic stem cells (hESCs) but its role in hESCs remains to be elucidated. In this study, knockout of IFITM1 mediated by CRISPR/Cas9 in hESCs did not affect self‐renewal, pluripotency, telomerase activity or telomeres. However expression of human endogenous retroviruses (HERVs) was higher than in wild‐type hESCs, and there was also a reduced level of trimethylation of histone H3 on lysine 9 at HERV loci. These data show that IFITM1 suppresses HERVs in hESCs by regulating epigenetic modifications.

Microarray-based Identification of Individual HERV Loci Expression: Application to Biomarker Discovery in Prostate Cancer

The prostate-specific antigen (PSA) is the main diagnostic biomarker for prostate cancer in clinical use, but it lacks specificity and sensitivity, particularly in low dosage values1​​. ‘How to use PSA' remains a current issue, either for diagnosis as a gray zone corresponding to a concentration in serum of 2.5-10 ng/ml which does not allow a clear differentiation to be made between cancer and noncancer2 or for patient follow-up as analysis of post-operative PSA kinetic parameters can pose considerable challenges for their practical application3,4. Alternatively, noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease, e.g. PCA3 in prostate cancer5,6 and to reveal uncharacterized aspects of tumor biology. Moreover, data from the ENCODE project published in 2012 showed that different RNA types cover about 62% of the genome. It also appears that the amount of transcriptional regulatory motifs is at least 4.5x higher than the one corresponding to protein-coding exons. Thus, long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) constitute a wide range of putative/candidate transcriptional regulatory sequences, as it is their primary function in infectious retroviruses. HERVs, which are spread throughout the human genome, originate from ancestral and independent infections within the germ line, followed by copy-paste propagation processes and leading to multicopy families occupying 8% of the human genome (note that exons span 2% of our genome). Some HERV loci still express proteins that have been associated with several pathologies including cancer7-10. We have designed a high-density microarray, in Affymetrix format, aiming to optimally characterize individual HERV loci expression, in order to better understand whether they can be active, if they drive ncRNA transcription or modulate coding gene expression. This tool has been applied in the prostate cancer field (Figure 1).

Microarray-based Identification of Individual HERV Loci Expression: Application to Biomarker Discovery in Prostate Cancer

The prostate-specific antigen (PSA) is the main diagnostic biomarker for prostate cancer in clinical use, but it lacks specificity and sensitivity, particularly in low dosage values1​​. ‘How to use PSA' remains a current issue, either for diagnosis as a gray zone corresponding to a concentration in serum of 2.5-10 ng/ml which does not allow a clear differentiation to be made between cancer and noncancer2 or for patient follow-up as analysis of post-operative PSA kinetic parameters can pose considerable challenges for their practical application3,4. Alternatively, noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease, e.g. PCA3 in prostate cancer5,6 and to reveal uncharacterized aspects of tumor biology. Moreover, data from the ENCODE project published in 2012 showed that different RNA types cover about 62% of the genome. It also appears that the amount of transcriptional regulatory motifs is at least 4.5x higher than the one corresponding to protein-coding exons. Thus, long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) constitute a wide range of putative/candidate transcriptional regulatory sequences, as it is their primary function in infectious retroviruses. HERVs, which are spread throughout the human genome, originate from ancestral and independent infections within the germ line, followed by copy-paste propagation processes and leading to multicopy families occupying 8% of the human genome (note that exons span 2% of our genome). Some HERV loci still express proteins that have been associated with several pathologies including cancer7-10. We have designed a high-density microarray, in Affymetrix format, aiming to optimally characterize individual HERV loci expression, in order to better understand whether they can be active, if they drive ncRNA transcription or modulate coding gene expression. This tool has been applied in the prostate cancer field (Figure 1).

Identification and association study with lung cancer for novel insertion polymorphisms of human endogenous retrovirus

Sequences of human endogenous retroviruses (HERVs) are members of the long terminal repeat (LTR) retrotransposon family. Although the expression of HERV has long been a topic of investigation, HERV-insertion polymorphisms are not well known, and a genetic association between HERV-insertion polymorphisms and cancer has never been reported. To identify novel HERV loci in the genome from cancer tissues, we carried out the inverse PCR method targeting a conserved LTR region of HML-2, which is the most recently acquired HERV group. Novel two insertions, HML-2_sLTR(1p13.2) and HML-2_sLTR(19q12), were identified as insertionally polymorphic solo LTRs. Furthermore, a significant prevalence of HML-2_sLTR(1p13.2) homozygosity was detected in female never-smoking patients aged 60 years and over who had lung adenocarcinoma [versus the other genotyping; odds ratio (OR): 1.97; 95% confidence interval (CI): 1.01-3.81]. In another cohort consisting of female never-smoking patients with lung adenocarcinoma, a prevalence of HML-2_sLTR(1p13.2) homozygosity tended to be high in patients aged 60 years and over (versus the other genotyping; OR: 2.03; 95% CI: 0.96-4.29), whereas a low prevalence of HML-2_sLTR(1p13.2) homozygosity was detected in patients <60 years old (versus the other genotyping; OR: 0.31; 95% CI: 0.11-0.94). Our results suggest that HML-2_sLTR(1p13.2) is involved with the susceptibility to lung adenocarcinoma in female never-smokers in an age-dependent manner and that other HERV polymorphisms related to human diseases might remain to be identified in the human genome.

The epigenetic control of HERV loci encoding for fusogenic envelope glycoproteins in trophoblast

Background Up to 8% of the human genome is of retroviral origin, comprising numerous human endogenous retroviruses (HERVs), remnants of germ line transmission of exogenous retroviruses during primate and human evolution. Some endogenous retroviruses within the human genome remain intact open reading frames with translational potential. Gag and Env genes are often expressed in various human tissues and provide superantigens or restriction factors interfering with exogenous infections. Several Env products display fusogenic properties and it is temting to implicite them in various human pathologies like cancer, multiple sclerosis, and muscular dystrophies. Of particular importace are the fusogenic Env proteins expressed in placenta, syncytin-1 and 2, encoded by ERVWE1 and ERV-FRD proviral loci, respectively. These fusogenic Envs contribute to the differentiation of multinucleated syncytiotrophoblast in placental chorionic villi. In non-placental tissues, however, expression of these fusogenic proteins has to be suppressed because of high risk for tissue integrity and cancer development. We have previously shown that the 5' long terminal repeat, the regulatory region of syncytin-1 is epigenetically suppressed by DNA methylation in non-placental tissues.

Expression patterns of transcribed human endogenous retrovirus HERV-K(HML-2) loci in human tissues and the need for a HERV Transcriptome Project

BackgroundA significant proportion of the human genome is comprised of human endogenous retroviruses (HERVs). HERV transcripts are found in every human tissue. Expression of proviruses of the HERV-K(HML-2) family has been associated with development of human tumors, in particular germ cell tumors (GCT). Very little is known about transcriptional activity of individual HML-2 loci in human tissues, though.ResultsBy employing private nucleotide differences between loci, we assigned ~1500 HML-2 cDNAs to individual HML-2 loci, identifying, in total, 23 transcriptionally active HML-2 proviruses. Several loci are active in various human tissue types. Transcription levels of some HML-2 loci appear higher than those of other loci. Several HML-2 Rec-encoding loci are expressed in GCT and non-GCT tissues. A provirus on chromosome 22q11.21 appears strongly upregulated in pathologic GCT tissues and may explain high HML-2 Gag protein levels in GCTs. Presence of Gag and Env antibodies in GCT patients is not correlated with activation of individual loci. HML-2 proviruses previously reported capable of forming an infectious HML-2 variant are transcriptionally active in germ cell tissue. Our study furthermore shows that Expressed Sequence Tag (EST) data are insufficient to describe transcriptional activity of HML-2 and other HERV loci in tissues of interest.ConclusionOur, to date, largest-scale study reveals in greater detail expression patterns of individual HML-2 loci in human tissues of clinical interest. Moreover, large-scale, specialized studies are indicated to better comprehend transcriptional activity and regulation of HERVs. We thus emphasize the need for a specialized HERV Transcriptome Project.

Expression analyses of human endogenous retroviruses (HERVs): tissue-specific and developmental stage-dependent expression of HERVs

The evolutional and biological roles of human endogenous retroviruses (HERVs) are less recognized compared to those of L1. In the present study, we focused on the transcriptional activity of HERVs in normal human tissues and found five HERV loci that are actively expressed in normal tissues. All but one showed tissue specificity of expression: one was expressed in stomach and small intestine and three were in placenta. We subsequently examined by TaqMan-based RT-PCR assays the temporal expression profiles of the three placenta-specific HERVs along with syncytin and syncytin 2 and observed three patterns. Syncytin and HERV-Fb showed almost constant expression through gestations. Syncytin 2 gradually decreased as pregnancy proceeded. In contrast, expression from the HERV-H/F and HERV-K(HML-6) loci increased remarkably in term placentas. Term placentas in general showed larger interindividual differences in HERV expression levels. Our results suggest that HERVs might have more diverse effects than currently thought.