The transcriptional landscape of endogenous retroelements delineates esophageal adenocarcinoma subtypes.

Most cancer types exhibit aberrant transcriptional activity, including derepression of retrotransposable elements (RTEs). However, the degree, specificity and potential consequences of RTE transcriptional activation may differ substantially among cancer types and subtypes. Representing one extreme of the spectrum, we characterize the transcriptional activity of RTEs in cohorts of esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) from the OCCAMS (Oesophageal Cancer Clinical and Molecular Stratification) consortium, and from TCGA (The Cancer Genome Atlas). We found exceptionally high RTE inclusion in the EAC transcriptome, driven primarily by transcription of genes incorporating intronic or adjacent RTEs, rather than by autonomous RTE transcription. Nevertheless, numerous chimeric transcripts straddling RTEs and genes, and transcripts from stand-alone RTEs, particularly KLF5- and SOX9-controlled HERVH proviruses, were overexpressed specifically in EAC. Notably, incomplete mRNA splicing and EAC-characteristic intronic RTE inclusion was mirrored by relative loss of the respective fully-spliced, functional mRNA isoforms, consistent with compromised cellular fitness. Defective RNA splicing was linked with strong transcriptional activation of a HERVH provirus on Chr Xp22.32 and defined EAC subtypes with distinct molecular features and prognosis. Our study defines distinguishable RTE transcriptional profiles of EAC, reflecting distinct underlying processes and prognosis, thus providing a framework for targeted studies.

Human endogenous retrovirus onco-exaptation counters cancer cell senescence through calbindin.

Increased levels and diversity of human endogenous retrovirus (HERV) transcription characterize most cancer types and are linked with disease outcomes. However, the underlying processes are incompletely understood. Here, we show that elevated transcription of HERVH proviruses predicted survival of lung squamous cell carcinoma (LUSC) and identified an isoform of CALB1, encoding calbindin, ectopically driven by an upstream HERVH provirus under the control of KLF5, as the mediator of this effect. HERVH-CALB1 expression was initiated in preinvasive lesions and associated with their progression. Calbindin loss in LUSC cell lines impaired in vitro and in vivo growth and triggered senescence, consistent with a protumor effect. However, calbindin also directly controlled the senescence-associated secretory phenotype (SASP), marked by secretion of CXCL8 and other neutrophil chemoattractants. In established carcinomas, CALB1-negative cancer cells became the dominant source of CXCL8, correlating with neutrophil infiltration and worse prognosis. Thus, HERVH-CALB1 expression in LUSC may display antagonistic pleiotropy, whereby the benefits of escaping senescence early during cancer initiation and clonal competition were offset by the prevention of SASP and protumor inflammation at later stages.

SARS-CoV-2-Host Chimeric RNA-Sequencing Reads Do Not Necessarily Arise From Virus Integration Into the Host DNA.

The human genome bears evidence of extensive invasion by retroviruses and other retroelements, as well as by diverse RNA and DNA viruses. High frequency of somatic integration of the RNA virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into the DNA of infected cells was recently suggested, based on a number of observations. One key observation was the presence of chimeric RNA-sequencing (RNA-seq) reads between SARS-CoV-2 RNA and RNA transcribed from human host DNA. Here, we examined the possible origin specifically of human-SARS-CoV-2 chimeric reads in RNA-seq libraries and provide alternative explanations for their origin. Chimeric reads were frequently detected also between SARS-CoV-2 RNA and RNA transcribed from mitochondrial DNA or episomal adenoviral DNA present in transfected cell lines, which was unlikely the result of SARS-CoV-2 integration. Furthermore, chimeric reads between SARS-CoV-2 RNA and RNA transcribed from nuclear DNA were highly enriched for host exonic, rather than intronic or intergenic sequences and often involved the same, highly expressed host genes. Although these findings do not rule out SARS-CoV-2 somatic integration, they nevertheless suggest that human-SARS-CoV-2 chimeric reads found in RNA-seq data may arise during library preparation and do not necessarily signify SARS-CoV-2 reverse transcription, integration in to host DNA and further transcription.

E3 ubiquitin ligase HECTD2 mediates melanoma progression and immune evasion.

The ubiquitin-proteasome system maintains protein homoeostasis, underpins the cell cycle, and is dysregulated in cancer. However, the role of individual E3 ubiquitin ligases, which mediate the final step in ubiquitin-mediated proteolysis, remains incompletely understood. Identified through screening for cancer-specific endogenous retroviral transcripts, we show that the little-studied E3 ubiquitin ligase HECTD2 exerts dominant control of tumour progression in melanoma. HECTD2 cell autonomously drives the proliferation of human and murine melanoma cells by accelerating the cell cycle. HECTD2 additionally regulates cancer cell production of immune mediators, initiating multiple immune suppressive pathways, which include the cyclooxygenase 2 (COX2) pathway. Accordingly, higher HECTD2 expression is associated with weaker anti-tumour immunity and unfavourable outcome of PD-1 blockade in human melanoma and counteracts immunity against a model tumour antigen in murine melanoma. This central, multifaceted role of HECTD2 in cancer cell-autonomous proliferation and in immune evasion may provide a single target for a multipronged therapy of melanoma.

A spontaneous genetically induced epiallele at a retrotransposon shapes host genome function.

Intracisternal A-particles (IAPs) are endogenous retroviruses (ERVs) responsible for most insertional mutations in the mouse. Full-length IAPs harbour genes flanked by long terminal repeats (LTRs). Here, we identify a solo LTR IAP variant (Iap5-1solo) recently formed in the inbred C57BL/6J mouse strain. In contrast to the C57BL/6J full-length IAP at this locus (Iap5-1full), Iap5-1solo lacks DNA methylation and H3K9 trimethylation. The distinct DNA methylation levels between the two alleles are established during preimplantation development, likely due to loss of KRAB zinc finger protein binding at the Iap5-1solo variant. Iap5-1solo methylation increases and becomes more variable in a hybrid genetic background yet is unresponsive to maternal dietary methyl supplementation. Differential epigenetic modification of the two variants is associated with metabolic differences and tissue-specific changes in adjacent gene expression. Our characterisation of Iap5-1 as a genetically induced epiallele with functional consequences establishes a new model to study transposable element repression and host-element co-evolution.

Our genome provides a complete set of genetic instructions for life. It begins by directing the growth and development of the embryo, and subsequently supports all the cells of the adult body in their daily routines. Yet approximately 10% of the DNA in mammalian genomes is made up of sequences originating from past retroviral infections, leaving a calling card in our genetic code. While these segments of retroviral DNA can no longer produce new infectious viruses, some of them retain the ability to copy themselves and jump into new parts of the genome. This can be problematic if they jump into and disrupt an important piece of genetic code. To protect against this, our bodies have evolved the ability to chemically strap down retroviral sequences by adding methyl groups to them and by modifying the proteins they are wrapped around. However, some of these endogenous retroviruses can dodge such so-called epigenetic modifications and disrupt genome function as a result. Studying a population of widely used inbred laboratory mice, Bertozzi et al. have identified a retroviral element that evades these epigenetic restraints. They discovered that some mice carry a full-length retroviral sequence while others have a shortened version of the same element. The shorter sequence lacked the repressive epigenetic marks found on the longer version, and this affected the expression of nearby genes. Moreover, the repressive marks could be partially restored by breeding the short-version mice with a distantly related mouse strain. Bertozzi et al. highlight an important issue for research using mouse models. Inbred laboratory mouse strains are assumed to have a fixed genetic code which allows scientists to conclude that any observed differences in their experiments are not a product of background genetic variation. However, this study emphasizes that this assumption is not guaranteed, and that hidden genetic diversity may be present in ostensibly genetically identical mice, with important implications for experimental outcomes. In addition, Bertozzi et al. provide a new mouse model for researchers to study the evolution and regulation of retroviral sequences and the impact of these processes on cell function.

Modulation of IL-6/STAT3 signaling axis in CD4+FOXP3- T cells represents a potential antitumor mechanism of azacitidine.

CD4+ T cells orchestrate immune responses and are actively engaged in shaping tumor immunity. Signal transducer and activator of transcription (STAT) signaling controls the epigenetic tuning of CD4+ T-cell differentiation and polarization, and perturbed STAT signaling networks in CD4+ T cells subvert antitumor immunity in malignancies. Azacitidine (AZA), the mainstay therapy for high-risk myelodysplastic syndromes (HR-MDS), affects CD4+ T-cell polarization and function, but whether this contributes to AZA efficacy is currently unknown. By using functional proteomic, transcriptomic, and mutational analyses in 73 HR-MDS patients undergoing AZA therapy, we demonstrate that responding patients exhibited a coordinated CD4+ T-cell immune response and downregulated the inflammatory cytokine signaling pathways in CD4+ T cells after AZA, in contrast to nonresponders who upregulated the same pathways. We further observed an AZA-mediated downregulation of intereukin-6 (IL-6)-induced STAT3 phosphorylation in CD4+FOXP3- conventional T cells (Tcons) that correlated independently with better response and survival, whereas it was also not associated with the mutation number and profile of the patients. The AZA-induced downregulation of IL-6/STAT3 axis in Tcons restored the STAT signaling architecture in CD4+ T-cell subsets, whereas STAT signaling networks remained disorganized in patients who upregulated IL-6/STAT3 activity in Tcons. Given the pivotal role of CD4+ T cells in adaptive immunity, our findings suggest that the downregulation of the IL-6/STAT3 pathway in Tcons potentially constitutes a previously unrecognized immune-mediated mechanism of action of AZA and sets the scene for developing rational strategies of AZA combinations with IL-6/STAT3 axis inhibitors.

Epigenetic therapy of myelodysplastic syndromes connects to cellular differentiation independently of endogenous retroelement derepression.

BACKGROUND: Myelodysplastic syndromes (MDS) and acute myeloid leukaemia (AML) are characterised by abnormal epigenetic repression and differentiation of bone marrow haematopoietic stem cells (HSCs). Drugs that reverse epigenetic repression, such as 5-azacytidine (5-AZA), induce haematological improvement in half of treated patients. Although the mechanisms underlying therapy success are not yet clear, induction of endogenous retroelements (EREs) has been hypothesised. METHODS: Using RNA sequencing (RNA-seq), we compared the transcription of EREs in bone marrow HSCs from a new cohort of MDS and chronic myelomonocytic leukaemia (CMML) patients before and after 5-AZA treatment with HSCs from healthy donors and AML patients. We further examined ERE transcription using the most comprehensive annotation of ERE-overlapping transcripts expressed in HSCs, generated here by de novo transcript assembly and supported by full-length RNA-seq. RESULTS: Consistent with prior reports, we found that treatment with 5-AZA increased the representation of ERE-derived RNA-seq reads in the transcriptome. However, such increases were comparable between treatment responses and failures. The extended view of HSC transcriptional diversity offered by de novo transcript assembly argued against 5-AZA-responsive EREs as determinants of the outcome of therapy. Instead, it uncovered pre-treatment expression and alternative splicing of developmentally regulated gene transcripts as predictors of the response of MDS and CMML patients to 5-AZA treatment. CONCLUSIONS: Our study identifies the developmentally regulated transcriptional signatures of protein-coding and non-coding genes, rather than EREs, as correlates of a favourable response of MDS and CMML patients to 5-AZA treatment and offers novel candidates for further evaluation.

ZFP57 regulation of transposable elements and gene expression within and beyond imprinted domains.

BACKGROUND: KRAB zinc finger proteins (KZFPs) represent one of the largest families of DNA-binding proteins in vertebrate genomes and appear to have evolved to silence transposable elements (TEs) including endogenous retroviruses through sequence-specific targeting of repressive chromatin states. ZFP57 is required to maintain the post-fertilization DNA methylation memory of parental origin at genomic imprints. Here we conduct RNA-seq and ChIP-seq analyses in normal and ZFP57 mutant mouse ES cells to understand the relative importance of ZFP57 at imprints, unique and repetitive regions of the genome. RESULTS: Over 80% of ZFP57 targets are TEs, however, ZFP57 is not essential for their repression. The remaining targets lie within unique imprinted and non-imprinted sequences. Though the loss of ZFP57 influences imprinted genes as expected, the majority of unique gene targets lose H3K9me3 with little effect on DNA methylation and very few exhibit alterations in expression. Comparison of ZFP57 mutants with DNA methyltransferase-deleted ES cells (TKO) identifies a remarkably similar pattern of H3K9me3 loss across the genome. These data define regions where H3K9me3 is secondary to DNA methylation and we propose that ZFP57 is the principal if not sole methylation-sensitive KZFP in mouse ES cells. Finally, we examine dynamics of DNA and H3K9 methylation during pre-implantation development and show that sites bound by ZFP57 in ES cells maintain DNA methylation and H3K9me3 at imprints and at non-imprinted regions on the maternally inherited chromosome throughout preimplantation development. CONCLUSION: Our analyses suggest the evolution of a rare DNA methylation-sensitive KZFP that is not essential for repeat silencing, but whose primary function is to maintain DNA methylation and repressive histone marks at germline-derived imprinting control regions.

Identification, Characterization, and Heritability of Murine Metastable Epialleles: Implications for Non-genetic Inheritance.

Generally repressed by epigenetic mechanisms, retrotransposons represent around 40% of the murine genome. At the Agouti viable yellow (Avy) locus, an endogenous retrovirus (ERV) of the intracisternal A particle (IAP) class retrotransposed upstream of the agouti coat-color locus, providing an alternative promoter that is variably DNA methylated in genetically identical individuals. This results in variable expressivity of coat color that is inherited transgenerationally. Here, a systematic genome-wide screen identifies multiple C57BL/6J murine IAPs with Avy epigenetic properties. Each exhibits a stable methylation state within an individual but varies between individuals. Only in rare instances do they act as promoters controlling adjacent gene expression. Their methylation state is locus-specific within an individual, and their flanking regions are enriched for CTCF. Variably methylated IAPs are reprogrammed after fertilization and re-established as variable loci in the next generation, indicating reconstruction of metastable epigenetic states and challenging the generalizability of non-genetic inheritance at these regions.

Interplay of cis and trans mechanisms driving transcription factor binding and gene expression evolution.

Noncoding regulatory variants play a central role in the genetics of human diseases and in evolution. Here we measure allele-specific transcription factor binding occupancy of three liver-specific transcription factors between crosses of two inbred mouse strains to elucidate the regulatory mechanisms underlying transcription factor binding variations in mammals. Our results highlight the pre-eminence of cis-acting variants on transcription factor occupancy divergence. Transcription factor binding differences linked to cis-acting variants generally exhibit additive inheritance, while those linked to trans-acting variants are most often dominantly inherited. Cis-acting variants lead to local coordination of transcription factor occupancies that decay with distance; distal coordination is also observed and may be modulated by long-range chromatin contacts. Our results reveal the regulatory mechanisms that interplay to drive transcription factor occupancy, chromatin state, and gene expression in complex mammalian cell states.