TEtranscripts: a package for including transposable elements in differential expression analysis of RNA-seq datasets

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programmes and repeats in pluripotent cells

H3K9me3 heterochromatin, established by lysine methyltransferases (KMTs) and compacted by heterochromatin protein 1 (HP1) isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3 heterochromatin stability is presently limited to individual domains and DNA repeats. Here we engineered Suv39h2-knockout mouse embryonic stem cells to degrade remaining two H3K9me3 KMTs within 1 hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A 'binary switch' governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMT depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening and exit from pluripotency within 12 h. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3 heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

Repetitive element transcript accumulation is associated with inflammaging in humans

Chronic, low-grade inflammation increases with aging, contributing to functional declines and diseases that reduce healthspan. Growing evidence suggests that transcripts from repetitive elements (RE) in the genome contribute to this "inflammaging" by stimulating innate immune activation, but evidence of RE-associated inflammation with aging in humans is limited. Here, we present transcriptomic and clinical data showing that RE transcript levels are positively related to gene expression of innate immune sensors, and to serum interleukin 6 (a marker of systemic inflammation), in a large group of middle-aged and older adults. We also: (1) use transcriptomics and whole-genome bisulfite (methylation) sequencing to show that many RE may be hypomethylated with aging, and that aerobic exercise, a healthspan-extending intervention, reduces RE transcript levels and increases RE methylation in older adults; and (2) extend our findings in a secondary dataset demonstrating age-related changes in RE chromatin accessibility. Collectively, our data support the idea that age-related RE transcript accumulation may play a role in inflammaging in humans, and that RE dysregulation with aging may be due in part to upstream epigenetic changes.

Spatial transcriptomic landscape unveils immunoglobin-associated senescence as a hallmark of aging

To systematically characterize the loss of tissue integrity and organ dysfunction resulting from aging, we produced an in-depth spatial transcriptomic profile of nine tissues in male mice during aging. We showed that senescence-sensitive spots (SSSs) colocalized with elevated entropy in organizational structure and that the aggregation of immunoglobulin-expressing cells is a characteristic feature of the microenvironment surrounding SSSs. Immunoglobulin G (IgG) accumulated across the aged tissues in both male and female mice, and a similar phenomenon was observed in human tissues, suggesting the potential of the abnormal elevation of immunoglobulins as an evolutionarily conserved feature in aging. Furthermore, we observed that IgG could induce a pro-senescent state in macrophages and microglia, thereby exacerbating tissue aging, and that targeted reduction of IgG mitigated aging across various tissues in male mice. This study provides a high-resolution spatial depiction of aging and indicates the pivotal role of immunoglobulin-associated senescence during the aging process.

YY1 is a transcriptional activator of the mouse LINE-1 Tf subfamily

Long interspersed element type 1 (LINE-1, L1) is an active autonomous transposable element in human and mouse genomes. L1 transcription is controlled by an internal RNA polymerase II promoter in the 5' untranslated region (5'UTR) of a full-length L1. It has been shown that transcription factor YY1 binds to a conserved sequence at the 5' end of the human L1 5'UTR and primarily dictates where transcription initiates. Putative YY1-binding motifs have been predicted in the 5'UTRs of two distinct mouse L1 subfamilies, Tf and Gf. Using site-directed mutagenesis, in vitro binding and gene knockdown assays, we experimentally tested the role of YY1 in mouse L1 transcription. Our results indicate that Tf, but not Gf subfamily, harbors functional YY1-binding sites in 5'UTR monomers and YY1 functions as a transcriptional activator for the mouse Tf subfamily. Activation of Tf transcription by YY1 during early embryogenesis is also supported by a reanalysis of published zygotic knockdown data. Furthermore, YY1-binding motifs are solely responsible for the synergistic interaction between Tf monomers, consistent with a model wherein distant monomers act as enhancers for mouse L1 transcription. The abundance of YY1-binding sites in Tf elements also raise important implications for gene regulation across the genome.

Targeting heterochromatin eliminates chronic myelomonocytic leukemia malignant stem cells through reactivation of retroelements and immune pathways

Chronic myelomonocytic leukemia (CMML) is a severe myeloid malignancy affecting the elderly, for which therapeutic options are limited. DNA hypomethylating agents (HMAs) provide transient responses, failing to eradicate the malignant clone. Hematopoietic stem cell (HSC) aging involves heterochromatin reorganization, evidenced by alterations in histone marks H3K9me2 and H3K9me3. These repressive marks together with DNA methylation are essential for suppressing transposable elements (TEs). In solid cancers, the antitumor efficacy of HMAs involves the derepression of TEs, mimicking a state of viral infection. In this study, we demonstrate a significant disorganization of heterochromatin in CMML HSCs and progenitors (HSPCs) characterized by an increase in the repressive mark H3K9me2, mainly at the level of TEs, and a repression of immune and age-associated transcripts. Combining HMAs with G9A/GLP H3K9me2 methyltransferase inhibitors reactivates these pathways, selectively targeting mutated cells while preserving wild-type HSCs, thus offering new therapeutic avenues for this severe myeloid malignancy.

Retrotransposons are co-opted to activate hematopoietic stem cells and erythropoiesis

Hematopoietic stem cells (HSCs) and erythropoiesis are activated during pregnancy and after bleeding by the derepression of retrotransposons, including endogenous retroviruses and long interspersed nuclear elements. Retrotransposon transcription activates the innate immune sensors cyclic guanosine 3',5'-monophosphate-adenosine 5'-monophosphate synthase (cGAS) and stimulator of interferon (IFN) genes (STING), which induce IFN and IFN-regulated genes in HSCs, increasing HSC division and erythropoiesis. Inhibition of reverse transcriptase or deficiency for cGAS or STING had little or no effect on hematopoiesis in nonpregnant mice but depleted HSCs and erythroid progenitors in pregnant mice, reducing red blood cell counts. Retrotransposons and IFN-regulated genes were also induced in mouse HSCs after serial bleeding and, in human HSCs, during pregnancy. Reverse transcriptase inhibitor use was associated with anemia in pregnant but not in nonpregnant people, suggesting conservation of these mechanisms from mice to humans.

Widespread transposable element dysregulation in human aging brains with Alzheimer's disease

INTRODUCTION

Transposable element (TE) dysregulation is associated with neuroinflammation in Alzheimer's disease (AD) brains. Yet, TE quantitative trait loci (teQTL) have not been well characterized in human aged brains with AD.

METHODS

We leveraged large-scale bulk and single-cell RNA sequencing, whole-genome sequencing (WGS), and xQTL from three human AD brain biobanks to characterize TE expression dysregulation and experimentally validate AD-associated TEs using CRISPR interference (CRISPRi) assays in human induced pluripotent stem cell (iPSC)-derived neurons.

RESULTS

We identified 26,188 genome-wide significant TE expression QTLs (teQTLs) in human aged brains. Subsequent colocalization analysis of teQTLs with AD genetic loci identified AD-associated teQTLs and linked locus TEs. Using CRISPRi assays, we pinpointed a neuron-specific suppressive role of the activated short interspersed nuclear element (SINE; chr11:47608036-47608220) on expression of C1QTNF4 via reducing neuroinflammation in human iPSC-derived neurons.

DISCUSSION

We identified widespread TE dysregulation in human AD brains and teQTLs offer a complementary analytic approach to identify likely AD risk genes.

HIGHLIGHTS

Widespread transposable element (TE) dysregulations are observed in human aging brains with degrees of neuropathology, apolipoprotein E (APOE) genotypes, and neuroinflammation in Alzheimer's disease (AD). A catalog of TE quantitative trait loci (teQTLs) in human aging brains was created using matched RNA sequencing and whole-genome sequencing data. CRISPR interference assays reveal that an upregulated intergenic TE from the MIR family (chr11: 47608036-47608220) suppresses expression of its nearest anti-inflammatory gene C1QTNF4 in human induced pluripotent stem cell-derived neurons.

Multiomic characterization of RNA microenvironments by oligonucleotide-mediated proximity-interactome mapping

RNA molecules form complex networks of molecular interactions that are central to their function and to cellular architecture. But these interaction networks are difficult to probe in situ. Here, we introduce Oligonucleotide-mediated proximity-interactome MAPping (O-MAP), a method for elucidating the biomolecules near an RNA of interest, within its native context. O-MAP uses RNA-fluorescence in situ hybridization-like oligonucleotide probes to deliver proximity-biotinylating enzymes to a target RNA in situ, enabling nearby molecules to be enriched by streptavidin pulldown. This induces exceptionally precise biotinylation that can be easily optimized and ported to new targets or sample types. Using the noncoding RNAs 47S, 7SK and Xist as models, we develop O-MAP workflows for discovering RNA-proximal proteins, transcripts and genomic loci, yielding a multiomic characterization of these RNAs' subcellular compartments and new regulatory interactions. O-MAP requires no genetic manipulation, uses exclusively off-the-shelf parts and requires orders of magnitude fewer cells than established methods, making it accessible to most laboratories.

Exogenous Janus Kinase 617 Codon Influences Small Noncoding RNAs and Gene Expression in Ba/F3 Cells

ABSTRACT

Myeloproliferative neoplasms (MPNs) are blood cancers caused by mutations that originate from hematopoietic stem cells. More than 50%-90% of MPN patients had a dominant negative valine (V) to phenylalanine (F) mutation at the Janus kinase 617 codon (JAK2V617F) within the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) signaling pathway; however, this mutation was also found in a high percentage of the general population, its penetrance varied, and its onset was shown to be polygenic. Consequently, it is still unknown what molecular mechanism underlies the MPN transformation produced by JAK2V617F. Patients with MPN have been shown to have dysregulation of noncoding RNAs, such as microRNA (miRNA) and PIWI-interacting RNA (piRNA), although there is not any concrete proof that JAK2V617F alone is responsible for the aberrant regulation of miRNA and piRNA. Human wild type versus V617F-mutated JAK2 are expressed in mouse Ba/F3 cells, and the expressed small and total RNAs were subjected to next generation sequencing analysis to determine the direct induction. Differentially expressed miRNAs, gene expression, and transcript and gene variations were found between exogenously expressed JAK2 and JAK2V617F in Ba/F3 cells. The differently expressed variations contained enriched transposable elements and piRNAs, indicating a rearranged epigenome. The results of the pathway analysis show that the transformation that self-validated the chosen sequencing target genes is impacted by the JAK-STAT pathway. The induction route is functionally conserved, according to exogenously produced miRNA and gene expression. These results may clarify how the JAK2V617F induces transformation.

Genetic and Developmental Divergence in the Neural Crest Program between Cichlid Fish Species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic program in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to interspecific morphological differences, such as craniofacial structures and pigmentation. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared with teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes-particularly those controlled by sox10s-are involved in generating morphological diversification between species and lays the groundwork for further investigations into the mechanisms underpinning vertebrate NC diversification.

Natural Diversity of Heat-Induced Transcription of Retrotransposons in Arabidopsis thaliana

Transposable elements (TEs) are major components of plant genomes, profoundly impacting the fitness of their hosts. However, technical bottlenecks have long hindered our mechanistic understanding of TEs. Using RNA-Seq and long-read sequencing with Oxford Nanopore Technologies' (ONT) direct cDNA sequencing, we analyzed the heat-induced transcription of TEs in three natural accessions of Arabidopsis thaliana (Cvi-0, Col-0, and Ler-1). In addition to the well-studied ONSEN retrotransposon family, we confirmed Copia-35 as a second heat-responsive retrotransposon family with particularly high activity in the relict accession Cvi-0. Our analysis revealed distinct expression patterns of individual TE copies and suggest different mechanisms regulating the GAG protein production in the ONSEN versus Copia-35 families. In addition, analogously to ONSEN, Copia-35 activation led to the upregulation of flanking genes such as APUM9 and potentially to the quantitative modulation of flowering time. ONT data allowed us to test the extent to which read-through formation is important in the regulation of adjacent genes. Unexpectedly, our results indicate that for both families, the upregulation of flanking genes is not predominantly directly initiated by transcription from their 3' long terminal repeats. These findings highlight the intraspecific expressional diversity linked to retrotransposon activation under stress.

IRescue: uncertainty-aware quantification of transposable elements expression at single cell level

Transposable elements (TEs) are mobile DNA repeats known to shape the evolution of eukaryotic genomes. In complex organisms, they exhibit tissue-specific transcription. However, understanding their role in cellular diversity across most tissues remains a challenge, when employing single-cell RNA sequencing (scRNA-seq), due to their widespread presence and genetic similarity. To address this, we present IRescue (Interspersed Repeats single-cell quantifier), a software capable of estimating the expression of TE subfamilies at the single-cell level. IRescue incorporates a unique UMI deduplication algorithm to rectify sequencing errors and employs an Expectation-Maximization procedure to effectively redistribute the counts of multi-mapping reads. Our study showcases the precision of IRescue through analysis of both simulated and real single cell and nuclei RNA-seq data from human colorectal cancer, brain, skin aging, and PBMCs during SARS-CoV-2 infection and recovery. By linking the expression patterns of TE signatures to specific conditions and biological contexts, we unveil insights into their potential roles in cellular heterogeneity and disease progression.

Transposable element small and long RNAs in aging brains and implications in Huntington's and Parkinson's disease

Transposable Elements (TEs) are implicated in aging and neurodegenerative disorders, but the impact of brain TE RNA dynamics on these phenomena is not fully understood. Therefore, we quantified TE RNA changes in aging post-mortem human and mouse brains and in the neurodegenerative disorders Huntington's Disease (HD) and Parkinson's Disease (PD). We tracked TE small RNAs (smRNAs) expression landscape to assess the relationship to the active processing from TE long RNAs (lnRNAs). Human brain transcriptomes from the BrainSpan Atlas displayed a significant shift of TE smRNA patterns at age 20 years, whereas aging mouse brains lacked any such marked change, despite clear shift in aging-associated mRNA levels. Human frontal cortex displayed pronounced sense TE smRNAs during aging with a negative relationship between the TE smRNAs and lnRNAs indicative of age associated regulatory effects. Our analysis revealed TE smRNAs dysregulation in HD, while PD showed a stronger impact on TE lnRNAs, potentially correlating with the early average age of death for HD relative to PD. Furthermore, TE-silencing factor TRIM28 was down-regulated only in aging human brains, possibly explaining the lack of substantial TE RNA changes in aging mouse brains. Our study suggests brain TE RNAs may serve as novel biomarkers of human brain aging and neurodegenerative disorders.

Targeting ZNF638 activates antiviral immune responses and potentiates immune checkpoint inhibition in glioblastoma

Viral mimicry refers to the activation of innate anti-viral immune responses due to the induction of endogenous retroelement (RE) expression. Viral mimicry has been previously described to augment anti-tumor immune responses and sensitize solid tumors to immunotherapy including colorectal cancer, melanoma, and clear renal cell carcinoma. Here, we found that targeting a novel, master epigenetic regulator, Zinc Finger Protein 638 (ZNF638), induces viral mimicry in glioblastoma (GBM) preclinical models and potentiates immune checkpoint inhibition (ICI). ZNF638 recruits the HUSH complex, which precipitates repressive H3K9me3 marks on endogenous REs. In GBM, ZNF638 is associated with marked locoregional immunosuppressive transcriptional signatures, reduced endogenous RE expression and poor immune cell infiltration (CD8 + T-cells, dendritic cells). ZNF638 knockdown decreased H3K9-trimethylation, increased cytosolic dsRNA and activated intracellular dsRNA-signaling cascades (RIG-I, MDA5 and IRF3). Furthermore, ZNF638 knockdown upregulated antiviral immune programs and significantly increased PD-L1 immune checkpoint expression in patient-derived GBM neurospheres and diverse murine models. Importantly, targeting ZNF638 sensitized mice to ICI in syngeneic murine orthotopic models through innate interferon signaling. This response was recapitulated in recurrent GBM (rGBM) samples with radiographic responses to checkpoint inhibition with widely increased expression of dsRNA, PD-L1 and perivascular CD8 cell infiltration, suggesting dsRNA-signaling may mediate response to immunotherapy. Finally, we showed that low ZNF638 expression was a biomarker of clinical response to ICI and improved survival in rGBM patients and melanoma patients. Our findings suggest that ZNF638 could serve as a target to potentiate immunotherapy in gliomas.

MATES: a deep learning-based model for locus-specific quantification of transposable elements in single cell

Transposable elements (TEs) are crucial for genetic diversity and gene regulation. Current single-cell quantification methods often align multi-mapping reads to either 'best-mapped' or 'random-mapped' locations and categorize them at the subfamily levels, overlooking the biological necessity for accurate, locus-specific TE quantification. Moreover, these existing methods are primarily designed for and focused on transcriptomics data, which restricts their adaptability to single-cell data of other modalities. To address these challenges, here we introduce MATES, a deep-learning approach that accurately allocates multi-mapping reads to specific loci of TEs, utilizing context from adjacent read alignments flanking the TE locus. When applied to diverse single-cell omics datasets, MATES shows improved performance over existing methods, enhancing the accuracy of TE quantification and aiding in the identification of marker TEs for identified cell populations. This development facilitates the exploration of single-cell heterogeneity and gene regulation through the lens of TEs, offering an effective transposon quantification tool for the single-cell genomics community.

cGAS deficient mice display premature aging associated with de-repression of LINE1 elements and inflammation

Aging-associated inflammation, or 'inflammaging" is a driver of multiple age-associated diseases. Cyclic GMP-AMP Synthase (cGAS) is a cytosolic DNA sensor that functions to activate interferon response upon detecting viral DNA in the cytoplasm. cGAS contributes to inflammaging by responding to endogenous signals such as damaged DNA or LINE1 (L1) cDNA which forms in aged cells. While cGAS knockout mice are viable their aging has not been examined. Unexpectedly, we found that cGAS knockout mice exhibit accelerated aging phenotype associated with induction of inflammation. Transcription of L1 elements was increased in both cGAS knockout mice and in cGAS siRNA knockdown cells associated with high levels of cytoplasmic L1 DNA and expression of ORF1 protein. Cells from cGAS knockout mice showed increased chromatin accessibility and decreased DNA methylation on L1 transposons. Stimulated emission depletion microscopy (STED) showed that cGAS forms nuclear condensates that co-localize with H3K9me3 heterochromatin marks, and H3K9me3 pattern is disrupted in cGAS knockout cells. Taken together these results suggest a previously undescribed role for cGAS in maintaining heterochromatin on transposable elements. We propose that loss of cGAS leads to loss of chromatin organization, de-repression of transposable elements and induction of inflammation resulting in accelerated aging.

ZIC2 and ZIC3 promote SWI/SNF recruitment to safeguard progression towards human primed pluripotency

The primed epiblast acts as a transitional stage between the relatively homogeneous naïve epiblast and the gastrulating embryo. Its formation entails coordinated changes in regulatory circuits driven by transcription factors and epigenetic modifications. Using a multi-omic approach in human embryonic stem cell models across the spectrum of peri-implantation development, we demonstrate that the transcription factors ZIC2 and ZIC3 have overlapping but essential roles in opening primed-specific enhancers. Together, they are essential to facilitate progression to and maintain primed pluripotency. ZIC2/3 accomplish this by recruiting SWI/SNF to chromatin and loss of ZIC2/3 or degradation of SWI/SNF both prevent enhancer activation. Loss of ZIC2/3 also results in transcriptome changes consistent with perturbed Polycomb activity and a shift towards the expression of genes linked to differentiation towards the mesendoderm. Additionally, we find an intriguing dependency on the transcriptional machinery for sustained recruitment of ZIC2/3 over a subset of primed-hESC specific enhancers. Taken together, ZIC2 and ZIC3 regulate highly dynamic lineage-specific enhancers and collectively act as key regulators of human primed pluripotency.

GNAS knockout potentiates HDAC3 inhibition through viral mimicry-related interferon responses in lymphoma

Despite selective HDAC3 inhibition showing promise in a subset of lymphomas with CREBBP mutations, wild-type tumors generally exhibit resistance. Here, using unbiased genome-wide CRISPR screening, we identify GNAS knockout (KO) as a sensitizer of resistant lymphoma cells to HDAC3 inhibition. Mechanistically, GNAS KO-induced sensitization is independent of the canonical G-protein activities but unexpectedly mediated by viral mimicry-related interferon (IFN) responses, characterized by TBK1 and IRF3 activation, double-stranded RNA formation, and transposable element (TE) expression. GNAS KO additionally synergizes with HDAC3 inhibition to enhance CD8+ T cell-induced cytotoxicity. Moreover, we observe in human lymphoma patients that low GNAS expression is associated with high baseline TE expression and upregulated IFN signaling and shares common disrupted biological activities with GNAS KO in histone modification, mRNA processing, and transcriptional regulation. Collectively, our findings establish an unprecedented link between HDAC3 inhibition and viral mimicry in lymphoma. We suggest low GNAS expression as a potential biomarker that reflects viral mimicry priming for enhanced response to HDAC3 inhibition in the clinical treatment of lymphoma, especially the CREBBP wild-type cases.

PICKLE and HISTONE DEACETYLASE6 coordinately regulate genes and transposable elements in Arabidopsis

Chromatin dynamics play essential roles in transcriptional regulation. The chromodomain helicase DNA-binding domain 3 chromatin remodeler PICKLE (PKL) and HISTONE DEACETYLASE6 (HDA6) are required for transcriptional gene silencing, but their coordinated function in gene repression requires further study. Through a genetic suppressor screen, we found that a point mutation at PKL could partially restore the developmental defects of a weak Polycomb repressive complex 1 (PRC1) mutant (ring1a-2 ring1b-3), in which RING1A expression is suppressed by a T-DNA insertion at the promoter. Compared to ring1a-2 ring1b-3, the expression of RING1A is increased, nucleosome occupancy is reduced, and the histone 3 lysine 9 acetylation (H3K9ac) level is increased at the RING1A locus in the pkl ring1a-2 ring1b-3 triple mutant. HDA6 interacts with PKL and represses RING1A expression similarly to PKL genetically and molecularly in the ring1a-2 ring1b-3 background. Furthermore, we show that PKL and HDA6 suppress the expression of a set of genes and transposable elements (TEs) by increasing nucleosome density and reducing H3K9ac. Genome-wide analysis indicated they possibly coordinately maintain DNA methylation as well. Our findings suggest that PKL and HDA6 function together to reduce H3K9ac and increase nucleosome occupancy, thereby facilitating gene/TE regulation in Arabidopsis (Arabidopsis thaliana).

Mini-heterochromatin domains constrain the cis-regulatory impact of SVA transposons in human brain development and disease

SVA (SINE (short interspersed nuclear element)-VNTR (variable number of tandem repeats)-Alu) retrotransposons remain active in humans and contribute to individual genetic variation. Polymorphic SVA alleles harbor gene regulatory potential and can cause genetic disease. However, how SVA insertions are controlled and functionally impact human disease is unknown. Here we dissect the epigenetic regulation and influence of SVAs in cellular models of X-linked dystonia parkinsonism (XDP), a neurodegenerative disorder caused by an SVA insertion at the TAF1 locus. We demonstrate that the KRAB zinc finger protein ZNF91 establishes H3K9me3 and DNA methylation over SVAs, including polymorphic alleles, in human neural progenitor cells. The resulting mini-heterochromatin domains attenuate the cis-regulatory impact of SVAs. This is critical for XDP pathology; removal of local heterochromatin severely aggravates the XDP molecular phenotype, resulting in increased TAF1 intron retention and reduced expression. Our results provide unique mechanistic insights into how human polymorphic transposon insertions are recognized and how their regulatory impact is constrained by an innate epigenetic defense system.

3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data

RNA sequencing is the gold-standard method to quantify transcriptomic changes between two conditions. The overwhelming majority of data analysis methods available are focused on polyadenylated RNA transcribed from single-copy genes and overlook transcripts from repeated sequences such as transposable elements (TEs). These self-autonomous genetic elements are increasingly studied, and specialized tools designed to handle multimapping sequencing reads are available. Transfer RNAs are transcribed by RNA polymerase III and are essential for protein translation. There is a need for integrated software that is able to analyze multiple types of RNA. Here, we present 3t-seq, a Snakemake pipeline for integrated differential expression analysis of transcripts from single-copy genes, TEs, and tRNA. 3t-seq produces an accessible report and easy-to-use results for downstream analysis starting from raw sequencing data and performing quality control, genome mapping, gene expression quantification, and statistical testing. It implements three methods to quantify TEs expression and one for tRNA genes. It provides an easy-to-configure method to manage software dependencies that lets the user focus on results. 3t-seq is released under MIT license and is available at https://github.com/boulardlab/3t-seq.

Genetic-epigenetic interplay in the determination of plant 3D genome organization

The 3D chromatin organization plays a major role in the control of gene expression. However, our comprehension of the governing principles behind nuclear organization remains incomplete. Particularly, the spatial segregation of loci with similar repressive transcriptional states in plants poses a significant yet poorly understood puzzle. In this study, employing a combination of genetics and advanced 3D genomics approaches, we demonstrated that a redistribution of facultative heterochromatin marks in regions usually occupied by constitutive heterochromatin marks disrupts the 3D genome compartmentalisation. This disturbance, in turn, triggers novel chromatin interactions between genic and transposable element (TE) regions. Interestingly, our results imply that epigenetic features, constrained by genetic factors, intricately mold the landscape of 3D genome organisation. This study sheds light on the profound genetic-epigenetic interplay that underlies the regulation of gene expression within the intricate framework of the 3D genome. Our findings highlight the complexity of the relationships between genetic determinants and epigenetic features in shaping the dynamic configuration of the 3D genome.

Investigation of chimeric transcripts derived from LINE-1 and Alu retrotransposons in cerebellar tissues of individuals with autism spectrum disorder (ASD)

LINE-1 and Alu retrotransposons are components of the human genome and have been implicated in many human diseases. These elements can influence human transcriptome plasticity in various mechanisms. Chimeric transcripts derived from LINE-1 and Alu can also impact the human transcriptome, such as exonization and post-transcriptional modification. However, its specific role in ASD neuropathology remains unclear, particularly in the cerebellum tissues. We performed RNA-sequencing of post-mortem cerebellum tissues from ASD and unaffected individuals for transposable elements profiling and chimeric transcript identification. The majority of free transcripts of transposable elements were not changed in the cerebellum tissues of ASD compared with unaffected individuals. Nevertheless, we observed that chimeric transcripts derived from LINE-1 and Alu were embedded in the transcripts of differentially expressed genes in the cerebellum of ASD, and these genes were related to developments and abnormalities of the cerebellum. In addition, the expression levels of these genes were correlated with the significantly decreased thickness of the molecular layer in the cerebellum of ASD. We also found that global methylation and expression of LINE-1 and Alu elements were not changed in ASD, but observed in the ASD sub-phenotypes. Our findings showed associations between transposable elements and cerebellar abnormalities in ASD, particularly in distinct phenotypic subgroups. Further investigations using appropriate models are warranted to elucidate the structural and functional implications of LINE-1 and Alu elements in ASD neuropathology.

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programs and repeats in pluripotent cells

H3K9me3-heterochromatin, established by lysine methyltransferases (KMTs) and compacted by HP1 isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3-heterochromatin stability is presently limited to individual domains and DNA repeats. We engineered Suv39h2 KO mouse embryonic stem cells to degrade remaining two H3K9me3-KMTs within one hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A "binary switch" governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMTs' depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening, and exit from pluripotency within 12 hr. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3-heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

Transcripts of repetitive DNA elements signal to block phagocytosis of hematopoietic stem cells

Macrophages maintain hematopoietic stem cell (HSC) quality by assessing cell surface Calreticulin (Calr), an "eat-me" signal induced by reactive oxygen species (ROS). Using zebrafish genetics, we identified Beta-2-microglobulin (B2m) as a crucial "don't eat-me" signal on blood stem cells. A chemical screen revealed inducers of surface Calr that promoted HSC proliferation without triggering ROS or macrophage clearance. Whole-genome CRISPR-Cas9 screening showed that Toll-like receptor 3 (Tlr3) signaling regulated b2m expression. Targeting b2m or tlr3 reduced the HSC clonality. Elevated B2m levels correlated with high expression of repetitive element (RE) transcripts. Overall, our data suggest that RE-associated double-stranded RNA could interact with TLR3 to stimulate surface expression of B2m on hematopoietic stem and progenitor cells. These findings suggest that the balance of Calr and B2m regulates macrophage-HSC interactions and defines hematopoietic clonality.

A Quantitative Computational Framework for Allopolyploid Single-Cell Data Integration and Core Gene Ranking in Development

Polyploidization drives regulatory and phenotypic innovation. How the merger of different genomes contributes to polyploid development is a fundamental issue in evolutionary developmental biology and breeding research. Clarifying this issue is challenging because of genome complexity and the difficulty in tracking stochastic subgenome divergence during development. Recent single-cell sequencing techniques enabled probing subgenome-divergent regulation in the context of cellular differentiation. However, analyzing single-cell data suffers from high error rates due to high dimensionality, noise, and sparsity, and the errors stack up in polyploid analysis due to the increased dimensionality of comparisons between subgenomes of each cell, hindering deeper mechanistic understandings. In this study, we develop a quantitative computational framework, called "pseudo-genome divergence quantification" (pgDQ), for quantifying and tracking subgenome divergence directly at the cellular level. Further comparing with cellular differentiation trajectories derived from single-cell RNA sequencing data allows for an examination of the relationship between subgenome divergence and the progression of development. pgDQ produces robust results and is insensitive to data dropout and noise, avoiding high error rates due to multiple comparisons of genes, cells, and subgenomes. A statistical diagnostic approach is proposed to identify genes that are central to subgenome divergence during development, which facilitates the integration of different data modalities, enabling the identification of factors and pathways that mediate subgenome-divergent activity during development. Case studies have demonstrated that applying pgDQ to single-cell and bulk tissue transcriptomic data promotes a systematic and deeper understanding of how dynamic subgenome divergence contributes to developmental trajectories in polyploid evolution.

Pathogen infection alters the gene expression landscape of transposable elements in Drosophila melanogaster

Transposable elements make up substantial proportions of eukaryotic genomes and many are thought to be remnants of ancient viral infections. Current research has begun to highlight the role transposable elements can play in the immune system response to infections. However, most of our knowledge about transposable element expression during infection is limited by the specific host and pathogen factors from each study, making it difficult to compare studies and develop broader patterns regarding the role of transposable elements during infection. Here, we use the tools and resources available in the model, Drosophila melanogaster, to analyze multiple gene expression datasets of flies subject to bacterial, fungal, and viral infections. We analyzed differences in pathogen species, host genotype, host tissue, and sex to understand how these factors impact transposable element expression during infection. Our results highlight both shared and unique transposable element expression patterns between pathogens and suggest a larger effect of pathogen factors over host factors for influencing transposable element expression.

Human cytomegalovirus harnesses host L1 retrotransposon for efficient replication

Genetic parasites, including viruses and transposons, exploit components from the host for their own replication. However, little is known about virus-transposon interactions within host cells. Here, we discover a strategy where human cytomegalovirus (HCMV) hijacks L1 retrotransposon encoded protein during its replication cycle. HCMV infection upregulates L1 expression by enhancing both the expression of L1-activating transcription factors, YY1 and RUNX3, and the chromatin accessibility of L1 promoter regions. Increased L1 expression, in turn, promotes HCMV replicative fitness. Affinity proteomics reveals UL44, HCMV DNA polymerase subunit, as the most abundant viral binding protein of the L1 ribonucleoprotein (RNP) complex. UL44 directly interacts with L1 ORF2p, inducing DNA damage responses in replicating HCMV compartments. While increased L1-induced mutagenesis is not observed in HCMV for genetic adaptation, the interplay between UL44 and ORF2p accelerates viral DNA replication by alleviating replication stress. Our findings shed light on how HCMV exploits host retrotransposons for enhanced viral fitness.

DNA methylation governs the sensitivity of repeats to restriction by the HUSH-MORC2 corepressor

The human silencing hub (HUSH) complex binds to transcripts of LINE-1 retrotransposons (L1s) and other genomic repeats, recruiting MORC2 and other effectors to remodel chromatin. How HUSH and MORC2 operate alongside DNA methylation, a central epigenetic regulator of repeat transcription, remains largely unknown. Here we interrogate this relationship in human neural progenitor cells (hNPCs), a somatic model of brain development that tolerates removal of DNA methyltransferase DNMT1. Upon loss of MORC2 or HUSH subunit TASOR in hNPCs, L1s remain silenced by robust promoter methylation. However, genome demethylation and activation of evolutionarily-young L1s attracts MORC2 binding, and simultaneous depletion of DNMT1 and MORC2 causes massive accumulation of L1 transcripts. We identify the same mechanistic hierarchy at pericentromeric α-satellites and clustered protocadherin genes, repetitive elements important for chromosome structure and neurodevelopment respectively. Our data delineate the epigenetic control of repeats in somatic cells, with implications for understanding the vital functions of HUSH-MORC2 in hypomethylated contexts throughout human development.

Combinatorial quantification of 5mC and 5hmC at individual CpG dyads and the transcriptome in single cells reveals modulators of DNA methylation maintenance fidelity

Inheritance of 5-methylcytosine from one cell generation to the next by DNA methyltransferase 1 (DNMT1) plays a key role in regulating cellular identity. While recent work has shown that the activity of DNMT1 is imprecise, it remains unclear how the fidelity of DNMT1 is tuned in different genomic and cell state contexts. Here we describe Dyad-seq, a method to quantify the genome-wide methylation status of cytosines at the resolution of individual CpG dinucleotides to find that the fidelity of DNMT1-mediated maintenance methylation is related to the local density of DNA methylation and the landscape of histone modifications. To gain deeper insights into methylation/demethylation turnover dynamics, we first extended Dyad-seq to quantify all combinations of 5-methylcytosine and 5-hydroxymethylcytosine at individual CpG dyads. Next, to understand how cell state transitions impact maintenance methylation, we scaled the method down to jointly profile genome-wide methylation levels, maintenance methylation fidelity and the transcriptome from single cells (scDyad&T-seq). Using scDyad&T-seq, we demonstrate that, while distinct cell states can substantially impact the activity of the maintenance methylation machinery, locally there exists an intrinsic relationship between DNA methylation density, histone modifications and DNMT1-mediated maintenance methylation fidelity that is independent of cell state.

PHF2-mediated H3K9me balance orchestrates heterochromatin stability and neural progenitor proliferation

Heterochromatin stability is crucial for progenitor proliferation during early neurogenesis. It relays on the maintenance of local hubs of H3K9me. However, understanding the formation of efficient localized levels of H3K9me remains limited. To address this question, we used neural stem cells to analyze the function of the H3K9me2 demethylase PHF2, which is crucial for progenitor proliferation. Through mass-spectroscopy and genome-wide assays, we show that PHF2 interacts with heterochromatin components and is enriched at pericentromeric heterochromatin (PcH) boundaries where it maintains transcriptional activity. This binding is essential for silencing the satellite repeats, preventing DNA damage and genome instability. PHF2's depletion increases the transcription of heterochromatic repeats, accompanied by a decrease in H3K9me3 levels and alterations in PcH organization. We further show that PHF2's PHD and catalytic domains are crucial for maintaining PcH stability, thereby safeguarding genome integrity. These results highlight the multifaceted nature of PHF2's functions in maintaining heterochromatin stability and regulating gene expression during neural development. Our study unravels the intricate relationship between heterochromatin stability and progenitor proliferation during mammalian neurogenesis.

ΔNp63-restricted viral mimicry response impedes cancer cell viability and remodels tumor microenvironment in esophageal squamous cell carcinoma

Tumor protein p63 isoform ΔNp63 plays roles in the squamous epithelium and squamous cell carcinomas (SCCs), including esophageal SCC (ESCC). By integrating data from cell lines and our latest patient-derived organoid cultures, derived xenograft models, and clinical sample transcriptomic analyses, we identified a novel and robust oncogenic role of ΔNp63 in ESCC. We showed that ΔNp63 maintains the repression of cancer cell endogenous retrotransposon expression and cellular double-stranded RNA sensing. These subsequently lead to a restricted cancer cell viral mimicry response and suppressed induction of tumor-suppressive type I interferon (IFN-I) signaling through the regulations of Signal transducer and activator of transcription 1, Interferon regulatory factor 1, and cGAS-STING pathway. The cancer cell ΔNp63/IFN-I signaling axis affects both the cancer cell and tumor-infiltrating immune cell (TIIC) compartments. In cancer cells, depletion of ΔNp63 resulted in reduced cell viability. ΔNp63 expression is negatively associated with the anticancer responses to viral mimicry booster treatments targeting cancer cells. In the tumor microenvironment, cancer cell TP63 expression negatively correlates with multiple TIIC signatures in ESCC clinical samples. ΔNp63 depletion leads to increased cancer cell antigen presentation molecule expression and enhanced recruitment and reprogramming of tumor-infiltrating myeloid cells. Similar IFN-I signaling and TIIC signature association with ΔNp63 were also observed in lung SCC. These results support the potential application of ΔNp63 as a therapeutic target and a biomarker to guide candidate anticancer treatments exploring viral mimicry responses.

Heat stress induced piRNA alterations in pachytene spermatocytes and round spermatids

BACKGROUND

Spermatogenesis is a temperature-sensitive process, and elevation in temperature hampers this process quickly and significantly. We studied the molecular effects of testicular heating on piRNAs and gene expression in rat testicular germ cells.

METHODS

We generated a cryptorchid rat model by displacing the testis from the scrotal sac (34 °C) to the abdominal area (37 °C) and sacrificed animals after 1 day, 3 days, and 5 days. Pachytene spermatocytes and round spermatids were purified using elutriation centrifugation and percoll gradient methods. We performed transcriptome sequencing in pachytene spermatocytes and round spermatids to identify differentially expressed piRNAs and their probable targets, i.e., TE transcripts and mRNAs.

RESULTS

As a result of heat stress, we observed significant upregulation of piRNAs and TE transcripts in testicular germ cells. In addition to this, piRNA biogenesis machinery and heat shock proteins (Hsp70 and Hsp90 family members) were upregulated. mRNAs have also been proposed as targets for piRNAs; therefore, we shortlisted certain piRNA-mRNA pairs with an inverse relationship of expression. We observed that in testicular heat stress, the heat shock proteins go hand-in-hand with the upregulation of piRNA biogenesis machinery. The dysregulation of piRNAs in heat-stressed germ cells, increased ping-pong activity, and disturbed expression of piRNA target transcripts suggest a connection between piRNAs, mRNAs, and TE transcripts.

CONCLUSIONS

In heat stress, piRNAs, piRNA machinery, and heat shock proteins are activated to deal with low levels of stress, which is followed by a rescue approach in prolonged stressaccompained by high TE activity to allow genetic mutations, perhaps for survival and adaptability.

Accurate allocation of multimapped reads enables regulatory element analysis at repeats

Transposable elements (TEs) and other repetitive regions have been shown to contain gene regulatory elements, including transcription factor binding sites. However, regulatory elements harbored by repeats have proven difficult to characterize using short-read sequencing assays such as ChIP-seq or ATAC-seq. Most regulatory genomics analysis pipelines discard "multimapped" reads that align equally well to multiple genomic locations. Because multimapped reads arise predominantly from repeats, current analysis pipelines fail to detect a substantial portion of regulatory events that occur in repetitive regions. To address this shortcoming, we developed Allo, a new approach to allocate multimapped reads in an efficient, accurate, and user-friendly manner. Allo combines probabilistic mapping of multimapped reads with a convolutional neural network that recognizes the read distribution features of potential peaks, offering enhanced accuracy in multimapping read assignment. Allo also provides read-level output in the form of a corrected alignment file, making it compatible with existing regulatory genomics analysis pipelines and downstream peak-finders. In a demonstration application on CTCF ChIP-seq data, we show that Allo results in the discovery of thousands of new CTCF peaks. Many of these peaks contain the expected cognate motif and/or serve as TAD boundaries. We additionally apply Allo to a diverse collection of ENCODE ChIP-seq data sets, resulting in multiple previously unidentified interactions between transcription factors and repetitive element families. Finally, we show that Allo may be particularly beneficial in identifying ChIP-seq peaks at centromeres, near segmentally duplicated genes, and in younger TEs, enabling new regulatory analyses in these regions.

ALS molecular subtypes are a combination of cellular, genetic, and pathological features learned by deep multiomics classifiers

Amyotrophic Lateral Sclerosis (ALS) is a complex syndrome with multiple genetic causes and wide variation in disease presentation. Despite this general heterogeneity, several common factors have been identified. For example, nearly all patients show pathological accumulations of phosphorylated TDP-43 protein in affected regions of the motor cortex and spinal cord. Moreover, large patient cohort studies have revealed that most patient samples can be grouped into a small number of ALS subtypes, as defined by their transcriptomic profiles. These ALS molecular subtypes can be grouped by whether postmortem motor cortex samples display signatures of: mitochondrial dysfunction and oxidative stress (ALS-Ox), microglial activation and neuroinflammation (ALS-Glia), or dense TDP-43 pathology and associated transposable element de-silencing (ALS-TE). In this study, we have built a deep layer ALS neural network classifier (DANcer) that has learned to accurately assign patient samples to these ALS subtypes, and which can be run on either bulk or single-cell datasets. Upon applying this classifier to an expanded ALS patient cohort from the NYGC ALS Consortium, we show that ALS Molecular Subtypes are robust across clinical centers, with no new subtypes appearing in a cohort that has quadrupled in size. Signatures from two of these molecular subtypes strongly correlate with disease duration: ALS-TE signatures in cortex and ALS-Glia signatures in spinal cord, revealing molecular correlates of clinical features. Finally, we use single nucleus RNA sequencing to reveal the cell type-specific contributions to ALS subtype, as determined by our single-cell classifier (scDANCer). Single-cell transcriptomes reveal that ALS molecular subtypes are recapitulated in neurons and glia, with both ALS-wide shared alterations in each cell type as well as ALS subtype-specific alterations. In summary, ALS molecular subtypes: (1) are robust across large cohorts of sporadic and familial ALS patient samples, (2) represent a combination of cellular, genetic, and pathological features, and (3) correlate with clinical features of ALS.

Abstract Figure

Endogenous retroviruses are dysregulated in ALS

Amyotrophic lateral sclerosis (ALS) is a universally fatal neurodegenerative disease with no cure. Human endogenous retroviruses (HERVs) have been implicated in its pathogenesis but their relevance to ALS is not fully understood. We examined bulk RNA-seq data from almost 2,000 ALS and unaffected control samples derived from the cortex and spinal cord. Using different methods of feature selection, including differential expression analysis and machine learning, we discovered that transcription of HERV-K loci 1q22 and 8p23.1 were significantly upregulated in the spinal cord of individuals with ALS. Additionally, we identified a subset of ALS patients with upregulated HERV-K expression in the cortex and spinal cord. We also found the expression of HERV-K loci 19q11 and 8p23.1 was correlated with protein coding genes previously implicated in ALS and dysregulated in ALS patients in this study. These results clarify the association of HERV-K and ALS and highlight specific genes in the pathobiology of late-stage ALS.

ChAHP2 and ChAHP control diverse retrotransposons by complementary activities

Retrotransposon control in mammals is an intricate process that is effectuated by a broad network of chromatin regulatory pathways. We previously discovered ChAHP, a protein complex with repressive activity against short interspersed element (SINE) retrotransposons that is composed of the transcription factor ADNP, chromatin remodeler CHD4, and HP1 proteins. Here we identify ChAHP2, a protein complex homologous to ChAHP, in which ADNP is replaced by ADNP2. ChAHP2 is predominantly targeted to endogenous retroviruses (ERVs) and long interspersed elements (LINEs) via HP1β-mediated binding of H3K9 trimethylated histones. We further demonstrate that ChAHP also binds these elements in a manner mechanistically equivalent to that of ChAHP2 and distinct from DNA sequence-specific recruitment at SINEs. Genetic ablation of ADNP2 alleviates ERV and LINE1 repression, which is synthetically exacerbated by additional depletion of ADNP. Together, our results reveal that the ChAHP and ChAHP2 complexes function to control both nonautonomous and autonomous retrotransposons by complementary activities, further adding to the complexity of mammalian transposon control.

Endogenous retroviruses mediate transcriptional rewiring in response to oncogenic signaling in colorectal cancer

Cancer cells exhibit rewired transcriptional regulatory networks that promote tumor growth and survival. However, the mechanisms underlying the formation of these pathological networks remain poorly understood. Through a pan-cancer epigenomic analysis, we found that primate-specific endogenous retroviruses (ERVs) are a rich source of enhancers displaying cancer-specific activity. In colorectal cancer and other epithelial tumors, oncogenic MAPK/AP1 signaling drives the activation of enhancers derived from the primate-specific ERV family LTR10. Functional studies in colorectal cancer cells revealed that LTR10 elements regulate tumor-specific expression of multiple genes associated with tumorigenesis, such as ATG12 and XRCC4. Within the human population, individual LTR10 elements exhibit germline and somatic structural variation resulting from a highly mutable internal tandem repeat region, which affects AP1 binding activity. Our findings reveal that ERV-derived enhancers contribute to transcriptional dysregulation in response to oncogenic signaling and shape the evolution of cancer-specific regulatory networks.

Manipulating mitochondrial reactive oxygen species alters survival in unexpected ways in a Drosophila Cdk5 model of neurodegeneration

Reactive oxygen species (ROS) are associated with aging and neurodegeneration, but the significance of this association remains obscure. Here, using a Drosophila Cdk5 model of age-related neurodegeneration, we probe this relationship in the pathologically relevant tissue, the brain, by quantifying three specific mitochondrial ROS and manipulating these redox species pharmacologically. Our goal is to ask whether pathology-associated changes in redox state are detrimental for survival, whether they may be beneficial responses to pathology, or whether they are covariates of pathology that do not alter viability. We find, surprisingly, that increasing mitochondrial H2O2 correlates with improved survival. We also find evidence that drugs that alter the mitochondrial glutathione redox potential modulate survival primarily through the compensatory effects they induce rather than through their direct effects on the final mitochondrial glutathione redox potential. We also find that the response to treatment with a redox-altering drug varies depending on the age and genotype of the individual receiving the drug as well as the duration of the treatment. These data have important implications for the design and interpretation of studies investigating the effect of redox state on health and disease as well as on efforts to modify the redox state to achieve therapeutic goals.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

Mutant IDH1 inhibition induces dsDNA sensing to activate tumor immunity

Isocitrate dehydrogenase 1 (IDH1) is the most commonly mutated metabolic gene across human cancers. Mutant IDH1 (mIDH1) generates the oncometabolite (R)-2-hydroxyglutarate, disrupting enzymes involved in epigenetics and other processes. A hallmark of IDH1-mutant solid tumors is T cell exclusion, whereas mIDH1 inhibition in preclinical models restores antitumor immunity. Here, we define a cell-autonomous mechanism of mIDH1-driven immune evasion. IDH1-mutant solid tumors show selective hypermethylation and silencing of the cytoplasmic double-stranded DNA (dsDNA) sensor CGAS, compromising innate immune signaling. mIDH1 inhibition restores DNA demethylation, derepressing CGAS and transposable element (TE) subclasses. dsDNA produced by TE-reverse transcriptase (TE-RT) activates cGAS, triggering viral mimicry and stimulating antitumor immunity. In summary, we demonstrate that mIDH1 epigenetically suppresses innate immunity and link endogenous RT activity to the mechanism of action of a US Food and Drug Administration-approved oncology drug.

The nuclear periphery confers repression on H3K9me2-marked genes and transposons to shape cell fate

Heterochromatic loci marked by histone H3 lysine 9 dimethylation (H3K9me2) are enriched at the nuclear periphery in metazoans, but the effect of spatial position on heterochromatin function has not been defined. Here, we remove three nuclear lamins and lamin B receptor (LBR) in mouse embryonic stem cells (mESCs) and show that heterochromatin detaches from the nuclear periphery. Mutant mESCs sustain naïve pluripotency and maintain H3K9me2 across the genome but cannot repress H3K9me2-marked genes or transposons. Further, mutant cells fail to differentiate into epiblast-like cells (EpiLCs), a transition that requires the expansion of H3K9me2 across the genome. Mutant EpiLCs can silence naïve pluripotency genes and activate epiblast-stage genes. However, H3K9me2 cannot repress markers of alternative fates, including primitive endoderm. We conclude that the nuclear periphery controls the spatial position, dynamic remodeling, and repressive capacity of H3K9me2-marked heterochromatin to shape cell fate decisions.

Inheritance of H3K9 methylation regulates genome architecture in Drosophila early embryos

Constitutive heterochromatin is essential for transcriptional silencing and genome integrity. The establishment of constitutive heterochromatin in early embryos and its role in early fruitfly development are unknown. Lysine 9 trimethylation of histone H3 (H3K9me3) and recruitment of its epigenetic reader, heterochromatin protein 1a (HP1a), are hallmarks of constitutive heterochromatin. Here, we show that H3K9me3 is transmitted from the maternal germline to the next generation. Maternally inherited H3K9me3, and the histone methyltransferases (HMT) depositing it, are required for the organization of constitutive heterochromatin: early embryos lacking H3K9 methylation display de-condensation of pericentromeric regions, centromere-centromere de-clustering, mitotic defects, and nuclear shape irregularities, resulting in embryo lethality. Unexpectedly, quantitative CUT&Tag and 4D microscopy measurements of HP1a coupled with biophysical modeling revealed that H3K9me2/3 is largely dispensable for HP1a recruitment. Instead, the main function of H3K9me2/3 at this developmental stage is to drive HP1a clustering and subsequent heterochromatin compaction. Our results show that HP1a binding to constitutive heterochromatin in the absence of H3K9me2/3 is not sufficient to promote proper embryo development and heterochromatin formation. The loss of H3K9 HMTs and H3K9 methylation alters genome organization and hinders embryonic development.

Technology to the rescue: how to uncover the role of transposable elements in preimplantation development

Transposable elements (TEs) are highly expressed in preimplantation development. Preimplantation development is the phase when the cells of the early embryo undergo the first cell fate choice and change from being totipotent to pluripotent. A range of studies have advanced our understanding of TEs in preimplantation, as well as their epigenetic regulation and functional roles. However, many questions remain about the implications of TE expression during early development. Challenges originate first due to the abundance of TEs in the genome, and second because of the limited cell numbers in preimplantation. Here we review the most recent technological advancements promising to shed light onto the role of TEs in preimplantation development. We explore novel avenues to identify genomic TE insertions and improve our understanding of the regulatory mechanisms and roles of TEs and their RNA and protein products during early development.

Combining EHMT and PARP Inhibition: A Strategy to Diminish Therapy-Resistant Ovarian Cancer Tumor Growth while Stimulating Immune Activation

Despite the success of poly-ADP-ribose polymerase inhibitors (PARPi) in the clinic, high rates of resistance to PARPi presents a challenge in the treatment of ovarian cancer, thus it is imperative to find therapeutic strategies to combat PARPi resistance. Here, we demonstrate that inhibition of epigenetic modifiers euchromatic histone lysine methyltransferases 1/2 (EHMT1/2) reduces the growth of multiple PARPi-resistant ovarian cancer cell lines and tumor growth in a PARPi-resistant mouse model of ovarian cancer. We found that combinatory EHMT and PARP inhibition increases immunostimulatory double-stranded RNA formation and elicits several immune signaling pathways in vitro. Using epigenomic profiling and transcriptomics, we found that EHMT2 is bound to transposable elements, and that EHMT inhibition leads to genome-wide epigenetic and transcriptional derepression of transposable elements. We validated EHMT-mediated activation of immune signaling and upregulation of transposable element transcripts in patient-derived, therapy-naïve, primary ovarian tumors, suggesting potential efficacy in PARPi-sensitive disease as well. Importantly, using multispectral immunohistochemistry, we discovered that combinatory therapy increased CD8 T-cell activity in the tumor microenvironment of the same patient-derived tissues. In a PARPi-resistant syngeneic murine model, EHMT and PARP inhibition combination inhibited tumor progression and increased Granzyme B+ cells in the tumor. Together, our results provide evidence that combinatory EHMT and PARP inhibition stimulates a cell autologous immune response in vitro, is an effective therapy to reduce PARPi-resistant ovarian tumor growth in vivo, and promotes antitumor immunity activity in the tumor microenvironment of patient-derived ex vivo tissues of ovarian cancer.

iPS cell generation-associated point mutations include many C > T substitutions via different cytosine modification mechanisms

Genomic aberrations are a critical impediment for the safe medical use of iPSCs and their origin and developmental mechanisms remain unknown. Here we find through WGS analysis of human and mouse iPSC lines that genomic mutations are de novo events and that, in addition to unmodified cytosine base prone to deamination, the DNA methylation sequence CpG represents a significant mutation-prone site. CGI and TSS regions show increased mutations in iPSCs and elevated mutations are observed in retrotransposons, especially in the AluY subfamily. Furthermore, increased cytosine to thymine mutations are observed in differentially methylated regions. These results indicate that in addition to deamination of cytosine, demethylation of methylated cytosine, which plays a central role in genome reprogramming, may act mutagenically during iPSC generation.

Genetic and developmental divergence in the neural crest programme between cichlid fish species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.

Walsh, Marazzi I, Kappei D, Guccione E, Jeyasekharan AD. A PRMT5-ZNF326 axis mediates innate immune activation upon replication stress

DNA replication stress (RS) is a widespread phenomenon in carcinogenesis, causing genomic instability and extensive chromatin alterations. DNA damage leads to activation of innate immune signaling, but little is known about transcriptional regulators mediating such signaling upon RS. Using a chemical screen, we identified protein arginine methyltransferase 5 (PRMT5) as a key mediator of RS-dependent induction of interferon-stimulated genes (ISGs). This response is also associated with reactivation of endogenous retroviruses (ERVs). Using quantitative mass spectrometry, we identify proteins with PRMT5-dependent symmetric dimethylarginine (SDMA) modification induced upon RS. Among these, we show that PRMT5 targets and modulates the activity of ZNF326, a zinc finger protein essential for ISG response. Our data demonstrate a role for PRMT5-mediated SDMA in the context of RS-induced transcriptional induction, affecting physiological homeostasis and cancer therapy.

Identification and characterization of a human MORC2 DNA binding region that is required for gene silencing

The eukaryotic microrchidia (MORC) protein family are DNA gyrase, Hsp90, histidine kinase, MutL (GHKL)-type ATPases involved in gene expression regulation and chromatin compaction. The molecular mechanisms underlying these activities are incompletely understood. Here we studied the full-length human MORC2 protein biochemically. We identified a DNA binding site in the C-terminus of the protein, and we observe that this region is heavily phosphorylated in cells. Phosphorylation of MORC2 reduces its affinity for DNA and appears to exclude the protein from the nucleus. We observe that DNA binding by MORC2 reduces its ATPase activity and that MORC2 can topologically entrap multiple DNA substrates between its N-terminal GHKL and C-terminal coiled coil 3 dimerization domains. Finally, we observe that the MORC2 C-terminal DNA binding region is required for gene silencing in cells. Together, our data provide a model to understand how MORC2 engages with DNA substrates to mediate gene silencing.

Overcoming Clinical Resistance to EZH2 Inhibition Using Rational Epigenetic Combination Therapy

Epigenetic dependencies have become evident in many cancers. On the basis of antagonism between BAF/SWI-SNF and PRC2 in SMARCB1-deficient sarcomas, we recently completed the clinical trial of the EZH2 inhibitor tazemetostat. However, the principles of tumor response to epigenetic therapy in general, and tazemetostat in particular, remain unknown. Using functional genomics and diverse experimental models, we define molecular mechanisms of tazemetostat resistance in SMARCB1-deficient tumors. We found distinct acquired mutations that converge on the RB1/E2F axis and decouple EZH2-dependent differentiation and cell-cycle control. This allows tumor cells to escape tazemetostat-induced G1 arrest, suggests a general mechanism for effective therapy, and provides prospective biomarkers for therapy stratification, including PRICKLE1. On the basis of this, we develop a combination strategy to circumvent tazemetostat resistance using bypass targeting of AURKB. This offers a paradigm for rational epigenetic combination therapy suitable for translation to clinical trials for epithelioid sarcomas, rhabdoid tumors, and other epigenetically dysregulated cancers.

SIGNIFICANCE

Genomic studies of patient epithelioid sarcomas and rhabdoid tumors identify mutations converging on a common pathway for response to EZH2 inhibition. Resistance mutations decouple drug-induced differentiation from cell-cycle control. We identify an epigenetic combination strategy to overcome resistance and improve durability of response, supporting its investigation in clinical trials. See related commentary by Paolini and Souroullas, p. 903. This article is featured in Selected Articles from This Issue, p. 897.