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

Shared and unique patterns of autonomous human endogenous retrovirus loci transcriptomes in CD14 + monocytes from individuals with physical trauma or infection with COVID-19

Since previous studies have suggested that the RNAs of human endogenous retrovirus (HERV) might be involved in regulating innate immunity, it is important to investigate the HERV transcriptome patterns in innate immune cell types such as CD14 + monocytes. Using single cell RNA-seq datasets from resting or stimulated PBMCs mapped to 3,220 known discrete autonomous proviral HERV loci, we found individual-specific variation in HERV transcriptomes between HERV loci in CD14 + monocytes. Analysis of paired datasets from the same individual that were cultured in vitro with LPS or without (i.e. control) revealed 36 HERV loci in CD14 + monocytes that were detected only after activation. To extend our analysis to in vivo activated CD14 + monocytes, we used two scRNA-seq datasets from studies that had demonstrated activation of circulating CD14 + monocytes in patients with physical trauma or patients hospitalized with COVID-19 infections. For direct comparison between the trauma and COVID-19 datasets, we first analyzed 1.625 billion sequence reads from a composite pangenome control of 21 normal individuals. Comparison of the sequence read depth of HERV loci in the trauma or COVID-19 samples to the pangenome control revealed that 39 loci in the COVID-19 and 11 HERV loci in the trauma samples were significantly different (Mann-Whitney U test), with 9 HERV loci shared between the COVID-19 and trauma datasets. The capacity to compare HERV loci transcriptome patterns in innate immune cells, like CD14 + monocytes, across different pathological conditions will lead to greater understanding of the physiological role of HERV expression in health and disease.

Human Papillomavirus Carcinogenicity and the Need of New Perspectives: Thoughts from a Retrospective Analysis on Human Papillomavirus Outcomes Conducted at the Hospital University of Bari, Apulia, Italy, between 2011 and 2022

BACKGROUND

The current manuscript's aim was to determine the human papillomavirus (HPV) genotype-specific prevalence and distribution among individuals, males, and females, of different ages in the region of Apulia, Italy, highlighting the possible variables involved in the carcinogenicity mechanism. In addition, we proposed two hypothetical models of HPV's molecular dynamics, intending to clarify the impact of prevention and therapeutic strategies, explicitly modeled by recent survey data.

METHODS

We presented clinical data from 9647 participants tested for either high-risk (HR) or low-risk (LR) HPV at the affiliated Bari Policlinic University Hospital of Bari from 2011 to 2022. HPV DNA detection was performed using nested-polymerase chain reaction (PCR) and multiplex real-time PCR assay. Statistical analysis showed significant associations for all genders and ages and both HR- and LR-HPV types. A major number of significant pairwise associations were detected for the higher-risk types and females and lower-risk types and males.

RESULTS

The overall prevalence of HPV was 50.5% (n-4.869) vs. 49.5% (n-4.778) of the study population, of which 74.4% (n-3621) were found to be HPV high-risk (HR-HPV) genotypes and 57.7% (n-2.807) low-risk HPV (LR-HPV) genotypes, of which males were 58% and females 49%; the three most prevalent HR-HPV genotypes were HPV 53 (n707-15%), 16 (n704-14%), and 31 (n589-12%), and for LR-HPV, they were 42 (19%), 6 (16%), and 54 (13%); 56% of patients screened for HPV were ≤ 30 years old, 53% were between 31 and 40 years old, 46% were 41-50 and 51-60 years old, and finally, 44% of subjects were >60 years old.

CONCLUSIONS

Our study provided comprehensive epidemiological data on HPV prevalence and genotype distribution among 9647 participants, which could serve as a significant reference for clinical practice, and it implied the necessity for more effective screening methods for HPV carcinogenesis covering the use of more specific molecular investigations. Although this is a predominantly descriptive and epidemiological study, the data obtained offer not only a fairly unique trend compared to other studies of different realities and latitudes but also lead us to focus on the HPV infection within two groups of young people and adults and hypothesize the possible involvement of dysbiosis, stem cells, and the retrotransposition mechanism.

Expression and high levels of insertional polymorphism of an endogenous gammaretrovirus lineage in dogs

Despite the absence of a confirmed exogenously replicating retrovirus in Canis lupus familiaris (C. familiaris), past retroviral infections are evident in the genomes of living animals via the presence of endogenous retroviruses (ERVs). Although gammaretrovirus-like transcripts and enzyme activities were previously reported to be present in canine leukemias and lymphomas, those findings were not further explored. Initial analysis of the C. familiaris reference genome revealed a minor subset of one ERV lineage, classified as CfERV-Fc1(a), or Fc1(a) here, with features characteristic of recent integration, including the presence of ORFs and identical or nearly identical LTRs. Our previous analysis of whole genome sequence data belonging to extant Canidae revealed a burst of past infections in Canis ancestors resulting in numerous young, polymorphic, and highly intact loci now segregating in dogs. Here, we demonstrate the expression of full-length Fc1(a) proviruses in tissues collected from healthy animals and from animals with cancer. We observed significantly higher expression in samples of dogs with various cancer diagnoses when compared to samples from healthy dogs. Genotyping of insertionally polymorphic Fc1(a) loci identified candidate expressed proviruses and delineated distributions over sample groups. Collectively, the data show that Fc1(a) proviruses retain biological activity in the domestic dog and provides a means to examine potential genetic links with disease states in this species.

Confounding factors in profiling of locus-specific human endogenous retrovirus (HERV) transcript signatures in primary T cells using multi-study-derived datasets

BACKGROUND

Human endogenous retroviruses (HERV) are repetitive sequence elements and a substantial part of the human genome. Their role in development has been well documented and there is now mounting evidence that dysregulated HERV expression also contributes to various human diseases. While research on HERV elements has in the past been hampered by their high sequence similarity, advanced sequencing technology and analytical tools have empowered the field. For the first time, we are now able to undertake locus-specific HERV analysis, deciphering expression patterns, regulatory networks and biological functions of these elements. To do so, we inevitable rely on omics datasets available through the public domain. However, technical parameters inevitably differ, making inter-study analysis challenging. We here address the issue of confounding factors for profiling locus-specific HERV transcriptomes using datasets from multiple sources.

METHODS

We collected RNAseq datasets of CD4 and CD8 primary T cells and extracted HERV expression profiles for 3220 elements, resembling most intact, near full-length proviruses. Looking at sequencing parameters and batch effects, we compared HERV signatures across datasets and determined permissive features for HERV expression analysis from multiple-source data.

RESULTS

We could demonstrate that considering sequencing parameters, sequencing-depth is most influential on HERV signature outcome. Sequencing samples deeper broadens the spectrum of expressed HERV elements. Sequencing mode and read length are secondary parameters. Nevertheless, we find that HERV signatures from smaller RNAseq datasets do reliably reveal most abundantly expressed HERV elements. Overall, HERV signatures between samples and studies overlap substantially, indicating a robust HERV transcript signature in CD4 and CD8 T cells. Moreover, we find that measures of batch effect reduction are critical to uncover genic and HERV expression differences between cell types. After doing so, differences in the HERV transcriptome between ontologically closely related CD4 and CD8 T cells became apparent.

CONCLUSION

In our systematic approach to determine sequencing and analysis parameters for detection of locus-specific HERV expression, we provide evidence that analysis of RNAseq datasets from multiple studies can aid confidence of biological findings. When generating de novo HERV expression datasets we recommend increased sequence depth ( > = 100 mio reads) compared to standard genic transcriptome pipelines. Finally, batch effect reduction measures need to be implemented to allow for differential expression analysis.

Endogenous Retroviruses as Modulators of Innate Immunity

Endogenous retroviruses (ERVs), or LTR retrotransposons, are a class of transposable elements that are highly represented in mammalian genomes. Human ERVs (HERVs) make up roughly 8.3% of the genome and over the course of evolution, HERV elements underwent positive selection and accrued mutations that rendered them non-infectious; thereby, the genome could co-opt them into constructive roles with important biological functions. In the past two decades, with the help of advances in sequencing technology, ERVs are increasingly considered to be important components of the innate immune response. While typically silenced, expression of HERVs can be induced in response to traumatic, toxic, or infection-related stress, leading to a buildup of viral transcripts and under certain circumstances, proteins, including functionally active reverse transcriptase and viral envelopes. The biological activity of HERVs in the context of the innate immune response can be based on the functional effect of four major viral components: (1) HERV LTRs, (2) HERV-derived RNAs, (3) HERV-derived RNA:DNA duplexes and cDNA, and (4) HERV-derived proteins and ribonucleoprotein complexes. In this review, we will discuss the implications of HERVs in all four contexts in relation to innate immunity and their association with various pathological disease states.

Alternative and aberrant splicing of human endogenous retroviruses in cancer. What about head and neck? —A mini review

Human endogenous retroviruses (HERVs) are transcribed in many cancer types, including head and neck cancer. Because of accumulating mutations at proviral loci over evolutionary time, HERVs are functionally defective and cannot complete their viral life cycle. Despite that, HERV transcripts, including full-length viral RNAs and viral RNAs spliced as expected at the conventional viral splice sites, can be detected in particular conditions, such as cancer. Interestingly, non-viral–related transcription, including aberrant, non-conventionally spliced RNAs, has been reported as well. The role of HERV transcription in cancer and its contribution to oncogenesis or progression are still debated. Nonetheless, HERVs may constitute a suitable cancer biomarker or a target for therapy. Thus, ongoing research aims both to clarify the basic mechanisms underlying HERV transcription in cancer and to exploit its potential toward clinical application. In this mini-review, we summarize the current knowledge, the most recent findings, and the future perspectives of research on HERV transcription and splicing, with particular focus on head and neck cancer.

Human Endogenous Retrovirus K Envelope in Spinal Fluid of Amyotrophic Lateral Sclerosis Is Toxic

OBJECTIVE

Human endogenous retroviruses have been implicated in neurodegenerative diseases including amyotrophic lateral sclerosis (ALS). Expression of human endogenous retrovirus K (HERV-K) subtype HML-2 envelope (Env) in human neuronal cultures and in transgenic mice results in neurotoxicity and neurodegeneration, and mice expressing HML-2 Env display behavioral and neuromuscular characteristics resembling ALS. This study aims to characterize the neurotoxic properties of HML-2 Env.

METHODS

Env neurotoxicity was detected in ALS cerebrospinal fluid and confirmed using recombinant Env protein in a cell-based assay and a mouse model. The mechanism of neurotoxicity was assessed with immunoprecipitation followed by mass spectrometry and Western blot, and by screening a panel of inhibitors.

RESULTS

We observed that recombinant HML-2 Env protein caused neurotoxicity resulting in neuronal cell death, retraction of neurites, and decreased neuronal electrical activity. Injection of the Env protein into the brains of mice also resulted in neuronal cell death. HML-2 Env protein was also found in the cerebrospinal fluid of patients with sporadic ALS. The neurotoxic properties of the Env and the cerebrospinal fluid could be rescued with the anti-Env antibody. The Env was found to bind to CD98HC complexed to β1 integrin on the neuronal cell surface. Using a panel of compounds to screen for their ability to block Env-induced neurotoxicity, we found that several compounds were protective and are linked to the β1 integrin pathway.

INTERPRETATION

HERV-K Env is released extracellularly in ALS and causes neurotoxicity via a novel mechanism. Present results pave the way for new treatment strategies in sporadic ALS. ANN NEUROL 2022;92:545-561.

Locus-Specific Characterization of Human Endogenous Retrovirus Expression in Prostate, Breast, and Colon Cancers

Human endogenous retroviruses (HERV) have been implicated in a variety of diseases including cancers. Recent research implicates HERVs in epigenetic gene regulation. Here we utilize a recently developed bioinformatics tool for identifying HERV expression at the locus-specific level to identify differential expression of HERVs in matched tumor-normal RNA-sequencing (RNA-seq) data from The Cancer Genome Atlas. Data from 52 prostate cancer, 111 breast cancer, and 24 colon cancer cases were analyzed. Locus-specific analysis identified active HERV elements and differentially expressed HERVs in prostate cancer, breast cancer, and colon cancer. In addition, differentially expressed host genes were identified across prostate, breast, and colon cancer datasets, respectively, including several involved in demethylation and antiviral response pathways, supporting previous findings regarding the pathogenic mechanisms of HERVs. A majority of differentially expressed HERVs intersected protein coding genes or lncRNAs in each dataset, and a subset of differentially expressed HERVs intersected differentially expressed genes in prostate, breast, and colon cancers, providing evidence towards regulatory function. Finally, patterns in HERV expression were identified in multiple cancer types, with 155 HERVs differentially expressed in all three cancer types. This analysis extends previous results identifying HERV transcription in cancer RNA-seq datasets to a locus-specific level, and in doing so provides a foundation for future studies investigating the functional role of HERV in cancers and identifies a number of novel targets for cancer biomarkers and immunotherapy. SIGNIFICANCE: Expressed human endogenous retroviruses are mapped at locus-specific resolution and linked to specific pathways to identify potential biomarkers and therapeutic targets in prostate, breast, and colon cancers.

SMARCB1 deletion in atypical teratoid rhabdoid tumors results in human endogenous retrovirus K (HML-2) expression

Atypical Teratoid Rhabdoid Tumor (AT/RT) is a rare pediatric central nervous system cancer often characterized by deletion or mutation of SMARCB1, a tumor suppressor gene. In this study, we found that SMARCB1 regulates Human Endogenous Retrovirus K (HERV-K, subtype HML-2) expression. HML-2 is a repetitive element scattered throughout the human genome, encoding several intact viral proteins that have been associated with stem cell maintenance and tumorigenesis. We found HML-2 env expression in both the intracellular and extracellular compartments in all AT/RT cell lines (n = 4) and in 95% of AT/RT patient tissues (n = 37) evaluated. SMARCB1 knock-down in neural stem cells (NSCs) led to an upregulation of HML-2 transcription. We found that SMARCB1 binds adjacent to the HML-2 promoter, repressing its transcription via chromatin immunoprecipitation; restoration of SMARCB1 expression in AT/RT cell lines significantly downregulated HML-2 expression. Further, targeted downregulation of HML-2 transcription via CRISPR-dCas9 coupled with suppressor proteins led to cellular dispersion, decreased proliferation, and cell death in vitro. HML-2 knock-down with shRNA, siRNA, and CRISPR-dCas9 significantly decreased Ras expression as measured by qRT-PCR, suggesting that HML-2 modulates MAPK/ERK signaling in AT/RT cells. Overexpression of NRAS was sufficient to restore cellular proliferation, and MYC, a transcription factor downstream of NRAS, was bound to the HERV-K LTR significantly more in the absence of SMARCB1 expression in AT/RT cells. We show a mechanism by which these undifferentiated tumors remain pluripotent, and we demonstrate that their formation is aided by aberrant HML-2 activation, which is dependent on SMARCB1 and its interaction with MYC.

Human Endogenous Retrovirus (HERV)-K env Gene Knockout Affects Tumorigenic Characteristics of nupr1 Gene in DLD-1 Colorectal Cancer Cells

Human endogenous retroviruses (HERVs) are suggested to be involved in the development of certain diseases, especially cancers. To elucidate the function of HERV-K Env protein in cancers, an HERV-K env gene knockout (KO) in DLD-1 colorectal cancer cell lines was generated using the CRISPR-Cas9 system. Transcriptome analysis of HERV-K env KO cells using next-generation sequencing (NGS) was performed to identify the key genes associated with the function of HERV-K Env protein. The proliferation of HERV-K env KO cells was significantly reduced in in vitro culture as well as in in vivo nude mouse model. Tumorigenic characteristics, including migration, invasion, and tumor colonization, were also significantly reduced in HERV-K env KO cells. Whereas, they were enhanced in HERV-K env over-expressing DLD-1 cells. The expression of nuclear protein-1 (NUPR1), an ER-stress response factor that plays an important role in cell proliferation, migration, and reactive oxygen species (ROS) generation in cancer cells, significantly reduced in HERV-K env KO cells. ROS levels and ROS-related gene expression was also significantly reduced in HERV-K env KO cells. Cells transfected with NUPR1 siRNA (small interfering RNA) exhibited the same phenotype as HERV-K env KO cells. These results suggest that the HERV-K env gene affects tumorigenic characteristics, including cell proliferation, migration, and tumor colonization through NUPR1 related pathway.

Human Endogenous Retroviruses in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most aggressive and deadly brain tumor. It is primarily diagnosed in the elderly and has a 5-year survival rate of less than 6% even with the most aggressive therapies. The lack of biomarkers has made the development of immunotherapy for GBM challenging. Human endogenous retroviruses (HERVs) are a group of viruses with long terminal repeat (LTR) elements, which are believed to be relics from ancient viral infections. Recent studies have found that those repetitive elements play important roles in regulating various biological processes. The differentially expressed LTR elements from HERVs are potential biomarkers for immunotherapy to treat GBM. However, the understanding of the LTR element expression in GBM is greatly lacking.

METHODS

We obtained 1077.4 GB of sequencing data from public databases. These data were generated from 111 GBM tissue studies, 30 GBM cell lines studies, and 45 normal brain tissues studies. We analyzed repetitive elements that were differentially expressed in GBM and normal brain samples.

RESULTS

We found that 48 LTR elements were differentially expressed (p-value < 0.05) between GBM and normal brain tissues, of which 46 were HERV elements. Among these 46 elements, 34 significantly changed HERVs belong to the ERV1 superfamily. Furthermore, 43 out of the 46 differentially expressed HERV elements were upregulated.

CONCLUSION

Our results indicate significant differential expression of many HERV LTR elements in GBM and normal brain tissues. Expression levels of these elements could be developed as biomarkers for GBM treatments.

Upregulation of human endogenous retrovirus-K (HML-2) mRNAs in hepatoblastoma: Identification of potential new immunotherapeutic targets and biomarkers

PURPOSE

Hepatoblastoma is the most common liver malignancy in children. In order to advance therapy against hepatoblastoma, novel immunologic targets and biomarkers are needed. Our purpose in this investigation is to examine hepatoblastoma transcriptomes for the expression of a class of genomic elements known as Human Endogenous Retrovirus (HERVs). HERVs are abundant in the human genome and are biologically active elements that have been associated with multiple malignancies and proposed as immunologic targets in a subset of tumors. A sub-family of HERVs, HERV-K(HML-2) (HERV-K), have been shown to be tightly regulated in fetal development, making investigation of these elements in pediatric tumors paramount.

METHODS

We first created a HERVK-FASTA file utilizing 91 previously described HML-2 proviruses. We then concatenated the file onto the GRCh38.95 cDNA library from Ensembl. We used this reference database to evaluate existing RNA-seq data from 10 hepatoblastoma tumors and 3 normal liver controls (GEO accession ID: GSE8977575). Quantification and differential proviral expression analysis between hepatoblastoma and normal liver controls was performed using the pseudo-alignment program Salmon and DESeq2, respectively.

RESULTS

HERV-K mRNA was expressed in hepatoblastoma from multiple proviral loci. All expressed HERV-K proviral loci were upregulated in hepatoblastoma compared to normal liver controls. Five HERV-K proviruses (1q21.3, 3q27.2, 7q22.2, 12q24.33 and 17p13.1) were significantly differentially expressed (p-adjusted value <0.05, |log2 fold change| > 1.5) across conditions. The provirus at 17p13.1 had an approximately 300-fold increased expression in hepatoblastoma as compared to normal liver. This was in part due to the near absence of HERV-K mRNA at the 17p13.1 locus in fully differentiated liver samples.

CONCLUSIONS

Our investigation demonstrates that HERV-K is expressed from multiple loci in hepatoblastoma and that the expression is increased for several proviruses compared to normal liver controls. Our results suggest that HERV-K mRNA expression may be useful as a biomarker in hepatoblastoma, given the large differential expression profiles in hepatoblastoma, with very low mRNA levels in liver control samples.

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

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

The role of endogenous retroviruses-K in human cancer

It is known that human endogenous retroviruses (HERVs) constitute almost 8% of the human genome. Although the expression of HERVs from the human genome is tightly regulated, different exogenous and endogenous factors could trigger HERV activation. Aberrant expression of different HERVs may potentially cause a variety of diseases such as neurological and autoimmune diseases as well as cancer. It is suggested that HERV-K can induce cancer through different mechanisms that are discussed. The interplay between some tumor viruses and HERV-K seems to be a key player in progression of viral-associated cancers because elevated levels of Rec and Np9 proteins are observed in several cancers. The frequent over expression of HERV proteins and some specific antibodies in cancer cells could be considered as suitable prognostic and therapeutic biomarkers in diagnosis and treatment of cancers. The expression of HERV proteins in cancers and development of immune responses against them may also be used as targets for cancer immunotherapy. Further studies are warranted to evaluate the role of HERVs in cancer formation and use of different HERV proteins in developing new diagnostic and therapeutic approaches for cancer treatments.

Human endogenous retrovirus HERV-K(HML-2) RNA causes neurodegeneration through Toll-like receptors

Although human endogenous retroviruses (HERVs) represent a substantial proportion of the human genome and some HERVs, such as HERV-K(HML-2), are reported to be involved in neurological disorders, little is known about their biological function. We report that RNA from an HERV-K(HML-2) envelope gene region binds to and activates human Toll-like receptor (TLR) 8, as well as murine Tlr7, expressed in neurons and microglia, thereby causing neurodegeneration. HERV-K(HML-2) RNA introduced into the cerebrospinal fluid (CSF) of either C57BL/6 wild-type mice or APPPS1 mice, a mouse model for Alzheimer's disease (AD), resulted in neurodegeneration and microglia accumulation. Tlr7-deficient mice were protected against neurodegenerative effects but were resensitized toward HERV-K(HML-2) RNA when neurons ectopically expressed murine Tlr7 or human TLR8. Transcriptome data sets of human AD brain samples revealed a distinct correlation of upregulated HERV-K(HML-2) and TLR8 RNA expression. HERV-K(HML-2) RNA was detectable more frequently in CSF from individuals with AD compared with controls. Our data establish HERV-K(HML-2) RNA as an endogenous ligand for species-specific TLRs 7/8 and imply a functional contribution of human endogenous retroviral transcripts to neurodegenerative processes, such as AD.

Characterising a human endogenous retrovirus(HERV)-derived tumour-associated antigen: enriched RNA-Seq analysis of HERV-K(HML-2) in mantle cell lymphoma cell lines

BACKGROUND

The cell-surface attachment protein (Env) of the HERV-K(HML-2) lineage of endogenous retroviruses is a potentially attractive tumour-associated antigen for anti-cancer immunotherapy. The human genome contains around 100 integrated copies (called proviruses or loci) of the HERV-K(HML-2) virus and we argue that it is important for therapy development to know which and how many of these contribute to protein expression, and how this varies across tissues. We measured relative provirus expression in HERV-K(HML-2), using enriched RNA-Seq analysis with both short- and long-read sequencing, in three Mantle Cell Lymphoma cell lines (JVM2, Granta519 and REC1). We also confirmed expression of the Env protein in two of our cell lines using Western blotting, and analysed provirus expression data from all other relevant published studies.

RESULTS

Firstly, in both our and other reanalysed studies, approximately 10% of the transcripts mapping to HERV-K(HML-2) came from Env-encoding proviruses. Secondly, in one cell line the majority of the protein expression appears to come from one provirus (12q14.1). Thirdly, we find a strong tissue-specific pattern of provirus expression.

CONCLUSIONS

A possible dependency of Env expression on a single provirus, combined with the earlier observation that this provirus is not present in all individuals and a general pattern of tissue-specific expression among proviruses, has serious implications for future HERV-K(HML-2)-targeted immunotherapy. Further research into HERV-K(HML-2) as a possible tumour-associated antigen in blood cancers requires a more targeted, proteome-based, screening protocol that will consider these polymorphisms within HERV-K(HML-2). We include a plan (and necessary alignments) for such work.

Association between ALS and retroviruses: evidence from bioinformatics analysis

BACKGROUND

Emerging evidence suggests retroviruses play a role in the pathophysiology of amyotrophic lateral sclerosis (ALS). Specifically, activation of ancient viral genes embedded in the human genome is theorized to lead to motor neuron degeneration. We explore whether connections exist between ALS and retroviruses through protein interaction networks (PIN) and pathway analysis, and consider the potential roles in drug target discovery. Protein database and pathway/network analytical software including Ingenuity Pathway BioProfiler, STRING, and CytoScape were utilized to identify overlapping protein interaction networks and extract core cluster (s) of retroviruses and ALS.

RESULTS

Topological and statistical analysis of the ALS-PIN and retrovirus-PIN identified a shared, essential protein network and a core cluster with significant connections with both networks. The identified core cluster has three interleukin molecules IL10, Il-6 and IL-1B, a central apoptosis regulator TP53, and several major transcription regulators including MAPK1, ANXA5, SQSTM1, SREBF2, and FADD. Pathway enrichment analysis showed that this core cluster is associated with the glucocorticoid receptor singling and neuroinflammation signaling pathways. For confirmation purposes, we applied the same methodology to the West Nile and Polio virus, which demonstrated trivial connectivity with ALS, supporting the unique connection between ALS and retroviruses.

CONCLUSIONS

Bioinformatics analysis provides evidence to support pathological links between ALS and retroviral activation. The neuroinflammation and apoptotic regulation pathways are specifically implicated. The continuation and further analysis of large scale genome studies may prove useful in exploring genes important in retroviral activation and ALS, which may help discover new drug targets.

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

Characterization of Human Endogenous Retrovirus (HERV) expression within the transcriptomic landscape using RNA-seq is complicated by uncertainty in fragment assignment because of sequence similarity. We present Telescope, a computational software tool that provides accurate estimation of transposable element expression (retrotranscriptome) resolved to specific genomic locations. Telescope directly addresses uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We demonstrate the utility of our approach through single locus analysis of HERV expression in 13 ENCODE cell types. When examined at this resolution, we find that the magnitude and breadth of the retrotranscriptome can be vastly different among cell types. Furthermore, our approach is robust to differences in sequencing technology and demonstrates that the retrotranscriptome has potential to be used for cell type identification. We compared our tool with other approaches for quantifying transposable element (TE) expression, and found that Telescope has the greatest resolution, as it estimates expression at specific TE insertions rather than at the TE subfamily level. Telescope performs highly accurate quantification of the retrotranscriptomic landscape in RNA-seq experiments, revealing a differential complexity in the transposable element biology of complex systems not previously observed. Telescope is available at https://github.com/mlbendall/telescope.

Comparative Expression Analysis of Human Endogenous Retrovirus Elements in Peripheral Blood of Children with Specific Language Impairment

Specific language impairment (SLI) is a psychiatric condition with a complex etiology and a substantial genetic basis that affects children's verbal communication abilities. In this study, we examined the expression of five different human endogenous retrovirus elements (HERVs) in a cohort of 25 children with SLI and 25 healthy children in the control group. Human endogenous retrovirus elements, a diverse group of repetitive DNA sequences, can potentially cause considerable genetic heterogeneity. They had been integrated in the genome of our ancestors throughout evolution and now consist of about 8.0% of the human genome. Several HERV loci are transcribed in various cell types. Their expression in peripheral blood and in the brain is altered in many neurological and psychiatric diseases. To date, HERV expression profiles have never been studied in patients with SLI. This study aimed to elucidate differentially regulated human endogenous retroelements in peripheral blood of children with SLI, in comparison with healthy controls, through quantitative reverse tran-scription-polymerase chain reaction (qRT-PCR) methodology. Our results show that two genes: HERV-K (HLM-2) gag and HERV-P env were expressed at lower levels in the blood samples from SLI children in comparison with those in the control group.

A human endogenous retrovirus encoded protease potentially cleaves numerous cellular proteins

Background

A considerable portion of the human genome derives from retroviruses inherited over millions of years. Human endogenous retroviruses (HERVs) are usually severely mutated, yet some coding-competent HERVs exist. The HERV-K(HML-2) group includes evolutionarily young proviruses that encode typical retroviral proteins. HERV-K(HML-2) has been implicated in various human diseases because transcription is often upregulated and some of its encoded proteins are known to affect cell biology. HERV-K(HML-2) Protease (Pro) has received little attention so far, although it is expressed in some disease contexts and other retroviral proteases are known to process cellular proteins.

Results

We set out to identify human cellular proteins that are substrates of HERV-K(HML-2) Pro employing a modified Terminal Amine Isotopic Labeling of Substrates (TAILS) procedure. Thousands of human proteins were identified by this assay as significantly processed by HERV-K(HML-2) Pro at both acidic and neutral pH. We confirmed cleavage of a majority of selected human proteins in vitro and in co-expression experiments in vivo. Sizes of processing products observed for some of the tested proteins coincided with product sizes predicted by TAILS. Processed proteins locate to various cellular compartments and participate in diverse, often disease-relevant cellular processes. A limited number of HERV-K(HML-2) reference and non-reference loci appears capable of encoding active Pro.

Conclusions

Our findings from an approach combining TAILS with experimental verification of candidate proteins in vitro and in cultured cells suggest that hundreds of cellular proteins are potential substrates of HERV-K(HML-2) Pro. It is therefore conceivable that even low-level expression of HERV-K(HML-2) Pro affects levels of a diverse array of proteins and thus has a functional impact on cell biology and possible relevance for human diseases. Further studies are indicated to elucidate effects of HERV-K(HML-2) Pro expression regarding human substrate proteins, cell biology, and disease. The latter also calls for studies on expression of specific HERV-K(HML-2) loci capable of encoding active Pro. Endogenous retrovirus-encoded Pro activity may also be relevant for disease development in species other than human.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0178-z) contains supplementary material, which is available to authorized users.

Human endogenous retrovirus-K (HML-2): a comprehensive review

The human genome contains a large number of retroviral elements acquired over the process of evolution, some of which are specific to primates. However, as many of these are defective or silenced through epigenetic changes, they were historically considered "junk DNA" and their potential role in human physiology or pathological circumstances have been poorly studied. The most recently acquired, human endogenous retrovirus-K (HERV-K), has multiple copies in the human genome and some of them have complete open reading frames that are transcribed and translated, especially in early embryogenesis. Phylogenetically, HERV-K is considered a supergroup of viruses. One of the subtypes, termed HML-2, seems to be the most active and hence, it is the best studied. Aberrant expression of HML-2 in adult tissues has been associated with certain types of cancer and with neurodegenerative diseases. This review discusses the discovery of these viruses, their classification, structure, regulation and potential for replication, physiological roles, and their involvement in disease pathogenesis. Finally, it presents different therapeutic approaches being considered to target these viruses.

Transcriptional profiling of HERV-K(HML-2) in amyotrophic lateral sclerosis and potential implications for expression of HML-2 proteins

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder. About 90% of ALS cases are without a known genetic cause. The human endogenous retrovirus multi-copy HERV-K(HML-2) group was recently reported to potentially contribute to neurodegeneration and disease pathogenesis in ALS because of transcriptional upregulation and toxic effects of HML-2 Envelope (Env) protein. Env and other proteins are encoded by some transcriptionally active HML-2 loci. However, more detailed information is required regarding which HML-2 loci are transcribed in ALS, which of their proteins are expressed, and differences between the disease and non-disease states.

METHODS

For brain and spinal cord tissue samples from ALS patients and controls, we identified transcribed HML-2 loci by generating and mapping HML-2-specific cDNA sequences. We predicted expression of HML-2 env gene-derived proteins based on the observed cDNA sequences. Furthermore, we determined overall HML-2 transcript levels by RT-qPCR and investigated presence of HML-2 Env protein in ALS and control tissue samples by Western blotting.

RESULTS

We identified 24 different transcribed HML-2 loci. Some of those loci are transcribed at relatively high levels. However, significant differences in HML-2 loci transcriptional activities were not seen when comparing ALS and controls. Likewise, overall HML-2 transcript levels, as determined by RT-qPCR, were not significantly different between ALS and controls. Indeed, we were unable to detect full-length HML-2 Env protein in ALS and control tissue samples despite reasonable sensitivity. Rather our analyses suggest that a number of HML-2 protein variants other than full-length Env may potentially be expressed in ALS patients.

CONCLUSIONS

Our results expand and refine recent publications on HERV-K(HML-2) and ALS. Some of our results are in conflict with recent findings and call for further specific analyses. Our profiling of HML-2 transcription in ALS opens up the possibility that HML-2 proteins other than canonical full-length Env may have to be considered when studying the role of HML-2 in ALS disease.

Differentiation-Dependent Regulation of Human Endogenous Retrovirus K Sequences and Neighboring Genes in Germ Cell Tumor Cells

Under physiological conditions, most human endogenous retroviruses (HERVs) are transcriptionally silent. However, re-activation of HERVs is observed under pathological conditions like inflammation or cancer. In addition to expression of HERV sequences, an impact of HERV-loci on expression of adjacent genes has been suggested as probably important patho-physiological mechanism. A candidate for such a gene is PRODH (proline dehydrogenase 1), which is located on chromosome 22 adjacent to HERVK-24. Germ cell tumors (GCTs) are known to express high level of HERVK sequences. In addition, non-seminomatous GCT are useful models to study HERV expression in the context of differentiation since they reflect aspects of cellular development during embryogenesis and usually contain different cell types. This is due to the embryonal carcinoma (EC) cells, which are the stem cell component of GCT. They are pluripotent, show high expression of pluripotency markers like OCT4 and LIN28A and can differentiate into either somatic derivatives (teratoma cells) or choriocarcinoma or yolk-sac tumor cells reflecting extra-embryonal differentiation. OCT4 is lost upon differentiation. We used GCT derived cell lines of varying differentiation stages to analyze expression of HERVK and PRODH. Differentiation status and cellular relationship of GCT cells was determined using microarray analysis and western blotting of the embryonic pluripotency markers OCT4 and LIN28A. The highest expression of HERVK was found in undifferentiated EC cells, which retain a stem cell phenotype and express both OCT4 and LIN28. In contrast, the lowest expression of HERVK was observed in somatic differentiated GCT cells which also lack OCT4 and LIN28A whereas GCT cells with differentiation characteristics of yolk-sac tumor expressed LIN28A but not OCT4 and showed intermediate level of HERVK. A similar pattern was found for PRODH. Differentiation of EC cells by siRNA mediated knock-down of OCT4 or treatment with differentiation inducing medium decreased expression of HERVK and PRODH. Treatment of differentiated GCT cells with 5′-azacytidine and trichostatin A increased expression of HERVK and PRODH, indicating that epigenetic mechanisms are responsible for altered expression of these genes. Our data suggest that HERVK expression is dependent on cellular differentiation stages regulated by epigenetic mechanisms, which can also affect expression of neighboring genes.

Downregulation of Human Endogenous Retrovirus Type K (HERV-K) Viral env RNA in Pancreatic Cancer Cells Decreases Cell Proliferation and Tumor Growth

Purpose: We investigated the role of the human endogenous retrovirus type K (HERV-K) envelope (env) gene in pancreatic cancer.Experimental Design: shRNA was employed to knockdown (KD) the expression of HERV-K in pancreatic cancer cells.Results: HERV-K env expression was detected in seven pancreatic cancer cell lines and in 80% of pancreatic cancer patient biopsies, but not in two normal pancreatic cell lines or uninvolved normal tissues. A new HERV-K splice variant was discovered in several pancreatic cancer cell lines. Reverse transcriptase activity and virus-like particles were observed in culture media supernatant obtained from Panc-1 and Panc-2 cells. HERV-K viral RNA levels and anti-HERV-K antibody titers were significantly higher in pancreatic cancer patient sera (N = 106) than in normal donor sera (N = 40). Importantly, the in vitro and in vivo growth rates of three pancreatic cancer cell lines were significantly reduced after HERV-K KD by shRNA targeting HERV-K env, and there was reduced metastasis to lung after treatment. RNA-Seq results revealed changes in gene expression after HERV-K env KD, including RAS and TP53. Furthermore, downregulation of HERV-K Env protein expression by shRNA also resulted in decreased expression of RAS, p-ERK, p-RSK, and p-AKT in several pancreatic cancer cells or tumors.Conclusions: These results demonstrate that HERV-K influences signal transduction via the RAS-ERK-RSK pathway in pancreatic cancer. Our data highlight the potentially important role of HERV-K in tumorigenesis and progression of pancreatic cancer, and indicate that HERV-K viral proteins may be attractive biomarkers and/or tumor-associated antigens, as well as potentially useful targets for detection, diagnosis, and immunotherapy of pancreatic cancer. Clin Cancer Res; 23(19); 5892-911. ©2017 AACR.

A comprehensive hybridization model allows whole HERV transcriptome profiling using high density microarray

Background

Human endogenous retroviruses (HERVs) have received much attention for their implications in the etiology of many human diseases and their profound effect on evolution. Notably, recent studies have highlighted associations between HERVs expression and cancers (Yu et al., Int J Mol Med 32, 2013), autoimmunity (Balada et al., Int Rev Immunol 29:351–370, 2010) and neurological (Christensen, J Neuroimmune Pharmacol 5:326–335, 2010) conditions. Their repetitive nature makes their study particularly challenging, where expression studies have largely focused on individual loci (De Parseval et al., J Virol 77:10414–10422, 2003) or general trends within families (Forsman et al., J Virol Methods 129:16–30, 2005; Seifarth et al., J Virol 79:341–352, 2005; Pichon et al., Nucleic Acids Res 34:e46, 2006).

Methods

To refine our understanding of HERVs activity, we introduce here a new microarray, HERV-V3. This work was made possible by the careful detection and annotation of genomic HERV/MaLR sequences as well as the development of a new hybridization model, allowing the optimization of probe performances and the control of cross-reactions.

Results

HERV-V3 offers an almost complete coverage of HERVs and their ancestors (mammalian apparent LTR-retrotransposons, MaLRs) at the locus level along with four other repertoires (active LINE-1 elements, lncRNA, a selection of 1559 human genes and common infectious viruses). We demonstrate that HERV-V3 analytical performances are comparable with commercial Affymetrix arrays, and that for a selection of tissue/pathological specific loci, the patterns of expression measured on HERV-V3 is consistent with those reported in the literature.

Conclusions

Given its large HERVs/MaLRs coverage and additional repertoires, HERV-V3 opens the door to multiple applications such as enhancers and alternative promoters identification, biomarkers identification as well as the characterization of genes and HERVs/MaLRs modulation caused by viral infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3669-7) contains supplementary material, which is available to authorized users.

Restricting retrotransposons: a review

Retrotransposons have generated about 40 % of the human genome. This review examines the strategies the cell has evolved to coexist with these genomic "parasites", focussing on the non-long terminal repeat retrotransposons of humans and mice. Some of the restriction factors for retrotransposition, including the APOBECs, MOV10, RNASEL, SAMHD1, TREX1, and ZAP, also limit replication of retroviruses, including HIV, and are part of the intrinsic immune system of the cell. Many of these proteins act in the cytoplasm to degrade retroelement RNA or inhibit its translation. Some factors act in the nucleus and involve DNA repair enzymes or epigenetic processes of DNA methylation and histone modification. RISC and piRNA pathway proteins protect the germline. Retrotransposon control is relaxed in some cell types, such as neurons in the brain, stem cells, and in certain types of disease and cancer, with implications for human health and disease. This review also considers potential pitfalls in interpreting retrotransposon-related data, as well as issues to consider for future research.

Human endogenous retrovirus-K contributes to motor neuron disease

The role of human endogenous retroviruses (HERVs) in disease pathogenesis is unclear. We show that HERV-K is activated in a subpopulation of patients with sporadic amyotrophic lateral sclerosis (ALS) and that its envelope (env) protein may contribute to neurodegeneration. The virus was expressed in cortical and spinal neurons of ALS patients, but not in neurons from control healthy individuals. Expression of HERV-K or its env protein in human neurons caused retraction and beading of neurites. Transgenic animals expressing the env gene developed progressive motor dysfunction accompanied by selective loss of volume of the motor cortex, decreased synaptic activity in pyramidal neurons, dendritic spine abnormalities, nucleolar dysfunction, and DNA damage. Injury to anterior horn cells in the spinal cord was manifested by muscle atrophy and pathological changes consistent with nerve fiber denervation and reinnervation. Expression of HERV-K was regulated by TAR (trans-activation responsive) DNA binding protein 43, which binds to the long terminal repeat region of the virus. Thus, HERV-K expression within neurons of patients with ALS may contribute to neurodegeneration and disease pathogenesis.

Viral symbiosis and the holobiontic nature of the human genome

The human genome is a holobiontic union of the mammalian nuclear genome, the mitochondrial genome and large numbers of endogenized retroviral genomes. This article defines and explores this symbiogenetic pattern of evolution, looking at the implications for human genetics, epigenetics, embryogenesis, physiology and the pathogenesis of inborn errors of metabolism and many other diseases.

Discovery of unfixed endogenous retrovirus insertions in diverse human populations

Endogenous retroviruses (ERVs) have contributed to more than 8% of the human genome. The majority of these elements lack function due to accumulated mutations or internal recombination resulting in a solitary (solo) LTR, although members of one group of human ERVs (HERVs), HERV-K, were recently active with members that remain nearly intact, a subset of which is present as insertionally polymorphic loci that include approximately full-length (2-LTR) and solo-LTR alleles in addition to the unoccupied site. Several 2-LTR insertions have intact reading frames in some or all genes that are expressed as functional proteins. These properties reflect the activity of HERV-K and suggest the existence of additional unique loci within humans. We sought to determine the extent to which other polymorphic insertions are present in humans, using sequenced genomes from the 1000 Genomes Project and a subset of the Human Genome Diversity Project panel. We report analysis of a total of 36 nonreference polymorphic HERV-K proviruses, including 19 newly reported loci, with insertion frequencies ranging from <0.0005 to >0.75 that varied by population. Targeted screening of individual loci identified three new unfixed 2-LTR proviruses within our set, including an intact provirus present at Xq21.33 in some individuals, with the potential for retained infectivity.

Human endogenous retroviruses and cancer

Human endogenous retroviruses (HERVs) are retroviruses that infected human genome millions of years ago and have persisted throughout human evolution. About 8% of our genome is composed of HERVs, most of which are nonfunctional because of epigenetic control or deactivating mutations. However, a correlation between HERVs and human cancer has been described and many tumors, such as melanoma, breast cancer, germ cell tumors, renal cancer or ovarian cancer, express HERV proteins, mainly HERV-K (HML6) and HERV-K (HML2). Although the causative role of HERVs in cancer is controversial, data from animal models demonstrated that endogenous retroviruses are potentially oncogenic. HERV protein expression in human cells generates an immune response by activating innate and adaptive immunities. Some HERV-derived peptides have antigenic properties. For example, HERV-K (HML-6) encodes the HER-K MEL peptide recognized by CD8+ lymphocytes. In addition, HERVs are two-edged immunomodulators. HERVs show immunosuppressive activity. The presence of genomic retroviral elements in host-cell cytosol may activate an interferon type I response. Therefore, targeting HERVs through cellular vaccines or immunomodulatory drugs combined with checkpoint inhibitors is attracting interest because they could be active in human tumors.

Human endogenous retrovirus HERV-K(HML-2) activity in prostate cancer is dominated by a few loci

BACKGROUND

Increased expression of human endogenous retroviruses, especially HERV-K(HML-2) proviruses, has recently been associated with prostate carcinoma progression. In particular, a HML-2 locus in chromosome 22q11.23 (H22q) is upregulated in many cases. We therefore aimed at delineating the extent and repertoire of HML-2 transcription in prostate cancer tissues and cell lines and to define the transcription pattern and biological effects of H22q.

METHODS

Sanger and high throughput amplicon sequencing was used to define the repertoire of expressed HML-2 in a selected set of samples. qRT-PCR was used to quantify expression of selected proviruses in an extended set of prostate cancer tissues. Transcription factor binding sites (TFBS) were compared bioinformatically using the Transfac database. Expression of H22q was further characterized by siRNA-mediated knockdown, 5' RACE mapping of transcriptional start sites (TSS) and identification of splice sites. Functional effects of H22q knockdown were investigated by viability and apoptosis assays.

RESULTS

In addition to H22q, a limited number of other proviruses were found expressed by sequencing. Of these, provirus ERVK-5 and to a lesser degree ERVK-15 were frequently upregulated in prostate cancer. In contrast, expression of ERVK-24, predominant in germ cell tumors, was not detectable in prostatic tissues. While HML-2 LTRs contain binding sites for the androgen receptor and cofactors, no consistent differences in transcription factor binding sites were found between expressed and non-expressed proviruses. The H22q locus contains two 5'-LTRs of which the upstream LTR is predominantly used in prostatic cells, with an imprecise TSS. Splicing of H22q transcripts is complex, generating, among others, a transcript with an Np9-like ORF. Knockdown of H22q did not significantly affect proliferation or apoptosis of prostate cancer cells.

CONCLUSIONS

Our findings further underline that HML-2 expression is commonly highly tissue-specific. In prostate cancer, a limited number of loci become activated, especially H22q and ERVK-5. As expressed and non-expressed proviruses do not differ significantly in TFBS, tissue- and tumor-specific expression may be governed primarily by chromatin context. Overexpression of HML-2 H22q is more likely consequence than cause of prostate cancer progression.

Viral metagenomics in drug-naïve, first-onset schizophrenia patients with prominent negative symptoms

Although several studies suggest a virus or (endogenous) retrovirus involvement at the time of onset of schizophrenia, the unequivocal identification of one or more infectious agents, by means of an undirected catch-all technique, has never been conducted. In this study VIDISCA, a virus discovery method, was used in combination with Roche-454 high-throughput sequencing as a tool to determine the possible presence of viruses (known or unknown) in blood of first-onset drugs-naïve schizophrenic patients with prominent negative symptoms. Two viruses (the Anellovirus Torque Teno virus and GB virus C) were detected. Both viruses are commonly found in healthy individuals and no clear link with disease was ever established. Viruses from the family Anelloviridae were also identified in the control population (4.8%). Besides, one patient sample was positive for human endogenous retroviruses type K (HML-2) RNA but no specific predominant strain was detected, instead 119 different variants were found. In conclusion, these findings indicate no evidence for viral or endogenous retroviral involvement in sera at the time of onset of schizophrenia.

What's the host and what's the microbe? The Marjory Stephenson Prize Lecture 2015

The interchange between retroviruses and their hosts is an intimate one because retroviruses integrate proviral DNA into host chromosomal DNA as an obligate step in the replication cycle. This has resulted in the occasional transduction of host genes into retroviral genomes as oncogenes, and also led to the integration of viral genomes into the host germ line that gives rise to endogenous retroviruses. I shall reflect on the evolutionary consequences of these events for virus and host. Then, I shall discuss the emergence of non-viral infections of host origin, namely, how malignant cells can give rise to eukaryotic single cell 'parasites' that colonize new hosts and how these in turn have been colonized by host mitochondria.

Np9, a cellular protein of retroviral ancestry restricted to human, chimpanzee and gorilla, binds and regulates ubiquitin ligase MDM2

Humans and primates are long-lived animals with long reproductive phases. One factor that appears to contribute to longevity and fertility in humans, as well as to cancer-free survival, is the transcription factor and tumor suppressor p53, controlled by its main negative regulator MDM2. However, p53 and MDM2 homologs are found throughout the metazoan kingdom from Trichoplacidae to Hominidae. Therefore the question arises, if p53/MDM2 contributes to the shaping of primate features, then through which mechanisms. Previous findings have indicated that the appearances of novel p53-regulated genes and wild-type p53 variants during primate evolution are important in this context. Here, we report on another mechanism of potential relevance. Human endogenous retrovirus K subgroup HML-2 (HERV-K(HML-2)) type 1 proviral sequences were formed in the genomes of the predecessors of contemporary Hominoidea and can be identified in the genomes of Nomascus leucogenys (gibbon) up to Homo sapiens. We previously reported on an alternative splicing event in HERV-K(HML-2) type 1 proviruses that can give rise to nuclear protein of 9 kDa (Np9). We document here the evolution of Np9-coding capacity in human, chimpanzee and gorilla, and show that the C-terminal half of Np9 binds directly to MDM2, through a domain of MDM2 that is known to be contacted by various cellular proteins in response to stress. Np9 can inhibit the MDM2 ubiquitin ligase activity toward p53 in the cell nucleus, and can support the transactivation of genes by p53. Our findings point to the possibility that endogenous retrovirus protein Np9 contributes to the regulation of the p53-MDM2 pathway specifically in humans, chimpanzees and gorillas.

Human endogenous retroviruses sustain complex and cooperative regulation of gene-containing loci and unannotated megabase-sized regions

Background

Evidence suggests that some human endogenous retroviruses and endogenous retrovirus-like repeats (here collectively ERVs) regulate the expression of neighboring genes in normal and disease states; e.g. the human globin locus is regulated by an ERV9 that coordinates long-range gene switching during hematopoiesis and activates also intergenic transcripts. While complex transcription regulation is associated with integration of certain exogenous retroviruses, comparable regulation sustained by ERVs is less understood.

Findings

We analyzed ERV transcription using ERV9 consensus sequences and publically available RNA-sequencing, chromatin immunoprecipitation with sequencing (ChIP-seq) and cap analysis gene expression (CAGE) data from ENCODE. We discovered previously undescribed and advanced transcription regulation mechanisms in several human reference cell lines. We show that regulation by ERVs involves long-ranging activations including complex RNA splicing patterns, and transcription of large unannotated regions ranging in size from several hundred kb to around 1 Mb. Moreover, regulation was found to be cooperatively sustained in some loci by multiple ERVs and also non-LTR repeats.

Conclusion

Our analyses show that endogenous retroviruses sustain advanced transcription regulation in human cell lines, which shows similarities to complex insertional mutagenesis effects exerted by exogenous retroviruses. By exposing previously undescribed regulation effects, this study should prove useful for understanding fundamental transcription mechanisms resulting from evolutionary acquisition of retroviral sequence in the human genome.

Electronic supplementary material

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

The Landscape of long noncoding RNA classification

Advances in the depth and quality of transcriptome sequencing have revealed many new classes of long noncoding RNAs (lncRNAs). lncRNA classification has mushroomed to accommodate these new findings, even though the real dimensions and complexity of the noncoding transcriptome remain unknown. Although evidence of functionality of specific lncRNAs continues to accumulate, conflicting, confusing, and overlapping terminology has fostered ambiguity and lack of clarity in the field in general. The lack of fundamental conceptual unambiguous classification framework results in a number of challenges in the annotation and interpretation of noncoding transcriptome data. It also might undermine integration of the new genomic methods and datasets in an effort to unravel the function of lncRNA. Here, we review existing lncRNA classifications, nomenclature, and terminology. Then, we describe the conceptual guidelines that have emerged for their classification and functional annotation based on expanding and more comprehensive use of large systems biology-based datasets.

Modulation of human endogenous retrovirus (HERV) transcription during persistent and de novo HIV-1 infection

BACKGROUND

The human genome contains multiple LTR elements including human endogenous retroviruses (HERVs) that together account for approximately 8-9% of the genomic DNA. At least 40 different HERV groups have been assigned to three major HERV classes on the basis of their homologies to exogenous retroviruses. Although most HERVs are silenced by a variety of genetic and epigenetic mechanisms, they may be reactivated by environmental stimuli such as exogenous viruses and thus may contribute to pathogenic conditions. The objective of this study was to perform an in-depth analysis of the influence of HIV-1 infection on HERV activity in different cell types.

RESULTS

A retrovirus-specific microarray that covers major HERV groups from all three classes was used to analyze HERV transcription patterns in three persistently HIV-1 infected cell lines of different cellular origins and in their uninfected counterparts. All three persistently infected cell lines showed increased transcription of multiple class I and II HERV groups. Up-regulated transcription of five HERV taxa (HERV-E, HERV-T, HERV-K (HML-10) and two ERV9 subgroups) was confirmed by quantitative reverse transcriptase PCR analysis and could be reversed by knock-down of HIV-1 expression with HIV-1-specific siRNAs. Cells infected de novo by HIV-1 showed stronger transcriptional up-regulation of the HERV-K (HML-2) group than persistently infected cells of the same origin. Analysis of transcripts from individual members of this group revealed up-regulation of predominantly two proviral loci (ERVK-7 and ERVK-15) on chromosomes 1q22 and 7q34 in persistently infected KE37.1 cells, as well as in de novo HIV-1 infected LC5 cells, while only one single HML-2 locus (ERV-K6) on chromosome 7p22.1 was activated in persistently infected LC5 cells.

CONCLUSIONS

Our results demonstrate that HIV-1 can alter HERV transcription patterns of infected cells and indicate a correlation between activation of HERV elements and the level of HIV-1 production. Moreover, our results suggest that the effects of HIV-1 on HERV activity may be far more extensive and complex than anticipated from initial studies with clinical material.

Lack of evidence that avian oncogenic viruses are infectious for humans: a review

Chickens may be infected with three different oncogenic viruses: avian leukosis virus (ALV), reticuloendotheliosis virus (REV), and Marek's disease herpesvirus (MDV). Several epidemiological studies have suggested a link between these viruses and different types of cancer in people working in poultry processing plants and with multiple sclerosis. In this article, we analyze the epidemiological evidence that these viruses are causative agents for human cancer, followed by description of the relevant key characteristics of ALV, REV, and MDV. Finally, we discuss the biological evidence or lack thereof that avian tumor viruses are involved in the etiology of human cancer and multiple sclerosis (MS). The recent primary epidemiologic articles that we reviewed as examples were only hypothesis-generating studies examining massive numbers of risk factors for associations with various imprecise, non-viral-specific outcomes. The studies lacked precise evidence of exposure to the relevant viruses and the statistical methods failed to adjust for the large risks of false-positive claims. ALV subgroups A-D and J have been eradicated in the United States from the pure lines down to the parent stocks by the breeder companies, which have greatly reduced the incidence of infection in layer flocks and broilers. As a consequence, potential exposure of humans to these viruses has greatly diminished. Infection of humans working in processing plants with ALV-A and ALV-B is unlikely, because broilers are generally resistant to infection with these two subgroups. Moreover, these viruses enter cells by specific receptors present on chicken, but not on mammalian, cells. Infection of mammalian cell cultures or animals with ALV-A, ALV-B, and ALV-J has not been reported. Moreover, humans vaccinated with exogenous or endogenous ALV-contaminated vaccines against yellow fever, measles, and mumps did not become antibody- or virus-positive for ALV. The risks for human infection with REV are similarly limited. First of all, REV also has been eradicated from pure lines down to parent stock by breeder companies in the United States. Broilers can still become infected with REV through infection with fowl pox virus containing REV. However, there is no indication that REV can infect human cells. Low levels of antibodies to ALV and REV in human sera have been reported by a few groups. Absorption of sera with chicken antigens reduced the antibody titers, and there was no clear association with contacts with poultry. Possible cross-reactions with human endogenous or exogenous retroviruses were not considered in these publications. MDV is typically associated with infection of chickens, and almost all experimental data show that MDV cannot infect mammalian cells or animals, including nonhuman primates. One study reports the presence of MDV gD DNA in human sera, but this finding could not be confirmed by another group. A Medline search of the term "gene expression in human cancers" was negative for publications with avian retroviruses or MDV. In conclusion, there is no indication that avian oncogenic viruses are involved in human cancer or MS or even able to infect and replicate in humans.

HERV-K(HML-2) rec and np9 transcripts not restricted to disease but present in many normal human tissues

Background

Human endogenous retroviruses of the HERV-K(HML-2) group have been associated with the development of tumor diseases. Various HERV-K(HML-2) loci encode retrovirus-like proteins, and expression of such proteins is upregulated in certain tumor types. HERV-K(HML-2)-encoded Rec and Np9 proteins interact with functionally important cellular proteins and may contribute to tumor development. Though, the biological role of HERV-K(HML-2) transcription and encoded proteins in health and disease is less understood. We therefore investigated transcription specifically of HERV-K(HML-2) rec and np9 mRNAs in a panel of normal human tissues.

Results

We obtained evidence for rec and np9 mRNA being present in all examined 16 normal tissue types. A total of 18 different HERV-K(HML-2) loci were identified as generating rec or np9 mRNA, among them loci not present in the human reference genome and several of the loci harboring open reading frames for Rec or Np9 proteins. Our analysis identified additional alternative splicing events of HERV-K(HML-2) transcripts, some of them encoding variant Rec/Np9 proteins. We also identified a second HERV-K(HML-2) locus formed by L1-mediated retrotransposition that is likewise transcribed in various human tissues.

Conclusions

HERV-K(HML-2) rec and np9 transcripts from different HERV-K(HML-2) loci appear to be present in various normal human tissues. It is conceivable that Rec and Np9 proteins and variants of those proteins are part of the proteome of normal human tissues and thus various cell types. Transcription of HERV-K(HML-2) may thus also have functional relevance in normal human cell physiology.

Electronic supplementary material

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

Differential expression of HERV-K (HML-2) proviruses in cells and virions of the teratocarcinoma cell line Tera-1

Human endogenous retrovirus (HERV-K (HML-2)) proviruses are among the few endogenous retroviral elements in the human genome that retain coding sequence. HML-2 expression has been widely associated with human disease states, including different types of cancers as well as with HIV-1 infection. Understanding of the potential impact of this expression requires that it be annotated at the proviral level. Here, we utilized the high throughput capabilities of next-generation sequencing to profile HML-2 expression at the level of individual proviruses and secreted virions in the teratocarcinoma cell line Tera-1. We identified well-defined expression patterns, with transcripts emanating primarily from two proviruses located on chromosome 22, only one of which was efficiently packaged. Interestingly, there was a preference for transcripts of recently integrated proviruses, over those from other highly expressed but older elements, to be packaged into virions. We also assessed the promoter competence of the 5’ long terminal repeats (LTRs) of expressed proviruses via a luciferase assay following transfection of Tera-1 cells. Consistent with the RNASeq results, we found that the activity of most LTRs corresponded to their transcript levels.

Increased seroreactivity to human T cell lymphoma/leukemia virus-related endogenous sequence-1 Gag peptides in patients with human T cell lymphoma/leukemia virus myelopathy

Previously, we had shown that although only 8% of patients with large granular lymphocytic leukemia (LGLL) were infected with human T cell lymphoma/leukemia virus (HTLV)-2, almost half had antibodies to HTLV Gag and Env peptides. Herein, we investigated whether this could be due to cross-reactive antibodies to two homologous peptides in the Gag protein of the endogenous retrovirus HTLV-related endogenous sequence-1 (HRES-1). In addition, we had previously shown that patients with HTLV neurodegenerative diseases had increased seroreactivity to homologous HERV-K10 endogenous retrovirus peptides. Hence, in this study we also examined whether these patients had increased seroreactivity to the aforementioned HRES-1 Gag peptides. Sera from 100 volunteer blood donors (VBD), 53 patients with LGLL, 74 subjects with HTLV-1 or 2 infection (58 nonmyelopathy and 16 myelopathy), and 83 patients with multiple sclerosis (MS) were evaluated. The HTLV-positive myelopathy (HAM) patients had a statistically increased prevalence of antibodies to both HRES-1 Gag peptides (81%) vs. the VBD (0%), LGLL patients (13%), and MS patients (1%), and the HTLV-positive nonmyelopathy subjects (21%). The data suggest that cross-reactivity to HRES-1 peptides could be involved in the pathogenesis of HAM. The difference between the VBD and LGLL patients was also statistically significant, also suggesting a possible association in a minority of patients.

The distribution of insertionally polymorphic endogenous retroviruses in breast cancer patients and cancer-free controls

Background

Integration of retroviral DNA into a germ cell can result in a provirus that is transmitted vertically to the host’s offspring. In humans, such endogenous retroviruses (HERVs) comprise >8% of the genome. The HERV-K(HML-2) proviruses consist of ~90 elements related to mouse mammary tumor virus, which causes breast cancer in mice. A subset of HERV-K(HML-2) proviruses has some or all genes intact, and even encodes functional proteins, though a replication competent copy has yet to be observed. More than 10% of HML-2 proviruses are human-specific, having integrated subsequent to the Homo-Pan divergence, and, of these, 11 are currently known to be polymorphic in integration site with variable frequencies among individuals. Increased expression of the most recent HML-2 proviruses has been observed in tissues and cell lines from several types of cancer, including breast cancer, for which expression may provide a meaningful marker of the disease.

Results

In this study, we performed a case–control analysis to investigate the possible relationship between the genome-wide presence of individual polymorphic HML-2 proviruses with the occurrence of breast cancer. For this purpose, we screened 50 genomic DNA samples from individuals diagnosed with breast cancer or without history of the disease (n = 25 per group) utilizing a combination of locus-specific PCR screening, in silico analysis of HML-2 content within the reference human genome sequence, and high-resolution genomic hybridization in semi-dried agarose. By implementing this strategy, we were able to analyze the distribution of both annotated and previously undescribed polymorphic HML-2 proviruses within our sample set, and to assess their possible association with disease outcome.

Conclusions

In a case–control analysis of 50 humans with regard to breast cancer diagnosis, we found no significant difference in the prevalence of proviruses between groups, suggesting common polymorphic HML-2 proviruses are not associated with breast cancer. Our findings indicate a higher level of putatively novel HML-2 sites within the population, providing support for additional recent insertion events, implying ongoing, yet rare, activities. These findings do not rule out either the possibility of involvement of such proviruses in a subset of breast cancers, or their possible utility as tissue-specific markers of disease.

Regulation of the human endogenous retrovirus K (HML-2) transcriptome by the HIV-1 Tat protein

Approximately 8% of the human genome is made up of endogenous retroviral sequences. As the HIV-1 Tat protein activates the overall expression of the human endogenous retrovirus type K (HERV-K) (HML-2), we used next-generation sequencing to determine which of the 91 currently annotated HERV-K (HML-2) proviruses are regulated by Tat. Transcriptome sequencing of total RNA isolated from Tat- and vehicle-treated peripheral blood lymphocytes from a healthy donor showed that Tat significantly activates expression of 26 unique HERV-K (HML-2) proviruses, silences 12, and does not significantly alter the expression of the remaining proviruses. Quantitative reverse transcription-PCR validation of the sequencing data was performed on Tat-treated PBLs of seven donors using provirus-specific primers and corroborated the results with a substantial degree of quantitative similarity.

IMPORTANCE

The expression of HERV-K (HML-2) is tightly regulated but becomes markedly increased following infection with HIV-1, in part due to the HIV-1 Tat protein. The findings reported here demonstrate the complexity of the genome-wide regulation of HERV-K (HML-2) expression by Tat. This work also demonstrates that although HERV-K (HML-2) proviruses in the human genome are highly similar in terms of DNA sequence, modulation of the expression of specific proviruses in a given biological situation can be ascertained using next-generation sequencing and bioinformatics analysis.

Unfixed endogenous retroviral insertions in the human population

One lineage of human endogenous retroviruses (HERVs), HERV-K(HML2), is upregulated in many cancers, some autoimmune/inflammatory diseases, and HIV-infected cells. Despite 3 decades of research, it is not known if these viruses play a causal role in disease, and there has been recent interest in whether they can be used as immunotherapy targets. Resolution of both these questions will be helped by an ability to distinguish between the effects of different integrated copies of the virus (loci). Research so far has concentrated on the 20 or so recently integrated loci that, with one exception, are in the human reference genome sequence. However, this viral lineage has been copying in the human population within the last million years, so some loci will inevitably be present in the human population but absent from the reference sequence. We therefore performed the first detailed search for such loci by mining whole-genome sequences generated by next-generation sequencing. We found a total of 17 loci, and the frequency of their presence ranged from only 2 of the 358 individuals examined to over 95% of them. On average, each individual had six loci that are not in the human reference genome sequence. Comparing the number of loci that we found to an expectation derived from a neutral population genetic model suggests that the lineage was copying until at least ∼250,000 years ago.

IMPORTANCE About 5% of the human genome sequence is composed of the remains of retroviruses that over millions of years have integrated into the chromosomes of egg and/or sperm precursor cells. There are indications that protein expression of these viruses is higher in some diseases, and we need to know (i) whether these viruses have a role in causing disease and (ii) whether they can be used as immunotherapy targets in some of them. Answering both questions requires a better understanding of how individuals differ in the viruses that they carry. We carried out the first careful search for new viruses in some of the many human genome sequences that are now available thanks to advances in sequencing technology. We also compared the number that we found to a theoretical expectation to see if it is likely that these viruses are still replicating in the human population today.

HIV-1 interacts with human endogenous retrovirus K (HML-2) envelopes derived from human primary lymphocytes

Human endogenous retroviruses (HERVs) are viruses that have colonized the germ line and spread through vertical passage. Only the more recently acquired HERVs, such as the HERV-K (HML-2) group, maintain coding open reading frames. Expression of HERV-Ks has been linked to different pathological conditions, including HIV infection, but our knowledge on which specific HERV-Ks are expressed in primary lymphocytes currently is very limited. To identify the most expressed HERV-Ks in an unbiased manner, we analyzed their expression patterns in peripheral blood lymphocytes using Pacific Biosciences (PacBio) single-molecule real-time (SMRT) sequencing. We observe that three HERV-Ks (KII, K102, and K18) constitute over 90% of the total HERV-K expression in primary human lymphocytes of five different donors. We also show experimentally that two of these HERV-K env sequences (K18 and K102) retain their ability to produce full-length and posttranslationally processed envelope proteins in cell culture. We show that HERV-K18 Env can be incorporated into HIV-1 but not simian immunodeficiency virus (SIV) particles. Moreover, HERV-K18 Env incorporation into HIV-1 virions is dependent on HIV-1 matrix. Taken together, we generated high-resolution HERV-K expression profiles specific for activated human lymphocytes. We found that one of the most abundantly expressed HERV-K envelopes not only makes a full-length protein but also specifically interacts with HIV-1. Our findings raise the possibility that these endogenous retroviral Env proteins could directly influence HIV-1 replication.

IMPORTANCE

Here, we report the HERV-K expression profile of primary lymphocytes from 5 different healthy donors. We used a novel deep-sequencing technology (PacBio SMRT) that produces the long reads necessary to discriminate the complexity of HERV-K expression. We find that primary lymphocytes express up to 32 different HERV-K envelopes, and that at least two of the most expressed Env proteins retain their ability to make a protein. Importantly, one of them, the envelope glycoprotein of HERV-K18, is incorporated into HIV-1 in an HIV matrix-specific fashion. The ramifications of such interactions are discussed, as the possibility of HIV-1 target tissue broadening and immune evasion are considered.

Regulation of Mouse Retroelement MuERV-L/MERVL Expression by REX1 and Epigenetic Control of Stem Cell Potency

About half of the mammalian genome is occupied by DNA sequences that originate from transposable elements. Retrotransposons can modulate gene expression in different ways and, particularly retrotransposon-derived long terminal repeats, profoundly shape expression of both surrounding and distant genomic loci. This is especially important in pre-implantation development, during which extensive reprograming of the genome takes place and cells pass through totipotent and pluripotent states. At this stage, the main mechanism responsible for retrotransposon silencing, i.e., DNA methylation, is inoperative. A particular retrotransposon called muERV-L/MERVL is expressed during pre-implantation stages and contributes to the plasticity of mouse embryonic stem cells. This review will focus on the role of MERVL-derived sequences as controlling elements of gene expression specific for pre-implantation development, two-cell stage-specific gene expression, and stem cell pluripotency, the epigenetic mechanisms that control their expression, and the contributions of the pluripotency marker REX1 and the related Yin Yang 1 family of transcription factors to this regulation process.