Endogenous retroelement expression in the gut microenvironment of people living with HIV-1.

BACKGROUND: Endogenous retroelements (EREs), including human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs), comprise almost half of the human genome. Our previous studies of the interferome in the gut suggest potential mechanisms regarding how IFNb may drive HIV-1 gut pathogenesis. As ERE activity is suggested to partake in type 1 immune responses and is incredibly sensitive to viral infections, we sought to elucidate underlying interactions between ERE expression and gut dynamics in people living with HIV-1 (PLWH). METHODS: ERE expression profiles from bulk RNA sequencing of colon biopsies and PBMC were compared between a cohort of PLWH not on antiretroviral therapy (ART) and uninfected controls. FINDINGS: 59 EREs were differentially expressed in the colon of PLWH when compared to uninfected controls (padj <0.05 and FC ≤ -1 or ≥ 1) [Wald's Test]. Of these 59, 12 EREs were downregulated in PLWH and 47 were upregulated. Colon expression of the ERE loci LTR19_12p13.31 and L1FLnI_1q23.1s showed significant correlations with certain gut immune cell subset frequencies in the colon. Furthermore L1FLnI_1q23.1s showed a significant upregulation in peripheral blood mononuclear cells (PBMCs) of PLWH when compared to uninfected controls suggesting a common mechanism of differential ERE expression in the colon and PBMC. INTERPRETATION: ERE activity has been largely understudied in genomic characterizations of human pathologies. We show that the activity of certain EREs in the colon of PLWH is deregulated, supporting our hypotheses that their underlying activity could function as (bio)markers and potential mediators of pathogenesis in HIV-1 reservoirs. FUNDING: US NIH grants NCI CA260691 (DFN) and NIAID UM1AI164559 (DFN).

Ribosomal profiling of human endogenous retroviruses in healthy tissues.

Human endogenous retroviruses (HERVs) are the germline embedded proviral fragments of ancient retroviral infections that make up roughly 8% of the human genome. Our understanding of HERVs in physiology primarily surrounds their non-coding functions, while their protein coding capacity remains virtually uncharacterized. Therefore, we applied the bioinformatic pipeline "hervQuant" to high-resolution ribosomal profiling of healthy tissues to provide a comprehensive overview of translationally active HERVs. We find that HERVs account for 0.1-0.4% of all translation in distinct tissue-specific profiles. Collectively, our study further supports claims that HERVs are actively translated throughout healthy tissues to provide sequences of retroviral origin to the human proteome.

Retroelement-Age Clocks: Epigenetic Age Captured by Human Endogenous Retrovirus and LINE-1 DNA methylation states.

Human endogenous retroviruses (HERVs), the remnants of ancient viral infections embedded within the human genome, and long interspersed nuclear elements 1 (LINE-1), a class of autonomous retrotransposons, are silenced by host epigenetic mechanisms including DNA methylation. The resurrection of particular retroelements has been linked to biological aging. Whether the DNA methylation states of locus specific HERVs and LINEs can be used as a biomarker of chronological age in humans remains unclear. We show that highly predictive epigenetic clocks of chronological age can be constructed from retroelement DNA methylation states in the immune system, across human tissues, and pan-mammalian species. We found retroelement epigenetic clocks were reversed during transient epigenetic reprogramming, accelerated in people living with HIV-1, responsive to antiretroviral therapy, and accurate in estimating long-term culture ages of human brain organoids. Our findings support the hypothesis of epigenetic dysregulation of retroelements as a potential contributor to the biological hallmarks of aging.

Locus specific human endogenous retroviruses reveal new lymphoma subtypes.

The heterogeneity of cancers are driven by diverse mechanisms underlying oncogenesis such as differential 'cell-of-origin' (COO) progenitors, mutagenesis, and viral infections. Classification of B-cell lymphomas have been defined by considering these characteristics. However, the expression and contribution of transposable elements (TEs) to B cell lymphoma oncogenesis or classification have been overlooked. We hypothesized that incorporating TE signatures would increase the resolution of B-cell identity during healthy and malignant conditions. Here, we present the first comprehensive, locus-specific characterization of TE expression in benign germinal center (GC) B-cells, diffuse large B-cell lymphoma (DLBCL), Epstein-Barr virus (EBV)-positive and EBV-negative Burkitt lymphoma (BL), and follicular lymphoma (FL). Our findings demonstrate unique human endogenous retrovirus (HERV) signatures in the GC and lymphoma subtypes whose activity can be used in combination with gene expression to define B-cell lineage in lymphoid malignancies, highlighting the potential of retrotranscriptomic analyses as a tool in lymphoma classification, diagnosis, and the identification of novel treatment groups.

Oncogenic Transformation Drives DNA Methylation Loss and Transcriptional Activation at Transposable Element Loci.

Transposable elements (TE) are typically silenced by DNA methylation and repressive histone modifications in differentiated healthy human tissues. However, TE expression increases in a wide range of cancers and is correlated with global hypomethylation of cancer genomes. We assessed expression and DNA methylation of TEs in fibroblast cells that were serially transduced with hTERT, SV40, and HRASR24C to immortalize and then transform them, modeling the different steps of the tumorigenesis process. RNA sequencing and whole-genome bisulfite sequencing were performed at each stage of transformation. TE expression significantly increased as cells progressed through transformation, with the largest increase in expression after the final stage of transformation, consistent with data from human tumors. The upregulated TEs were dominated by endogenous retroviruses [long terminal repeats (LTR)]. Most differentially methylated regions (DMR) in all stages were hypomethylated, with the greatest hypomethylation in the final stage of transformation. A majority of the DMRs overlapped TEs from the RepeatMasker database, indicating that TEs are preferentially demethylated. Many hypomethylated TEs displayed a concordant increase in expression. Demethylation began during immortalization and continued into transformation, while upregulation of TE transcription occurred in transformation. Numerous LTR elements upregulated in the model were also identified in The Cancer Genome Atlas datasets of breast, colon, and prostate cancer. Overall, these findings indicate that TEs, specifically endogenous retroviruses, are demethylated and transcribed during transformation. SIGNIFICANCE: Analysis of epigenetic and transcriptional changes in a transformation model reveals that transposable element expression and methylation are dysregulated during oncogenic transformation.

Endogenous Reverse Transcriptase Inhibition Attenuates TLR5-Mediated Inflammation.

Transposable elements (TEs) are mobile genomic sequences that encompass roughly 50% of the human genome. Class 1 TEs, or "retrotransposons," mobilize through the production of an RNA intermediate that is then reverse transcribed to form complementary DNA (cDNA) molecules capable of genomic reinsertion. While TEs are traditionally silenced to maintain genomic integrity, the recognition of immunostimulatory cues, such as those provided by microorganisms, drastically alters host transcription to induce the differential expression of TEs. Emerging evidence demonstrates that the inducible production of TE cDNA is not an inert phenomenon but instead has been coopted by host immunity to facilitate cross talk between host and constituents of the microbiota by agonizing intrinsic antiviral receptors. Here, we demonstrate that immunostimulation of toll-like receptor 4 (TLR4) with lipopolysaccharide (LPS) and TLR5 with bacterial flagella (FLA) alters the expression of retrotransposons, such as human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs). Next, we demonstrate that reverse transcriptase inhibitor (RTi) delivery ameliorates the acute production of the proinflammatory cytokine "tumor necrosis factor alpha" (TNF-α) in response to FLA in a monocytic cell line (THP-1). Collectively, our findings demonstrate that TLR5-mediated cross talk between the host and microbiota is partially dependent on the reverse transcription (RT) of retrotransposons. IMPORTANCE The microbiota is a potent reservoir of immunostimulatory and immunosuppressive motifs that fundamentally shape host immunity. Despite broad associations between microbial composition and host immunity, the mechanisms underlying host microbiota-induced immunoregulation remain poorly defined. Here, we demonstrate a novel mechanism by which motifs overabundant during dysbiotic conditions influence host immunity through the upregulation of endogenous RT to produce motifs that agonize antiviral receptors.

A single-cell transposable element atlas of human cell identity.

Single cell RNA sequencing (scRNA-seq) is revolutionizing the study of complex biological systems. However, most sequencing studies overlook the contribution of transposable element (TE) expression to the transcriptome. In both scRNA-seq and bulk tissue RNA sequencing (RNA-seq), quantification of TE expression is challenging due to repetitive sequence content and poorly characterized TE gene models. Here, we developed a tool and analysis pipeline for Single cell Transposable Element Locus Level Analysis of scRNA Sequencing (Stellarscope) that reassigns multi-mapped reads to specific genomic loci using an expectation-maximization algorithm. Using Stellarscope, we built an atlas of TE expression in human PBMCs. We found that locus-specific TEs delineate cell types and define new cell subsets not identified by standard mRNA expression profiles. Altogether, this study provides comprehensive insights into the influence of transposable elements in human biology.

Transposable elements may enhance antiviral resistance in HIV-1 elite controllers.

Less than 0.5% of people living with HIV-1 are elite controllers (ECs) - individuals who have a replication-competent viral reservoir in their CD4+ T cells but maintain undetectable plasma viremia without the help of antiretroviral therapy. While the EC CD4+ T cell transcriptome has been investigated for gene expression signatures associated with disease progression (or, in this case, a lack thereof), the expression and regulatory activity of transposable elements (TEs) in ECs has not been explored. Yet previous studies have established that TEs can directly impact the immune response to pathogens, including HIV-1. Thus, we hypothesize that the regulatory activities of TEs could contribute to the natural resistance of ECs against HIV-1. We perform a TE-centric analysis of previously published multi-omics data derived from EC individuals and other populations. We find that the CD4+ T cell transcriptome and retrotranscriptome of ECs are distinct from healthy controls, treated patients, and viremic progressors. However, there is a substantial level of transcriptomic heterogeneity among ECs. We categorize individuals with distinct chromatin accessibility and expression profiles into four clusters within the EC group, each possessing unique repertoires of TEs and antiviral factors. Notably, several TE families with known immuno-regulatory activity are differentially expressed among ECs. Their transcript levels in ECs positively correlate with their chromatin accessibility and negatively correlate with the expression of their KRAB zinc-finger (KZNF) repressors. This coordinated variation is seen at the level of individual TE loci likely acting or, in some cases, known to act as cis-regulatory elements for nearby genes involved in the immune response and HIV-1 restriction. Based on these results, we propose that the EC phenotype is driven in part by the reduced availability of specific KZNF proteins to repress TE-derived cis-regulatory elements for antiviral genes, thereby heightening their basal level of resistance to HIV-1 infection. Our study reveals considerable heterogeneity in the CD4+ T cell transcriptome of ECs, including variable expression of TEs and their KZNF controllers, that must be taken into consideration to decipher the mechanisms enabling HIV-1 control.

A field guide to endogenous retrovirus regulatory networks.

Germ cells are subject to exogenous retrovirus infections occasionally resulting in the genomic integration of retroviral gene sequences. These endogenized retroviruses (ERVs) are found throughout mammalian genomes. Initially thought to be inert, it is now appreciated that ERVs have often been co-opted for complex physiological processes. However, unregulated ERV transposition and expression are a threat to cellular fitness and genomic integrity, and so mammalian cells must control ERVs through pre- and post-transcriptional mechanisms. Here, we provide a field guide to the molecular machinery that identifies and silences ERVs.

How human endogenous retroviruses interact with the microbiota in health and disease.

The microbiota is a collective of microorganisms whose composition is intimately linked with human health and disease. Emerging evidence demonstrates that endogenous retroviruses facilitate crosstalk between the host and microbiota to fundamentally shape immunity.

Specific human endogenous retroviruses predict metastatic potential in uveal melanoma.

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

Naso-oropharyngeal microbiome from breast cancer patients diagnosed with COVID-19.

Due to immunosuppressive cancer therapies, cancer patients diagnosed with COVID-19 have a higher chance of developing severe symptoms and present a higher mortality rate in comparison to the general population. Here we show a comparative analysis of the microbiome from naso-oropharyngeal samples of breast cancer patients with respect to SARS-CoV-2 status and identified bacteria associated with symptom severity. Total DNA of naso-oropharyngeal swabs from 74 women with or without breast cancer, positive or negative for SARS-CoV-2 were PCR-amplified for 16S-rDNA V3 and V4 regions and submitted to massive parallel sequencing. Sequencing data were analyzed with QIIME2 and taxonomic identification was performed using the q2-feature-classifier QIIME2 plugin, the Greengenes Database, and amplicon sequence variants (ASV) analysis. A total of 486 different bacteria were identified. No difference was found in taxa diversity between sample groups. Cluster analysis did not group the samples concerning SARS-CoV-2 status, breast cancer diagnosis, or symptom severity. Three taxa (Pseudomonas, Moraxella, and Klebsiella,) showed to be overrepresented in women with breast cancer and positive for SARS-CoV-2 when compared to the other women groups, and five bacterial groups were associated with COVID-19 severity among breast cancer patients: Staphylococcus, Staphylococcus epidermidis, Scardovia, Parasegitibacter luogiensis, and Thermomonas. The presence of Staphylococcus in COVID-19 breast cancer patients may possibly be a consequence of nosocomial infection.

SARS-CoV-2 infection mediates differential expression of human endogenous retroviruses and long interspersed nuclear elements.

SARS-CoV-2 promotes an imbalanced host response that underlies the development and severity of COVID-19. Infections with viruses are known to modulate transposable elements (TEs), which can exert downstream effects by modulating host gene expression, innate immune sensing, or activities encoded by their protein products. We investigated the impact of SARS-CoV-2 infection on TE expression using RNA-Seq data from cell lines and from primary patient samples. Using a bioinformatics tool, Telescope, we showed that SARS-CoV-2 infection led to upregulation or downregulation of TE transcripts, a subset of which differed from cells infected with SARS, Middle East respiratory syndrome coronavirus (MERS-CoV or MERS), influenza A virus (IAV), respiratory syncytial virus (RSV), and human parainfluenza virus type 3 (HPIV3). Differential expression of key retroelements specifically identified distinct virus families, such as Coronaviridae, with unique retroelement expression subdividing viral species. Analysis of ChIP-Seq data showed that TEs differentially expressed in SARS-CoV-2 infection were enriched for binding sites for transcription factors involved in immune responses and for pioneer transcription factors. In samples from patients with COVID-19, there was significant TE overexpression in bronchoalveolar lavage fluid and downregulation in PBMCs. Thus, although the host gene transcriptome is altered by infection with SARS-CoV-2, the retrotranscriptome may contain the most distinctive features of the cellular response to SARS-CoV-2 infection.

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.

Expression of Human Endogenous Retroviruses in Systemic Lupus Erythematosus: Multiomic Integration With Gene Expression.

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies predominantly to nuclear material. Many aspects of disease pathology are mediated by the deposition of nucleic acid containing immune complexes, which also induce the type 1interferon response, a characteristic feature of SLE. Notably, SLE is remarkably heterogeneous, with a variety of organs involved in different individuals, who also show variation in disease severity related to their ancestries. Here, we probed one potential contribution to disease heterogeneity as well as a possible source of immunoreactive nucleic acids by exploring the expression of human endogenous retroviruses (HERVs). We investigated the expression of HERVs in SLE and their potential relationship to SLE features and the expression of biochemical pathways, including the interferon gene signature (IGS). Towards this goal, we analyzed available and new RNA-Seq data from two independent whole blood studies using Telescope. We identified 481 locus specific HERV encoding regions that are differentially expressed between case and control individuals with only 14% overlap of differentially expressed HERVs between these two datasets. We identified significant differences between differentially expressed HERVs and non-differentially expressed HERVs between the two datasets. We also characterized the host differentially expressed genes and tested their association with the differentially expressed HERVs. We found that differentially expressed HERVs were significantly more physically proximal to host differentially expressed genes than non-differentially expressed HERVs. Finally, we capitalized on locus specific resolution of HERV mapping to identify key molecular pathways impacted by differential HERV expression in people with SLE.

HAPHPIPE: Haplotype Reconstruction and Phylodynamics for Deep Sequencing of Intrahost Viral Populations.

Deep sequencing of viral populations using next-generation sequencing (NGS) offers opportunities to understand and investigate evolution, transmission dynamics, and population genetics. Currently, the standard practice for processing NGS data to study viral populations is to summarize all the observed sequences from a sample as a single consensus sequence, thus discarding valuable information about the intrahost viral molecular epidemiology. Furthermore, existing analytical pipelines may only analyze genomic regions involved in drug resistance, thus are not suited for full viral genome analysis. Here, we present HAPHPIPE, a HAplotype and PHylodynamics PIPEline for genome-wide assembly of viral consensus sequences and haplotypes. The HAPHPIPE protocol includes modules for quality trimming, error correction, de novo assembly, alignment, and haplotype reconstruction. The resulting consensus sequences, haplotypes, and alignments can be further analyzed using a variety of phylogenetic and population genetic software. HAPHPIPE is designed to provide users with a single pipeline to rapidly analyze sequences from viral populations generated from NGS platforms and provide quality output properly formatted for downstream evolutionary analyses.

Human Endogenous Retrovirus Expression Is Upregulated in the Breast Cancer Microenvironment of HIV Infected Women: A Pilot Study.

In people living with HIV (PLWH), chronic inflammation can lead to cancer initiation and progression, besides driving a dysregulated and diminished immune responsiveness. HIV infection also leads to increased transcription of Human Endogenous Retroviruses (HERVs), which could increase an inflammatory environment and create a tumor growth suppressive environment with high expression of pro-inflammatory cytokines. In order to determine the impact of HIV infection to HERV expression on the breast cancer microenvironment, we sequenced total RNA from formalin-fixed paraffin-embedded (FFPE) breast cancer samples of women HIV-negative and HIV-positive for transcriptome and retrotranscriptome analyses. We performed RNA extraction from FFPE samples, library preparation and total RNA sequencing (RNA-seq). The RNA-seq analysis shows 185 differentially expressed genes: 181 host genes (178 upregulated and three downregulated) and four upregulated HERV transcripts in HIV-positive samples. We also explored the impact of HERV expression in its neighboring breast cancer development genes (BRCA1, CCND1, NBS1/NBN, RAD50, KRAS, PI3K/PIK3CA) and in long non-coding RNA expression (AC060780.1, also known as RP11-242D8.1). We found a significant positive association of HERV expression with RAD50 and with AC060780.1, which suggest a possible role of HERV in regulating breast cancer genes from PLWH with breast cancer. In addition, we found immune system, extracellular matrix organization and metabolic signaling genes upregulated in HIV-positive breast cancer. In conclusion, our findings provide evidence of transcriptional and retrotranscriptional changes in breast cancer from PLWH compared to non-HIV breast cancer, including dysregulation of HERVs, suggesting an indirect effect of the virus on the breast cancer microenvironment.

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.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses.

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

Author Correction: A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.