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

PUS10-induced tRNA fragmentation impacts retrotransposon-driven inflammation

Pseudouridine synthases (PUSs) catalyze the isomerization of uridine (U)-to-pseudouridine (Ψ) and have emerging roles in development and disease. How PUSs adapt gene expression under stress remains mostly unexplored. We identify an unconventional role for the Ψ "writer" PUS10 impacting intracellular innate immunity. Using Pus10 knockout mice, we uncover cell-intrinsic upregulation of interferon (IFN) signaling, conferring resistance to inflammation in vivo. Pus10 loss alters tRNA-derived small RNAs (tdRs) abundance, perturbing translation and endogenous retroelements expression. These alterations promote proinflammatory RNA-DNA hybrids accumulation, potentially activating cyclic GMP-AMP synthase (cGAS)-stimulator of interferon gene (STING). Supplementation with selected tdR pools partly rescues these effects through interactions with RNA processing factors that modulate immune responses, revealing a regulatory circuit that counteracts cell-intrinsic inflammation. By extension, we define a PUS10-specific molecular fingerprint linking its dysregulation to human autoimmune disorders, including inflammatory bowel diseases. Collectively, these findings establish PUS10 as a viral mimicry modulator, with broad implications for innate immune homeostasis and autoimmunity.

Limited impact of the siRNA pathway on transposable element expression in Aedes aegypti

BACKGROUND

Transposable elements (TEs) are DNA sequences that can change their position within a genome. In insects, small RNA pathways are central to the transcriptional and post-transcriptional regulation of TE expression. The Piwi-interacting RNA (piRNA) pathway is particularly important in germline tissues, where it silences TE transcripts via small RNAs of 24-30 nucleotides (nt) in length produced from genomic precursor transcripts as well as through a "ping-pong" amplification cycle. The small interfering RNA (siRNA) pathway helps restrict TE expression in somatic tissues via 21nt small RNAs produced from double-stranded RNA by the endonuclease Dicer2, which guide an RNA-induced silencing complex to degrade complementary RNAs. However, much of this knowledge comes from studies of the model insect Drosophila melanogaster. In the mosquito Aedes aegypti, a medically significant vector species, the siRNA pathway has mainly been investigated in connection with its antiviral role, leaving open whether it also regulates TE expression.

RESULTS

We investigated the expression of TEs and small RNAs in both somatic and gonadal tissues of a Dicer2 mutant line of Ae. aegypti and its wild-type counterpart. Our results show a modified pattern of TE expression and a decrease in TE-derived 21nt RNAs in the Dicer2 mutant, but no major shift of TE transcript abundance. The lack of a functional siRNA pathway also causes perturbations in piRNA ping-pong signatures and the expression of certain piRNA-associated genes, but without clear evidence for compensation by increased piRNA pathway activity.

CONCLUSIONS

The mosquito Ae. aegypti produces siRNAs derived from TEs but these lack a critical role in the regulation of TE expression both in somatic and in gonadal tissues.

DUX4 activates common and context-specific intergenic transcripts and isoforms

DUX4 regulates the expression of genic and nongenic elements and modulates chromatin accessibility during zygotic genome activation in cleavage stage embryos. Its misexpression in skeletal muscle causes facioscapulohumeral dystrophy (FSHD). By leveraging full-length RNA isoform sequencing with short-read RNA sequencing of DUX4-inducible myoblasts, we elucidate an isoform-resolved transcriptome featuring numerous unannotated isoforms from known loci and novel intergenic loci. While DUX4 activates similar programs in early embryos and FSHD muscle, the isoform usage of known DUX4 targets is notably distinct between the two contexts. DUX4 also activates hundreds of previously unannotated intergenic loci dominated by repetitive elements. The transcriptional and epigenetic profiles of these loci in myogenic and embryonic contexts indicate that the usage of DUX4-binding sites at these intergenic loci is influenced by the cellular environment. These findings demonstrate that DUX4 induces context-specific transcriptomic programs, enriching our understanding of DUX4-induced muscle pathology.

Proteogenomic reprogramming to a functional human blastomere-like stem cell state via a PARP-DUX4 regulatory axis

Here, we show that conventional human pluripotent stem cells cultured with non-specific tankyrase-PARP1-inhibited conditions underwent proteogenomic reprogramming to functional blastomere-like tankyrase/PARP inhibitor-regulated naive stem cells (TIRN-SC). TIRN-SCs concurrently expressed hundreds of pioneer factors in hybrid 2C-8C-morula-ICM programs that were augmented by induced expression of DUX4. Injection of TIRN-SCs into 8C-staged murine embryos equipotently differentiated human cells to the extra-embryonic and embryonic compartments of chimeric blastocysts and fetuses. Ectopic expression of murine-E-Cadherin in TIRN-SCs further enhanced interspecific chimeric tissue targeting. TIRN-SC-derived trophoblast stem cells efficiently generated placental chimeras. Proteome-ubiquitinome analyses revealed increased TNKS and reduced PARP1 levels and an ADP-ribosylation-deficient, hyper-ubiquitinated proteome that impacted expression of both tankyrase and PARP1 substrates. ChIP-seq of NANOG-SOX2-OCT4 and PARP1 (NSOP) revealed genome-wide NSOP co-binding at DUX4-accessible enhancers of embryonic lineage factors; suggesting a DUX4-NSOP axis regulated TIRN-SC lineage plasticity. TIRN-SCs may serve as valuable models for studying the proteogenomic regulation of pre-lineage human embryogenesis. VIDEO ABSTRACT.

Transposable elements as genome regulators in normal and malignant haematopoiesis

Transposable elements (TEs) constitute over half of the human genome and have played a profound role in genome evolution. While most TEs have lost the ability to transpose, many retain functional elements that serve as drivers of genome innovation, including the emergence of novel genes and regulatory elements. Recent advances in experimental and bioinformatic methods have provided new insights into their roles in human biology, both in health and disease. In this review, we discuss the multifaceted roles of TEs in haematopoiesis, highlighting their contributions to both normal and pathological contexts. TEs influence gene regulation by reshaping gene-regulatory networks, modulating transcriptional activity, and creating novel regulatory elements. These activities play key roles in maintaining normal haematopoietic processes and supporting cellular regeneration. However, in haematological malignancies, TE reactivation can disrupt genomic integrity, induce structural variations, and dysregulate transcriptional programmes, thereby driving oncogenesis. By examining the impact of TE activity on genome regulation and variation, we highlight their pivotal roles in both normal haematopoietic processes and haematological cancers.

Gag proteins encoded by endogenous retroviruses are required for zebrafish development

Transposable elements (TEs) make up the bulk of eukaryotic genomes and examples abound of TE-derived sequences repurposed for organismal function. The process by which TEs become coopted remains obscure because most cases involve ancient, transpositionally inactive elements. Reports of active TEs serving beneficial functions are scarce and often contentious due to difficulties in manipulating repetitive sequences. Here, we show that recently active TEs in zebrafish encode products critical for embryonic development. Knockdown and rescue experiments demonstrate that the endogenous retrovirus family BHIKHARI-1 (Bik-1) encodes a Gag protein essential for mesoderm development. Mechanistically, Bik-1 Gag associates with the cell membrane, and its ectopic expression in chicken embryos alters cell migration. Similarly, depletion of BHIKHARI-2 Gag, a relative of Bik-1, causes defects in neural crest development in zebrafish. We propose an "addiction" model to explain how active TEs can be integrated into conserved developmental processes.

Intermittent fasting boosts sexual behavior by limiting the central availability of tryptophan and serotonin

Aging affects reproductive capabilities in males through physiological and behavioral alterations, including endocrine changes and decreased libido. In this study, we investigated the influence of intermittent fasting (IF) on these aging-related declines, using male C57BL/6J mice. Our findings revealed that IF significantly preserved reproductive success in aged mice, not by improving traditional reproductive metrics such as sperm quality or endocrine functions but by enhancing mating behavior. This behavioral improvement was attributed to IF's ability to counter age-dependent increases in serotonergic inhibition, primarily through the decreased supply of the serotonin precursor tryptophan from the periphery to the brain. Our research underscores the potential of dietary interventions like IF in mitigating age-associated declines in male reproductive health and suggests a novel approach to managing conditions related to reduced sexual desire, highlighting the complex interplay between diet, metabolism, and reproductive behavior.

H3K27me3 and the PRC1-H2AK119ub pathway cooperatively maintain heterochromatin and transcriptional silencing after the loss of H3K9 methylation

BACKGROUND

Heterochromatin is a fundamental component of eukaryotic chromosome architecture, crucial for genome stability and cell type-specific gene regulation. In mammalian nuclei, heterochromatin forms condensed B compartments, distinct from the transcriptionally active euchromatic A compartments. Histone H3 lysine 9 and lysine 27 trimethylation (H3K9me3 and H3K27me3) are two major epigenetic modifications that enrich constitutive and facultative heterochromatin, respectively. Previously, we found that the redistribution of H3K27me3 following the loss of H3K9 methylation contributes to heterochromatin maintenance, while the simultaneous loss of both H3K27me3 and H3K9 methylation induces heterochromatin decondensation in mouse embryonic fibroblasts. However, the spatial positioning of B compartments largely persists, suggesting additional mechanisms are involved.

RESULTS

In this study, we investigated the role of H2AK119 monoubiquitylation (uH2A), a repressive chromatin mark deposited by Polycomb Repressive Complex 1 (PRC1), in maintaining heterochromatin structure following the loss of H3K9 and H3K27 methylation. We observed that uH2A and H3K27me3 are independently enriched in B compartments after H3K9 methylation loss. Despite the absence of H3K9me3 and H3K27me3, uH2A remained localized and contributed to heterochromatin retention. These results suggest that PRC1-mediated uH2A functions independently and cooperatively with H3K27me3 to maintain heterochromatin organization originally created by H3K9 methylation.

CONCLUSION

Our findings highlight a compensatory role for uH2A in preserving heterochromatin structure after the loss of other repressive chromatin modifications. The PRC1-uH2A pathway plays a critical role in maintaining the integrity of B compartments and suggests that heterochromatin architecture is supported by a network of redundant epigenetic mechanisms in mammalian cells.

α-Ketoglutarate inhibits the pluripotent-to-totipotent state transition in stem cells

In early mouse embryogenesis, the distinct enrichment of α-ketoglutarate (αKG) in blastocysts and L-2-hydroxyglutarate (L-2HG) in 2-cell (2C) embryos serves as a key metabolic signature. While elevated L-2HG levels inhibit the resolution of totipotency during the transition from the 2C stage to the blastocyst, the role of αKG remains elusive. Mouse embryonic stem cells (mESCs) cultured in vitro naturally harbor a subpopulation that transitions dynamically into a 2C-like totipotent state, providing a convenient model to investigate the role of αKG in totipotency reprogramming. This study demonstrates that αKG significantly inhibits the pluripotency to totipotency transition through upregulating ten-eleven translocation (TET) DNA hydroxylases. We further show that reducing endogenous αKG levels via glutamine withdrawal or inhibiting αKG-dependent dioxygenases by blocking succinate dehydrogenase (SDH) markedly enhances the induction of 2C-like cells (2CLCs). Finally, leveraging the potent SDH inhibitor dimethyl malonate (DMM), we have developed a highly efficient protocol for 2CLC induction, producing cells that transcriptionally resemble mid-to-late 2C embryos. Our findings deepen the understanding of the metabolic regulation of totipotency and provide a previously undescribed approach for capturing totipotent-like stem cells in vitro.

Transposable Elements Drive Regulatory and Functional Innovation of F-box Genes

Protein domains of transposable elements (TEs) and viruses increase the protein diversity of host genomes by recombining with other protein domains. By screening 10 million eukaryotic proteins, we identified several domains that define multicopy gene families and frequently co-occur with TE/viral domains. Among these, a Tc1/Mariner transposase helix-turn-helix (HTH) domain was captured by F-box genes in the Caenorhabditis genus, creating a new class of F-box genes. For specific members of this class, like fbxa-215, we found that the HTH domain is required for diverse processes including germ granule localization, fertility, and thermotolerance. Furthermore, we provide evidence that Heat Shock Factor 1 (HSF-1) mediates the transcriptional integration of fbxa-215 into the heat shock response by binding to Helitron TEs directly upstream of the fbxa-215 locus. The interactome of HTH-bearing F-box factors suggests roles in post-translational regulation and proteostasis, consistent with established functions of F-box proteins. Based on AlphaFold2 multimer proteome-wide screens, we propose that the HTH domain may diversify the repertoire of protein substrates that F-box factors regulate post-translationally. We also describe an independent capture of a TE domain by F-box genes in zebrafish. In conclusion, we identify two independent TE domain captures by F-box genes in eukaryotes and provide insights into how these novel proteins are integrated within host gene regulatory networks.

Paramutation-Like Behavior of Genic piRNA-Producing Loci in Drosophila virilis

Piwi-interacting RNAs (piRNAs) play a crucial role in silencing transposable elements (TEs) in the germ cells of Metazoa by acting as sequence-specific guides. Originating from distinct genomic loci, called piRNA clusters, piRNA can trigger an epigenetic conversion of TE insertions into piRNA clusters by means of a paramutation-like process. However, the variability in piRNA clusters' capacity to induce such conversions remains poorly understood. Here, we investigated two Drosophila virilis strains with differing capacities to produce piRNAs from the subtelomeric RhoGEF3 and Adar gene loci. We found that active piRNA generation correlates with high levels of the heterochromatic mark histone 3 lysine 9 trimethylation (H3K9me3) over genomic regions that give rise to piRNAs. Importantly, the maternal transmission of piRNAs drives their production in the progeny, even from homologous loci previously inactive in piRNA biogenesis. The RhoGEF3 locus, once epigenetically converted, maintained enhanced piRNA production in subsequent generations lacking the original allele carrying the active piRNA cluster. In contrast, piRNA expression from the converted Adar locus was lost in offspring lacking the inducer allele. The present findings suggest that the paramutation-like behavior of piRNA clusters may be influenced not only by piRNAs but also by structural features and the chromatin environment in the proximity to telomeres, providing new insights into the epigenetic regulation of the Drosophila genome.

DUX-mediated configuration of p300/CBP drives minor zygotic genome activation independent of its catalytic activity

Maternal-deposited factors initiate zygotic genome activation (ZGA), driving the maternal-to-zygotic transition; however, the coordination between maternal coactivators and transcription factors (TFs) in this process remains unclear. In this study, by profiling the dynamic landscape of p300 during mouse ZGA, we reveal its role in promoting RNA polymerase II (Pol II) pre-configuration at ZGA gene regions and sequentially establishing enhancer activity and regulatory networks. Moreover, p300/CBP-catalyzed acetylation drives Pol II elongation and minor ZGA gene expression by inducing pivotal TFs such as Dux. Remarkably, the supplementation of exogenous Dux rescues ZGA failure and developmental defects caused by the loss of p300/CBP acetylation. DUX functions as a pioneer factor, guiding p300 and Pol II to minor ZGA gene regions and activating them in a manner dependent on the non-catalytic functions of p300/CBP. Together, our findings reveal a mutual dependency between p300/CBP and DUX, highlighting their coordinated role in regulating minor ZGA activation.

The establishment of nuclear organization in mouse embryos is orchestrated by multiple epigenetic pathways

The folding of the genome in the 3D nuclear space is fundamental for regulating all DNA-related processes. The association of the genome with the nuclear lamina into lamina-associated domains (LADs) represents the earliest feature of nuclear organization during development. Here, we performed a gain-of-function screen in mouse embryos to obtain mechanistic insights. We find that perturbations impacting histone H3 modifications, heterochromatin, and histone content are crucial for the establishment of nuclear architecture in zygotes and/or 2-cell-stage embryos. Notably, some perturbations exerted differential effects on zygotes versus 2-cell-stage embryos. Moreover, embryos with disrupted LADs can rebuild nuclear architecture at the 2-cell stage, indicating that the initial establishment of LADs in zygotes might be dispensable for early development. Our findings provide valuable insights into the functional interplay between chromatin and structural components of the nucleus that guide genome-lamina interactions during the earliest developmental stages.

Transposable element dynamics in Xenopus laevis embryogenesis: a tale of two coexisting subgenomes

The African clawed frog Xenopus laevis has an allotetraploid genome consisting of two subgenomes referred as L relating to the Long chromosomes and S relating to the Short chromosomes. While the L subgenome presents conserved synteny with X. tropicalis chromosomes, the S subgenome has undergone rearrangements and deletions leading to differences in gene and transposable element (TE) content between the two subgenomes. The asymmetry in the evolution of the two subgenomes is also detectable in gene expression levels and TE mobility. TEs, also known as "jumping genes", are mobile genetic elements having a key role in genome evolution and gene regulation. However, due to their potential deleterious effects, TEs are controlled by host defense mechanisms such as the nucleosome remodeling and deacetylase (NuRD) complex and the Argonaute proteins that mainly modify the heterochromatin environment. In embryogenesis, TEs can escape the silencing mechanisms during the maternal-to-zygotic transition when a transcriptionally permissive environment is created. Moreover, further evidence highlighted that the reactivation of TEs during early developmental stages is not the result of this genome-wide reorganization of chromatin but it is class and stage-specific, suggesting a precise regulation. In line with these premises, we explored the impact of TE transcriptional contribution in six developmental stages of X. laevis. Overall, the expression pattern referred to the entire set of transcribed TEs was constant across the six developmental stages and in line with their abundance in the genome. However, focusing on subgenome-specific TEs, our analyses revealed a distinctive transcriptional pattern dominated by LTR retroelements in the L subgenome and LINE retroelements in the S subgenome attributable to young copies. Interestingly, genes encoding proteins involved in maintaining the repressive chromatin environment were active in both subgenomes highlighting that TE controlling systems were active in X. laevis embryogenesis and evolved symmetrically.

Chronic Viral Mimicry Induction following p53 Loss Promotes Immune Evasion

SIGNIFICANCE

Our landmark discovery of viral mimicry characterized repetitive elements as immunogenic stimuli that cull cancer cells. If expressed repetitive elements cull cancer cells, why does every human cancer express repetitive elements? Our report offers an exciting advancement toward understanding this paradox and how to exploit this mechanism for cancer interception. See related commentary by Murayama and Cañadas, p. 670.

Environmentally induced variation in sperm sRNAs is linked to gene expression and transposable elements in zebrafish offspring

Environmental factors affect not only paternal condition but may translate into the following generations where sperm-mediated small RNAs (sRNAs) can contribute to the transmission of paternal effects. sRNAs play a key role in the male germ line in genome maintenance and repair, and particularly in response to environmental stress and the resulting increase in transposable element (TE) activity. Here, we investigated how the social environment (high competition, low competition) of male zebrafish Danio rerio affects sRNAs in sperm and how these are linked to gene expression and TE activity in their offspring. In a first experiment, we collected sperm samples after exposing males to each social environment for 2 weeks to test for differentially expressed sperm micro- (miRNA) and piwi-interacting RNAs (piRNA). In a separate experiment, we performed in vitro fertilisations after one 2-week period using a split-clutch design to control for maternal effects and collected embryos at 24 h to test for differentially expressed genes and TEs. We developed new computational prediction tools to link sperm sRNAs with differentially expressed TEs and genes in the embryos. Our results support the idea that the molecular stress response in the male germ line has significant down-stream effects on the molecular pathways, and we provide a direct link between sRNAs, TEs and gene expression.

The transcription factor SRF regulates MERVL retrotransposons and gene expression during zygotic genome activation

The regulatory circuitry of cell-specific transcriptional programs is thought to be influenced by transposable elements (TEs), whereby TEs serve as raw material for the diversification and genome-wide distribution of genetic elements that contain cis-regulatory activity. However, the transcriptional activators of TEs in relevant physiological contexts are largely unknown. Here, we undertook an evolutionary approach to identify regulators of two main families of MERVL, a major regulator of transcription during early mouse development. Using a combination of phyloregulatory, transcriptomic, and loss-of-function approaches, we demonstrate that SRF is a novel regulator of MERVL and embryonic transcription during zygotic genome activation. By resolving the phylogenetic history of two major MERVL families, we delineate the evolutionary acquisition of SRF and DUX binding sites and show that the acquisition of the SRF site precedes that of DUX. SRF contributes to embryonic transcription through the regulation of MERVLs, which in turn serve as promoters for host genes. Our work identifies new transcriptional regulators and TEs that shape the gene expression programs in early embryos and highlights the process of TE domestication via the sequential acquisition of transcription factor binding sites and coevolution with the host.

Chromatin accessibility landscape of mouse early embryos revealed by single-cell NanoATAC-seq2

In mammals, fertilized eggs undergo genome-wide epigenetic reprogramming to generate the organism. However, our understanding of epigenetic dynamics during preimplantation development at single-cell resolution remains incomplete. Here, we developed scNanoATAC-seq2, a single-cell assay for transposase-accessible chromatin using long-read sequencing for scarce samples. We present a detailed chromatin accessibility landscape of mouse preimplantation development, revealing distinct chromatin signatures in the epiblast, primitive endoderm, and trophectoderm during lineage segregation. Differences between zygotes and two-cell embryos highlight reprogramming in chromatin accessibility during the maternal-to-zygotic transition. Single-cell long-read sequencing enables in-depth analysis of chromatin accessibility in noncanonical imprinting, imprinted X chromosome inactivation, and low-mappability genomic regions, such as repetitive elements and paralogs. Our data provide insights into chromatin dynamics during mammalian preimplantation development and lineage differentiation.

Survey of gene, lncRNA and transposon transcription patterns in four mouse organs highlights shared and organ-specific sex-biased regulation

BACKGROUND

Sex-biased gene regulation is the basis of sexual dimorphism in phenotypes and has been studied across different cell types and different developmental stages. However, sex-biased expression of transposable elements (TEs), which represent nearly half of the mammalian genome and have the potential of influencing genome integrity and regulation, remains underexplored.

RESULTS

We report a survey of gene, lncRNA, and TE expression in four organs from mice with different combinations of gonadal and genetic sex. The data show remarkable variability among organs with respect to the impact of gonadal sex on transcription with the strongest effects observed in the liver. In contrast, the X-chromosome dosage alone had a modest influence on sex-biased transcription across organs, albeit interaction between X-dosage and gonadal sex cannot be ruled out. The presence of the Y-chromosome influences TE, but not gene or lncRNA, expression in the liver. Notably, 90% of sex-biased TEs (sDETEs) reside in clusters. Moreover, 54% of these clusters overlap or reside less than 100 kb from sex-biased genes or lncRNAs, share the same sex bias, and also have higher expression levels than sDETE clusters that do not co-localize with other types of sex-biased transcripts. We test the heterochromatic sink hypothesis that predicts higher expression of TEs in XX individuals finding no evidence to support it.

CONCLUSIONS

Our data show that sex-biased expression of TEs varies among organs with the highest numbers of sDETEs found in the liver following trends observed for genes and lncRNAs. It is enhanced by proximity to other types of sex-biased transcripts.

LocusMasterTE: integrating long-read RNA sequencing improves locus-specific quantification of transposable element expression

Transposable elements (TEs) can influence human diseases by disrupting genome integrity, yet their quantification has been challenging due to the repetitive nature of these sequences across the genome. We develop LocusMasterTE, a method that integrates long-read with short-read RNA-seq to increase the accuracy of TE expression quantification. By incorporating fractional transcript per million values from long-read sequencing data into an expectation-maximization algorithm, LocusMasterTE reassigns multi-mapped reads, enhancing accuracy in short-read-based TE quantification. We validate the method with simulated and human datasets. LocusMasterTE may give new insights into TE functions through precise quantification.

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

Amyotrophic lateral sclerosis (ALS) is a complex syndrome with multiple genetic causes and wide variation in disease presentation. Despite this heterogeneity, large-scale genomics studies revealed that ALS postmortem samples can be grouped into a small number of subtypes, defined by transcriptomic signatures of mitochondrial dysfunction and oxidative stress (ALS-Ox), microglial activation and neuroinflammation (ALS-Glia), or TDP-43 pathology and associated transposable elements (ALS-TE). In this study, we present a deep ALS neural net classifier (DANCer) for ALS molecular subtypes. Applying DANCer to an expanded cohort from the NYGC ALS Consortium highlights two subtypes that strongly correlate with disease duration: ALS-TE in cortex and ALS-Glia in spinal cord. Finally, single-nucleus transcriptomes demonstrate that ALS subtypes are recapitulated in neurons and glia, with both ALS-wide and subtype-specific alterations in all cell types. In summary, ALS molecular subtypes represent a combination of cellular and pathological features that correlate with clinical features of ALS.

Multi-omics analysis reveals critical cis-regulatory roles of transposable elements in livestock genomes

Transposable elements (TEs) are important sources of genetic and regulatory variation, yet their functional roles in domesticated animals remain insufficiently explored. To address this gap, we comprehensively annotated TE types, ages, and distributions in the genomes of pig (Sus scrofa), cattle (Bos taurus), and chicken (Gallus gallus). Our analysis revealed species-specific patterns in TE abundance, amplification, and activity in modern genomes. By integrating transcriptomic and epigenomic data, we explored the impact of specific TE types on cis-regulatory elements (CREs) and constructed a TE expression atlas across five tissues in all three species. Our findings underscored the critical roles of tissue-specific TE expression and chromatin accessibility in regulating tissue-specific biological processes. Most notably, we developed a computational framework to uncover TE-mediated gene regulatory networks (TE-GRNs). Our findings provide valuable insights into the regulatory functions of TEs in livestock and offer a robust approach for studying TE-GRNs in diverse biological contexts.

An evolutionarily conserved constellation of functional cis-elements programs the virus-responsive fate of the human (epi)genome

Human health depends on perplexing defensive cellular responses against microbial pathogens like Viruses. Despite the major effort undertaken, the (epi)genomic mechanisms that human cells utilize to tailor defensive gene expression programs against microbial attacks have remained inadequately understood, mainly due to a significant lack of recording of the in vivo functional cis-regulatory modules (CRMs) of the human genome. Here, we introduce the virus-responsive fate of the human (epi)genome as characterized in naïve and infected cells by functional genomics, computational biology, DNA evolution, and DNA Grammar and Syntax investigations. We discovered that multitudes of novel functional virus-responsive CRMs (vrCRMs) compose typical enhancers (tEs), super-enhancers (SEs), repetitive-DNA enhancers (rDEs), and stand-alone functional genomic stretches that grant human cells regulatory underpinnings for layering basal immunity and eliminating illogical/harmful defensive responses under homeostasis, yet stimulating virus-responsive genes and transposable elements (TEs) upon infection. Moreover, extensive epigenomic reprogramming of previously unknown SE landscapes marks the transition from naïve to antiviral human cell states and involves the functions of the antimicrobial transcription factors (TFs), including interferon response factor 3 (IRF3) and nuclear factor-κB (NF-κB), as well as coactivators and transcriptional apparatus, along with intensive modifications/alterations in histone marks and chromatin accessibility. Considering the polyphyletic evolutionary fingerprints of the composite DNA sequences of the vrCRMs assessed by TFs-STARR-seq, ranging from the animal to microbial kingdoms, the conserved features of antimicrobial TFs and chromatin complexes, and their pluripotent stimulus-induced activation, these findings shed light on how mammalian (epi)genomes evolved their functions to interpret the exogenous stress inflicted and program defensive transcriptional responses against microbial agents. Crucially, many known human short variants, e.g. single-nucleotide polymorphisms (SNPs), insertions, deletions etc., and quantitative trait loci (QTLs) linked to autoimmune diseases, such as multiple sclerosis (MS), systemic lupus erythematosus (SLE), Crohn's disease (CD) etc., were mapped within or vastly proximal (±2.5 kb) to the novel in vivo functional SEs and vrCRMs discovered, thus underscoring the impact of their (mal)functions on human physiology and disease development. Hence, we delved into the virus-responsive fate of the human (epi)genome and illuminated its architecture, function, evolutionary origins, and its significance for cellular homeostasis. These results allow us to chart the "Human hyper-Atlas of virus-infection", an integrated "molecular in silico" encyclopedia situated in the UCSC Genome Browser that benefits our mechanistic understanding of human infectious/(auto)immune diseases development and can facilitate the generation of in vivo preclinical animal models, drug design, and evolution of therapeutic applications.

The domesticated transposon protein L1TD1 associates with its ancestor L1 ORF1p to promote LINE-1 retrotransposition

Repression of retrotransposition is crucial for the successful fitness of a mammalian organism. The domesticated transposon protein L1TD1, derived from LINE-1 (L1) ORF1p, is an RNA-binding protein that is expressed only in some cancers and early embryogenesis. In human embryonic stem cells, it is found to be essential for maintaining pluripotency. In cancer, L1TD1 expression is highly correlative with malignancy progression and as such considered a potential prognostic factor for tumors. However, its molecular role in cancer remains largely unknown. Our findings reveal that DNA hypomethylation induces the expression of L1TD1 in HAP1 human tumor cells. L1TD1 depletion significantly modulates both the proteome and transcriptome and thereby reduces cell viability. Notably, L1TD1 associates with L1 transcripts and interacts with L1 ORF1p protein, thereby facilitating L1 retrotransposition. Our data suggest that L1TD1 collaborates with its ancestral L1 ORF1p as an RNA chaperone, ensuring the efficient retrotransposition of L1 retrotransposons, rather than directly impacting the abundance of L1TD1 targets. In this way, L1TD1 might have an important role not only during early development but also in tumorigenesis.

Multi-tissue transcriptomic characterization of endogenous retrovirus-derived transcripts in Capra hircus

Background

Transposable elements (TEs, or transposons) are repetitive genomic sequences, accounting for half of a mammal genome. Most TEs are transcriptionally silenced, whereas some TEs, especially endogenous retroviruses (ERVs, long terminal repeat retrotransposons), are physiologically expressed in certain conditions. However, the expression pattern of TEs in those less studied species, like goat (Capra hircus), remains unclear. To obtain an overview of the genomic and transcriptomic features of TEs and ERVs in goat, an important farm species, we herein analyzed transcriptomes of ten C. hircus tissues and cells under various physiological and pathological conditions.

Method

Distribution of classes, families, and subfamilies of TEs in the C. hircus genome were systematically annotated. The expression patterns of TE-derived transcripts in multiple tissues were investigated at subfamily and location levels. Differential expression of ERV-derived reads was measured under various physiological and pathological conditions, such as embryo development and virus infection challenges. Co-expression between ERV-reads and their proximal genes was also explored to decipher the expression regulation of ERV-derived transcripts.

Results

There are around 800 TE subfamilies in the goat genome, accounting for 49.1% of the goat genome sequence. TE-derived reads account for 10% of the transcriptome and their abundance are comparable in various goat tissues, while expression of ERVs are variable among tissues. We further characterized expression pattern of ERV reads in various tissues. Differential expression analysis showed that ERVs are highly active in 16-cell embryos, when the genome of the zygote begins to transcribe its own genes. We also recognized numerous activated ERV reads in response to RNA virus infection in lung, spleen, caecum, and immune cells. CapAeg_1.233:ERVK in chromosome 1 and 17 are dysregulated under endometrium development and infection conditions. They showed strong co-expression with their proximal gene OAS1 and TMPRSS2, indicating the impact of activated proximal gene expression on nearby ERVs.

Conclusion

We generated ERV transcriptomes across goat tissues, and identified ERVs activated in response to different physiological and pathological conditions.

Short-term and long-term high-fat diet promote metabolic disorder through reprogramming mRNA m6A in white adipose tissue by gut microbiota

BACKGROUND

Although short-term high-fat diet (S-HFD) and long-term high-fat diet (L-HFD) induce metabolic disorder, the underlying epigenetic mechanism is still unclear.

RESULTS

Here, we found that both 4 days of S-HFD and 10 weeks of L-HFD increased mRNA m6A level in epididymal white adipose tissue (eWAT) and impaired metabolic health. Interestingly, S-HFD activated transposable elements (TEs), especially endogenous retroviruses (ERVs) in eWAT, while L-HFD activated long interspersed elements (LINEs). Subsequently, we demonstrated that both S-HFD and L-HFD increased m6A level of Ehmt2 and decreased EHMT2 protein expression and H3K9me2 level, accounting for activation of ERVs and LINEs. Overexpression of EHMT2 in eWAT or inhibition of ERVs and LINEs by antiviral therapy improved metabolic health under HFD feeding. Notably, we found that both short-term and long-term HFD feeding increased Fimicutes/Bacteroidota ratio and decreased the gut microbiome health index. Fecal microbiota transplantation (FMT) experiments demonstrated that gut microbiota from S-HFD and L-HFD was responsible for increased m6A level in eWAT, resulting in glucose intolerance and insulin insensitivity. Furthermore, we identified that both S-HFD and L-HFD increased the abundance of the gut microbial metabolite homogentisic acid (HGA), and HGA level was positively correlated with unclassified_f__Lachnospiraceae which was both increased in S-HFD and L-HFD feeding mice. Administration of HGA increased the m6A level of Ehmt2 and decreased the EHMT2 protein expression and H3K9me2 level in eWAT, leading to metabolic disorder in mice.

CONCLUSIONS

Together, this study reveals a novel mechanism that S-HFD and L-HFD induce metabolism disorder through gut microbiota-HGA-m6A-Ehmt2-ERV/LINE signaling. These findings may provide a novel insight for prevention and treatment of metabolism disorder upon short-term or long-term dietary fat intake. Video Abstract.

Activating antiviral immune responses potentiates immune checkpoint inhibition in glioblastoma models

Viral mimicry refers to the activation of innate antiviral immune responses due to the induction of endogenous retroelements (REs). Viral mimicry augments antitumor immune responses and sensitizes solid tumors to immunotherapy. Here, we found that targeting what we believe to be a novel, master epigenetic regulator, Zinc Finger Protein 638 (ZNF638), induces viral mimicry in glioblastoma (GBM) preclinical models and potentiates immune checkpoint inhibition (ICI). ZNF638 recruits the HUSH complex, which precipitates repressive H3K9me3 marks on endogenous REs. In GBM, ZNF638 is associated with marked locoregional immunosuppressive transcriptional signatures, reduced endogenous RE expression, and poor immune cell infiltration. Targeting ZNF638 decreased H3K9 trimethylation, increased REs, and activated intracellular dsRNA signaling cascades. Furthermore, ZNF638 knockdown upregulated antiviral immune programs and significantly increased PD-L1 immune checkpoint expression in diverse GBM models. Importantly, targeting ZNF638 sensitized mice to ICI in syngeneic murine orthotopic models through innate IFN signaling. This response was recapitulated in recurrent GBM (rGBM) samples with radiographic responses to checkpoint inhibition with widely increased expression of dsRNA, PD-L1, and perivascular CD8 cell infiltration, suggesting that dsRNA signaling may mediate response to immunotherapy. Finally, low ZNF638 expression was a biomarker of clinical response to ICI and improved survival in patients with rGBM and patients with melanoma. Our findings suggest that ZNF638 could serve as a target to potentiate immunotherapy in gliomas.

Contrasting distributions and expression characteristics of transcribing repeats in Setaria viridis

Repetitive DNA contributes significantly to plant genome size, adaptation, and evolution. However, little is understood about the transcription of repeats. This is addressed here in the plant green foxtail millet (Setaria viridis). First, we used RepeatExplorer2 to calculate the genome proportion (GP) of all repeat types and compared the GP of long terminal repeat (LTR) retroelements against annotated complete and incomplete LTR retroelements (Ty1/copia and Ty3/gypsy) identified by DANTE in a whole genome assembly. We show that DANTE-identified LTR retroelements can comprise ∼0.75% of the inflorescence poly-A transcriptome and ∼0.24% of the stem ribo-depleted transcriptome. In the RNA libraries from inflorescence tissue, both LTR retroelements and DNA transposons identified by RepeatExplorer2 were highly abundant, where they may be taking advantage of the reduced epigenetic silencing in the germ line to amplify. Typically, there was a higher representation of DANTE-identified LTR retroelements in the transcriptome than RepeatExplorer2-identified LTR retroelements, potentially reflecting the transcription of elements that have insufficient genomic copy numbers to be detected by RepeatExplorer2. In contrast, for ribo-depleted libraries of stem tissues, the reverse was observed, with a higher transcriptome representation of RepeatExplorer2-identified LTR retroelements. For RepeatExplorer2-identified repeats, we show that the GP of most Ty1/copia and Ty3/gypsy families were positively correlated with their transcript proportion. In addition, guanine- and cytosine-rich repeats with high sequence similarity were also the most abundant in the transcriptome, and these likely represent young elements that are most capable of amplification due to their ability to evade epigenetic silencing.

Transposable Element Landscape in the Monotypic Species Barthea barthei (Hance) Krass (Melastomataceae) and Its Role in Ecological Adaptation

Transposable elements (TEs) are crucial for genome evolution and ecological adaptation, but their dynamics in non-model plants are poorly understood. Using genomic, transcriptomic, and population genomic approaches, we analyzed the TE landscape of Barthea barthei (Melastomataceae), a species distributed across tropical and subtropical southern China. We identified 64,866 TE copies (16.76% of a 235 Mb genome), dominated by Ty3/Gypsy retrotransposons (8.82%) and DNA/Mutator elements (2.7%). A genome-wide analysis revealed 13 TE islands enriched in genes related to photosynthesis, tryptophan metabolism, and stress response. We identified 3859 high-confidence TE insertion polymorphisms (TIPs), including 29 fixed insertions between red and white flower ecotypes, affecting genes involved in cell wall modification, stress response, and secondary metabolism. A transcriptome analysis of the flower buds identified 343 differentially expressed TEs between the ecotypes, 30 of which were near or within differentially expressed genes. The non-random distribution (primarily within 5 kb of genes) and association with adaptive traits suggest a significant role in B. barthei's successful colonization of diverse habitats. Our findings provide insights into how TEs contribute to plant genome evolution and ecological adaptation in tropical forests, particularly through their influence on regulatory networks governing stress response and development.

An atlas of transcription initiation reveals regulatory principles of gene and transposable element expression in early mammalian development

Transcriptional activation of the embryonic genome (EGA) is a major developmental landmark enabling the embryo to become independent from maternal control. The magnitude and control of transcriptional reprogramming during this event across mammals remains poorly understood. Here, we developed Smart-seq+5' for high sensitivity, full-length transcript coverage and simultaneous capture of 5' transcript information from single cells and single embryos. Using Smart-seq+5', we profiled 34 developmental stages in 5 mammalian species and provide an extensive characterization of the transcriptional repertoire of early development before, during, and after EGA. We demonstrate widespread transposable element (TE)-driven transcription across species, including, remarkably, of DNA transposons. We identify 19,657 TE-driven genic transcripts, suggesting extensive TE co-option in early development over evolutionary timescales. TEs display similar expression dynamics across species and species-specific patterns, suggesting shared and divergent regulation. Our work provides a powerful resource for understanding transcriptional regulation of mammalian development.

Multi-ancestry GWAS reveals loci linked to human variation in LINE-1- and Alu-insertion numbers

LINE-1 (L1) and Alu are two families of transposable elements (TEs) occupying ~17% and ~11% of the human genome, respectively. Though only a small fraction of L1 copies is able to produce the machinery to mobilize autonomously, Alu and degenerate L1s can hijack their functional machinery and mobilize in trans. The expression and subsequent mobilization of L1 and Alu can exert pathological effects on their hosts. These features have made them promising focus subjects in studies of aging where they can become active. However, mechanisms regulating TE activity are incompletely characterized, especially in diverse human populations. To address these gaps, we leveraged genomic data from the 1000 Genomes Project to carry out a trans-ethnic GWAS of L1/Alu insertion singletons. These are rare, recently acquired insertions observed in only one person and which we used as proxies for variation in L1/Alu insertion numbers. Our approach identified SNVs in genomic regions containing genes with potential and known TE regulatory properties, and it enriched for SNVs in regions containing known regulators of L1 expression. Moreover, we identified reference TE copies and structural variants that associated with L1/Alu singletons, suggesting their potential contribution to TE insertion number variation. Finally, a transcriptional analysis of lymphoblastoid cells highlighted potential cell cycle alterations in a subset of samples harboring L1/Alu singletons. Collectively, our results suggest that known TE regulatory mechanisms may be active in diverse human populations, expand the list of loci implicated in TE insertion number variability, and reinforce links between TEs and disease.

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

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

Heterochromatin fidelity is a therapeutic vulnerability in lymphoma and other human cancers

Genes involved in the regulation of chromatin structure are frequently disrupted in cancer, contributing to an aberrant transcriptome and phenotypic plasticity. Yet, therapeutics targeting mutant forms of chromatin-modifying enzymes have yielded only modest clinical utility, underscoring the difficulty of targeting the epigenomic underpinnings of aberrant gene regulatory networks. Here, we sought to identify novel epigenetic vulnerabilities in diffuse large B-cell lymphoma (DLBCL). Through phenotypic screens and biochemical analysis, we demonstrated that inhibition of the H3K9 demethylases KDM4A and KDM4C elicits potent, subtype-agnostic cytotoxicity by antagonizing transcriptional networks associated with B-cell identity and epigenetically rewiring heterochromatin. KDM4 demethylases associated with the KRAB zinc finger ZNF587, and their enzymatic inhibition led to DNA replication stress and DNA damage-einduced cGAS-STING activation. Broad surveys of transcriptional data from patients also revealed KDM4 family dysregulation in several other cancer types. To explore this potential therapeutic avenue, we performed high-throughput small molecule screens with H3K9me3 nucleosome substrates and identified novel KDM4 demethylase inhibitors. AI-guided protein-ligand binding predictions suggested diverse modes of action for various small molecule hits. Our findings underscore the relevance of targeting fundamental transcriptional and epigenetic mechanisms for anti-cancer therapy.

Capsid transfer of the retrotransposon Copia controls structural synaptic plasticity in Drosophila

Transposons are parasitic genome elements that can also serve as raw material for the evolution of new cellular functions. However, how retrotransposons are selected and domesticated by host organisms to modulate synaptic plasticity remains largely unknown. Here, we show that the Ty1 retrotransposon Copia forms virus-like capsids in vivo and transfers between cells. Copia is enriched at the Drosophila neuromuscular junction (NMJ) and transported across synapses, and disrupting its expression promotes both synapse development and structural synaptic plasticity. We show that proper synaptic plasticity is maintained in Drosophila by the balance of Copia and the Arc1 (activity-regulated cytoskeleton-associated protein) homolog. High-resolution cryogenic-electron microscopy imaging shows that the structure of the Copia capsid has a large capacity and pores like retroviruses but is distinct from domesticated capsids such as dArc1. Our results suggest a fully functional transposon mediates synaptic plasticity, possibly representing an early stage of domestication of a retrotransposon.

FishPi: a bioinformatic prediction tool to link piRNA and transposable elements

BACKGROUND

Piwi-interacting RNAs (piRNA)s are non-coding small RNAs that post-transcriptionally affect gene expression and regulation. Through complementary seed region binding with transposable elements (TEs), piRNAs protect the genome from transposition. A tool to link piRNAs with complementary TE targets will improve our understanding of the role of piRNAs in genome maintenance and gene regulation. Existing tools such as TEsmall can process sRNA-seq datasets to produce differentially expressed piRNAs, and piRScan developed for nematodes can link piRNAs and TEs but it requires knowledge about the target region of interest and works backwards.

RESULTS

We developed FishPi to predict the pairings between piRNA and TEs for available genomes from zebrafish, medaka and tilapia, with full user customisation of parameters including orientation of piRNA, mismatches in the piRNA seed binding to TE and scored output lists of piRNA-TE matches. FishPi works with individual piRNAs or a list of piRNA sequences in fasta format. The software focuses on the piRNA-TE seed region and analyses reference TEs for piRNA complementarity. TE type is examined, counted and stored to a dictionary, with genomic loci recorded. Any updates to piRNA-TE binding rules can easily be incorporated by changing the seed-region options in the graphic user-interface. FishPi provides a graphic interface using tkinter for the user to input piRNA sequences to generate comprehensive reports on piRNA-TE interactions. FishPi can easily be adapted to genomes from other species and taxa opening the interpretation of piRNA functionality to a wide community.

CONCLUSIONS

Users will gain insight into genome mobility and FishPi will help further our understanding of the biological role of piRNAs and their interaction with TEs in a similar way that public databases have improved the access to and the understanding of the role of small RNAs.

Characterization of novel human endogenous retrovirus structures on chromosomes 6 and 7

Human endogenous retroviruses (HERV) represent nearly 8% of the human genome. Of these, HERV-K subtype HML-2 is a transposable element that plays a critical role in embryonic development and in the pathogenesis of several diseases. Quantification and characterization of these multiple HML-2 insertions in the human chromosome has been challenging due to their size, sequence homology with each other, and their repetitive nature. We examined a cohort of 222 individuals for HML-2 proviruses 6q14.1 and 7p22.1a, two loci that are capable of producing full-length viral proteins and have been previously implicated in several cancers, autoimmune disorders and neurodegenerative diseases, using long-read DNA sequencing. While the reference genome for both regions suggests these two loci are structurally dissimilar, we found that for both loci about 5% of individuals have a unique tandem repeat-like sequence (three long terminal repeat sequences sandwiching two internal, potentially protein coding sequences), while most individuals have a standard proviral structure (one internal region sandwiched by two long terminal repeats). Moreover, both proviruses can make full-length, or nearly full-length, HERV-K proteins in multiple transcription orientations. The amino acid sequences from different loci in the same transcriptional orientation share sequence homology with each other. These results demonstrate a clear, previously unreported, relationship between HML-2 loci 6q14.1 and 7p22.1a and highlight the utility of long-read sequencing to study repetitive elements. Future studies need to determine if these polymorphisms determine genetic susceptibility to diseases that are associated with them.

Multi-ancestry GWAS reveals loci linked to human variation in LINE-1- and Alu-insertion numbers

LINE-1 (L1) and Alu are two families of transposable elements (TEs) occupying ~17% and ~11% of the human genome, respectively. Though only a small fraction of L1 copies is able to produce the machinery to mobilize autonomously, Alu and degenerate L1s can hijack their functional machinery and mobilize in trans. The expression and subsequent mobilization of L1 and Alu can exert pathological effects on their hosts. These features have made them promising focus subjects in studies of aging where they can become active. However, mechanisms regulating TE activity are incompletely characterized, especially in diverse human populations. To address these gaps, we leveraged genomic data from the 1000 Genomes Project to carry out a trans-ethnic GWAS of L1/Alu insertion singletons. These are rare, recently acquired insertions observed in only one person and which we used as proxies for variation in L1/Alu insertion numbers. Our approach identified SNVs in genomic regions containing genes with potential and known TE regulatory properties, and it enriched for SNVs in regions containing known regulators of L1 expression. Moreover, we identified reference TE copies and structural variants that associated with L1/Alu singletons, suggesting their potential contribution to TE insertion number variation. Finally, a transcriptional analysis of lymphoblastoid cells highlighted potential cell cycle alterations in a subset of samples harboring L1/Alu singletons. Collectively, our results suggest that known TE regulatory mechanisms may be active in diverse human populations, expand the list of loci implicated in TE insertion number variability, and reinforce links between TEs and disease.

Transposable Elements are Dysregulated in Brains of Individuals with Major Depressive Disorder

Transposable elements (TEs) are repetitive DNA sequences capable of being transcribed and reintegrated, or transposed, into distinct loci throughout the genome. While thought to be largely transcriptionally silenced in brain, TE transcription is increasingly recognized as dynamic and involved in human health and disease states, including in disorders of the brain. In this study, we annotated TE transcripts in publicly available RNA sequencing (RNAseq) of postmortem human brain tissue to investigate the expression profile of TE transcripts in individuals with Major Depressive Disorder (MDD) compared to healthy controls. Our findings reveal a robust and uniform downregulation of TE transcript expression in the brains of subjects with MDD relative to controls, this occurs most prominently in the orbitofrontal cortex (OFC) brain region, and MDD differentially impacts this TE expression by age and sex. This work points to the aberrant transcription of cortical TEs as a potentially overlooked molecular signature of MDD.

ERH Enables Early Embryonic Differentiation and Overlays H3K9me3 Heterochromatin on a Cryptic Pluripotency H3K9me3 Landscape in Somatic Cells

Enhancer of Rudimentary Homolog (ERH) is an evolutionarily conserved protein originally characterized in fission yeast 1 and recently shown to maintain H3K9me3 in human fibroblasts 2 . Here, we find that ERH depletion in fibroblasts reverts the H3K9me3 landscape to an embryonic stem cell (ESC) state and enables activation of naïve and pluripotency genes and transposable elements during induced pluripotent stem cell (iPSC) reprogramming. We find that ERH similarly represses totipotent and alternative lineage programs during mouse preimplantation development and is required for proper segregation of the inner cell mass and trophectoderm cell lineages. During human ESC differentiation into germ layer lineages, ERH silences naïve and pluripotency genes, transposable elements, and alternative lineage somatic genes. As in fission yeast, we find that mammalian ERH interacts with RNA-binding proteins to engage and repress its chromatin targets. Our findings reveal a fundamental role for ERH in cell fate specification via the initiation and maintenance of early developmental gene repression.

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

BACKGROUND

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. We set out to bridge this gap and to understand the interactions between TEs and their cichlid hosts.

RESULTS

Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed strong conservation of TE silencing factors in cichlids, and 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 targets 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.

CONCLUSIONS

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.

SMARCA4 regulates the NK-mediated killing of senescent cells

Induction of senescence by chemotherapeutic agents arrests cancer cells and activates immune surveillance responses to contribute to therapy outcomes. In this investigation, we searched for ways to enhance the NK-mediated elimination of senescent cells. We used a staggered screen approach, first identifying siRNAs potentiating the secretion of immunomodulatory cytokines to later test for their ability to enhance NK-mediated killing of senescent cells. We identified that genetic or pharmacological inhibition of SMARCA4 enhanced senescent cell elimination by NK cells. SMARCA4 expression is elevated during senescence and its inhibition derepresses repetitive elements, inducing the SASP via activation of cGAS/STING and MAVS/MDA5 pathways. Moreover, a PROTAC targeting SMARCA4 synergized with cisplatin to increase the infiltration of CD8 T cells and mature, activated NK cells in an immunocompetent model of ovarian cancer. Our results indicate that SMARCA4 inhibitors enhance NK-mediated surveillance of senescent cells and may represent senotherapeutic interventions for ovarian cancer.

RNA binding by Periphilin plays an essential role in initiating silencing by the HUSH complex

The human silencing hub (HUSH) complex is a transcription-dependent, epigenetic repressor complex that provides a genome-wide immunosurveillance system for the recognition and silencing of newly-integrated retroelements. The core HUSH complex of TASOR, MPP8 and Periphilin, represses these retroelements through SETDB1-mediated H3K9me3 deposition and MORC2-dependent chromatin compaction. HUSH-dependent silencing is RNA-mediated, yet no HUSH component contains a recognised RNA-binding domain. Here we used an unbiased approach to identify which HUSH component was able to bind RNA and determine whether RNA-binding was essential for HUSH function. We identify Periphilin as the major RNA-binding component of the HUSH complex and show that Periphilin's N-terminal domain is essential for both RNA binding and HUSH function. Periphilin binding to RNA was independent of its interaction with TASOR or MPP8, as its N-terminal domain was sufficient for RNA targeting. The artificial tethering of Periphilin to a HUSH-insensitive, nascent transcript, enabled the HUSH-dependent silencing of the transcript. This tethering of Periphilin allowed the RNA-binding region of Periphilin to be removed such that only its C-terminal domain was required for oligomerisation and interaction with TASOR. We therefore show that Periphilin is the predominant RNA-binding protein of the HUSH complex and this RNA-binding is essential for HUSH activity.

Proteolethargy is a pathogenic mechanism in chronic disease

The pathogenic mechanisms of many diseases are well understood at the molecular level, but there are prevalent syndromes associated with pathogenic signaling, such as diabetes and chronic inflammation, where our understanding is more limited. Here, we report that pathogenic signaling suppresses the mobility of a spectrum of proteins that play essential roles in cellular functions known to be dysregulated in these chronic diseases. The reduced protein mobility, which we call proteolethargy, was linked to cysteine residues in the affected proteins and signaling-related increases in excess reactive oxygen species. Diverse pathogenic stimuli, including hyperglycemia, dyslipidemia, and inflammation, produce similar reduced protein mobility phenotypes. We propose that proteolethargy is an overlooked cellular mechanism that may account for various pathogenic features of diverse chronic diseases.

Multi-tissue characterization of the constitutive heterochromatin proteome in Drosophila identifies a link between satellite DNA organization and transposon repression

Noncoding satellite DNA repeats are abundant at the pericentromeric heterochromatin of eukaryotic chromosomes. During interphase, sequence-specific DNA-binding proteins cluster these repeats from multiple chromosomes into nuclear foci known as chromocenters. Despite the pivotal role of chromocenters in cellular processes like genome encapsulation and gene repression, the associated proteins remain incompletely characterized. Here, we use 2 satellite DNA-binding proteins, D1 and Prod, as baits to characterize the chromocenter-associated proteome in Drosophila embryos, ovaries, and testes through quantitative mass spectrometry. We identify D1- and Prod-associated proteins, including known heterochromatin proteins as well as proteins previously unlinked to satellite DNA or chromocenters, thereby laying the foundation for a comprehensive understanding of cellular functions enabled by satellite DNA repeats and their associated proteins. Interestingly, we find that multiple components of the transposon-silencing piRNA pathway are associated with D1 and Prod in embryos. Using genetics, transcriptomics, and small RNA profiling, we show that flies lacking D1 during embryogenesis exhibit transposon expression and gonadal atrophy as adults. We further demonstrate that this gonadal atrophy can be rescued by mutating the checkpoint kinase, Chk2, which mediates germ cell arrest in response to transposon mobilization. Thus, we reveal that a satellite DNA-binding protein functions during embryogenesis to silence transposons, in a manner that is heritable across later stages of development.

The chromatin remodeling factor OsINO80 promotes H3K27me3 and H3K9me2 deposition and maintains TE silencing in rice

The INO80 chromatin remodeling complex plays a critical role in shaping the dynamic chromatin environment. The diverse functions of the evolutionarily conserved INO80 complex have been widely reported. However, the role of INO80 in modulating the histone variant H2A.Z is controversial. Moreover, whether INO80 helps regulate heterochromatin remains unknown. Here, we characterize the regulatory effects of OsINO80 on protein-coding genes and transposable elements (TEs) in rice. Upon OsINO80 overexpression in rice, we found three types of OsINO80-occupied regions with different chromatin signatures: type I (enriched with H2A.Z), type II (enriched with H3K9me2), and type III (deficient in H2A.Z/H3K9me2). Loss of OsINO80 results in a decrease in H3K27me3, but not H2A.Z, at type I regions as well as a decrease in H3K9me2 at type II regions, which correlates with TE activation and transposition. Our findings reveal that OsINO80 facilitates H3K27me3 establishment, promotes H3K9me2 deposition, and maintains TE silencing.

Coordinated repression of totipotency-associated gene loci by histone methyltransferase EHMT2 through binding to LINE-1 regulatory elements

Mouse embryonic stem cells (mESCs) and other naïve pluripotent stem cells can reverse typical developmental trajectories and, at low frequency, de-differentiate into 2-cell-like cells (2CLCs) that resemble the mammalian embryo during zygotic genome activation (ZGA). This affords the opportunity to reveal molecular principles that govern the pre-implantation stages of mammalian development. We leveraged a multipurpose allele for acute protein depletion and efficient immunoprecipitation to dissect the molecular functions of the chromatin repressor EHMT2, a candidate antagonist of the mESC-to-2CLC transition. This allowed us to define categories of EHMT2 target genes characterized by distinct modes of EHMT2 chromatin engagement and repression. Most notably, EHMT2 directly represses large clusters of co-regulated gene loci that comprise a significant fraction of the 2CLC-specific transcriptome by initiating H3K9me2 spreading from distal LINE-1 elements. EHMT2 counteracts the recruitment of the activator DPPA2/4 to promoter-proximal endogenous retroviral elements (ERVs) at 2CLC genes. EHMT2 depletion elevates the expression of ZGA-associated transcripts in 2CLCs and synergizes with spliceosome inhibition and retinoic acid signaling in facilitating the mESC-to-2CLC transition. In contrast to ZGA-associated genes, repression of germ layer-associated transcripts by EHMT2 occurs outside of gene clusters in collaboration with ZFP462 and entails binding to non-repeat enhancers. Our observations show that EHMT2 attenuates the bidirectional differentiation potential of mouse pluripotent stem cells and define molecular modes for locus-specific transcriptional repression by this essential histone methyltransferase.

TPR is required for cytoplasmic chromatin fragment formation during senescence

During oncogene-induced senescence there are striking changes in the organisation of heterochromatin in the nucleus. This is accompanied by activation of a pro-inflammatory gene expression programme - the senescence-associated secretory phenotype (SASP) - driven by transcription factors such as NF-κB. The relationship between heterochromatin re-organisation and the SASP has been unclear. Here, we show that TPR, a protein of the nuclear pore complex basket required for heterochromatin re-organisation during senescence, is also required for the very early activation of NF-κB signalling during the stress-response phase of oncogene-induced senescence. This is prior to activation of the SASP and occurs without affecting NF-κB nuclear import. We show that TPR is required for the activation of innate immune signalling at these early stages of senescence and we link this to the formation of heterochromatin-enriched cytoplasmic chromatin fragments thought to bleb off from the nuclear periphery. We show that HMGA1 is also required for cytoplasmic chromatin fragment formation. Together these data suggest that re-organisation of heterochromatin is involved in altered structural integrity of the nuclear periphery during senescence, and that this can lead to activation of cytoplasmic nucleic acid sensing, NF-κB signalling, and activation of the SASP.

Human immunodeficiency virus-1 induces host genomic R-loops and preferentially integrates its genome near the R-loop regions

Although HIV-1 integration sites favor active transcription units in the human genome, high-resolution analysis of individual HIV-1 integration sites has shown that the virus can integrate into a variety of host genomic locations, including non-genic regions. The invisible infection by HIV-1 integrating into non-genic regions, challenging the traditional understanding of HIV-1 integration site selection, is more problematic because they are selected for preservation in the host genome during prolonged antiretroviral therapies. Here, we showed that HIV-1 integrates its viral genome into the vicinity of R-loops, a genomic structure composed of DNA-RNA hybrids. VSV-G-pseudotyped HIV-1 infection initiates the formation of R-loops in both genic and non-genic regions of the host genome and preferentially integrates into R-loop-rich regions. Using a HeLa cell model that can independently control transcriptional activity and R-loop formation, we demonstrated that the exogenous formation of R-loops directs HIV-1 integration-targeting sites. We also found that HIV-1 integrase proteins physically bind to the host genomic R-loops. These findings provide novel insights into the mechanisms underlying retroviral integration and the new strategies for antiretroviral therapy against HIV-1 latent infection.

Separable Cell Cycle Arrest and Immune Response Elicited through Pharmacological CDK4/6 and MEK Inhibition in RASmut Disease Models

The combination of CDK4/6 and MEK inhibition as a therapeutic strategy has shown promise in various cancer models, particularly in those harboring RAS mutations. An initial high-throughput drug screen identified high synergy between the CDK4/6 inhibitor palbociclib and the MEK inhibitor trametinib when used in combination in soft tissue sarcomas. In RAS mutant models, combination treatment with palbociclib and trametinib induced significant G1 cell cycle arrest, resulting in a marked reduction in cell proliferation and growth. CRISPR-mediated RB1 depletion resulted in a decreased response to CDK4/6 and MEK inhibition, which was validated in both cell culture and xenograft models. Beyond its cell cycle inhibitory effects, pathway enrichment analysis revealed the robust activation of interferon pathways upon CDK4/6 and MEK inhibition. This induction of gene expression was associated with the upregulation of retroviral elements. The TANK-binding kinase 1 inhibitor GSK8612 selectively blocked the induction of interferon-related genes induced by palbociclib and trametinib treatment and highlighted the separable epigenetic responses elicited by combined CDK4/6 and MEK inhibition. Together, these findings provide key mechanistic insights into the therapeutic potential of CDK4/6 and MEK inhibition in soft tissue sarcomas.

Activation characteristics of Ty3-retrotransposons after spaceflight and genetic stability of insertion sites in rice progeny

Introduction

The space environment is mutagenic and may induce genomic and phenotypic variations. Exploring the changes in transposon activity in the rice genome under space radiation is of great significance.

Methods

To analyze the activation characteristics of Ty3-retrotransposons and genetic stability of insertion sites in rice progeny after spaceflight, seeds of Nipponbare, DN416, and DN423 were exposed on board the SJ-10 recoverable satellite for 12.5 days. The differential methylation and transcription levels of Ty3-retrotransposons in the genome of Nipponbare's F0 generation after spaceflight, as well as the genetic stability of Ty3-retrotransposon insertion sites in DN416 and DN423 from F3 to F5 generations, was analyzed.

Results

The study found that the retrotransposons of ancient and young transposon families underwent demethylation from the tillering to heading stages of Nipponbare plants, which were F0 generation of space-exposed seeds, when the Nipponbare seeds were hit by single space high charge and energy (HZE) particles with LET ≥ 100 keV/μm. the transcription levels significantly increased in ancient transposon families (osr30, osr40, and rire10) and young transposon families (dagul, rn215-125, osr37, RLG_15, osr34, rire8, rire3, rire2, and hopi) (p ≤ 0.05) when LET > 100 keV/μm. Furthermore, the young Ty3-retrotransposons, which included the hopi, squiq, dasheng, rire2, rire3, rire8, osr34, rn_215-125, dagul, and RLG_15 families, underwent 1 to 8 transpositions in the F3 to F5 of DN416 and DN423 mutants, and some of these transposon insertion sites were stably inherited.

Discussion

The research holds great significance for understanding the activation characteristics of Ty3-retrotransposons in the rice genome induced by space radiation and the genetic characteristics of transposon insertion sites in its progeny.