SQuIRE reveals locus-specific regulation of interspersed repeat expression

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.

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.

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.

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

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.

Epigenetic therapy potentiates transposable element transcription to create tumor-enriched antigens in glioblastoma cells

Inhibiting epigenetic modulators can transcriptionally reactivate transposable elements (TEs). These TE transcripts often generate unique peptides that can serve as immunogenic antigens for immunotherapy. Here, we ask whether TEs activated by epigenetic therapy could appreciably increase the antigen repertoire in glioblastoma, an aggressive brain cancer with low mutation and neoantigen burden. We treated patient-derived primary glioblastoma stem cell lines, an astrocyte cell line and primary fibroblast cell lines with epigenetic drugs, and identified treatment-induced, TE-derived transcripts that are preferentially expressed in cancer cells. We verified that these transcripts could produce human leukocyte antigen class I-presented antigens using liquid chromatography with tandem mass spectrometry pulldown experiments. Importantly, many TEs were also transcribed, even in proliferating nontumor cell lines, after epigenetic therapy, which suggests that targeted strategies like CRISPR-mediated activation could minimize potential side effects of activating unwanted genomic regions. The results highlight both the need for caution and the promise of future translational efforts in harnessing treatment-induced TE-derived antigens for targeted immunotherapy.

LINE-1 RNA triggers matrix formation in bone cells via a PKR-mediated inflammatory response

Transposable elements (TEs) are mobile genetic modules of viral derivation that have been co-opted to become modulators of mammalian gene expression. TEs are a major source of endogenous dsRNAs, signaling molecules able to coordinate inflammatory responses in various physiological processes. Here, we provide evidence for a positive involvement of TEs in inflammation-driven bone repair and mineralization. In newly fractured mice bone, we observed an early transient upregulation of repeats occurring concurrently with the initiation of the inflammatory stage. In human bone biopsies, analysis revealed a significant correlation between repeats expression, mechanical stress and bone mineral density. We investigated a potential link between LINE-1 (L1) expression and bone mineralization by delivering a synthetic L1 RNA to osteoporotic patient-derived mesenchymal stem cells and observed a dsRNA-triggered protein kinase (PKR)-mediated stress response that led to strongly increased mineralization. This response was associated with a strong and transient inflammation, accompanied by a global translation attenuation induced by eIF2α phosphorylation. We demonstrated that L1 transfection reshaped the secretory profile of osteoblasts, triggering a paracrine activity that stimulated the mineralization of recipient cells.

Dissecting transposable elements and endogenous retroviruses upregulation by HDAC inhibitors in leiomyosarcoma cells: Implications for the interferon response

Transposable elements (TEs) are of interest as immunomodulators for cancer therapies. TEs can fold into dsRNAs that trigger the interferon response. Here, we investigated the effect of different HDAC inhibitors (HDACIs) on the expression of TEs in leiomyosarcoma cells. Our data show that endogenous retroviruses (ERVs), especially ERV1 elements, are upregulated after treatment with HDAC1/2/3-specific inhibitors. Surprisingly, the interferon response was not activated. We observed an increase in A-to-I editing of upregulated ERV1. This could have an impact on the stability of dsRNAs and the activation of the interferon response. We also found that H3K27ac levels are increased in the LTR12 subfamilies, which could be regulatory elements controlling the expression of proapoptotic genes such as TNFRSF10B. In summary, we provide a detailed characterization of TEs modulation in response to HDACIs and suggest the use of HDACIs in combination with ADAR inhibitors to induce cell death and support immunotherapy in cancer.

H3K9me3 loss and ERVs activation as hallmarks for osteoarthritis progression and knee joint aging

OBJECTIVE

This study aims to link aberrant endogenous retroviruses (ERVs) activation and osteoarthritis (OA) progression by comparing the chromatin accessibility and transcriptomic landscapes of diseased or intact joint tissues of OA patients.

METHOD

We performed ERVs-centric analysis on published ATAC-seq and RNA-seq data from OA patients' cartilage tissues. Here, we compared the outer region of the lateral tibial plateau, representing intact cartilage, to the inner region of the medial tibial plateau, representing damaged cartilage. In addition, cartilage tissue sections from OA patients and post-traumatic OA mouse models were assayed for global H3K9me3 abundance through immunohistochemistry staining.

RESULTS

Chromatin accessibility and transcription of ERVs, particularly from evolutionarily "intermediate age" ERVs families (ERV1 and ERVL), were enriched and elevated in OA cartilage. This integrative analysis suggests that H3K9me3-related heterochromatin loss might be mechanistically connected to ERV activation in OA tissue. We further verified that global H3K9me3 levels were reduced in diseased cartilage relative to intact tissue in OA patients and injury-induced OA mice.

CONCLUSION

The findings suggest a compelling hypothesis that the loss of H3K9me3, either due to aging or cellular stressors, may lead to ERVs reactivation that contributes to tissue inflammation and OA progression. This study unveils the intricate relationship between epigenetic alterations, ERVs activation, and OA, paving the way for potential therapeutic interventions targeting these pathogenic mechanisms.

The transcription factor PAX5 activates human LINE1 retrotransposons to induce cellular senescence

As a hallmark of senescent cells, the derepression of Long Interspersed Elements 1 (LINE1) transcription results in accumulated LINE1 cDNA, which triggers the secretion of the senescence-associated secretory phenotype (SASP) and paracrine senescence in a cGAS-STING pathway-dependent manner. However, transcription factors that govern senescence-associated LINE1 reactivation remain ill-defined. Here, we predict several transcription factors that bind to human LINE1 elements to regulate their transcription by analyzing the conserved binding motifs in the 5'-untranslated regions (UTR) of the commonly upregulated LINE1 elements in different types of senescent cells. Further analysis reveals that PAX5 directly binds to LINE1 5'-UTR and the binding is enhanced in senescent cells. The enrichment of PAX5 at the 5'-UTR promotes cellular senescence and SASP by activating LINE1. We also demonstrate that the longevity gene SIRT6 suppresses PAX5 transcription by directly binding to the PAX5 promoter, and overexpressing PAX5 abrogates the suppressive effect of SIRT6 on stress-dependent cellular senescence. Our work suggests that PAX5 could serve as a potential target for drug development aiming to suppress LINE1 activation and treat senescence-associated diseases.

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.

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.

Transcriptome analysis of cynomolgus macaques throughout their lifespan reveals age-related immune patterns

Despite the different perspectives by diverse research sectors spanning several decades, aging research remains uncharted territory for human beings. Therefore, we investigated the transcriptomic characteristics of eight male healthy cynomolgus macaques, and the annual sampling was designed with two individuals in four age groups. As a laboratory animal, the macaques were meticulously shielded from all environmental factors except aging. The results showed recent findings of certain immune response and the age-associated network of primate immunity. Three important aging patterns were identified and each gene clusters represented a different immune response. The increased expression pattern was predominantly associated with innate immune cells, such as Neutrophils and NK cells, causing chronic inflammation with aging whereas the other two decreased patterns were associated with adaptive immunity, especially "B cell activation" affecting antibody diversity of aging. Furthermore, the hub gene network of the patterns reflected transcriptomic age and correlated with human illness status, aiding in future human disease prediction. Our macaque transcriptome profiling results offer systematic insights into the age-related immunological features of primates.

An integrative epigenome-based strategy for unbiased functional profiling of clinical kinase inhibitors

More than 500 kinases are implicated in the control of most cellular process in mammals, and deregulation of their activity is linked to cancer and inflammatory disorders. 80 clinical kinase inhibitors (CKIs) have been approved for clinical use and hundreds are in various stages of development. However, CKIs inhibit other kinases in addition to the intended target(s), causing both enhanced clinical effects and undesired side effects that are only partially predictable based on in vitro selectivity profiling. Here, we report an integrative approach grounded on the use of chromatin modifications as unbiased, information-rich readouts of the functional effects of CKIs on macrophage activation. This approach exceeded the performance of transcriptome-based approaches and allowed us to identify similarities and differences among CKIs with identical intended targets, to recognize novel CKI specificities and to pinpoint CKIs that may be repurposed to control inflammation, thus supporting the utility of this strategy to improve selection and use of CKIs in clinical settings.

Ocrelizumab B cell depletion has no effect on HERV RNA expression in PBMC in MS patients

BACKGROUND

Epstein barr virus (EBV) infection of B cells is now understood to be one of the triggering events for the development of Multiple Sclerosis (MS), a progressive immune-mediated disease of the central nervous system. EBV infection is also linked to expression of human endogenous retroviruses (HERVs) of the HERV-W group, a further risk factor for the development of MS. Ocrelizumab is a high-potency disease-modifying treatment (DMT) for MS, which depletes B cells by targeting CD20.

OBJECTIVES

We studied the effects of ocrelizumab on gene expression in peripheral blood mononuclear cells (PBMC) from paired samples from 20 patients taken prior to and 6 months after beginning ocrelizumab therapy. We hypothesised that EBV and HERV-W loads would be lower in post-treatment samples.

METHODS

Samples were collected in Paxgene tubes, subject to RNA extraction and Illumina paired end short read mRNA sequencing with mapping of sequence reads to the human genome using Salmon and differential gene expression compared with DeSeq2. Mapping was also performed separately to the HERV-D database of HERV sequences and the EBV reference sequence.

RESULTS

Patient samples were more strongly clustered by individual rather than disease type (relapsing/remitting or primary progressive), treatment (pre and post), age, or sex. Fourteen genes, all clearly linked to B cell function were significantly down regulated in the post treatment samples. Interestingly only one pre-treatment sample had detectable EBV RNA and there were no significant differences in HERV expression (of any group) between pre- and post-treatment samples.

CONCLUSIONS

While EBV and HERV expression are clearly linked to triggering MS pathogenesis, it does not appear that high level expression of these viruses is a part of the ongoing disease process or that changes in virus load are associated with ocrelizumab treatment.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks. We additionally observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Finally, we demonstrated that loss of ATRX in a mesenchymal malignancy, undifferentiated pleomorphic sarcoma, results in similar epigenetic disruption and de-repression of transposable elements. Together, our results reveal a role for ATRX in maintaining epigenetic states and transcriptional repression in mesenchymal progenitors and tumor cells and in preventing aberrant differentiation in the progenitor context.

Pivotal role for long noncoding RNAs in zygotic genome activation in mice

Vertebrate life begins with fertilization, and then the zygote genome is activated after transient silencing, a process termed zygotic genome activation (ZGA). Despite its fundamental role in totipotency and the initiation of life, the precise mechanism underlying ZGA initiation remains unclear. The existence of minor ZGA implies the possible critical role of noncoding RNAs in the initiation of ZGA. Here, we delineate the expression profile of long noncoding RNAs (lncRNAs) in early mouse embryonic development and elucidate their critical role in minor ZGA. Compared with protein-coding genes (PCGs), lncRNAs exhibit a stronger correlation with minor ZGA. Distinct H3K9me3 profiles can be observed between lncRNA genes and PCGs, and the enrichment of H3K9me3 before ZGA might explain the suspended expression of major ZGA-related PCGs despite possessing PolII pre-configuration. Furthermore, we identified the presence of PolII-enriched MuERV-L around the transcriptional start site of minor ZGA-related lncRNAs, and these repeats are responsible for the activation of minor ZGA-related lncRNAs and subsequent embryo development. Our study suggests that MuERV-L mediates minor ZGA lncRNA activation as a critical driver between epigenetic reprogramming triggered by fertilization and the embryo developmental program, thus providing clues for understanding the regulatory mechanism of totipotency and establishing bona fide totipotent stem cells.

Gonadal transcriptome sequencing reveals sexual dimorphism in expression profiling of sex-related genes in Asian arowana (Scleropages formosus)

Asia arowana (Scleropages formosus) is an ornamental fish with high economic value, while its sex determination mechanism is still poorly understood. By far, no morphological evidence or molecular marker has been developed for effective distinguishment of genders, which poses a critical challenge to our captive breeding efforts. In this study, we sequenced gonadal transcriptomes of adult Asian arowanas and revealed differential expression profiling of sex-related genes. Based on the comparative transcriptomics analysis of testes (n = 3) and ovaries (n = 3), we identified a total of 8,872 differentially expressed genes (DEGs) and 18,490 differentially expressed transposable elements (TEs) between male and female individuals. Interestingly, the expression of TEs usually has been more significantly testis-biased than related coding genes. As expected, several genes related to females (such as foxl2 and cyp19a1a) are significantly transcribed in the ovary, and some genes related to male gonad development (such as dmrt1, gsdf and amh) are highly expressed in the testis. This sexual dimorphism is valuable for ascertaining the differential expression patterns of sex-related genes and enriching the genetic resources of this economically important species. These valuable genetic materials thereby provide instructive references for gender identification and one-to-one breeding practices so as to expand fish numbers for a rapid elevation of economic value.

Jump-starting life: balancing transposable element co-option and genome integrity in the developing mammalian embryo

Remnants of transposable elements (TEs) are widely expressed throughout mammalian embryo development. Originally infesting our genomes as selfish elements and acting as a source of genome instability, several of these elements have been co-opted as part of a complex system of genome regulation. Many TEs have lost transposition ability and their transcriptional potential has been tampered as a result of interactions with the host throughout evolutionary time. It has been proposed that TEs have been ultimately repurposed to function as gene regulatory hubs scattered throughout our genomes. In the early embryo in particular, TEs find a perfect environment of naïve chromatin to escape transcriptional repression by the host. As a consequence, it is thought that hosts found ways to co-opt TE sequences to regulate large-scale changes in chromatin and transcription state of their genomes. In this review, we discuss several examples of TEs expressed during embryo development, their potential for co-option in genome regulation and the evolutionary pressures on TEs and on our genomes.

Dynamic co-evolution of transposable elements and the piRNA pathway in African cichlid fishes

East African cichlid fishes have diversified in an explosive fashion, but the (epi)genetic basis of the phenotypic diversity of these fishes remains largely unknown. Although transposable elements (TEs) have been associated with phenotypic variation in cichlids, little is known about their transcriptional activity and epigenetic silencing. Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed an expansion of piwil1 genes in Lake Malawi cichlids, likely driven by PiggyBac TEs. The expanded piwil1 copies have signatures of positive selection and retain amino acid residues essential for catalytic activity. Furthermore, the gonads of African cichlids express a Piwi-interacting RNA (piRNA) pathway that target TEs. We define the genomic sites of piRNA production in African cichlids and find divergence in closely related species, in line with fast evolution of piRNA-producing loci. Our findings suggest dynamic co-evolution of TEs and host silencing pathways in the African cichlid radiations. We propose that this co-evolution has contributed to cichlid genomic diversity.

Regulation and function of transposable elements in cancer genomes

Over half of human genomic DNA is composed of repetitive sequences generated throughout evolution by prolific mobile genetic parasites called transposable elements (TEs). Long disregarded as "junk" or "selfish" DNA, TEs are increasingly recognized as formative elements in genome evolution, wired intimately into the structure and function of the human genome. Advances in sequencing technologies and computational methods have ushered in an era of unprecedented insight into how TE activity impacts human biology in health and disease. Here we discuss the current views on how TEs have shaped the regulatory landscape of the human genome, how TE activity is implicated in human cancers, and how recent findings motivate novel strategies to leverage TE activity for improved cancer therapy. Given the crucial role of methodological advances in TE biology, we pair our conceptual discussions with an in-depth review of the inherent technical challenges in studying repeats, specifically related to structural variation, expression analyses, and chromatin regulation. Lastly, we provide a catalog of existing and emerging assays and bioinformatic software that altogether are enabling the most sophisticated and comprehensive investigations yet into the regulation and function of interspersed repeats in cancer genomes.

Select EZH2 inhibitors enhance viral mimicry effects of DNMT inhibition through a mechanism involving NFAT:AP-1 signaling

DNA methyltransferase inhibitor (DNMTi) efficacy in solid tumors is limited. Colon cancer cells exposed to DNMTi accumulate lysine-27 trimethylation on histone H3 (H3K27me3). We propose this Enhancer of Zeste Homolog 2 (EZH2)-dependent repressive modification limits DNMTi efficacy. Here, we show that low-dose DNMTi treatment sensitizes colon cancer cells to select EZH2 inhibitors (EZH2is). Integrative epigenomic analysis reveals that DNMTi-induced H3K27me3 accumulates at genomic regions poised with EZH2. Notably, combined EZH2i and DNMTi alters the epigenomic landscape to transcriptionally up-regulate the calcium-induced nuclear factor of activated T cells (NFAT):activating protein 1 (AP-1) signaling pathway. Blocking this pathway limits transcriptional activating effects of these drugs, including transposable element and innate immune response gene expression involved in viral defense. Analysis of primary human colon cancer specimens reveals positive correlations between DNMTi-, innate immune response-, and calcium signaling-associated transcription profiles. Collectively, we show that compensatory EZH2 activity limits DNMTi efficacy in colon cancer and link NFAT:AP-1 signaling to epigenetic therapy-induced viral mimicry.

Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

Activation of human endogenous retroviruses and its physiological consequences

Human endogenous retroviruses (HERVs) are abundant sequences that persist within the human genome as remnants of ancient retroviral infections. These sequences became fixed and accumulate mutations or deletions over time. HERVs have affected human evolution and physiology by providing a unique repertoire of coding and non-coding sequences to the genome. In healthy individuals, HERVs participate in immune responses, formation of syncytiotrophoblasts and cell-fate specification. In this Review, we discuss how endogenized retroviral motifs and regulatory sequences have been co-opted into human physiology and how they are tightly regulated. Infections and mutations can derail this regulation, leading to differential HERV expression, which may contribute to pathologies including neurodegeneration, pathological inflammation and oncogenesis. Emerging evidence demonstrates that HERVs are crucial to human health and represent an understudied facet of many diseases, and we therefore argue that investigating their fundamental properties could improve existing therapies and help develop novel therapeutic strategies.

XRN1 deletion induces PKR-dependent cell lethality in interferon-activated cancer cells

Emerging data suggest that induction of viral mimicry responses through activation of double-stranded RNA (dsRNA) sensors in cancer cells is a promising therapeutic strategy. One approach to induce viral mimicry is to target molecular regulators of dsRNA sensing pathways. Here, we show that the exoribonuclease XRN1 is a negative regulator of the dsRNA sensor protein kinase R (PKR) in cancer cells with high interferon-stimulated gene expression. XRN1 deletion causes PKR pathway activation and consequent cancer cell lethality. Disruption of interferon signaling with the JAK1/2 inhibitor ruxolitinib can decrease cellular PKR levels and rescue sensitivity to XRN1 deletion. Conversely, interferon-β stimulation can increase PKR levels and induce sensitivity to XRN1 inactivation. Lastly, XRN1 deletion causes accumulation of endogenous complementary sense/anti-sense RNAs, which may represent candidate PKR ligands. Our data demonstrate how XRN1 regulates PKR and how this interaction creates a vulnerability in cancer cells with an activated interferon cell state.

Towards targeting transposable elements for cancer therapy

Transposable elements (TEs) represent almost half of the human genome. Historically deemed 'junk DNA', recent technological advancements have stimulated a wave of research into the functional impact of TEs on gene-regulatory networks in evolution and development, as well as in diseases including cancer. The genetic and epigenetic evolution of cancer involves the exploitation of TEs, whereby TEs contribute directly to cancer-specific gene activities. This Review provides a perspective on the role of TEs in cancer as being a 'double-edged sword', both promoting cancer evolution and representing a vulnerability that could be exploited in cancer therapy. We discuss how TEs affect transcriptome regulation and other cellular processes in cancer. We highlight the potential of TEs as therapeutic targets for cancer. We also summarize technical hurdles in the characterization of TEs with genomic assays. Last, we outline open questions and exciting future research avenues.

Significant Variations in Double-Stranded RNA Levels in Cultured Skin Cells

Endogenous double-stranded RNA has emerged as a potent stimulator of innate immunity. Under physiological conditions, endogenous dsRNA is maintained in the cell nucleus or the mitochondria; however, if protective mechanisms are breached, it leaches into the cytoplasm and triggers immune signaling pathways. Ectopic activation of innate immune pathways is associated with various diseases and senescence and can trigger apoptosis. Hereby, the level of cytoplasmic dsRNA is crucial. We have enriched dsRNA from two melanoma cell lines and primary dermal fibroblasts, including a competing probe, and analyzed the dsRNA transcriptome using RNA sequencing. There was a striking difference in read counts between the cell lines and the primary cells, and the effect was confirmed by northern blotting and immunocytochemistry. Both mitochondria (10-20%) and nuclear transcription (80-90%) contributed significantly to the dsRNA transcriptome. The mitochondrial contribution was lower in the cancer cells compared to fibroblasts. The expression of different transposable element families was comparable, suggesting a general up-regulation of transposable element expression rather than stimulation of a specific sub-family. Sequencing of the input control revealed minor differences in dsRNA processing pathways with an upregulation of oligoadenylate synthase and RNP125 that negatively regulates the dsRNA sensors RIG1 and MDA5. Moreover, RT-qPCR, Western blotting, and immunocytochemistry confirmed the relatively minor adaptations to the hugely different dsRNA levels. As a consequence, these transformed cell lines are potentially less tolerant to interventions that increase the formation of endogenous dsRNA.

Transposable elements acquire time- and sex-specific transcriptional and epigenetic signatures along mouse fetal gonad development

Gonadal sex determination in mice is a complex and dynamic process, which is crucial for the development of functional reproductive organs. The expression of genes involved in this process is regulated by a variety of genetic and epigenetic mechanisms. Recently, there has been increasing evidence that transposable elements (TEs), which are a class of mobile genetic elements, play a significant role in regulating gene expression during embryogenesis and organ development. In this study, we aimed to investigate the involvement of TEs in the regulation of gene expression during mouse embryonic gonadal development. Through bioinformatics analysis, we aimed to identify and characterize specific TEs that operate as regulatory elements for sex-specific genes, as well as their potential mechanisms of regulation. We identified TE loci expressed in a time- and sex-specific manner along fetal gonad development that correlate positively and negatively with nearby gene expression, suggesting that their expression is integrated to the gonadal regulatory network. Moreover, chromatin accessibility and histone post-transcriptional modification analyses in differentiating supporting cells revealed that TEs are acquiring a sex-specific signature for promoter-, enhancer-, and silencer-like elements, with some of them being proximal to critical sex-determining genes. Altogether, our study introduces TEs as the new potential players in the gene regulatory network that controls gonadal development in mammals.

LINE-1 retrotransposition and its deregulation in cancers: implications for therapeutic opportunities

Long interspersed element 1 (LINE-1) is the only protein-coding transposon that is active in humans. LINE-1 propagates in the genome using RNA intermediates via retrotransposition. This activity has resulted in LINE-1 sequences occupying approximately one-fifth of our genome. Although most copies of LINE-1 are immobile, ∼100 copies are retrotransposition-competent. Retrotransposition is normally limited via epigenetic silencing, DNA repair, and other host defense mechanisms. In contrast, LINE-1 overexpression and retrotransposition are hallmarks of cancers. Here, we review mechanisms of LINE-1 regulation and how LINE-1 may promote genetic heterogeneity in tumors. Finally, we discuss therapeutic strategies to exploit LINE-1 biology in cancers.

An Analysis of a Transposable Element Expression Atlas during 27 Developmental Stages in Porcine Skeletal Muscle: Unveiling Molecular Insights into Pork Production Traits

The development and growth of porcine skeletal muscle determine pork quality and yield. While genetic regulation of porcine skeletal muscle development has been extensively studied using various omics data, the role of transposable elements (TEs) in this context has been less explored. To bridge this gap, we constructed a comprehensive atlas of TE expression throughout the developmental stages of porcine skeletal muscle. This was achieved by integrating porcine TE genomic coordinates with whole-transcriptome RNA-Seq data from 27 developmental stages. We discovered that in pig skeletal muscle, active Tes are closely associated with active epigenomic marks, including low levels of DNA methylation, high levels of chromatin accessibility, and active histone modifications. Moreover, these TEs include 6074 self-expressed TEs that are significantly enriched in terms of muscle cell development and myofibril assembly. Using the TE expression data, we conducted a weighted gene co-expression network analysis (WGCNA) and identified a module that is significantly associated with muscle tissue development as well as genome-wide association studies (GWAS) of the signals of pig meat and carcass traits. Within this module, we constructed a TE-mediated gene regulatory network by adopting a unique multi-omics integration approach. This network highlighted several established candidate genes associated with muscle-relevant traits, including HES6, CHRNG, ACTC1, CHRND, MAMSTR, and PER2, as well as novel genes like ENSSSCG00000005518, ENSSSCG00000033601, and PIEZO2. These novel genes hold promise for regulating muscle-related traits in pigs. In summary, our research not only enhances the TE-centered dissection of the genetic basis underlying pork production traits, but also offers a general approach for constructing TE-mediated regulatory networks to study complex traits or diseases.

Selective control of transposable element expression during T cell exhaustion and anti-PD-1 treatment

In chronic infections and cancer, T cells are exposed to prolonged antigen stimulation, resulting in loss of function (or exhaustion) and impairment of effective immunological protection. Exhausted T cells are heterogeneous and include early progenitors (Tpex) and terminally exhausted cells (Tex). Here, we used bulk and single-cell transcriptomics to analyze expression of transposable elements (TEs) in subpopulations of mouse and human CD8+ tumor-infiltrating T lymphocytes (TILs). We show that in mice, members of the virus-like murine VL30 TE family (mostly intact, evolutionary young ERV1s) are strongly repressed in terminally exhausted CD8+ T cells in both tumor and viral models of exhaustion. Tpex expression of these VL30s, which are mainly intergenic and transcribed independently of their closest gene neighbors, was driven by Fli1, a transcription factor involved in progression from Tpex to Tex. Immune checkpoint blockade (ICB) in both mice and patients with cancer increased TE expression (including VL30 in mice), demonstrating that TEs may be applicable as ICB response biomarkers. We conclude that expression of TEs is tightly regulated in TILs during establishment of exhaustion and reprogramming by ICB. Analyses of TE expression on single cells and bulk populations open opportunities for understanding immune cell identity and heterogeneity, as well as for defining cellular gene expression signatures and disease biomarkers.

Vitamin C activates young LINE-1 elements in mouse embryonic stem cells via H3K9me3 demethylation

BACKGROUND

Vitamin C (vitC) enhances the activity of 2-oxoglutarate-dependent dioxygenases, including TET enzymes, which catalyse DNA demethylation, and Jumonji-domain histone demethylases. The epigenetic remodelling promoted by vitC improves the efficiency of induced pluripotent stem cell derivation, and is required to attain a ground-state of pluripotency in embryonic stem cells (ESCs) that closely mimics the inner cell mass of the early blastocyst. However, genome-wide DNA and histone demethylation can lead to upregulation of transposable elements (TEs), and it is not known how vitC addition in culture media affects TE expression in pluripotent stem cells.

RESULTS

Here we show that vitC increases the expression of several TE families, including evolutionarily young LINE-1 (L1) elements, in mouse ESCs. We find that TET activity is dispensable for L1 upregulation, and that instead it occurs largely as a result of H3K9me3 loss mediated by KDM4A/C histone demethylases. Despite increased L1 levels, we did not detect increased somatic insertion rates in vitC-treated cells. Notably, treatment of human ESCs with vitC also increases L1 protein levels, albeit through a distinct, post-transcriptional mechanism.

CONCLUSION

VitC directly modulates the expression of mouse L1s and other TEs through epigenetic mechanisms, with potential for downstream effects related to the multiple emerging roles of L1s in cellular function.

capTEs enables locus-specific dissection of transcriptional outputs from reference and nonreference transposable elements

Transposable elements (TEs) serve as both insertional mutagens and regulatory elements in cells, and their aberrant activity is increasingly being revealed to contribute to diseases and cancers. However, measuring the transcriptional consequences of nonreference and young TEs at individual loci remains challenging with current methods, primarily due to technical limitations, including short read lengths generated and insufficient coverage in target regions. Here, we introduce a long-read targeted RNA sequencing method, Cas9-assisted profiling TE expression sequencing (capTEs), for quantitative analysis of transcriptional outputs for individual TEs, including transcribed nonreference insertions, noncanonical transcripts from various transcription patterns and their correlations with expression changes in related genes. This method selectively identified TE-containing transcripts and outputted data with up to 90% TE reads, maintaining a comparable data yield to whole-transcriptome sequencing. We applied capTEs to human cancer cells and found that internal and inserted Alu elements may employ distinct regulatory mechanisms to upregulate gene expression. We expect that capTEs will be a critical tool for advancing our understanding of the biological functions of individual TEs at the locus level, revealing their roles as both mutagens and regulators in biological and pathogenic processes.

Transposable elements as tissue-specific enhancers in cancers of endodermal lineage

Transposable elements (TE) are repetitive genomic elements that harbor binding sites for human transcription factors (TF). A regulatory role for TEs has been suggested in embryonal development and diseases such as cancer but systematic investigation of their functions has been limited by their widespread silencing in the genome. Here, we utilize unbiased massively parallel reporter assay data using a whole human genome library to identify TEs with functional enhancer activity in two human cancer types of endodermal lineage, colorectal and liver cancers. We show that the identified TE enhancers are characterized by genomic features associated with active enhancers, such as epigenetic marks and TF binding. Importantly, we identify distinct TE subfamilies that function as tissue-specific enhancers, namely MER11- and LTR12-elements in colon and liver cancers, respectively. These elements are bound by distinct TFs in each cell type, and they have predicted associations to differentially expressed genes. In conclusion, these data demonstrate how different cancer types can utilize distinct TEs as tissue-specific enhancers, paving the way for comprehensive understanding of the role of TEs as bona fide enhancers in the cancer genomes.

The Regulation and Immune Signature of Retrotransposons in Cancer

Advances in sequencing technologies and the bioinformatic analysis of big data facilitate the study of jumping genes' activity in the human genome in cancer from a broad perspective. Retrotransposons, which move from one genomic site to another by a copy-and-paste mechanism, are regulated by various molecular pathways that may be disrupted during tumorigenesis. Active retrotransposons can stimulate type I IFN responses. Although accumulated evidence suggests that retrotransposons can induce inflammation, the research investigating the exact mechanism of triggering these responses is ongoing. Understanding these mechanisms could improve the therapeutic management of cancer through the use of retrotransposon-induced inflammation as a tool to instigate immune responses to tumors.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in several cancer types including sarcomas, which are mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors (Pparγ and Cebpα) and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks at putative enhancer elements and promoters. Finally, we observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Our results demonstrate that ATRX functions to buffer against differentiation in mesenchymal progenitor cells, which has implications for understanding ATRX loss of function in sarcomas.

XRN1 deletion induces PKR-dependent cell lethality in interferon-activated cancer cells

Emerging data suggest that induction of viral mimicry responses through activation of double-stranded RNA (dsRNA) sensors in cancer cells is a promising therapeutic strategy. One approach to induce viral mimicry is to target molecular regulators of dsRNA sensing pathways. Here, we show that the exoribonuclease XRN1 is a negative regulator of the dsRNA sensor protein kinase R (PKR) in cancer cells with high interferon-stimulated gene (ISG) expression. XRN1 deletion causes PKR activation and consequent cancer cell lethality. Disruption of interferon signaling with the JAK1/2 inhibitor ruxolitinib can decrease cellular PKR levels and rescue sensitivity to XRN1 deletion. Conversely, interferon-β stimulation can increase PKR levels and induce sensitivity to XRN1 inactivation. Lastly, XRN1 deletion causes accumulation of endogenous complementary sense/anti-sense RNAs, which may represent candidate PKR ligands. Our data demonstrate how XRN1 regulates PKR and nominate XRN1 as a potential therapeutic target in cancer cells with an activated interferon cell state.

Lack of Stim2 Affects Vision-Dependent Behavior and Sensitivity to Hypoxia

Stromal interaction molecules (STIMs) are endoplasmic reticulum-resident proteins that regulate Ca2+ homeostasis and signaling by store-operated calcium entry (SOCE). The different properties and functions of STIM1 and STIM2 have been described mostly based on work in vitro. STIM2 knockout mice do not survive until adulthood. Therefore, we generated and characterized stim2a and stim2b double-knockout zebrafish. The (stim2a;stim2b)-/- zebrafish did not have any apparent morphological phenotype. However, RNA sequencing revealed 1424 differentially expressed genes. One of the most upregulated genes was annexin A3a, which is a marker of activated microglia. This corresponded well to an increase in Neutral Red staining in the in vivo imaging of the (stim2a;stim2b)-/- zebrafish brain. The lack of Stim2 decreased zebrafish survival under low oxygen conditions. Behavioral tests, such as the visual-motor response test and dark-light preference test, indicated that (stim2a;stim2b)-/- larvae might have problems with vision. This was consistent with the downregulation of many genes that are related to light perception. The periodic acid-Schiff staining of retina sections from adult zebrafish revealed alterations of the stratum pigmentosum, suggesting the involvement of a Stim2-dependent process in visual perception. Altogether, these data reveal new functions for Stim2 in the nervous system.

Loss of epigenetic suppression of retrotransposons with oncogenic potential in aging mammary luminal epithelial cells

A primary function of DNA methylation in mammalian genomes is to repress transposable elements (TEs). The widespread methylation loss that is commonly observed in cancer cells results in the loss of epigenetic repression of TEs. The aging process is similarly characterized by changes to the methylome. However, the impact of these epigenomic alterations on TE silencing and the functional consequences of this have remained unclear. To assess the epigenetic regulation of TEs in aging, we profiled DNA methylation in human mammary luminal epithelial cells (LEps)-a key cell lineage implicated in age-related breast cancers-from younger and older women. We report here that several TE subfamilies function as regulatory elements in normal LEps, and a subset of these display consistent methylation changes with age. Methylation changes at these TEs occurred at lineage-specific transcription factor binding sites, consistent with loss of lineage specificity. Whereas TEs mainly showed methylation loss, CpG islands (CGIs) that are targets of the Polycomb repressive complex 2 (PRC2) show a gain of methylation in aging cells. Many TEs with methylation loss in aging LEps have evidence of regulatory activity in breast cancer samples. We furthermore show that methylation changes at TEs impact the regulation of genes associated with luminal breast cancers. These results indicate that aging leads to DNA methylation changes at TEs that undermine the maintenance of lineage specificity, potentially increasing susceptibility to breast cancer.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

Immunogenicity in renal cell carcinoma: shifting focus to alternative sources of tumour-specific antigens

Renal cell carcinoma (RCC) comprises a group of malignancies arising from the kidney with unique tumour-specific antigen (TSA) signatures that can trigger cytotoxic immunity. Two classes of TSAs are now considered potential drivers of immunogenicity in RCC: small-scale insertions and deletions (INDELs) that result in coding frameshift mutations, and activation of human endogenous retroviruses. The presence of neoantigen-specific T cells is a hallmark of solid tumours with a high mutagenic burden, which typically have abundant TSAs owing to non-synonymous single nucleotide variations within the genome. However, RCC exhibits high cytotoxic T cell reactivity despite only having an intermediate non-synonymous single nucleotide variation mutational burden. Instead, RCC tumours have a high pan-cancer proportion of INDEL frameshift mutations, and coding frameshift INDELs are associated with high immunogenicity. Moreover, cytotoxic T cells in RCC subtypes seem to recognize tumour-specific endogenous retrovirus epitopes, whose presence is associated with clinical responses to immune checkpoint blockade therapy. Here, we review the distinct molecular landscapes in RCC that promote immunogenic responses, discuss clinical opportunities for discovery of biomarkers that can inform therapeutic immune checkpoint blockade strategies, and identify gaps in knowledge for future investigations.

Exploratory analysis of L1 retrotransposons expression in autism

BACKGROUND

Autism spectrum disorder (ASD) is a set of highly heterogeneous neurodevelopmental diseases whose genetic etiology is not completely understood. Several investigations have relied on transcriptome analysis from peripheral tissues to dissect ASD into homogenous molecular phenotypes. Recently, analysis of changes in gene expression from postmortem brain tissues has identified sets of genes that are involved in pathways previously associated with ASD etiology. In addition to protein-coding transcripts, the human transcriptome is composed by a large set of non-coding RNAs and transposable elements (TEs). Advancements in sequencing technologies have proven that TEs can be transcribed in a regulated fashion, and their dysregulation might have a role in brain diseases.

METHODS

We exploited published datasets comprising RNA-seq data from (1) postmortem brain of ASD subjects, (2) in vitro cell cultures where ten different ASD-relevant genes were knocked out and (3) blood of discordant siblings. We measured the expression levels of evolutionarily young full-length transposable L1 elements and characterized the genomic location of deregulated L1s assessing their potential impact on the transcription of ASD-relevant genes. We analyzed every sample independently, avoiding to pool together the disease subjects to unmask the heterogeneity of the molecular phenotypes.

RESULTS

We detected a strong upregulation of intronic full-length L1s in a subset of postmortem brain samples and in in vitro differentiated neurons from iPSC knocked out for ATRX. L1 upregulation correlated with an high number of deregulated genes and retained introns. In the anterior cingulate cortex of one subject, a small number of significantly upregulated L1s overlapped with ASD-relevant genes that were significantly downregulated, suggesting the possible existence of a negative effect of L1 transcription on host transcripts.

LIMITATIONS

Our analyses must be considered exploratory and will need to be validated in bigger cohorts. The main limitation is given by the small sample size and by the lack of replicates for postmortem brain samples. Measuring the transcription of locus-specific TEs is complicated by the repetitive nature of their sequence, which reduces the accuracy in mapping sequencing reads to the correct genomic locus.

CONCLUSIONS

L1 upregulation in ASD appears to be limited to a subset of subjects that are also characterized by a general deregulation of the expression of canonical genes and an increase in intron retention. In some samples from the anterior cingulate cortex, L1s upregulation seems to directly impair the expression of some ASD-relevant genes by a still unknown mechanism. L1s upregulation may therefore identify a group of ASD subjects with common molecular features and helps stratifying individuals for novel strategies of therapeutic intervention.

One-carbon metabolism is required for epigenetic stability in the mouse placenta

One-carbon metabolism, including the folate cycle, has a crucial role in fetal development though its molecular function is complex and unclear. The hypomorphic Mtrr ^gt allele is known to disrupt one-carbon metabolism, and thus methyl group availability, leading to several developmental phenotypes (e.g., neural tube closure defects, fetal growth anomalies). Remarkably, previous studies showed that some of the phenotypes were transgenerationally inherited. Here, we explored the genome-wide epigenetic impact of one-carbon metabolism in placentas associated with fetal growth phenotypes and determined whether specific DNA methylation changes were inherited. Firstly, methylome analysis of Mtrr ^gt/gt homozygous placentas revealed genome-wide epigenetic instability. Several differentially methylated regions (DMRs) were identified including at the Cxcl1 gene promoter and at the En2 gene locus, which may have phenotypic implications. Importantly, we discovered hypomethylation and ectopic expression of a subset of ERV elements throughout the genome of Mtrr ^gt/gt placentas with broad implications for genomic stability. Next, we determined that known spermatozoan DMRs in Mtrr ^gt/gt males were reprogrammed in the placenta with little evidence of direct or transgenerational germline DMR inheritance. However, some spermatozoan DMRs were associated with placental gene misexpression despite normalisation of DNA methylation, suggesting the inheritance of an alternative epigenetic mechanism. Integration of published wildtype histone ChIP-seq datasets with Mtrr ^gt/gt spermatozoan methylome and placental transcriptome datasets point towards H3K4me3 deposition at key loci. These data suggest that histone modifications might play a role in epigenetic inheritance in this context. Overall, this study sheds light on the mechanistic complexities of one-carbon metabolism in development and epigenetic inheritance.

In silico characterisation of minor wave genes and LINE-1s transcriptional dynamics at murine zygotic genome activation

Introduction: In mouse, the zygotic genome activation (ZGA) is coordinated by MERVL elements, a class of LTR retrotransposons. In addition to MERVL, another class of retrotransposons, LINE-1 elements, recently came under the spotlight as key regulators of murine ZGA. In particular, LINE-1 transcripts seem to be required to switch-off the transcriptional program started by MERVL sequences, suggesting an antagonistic interplay between LINE-1 and MERVL pathways. Methods: To better investigate the activities of LINE-1 and MERVL elements at ZGA, we integrated publicly available transcriptomics (RNA-seq), chromatin accessibility (ATAC-seq) and Pol-II binding (Stacc-seq) datasets and characterised the transcriptional and epigenetic dynamics of such elements during murine ZGA. Results: We identified two likely distinct transcriptional activities characterising the murine zygotic genome at ZGA onset. On the one hand, our results confirmed that ZGA minor wave genes are preferentially transcribed from MERVL-rich and gene-dense genomic compartments, such as gene clusters. On the other hand, we identified a set of evolutionary young and likely transcriptionally autonomous LINE-1s located in intergenic and gene-poor regions showing, at the same stage, features such as open chromatin and RNA Pol II binding suggesting them to be, at least, poised for transcription. Discussion: These results suggest that, across evolution, transcription of two different classes of transposable elements, MERVLs and LINE-1s, have likely been confined in genic and intergenic regions respectively in order to maintain and regulate two successive transcriptional programs at ZGA.

Cancer cells co-evolve with retrotransposons to mitigate viral mimicry

Overexpression of repetitive elements is an emerging hallmark of human cancers ¹ . Diverse repeats can mimic viruses by replicating within the cancer genome through retrotransposition, or presenting pathogen-associated molecular patterns (PAMPs) to the pattern recognition receptors (PRRs) of the innate immune system ^2-5 . Yet, how specific repeats affect tumor evolution and shape the tumor immune microenvironment (TME) in a pro- or anti-tumorigenic manner remains poorly defined. Here, we integrate whole genome and total transcriptome data from a unique autopsy cohort of multiregional samples collected in pancreatic ductal adenocarcinoma (PDAC) patients, into a comprehensive evolutionary analysis. We find that more recently evolved S hort I nterspersed N uclear E lements (SINE), a family of retrotransposable repeats, are more likely to form immunostimulatory double-strand RNAs (dsRNAs). Consequently, younger SINEs are strongly co-regulated with RIG-I like receptor associated type-I interferon genes but anti-correlated with pro-tumorigenic macrophage infiltration. We discover that immunostimulatory SINE expression in tumors is regulated by either L ong I nterspersed N uclear E lements 1 (LINE1/L1) mobility or ADAR1 activity in a TP53 mutation dependent manner. Moreover, L1 retrotransposition activity tracks with tumor evolution and is associated with TP53 mutation status. Altogether, our results suggest pancreatic tumors actively evolve to modulate immunogenic SINE stress and induce pro-tumorigenic inflammation. Our integrative, evolutionary analysis therefore illustrates, for the first time, how dark matter genomic repeats enable tumors to co-evolve with the TME by actively regulating viral mimicry to their selective advantage.

Low impact of polyploidization on the transcriptome of synthetic allohexaploid wheat

BACKGROUND

Bread wheat is a recent allohexaploid (genomic constitution AABBDD) that emerged through a hybridization between tetraploid Triticum turgidum (AABB) and diploid Aegilops tauschii (DD) less than 10,000 years ago. The hexaploidization can be re-created artificially, producing synthetic wheat that has been used to study immediate genomic responses to polyploidization. The scale of the consequences of polyploidization, and their mechanism of establishment, remain uncertain.

RESULTS

Here we sampled several synthetic wheats from alternative parental genotypes and reciprocal crosses, and examined transcriptomes from two different tissues and successive generations. We did not detect any massive reprogramming in gene expression, with only around 1% of expressed genes showing significant differences compared to their lower-ploidy parents. Most of this differential expression is located on the D subgenome, without consistency in the direction of the expression change. Homoeolog expression bias in synthetic wheat is similar to the pattern observed in the parents. Both differential expression and homoeolog bias are tissue-specific. While up to three families of transposable elements became upregulated in wheat synthetics, their position and distance are not significantly associated with expression changes in proximal genes.

DISCUSSION

While only a few genes change their expression pattern after polyploidization, they can be involved in agronomically important pathways. Alternative parental combinations can lead to opposite changes on the same subset of D-located genes, which is relevant for harnessing new diversity in wheat breeding. Tissue specificity of the polyploidization-triggered expression changes indicates the remodelling of transcriptomes in synthetic wheat is plastic and likely caused by regulome interactions rather than permanent changes. We discuss the pitfalls of transcriptomic comparisons across ploidy levels that can inflate the de-regulation signal.

CONCLUSIONS

Transcriptomic response to polyploidization in synthetic AABBDD wheat is modest and much lower than some previous estimates. Homoeolog expression bias in wheat allohexaploids is mostly attributed to parental legacy, with polyploidy having a mild balancing effect.

Regulation of human trophoblast gene expression by endogenous retroviruses

The placenta is a fast-evolving organ with large morphological and histological differences across eutherians, but the genetic changes driving placental evolution have not been fully elucidated. Transposable elements, through their capacity to quickly generate genetic variation and affect host gene regulation, may have helped to define species-specific trophoblast gene expression programs. Here we assess the contribution of transposable elements to human trophoblast gene expression as enhancers or promoters. Using epigenomic data from primary human trophoblast and trophoblast stem-cell lines, we identified multiple endogenous retrovirus families with regulatory potential that lie close to genes with preferential expression in trophoblast. These largely primate-specific elements are associated with inter-species gene expression differences and are bound by transcription factors with key roles in placental development. Using genetic editing, we demonstrate that several elements act as transcriptional enhancers of important placental genes, such as CSF1R and PSG5. We also identify an LTR10A element that regulates ENG expression, affecting secretion of soluble endoglin, with potential implications for preeclampsia. Our data show that transposons have made important contributions to human trophoblast gene regulation, and suggest that their activity may affect pregnancy outcomes.

RUNX1 is required in granulocyte-monocyte progenitors to attenuate inflammatory cytokine production by neutrophils

The transcription factor RUNX1 is mutated in familial platelet disorder with associated myeloid malignancies (FPDMM) and in sporadic myelodysplastic syndrome and leukemia. RUNX1 regulates inflammation in multiple cell types. Here we show that RUNX1 is required in granulocyte-monocyte progenitors (GMPs) to restrict the inflammatory response of neutrophils to toll-like receptor 4 (TLR4) signaling. Loss of RUNX1 in GMPs increased the TLR4 coreceptor CD14 on neutrophils, which contributed to neutrophils’ increased inflammatory cytokine production in response to the TLR4 ligand lipopolysaccharide. RUNX1 loss increased the chromatin accessibility of retrotransposons in GMPs and neutrophils and induced a type I interferon signature characterized by enriched footprints for signal transducer and activator of transcription (STAT1::STAT2) and interferon regulatory factors (IRF) in opened chromatin, and increased expression of interferon-stimulated genes. The overproduction of inflammatory cytokines by neutrophils was reversed by inhibitors of type I IFN signaling. We conclude that RUNX1 restrains the chromatin accessibility of retrotransposons in GMPs and neutrophils, and that loss of RUNX1 increases proinflammatory cytokine production by elevating tonic type I interferon signaling.

Epigenetic resetting in the human germ line entails histone modification remodeling

Epigenetic resetting in the mammalian germ line entails acute DNA demethylation, which lays the foundation for gametogenesis, totipotency, and embryonic development. We characterize the epigenome of hypomethylated human primordial germ cells (hPGCs) to reveal mechanisms preventing the widespread derepression of genes and transposable elements (TEs). Along with the loss of DNA methylation, we show that hPGCs exhibit a profound reduction of repressive histone modifications resulting in diminished heterochromatic signatures at most genes and TEs and the acquisition of a neutral or paused epigenetic state without transcriptional activation. Efficient maintenance of a heterochromatic state is limited to a subset of genomic loci, such as evolutionarily young TEs and some developmental genes, which require H3K9me3 and H3K27me3, respectively, for efficient transcriptional repression. Accordingly, transcriptional repression in hPGCs presents an exemplary balanced system relying on local maintenance of heterochromatic features and a lack of inductive cues.

Molecular subtypes of ALS are associated with differences in patient prognosis

Amyotrophic Lateral Sclerosis (ALS) is a neurodegenerative disease with poorly understood clinical heterogeneity, underscored by significant differences in patient age at onset, symptom progression, therapeutic response, disease duration, and comorbidity presentation. We perform a patient stratification analysis to better understand the variability in ALS pathology, utilizing postmortem frontal and motor cortex transcriptomes derived from 208 patients. Building on the emerging role of transposable element (TE) expression in ALS, we consider locus-specific TEs as distinct molecular features during stratification. Here, we identify three unique molecular subtypes in this ALS cohort, with significant differences in patient survival. These results suggest independent disease mechanisms drive some of the clinical heterogeneity in ALS.