HERVd: database of human endogenous retroviruses

VITALdb: to select the best viroinformatics tools for a desired virus or application

The recent pandemics of viral diseases, COVID-19/mpox (humans) and lumpy skin disease (cattle), have kept us glued to viral research. These pandemics along with the recent human metapneumovirus outbreak have exposed the urgency for early diagnosis of viral infections, vaccine development, and discovery of novel antiviral drugs and therapeutics. To support this, there is an armamentarium of virus-specific computational tools that are currently available. VITALdb (VIroinformatics Tools and ALgorithms database) is a resource of ~360 viroinformatics tools encompassing all major viruses (SARS-CoV-2, influenza virus, human immunodeficiency virus, papillomavirus, herpes simplex virus, hepatitis virus, dengue virus, Ebola virus, Zika virus, etc.) and several diverse applications [structural and functional annotation, antiviral peptides development, subspecies characterization, recognition of viral recombination, inhibitors identification, phylogenetic analysis, virus-host prediction, viral metagenomics, detection of mutation(s), primer designing, etc.]. Resources, tools, and other utilities mentioned in this article will not only facilitate further developments in the realm of viroinformatics but also provide tremendous fillip to translate fundamental knowledge into applied research. Most importantly, VITALdb is an inevitable tool for selecting the best tool(s) to carry out a desired task and hence will prove to be a vital database (VITALdb) for the scientific community. Database URL: https://compbio.iitr.ac.in/vitaldb.

Viral influencers: deciphering the role of endogenous retroviral LTR12 repeats in cellular gene expression

The human genome is like a museum of ancient retroviral infections. It contains a large number of endogenous retroviruses (ERVs) that bear witness to past integration events. About 5,000 of them are so-called long terminal repeat 12 (LTR12) elements. Compared with 20,000 human genes, this is a remarkable number. Although LTR12 elements can act as promoters or enhancers of cellular genes, the function of most of these retroviral elements has remained unclear. In our mini-review, we show that different LTR12 elements share many similarities, including common transcription factor binding sites. Furthermore, we summarize novel insights into the epigenetic mechanisms governing their silencing and activation. Specific examples of genes and pathways that are regulated by LTR12 loci are used to illustrate the regulatory network built by these repetitive elements. A particular focus is on their role in the regulation of antiviral immune responses, tumor cell proliferation, and senescence. Finally, we describe how a targeted activation of this fascinating ERV family could be used for diagnostic or therapeutic purposes.

Endogenous viral elements: insights into data availability and accessibility

Endogenous viral elements (EVEs) are remnants of viral genetic material endogenized into the host genome. They have, in the last decades, attracted attention for their role as potential contributors to pathogenesis, drivers of selective advantage for the host, and genomic remnants of ancient viruses. EVEs have a nuanced and complex influence on both host health and evolution, and can offer insights on the deep evolutionary history of viruses. As an emerging field of research, several factors limit a comprehensive understanding of EVEs: they are currently underestimated and periodically overlooked in studies of the host genome, transcriptome, and virome. The absence of standardized guidelines for ensuring EVE-related data availability and accessibility following the FAIR ('findable, accessible, interoperable, and reusable') principles obstructs our ability to gather and connect information. Here, we discuss challenges to the availability and accessibility of EVE-related data and propose potential solutions. We identified the biological and research focus imbalance between different types of EVEs, and their overall biological complexity as genomic loci with viral ancestry, as potential challenges that can be addressed with the development of a user-oriented identification tool. In addition, reports of EVE identification are scattered between different subfields under different keywords, and EVE sequences and associated data are not properly gathered in databases. While developing an open and dedicated database might be ideal, targeted improvements of generalist databases might provide a pragmatic solution to EVE data and metadata accessibility. The implementation of these solutions, as well as the collective effort by the EVE scientific community in discussing and setting guidelines, is now drastically needed to lead the development of EVE research and offer insights into host-virus interactions and their evolutionary history.

Potential Contributions of Human Endogenous Retroviruses in Innate Immune Memory

The phenomenon wherein innate immune cells adopt long-term inflammatory phenotypes following the first stimuli is named trained immunity and can improve host defense against infections. Transcriptional and epigenetic reprogramming are critical mechanisms of trained immunity; however, the regulatory networks are not entirely clear at present. The human endogenous retroviruses (HERVs) provide large amounts of transcriptional regulators in the regulatory pathways. In this study, we analyzed published large omics data to explore the roles of such "dark matter" of the human genome in trained and tolerant macrophages. We collected 80 RNA sequencing data and 62 sequencing data to detect histone modifications and active regulatory regions from nine published studies on trained and tolerant macrophages. By analyzing the characteristics of transcription and epigenetic modification of HERVs, as well as their association with gene expression, we found that 15.3% of HERVs were transcribed nonrandomly from noncoding regions and enriched in specific HERV families and specific chromosomes, such as chromosomes 11, 15, 17, and 19, and they were highly related with the expression of adjacent genes. We found that 295 differentially expressed HERVs are located in 50-kbp flanking regions of 142 differentially expressed genes. We found epigenetic changes of these HERVs and that overlap with predicted enhancers and identified 35 enhancer-like HERVs. The related genes were highly involved in the activation and inflammatory responses, such as the TLR pathway. Other pathways including phosphoinositide signaling and transport of folate and K+ might be also related with trained immunity, which require further study. These results demonstrated that HERVs might play important roles in trained immunity.

Ocrelizumab B cell depletion has no effect on HERV RNA expression in PBMC in MS patients

BACKGROUND

Epstein barr virus (EBV) infection of B cells is now understood to be one of the triggering events for the development of Multiple Sclerosis (MS), a progressive immune-mediated disease of the central nervous system. EBV infection is also linked to expression of human endogenous retroviruses (HERVs) of the HERV-W group, a further risk factor for the development of MS. Ocrelizumab is a high-potency disease-modifying treatment (DMT) for MS, which depletes B cells by targeting CD20.

OBJECTIVES

We studied the effects of ocrelizumab on gene expression in peripheral blood mononuclear cells (PBMC) from paired samples from 20 patients taken prior to and 6 months after beginning ocrelizumab therapy. We hypothesised that EBV and HERV-W loads would be lower in post-treatment samples.

METHODS

Samples were collected in Paxgene tubes, subject to RNA extraction and Illumina paired end short read mRNA sequencing with mapping of sequence reads to the human genome using Salmon and differential gene expression compared with DeSeq2. Mapping was also performed separately to the HERV-D database of HERV sequences and the EBV reference sequence.

RESULTS

Patient samples were more strongly clustered by individual rather than disease type (relapsing/remitting or primary progressive), treatment (pre and post), age, or sex. Fourteen genes, all clearly linked to B cell function were significantly down regulated in the post treatment samples. Interestingly only one pre-treatment sample had detectable EBV RNA and there were no significant differences in HERV expression (of any group) between pre- and post-treatment samples.

CONCLUSIONS

While EBV and HERV expression are clearly linked to triggering MS pathogenesis, it does not appear that high level expression of these viruses is a part of the ongoing disease process or that changes in virus load are associated with ocrelizumab treatment.

Endogenous retroviral solo-LTRs in human genome

Human endogenous retroviruses (HERVs) are derived from the infection and integration of exogenetic retroviruses. HERVs account for 8% of human genome, and the majority of HERVs are solitary LTRs (solo-LTRs) due to homologous recombination. Multiple findings have showed that solo-LTRs could provide an enormous reservoir of transcriptional regulatory sequences involved in diverse biological processes, especially carcinogenesis and cancer development. The link between solo-LTRs and human diseases still remains poorly understood. This review focuses on the regulatory modules of solo-LTRs, which contribute greatly to the diversification and evolution of human genes. More importantly, although inactivating mutations, insertions and deletions have been identified in solo-LTRs, the inherited regulatory elements of solo-LTRs initiate the expression of chimeric lncRNA transcripts, which have been reported to play crucial roles in human health and disease. These findings provide valuable insights into the evolutionary and functional mechanisms underlying the presence of HERVs in human genome. Taken together, in this review, we will present evidences showing the regulatory and encoding capacity of solo-LTRs as well as the significant impact on various aspects of human biology.

HervD Atlas: a curated knowledgebase of associations between human endogenous retroviruses and diseases

Human endogenous retroviruses (HERVs), as remnants of ancient exogenous retrovirus infected and integrated into germ cells, comprise ∼8% of the human genome. These HERVs have been implicated in numerous diseases, and extensive research has been conducted to uncover their specific roles. Despite these efforts, a comprehensive source of HERV-disease association still needs to be added. To address this gap, we introduce the HervD Atlas (https://ngdc.cncb.ac.cn/hervd/), an integrated knowledgebase of HERV-disease associations manually curated from all related published literature. In the current version, HervD Atlas collects 60 726 HERV-disease associations from 254 publications (out of 4692 screened literature), covering 21 790 HERVs (21 049 HERV-Terms and 741 HERV-Elements) belonging to six types, 149 diseases and 610 related/affected genes. Notably, an interactive knowledge graph that systematically integrates all the HERV-disease associations and corresponding affected genes into a comprehensive network provides a powerful tool to uncover and deduce the complex interplay between HERVs and diseases. The HervD Atlas also features a user-friendly web interface that allows efficient browsing, searching, and downloading of all association information, research metadata, and annotation information. Overall, the HervD Atlas is an essential resource for comprehensive, up-to-date knowledge on HERV-disease research, potentially facilitating the development of novel HERV-associated diagnostic and therapeutic strategies.

VariantSurvival: a tool to identify genotype-treatment response

Motivation: For a number of neurological diseases, such as Alzheimer's disease, amyotrophic lateral sclerosis, and many others, certain genes are known to be involved in the disease mechanism. A common question is whether a structural variant in any such gene may be related to drug response in clinical trials and how this relationship can contribute to the lifecycle of drug development. Results: To this end, we introduce VariantSurvival, a tool that identifies changes in survival relative to structural variants within target genes. VariantSurvival matches annotated structural variants with genes that are clinically relevant to neurological diseases. A Cox regression model determines the change in survival between the placebo and clinical trial groups with respect to the number of structural variants in the drug target genes. We demonstrate the functionality of our approach with the exemplary case of the SETX gene. VariantSurvival has a user-friendly and lightweight graphical user interface built on the shiny web application package.

Navigating the Landscape: A Comprehensive Review of Current Virus Databases

Viruses are abundant and diverse entities that have important roles in public health, ecology, and agriculture. The identification and surveillance of viruses rely on an understanding of their genome organization, sequences, and replication strategy. Despite technological advancements in sequencing methods, our current understanding of virus diversity remains incomplete, highlighting the need to explore undiscovered viruses. Virus databases play a crucial role in providing access to sequences, annotations and other metadata, and analysis tools for studying viruses. However, there has not been a comprehensive review of virus databases in the last five years. This study aimed to fill this gap by identifying 24 active virus databases and included an extensive evaluation of their content, functionality and compliance with the FAIR principles. In this study, we thoroughly assessed the search capabilities of five database catalogs, which serve as comprehensive repositories housing a diverse array of databases and offering essential metadata. Moreover, we conducted a comprehensive review of different types of errors, encompassing taxonomy, names, missing information, sequences, sequence orientation, and chimeric sequences, with the intention of empowering users to effectively tackle these challenges. We expect this review to aid users in selecting suitable virus databases and other resources, and to help databases in error management and improve their adherence to the FAIR principles. The databases listed here represent the current knowledge of viruses and will help aid users find databases of interest based on content, functionality, and scope. The use of virus databases is integral to gaining new insights into the biology, evolution, and transmission of viruses, and developing new strategies to manage virus outbreaks and preserve global health.

Human endogenous retrovirus K contributes to a stem cell niche in glioblastoma

Human endogenous retroviruses (HERVs) are ancestral viral relics that constitute nearly 8% of the human genome. Although normally silenced, the most recently integrated provirus HERV-K (HML-2) can be reactivated in certain cancers. Here, we report pathological expression of HML-2 in malignant gliomas in both cerebrospinal fluid and tumor tissue that was associated with a cancer stem cell phenotype and poor outcomes. Using single-cell RNA-Seq, we identified glioblastoma cellular populations with elevated HML-2 transcripts in neural progenitor-like cells (NPC-like) that drive cellular plasticity. Using CRISPR interference, we demonstrate that HML-2 critically maintained glioblastoma stemness and tumorigenesis in both glioblastoma neurospheres and intracranial orthotopic murine models. Additionally, we demonstrate that HML-2 critically regulated embryonic stem cell programs in NPC-derived astroglia and altered their 3D cellular morphology by activating the nuclear transcription factor OCT4, which binds to an HML-2-specific long-terminal repeat (LTR5Hs). Moreover, we discovered that some glioblastoma cells formed immature retroviral virions, and inhibiting HML-2 expression with antiretroviral drugs reduced reverse transcriptase activity in the extracellular compartment, tumor viability, and pluripotency. Our results suggest that HML-2 fundamentally contributes to the glioblastoma stem cell niche. Because persistence of glioblastoma stem cells is considered responsible for treatment resistance and recurrence, HML-2 may serve as a unique therapeutic target.

A novel colony-stimulating factor 1 (CSF1) translocation involving human endogenous retroviral element in a tenosynovial giant cell tumor

Tenosynovial giant cell tumors (TSGCTs) are rare tumors arising in tendons or the synoviae of joints and bursae. The localized type is benign while the diffuse type shows expansive growth leading to greater morbidity and is therefore considered locally aggressive. Typical recurrent chromosomal aberrations are found in the majority of TSCGT and the CSF1 gene is frequently involved. In this article, we describe a newly identified gene fusion mediated by an inversion in a case of diffuse TSGCT. Multicolor-fluorescence in situ hybridization (FISH) molecular karyotyping identified a pericentric inversion of chromosome 1 in 7 out of 17 analyzed cells 46,XX,inv(1)(p13.3q24.3) [7]/46,XX [10], and with interphase FISH the involvement the CSF1 locus was detected. After performing transcriptome sequencing analysis for fusion detection, only one out of five fusion gene algorithms detected a fusion involving the CSF1 gene product. The resulting chimera fuses a sequence from a human endogenous retrovirus (HERV) gene to CSF1 Exon 6 on chromosome 1, abrogating the regulatory element of the 3' untranslated region of the CSF1 gene. This new translocation involving Exon 6 of the CSF1 gene fused to 1q24.1, supports the hypothesis that a mutated CSF1 protein is likely to play a vital role in the pathogenesis of TSGCT. The role of the HERV partner identified as a translocation partner, however, remains unclear. Our data add to the complexity of involved translocation partners in TSGCT and point to the potential difficulty of identifying fusion partners in tumor diagnostics using transcriptome sequencing when HERV or other repeat elements are involved.

Contribution of Retrotransposons to the Pathogenesis of Type 1 Diabetes and Challenges in Analysis Methods

Type 1 diabetes (T1D) is one of the most common chronic diseases of the endocrine system, associated with several life-threatening comorbidities. While the etiopathogenesis of T1D remains elusive, a combination of genetic susceptibility and environmental factors, such as microbial infections, are thought to be involved in the development of the disease. The prime model for studying the genetic component of T1D predisposition encompasses polymorphisms within the HLA (human leukocyte antigen) region responsible for the specificity of antigen presentation to lymphocytes. Apart from polymorphisms, genomic reorganization caused by repeat elements and endogenous viral elements (EVEs) might be involved in T1D predisposition. Such elements are human endogenous retroviruses (HERVs) and non-long terminal repeat (non-LTR) retrotransposons, including long and short interspersed nuclear elements (LINEs and SINEs). In line with their parasitic origin and selfish behaviour, retrotransposon-imposed gene regulation is a major source of genetic variation and instability in the human genome, and may represent the missing link between genetic susceptibility and environmental factors long thought to contribute to T1D onset. Autoreactive immune cell subtypes with differentially expressed retrotransposons can be identified with single-cell transcriptomics, and personalized assembled genomes can be constructed, which can then serve as a reference for predicting retrotransposon integration/restriction sites. Here we review what is known to date about retrotransposons, we discuss the involvement of viruses and retrotransposons in T1D predisposition, and finally we consider challenges in retrotransposons analysis methods.

Endogenous Retrovirus RNA Expression Differences between Race, Stage and HPV Status Offer Improved Prognostication among Women with Cervical Cancer

Endogenous human retroviruses (ERVs) are remnants of exogenous retroviruses that have integrated into the human genome. Using publicly available RNA-seq data from 63 cervical cancer patients, we investigated the expression of ERVs in cervical cancers. Four aspects of cervical cancer were investigated: patient ancestral background, tumor HPV type, tumor stage and patient survival. Between the racial subgroups, 74 ERVs were significantly differentially expressed, with Black Americans having 30 upregulated and 44 downregulated (including MER21C, HERV9-int, and HERVH-int) ERVs when compared to White Americans. We found that 3313 ERVs were differentially expressed between HPV subgroups, including MER41A, HERVH-int and HERVK9. There were 28 downregulated (including MLT1D and HERVH-int) and 61 upregulated (including MER41A) ERVs in locally advanced-stage compared to early-stage samples. Tissue microarrays of cervical cancer patients were used to investigate the protein expression of ERVs with protein coding potential (i.e., HERVK and ERV3). Significant differences in protein expression of ERV3 (p = 0.000905) were observed between early-stage and locally advanced-stage tumors. No significant differential expression at the protein level was found for HERVK7 (p = 0.243). We also investigated a prognostic model, supplementing a baseline prediction model using FIGO stage, age and HPV positivity with ERVs data. The expression levels of all ERVs in the HERVd were input into a Lasso-Cox proportional hazards model, developing a predictive 67-ERV panel. When ERVs expression levels were supplemented with the clinical data, a significant increase in prognostic power (p = 9.433 × 10-15) relative to that obtained with the clinical parameters alone (p = 0.06027) was observed. In summary, ERV RNA expression in cervical cancer tumors is significantly different among racial cohorts, HPV subgroups and disease stages. The combination of the expression of certain ERVs in cervical cancers with clinical factors significantly improved prognostication compared to clinical factors alone; therefore, ERVs may serve as future prognostic biomarkers and therapeutic targets. Novelty and Impact: When endogenous retroviral (ERV) expression signatures were combined with currently employed clinical prognosticators of relapse of cervical cancer, the combination outperformed prediction models based on clinical prognosticators alone. ERV expression signatures in tumor biopsies may therefore be useful to help identify patients at greater risk of recurrence. The novel ERV expression signatures or adjacent genes possibly impacted by ERV expression described here may also be targets for the development of future therapeutic interventions.

The landscape of hervRNAs transcribed from human endogenous retroviruses across human body sites

BACKGROUND

Human endogenous retroviruses (HERVs), the remnants of ancient retroviruses, account for 8% of the human genome, but most have lost their transcriptional abilities under physiological conditions. However, mounting evidence shows that several expressed HERVs do exert biological functions. Here, we systematically characterize physiologically expressed HERVs and examine whether they may give insight into the molecular fundamentals of human development and disease.

RESULTS

We systematically identify 13,889 expressed HERVs across normal body sites and demonstrate that they are expressed in body site-specific patterns and also by sex, ethnicity, and age. Analyzing cis-ERV-related quantitative trait loci, we find that 5435 hervRNAs are regulated by genetic variants. Combining this with a genome-wide association study, we elucidate that the dysregulation of expressed HERVs might be associated with various complex diseases, particularly neurodegenerative and psychiatric diseases. We further find that physiologically activated hervRNAs are associated with histone modifications rather than DNA demethylation.

CONCLUSIONS

Our results present a locus-specific landscape of physiologically expressed hervRNAs, which represent a hidden layer of genetic architecture in development and disease.

Influenza A Virus Infection Reactivates Human Endogenous Retroviruses Associated with Modulation of Antiviral Immunity

Human endogenous retrovirus (HERVs), normally silenced by methylation or mutations, can be reactivated by multiple environmental factors, including infections with exogenous viruses. In this work, we investigated the transcriptional activity of HERVs in human A549 cells infected by two wild-type (PR8M, SC35M) and one mutated (SC35MΔNS1) strains of Influenza A virus (IAVs). We found that the majority of differentially expressed HERVs (DEHERVS) and genes (DEGs) were up-regulated in the infected cells, with the most significantly enriched biological processes associated with the genes differentially expressed exclusively in SC35MΔNS1 being linked to the immune system. Most DEHERVs in PR8M and SC35M are mammalian apparent LTR retrotransposons, while in SC35MΔNS1, more HERV loci from the HERVW9 group were differentially expressed. Furthermore, up-regulated pairs of HERVs and genes in close chromosomal proximity to each other tended to be associated with immune responses, which implies that specific HERV groups might have the potential to trigger specific gene networks and influence host immunological pathways.

Switching Sides: How Endogenous Retroviruses Protect Us from Viral Infections

Long disregarded as junk DNA or genomic dark matter, endogenous retroviruses (ERVs) have turned out to represent important components of the antiviral immune response. These remnants of once-infectious retroviruses not only regulate cellular immune activation, but may even directly target invading viral pathogens. In this Gem, we summarize mechanisms by which retroviral fossils protect us from viral infections. One focus will be on recent advances in the role of ERVs as regulators of antiviral gene expression.

Transcriptional analysis of multiple ovarian cancer cohorts reveals prognostic and immunomodulatory consequences of ERV expression

BACKGROUND

Endogenous retroviruses (ERVs) play a role in a variety of biological processes, including embryogenesis and cancer. DNA methyltransferase inhibitors (DNMTi)-induced ERV expression triggers interferon responses in ovarian cancer cells via the viral sensing machinery. Baseline expression of ERVs also occurs in cancer cells, though this process is poorly understood and previously unexplored in epithelial ovarian cancer (EOC). Here, the prognostic and immunomodulatory consequences of baseline ERV expression was assessed in EOC.

METHODS

ERV expression was assessed using EOC transcriptional data from The Cancer Genome Atlas (TCGA) and from an independent cohort (Hammersmith Hospital, HH), as well as from untreated or DNMTi-treated EOC cell lines. Least absolute shrinkage and selection operator (LASSO) logistic regression defined an ERV expression score to predict patient prognosis. Immunohistochemistry (IHC) was conducted on the HH cohort. Combination of DNMTi treatment with γδ T cells was tested in vitro, using EOC cell lines and patient-derived tumor cells.

RESULTS

ERV expression was found to define clinically relevant subsets of EOC patients. An ERV prognostic score was successfully generated in TCGA and validated in the independent cohort. In EOC patients from this cohort, a high ERV score was associated with better survival (log-rank p=0.0009) and correlated with infiltration of CD8+PD1+T cells (r=0.46, p=0.0001). In the TCGA dataset, a higher ERV score was found in BRCA1/2 mutant tumors, compared to wild type (p=0.015), while a lower ERV score was found in CCNE1 amplified tumors, compared to wild type (p=0.019). In vitro, baseline ERV expression dictates the level of ERV induction in response to DNMTi. Manipulation of an ERV expression threshold by DNMTi resulted in improved EOC cell killing by cytotoxic immune cells.

CONCLUSIONS

These findings uncover the potential for baseline ERV expression to robustly inform EOC patient prognosis, influence tumor immune infiltration and affect antitumor immunity.

HIV-1 infection activates endogenous retroviral promoters regulating antiviral gene expression

Although endogenous retroviruses (ERVs) are known to harbor cis-regulatory elements, their role in modulating cellular immune responses remains poorly understood. Using an RNA-seq approach, we show that several members of the ERV9 lineage, particularly LTR12C elements, are activated upon HIV-1 infection of primary CD4⁺ T cells. Intriguingly, HIV-1-induced ERVs harboring transcription start sites are primarily found in the vicinity of immunity genes. For example, HIV-1 infection activates LTR12C elements upstream of the interferon-inducible genes GBP2 and GBP5 that encode for broad-spectrum antiviral factors. Reporter assays demonstrated that these LTR12C elements drive gene expression in primary CD4⁺ T cells. In line with this, HIV-1 infection triggered the expression of a unique GBP2 transcript variant by activating a cryptic transcription start site within LTR12C. Furthermore, stimulation with HIV-1-induced cytokines increased GBP2 and GBP5 expression in human cells, but not in macaque cells that naturally lack the GBP5 gene and the LTR12C element upstream of GBP2. Finally, our findings suggest that GBP2 and GBP5 have already been active against ancient viral pathogens as they suppress the maturation of the extinct retrovirus HERV-K (HML-2). In summary, our findings uncover how human cells can exploit remnants of once-infectious retroviruses to regulate antiviral gene expression.

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

Background

Endogenous retroviruses (ERVs) are remnants of ancient retroviral infections of mammalian germline cells. A large proportion of ERVs lose their open reading frames (ORFs), while others retain them and become exapted by the host species. However, it remains unclear what proportion of ERVs possess ORFs (ERV-ORFs), become transcribed, and serve as candidates for co-opted genes.

Results

We investigated characteristics of 176,401 ERV-ORFs containing retroviral-like protein domains (gag, pro, pol, and env) in 19 mammalian genomes. The fractions of ERVs possessing ORFs were overall small (~ 0.15%) although they varied depending on domain types as well as species. The observed divergence of ERV-ORF from their consensus sequences showed bimodal distributions, suggesting that a large proportion of ERV-ORFs either recently, or anciently, inserted themselves into mammalian genomes. Alternatively, very few ERVs lacking ORFs were found to exhibit similar divergence patterns. To identify candidates for ERV-derived genes, we estimated the ratio of non-synonymous to synonymous substitution rates (dN/dS) for ERV-ORFs in human and non-human mammalian pairs, and found that approximately 42% of the ERV-ORFs showed dN/dS < 1. Further, using functional genomics data including transcriptome sequencing, we determined that approximately 9.7% of these selected ERV-ORFs exhibited transcriptional potential.

Conclusions

These results suggest that purifying selection operates on a certain portion of ERV-ORFs, some of which may correspond to uncharacterized functional genes hidden within mammalian genomes. Together, our analyses suggest that more ERV-ORFs may be co-opted in a host-species specific manner than we currently know, which are likely to have contributed to mammalian evolution and diversification.

Transcription profile of human endogenous retroviruses in response to dengue virus serotype 2 infection

Human endogenous retroviruses (HERVs), the remains of retroviruses infection in our ancestors' germline cell over millions of years, take up about 8% of the human genome in total. HERV transcription has been detected in various cancers and diseases. However, the interaction between HERV expression and viral infection has not been fully elucidated. Here, we provided the first transcriptional profile of HERVs in dengue virus serotype 2 (DENV-2) infected A549 cells by using high-throughput RNA sequencing. The results showed that a number of HERVs and human genes were significantly differentially expressed in response to DENV-2 infection. Further bioinformatic analyses indicated a correlation between HERVs and human genes. In particular, the genes near the HERVs activated by dengue infection were dominantly enriched in the antiviral immune response. Taken together, our findings suggest that activated HERVs may be involved in the cellular immune response to viral infection by coexpressing with nearby host genes.

Endogenous Retroviruses Activity as a Molecular Signature of Neurodevelopmental Disorders

Human endogenous retroviruses (HERVs) are genetic elements resulting from relics of ancestral infection of germline cells, now recognized as cofactors in the etiology of several complex diseases. Here we present a review of findings supporting the role of the abnormal HERVs activity in neurodevelopmental disorders. The derailment of brain development underlies numerous neuropsychiatric conditions, likely starting during prenatal life and carrying on during subsequent maturation of the brain. Autism spectrum disorders, attention deficit hyperactivity disorders, and schizophrenia are neurodevelopmental disorders that arise clinically during early childhood or adolescence, currently attributed to the interplay among genetic vulnerability, environmental risk factors, and maternal immune activation. The role of HERVs in human embryogenesis, their intrinsic responsiveness to external stimuli, and the interaction with the immune system support the involvement of HERVs in the derailed neurodevelopmental process. Although definitive proofs that HERVs are involved in neurobehavioral alterations are still lacking, both preclinical models and human studies indicate that the abnormal expression of ERVs could represent a neurodevelopmental disorders-associated biological trait in affected individuals and their parents.

ERVmap analysis reveals genome-wide transcription of human endogenous retroviruses

Endogenous retrovirus (ERV) sequences make up a large fraction of our genome, yet little is understood about their function and biological relevance. Deep-sequencing data contain valuable information on a genome-wide scale. Yet, due to their highly repetitive nature, analysis of ERVs has been computationally challenging. We describe a bioinformatics tool called ERVmap to analyze transcription of unique sets of human ERVs in a range of cell types in health and disease settings. Our open-source code and accompanied web tool should facilitate researchers in all fields to study the expression patterns of ERVs in sequencing data and should lead to significant advancement in understanding the biological relevance of ERVs in health and disease.

Endogenous retroviruses (ERVs) are integrated retroviral elements that make up 8% of the human genome. However, the impact of ERVs on human health and disease is not well understood. While select ERVs have been implicated in diseases, including autoimmune disease and cancer, the lack of tools to analyze genome-wide, locus-specific expression of proviral autonomous ERVs has hampered the progress in the field. Here we describe a method called ERVmap, consisting of an annotated database of 3,220 human proviral ERVs and a pipeline that allows for locus-specific genome-wide identification of proviral ERVs that are transcribed based on RNA-sequencing data, and provide examples of the utility of this tool. Using ERVmap, we revealed cell-type–specific ERV expression patterns in commonly used cell lines as well as in primary cells. We identified 124 unique ERV loci that are significantly elevated in the peripheral blood mononuclear cells of patients with systemic lupus erythematosus that represent an IFN-independent signature. Finally, we identified additional tumor-associated ERVs that correlate with cytolytic activity represented by granzyme and perforin expression in breast cancer tissue samples. The open-source code of ERVmap and the accompanied web tool are made publicly available to quantify proviral ERVs in RNA-sequencing data with ease. Use of ERVmap across a range of diseases and experimental conditions has the potential to uncover novel disease-associated antigens and effectors involved in human health that is currently missed by focusing on protein-coding sequences.

Endogenous retroviral signatures predict immunotherapy response in clear cell renal cell carcinoma

Human endogenous retroviruses (hERVs) are remnants of exogenous retroviruses that have integrated into the genome throughout evolution. We developed a computational workflow, hervQuant, which identified more than 3,000 transcriptionally active hERVs within The Cancer Genome Atlas (TCGA) pan-cancer RNA-Seq database. hERV expression was associated with clinical prognosis in several tumor types, most significantly clear cell renal cell carcinoma (ccRCC). We explored two mechanisms by which hERV expression may influence the tumor immune microenvironment in ccRCC: (i) RIG-I-like signaling and (ii) retroviral antigen activation of adaptive immunity. We demonstrated the ability of hERV signatures associated with these immune mechanisms to predict patient survival in ccRCC, independent of clinical staging and molecular subtyping. We identified potential tumor-specific hERV epitopes with evidence of translational activity through the use of a ccRCC ribosome profiling (Ribo-Seq) dataset, validated their ability to bind HLA in vitro, and identified the presence of MHC tetramer-positive T cells against predicted epitopes. hERV sequences identified through this screening approach were significantly more highly expressed in ccRCC tumors responsive to treatment with programmed death receptor 1 (PD-1) inhibition. hervQuant provides insights into the role of hERVs within the tumor immune microenvironment, as well as evidence that hERV expression could serve as a biomarker for patient prognosis and response to immunotherapy.

Transcriptional Modulation of Human Endogenous Retroviruses in Primary CD4+ T Cells Following Vorinostat Treatment

The greatest obstacle to a cure for HIV is the provirus that integrates into the genome of the infected cell and persists despite antiretroviral therapy. A "shock and kill" approach has been proposed as a strategy for an HIV cure whereby drugs and compounds referred to as latency-reversing agents (LRAs) are used to "shock" the silent provirus into active replication to permit "killing" by virus-induced pathology or immune recognition. The LRA most utilized to date in clinical trials has been the histone deacetylase (HDAC) inhibitor-vorinostat. Potentially, pathological off-target effects of vorinostat may result from the activation of human endogenous retroviruses (HERVs), which share common ancestry with exogenous retroviruses including HIV. To explore the effects of HDAC inhibition on HERV transcription, an unbiased pharmacogenomics approach (total RNA-Seq) was used to evaluate HERV expression following the exposure of primary CD4+ T cells to a high dose of vorinostat. Over 2,000 individual HERV elements were found to be significantly modulated by vorinostat, whereby elements belonging to the ERVL family (e.g., LTR16C and LTR33) were predominantly downregulated, in contrast to LTR12 elements of the HERV-9 family, which exhibited the greatest signal, with the upregulation of 140 distinct elements. The modulation of three different LTR12 elements by vorinostat was confirmed by droplet digital PCR along a dose-response curve. The monitoring of LTR12 expression during clinical trials with vorinostat may be indicated to assess the impact of this HERV on the human genome and host immunity.

Analysis of Human Endogenous Retrovirus Expression in Multiple Sclerosis Plaques

BACKGROUND

It has been suggested that Human endogenous retroviruses (HERVs) are associated with multiple sclerosis (MS) pathogenesis. The objective of this study was to broadly evaluate the expression of HERV core (GAG) and envelope (ENV) genes in diseased brain white matter samples from MS patients compared to normal controls.

METHODS

Twenty-eight HERV GAG and 88 ENV gene sequences were retrieved, classified by phylogeny, and grouped into clades. Consensus qPCR primers were designed for each clade, and quantitative PCR was performed on 33 MS and 9 normal control frozen brain samples. MS samples included chronic progressive (n=5), primary progressive (n=4), secondary progressive (n=14), relapsing remitting (n=3) and unclassified confirmed MS cases (n=7). The levels of GAG and ENV RNA within each of the samples were quantitated and normalized using the neuronal reference gene RPL19. Expression differences were analyzed for MS vs control.

RESULTS

Expression of GAG clades 1A, 3B, and 3C mapping to HERV-E and HERV-K were significantly increased compared to controls, while GAG clade 3A expression was decreased. Expression of HERV ENV clades 2, 3A, 3B, mapping to RTVL, HERV-E and HERV-K and MSRV (HERV-W), were significantly increased in the MS group. However, the relative expression differences between the MS and control groups were small, differing less than 1.5-fold.

CONCLUSION

Expression of GAG and ENV mapping to HERV-E, RTVL and HERV-K10 families were significantly increased in the MS group. However, the relative expression differences between the MS and control groups were small, differing less than 1.5-fold. These results indicate that the expression of HERV GAG and ENV regions do not differ greatly between MS and controls in these frozen brain samples.

How to tame an endogenous retrovirus: HERVH and the evolution of human pluripotency

HERVH is one of the most successful endogenous retrovirus in the human genome. Relative to other endogenous retroviruses, slower degradation of HERVH internal sequences indicates their potential relevance for the host. HERVH is transcriptionally active during human preimplantation embryogenesis. In this review, we focus on the role of HERVH in regulating human pluripotency. The HERVH-mediated pluripotency network has been evolved recently in primates. Nevertheless, it became an essential feature of human pluripotency. We discuss how HERVH modulates the human pluripotency network by providing alternative transcription factor binding sites, functioning as a long-range enhancer, and as being a major source for pluripotency specific long non-coding RNAs and chimeric transcripts.

EnHERV: Enrichment analysis of specific human endogenous retrovirus patterns and their neighboring genes

Human endogenous retroviruses (HERVs) are flanked by long terminal repeats (LTRs), which contain the regulation part of the retrovirus. Remaining HERVs constitute 7% to 8% of the present day human genome, and most have been identified as solo LTRs. The HERV sequences have been associated with several molecular functions as well as certain diseases in human, but their roles in human diseases are yet to be established. We designed EnHERV to make accessible the identified endogenous retrovirus repetitive sequences from Repbase Update (a database of eukaryotic repetitive elements) that are present in the human genome. Defragmentation process was done to improve the RepeatMasker annotation output. The defragmented elements were used as core database in EnHERV. EnHERV is available at http://sysbio.chula.ac.th/enherv and can be searched using either gene lists of user interest or HERV characteristics. Besides the search function, EnHERV also provides an enrichment analysis function that allows users to perform enrichment analysis between selected HERV characteristics and user-input gene lists, especially genes with the expression profile of a certain disease. EnHERV will facilitate exploratory studies of specific HERV characteristics that control gene expression patterns related to various disease conditions. Here we analyzed 25 selected HERV groups/names from all four HERV superfamilies, using the sense and anti-sense directions of the HERV and gene expression profiles from 49 specific tissue and disease conditions. We found that intragenic HERVs were associated with down-regulated genes in most cancer conditions and in psoriatic skin tissues and associated with up-regulated genes in immune cells particularly from systemic lupus erythematosus (SLE) patients. EnHERV allowed the analysis of how different types of LTRs were differentially associated with specific gene expression profiles in particular disease conditions for further studies into their mechanisms and functions.

Phylogenetic Analysis Reveals That ERVs "Die Young" but HERV-H Is Unusually Conserved

About 8% of the human genome is made up of endogenous retroviruses (ERVs). Though most human endogenous retroviruses (HERVs) are thought to be irrelevant to our biology notable exceptions include members of the HERV-H family that are necessary for the correct functioning of stem cells. ERVs are commonly found in two forms, the full-length proviral form, and the more numerous solo-LTR form, thought to result from homologous recombination events. Here we introduce a phylogenetic framework to study ERV insertion and solo-LTR formation. We then apply the framework to site patterns sampled from a set of long alignments covering six primate genomes. Studying six categories of ERVs we quantitatively recapitulate patterns of insertional activity that are usually described in qualitative terms in the literature. A slowdown in most ERV groups is observed but we suggest that HERV-K activity may have increased in humans since they diverged from chimpanzees. We find that the rate of solo-LTR formation decreases rapidly as a function of ERV age and that an age dependent model of solo-LTR formation describes the history of ERVs more accurately than the commonly used exponential decay model. We also demonstrate that HERV-H loci are markedly less likely to form solo-LTRs than ERVs from other families. We conclude that the slower dynamics of HERV-H suggest a host role for the internal regions of these exapted elements and posit that in future it will be possible to use the relationship between full-length proviruses and solo-LTRs to help identify large scale co-options in distant vertebrate genomes.

Animal genomes contain ancient pathogens known as endogenous retroviruses (ERVs). Though the widespread abundance of ERVs is due to their ability to self replicate, some ERVs are known to have become important to host processes including placentation, and in the case of HERV-H, the functioning of human stem cells. In our study we place the insertion and deletion activity of primate ERV families in direct quantitative comparison. In particular, we show that ERV deletion is an age dependent process, so that as an ERV ages it becomes less likely to be deleted at any given instant. We also find that ERVs from the HERV-H family are unusually slowly deleted, an interesting result that suggests that the exaptation of HERV-H may have involved internal regions of the virus and not just its terminal promoters. Assuming the behaviour of primate ERVs is not unusual, our study suggests that future bioinformatics screening for ERVs with slow deletion dynamics could help identify large-scale exaptations in distant species. Furthermore, as we demonstrate that ERVs are deleted rapidly, we think that such screening could be performed using ratios of conserved to deleted elements and could therefore be applied to single genomes.

A Novel, Highly Selective RT-QPCR Method for Quantification of MSRV Using PNA Clamping Syncytin-1 (ERVWE1)

HERV-W is a multi-locus family of human endogenous retroviruses (HERVs) that has been found to play an important role in human physiology and pathology. Two particular members of HERV-W family are of special interests: ERVWE1 (coding syncytin-1, which is a glycoprotein essential in the formation of the placenta) and MSRV (multiple sclerosis-associated retrovirus that is thought to play a significant role in human pathology as a result of its increased expression in the brain tissue and blood cells derived from patients with multiple sclerosis (MS)). Both ERVWE1 and MSRV mRNA share high level of similarity and hence a method that allows to exclusively quantify the MSRV expression in clinical samples would be desirable. We developed a quantitative polymerase chain reaction (QPCR) technique for the detection and quantification of the multiple sclerosis-associated retrovirus. The assay utilises fluorescently labelled oligonucleotide probe, which is complementary to the conservative fragment of MSRV env gene and a peptide nucleic acid (PNA) probe, fully complementary to the ERVWE1 sequence fragment that efficiently blocks the polymerase action on ERVWE1 templates. The PNA molecule, if used parallel with hydrolysis probe in QPCR analysis, greatly facilitates the detection efficiency of MSRV even if ERVWE1 is present abundantly in respect to MSRV in the analysed sample. We achieved a wide and measurable range from 1 × 10 e² to 1 × 10 e⁸ copies/reaction; the linearity of the technique was maintained even at the low MSRV level of 1 % in respect to ERVWE1. Using our newly developed method we confirmed that the expression of MSRV takes place in normal human astrocytes and in human umbilical vein endothelial cells in vitro. We also found that the stimulation of human monocytes did not influence the specific expression of MSRV but it caused changes in mRNA level of distinct HERV-W templates.

Differential expression analysis of human endogenous retroviruses based on ENCODE RNA-seq data

Background

Human endogenous retroviruses (HERVs) are flanked by long terminal repeats (LTRs), which possess promoter activity and can therefore influence the expression of neighboring genes. HERV involvement in different types of cancer has already been thoroughly documented. However, so far there has been no systematic study of HERV expression patterns in a multitude of cell types in health and disease. In particular, the publication of the comprehensive ENCODE dataset has already facilitated many gene expression studies, but none so far focusing exclusively on HERVs.

Results

We present a comprehensive differential analysis of HERV expression based on ENCODE Tier 1 and Tier 2 RNA-seq data produced by Cold Spring Harbor Laboratories and the California Institute of Technology. This analysis was conducted for individual HERV loci and for entire HERV families in twelve different cell lines, of which six correspond to the normal condition and the other six represent cancer cell types. Although the principal component analysis revealed that the two groups of cells show distinguishable expression patterns, we were not able to link these differences to one or multiple particular HERV families. Two samples exhibit expression patterns, which are not similar to the corresponding cell lines of the other producing lab. Instead they show signs of cancer formation and expression of the pluripotency marker HERVH, despite being classified as a normal cell line and a differentiated cell, respectively.

Conclusions

Our study demonstrates that ENCODE data are generally comparable between the different contributing labs and that the analysis of HERV elements can provide novel insights into differentiation and disease state of a cell that are easily overlooked when focusing on protein-coding genes. Our findings hint at a change in HERV expression during cancerogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0146-5) contains supplementary material, which is available to authorized users.

Endogenous viruses: Connecting recent and ancient viral evolution

The rapid rates of viral evolution allow us to reconstruct the recent history of viruses in great detail. This feature, however, also results in rapid erosion of evolutionary signal within viral molecular data, impeding studies of their deep history. Thus, the further back in time, the less accurate the inference becomes. Furthermore, reconstructing complex histories of transmission can be challenging, especially where extinct viral lineages are concerned. This problem has been partially solved by the discovery of viruses embedded in host genomes, known as endogenous viral elements (EVEs). Some of these endogenous viruses are derived from ancient relatives of extant viruses, allowing us to better examine ancient viral host range, geographical distribution and transmission routes. Moreover, our knowledge of viral evolutionary timescales and rate dynamics has also been greatly improved by their discovery, thereby bridging the gap between recent and ancient viral evolution.

Unraveling the web of viroinformatics: computational tools and databases in virus research

The beginning of the second century of research in the field of virology (the first virus was discovered in 1898) was marked by its amalgamation with bioinformatics, resulting in the birth of a new domain--viroinformatics. The availability of more than 100 Web servers and databases embracing all or specific viruses (for example, dengue virus, influenza virus, hepatitis virus, human immunodeficiency virus [HIV], hemorrhagic fever virus [HFV], human papillomavirus [HPV], West Nile virus, etc.) as well as distinct applications (comparative/diversity analysis, viral recombination, small interfering RNA [siRNA]/short hairpin RNA [shRNA]/microRNA [miRNA] studies, RNA folding, protein-protein interaction, structural analysis, and phylotyping and genotyping) will definitely aid the development of effective drugs and vaccines. However, information about their access and utility is not available at any single source or on any single platform. Therefore, a compendium of various computational tools and resources dedicated specifically to virology is presented in this article.

Human endogenous retrovirus-K(II) envelope induction protects neurons during HIV/AIDS

Human endogenous retroviruses (HERVs) are differentially expressed depending on the cell type and physiological circumstances. HERV-K has been implicated in the pathogenesis of several diseases although the functional consequences of its expression remain unknown. Human immunodeficiency virus (HIV) infection causes neuroinflammation with neuronal damage and death. Herein, we investigated HERV-K(II)/(HML-2) envelope (Env) expression and its actions in the brain during HIV/AIDS. HERV-K(II) Env expression was assessed in healthy brain tissues, autopsied HIV HIV− infected (HIV+) and uninfected (HIV−) brains and in neural cell cultures by real time RT-PCR, massively parallel (deep) sequencing, immunoblotting and immunohistochemistry. Neuronal and neural stem cells expressing HERV-K(II) Env were analyzed in assays of host responses including cellular viability, immune responses and neurobehavioral outcomes. Deep sequencing of human brain transcriptomes disclosed that RNA sequences encoded by HERV-K were among the most abundant HERV sequences detected in human brain. Comparison of different cell types revealed that HERV-K(II) env RNA abundance was highest in cultured human neurons but was suppressed by epidermal growth factor exposure. HERV-K(II) Env immunoreactivity was increased in the cerebral cortex from persons with HIV/AIDS, principally localized in neurons. Human neuronal cells transfected with HERV-K(II) Env exhibited increased NGF and BDNF expression. Expression of HERV-K(II) Env in neuronal cells increased cellular viability and prevented neurotoxicity mediated by HIV-1 Vpr. Intracerebral delivery of HERV-K(II) Env expressed by neural stem cells suppressed TNF-α expression and microglial activation while also improving neurobehavioral deficits in vpr/RAG1^−/− mice. HERV-K(II) Env was highly expressed in human neurons, especially during HIV/AIDS, but in addition exerted neuroprotective effects. These findings imply that HERV gene products might exert adaptive effects in circumstances of pathophysiological stress, perhaps underlying the conservation of HERVs within the human genome.

Porcine endogenous retrovirus (PERV) infection of HEK-293 cell line alters expression of human endogenous retrovirus (HERV-W) sequences

The risk of infections of human recipients after xenotransplantations is now mainly represented by porcine endogenous retroviruses (PERVs) as these particles are part of the porcine genome. As in all vertebrates, human genome harbours its own numerous genetic sequences of retroviral origin; it is estimated that they comprise about 8 % of the human genome. Because some of them play an important role in human physiology, it is valuable to estimate whether the presence of PERVs in human cells influences homeostasis of the human endogenous retrovirus (HERV) expression pattern. The aim of the study was to evaluate whether the expression profile of HERV-W genes changes after infection of cells by porcine endogenous retroviruses. In the experimental settings, human embryonic kidney cell line (HEK-293) was infected by PERV particles and cultivated up to 22th passage after infection. HERV-W gag, pol and env, as well as env from locus 7q21.2 gene expression was monitored by means of realtime reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and Western blot techniques. We found that the expression level of HERV-W genes differs in PERV-infected HEK-293 cell cultures in comparison with that from non-infected cultures. Relative HERV-W gene expression also differed significantly between particular passages (P < 0.05). Moreover, we have noticed a high correlation between the HERV-W Env(7q21.2) mRNA and protein level (Spearman rank r = 0.65; P < 0.05) during the course of the experiment. As previously hypothesized, human genomic sequences of retroviral origin may be changed by the presence of porcine endogenous retroviruses.

Association of endogenous retroviruses and long terminal repeats with human disorders

Since the human genome sequences became available in 2001, our knowledge about the human transposable elements which comprise ∼40% of the total nucleotides has been expanding. Non-long terminal repeat (non-LTR) retrotransposons are actively transposing in the present-day human genome, and have been found to cause ∼100 identified clinical cases of varied disorders. In contrast, almost all of the human endogenous retroviruses (HERVs) originating from ancient infectious retroviruses lost their infectivity and transposing activity at various times before the human-chimpanzee speciation (∼6 million years ago), and no known HERV is presently infectious. Insertion of HERVs and mammalian apparent LTR retrotransposons (MaLRs) into the chromosomal DNA influenced a number of host genes in various modes during human evolution. Apart from the aspect of genome evolution, HERVs and solitary LTRs being suppressed in normal biological processes can potentially act as extra transcriptional apparatuses of cellular genes by re-activation in individuals. There has been a reasonable prediction that aberrant LTR activation could trigger malignant disorders and autoimmune responses if epigenetic changes including DNA hypomethylation occur in somatic cells. Evidence supporting this hypothesis has begun to emerge only recently: a MaLR family LTR activation in the pathogenesis of Hodgkin’s lymphoma and a HERV-E antigen expression in an anti-renal cell carcinoma immune response. This mini review addresses the impacts of the remnant-form LTR retrotransposons on human pathogenesis.

Microarray-based sketches of the HERV transcriptome landscape

Human endogenous retroviruses (HERVs) are spread throughout the genome and their long terminal repeats (LTRs) constitute a wide collection of putative regulatory sequences. Phylogenetic similarities and the profusion of integration sites, two inherent characteristics of transposable elements, make it difficult to study individual locus expression in a large-scale approach, and historically apart from some placental and testis-regulated elements, it was generally accepted that HERVs are silent due to epigenetic control. Herein, we have introduced a generic method aiming to optimally characterize individual loci associated with 25-mer probes by minimizing cross-hybridization risks. We therefore set up a microarray dedicated to a collection of 5,573 HERVs that can reasonably be assigned to a unique genomic position. We obtained a first view of the HERV transcriptome by using a composite panel of 40 normal and 39 tumor samples. The experiment showed that almost one third of the HERV repertoire is indeed transcribed. The HERV transcriptome follows tropism rules, is sensitive to the state of differentiation and, unexpectedly, seems not to correlate with the age of the HERV families. The probeset definition within the U3 and U5 regions was used to assign a function to some LTRs (i.e. promoter or polyA) and revealed that (i) autonomous active LTRs are broadly subjected to operational determinism (ii) the cellular gene density is substantially higher in the surrounding environment of active LTRs compared to silent LTRs and (iii) the configuration of neighboring cellular genes differs between active and silent LTRs, showing an approximately 8 kb zone upstream of promoter LTRs characterized by a drastic reduction in sense cellular genes. These gathered observations are discussed in terms of virus/host adaptive strategies, and together with the methods and tools developed for this purpose, this work paves the way for further HERV transcriptome projects.

A recurrent translocation is mediated by homologous recombination between HERV-H elements

BACKGROUND

Chromosome rearrangements are caused by many mutational mechanisms; of these, recurrent rearrangements can be particularly informative for teasing apart DNA sequence-specific factors. Some recurrent translocations are mediated by homologous recombination between large blocks of segmental duplications on different chromosomes. Here we describe a recurrent unbalanced translocation casued by recombination between shorter homologous regions on chromosomes 4 and 18 in two unrelated children with intellectual disability.

RESULTS

Array CGH resolved the breakpoints of the 6.97-Megabase (Mb) loss of 18q and the 7.30-Mb gain of 4q. Sequencing across the translocation breakpoints revealed that both translocations occurred between 92%-identical human endogenous retrovirus (HERV) elements in the same orientation on chromosomes 4 and 18. In addition, we find sequence variation in the chromosome 4 HERV that makes one allele more like the chromosome 18 HERV.

CONCLUSIONS

Homologous recombination between HERVs on the same chromosome is known to cause chromosome deletions, but this is the first report of interchromosomal HERV-HERV recombination leading to a translocation. It is possible that normal sequence variation in substrates of non-allelic homologous recombination (NAHR) affects the alignment of recombining segments and influences the propensity to chromosome rearrangement.

Bioinformatics and genomic analysis of transposable elements in eukaryotic genomes

A major portion of most eukaryotic genomes are transposable elements (TEs). During evolution, TEs have introduced profound changes to genome size, structure, and function. As integral parts of genomes, the dynamic presence of TEs will continue to be a major force in reshaping genomes. Early computational analyses of TEs in genome sequences focused on filtering out "junk" sequences to facilitate gene annotation. When the high abundance and diversity of TEs in eukaryotic genomes were recognized, these early efforts transformed into the systematic genome-wide categorization and classification of TEs. The availability of genomic sequence data reversed the classical genetic approaches to discovering new TE families and superfamilies. Curated TE databases and their accurate annotation of genome sequences in turn facilitated the studies on TEs in a number of frontiers including: (1) TE-mediated changes of genome size and structure, (2) the influence of TEs on genome and gene functions, (3) TE regulation by host, (4) the evolution of TEs and their population dynamics, and (5) genomic scale studies of TE activity. Bioinformatics and genomic approaches have become an integral part of large-scale studies on TEs to extract information with pure in silico analyses or to assist wet lab experimental studies. The current revolution in genome sequencing technology facilitates further progress in the existing frontiers of research and emergence of new initiatives. The rapid generation of large-sequence datasets at record low costs on a routine basis is challenging the computing industry on storage capacity and manipulation speed and the bioinformatics community for improvement in algorithms and their implementations.

Human endogenous retrovirus K106 (HERV-K106) was infectious after the emergence of anatomically modern humans

HERV-K113 and HERV-K115 have been considered to be among the youngest HERVs because they are the only known full-length proviruses that are insertionally polymorphic and maintain the open reading frames of their coding genes. However, recent data suggest that HERV-K113 is at least 800,000 years old, and HERV-K115 even older. A systematic study of HERV-K HML2 members to identify HERVs that may have infected the human genome in the more recent evolutionary past is lacking. Therefore, we sought to determine how recently HERVs were exogenous and infectious by examining sequence variation in the long terminal repeat (LTR) regions of all full-length HERV-K loci. We used the traditional method of inter-LTR comparison to analyze all full length HERV-Ks and determined that two insertions, HERV-K106 and HERV-K116 have no differences between their 5′ and 3′ LTR sequences, suggesting that these insertions were endogenized in the recent evolutionary past. Among these insertions with no sequence differences between their LTR regions, HERV-K106 had the most intact viral sequence structure. Coalescent analysis of HERV-K106 3′ LTR sequences representing 51 ethnically diverse individuals suggests that HERV-K106 integrated into the human germ line approximately 150,000 years ago, after the emergence of anatomically modern humans.

Effects of retroviruses on host genome function

For millions of years, retroviral infections have challenged vertebrates, occasionally leading to germline integration and inheritance as ERVs, genetic parasites whose remnants today constitute some 7% to 8% of the human genome. Although they have had significant evolutionary side effects, it is useful to view ERVs as fossil representatives of retroviruses extant at the time of their insertion into the germline and not as direct players in the evolutionary process itself. Expression of particular ERVs is associated with several positive physiological functions as well as certain diseases, although their roles in human disease as etiological agents, possible contributing factors, or disease markers-well demonstrated in animal models-remain to be established. Here we discuss ERV contributions to host genome structure and function, including their ability to mediate recombination, and physiological effects on the host transcriptome resulting from their integration, expression, and other events.

Human endogenous retrovirus K (HML-2) elements in the plasma of people with lymphoma and breast cancer

Actively replicating endogenous retroviruses entered the human genome millions of years ago and became a stable part of the inherited genetic material. They subsequently acquired multiple mutations, leading to the assumption that these viruses no longer replicate. However, certain human tumor cell lines have been shown to release endogenous retroviral particles. Here we show that RNA from human endogenous retrovirus K (HERV-K) (HML-2), a relatively recent entrant into the human genome, can be found in very high titers in the plasma of patients with lymphomas and breast cancer as measured by either reverse transcriptase PCR or nucleic acid sequence-based amplification. Further, these titers drop dramatically with cancer treatment. We also demonstrate the presence of reverse transcriptase and viral RNA in plasma fractions that contain both immature and correctly processed HERV-K (HML-2) Gag and envelope proteins. Finally, using immunoelectron microscopy, we show the presence of HERV-K (HML-2) virus-like particles in the plasma of lymphoma patients. Taken together, these findings demonstrate that elements of the endogenous retrovirus HERV-K (HML-2) can be found in the blood of modern-day humans with certain cancers.

A paradigm for virus-host coevolution: sequential counter-adaptations between endogenous and exogenous retroviruses

Endogenous retroviruses (ERVs) are remnants of ancient retroviral infections of the host germline transmitted vertically from generation to generation. It is hypothesized that some ERVs are used by the host as restriction factors to block the infection of pathogenic retroviruses. Indeed, some ERVs efficiently interfere with the replication of related exogenous retroviruses. However, data suggesting that these mechanisms have influenced the coevolution of endogenous and/or exogenous retroviruses and their hosts have been more difficult to obtain. Sheep are an interesting model system to study retrovirus-host coevolution because of the coexistence in this animal species of two exogenous (i.e., horizontally transmitted) oncogenic retroviruses, Jaagsiekte sheep retrovirus and Enzootic nasal tumor virus, with highly related and biologically active endogenous retroviruses (enJSRVs). Here, we isolated and characterized the evolutionary history and molecular virology of 27 enJSRV proviruses. enJSRVs have been integrating in the host genome for the last 5–7 million y. Two enJSRV proviruses (enJS56A1 and enJSRV-20), which entered the host genome within the last 3 million y (before and during speciation within the genus Ovis), acquired in two temporally distinct events a defective Gag polyprotein resulting in a transdominant phenotype able to block late replication steps of related exogenous retroviruses. Both transdominant proviruses became fixed in the host genome before or around sheep domestication (∼ 9,000 y ago). Interestingly, a provirus escaping the transdominant enJSRVs has emerged very recently, most likely within the last 200 y. Thus, we determined sequentially distinct events during evolution that are indicative of an evolutionary antagonism between endogenous and exogenous retroviruses. This study strongly suggests that endogenization and selection of ERVs acting as restriction factors is a mechanism used by the host to fight retroviral infections.

The genome of all vertebrates is heavily colonized by “endogenous” retroviruses (ERVs). ERVs derive from retrovirus infections of the germ cells of the host during evolution, leading to permanent integration of the viral genome into the host DNA. Because ERVs are integrated in the host genome, they are transmitted to subsequent generations like any other host gene. The function of endogenous retroviruses is not completely clear, but some ERVs can block the replication cycle of horizontally transmitted “exogenous” pathogenic retroviruses. These observations lead to the hypothesis that ERVs have protected the host during evolution against incoming pathogenic retroviruses. Here, by characterizing the evolutionary history and molecular virology of a particular group of endogenous betaretroviruses of sheep (enJSRVs) we show a fascinating series of events unveiling the endless struggle between host and retroviruses. In particular, we discovered that: (i) two enJSRV loci that entered the host genome before speciation within the genus Ovis (∼ 3 million y ago) acquired, after their integration, a mutated defective viral protein capable of blocking exogenous related retroviruses; (ii) both these transdominant enJSRV loci became fixed in the host genome before or around sheep domestication (∼ 10,000 y ago); (iii) the invasion of the sheep genome by ERVs of the JSRV/enJSRVs group is still in progress; and (iv) new viruses have recently emerged (less than 200 y ago) that can escape the transdominant enJSRV loci. This study strongly suggests that endogenization and selection of ERVs acting as restriction factors is a mechanism used by the host to fight retroviral infections.

Rate of recombinational deletion among human endogenous retroviruses

The fate of most human endogenous retroviruses (HERVs) has been to undergo recombinational deletion. This process involves homologous recombination between the flanking long terminal repeats (LTRs) of a full-length element, leaving a relic structure in the genome termed a solo LTR. We examined loci in one family, HERV-K(HML2), and found that the deletion rate decreased markedly with age: the rate among recently integrated loci was almost 200-fold higher than that among loci whose insertion predated the divergence of humans and chimpanzees (8 x 10(-5) and 4 x 10(-7) recombinational deletion events per locus per generation, respectively). One hypothesis for this finding is that increasing mutational divergence between the flanking LTRs reduces the probability of homologous recombination and thus the rate of solo LTR formation. Consistent with this idea, we were able to replicate the observed rates by a simulation in which the probability of recombinational deletion was reduced 10-fold by a single mutation and 100-fold by any additional mutations. We also discuss the evidence for other factors that may influence the relationship between locus age and the rate of deletion, for example, host recombination rates and selection, and highlight the consequences of recombinational deletion for dating recent HERV integrations.