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

LOXL2-mediated chromatin compaction is required to maintain the oncogenic properties of triple-negative breast cancer cells

Oxidation of histone H3 at lysine 4 (H3K4ox) is catalyzed by lysyl oxidase homolog 2 (LOXL2). This histone modification is enriched in heterochromatin in triple-negative breast cancer (TNBC) cells and has been linked to the maintenance of compacted chromatin. However, the molecular mechanism underlying this maintenance is still unknown. Here, we show that LOXL2 interacts with RuvB-Like 1 (RUVBL1), RuvB-Like 2 (RUVBL2), Actin-like protein 6A (ACTL6A), and DNA methyltransferase 1associated protein 1 (DMAP1), a complex involved in the incorporation of the histone variant H2A.Z. Our experiments indicate that this interaction and the active form of RUVBL2 are required to maintain LOXL2-dependent chromatin compaction. Genome-wide experiments showed that H2A.Z, RUVBL2, and H3K4ox colocalize in heterochromatin regions. In the absence of LOXL2 or RUVBL2, global levels of the heterochromatin histone mark H3K9me3 were strongly reduced, and the ATAC-seq signal in the H3K9me3 regions was increased. Finally, we observed that the interplay between these series of events is required to maintain H3K4ox-enriched heterochromatin regions, which in turn is key for maintaining the oncogenic properties of the TNBC cell line tested (MDA-MB-231).

Effect of a FOXO1 inhibitor on trophoblast differentiation from human pluripotent stem cells and ERV-associated gene expression

Introduction

In human placental development, the trophectoderm (TE) appears in blastocysts on day 5 post-fertilization and develops after implantation into three types of trophoblast lineages: cytotrophoblast (CT), syncytiotrophoblast (ST), and extravillous trophoblast (EVT). CDX2/Cdx2 is expressed in the TE, and Cdx2 expression is upregulated by knockdown of Foxo1 in mouse ESCs. However, the significance of FOXO1 in trophoblast lineage differentiation during the early developmental period remains unclear. In this study, we examined the effect of FOXO1 inhibition on the differentiation of naive human induced pluripotent stem cells (iPSCs) into TE and trophoblast lineages.

Methods

We induced TE differentiation from naive iPSCs in the presence or absence of a FOXO1 inhibitor, and the resulting cells were subjected to trophoblast differentiation procedures without the FOXO1 inhibitor. The cells obtained in these processes were assessed for morphology, gene expression, and hCG secretion using phase-contrast microscopy, reverse transcription polymerase chain reaction (RT-PCR), quantitative RT-PCR (RT-qPCR), RNA-seq, immunochromatography, and a chemiluminescent enzyme immunoassay.

Results

In the induction of trophoblast differentiation from naive iPSCs, treatment with a FOXO1 inhibitor resulted in the enhanced expression of TE markers, CDX2 and HAND1, but conversely decreased the expression of ST markers, such as ERVW1 (Syncytin-1) and GCM1, and an EVT marker, HLA-G. The proportion of cells positive for an early TE marker TACSTD2 and negative for a late TE marker ENPEP was higher in FOXO1 inhibitor-treated cells than in non-treated cells. The expressions of ERVW1 (Syncytin-1), ERVFRD-1 (Syncytin-2), and other endogenous retrovirus (ERV)-associated genes that have been reported to be expressed in trophoblasts were suppressed in the cells obtained by differentiating the TE cells treated with FOXO1 inhibitor.

Conclusions

Treatment with a FOXO1 inhibitor during TE induction from naive iPSCs promotes early TE differentiation but hinders the progression of differentiation into ST and EVT. The suppression of ERV-associated genes may be involved in this process.

An eQTL-based approach reveals candidate regulators of LINE-1 RNA levels in lymphoblastoid cells

Long interspersed element 1 (LINE-1; L1) are a family of transposons that occupy ~17% of the human genome. Though a small number of L1 copies remain capable of autonomous transposition, the overwhelming majority of copies are degenerate and immobile. Nevertheless, both mobile and immobile L1s can exert pleiotropic effects (promoting genome instability, inflammation, or cellular senescence) on their hosts, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSD17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for regulators of transposon RNA levels.

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.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

Lifespan analysis of repeat expression reveals age-dependent upregulation of HERV-K in the neurotypical human brain

DNA repetitive sequences (or repeats) comprise over 50% of the human genome and have a crucial regulatory role, specifically regulating transcription machinery. The human brain is the tissue with the highest detectable repeat expression and dysregulations on the repeat activity are related to several neurological and neurodegenerative disorders, as repeat-derived products can stimulate a pro-inflammatory response. Even so, it is unclear how repeat expression acts on the aging neurotypical brain. Here, we leverage a large postmortem transcriptome cohort spanning the human lifespan to assess global repeat expression in the neurotypical brain. We identified 21,696 differentially expressed repeats (DERs) that varied across seven age bins (Prenatal; 0-15; 16-29; 30-39; 40-49; 50-59; 60+) across the caudate nucleus (n=271), dorsolateral prefrontal cortex (n=304), and hippocampus (n=310). Interestingly, we found that long interspersed nuclear elements and long terminal repeats (LTRs) DERs were the most abundant repeat families when comparing infants to early adolescence (0-15) with older adults (60+). Of these differentially regulated LTRs, we identified 17 shared across all brain regions, including increased expression of HERV-K-int in older adult brains (60+). Co-expression analysis from each of the three brain regions also showed repeats from the HERV subfamily were intramodular hubs in its subnetworks. While we do not observe a strong global relationship between repeat expression and age, we identified HERV-K as a repeat signature associated with the aging neurotypical brain. Our study is the first global assessment of repeat expression in the neurotypical brain.

Navigating the landscape of epitranscriptomics and host immunity

RNA modifications, also termed epitranscriptomic marks, encompass chemical alterations to individual nucleotides, including processes such as methylation and editing. These marks contribute to a wide range of biological processes, many of which are related to host immune system defense. The functions of immune-related RNA modifications can be categorized into three main groups: regulation of immunogenic RNAs, control of genes involved in innate immune response, and facilitation of adaptive immunity. Here, we provide an overview of recent research findings that elucidate the contributions of RNA modifications to each of these processes. We also discuss relevant methods for genome-wide identification of RNA modifications and their immunogenic substrates. Finally, we highlight recent advances in cancer immunotherapies that aim to reduce cancer cell viability by targeting the enzymes responsible for RNA modifications. Our presentation of these dynamic research avenues sets the stage for future investigations in this field.

Strain-specific evolution and host-specific regulation of transposable elements in the model plant symbiont Rhizophagus irregularis

Transposable elements (TEs) are repetitive DNA that can create genome structure and regulation variability. The genome of Rhizophagus irregularis, a widely studied arbuscular mycorrhizal fungus (AMF), comprises ∼50% repetitive sequences that include TEs. Despite their abundance, two-thirds of TEs remain unclassified, and their regulation among AMF life stages remains unknown. Here, we aimed to improve our understanding of TE diversity and regulation in this model species by curating repeat datasets obtained from chromosome-level assemblies and by investigating their expression across multiple conditions. Our analyses uncovered new TE superfamilies and families in this model symbiont and revealed significant differences in how these sequences evolve both within and between R. irregularis strains. With this curated TE annotation, we also found that the number of upregulated TE families in colonized roots is 4 times higher than in the extraradical mycelium, and their overall expression differs depending on the plant host. This work provides a fine-scale view of TE diversity and evolution in model plant symbionts and highlights their transcriptional dynamism and specificity during host-microbe interactions. We also provide Hidden Markov Model profiles of TE domains for future manual curation of uncharacterized sequences (https://github.com/jordana-olive/TE-manual-curation/tree/main).

DDR2 signaling and mechanosensing orchestrate neuroblastoma cell fate through different transcriptome mechanisms

The extracellular matrix (ECM) regulates carcinogenesis by interacting with cancer cells via cell surface receptors. Discoidin Domain Receptor 2 (DDR2) is a collagen-activated receptor implicated in cell survival, growth, and differentiation. Dysregulated DDR2 expression has been identified in various cancer types, making it as a promising therapeutic target. Additionally, cancer cells exhibit mechanosensing abilities, detecting changes in ECM stiffness, which is particularly important for carcinogenesis given the observed ECM stiffening in numerous cancer types. Despite these, whether collagen-activated DDR2 signaling and ECM stiffness-induced mechanosensing exert similar effects on cancer cell behavior and whether they operate through analogous mechanisms remain elusive. To address these questions, we performed bulk RNA sequencing (RNA-seq) on human SH-SY5Y neuroblastoma cells cultured on collagen-coated substrates. Our results show that DDR2 downregulation induces significant changes in the cell transcriptome, with changes in expression of 15% of the genome, specifically affecting the genes associated with cell division and differentiation. We validated the RNA-seq results by showing that DDR2 knockdown redirects the cell fate from proliferation to senescence. Like DDR2 knockdown, increasing substrate stiffness diminishes cell proliferation. Surprisingly, RNA-seq indicates that substrate stiffness has no detectable effect on the transcriptome. Furthermore, DDR2 knockdown influences cellular responses to substrate stiffness changes, highlighting a crosstalk between these two ECM-induced signaling pathways. Based on our results, we propose that the ECM could activate DDR2 signaling and mechanosensing in cancer cells to orchestrate their cell fate through distinct mechanisms, with or without involving gene expression, thus providing novel mechanistic insights into cancer progression.

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.

Identification and characterization of endogenous retroviruses upon SARS-CoV-2 infection

Endogenous retroviruses (ERVs) derived from the long terminal repeat (LTR) family of transposons constitute a significant portion of the mammalian genome, with origins tracing back to ancient viral infections. Despite comprising approximately 8% of the human genome, the specific role of ERVs in the pathogenesis of COVID-19 remains unclear. In this study, we conducted a genome-wide identification of ERVs in human peripheral blood mononuclear cells (hPBMCs) and primary lung epithelial cells from monkeys and mice, both infected and uninfected with SARS-CoV-2. We identified 405, 283, and 206 significantly up-regulated transposable elements (TEs) in hPBMCs, monkeys, and mice, respectively. This included 254, 119, 68, and 28 ERVs found in hPBMCs from severe and mild COVID-19 patients, monkeys, and transgenic mice expressing the human ACE2 receptor (hACE2) and infected with SARS-CoV-2. Furthermore, analysis using the Genomic Regions Enrichment of Annotations Tool (GREAT) revealed certain parental genomic sequences of these up-regulated ERVs in COVID-19 patients may be involved in various biological processes, including histone modification and viral replication. Of particular interest, we identified 210 ERVs specifically up-regulated in the severe COVID-19 group. The genes associated with these differentially expressed ERVs were enriched in processes such as immune response activation and histone modification. HERV1_I-int: ERV1:LTR and LTR7Y: ERV1:LTR were highlighted as potential biomarkers for evaluating the severity of COVID-19. Additionally, validation of our findings using RT-qPCR in Bone Marrow-Derived Macrophages (BMDMs) from mice infected by HSV-1 and VSV provided further support to our results. This study offers insights into the expression patterns and potential roles of ERVs following viral infection, providing a valuable resource for future studies on ERVs and their interaction with SARS-CoV-2.

Induction of Viral Mimicry Upon Loss of DHX9 and ADAR1 in Breast Cancer Cells

Detection of viral double-stranded RNA (dsRNA) is an important component of innate immunity. However, many endogenous RNAs containing double-stranded regions can be misrecognized and activate innate immunity. The IFN-inducible ADAR1-p150 suppresses dsRNA sensing, an essential function for adenosine deaminase acting on RNA 1 (ADAR1) in many cancers, including breast. Although ADAR1-p150 has been well established in this role, the functions of the constitutively expressed ADAR1-p110 isoform are less understood. We used proximity labeling to identify putative ADAR1-p110-interacting proteins in breast cancer cell lines. Of the proteins identified, the RNA helicase DHX9 was of particular interest. Knockdown of DHX9 in ADAR1-dependent cell lines caused cell death and activation of the dsRNA sensor PKR. In ADAR1-independent cell lines, combined knockdown of DHX9 and ADAR1, but neither alone, caused activation of multiple dsRNA sensing pathways leading to a viral mimicry phenotype. Together, these results reveal an important role for DHX9 in suppressing dsRNA sensing by multiple pathways.

SIGNIFICANCE

These findings implicate DHX9 as a suppressor of dsRNA sensing. In some cell lines, loss of DHX9 alone is sufficient to cause activation of dsRNA sensing pathways, while in other cell lines DHX9 functions redundantly with ADAR1 to suppress pathway activation.

A concerted increase in readthrough and intron retention drives transposon expression during aging and senescence

Aging and senescence are characterized by pervasive transcriptional dysfunction, including increased expression of transposons and introns. Our aim was to elucidate mechanisms behind this increased expression. Most transposons are found within genes and introns, with a large minority being close to genes. This raises the possibility that transcriptional readthrough and intron retention are responsible for age-related changes in transposon expression rather than expression of autonomous transposons. To test this, we compiled public RNA-seq datasets from aged human fibroblasts, replicative and drug-induced senescence in human cells, and RNA-seq from aging mice and senescent mouse cells. Indeed, our reanalysis revealed a correlation between transposons expression, intron retention, and transcriptional readthrough across samples and within samples. Both intron retention and readthrough increased with aging or cellular senescence and these transcriptional defects were more pronounced in human samples as compared to those of mice. In support of a causal connection between readthrough and transposon expression, analysis of models showing induced transcriptional readthrough confirmed that they also show elevated transposon expression. Taken together, our data suggest that elevated transposon reads during aging seen in various RNA-seq dataset are concomitant with multiple transcriptional defects. Intron retention and transcriptional readthrough are the most likely explanation for the expression of transposable elements that lack a functional promoter.

Lysine-36 of Drosophila histone H3.3 supports adult longevity

Aging is a multifactorial process that disturbs homeostasis, increases disease susceptibility, and ultimately results in death. Although the definitive set of molecular mechanisms responsible for aging remain to be discovered, epigenetic change over time is proving to be a promising piece of the puzzle. Several post-translational histone modifications have been linked to the maintenance of longevity. Here, we focus on lysine-36 of the replication-independent histone protein, H3.3 (H3.3K36). To interrogate the role of this residue in Drosophila developmental gene regulation, we generated a lysine-to-arginine mutant that blocks the activity of its cognate-modifying enzymes. We found that an H3.3BK36R mutation causes a significant reduction in adult lifespan, accompanied by dysregulation of the genomic and transcriptomic architecture. Transgenic co-expression of wild-type H3.3B completely rescues the longevity defect. Because H3.3 is known to accumulate in nondividing tissues, we carried out transcriptome profiling of young vs aged adult fly heads. The data show that loss of H3.3K36 results in age-dependent misexpression of NF-κB and other innate immune target genes, as well as defects in silencing of heterochromatin. We propose H3.3K36 maintains the postmitotic epigenomic landscape, supporting longevity by regulating both pericentric and telomeric retrotransposons and by suppressing aberrant immune signaling.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

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.

The homeobox transcription factor DUXBL controls exit from totipotency

In mice, exit from the totipotent two-cell (2C) stage embryo requires silencing of the 2C-associated transcriptional program. However, the molecular mechanisms involved in this process remain poorly understood. Here we demonstrate that the 2C-specific transcription factor double homeobox protein (DUX) mediates an essential negative feedback loop by inducing the expression of DUXBL to promote this silencing. We show that DUXBL gains accessibility to DUX-bound regions specifically upon DUX expression. Furthermore, we determine that DUXBL interacts with TRIM24 and TRIM33, members of the TRIM superfamily involved in gene silencing, and colocalizes with them in nuclear foci upon DUX expression. Importantly, DUXBL overexpression impairs 2C-associated transcription, whereas Duxbl inactivation in mouse embryonic stem cells increases DUX-dependent induction of the 2C-transcriptional program. Consequently, DUXBL deficiency in embryos results in sustained expression of 2C-associated transcripts leading to early developmental arrest. Our study identifies DUXBL as an essential regulator of totipotency exit enabling the first divergence of cell fates.

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.

MANIPULATING MITOCHONDRIAL REACTIVE OXYGEN SPECIES ALTERS SURVIVAL IN UNEXPECTED WAYS IN A DROSOPHILA 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 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, or whether they are simply 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 per se. We also find that the response to treatment with a redox-altering drug varies dramatically depending on the age at which the drug is administered, the duration of the treatment, and the genotype of the individual receiving the drug. 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.

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Activation of endogenous retroviruses and induction of viral mimicry by MEK1/2 inhibition in pancreatic cancer

While pancreatic ductal adenocarcinomas (PDACs) are addicted to KRAS-activating mutations, inhibitors of downstream KRAS effectors, such as the MEK1/2 kinase inhibitor trametinib, are devoid of therapeutic effects. However, the extensive rewiring of regulatory circuits driven by the attenuation of the KRAS pathway may induce vulnerabilities of therapeutic relevance. An in-depth molecular analysis of the transcriptional and epigenomic alterations occurring in PDAC cells in the initial hours after MEK1/2 inhibition by trametinib unveiled the induction of endogenous retroviruses (ERVs) escaping epigenetic silencing, leading to the production of double-stranded RNAs and the increased expression of interferon (IFN) genes. We tracked ERV activation to the early induction of the transcription factor ELF3, which extensively bound and activated nonsilenced retroelements and synergized with IRF1 (interferon regulatory factor 1) in the activation of IFNs and IFN-stimulated genes. Trametinib-induced viral mimicry in PDAC may be exploited in the rational design of combination therapies in immuno-oncology.

Systematic evaluation of retroviral LTRs as cis-regulatory elements in mouse embryos

In mammals, many retrotransposons are de-repressed during zygotic genome activation (ZGA). However, their functions in early development remain elusive largely due to the challenge to simultaneously manipulate thousands of retrotransposon insertions in embryos. Here, we applied CRISPR interference (CRISPRi) to perturb the long terminal repeat (LTR) MT2_Mm, a well-known ZGA and totipotency marker that exists in ∼2,667 insertions throughout the mouse genome. CRISPRi robustly perturbed 2,485 (∼93%) MT2_Mm insertions and 1,090 (∼55%) insertions of the closely related MT2C_Mm in 2-cell embryos. Remarkably, such perturbation caused downregulation of hundreds of ZGA genes and embryonic arrest mostly at the morula stage. Mechanistically, MT2 LTRs are globally enriched for open chromatin and H3K27ac and function as promoters/enhancers downstream of OBOX/DUX proteins. Thus, we not only provide direct evidence to support the functional importance of MT2 activation in development but also systematically define cis-regulatory function of MT2 in embryos by integrating functional perturbation and multi-omic analyses.

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

Chronically high stress hormone levels dysregulate sperm long noncoding RNAs and their embryonic microinjection alters development and affective behaviours

Previous studies on paternal epigenetic inheritance have shown that sperm RNAs play a role in this type of inheritance. The microinjection of sperm small noncoding RNAs into fertilised mouse oocytes induces reprogramming of the early embryo, which is thought to be responsible for the differences observed in adult phenotype. While sperm long noncoding RNAs (lncRNAs) have also been investigated in a previous study, their microinjection into fertilised oocytes did not yield conclusive results regarding their role in modulating brain development and adult behavioural phenotypes. Therefore, in the current study we sought to investigate this further. We used our previously established paternal corticosterone (stress hormone) model to assess sperm lncRNA expression using CaptureSeq, a sequencing technique that is more sensitive than the ones used in other studies in the field. Paternal corticosterone exposure led to dysregulation of sperm long noncoding RNA expression, which encompassed lncRNAs, circular RNAs and transposable element transcripts. Although they have limited functional annotation, bioinformatic approaches indicated the potential of these lncRNAs in regulating brain development and function. We then separated and isolated the sperm lncRNAs and performed microinjections into fertilised oocytes, to generate embryos with modulated lncRNA populations. We observed that the resulting adult offspring had lower body weight and altered anxiety and affective behavioural responses, demonstrating roles for lncRNAs in modulating development and brain function. This study provides novel insights into the roles of lncRNAs in epigenetic inheritance, including impacts on brain development and behaviours of relevance to affective disorders.

Genomic profiling of six human somatic histone H1 variants denotes that H1X accumulates at recently incorporated transposable elements

Histone H1, a vital component in chromatin structure, binds to linker DNA and regulates nuclear processes. We have investigated the distribution of histone H1 variants in a breast cancer cell line using ChIP-Seq. Two major groups of variants are identified: H1.2, H1.3, H1.5 and H1.0 are abundant in low GC regions (B compartment), while H1.4 and H1X preferentially localize in high GC regions (A compartment). Examining their abundance within transposable elements (TEs) reveals that H1X and H1.4 are enriched in recently-incorporated TEs (SVA and SINE-Alu), while H1.0/H1.2/H1.3/H1.5 are more abundant in older elements. Notably, H1X is particularly enriched in SVA families, while H1.4 shows the highest abundance in young AluY elements. Although low GC variants are generally enriched in LINE, LTR and DNA repeats, H1X and H1.4 are also abundant in a subset of recent LINE-L1 and LTR repeats. H1X enrichment at SVA and Alu is consistent across multiple cell lines. Further, H1X depletion leads to TE derepression, suggesting its role in maintaining TE repression. Overall, this study provides novel insights into the differential distribution of histone H1 variants among repetitive elements, highlighting the potential involvement of H1X in repressing TEs recently incorporated within the human genome.

Inhibition of Aurora Kinase Induces Endogenous Retroelements to Induce a Type I/III IFN Response via RIG-I

Type I IFN signaling is a crucial component of antiviral immunity that has been linked to promoting the efficacy of some chemotherapeutic drugs. We developed a reporter system in HCT116 cells that detects activation of the endogenous IFI27 locus, an IFN target gene. We screened a library of annotated compounds in these cells and discovered Aurora kinase inhibitors (AURKi) as strong hits. Type I IFN signaling was found to be the most enriched gene signature after AURKi treatment in HCT116, and this signature was also strongly enriched in other colorectal cancer cell lines. The ability of AURKi to activate IFN in HCT116 was dependent on MAVS and RIG-I, but independent of STING, whose signaling is deficient in these cells. MAVS dependence was recapitulated in other colorectal cancer lines with STING pathway deficiency, whereas in cells with intact STING signaling, the STING pathway was required for IFN induction by AURKi. AURKis were found to induce expression of endogenous retroviruses (ERV). These ERVs were distinct from those induced by the DNA methyltransferase inhibitors (DNMTi), which can induce IFN signaling via ERV induction, suggesting a novel mechanism of action. The antitumor effect of alisertib in mice was accompanied by an induction of IFN expression in HCT116 or CT26 tumors. CT26 tumor growth inhibition by alisertib was absent in NSG mice versus wildtype (WT) mice, and tumors from WT mice with alisertib treatment showed increased in CD8+ T-cell infiltration, suggesting that antitumor efficacy of AURKi depends, at least in part, on an intact immune response.

SIGNIFICANCE

Some cancers deactivate STING signaling to avoid consequences of DNA damage from aberrant cell division. The surprising activation of MAVS/RIG-I signaling by AURKi might represent a vulnerability in STING signaling deficient cancers.

Locus-level L1 DNA methylation profiling reveals the epigenetic and transcriptional interplay between L1s and their integration sites

Long interspersed element 1 (L1) retrotransposons are implicated in human disease and evolution. Their global activity is repressed by DNA methylation, but deciphering the regulation of individual copies has been challenging. Here, we combine short- and long-read sequencing to unveil L1 methylation heterogeneity across cell types, families, and individual loci and elucidate key principles involved. We find that the youngest primate L1 families are specifically hypomethylated in pluripotent stem cells and the placenta but not in most tumors. Locally, intronic L1 methylation is intimately associated with gene transcription. Conversely, the L1 methylation state can propagate to the proximal region up to 300 bp. This phenomenon is accompanied by the binding of specific transcription factors, which drive the expression of L1 and chimeric transcripts. Finally, L1 hypomethylation alone is typically insufficient to trigger L1 expression due to redundant silencing pathways. Our results illuminate the epigenetic and transcriptional interplay between retrotransposons and their host genome.

Recent reactivation of a pathogenicity-associated transposable element is associated with major chromosomal rearrangements in a fungal wheat pathogen

Transposable elements (TEs) are key drivers of genomic variation contributing to recent adaptation in most species. Yet, the evolutionary origins and insertion dynamics within species remain poorly understood. We recapitulate the spread of the pathogenicity-associated Styx element across five species that last diverged ∼11 000 years ago. We show that the element likely originated in the Zymoseptoria fungal pathogen genus and underwent multiple independent reactivation events. Using a global 900-genome panel of the wheat pathogen Zymoseptoria tritici, we assess Styx copy number variation and identify renewed transposition activity in Oceania and South America. We show that the element can mobilize to create additional Styx copies in a four-generation pedigree. Importantly, we find that new copies of the element are not affected by genomic defenses suggesting minimal control against the element. Styx copies are preferentially located in recombination breakpoints and likely triggered multiple types of large chromosomal rearrangements. Taken together, we establish the origin, diversification and reactivation of a highly active TE with likely major consequences for chromosomal integrity and the expression of disease.

Autonomous transposons tune their sequences to ensure somatic suppression

Transposable elements (TEs) are a major constituent of human genes, occupying approximately half of the intronic space. During pre-messenger RNA synthesis, intronic TEs are transcribed along with their host genes but rarely contribute to the final mRNA product because they are spliced out together with the intron and rapidly degraded. Paradoxically, TEs are an abundant source of RNA-processing signals through which they can create new introns1, and also functional2 or non-functional chimeric transcripts3. The rarity of these events implies the existence of a resilient splicing code that is able to suppress TE exonization without compromising host pre-mRNA processing. Here we show that SAFB proteins protect genome integrity by preventing retrotransposition of L1 elements while maintaining splicing integrity, via prevention of the exonization of previously integrated TEs. This unique dual role is possible because of L1's conserved adenosine-rich coding sequences that are bound by SAFB proteins. The suppressive activity of SAFB extends to tissue-specific, giant protein-coding cassette exons, nested genes and Tigger DNA transposons. Moreover, SAFB also suppresses LTR/ERV elements in species in which they are still active, such as mice and flies. A significant subset of splicing events suppressed by SAFB in somatic cells are activated in the testis, coinciding with low SAFB expression in postmeiotic spermatids. Reminiscent of the division of labour between innate and adaptive immune systems that fight external pathogens, our results uncover SAFB proteins as an RNA-based, pattern-guided, non-adaptive defence system against TEs in the soma, complementing the RNA-based, adaptive Piwi-interacting RNA pathway of the germline.

Single-cell multi-omics profiling of human preimplantation embryos identifies cytoskeletal defects during embryonic arrest

Human in vitro fertilized embryos exhibit low developmental capabilities, and the mechanisms that underlie embryonic arrest remain unclear. Here using a single-cell multi-omics sequencing approach, we simultaneously analysed alterations in the transcriptome, chromatin accessibility and the DNA methylome in human embryonic arrest due to unexplained reasons. Arrested embryos displayed transcriptome disorders, including a distorted microtubule cytoskeleton, increased genomic instability and impaired glycolysis, which were coordinated with multiple epigenetic reprogramming defects. We identified Aurora A kinase (AURKA) repression as a cause of embryonic arrest. Mechanistically, arrested embryos induced through AURKA inhibition resembled the reprogramming abnormalities of natural embryonic arrest in terms of the transcriptome, the DNA methylome, chromatin accessibility and H3K4me3 modifications. Mitosis-independent sequential activation of the zygotic genome in arrested embryos showed that YY1 contributed to human major zygotic genome activation. Collectively, our study decodes the reprogramming abnormalities and mechanisms of human embryonic arrest and the key regulators of zygotic genome activation.

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.

A zinc finger transcription factor enables social behaviors while controlling transposable elements and immune response in prefrontal cortex

The neurobiological origins of social behaviors are incompletely understood. Here we utilized synthetic biology approaches to reprogram the function of ZFP189, a transcription factor whose expression and function in rodent prefrontal cortex was previously demonstrated to be protective against stress-induced social deficits. We created novel synthetic ZFP189 transcription factors including ZFP189VPR, which activates the transcription of target genes and therefore exerts opposite functional control from the endogenous, transcriptionally repressive ZFP189WT. Following viral delivery of these synthetic ZFP189 transcription factors to mouse prefrontal cortex, we observe that ZFP189-mediated transcriptional control promotes mature dendritic spine morphology on transduced pyramidal neurons. Interestingly, inversion of ZFP189-mediated transcription in this brain area, achieved by viral delivery of synthetic ZFP189VPR, precipitates social behavioral deficits in terms of social interaction, motivation, and the cognition necessary for the maintenance of social hierarchy, without other observable behavioral deficits. RNA sequencing of virally manipulated prefrontal cortex tissues reveals that ZFP189 transcription factors of opposing regulatory function (ZFP189WT versus ZFP189VPR) have opposite influence on the expression of genetic transposable elements as well as genes that participate in adaptive immune functions. Collectively, this work reveals that ZFP189 function in the prefrontal cortex coordinates structural and transcriptional neuroadaptations necessary for complex social behaviors while regulating transposable element-rich regions of DNA and the expression of immune-related genes. Given the evidence for a co-evolution of social behavior and the brain immune response, we posit that ZFP189 may have evolved to augment brain transposon-associated immune function as a way of enhancing an animal's capacity for functioning in social groups.

TDP-43 chronic deficiency leads to dysregulation of transposable elements and gene expression by affecting R-loop and 5hmC crosstalk

TDP-43 is an RNA/DNA-binding protein that forms aggregates in various brain disorders. TDP-43 engages in many aspects of RNA metabolism, but its molecular roles in regulating genes and transposable elements (TEs) have not been extensively explored. Chronic TDP-43 knockdown impairs cell proliferation and cellular responses to DNA damage. At the molecular level, TDP-43 chronic deficiency affects gene expression either locally or distally by concomitantly altering the crosstalk between R-loops and 5-hydroxymethylcytosine (5hmC) in gene bodies and long-range enhancer/promoter interactions. Furthermore, TDP-43 knockdown induces substantial disease-relevant TE activation by influencing their R-loop and 5hmC homeostasis in a locus-specific manner. Together, our findings highlight the genomic roles of TDP-43 in modulating R-loop-5hmC coordination in coding genes, distal regulatory elements, and TEs, presenting a general and broad molecular mechanism underlying the contributions of proteinopathies to the etiology of neurodegenerative disorders.

Pig H3K4me3, H3K27ac, and gene expression profiles reveal reproductive tissue-specific activity of transposable elements

Regulatory sequences and transposable elements (TEs) account for a large proportion of the genomic sequences of species; however, their roles in gene transcription, especially tissue-specific expression, remain largely unknown. Pigs serve as an excellent animal model for studying genomic sequence biology due to the extensive diversity among their wild and domesticated populations. Here, we conducted an integrated analysis using H3K27ac ChIP-seq, H3K4me3 ChIP-seq, and RNA-seq data from 10 different tissues of seven fetuses and eight closely related adult pigs. We aimed to annotate the regulatory elements and TEs to elucidate their associations with histone modifications and mRNA expression across different tissues and developmental stages. Based on correlation analysis between mRNA expression and H3K27ac and H3K4me3 peak activity, results indicated that H3K27ac exhibited stronger associations with gene expression than H3K4me3. Furthermore, 1.45% of TEs overlapped with either the H3K27ac or H3K4me3 peaks, with the majority displaying tissue-specific activity. Notably, a TE subfamily (LTR4C_SS), containing binding motifs for SIX1 and SIX4, showed specific enrichment in the H3K27ac peaks of the adult and fetal ovaries. RNA-seq analysis also revealed widespread expression of TEs in the exons or promoters of genes, including 4 688 TE-containing transcripts with distinct development stage-specific and tissue-specific expression. Of note, 1 967 TE-containing transcripts were enriched in the testes. We identified a long terminal repeat (LTR), MLT1F1, acting as a testis-specific alternative promoter in SRPK2 (a cell cycle-related protein kinase) in our pig dataset. This element was also conserved in humans and mice, suggesting either an ancient integration of TEs in genes specifically expressed in the testes or parallel evolutionary patterns. Collectively, our findings demonstrate that TEs are deeply embedded in the genome and exhibit important tissue-specific biological functions, particularly in the reproductive organs.

Auto-suppression of Tet dioxygenases protects the mouse oocyte genome from oxidative demethylation

DNA cytosine methylation plays a vital role in repressing retrotransposons, and such derepression is linked with developmental failure, tumorigenesis and aging. DNA methylation patterns are formed by precisely regulated actions of DNA methylation writers (DNA methyltransferases) and erasers (TET, ten-eleven translocation dioxygenases). However, the mechanisms underlying target-specific oxidation of 5mC by TET dioxygenases remain largely unexplored. Here we show that a large low-complexity domain (LCD), located in the catalytic part of Tet enzymes, negatively regulates the dioxygenase activity. Recombinant Tet3 lacking LCD is shown to be hyperactive in converting 5mC into oxidized species in vitro. Endogenous expression of the hyperactive Tet3 mutant in mouse oocytes results in genome-wide 5mC oxidation. Notably, the occurrence of aberrant 5mC oxidation correlates with a consequent loss of the repressive histone mark H3K9me3 at ERVK retrotransposons. The erosion of both 5mC and H3K9me3 causes ERVK derepression along with upregulation of their neighboring genes, potentially leading to the impairment of oocyte development. These findings suggest that Tet dioxygenases use an intrinsic auto-regulatory mechanism to tightly regulate their enzymatic activity, thus achieving spatiotemporal specificity of methylome reprogramming, and highlight the importance of methylome integrity for development.

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.

Dynamic chromatin remodeling in cycling human endometrium at single-cell level

Estrogen-dependent proliferation followed by progesterone-dependent differentiation of the endometrium culminates in a short implantation window. We performed single-cell assay for transposase-accessible chromatin with sequencing on endometrial samples obtained across the menstrual cycle to investigate the regulation of temporal gene networks that control embryo implantation. We identify uniquely accessible chromatin regions in all major cellular constituents of the endometrium, delineate temporal patterns of coordinated chromatin remodeling in epithelial and stromal cells, and gain mechanistic insights into the emergence of a receptive state through integrated analysis of enriched transcription factor (TF) binding sites in dynamic chromatin regions, chromatin immunoprecipitation sequencing analyses, and gene expression data. We demonstrate that the implantation window coincides with pervasive cooption of transposable elements (TEs) into the regulatory chromatin landscape of decidualizing cells and expression of TE-derived transcripts in a spatially defined manner. Our data constitute a comprehensive map of the chromatin changes that control TF activities in a cycling endometrium at cellular resolution.

Lysine-36 of Drosophila histone H3.3 supports adult longevity

Aging is a multifactorial process that disturbs homeostasis, increases disease susceptibility, and ultimately results in death. Although the definitive set of molecular mechanisms responsible for aging remain to be discovered, epigenetic change over time is proving to be a promising piece of the puzzle. Several posttranslational histone modifications (PTMs) have been linked to the maintenance of longevity. Here, we focus on lysine-36 of the replication-independent histone protein, H3.3 (H3.3K36). To interrogate the role of this residue in Drosophila developmental gene regulation, we generated a lysine to arginine mutant that blocks the activity of its cognate modifying enzymes. We found that an H3.3BK36R mutation causes a significant reduction in adult lifespan, accompanied by dysregulation of the genomic and transcriptomic architecture. Transgenic co-expression of wild-type H3.3B completely rescues the longevity defect. Because H3.3 is known to accumulate in non-dividing tissues, we carried out transcriptome profiling of young vs aged adult fly heads. The data show that loss of H3.3K36 results in age-dependent misexpression of NF-κB and other innate immune target genes, as well as defects in silencing of heterochromatin. We propose H3.3K36 maintains the postmitotic epigenomic landscape, supporting longevity by regulating both pericentric and telomeric retrotransposons and by suppressing aberrant immune signaling.

Herpesviruses mimic zygotic genome activation to promote viral replication

DUX4 is a germline transcription factor and a master regulator of zygotic genome activation. During early embryogenesis, DUX4 is crucial for maternal to zygotic transition at the 2-8-cell stage in order to overcome silencing of genes and enable transcription from the zygotic genome. In adult somatic cells, DUX4 expression is silenced and its activation in adult muscle cells causes the genetic disorder Facioscapulohumeral Muscular Dystrophy (FSHD). Here we show that herpesviruses from alpha-, beta- and gamma-herpesvirus subfamilies as well as papillomaviruses actively induce DUX4 expression to promote viral transcription and replication. We demonstrate that HSV-1 immediate early proteins directly induce expression of DUX4 and its target genes including endogenous retroelements, which mimics zygotic genome activation. We further show that DUX4 directly binds to the viral genome and promotes viral transcription. DUX4 is functionally required for herpesvirus infection, since genetic depletion of DUX4 by CRISPR/Cas9 abrogates viral replication. Our results show that herpesviruses induce DUX4 expression and its downstream germline-specific genes and retroelements, thus mimicking an early embryonic-like transcriptional program that prevents epigenetic silencing of the viral genome and facilitates herpesviral gene expression.

Cellular RNA and DNA sensing pathways are essential for the dose-dependent response of human monocytes to ionizing radiation

Circulating monocytes are important players of the inflammatory response to ionizing radiation (IR). These IR-resistant immune cells migrate to radiation-damaged tissues and differentiate into macrophages that phagocytize dying cells, but also facilitate inflammation. Besides the effect of damage-associated molecular patterns, released from irradiated tissues, the inflammatory activation of monocytes and macrophages is largely dependent on IR-induced DNA damage and aberrant transcriptional activity, which may facilitate expression of type I interferons (IFN-I) and numerous inflammation-related genes. We analyzed the accumulation of dsRNA, dsDNA fragments, and RNA:DNA hybrids in the context of induction of RNA-triggered MAVS-mediated and DNA-triggered STING-mediated signaling pathways, in primary human monocytes and a monocytic cell line, THP1, in response to various doses of gamma IR. We found that exposure to lower doses (<7.5 Gy) led to the accumulation of dsRNA, along with dsDNA and RNA:DNA hybrids and activated both MAVS and STING pathway-induced gene expression and signaling activity of IFN-I. Higher doses of IR resulted in the reduced dsRNA level, degradation of RNA-sensing mediators involved in MAVS signaling and coincided with an increased accumulation of dsDNA and RNA:DNA hybrids that correlated with elevated STING signaling and NF-κB-dependent gene expression. While both pathways activate IFN-I expression, using MAVS- and STING-knockout THP1 cells, we identified differences in the spectra of interferon-stimulated genes (ISGs) that are associated with each specific signaling pathway and outlined a large group of STING signaling-associated genes. Using the RNAi technique, we found that increasing the dose of IR activates STING signaling through the DNA sensor cGAS, along with suppression of the DDX41 helicase, which is known to reduce the accumulation of RNA:DNA hybrids and thereby limit cGAS/STING signaling activity. Together, these results indicate that depending on the applied dose, IR leads to the activation of either dsRNA-induced MAVS signaling, which predominantly leads to the expression of both pro- and anti-inflammatory markers, or dsDNA-induced STING signaling that contributes to pro-inflammatory activation of the cells. While RNA:DNA hybrids boost both MAVS- and STING-mediated signaling pathways, these structures being accumulated upon high IR doses promote type I interferon expression and appear to be potent enhancers of radiation dose-dependent pro-inflammatory activation of monocytes.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Socrates: A Novel N-Ethyl-N-nitrosourea-Induced Mouse Mutant with Audiogenic Epilepsy

Epilepsy is one of the common neurological diseases that affects not only adults but also infants and children. Because epilepsy has been studied for a long time, there are several pharmacologically effective anticonvulsants, which, however, are not suitable as therapy for all patients. The genesis of epilepsy has been extensively investigated in terms of its occurrence after injury and as a concomitant disease with various brain diseases, such as tumors, ischemic events, etc. However, in the last decades, there are multiple reports that both genetic and epigenetic factors play an important role in epileptogenesis. Therefore, there is a need for further identification of genes and loci that can be associated with higher susceptibility to epileptic seizures. Use of mouse knockout models of epileptogenesis is very informative, but it has its limitations. One of them is due to the fact that complete deletion of a gene is not, in many cases, similar to human epilepsy-associated syndromes. Another approach to generating mouse models of epilepsy is N-Ethyl-N-nitrosourea (ENU)-directed mutagenesis. Recently, using this approach, we generated a novel mouse strain, soc (socrates, formerly s8-3), with epileptiform activity. Using molecular biology methods, calcium neuroimaging, and immunocytochemistry, we were able to characterize the strain. Neurons isolated from soc mutant brains retain the ability to differentiate in vitro and form a network. However, soc mutant neurons are characterized by increased spontaneous excitation activity. They also demonstrate a high degree of Ca2+ activity compared to WT neurons. Additionally, they show increased expression of NMDA receptors, decreased expression of the Ca2+-conducting GluA2 subunit of AMPA receptors, suppressed expression of phosphoinositol 3-kinase, and BK channels of the cytoplasmic membrane involved in protection against epileptogenesis. During embryonic and postnatal development, the expression of several genes encoding ion channels is downregulated in vivo, as well. Our data indicate that soc mutation causes a disruption of the excitation-inhibition balance in the brain, and it can serve as a mouse model of epilepsy.

An eQTL-based Approach Reveals Candidate Regulators of LINE-1 RNA Levels in Lymphoblastoid Cells

Long interspersed element 1 (L1) are a family of autonomous, actively mobile transposons that occupy ~17% of the human genome. A number of pleiotropic effects induced by L1 (promoting genome instability, inflammation, or cellular senescence) have been observed, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSDB17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for transposon transcriptional regulators.

Single-cell transcriptomics of human traumatic brain injury reveals activation of endogenous retroviruses in oligodendroglia

Traumatic brain injury (TBI) is a leading cause of chronic brain impairment and results in a robust, but poorly understood, neuroinflammatory response that contributes to the long-term pathology. We used single-nuclei RNA sequencing (snRNA-seq) to study transcriptomic changes in different cell populations in human brain tissue obtained acutely after severe, life-threatening TBI. This revealed a unique transcriptional response in oligodendrocyte precursors and mature oligodendrocytes, including the activation of a robust innate immune response, indicating an important role for oligodendroglia in the initiation of neuroinflammation. The activation of an innate immune response correlated with transcriptional upregulation of endogenous retroviruses in oligodendroglia. This observation was causally linked in vitro using human glial progenitors, implicating these ancient viral sequences in human neuroinflammation. In summary, this work provides insight into the initiating events of the neuroinflammatory response in TBI, which has therapeutic implications.

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.

Identifying new players in structural synaptic plasticity through dArc1 interrogation

The formation, expansion, and pruning of synapses, known as structural synaptic plasticity, is needed for learning and memory, and perturbation of plasticity is associated with many neurological disorders and diseases. Previously, we observed that the Drosophila homolog of Activity-regulated cytoskeleton-associated protein (dArc1), forms a capsid-like structure, associates with its own mRNA, and is transported across synapses. We demonstrated that this transfer is needed for structural synaptic plasticity. To identify mRNAs that are modified by dArc1 in presynaptic neuron and postsynaptic muscle, we disrupted the expression of dArc1 and performed genomic analysis with deep sequencing. We found that dArc1 affects the expression of genes involved in metabolism, phagocytosis, and RNA-splicing. Through immunoprecipitation we also identified potential mRNA cargos of dArc1 capsids. This study suggests that dArc1 acts as a master regulator of plasticity by affecting several distinct and highly conserved cellular processes.

Expression analysis suggests that DNMT3L is required for oocyte de novo DNA methylation only in Muridae and Cricetidae rodents

BACKGROUND

During early mammalian development, DNA methylation undergoes two waves of reprogramming, enabling transitions between somatic cells, oocyte and embryo. The first wave of de novo DNA methylation establishment occurs in oocytes. Its molecular mechanisms have been studied in mouse, a classical mammalian model. Current model describes DNA methyltransferase 3A (DNMT3A) and its cofactor DNMT3L as two essential factors for oocyte DNA methylation-the ablation of either leads to nearly complete abrogation of DNA methylation. However, DNMT3L is not expressed in human oocytes, suggesting that the mechanism uncovered in mouse is not universal across mammals.

RESULTS

We analysed available RNA-seq data sets from oocytes of multiple mammals, including our novel data sets of several rodent species, and revealed that Dnmt3l is expressed only in the oocytes of mouse, rat and golden hamster, and at a low level in guinea pigs. We identified a specific promoter sequence recognised by an oocyte transcription factor complex associated with strong Dnmt3l activity and demonstrated that it emerged in the rodent clade Eumuroida, comprising the families Muridae (mice, rats, gerbils) and Cricetidae (hamsters). In addition, an evolutionarily novel promoter emerged in the guinea pig, driving weak Dnmt3l expression, likely without functional relevance. Therefore, Dnmt3l is expressed and consequently plays a role in oocyte de novo DNA methylation only in a small number of rodent species, instead of being an essential pan-mammalian factor. In contrast to somatic cells, where catalytically inactive DNMT3B interacts with DNMT3A, forming a heterotetramer, we did not find evidence for the expression of such inactive Dnmt3b isoforms in the oocytes of the tested species.

CONCLUSIONS

The analysis of RNA-seq data and genomic sequences revealed that DNMT3L is likely to play a role in oocytes de novo DNA methylation only in mice, rats, gerbils and hamsters. The mechanism governing de novo DNA methylation in the oocytes of most mammalian species, including humans, occurs through a yet unknown mechanism that differs from the current model discovered in mouse.

LINE-1 retrotransposons drive human neuronal transcriptome complexity and functional diversification

The genetic mechanisms underlying the expansion in size and complexity of the human brain remain poorly understood. Long interspersed nuclear element-1 (L1) retrotransposons are a source of divergent genetic information in hominoid genomes, but their importance in physiological functions and their contribution to human brain evolution are largely unknown. Using multiomics profiling, we here demonstrate that L1 promoters are dynamically active in the developing and the adult human brain. L1s generate hundreds of developmentally regulated and cell type-specific transcripts, many that are co-opted as chimeric transcripts or regulatory RNAs. One L1-derived long noncoding RNA, LINC01876, is a human-specific transcript expressed exclusively during brain development. CRISPR interference silencing of LINC01876 results in reduced size of cerebral organoids and premature differentiation of neural progenitors, implicating L1s in human-specific developmental processes. In summary, our results demonstrate that L1-derived transcripts provide a previously undescribed layer of primate- and human-specific transcriptome complexity that contributes to the functional diversification of the human brain.

Genome-wide perturbations of A-to-I RNA editing dysregulated circular RNAs promoting the development of cervical cancer

Cervical cancer, the second most common female malignant tumor, seriously threatens women's health and lives. Despite the availability of the HPV vaccine, effective treatment options for cervical cancer are still lacking. New research perspectives now clarify that RNA editing dysregulation and changes in circRNA expression are jointly involved in disease pathogenesis, so molecular changes associated with circRNA and RNA editing may provide clues for the development of new therapeutic strategies for cervical cancer. In this study, we designed a series of pipelines to identify and analyze dysregulated RNA editing events in circRNAs. Our findings indicate a decrease in A-to-I RNA editing levels in cervical cancer compared to normal tissues, and editing may influence the back-splicing process of circRNAs through structural modifications of Alu elements. Moreover, our research reveals that RNA editing could modulate circRNA biogenesis by influencing RNA binding protein (RBP) binding on a transcriptome-wide scale, as well as influence the expression and coding potential of circRNAs. Importantly, we identified three RNA editing sites that could serve as potential biomarkers. In summary, our study presents a comprehensive landscape of RNA editing perturbations in circRNAs, providing new insights into the complex relationship between RNA editing and circRNA dysregulation in cervical cancer.