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.

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.

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.

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.

Retro-age: A unique epigenetic biomarker of aging captured by DNA methylation states of retroelements.

Reactivation of retroelements in the human genome has been linked to aging. However, whether the epigenetic state of specific retroelements can predict chronological age remains unknown. We provide evidence that locus-specific retroelement DNA methylation can be used to create retroelement-based epigenetic clocks that accurately measure chronological age in the immune system, across human tissues, and pan-mammalian species. We also developed a highly accurate retroelement epigenetic clock compatible with EPICv.2.0 data that was constructed from CpGs that did not overlap with existing first- and second-generation epigenetic clocks, suggesting a unique signal for epigenetic clocks not previously captured. We found retroelement-based epigenetic clocks were reversed during transient epigenetic reprogramming, accelerated in people living with HIV-1, and responsive to antiretroviral therapy. Our findings highlight the utility of retroelement-based biomarkers of aging and support a renewed emphasis on the role of retroelements in geroscience.

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).

Integrating human endogenous retroviruses into transcriptome-wide association studies highlights novel risk factors for major psychiatric conditions.

Human endogenous retroviruses (HERVs) are repetitive elements previously implicated in major psychiatric conditions, but their role in aetiology remains unclear. Here, we perform specialised transcriptome-wide association studies that consider HERV expression quantified to precise genomic locations, using RNA sequencing and genetic data from 792 post-mortem brain samples. In Europeans, we identify 1238 HERVs with expression regulated in cis, of which 26 represent expression signals associated with psychiatric disorders, with ten being conditionally independent from neighbouring expression signals. Of these, five are additionally significant in fine-mapping analyses and thus are considered high confidence risk HERVs. These include two HERV expression signatures specific to schizophrenia risk, one shared between schizophrenia and bipolar disorder, and one specific to major depressive disorder. No robust signatures are identified for autism spectrum conditions or attention deficit hyperactivity disorder in Europeans, or for any psychiatric trait in other ancestries, although this is likely a result of relatively limited statistical power. Ultimately, our study highlights extensive HERV expression and regulation in the adult cortex, including in association with psychiatric disorder risk, therefore providing a rationale for exploring neurological HERV expression in complex neuropsychiatric traits.

Exploring the roles of DNA methylation in the metal-reducing bacterium Shewanella oneidensis MR-1.

We performed whole-genome analyses of DNA methylation in Shewanella oneidensis MR-1 to examine its possible role in regulating gene expression and other cellular processes. Single-molecule real-time (SMRT) sequencing revealed extensive methylation of adenine (N6mA) throughout the genome. These methylated bases were located in five sequence motifs, including three novel targets for type I restriction/modification enzymes. The sequence motifs targeted by putative methyltranferases were determined via SMRT sequencing of gene knockout mutants. In addition, we found that S. oneidensis MR-1 cultures grown under various culture conditions displayed different DNA methylation patterns. However, the small number of differentially methylated sites could not be directly linked to the much larger number of differentially expressed genes under these conditions, suggesting that DNA methylation is not a major regulator of gene expression in S. oneidensis MR-1. The enrichment of methylated GATC motifs in the origin of replication indicates that DNA methylation may regulate genome replication in a manner similar to that seen in Escherichia coli. Furthermore, comparative analyses suggest that many Gammaproteobacteria, including all members of the Shewanellaceae family, may also utilize DNA methylation to regulate genome replication.

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.

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.

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.

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.

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.

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.

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.

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.

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.