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

Deficiency for SAMHD1 activates MDA5 in a cGAS/STING-dependent manner

Defects in nucleic acid metabolizing enzymes can lead to spontaneous but selective activation of either cGAS/STING or RIG-like receptor (RLR) signaling, causing type I interferon-driven inflammatory diseases. In these pathophysiological conditions, activation of the DNA sensor cGAS and IFN production are linked to spontaneous DNA damage. Physiological, or tonic, IFN signaling on the other hand is essential to functionally prime nucleic acid sensing pathways. Here, we show that low-level chronic DNA damage in mice lacking the Aicardi-Goutières syndrome gene SAMHD1 reduced tumor-free survival when crossed to a p53-deficient, but not to a DNA mismatch repair-deficient background. Increased DNA damage did not result in higher levels of type I interferon. Instead, we found that the chronic interferon response in SAMHD1-deficient mice was driven by the MDA5/MAVS pathway but required functional priming through the cGAS/STING pathway. Our work positions cGAS/STING upstream of tonic IFN signaling in Samhd1-deficient mice and highlights an important role of the pathway in physiological and pathophysiological innate immune priming.

The HUSH complex controls brain architecture and protocadherin fidelity

The HUSH (human silencing hub) complex contains the H3K9me3 binding protein M-phase phosphoprotein 8 (MPP8) and recruits the histone methyltransferase SETDB1 as well as Microrchidia CW-type zinc finger protein 2 (MORC2). Functional and mechanistic studies of the HUSH complex have hitherto been centered around SETDB1 while the in vivo functions of MPP8 and MORC2 remain elusive. Here, we show that genetic inactivation of Mphosph8 or Morc2a in the nervous system of mice leads to increased brain size, altered brain architecture, and behavioral changes. Mechanistically, in both mouse brains and human cerebral organoids, MPP8 and MORC2 suppress the repetitive-like protocadherin gene cluster in an H3K9me3-dependent manner. Our data identify MPP8 and MORC2, previously linked to silencing of repetitive elements via the HUSH complex, as key epigenetic regulators of protocadherin expression in the nervous system and thereby brain development and neuronal individuality in mice and humans.

The landscape of hervRNAs transcribed from human endogenous retroviruses across human body sites

BACKGROUND

Human endogenous retroviruses (HERVs), the remnants of ancient retroviruses, account for 8% of the human genome, but most have lost their transcriptional abilities under physiological conditions. However, mounting evidence shows that several expressed HERVs do exert biological functions. Here, we systematically characterize physiologically expressed HERVs and examine whether they may give insight into the molecular fundamentals of human development and disease.

RESULTS

We systematically identify 13,889 expressed HERVs across normal body sites and demonstrate that they are expressed in body site-specific patterns and also by sex, ethnicity, and age. Analyzing cis-ERV-related quantitative trait loci, we find that 5435 hervRNAs are regulated by genetic variants. Combining this with a genome-wide association study, we elucidate that the dysregulation of expressed HERVs might be associated with various complex diseases, particularly neurodegenerative and psychiatric diseases. We further find that physiologically activated hervRNAs are associated with histone modifications rather than DNA demethylation.

CONCLUSIONS

Our results present a locus-specific landscape of physiologically expressed hervRNAs, which represent a hidden layer of genetic architecture in development and disease.

Single-cell RNA sequencing highlights the functional role of human endogenous retroviruses in gallbladder cancer

Background

Gallbladder cancer (GBC), the most common malignancy of the biliary tract, shows late diagnosis and low survival rate and requires continued search for new diagnostic biomarkers and therapeutic targets. Human endogenous retroviruses (HERVs) are specifically prone to be reactivated in diverse cancers and are implicated in cancer progression and immunotherapy.

Methods

Single-cell RNA sequencing was performed on tumor tissues and paired adjacent tissues from 4 GBC patients. Dual-luciferase reporter assay was applied to measure enhancer activity of HERV sequences.

Findings

We dissected the cellular diversity and described the HERV transcriptomic landscape for GBC. We found that HERVs were transcribed in a cell type-specific manner and different HERV families were associated with diverse biological effects. HERVs could function as enhancers, presumably causing altered expression of neighboring genes. The transcription level of HERVH was gradually elevated with the malignant transformation of epithelial cells, suggesting HERVH may be a potential early diagnostic biomarker of GBC. HHLA2, a newly emerging immune checkpoint, was derived by HERVH, exhibited an expressional correlation with HERVH, and was identified as a promising target for immunotherapy.

Interpretation

Exploring the transcriptional landscape and potential functional impact of HERVs highlights the important role of HERVs in GBC and provides a fresh perspective on managing GBC.

Funding

This study was supported by the National Natural Science Foundation of China (31970176, 81972256) and the research grants from the Innovation Capacity Building Project of Jiangsu province (BM2020019).

A satellite DNA array barcodes chromosome 7 and regulates totipotency via ZFP819

Mammalian genomes are a battleground for genetic conflict between repetitive elements and KRAB-zinc finger proteins (KZFPs). We asked whether KZFPs can regulate cell fate by using ZFP819, which targets a satellite DNA array, ZP3AR. ZP3AR coats megabase regions of chromosome 7 encompassing genes encoding ZSCAN4, a master transcription factor of totipotency. Depleting ZFP819 in mouse embryonic stem cells (mESCs) causes them to transition to a 2-cell (2C)-like state, whereby the ZP3AR array switches from a poised to an active enhancer state. This is accompanied by a global erosion of heterochromatin roadblocks, which we link to decreased SETDB1 stability. These events result in transcription of active LINE-1 elements and impaired differentiation. In summary, ZFP819 and TRIM28 partner up to close chromatin across Zscan4, to promote exit from totipotency. We propose that satellite DNAs may control developmental fate transitions by barcoding and switching off master transcription factor genes.

Transcription start site signal profiling improves transposable element RNA expression analysis at locus-level

The transcriptional activity of Transposable Elements (TEs) has been involved in numerous pathological processes, including neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal lobar degeneration. The TE expression analysis from short-read sequencing technologies is, however, challenging due to the multitude of similar sequences derived from singular TEs subfamilies and the exaptation of TEs within longer coding or non-coding RNAs. Specialised tools have been developed to quantify the expression of TEs that either relies on probabilistic re-distribution of multimapper count fractions or allow for discarding multimappers altogether. Until now, the benchmarking across those tools was largely limited to aggregated expression estimates over whole TEs subfamilies. Here, we compared the performance of recently published tools (SQuIRE, TElocal, SalmonTE) with simplistic quantification strategies (featureCounts in unique, fraction and random modes) at the individual loci level. Using simulated datasets, we examined the false discovery rate and the primary driver of those false positive hits in the optimal quantification strategy. Our findings suggest a high false discovery number that exceeds the total number of correctly recovered active loci for all the quantification strategies, including the best performing tool TElocal. As a remedy, filtering based on the minimum number of read counts or baseMean expression improves the F1 score and decreases the number of false positives. Finally, we demonstrate that additional profiling of Transcription Start Site mapping statistics (using a k-means clustering approach) significantly improves the performance of TElocal while reporting a reliable set of detected and differentially expressed TEs in human simulated RNA-seq data.

SoloTE for improved analysis of transposable elements in single-cell RNA-Seq data using locus-specific expression

Transposable Elements (TEs) contribute to the repetitive fraction in almost every eukaryotic genome known to date, and their transcriptional activation can influence the expression of neighboring genes in healthy and disease states. Single cell RNA-Seq (scRNA-Seq) is a technical advance that allows the study of gene expression on a cell-by-cell basis. Although a current computational approach is available for the single cell analysis of TE expression, it omits their genomic location. Here we show SoloTE, a pipeline that outperforms the previous approach in terms of computational resources and by allowing the inclusion of locus-specific TE activity in scRNA-Seq expression matrixes. We then apply SoloTE to several datasets to reveal the repertoire of TEs that become transcriptionally active in different cell groups, and based on their genomic location, we predict their potential impact on gene expression. As our tool takes as input the resulting files from standard scRNA-Seq processing pipelines, we expect it to be widely adopted in single cell studies to help researchers discover patterns of cellular diversity associated with TE expression.

Exploration of nuclear body-enhanced sumoylation reveals that PML represses 2-cell features of embryonic stem cells

Membrane-less organelles are condensates formed by phase separation whose functions often remain enigmatic. Upon oxidative stress, PML scaffolds Nuclear Bodies (NBs) to regulate senescence or metabolic adaptation. PML NBs recruit many partner proteins, but the actual biochemical mechanism underlying their pleiotropic functions remains elusive. Similarly, PML role in embryonic stem cell (ESC) and retro-element biology is unsettled. Here we demonstrate that PML is essential for oxidative stress-driven partner SUMO2/3 conjugation in mouse ESCs (mESCs) or leukemia, a process often followed by their poly-ubiquitination and degradation. Functionally, PML is required for stress responses in mESCs. Differential proteomics unravel the KAP1 complex as a PML NB-dependent SUMO2-target in arsenic-treated APL mice or mESCs. PML-driven KAP1 sumoylation enables activation of this key epigenetic repressor implicated in retro-element silencing. Accordingly, Pml^-/- mESCs re-express transposable elements and display 2-Cell-Like features, the latter enforced by PML-controlled SUMO2-conjugation of DPPA2. Thus, PML orchestrates mESC state by coordinating SUMO2-conjugation of different transcriptional regulators, raising new hypotheses about PML roles in cancer.

Transcriptional dynamics of transposable elements in the type I IFN response in Myotis lucifugus cells

Background

Bats are a major reservoir of zoonotic viruses, and there has been growing interest in characterizing bat-specific features of innate immunity and inflammation. Recent studies have revealed bat-specific adaptations affecting interferon (IFN) signaling and IFN-stimulated genes (ISGs), but we still have a limited understanding of the genetic mechanisms that have shaped the evolution of bat immunity. Here we investigated the transcriptional and epigenetic dynamics of transposable elements (TEs) during the type I IFN response in little brown bat (Myotis lucifugus) primary embryonic fibroblast cells, using RNA-seq and CUT&RUN.

Results

We found multiple bat-specific TEs that undergo both locus-specific and family-level transcriptional induction in response to IFN. Our transcriptome reassembly identified multiple ISGs that have acquired novel exons from bat-specific TEs, including NLRC5, SLNF5 and a previously unannotated isoform of the IFITM2 gene. We also identified examples of TE-derived regulatory elements, but did not find strong evidence supporting genome-wide epigenetic activation of TEs in response to IFN.

Conclusion

Collectively, our study uncovers numerous TE-derived transcripts, proteins, and alternative isoforms that are induced by IFN in Myotis lucifugus cells, highlighting candidate loci that may contribute to bat-specific immune function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00277-z.

Mammalian genome innovation through transposon domestication

Since the discovery of transposons, their sheer abundance in host genomes has puzzled many. While historically viewed as largely harmless 'parasitic' DNAs during evolution, transposons are not a mere record of ancient genome invasion. Instead, nearly every element of transposon biology has been integrated into host biology. Here we review how host genome sequences introduced by transposon activities provide raw material for genome innovation and document the distinct evolutionary path of each species.

Detailed Analysis of Dorsal-Ventral Gradients of Gene Expression in the Hippocampus of Adult Rats

We performed RNA sequencing of the dorsal and ventral parts of the hippocampus and compared it with previously published data to determine the differences in the dorsoventral gradients of gene expression that may result from biological or technical variability. Our data suggest that the dorsal and ventral parts of the hippocampus differ in the expression of genes related to signaling pathways mediated by classical neurotransmitters (glutamate, GABA, monoamines, etc.) as well as peptide and Wnt ligands. These hippocampal parts also diverge in the expression of axon-guiding molecules (both receptors and ligands) and splice isoforms of genes associated with intercellular signaling and cell adhesion. Furthermore, analysis of differential expressions of genes specific for astrocytes, microglia, oligodendrocytes, and vascular cells suggests that non-neuronal cells may also differ in the characteristics between hippocampal parts. Analysis of expression of transposable elements showed that depletion of ribosomal RNA strongly increased the representation of transposable elements in the RNA libraries and helped to detect a weak predominance of expression of these elements in the ventral hippocampus. Our data revealed new molecular dimensions of functional differences between the dorsal and ventral hippocampus and points to possible cascades that may be involved in the longitudinal organization of the hippocampus.

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

Transcriptomic diversity in human medullary thymic epithelial cells

The induction of central T cell tolerance in the thymus depends on the presentation of peripheral self-epitopes by medullary thymic epithelial cells (mTECs). This promiscuous gene expression (pGE) drives mTEC transcriptomic diversity, with non-canonical transcript initiation, alternative splicing, and expression of endogenous retroelements (EREs) representing important but incompletely understood contributors. Here we map the expression of genome-wide transcripts in immature and mature human mTECs using high-throughput 5' cap and RNA sequencing. Both mTEC populations show high splicing entropy, potentially driven by the expression of peripheral splicing factors. During mTEC maturation, rates of global transcript mis-initiation increase and EREs enriched in long terminal repeat retrotransposons are up-regulated, the latter often found in proximity to differentially expressed genes. As a resource, we provide an interactive public interface for exploring mTEC transcriptomic diversity. Our findings therefore help construct a map of transcriptomic diversity in the healthy human thymus and may ultimately facilitate the identification of those epitopes which contribute to autoimmunity and immune recognition of tumor antigens.

TET1 regulates gene expression and repression of endogenous retroviruses independent of DNA demethylation

DNA methylation (5-methylcytosine (5mC)) is critical for genome stability and transcriptional regulation in mammals. The discovery that ten-eleven translocation (TET) proteins catalyze the oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) revolutionized our perspective on the complexity and regulation of DNA modifications. However, to what extent the regulatory functions of TET1 can be attributed to its catalytic activity remains unclear. Here, we use genome engineering and quantitative multi-omics approaches to dissect the precise catalytic vs. non-catalytic functions of TET1 in murine embryonic stem cells (mESCs). Our study identifies TET1 as an essential interaction hub for multiple chromatin modifying complexes and a global regulator of histone modifications. Strikingly, we find that the majority of transcriptional regulation depends on non-catalytic functions of TET1. In particular, we show that TET1 is critical for the establishment of H3K9me3 and H4K20me3 at endogenous retroviral elements (ERVs) and their silencing that is independent of its canonical role in DNA demethylation. Furthermore, we provide evidence that this repression of ERVs depends on the interaction between TET1 and SIN3A. In summary, we demonstrate that the non-catalytic functions of TET1 are critical for regulation of gene expression and the silencing of endogenous retroviruses in mESCs.

TEspeX: consensus-specific quantification of transposable element expression preventing biases from exonized fragments

SUMMARY

Transposable elements (TEs) play key roles in crucial biological pathways. Therefore, several tools enabling the quantification of their expression were recently developed. However, many of the existing tools lack the capability to distinguish between the transcription of autonomously expressed TEs and TE fragments embedded in canonical coding/non-coding non-TE transcripts. Consequently, an apparent change in the expression of a given TE may simply reflect the variation in the expression of the transcripts containing TE-derived sequences. To overcome this issue, we have developed TEspeX, a pipeline for the quantification of TE expression at the consensus level. TEspeX uses Illumina RNA-seq short reads to quantify TE expression avoiding counting reads deriving from inactive TE fragments embedded in canonical transcripts.

AVAILABILITY AND IMPLEMENTATION

The tool is implemented in python3, distributed under the GNU General Public License (GPL) and available on Github at https://github.com/fansalon/TEspeX (Zenodo URL: https://doi.org/10.5281/zenodo.6800331).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Somatic retrotransposition in the developing rhesus macaque brain

The retrotransposon LINE-1 (L1) is central to the recent evolutionary history of the human genome and continues to drive genetic diversity and germline pathogenesis. However, the spatiotemporal extent and biological significance of somatic L1 activity are poorly defined and are virtually unexplored in other primates. From a single L1 lineage active at the divergence of apes and Old World monkeys, successive L1 subfamilies have emerged in each descendant primate germline. As revealed by case studies, the presently active human L1 subfamily can also mobilize during embryonic and brain development in vivo. It is unknown whether nonhuman primate L1s can similarly generate somatic insertions in the brain. Here we applied approximately 40× single-cell whole-genome sequencing (scWGS), as well as retrotransposon capture sequencing (RC-seq), to 20 hippocampal neurons from two rhesus macaques (Macaca mulatta). In one animal, we detected and PCR-validated a somatic L1 insertion that generated target site duplications, carried a short 5' transduction, and was present in ∼7% of hippocampal neurons but absent from cerebellum and nonbrain tissues. The corresponding donor L1 allele was exceptionally mobile in vitro and was embedded in PRDM4, a gene expressed throughout development and in neural stem cells. Nanopore long-read methylome and RNA-seq transcriptome analyses indicated young retrotransposon subfamily activation in the early embryo, followed by repression in adult tissues. These data highlight endogenous macaque L1 retrotransposition potential, provide prototypical evidence of L1-mediated somatic mosaicism in a nonhuman primate, and allude to L1 mobility in the brain over the past 30 million years of human evolution.

Epigenomic analysis reveals prevalent contribution of transposable elements to cis-regulatory elements, tissue-specific expression, and alternative promoters in zebrafish

Transposable elements (TEs) encode regulatory elements that impact gene expression in multiple species, yet a comprehensive analysis of zebrafish TEs in the context of gene regulation is lacking. Here, we systematically investigate the epigenomic and transcriptomic landscape of TEs across 11 adult zebrafish tissues using multidimensional sequencing data. We find that TEs contribute substantially to a diverse array of regulatory elements in the zebrafish genome and that 37% of TEs are positioned in active regulatory states in adult zebrafish tissues. We identify TE subfamilies enriched in highly specific regulatory elements among different tissues. We use transcript assembly to discover TE-derived transcriptional units expressed across tissues. Finally, we show that novel TE-derived promoters can initiate tissue-specific transcription of alternate gene isoforms. This work provides a comprehensive profile of TE activity across normal zebrafish tissues, shedding light on mechanisms underlying the regulation of gene expression in this widely used model organism.

Retroviral Elements in Pathophysiology and as Therapeutic Targets for Amyotrophic Lateral Sclerosis

The study of the role of retroviruses in amyotrophic lateral sclerosis (ALS) dates back to the 1960s shortly after transposable elements themselves were first discovered. It was quickly realized that in wild mice both horizontal and vertical transmissions of retroviral elements were key to the development of an ALS-like syndrome leading to the postulate that endogenous retroviruses (ERVs) contribute significantly to the pathogenicity of this disease. Subsequent studies identified retroviral reverse transcriptase activity in brains of individuals with ALS from Guam. However, except for a single study from the former Soviet Union, ALS could not be transmitted to rhesus macaques. The discovery of an ALS-like syndrome in human immunodeficiency virus (HIV) and human T cell leukemia virus infected individuals led to renewed interest in the field and reverse transcriptase activity was found in the blood and cerebrospinal fluid of individuals with sporadic ALS. However, exogenous retroviruses could not be found in individuals with ALS which further reinforced the possibility of involvement of a human ERV (HERV). The first demonstration of the involvement of a HERV was the discovery of the activation of human endogenous retrovirus-K subtype HML-2 in the brains of individuals with ALS. The envelope protein of HML-2 is neurotoxic and transgenic animals expressing the envelope protein develop an ALS-like syndrome. Activation of HML-2 occurs in the context of generalized transposable element activation and is not specific for ALS. Individuals with HIV-associated ALS show a remarkable response to antiretroviral therapy; however, antiretroviral trials in ALS down-regulate HML-2 without ameliorating the disease. This highlights the need for specific drugs to be developed against HML-2 as a novel therapeutic target for ALS. Other approaches might include antisense oligonucleotides, shRNA targeted against the envelope gene or antibodies that can target the extracellular envelope protein. Future clinical trials in ALS should consider combination therapies to control these ERVs.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Recapitulating early human development with 8C-like cells

In human embryos, major zygotic genome activation (ZGA) initiates at the eight-cell (8C) stage. Abnormal ZGA leads to developmental defects and even contributes to the failure of human blastocyst formation or implantation. An in vitro cell model mimicking human 8C blastomeres would be invaluable to understanding the mechanisms regulating key biological events during early human development. Using the non-canonical promoter of LEUTX that putatively regulates human ZGA, we developed an 8C::mCherry reporter, which specifically marks the 8C state, to isolate rare 8C-like cells (8CLCs) from human preimplantation epiblast-like stem cells. The 8CLCs express a panel of human ZGA genes and have a unique transcriptome resembling that of the human 8C embryo. Using the 8C::mCherry reporter, we further optimize the chemical-based culture condition to increase and maintain the 8CLC population. Functionally, 8CLCs can self-organize to form blastocyst-like structures. The discovery and maintenance of 8CLCs provide an opportunity to recapitulate early human development.

piRNA-like small RNAs target transposable elements in a Clade IV parasitic nematode

The small RNA (sRNA) pathways identified in the model organism Caenorhabditis elegans are not widely conserved across nematodes. For example, the PIWI pathway and PIWI-interacting RNAs (piRNAs) are involved in regulating and silencing transposable elements (TE) in most animals but have been lost in nematodes outside of the C. elegans group (Clade V), and little is known about how nematodes regulate TEs in the absence of the PIWI pathway. Here, we investigated the role of sRNAs in the Clade IV parasitic nematode Strongyloides ratti by comparing two genetically identical adult stages (the parasitic female and free-living female). We identified putative small-interfering RNAs, microRNAs and tRNA-derived sRNA fragments that are differentially expressed between the two adult stages. Two classes of sRNAs were predicted to regulate TE activity including (i) a parasite-associated class of 21-22 nt long sRNAs with a 5' uridine (21-22Us) and a 5' monophosphate, and (ii) 27 nt long sRNAs with a 5' guanine/adenine (27GAs) and a 5' modification. The 21-22Us show striking resemblance to the 21U PIWI-interacting RNAs found in C. elegans, including an AT rich upstream sequence, overlapping loci and physical clustering in the genome. Overall, we have shown that an alternative class of sRNAs compensate for the loss of piRNAs and regulate TE activity in nematodes outside of Clade V.

Single-cell RNA-seq-based proteogenomics identifies glioblastoma-specific transposable elements encoding HLA-I-presented peptides

We analyze transposable elements (TEs) in glioblastoma (GBM) patients using a proteogenomic pipeline that combines single-cell transcriptomics, bulk RNA sequencing (RNA-seq) samples from tumors and healthy-tissue cohorts, and immunopeptidomic samples. We thus identify 370 human leukocyte antigen (HLA)-I-bound peptides encoded by TEs differentially expressed in GBM. Some of the peptides are encoded by repeat sequences from intact open reading frames (ORFs) present in up to several hundred TEs from recent long interspersed nuclear element (LINE)-1, long terminal repeat (LTR), and SVA subfamilies. Other HLA-I-bound peptides are encoded by single copies of TEs from old subfamilies that are expressed recurrently in GBM tumors and not expressed, or very infrequently and at low levels, in healthy tissues (including brain). These peptide-coding, GBM-specific, highly recurrent TEs represent potential tumor-specific targets for cancer immunotherapies.

The lysine deacetylase activity of histone deacetylases 1 and 2 is required to safeguard zygotic genome activation in mice and cattle

The genome is transcriptionally inert at fertilization and must be activated through a remarkable developmental process called zygotic genome activation (ZGA). Epigenetic reprogramming contributes significantly to the dynamic gene expression during ZGA; however, the mechanism has yet to be resolved. Here, we find histone deacetylases 1 and 2 (HDAC1/2) can regulate ZGA through lysine deacetylase activity. Notably, in mouse embryos, overexpression of a HDAC1/2 dominant-negative mutant leads to developmental arrest at the two-cell stage. RNA-seq reveals that 64% of downregulated genes are ZGA genes and 49% of upregulated genes are developmental genes. Inhibition of the deacetylase activity of HDAC1/2 causes a failure of histone deacetylation at multiple sites, including H4K5, H4K16, H3K14, H3K18 and H3K27. ChIP-seq analysis exhibits an increase and decrease of H3K27ac enrichment at promoters of up- and downregulated genes, respectively. Moreover, HDAC1 mutants prohibit the removal of H3K4me3 by impeding expression of Kdm5 genes. Importantly, the developmental block can be greatly rescued by Kdm5b injection and by partially correcting the expression of the majority of dysregulated genes. Similar functional significance of HDAC1/2 is conserved in bovine embryos. Overall, we propose that HDAC1/2 are indispensable for ZGA by creating correct transcriptional repressive and active states in mouse and bovine embryos.

Snrpb is required in murine neural crest cells for proper splicing and craniofacial morphogenesis

Heterozygous mutations in SNRPB, an essential core component of the five small ribonucleoprotein particles of the spliceosome, are responsible for cerebrocostomandibular syndrome (CCMS). We show that Snrpb heterozygous mouse embryos arrest shortly after implantation. Additionally, heterozygous deletion of Snrpb in the developing brain and neural crest cells models craniofacial malformations found in CCMS, and results in death shortly after birth. RNAseq analysis of mutant heads prior to morphological defects revealed increased exon skipping and intron retention in association with increased 5′ splice site strength. We found increased exon skipping in negative regulators of the P53 pathway, along with increased levels of nuclear P53 and P53 target genes. However, removing Trp53 in Snrpb heterozygous mutant neural crest cells did not completely rescue craniofacial development. We also found a small but significant increase in exon skipping of several transcripts required for head and midface development, including Smad2 and Rere. Furthermore, mutant embryos exhibited ectopic or missing expression of Fgf8 and Shh, which are required to coordinate face and brain development. Thus, we propose that mis-splicing of transcripts that regulate P53 activity and craniofacial-specific genes contributes to craniofacial malformations.

This article has an associated First Person interview with the first author of the paper.

Summary: We report the first mouse model for cerebrocostomandibular syndrome, showing that mis-splicing of transcripts that regulate P53 activity and craniofacial-specific genes contributes to craniofacial malformations.

Functional editing of endogenous genes through rapid selection of cell pools (Rapid generation of endogenously tagged genes in Drosophila ovarian somatic sheath cells)

The combination of genome-editing and epitope tagging provides a powerful strategy to study proteins with high affinity and specificity while preserving their physiological expression patterns. However, stably modifying endogenous genes in cells that do not allow for clonal selection has been challenging. Here, we present a simple and fast strategy to generate stable, endogenously tagged alleles in a non-transformed cell culture model. At the example of piwi in Drosophila ovarian somatic sheath cells, we show that this strategy enables the generation of an N-terminally tagged protein that emulates the expression level and subcellular localization of the wild type protein and forms functional Piwi-piRNA complexes. We further present a concise workflow to establish endogenously N-terminally and C-terminally tagged proteins, and knockout alleles through rapid selection of cell pools in fly and human models.

ExplorATE: A new pipeline to explore active transposable elements from RNA-seq data

MOTIVATION

Transposable elements (TEs) are ubiquitous in genomes and many remain active. TEs comprise an important fraction of the transcriptomes with potential effects on the host genome, either by generating deleterious mutations or promoting evolutionary novelties. However, their functional study is limited by the difficulty in their identification and quantification, particularly in non-model organisms.

RESULTS

We developed a new pipeline (ExplorATE or Explore Active Transposable Elements) implemented in R and bash that allows the quantification of active TEs in both model and non-model organisms. ExplorATE creates TE-specific indexes and uses the Selective Alignment (SA) to filter out co-transcribed transposons within genes based on alignment scores. Moreover, our software incorporates a Wicker-like criteria to refine a set of target TEs and avoid spurious mapping. Based on simulated and real data, we show that the SA strategy adopted by ExplorATE achieved better estimates of non-co-transcribed elements than other available alignment-based or mapping-based software. ExplorATE results showed high congruence with alignment-based tools with and without a reference genome, yet ExplorATE required less execution time. Likewise, ExplorATE expands and complements most previous TE analyses by incorporating the co-transcription and multi-mapping effects during quantification, and provides a seamless integration with other downstream tools within the R environment.

AVAILABILITY

Source code is available at https://github.com/FemeniasM/ExplorATEproject and https://github.com/FemeniasM/ExplorATE_shell_script.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comprehensive analysis of both long and short read transcriptomes of a clonal and a seed-propagated model species reveal the prerequisites for transcriptional activation of autonomous and non-autonomous transposons in plants

Background

Transposable element (TE) transcription is a precursor to its mobilisation in host genomes. However, the characteristics of expressed TE loci, the identification of self-competent transposon loci contributing to new insertions, and the genomic conditions permitting their mobilisation remain largely unknown.

Results

Using Vitis vinifera embryogenic callus, we explored the impact of biotic stressors on transposon transcription through the exposure of the callus to live cultures of an endemic grapevine yeast, Hanseniaspora uvarum. We found that only 1.7–2.5% of total annotated TE loci were transcribed, of which 5–10% of these were full-length, and the expressed TE loci exhibited a strong location bias towards expressed genes. These trends in transposon transcription were also observed in RNA-seq data from Arabidopsis thaliana wild-type plants but not in epigenetically compromised Arabidopsis ddm1 mutants. Moreover, differentially expressed TE loci in the grapevine tended to share expression patterns with co-localised differentially expressed genes. Utilising nanopore cDNA sequencing, we found a strong correlation between the inclusion of intronic TEs in gene transcripts and the presence of premature termination codons in these transcripts. Finally, we identified low levels of full-length transcripts deriving from structurally intact TE loci in the grapevine model.

Conclusion

Our observations in two disparate plant models representing clonally and seed propagated plant species reveal a closely connected transcriptional relationship between TEs and co-localised genes, particularly when epigenetic silencing is not compromised. We found that the stress treatment alone was insufficient to induce large-scale full-length transcription from structurally intact TE loci, a necessity for non-autonomous and autonomous mobilisation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00271-5.

Transposable element accumulation drives size differences among polymorphic Y Chromosomes in Drosophila

Y Chromosomes of many species are gene poor and show low levels of nucleotide variation, yet often display high amounts of structural diversity. Dobzhansky cataloged several morphologically distinct Y Chromosomes in Drosophila pseudoobscura that differ in size and shape, but the molecular causes of their dramatic size differences are unclear. Here we use cytogenetics and long-read sequencing to study the sequence content of polymorphic Y Chromosomes in D. pseudoobscura We show that Y Chromosomes differ almost 2-fold in size, ranging from 30 to 60 Mb. Most of this size difference is caused by a handful of active transposable elements (TEs) that have recently expanded on the largest Y Chromosome, with different elements being responsible for Y expansion on differently sized D. pseudoobscura Y's. We show that Y Chromosomes differ in their heterochromatin enrichment, expression of Y-enriched TEs, and also influence expression of dozens of autosomal and X-linked genes. The same helitron element that showed the most drastic amplification on the largest Y in D. pseudoobscura independently amplified on a polymorphic large Y Chromosome in D. affinis, suggesting that some TEs are inherently more prone to become deregulated on Y Chromosomes.

Locus-specific chromatin profiling of evolutionarily young transposable elements

Despite a vast expansion in the availability of epigenomic data, our knowledge of the chromatin landscape at interspersed repeats remains highly limited by difficulties in mapping short-read sequencing data to these regions. In particular, little is known about the locus-specific regulation of evolutionarily young transposable elements (TEs), which have been implicated in genome stability, gene regulation and innate immunity in a variety of developmental and disease contexts. Here we propose an approach for generating locus-specific protein-DNA binding profiles at interspersed repeats, which leverages information on the spatial proximity between repetitive and non-repetitive genomic regions. We demonstrate that the combination of HiChIP and a newly developed mapping tool (PAtChER) yields accurate protein enrichment profiles at individual repetitive loci. Using this approach, we reveal previously unappreciated variation in the epigenetic profiles of young TE loci in mouse and human cells. Insights gained using our method will be invaluable for dissecting the molecular determinants of TE regulation and their impact on the genome.

scTEM-seq: Single-cell analysis of transposable element methylation to link global epigenetic heterogeneity with transcriptional programs

Global changes in DNA methylation are observed in development and disease, and single-cell analyses are highlighting the heterogeneous regulation of these processes. However, technical challenges associated with single-cell analysis of DNA methylation limit these studies. We present single-cell transposable element methylation sequencing (scTEM-seq) for cost-effective estimation of average DNA methylation levels. By targeting high-copy SINE Alu elements, we achieve amplicon bisulphite sequencing with thousands of loci covered in each scTEM-seq library. Parallel transcriptome analysis is also performed to link global DNA methylation estimates with gene expression. We apply scTEM-seq to KG1a acute myeloid leukaemia (AML) cells, and primary AML cells. Our method reveals global DNA methylation heterogeneity induced by decitabine treatment of KG1a cells associated with altered expression of immune process genes. We also compare global DNA methylation estimates to expression of transposable elements and find a predominance of negative correlations. Finally, we observe co-ordinated upregulation of many transposable elements in a sub-set of decitabine treated cells. By linking global DNA methylation heterogeneity with transcription, scTEM-seq will refine our understanding of epigenetic regulation in cancer and beyond.

The PIWI/piRNA response is relaxed in a rodent that lacks mobilizing transposable elements

Transposable elements (TEs) are genomic parasites that can propagate throughout host genomes. Mammalian genomes are typically dominated by LINE retrotransposons and their associated SINEs, and germline mobilization is a challenge to genome integrity. There are defenses against TE proliferation and the PIWI/piRNA defense is among the most well understood. However, the PIWI/piRNA system has been investigated largely in animals with actively mobilizing TEs and it is unclear how the PIWI/piRNA system functions in the absence of mobilizing TEs. The 13-lined ground squirrel provides the opportunity to examine PIWI/piRNA and TE dynamics within the context of minimal, and possibly nonexistent, TE accumulation. To do so, we compared the PIWI/piRNA dynamics in squirrels to observations from the rabbit and mouse. Despite a lack of young insertions in squirrels, TEs were still actively transcribed at higher levels compared to mouse and rabbit. All three Piwi genes were not expressed, prior to P8 in squirrel testis, and there was little TE expression change with the onset of Piwi expression. We also demonstrated there was not a major expression change in the young squirrel LINE families in the transition from juvenile to adult testis in contrast to young mouse and rabbit LINE families. These observations lead us to conclude that PIWI suppression, was weaker for squirrel LINEs and SINEs and did not strongly reduce their transcription. We speculate that, although the PIWI/piRNA system is adaptable to novel TE threats, transcripts from TEs that are no longer threatening receive less attention from PIWI proteins.

Chromatin modifier HUSH co-operates with RNA decay factor NEXT to restrict transposable element expression

Transposable elements (TEs) are widespread genetic parasites known to be kept under tight transcriptional control. Here, we describe a functional connection between the mouse-orthologous “nuclear exosome targeting” (NEXT) and “human silencing hub” (HUSH) complexes, involved in nuclear RNA decay and the epigenetic silencing of TEs, respectively. Knocking out the NEXT component ZCCHC8 in embryonic stem cells results in elevated TE RNA levels. We identify a physical interaction between ZCCHC8 and the MPP8 protein of HUSH and establish that HUSH recruits NEXT to chromatin at MPP8-bound TE loci. However, while NEXT and HUSH both dampen TE RNA expression, their activities predominantly affect shorter non-polyadenylated and full-length polyadenylated transcripts, respectively. Indeed, our data suggest that the repressive action of HUSH promotes a condition favoring NEXT RNA decay activity. In this way, transcriptional and post-transcriptional machineries synergize to suppress the genotoxic potential of TE RNAs.

Garland et al. report a physical and functional connection between the NEXT complex, involved in RNA decay, and the HUSH complex, involved in chromatin regulation. Together, NEXT and HUSH cooperate to control transposable element (TE) RNA expression in embryonic stem cells, suppressing pA⁻ and pA⁺ transcripts, respectively.

H3K9 dimethylation safeguards cancer cells against activation of the interferon pathway

Activation of interferon genes constitutes an important anticancer pathway able to restrict proliferation of cancer cells. Here, we demonstrate that the H3K9me3 histone methyltransferase (HMT) suppressor of variegation 3-9 homolog 1 (SUV39H1) is required for the proliferation of acute myeloid leukemia (AML) and find that its loss leads to activation of the interferon pathway. Mechanistically, we show that this occurs via destabilization of a complex composed of SUV39H1 and the two H3K9me2 HMTs, G9A and GLP. Indeed, loss of H3K9me2 correlated with the activation of key interferon pathway genes, and interference with the activities of G9A/GLP largely phenocopied loss of SUV39H1. Last, we demonstrate that inhibition of G9A/GLP synergized with DNA demethylating agents and that SUV39H1 constitutes a potential biomarker for the response to hypomethylation treatment. Collectively, we uncovered a clinically relevant role for H3K9me2 in safeguarding cancer cells against activation of the interferon pathway.

LncRNA Biomarkers of Inflammation and Cancer

Long noncoding RNAs (lncRNAs) are promising candidates as biomarkers of inflammation and cancer. LncRNAs have several properties that make them well-suited as molecular markers of disease: (1) many lncRNAs are expressed in a tissue-specific manner, (2) distinct lncRNAs are upregulated based on different inflammatory or oncogenic stimuli, (3) lncRNAs released from cells are packaged and protected in extracellular vesicles, and (4) circulating lncRNAs in the blood are detectable using various RNA sequencing approaches. Here we focus on the potential for lncRNA biomarkers to detect inflammation and cancer, highlighting key biological, technological, and analytical considerations that will help advance the development of lncRNA-based liquid biopsies.

ATF7ip Targets Transposable Elements for H3K9me3 Deposition to Modify CD8+ T Cell Effector and Memory Responses

CD8+ T cells are critical for the immune response to pathogens and tumors, and CD8+ T cell memory protects against repeat infections. In this study, we identify the activating transcription factor 7 interacting protein (ATF7ip) as a critical regulator of CD8+ T cell immune responses. Mice with a T cell-specific deletion of ATF7ip have a CD8+ T cell-intrinsic enhancement of Il7r expression and Il2 expression leading to enhanced effector and memory responses. Chromatin immunoprecipitation sequencing studies identified ATF7ip as a repressor of Il7r and Il2 gene expression through the deposition of the repressive histone mark H3K9me3 at the Il7r gene and Il2-Il21 intergenic region. Interestingly, ATF7ip targeted transposable elements for H3K9me3 deposition at both the IL7r locus and the Il2-Il21 intergenic region, indicating that ATF7ip silencing of transposable elements is important for regulating CD8+ T cell function. These results demonstrate a new epigenetic pathway by which IL-7R and IL-2 production are constrained in CD8+ T cells, and this may open up new avenues for modulating their production.

Dicer promotes genome stability via the bromodomain transcriptional co-activator BRD4

RNA interference is required for post-transcriptional silencing, but also has additional roles in transcriptional silencing of centromeres and genome stability. However, these roles have been controversial in mammals. Strikingly, we found that Dicer-deficient embryonic stem cells have strong proliferation and chromosome segregation defects as well as increased transcription of centromeric satellite repeats, which triggers the interferon response. We conducted a CRISPR-Cas9 genetic screen to restore viability and identified transcriptional activators, histone H3K9 methyltransferases, and chromosome segregation factors as suppressors, resembling Dicer suppressors identified in independent screens in fission yeast. The strongest suppressors were mutations in the transcriptional co-activator Brd4, which reversed the strand-specific transcription of major satellite repeats suppressing the interferon response, and in the histone acetyltransferase Elp3. We show that identical mutations in the second bromodomain of Brd4 rescue Dicer-dependent silencing and chromosome segregation defects in both mammalian cells and fission yeast. This remarkable conservation demonstrates that RNA interference has an ancient role in transcriptional silencing and in particular of satellite repeats, which is essential for cell cycle progression and proper chromosome segregation. Our results have pharmacological implications for cancer and autoimmune diseases characterized by unregulated transcription of satellite repeats.

While RNA interference is conserved across species, small RNA pathways are very diverse. In this study, Gutbrod et al. find that non-canonical roles of Dicer in genome stability are in fact deeply conserved from yeast to humans.

Pathogenic variants in the human m6A reader YTHDC2 are associated with primary ovarian insufficiency

Primary ovarian insufficiency (POI) affects 1% of women and carries significant medical and psychosocial sequelae. Approximately 10% of POI has a defined genetic cause, with most implicated genes relating to biological processes involved in early fetal ovary development and function. Recently, Ythdc2, an RNA helicase and N6-methyladenosine (m6a) reader, has emerged as a novel regulator of meiosis in mice. Here, we describe homozygous pathogenic variants in YTHDC2 in three women with early-onset POI from two families: c. 2567C>G, p.P856R in the helicase-associated (HA2) domain; and c.1129G>T, p.E377*. We demonstrate that YTHDC2 is expressed in the developing human fetal ovary and is upregulated in meiotic germ cells, together with related meiosis-associated factors. The p.P856R variant results in a less flexible protein that likely disrupts downstream conformational kinetics of the HA2 domain, whereas the p.E377* variant truncates the helicase core. Taken together, our results reveal that YTHDC2 is a key new regulator of meiosis in humans and pathogenic variants within this gene are associated with POI.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Human reproduction is regulated by retrotransposons derived from ancient Hominidae-specific viral infections

Germ cells are essential to pass DNA from one generation to the next. In human reproduction, germ cell development begins with the specification of primordial germ cells (PGCs) and a failure to specify PGCs leads to human infertility. Recent studies have revealed that the transcription factor network required for PGC specification has diverged in mammals, and this has a significant impact on our understanding of human reproduction. Here, we reveal that the Hominidae-specific Transposable Elements (TEs) LTR5Hs, may serve as TEENhancers (TE Embedded eNhancers) to facilitate PGC specification. LTR5Hs TEENhancers become transcriptionally active during PGC specification both in vivo and in vitro with epigenetic reprogramming leading to increased chromatin accessibility, localized DNA demethylation, enrichment of H3K27ac, and occupation of key hPGC transcription factors. Inactivation of LTR5Hs TEENhancers with KRAB mediated CRISPRi has a significant impact on germ cell specification. In summary, our data reveals the essential role of Hominidae-specific LTR5Hs TEENhancers in human germ cell development.

The transcription factor network required for primordial germ cell (PGC) specification is known to diverge in mammals. Here the authors show that hominidae-specific transposable element (TE) LTR5Hs becomes transcriptionally active during PGC specification, and LTR5Hs inactivation abrogates human PGC specification

Genome-wide perturbations of Alu expression and Alu-associated post-transcriptional regulations distinguish oligodendroglioma from other gliomas

Alu is a primate-specific repeat element in the human genome and has been increasingly appreciated as a regulatory element in many biological processes. But the appreciation of Alu has been limited in tumorigenesis, especially for brain tumor. To investigate the relevance of Alu to the gliomagenesis, we studied Alu element-associated post-transcriptional processes and the RNA expression of the element by performing RNA-seq for a total of 41 pairs of neurotypical and diverse glioma brain tissues. We find that A-to-I editing and circular RNA levels, as well as Alu RNA expression, are decreased overall in gliomas, compared to normal tissue. Interestingly, grade 2 oligodendrogliomas are least affected in A-to-I editing and circular RNA levels among gliomas, whereas they have a higher proportion of down-regulated Alu subfamilies, compared to the other gliomas. These findings collectively imply a unique pattern of Alu-associated transcriptomes in grade 2 oligodendroglioma, providing an insight to gliomagenesis from the perspective of an evolutionary genetic element.

The insertion and dysregulation of transposable elements in osteosarcoma and their association with patient event-free survival

The dysregulation of transposable elements (TEs) has been explored in a variety of cancers. However, TE activities in osteosarcoma (OS) have not been extensively studied yet. By integrative analysis of RNA-seq, whole-genome sequencing (WGS), and methylation data, we showed aberrant TE activities associated with dysregulations of TEs in OS tumors. Specifically, expression levels of LINE-1 and Alu of different evolutionary ages, as well as subfamilies of SVA and HERV-K, were significantly up-regulated in OS tumors, accompanied by enhanced DNA repair responses. We verified the characteristics of LINE-1 mediated TE insertions, including target site duplication (TSD) length (centered around 15 bp) and preferential insertions into intergenic and AT-rich regions as well as intronic regions of longer genes. By filtering polymorphic TE insertions reported in 1000 genome project (1KGP), besides 148 tumor-specific somatic TE insertions, we found most OS patient-specific TE insertions (3175 out of 3326) are germline insertions, which are associated with genes involved in neuronal processes or with transcription factors important for cancer development. In addition to 68 TE-affected cancer genes, we found recurrent germline TE insertions in 72 non-cancer genes with high frequencies among patients. We also found that +/− 500 bps flanking regions of transcription start sites (TSS) of LINE-1 (young) and Alu showed lower methylation levels in OS tumor samples than controls. Interestingly, by incorporating patient clinical data and focusing on TE activities in OS tumors, our data analysis suggested that higher TE insertions in OS tumors are associated with a longer event-free survival time.

TASOR epigenetic repressor cooperates with a CNOT1 RNA degradation pathway to repress HIV

The Human Silencing Hub (HUSH) complex constituted of TASOR, MPP8 and Periphilin recruits the histone methyl-transferase SETDB1 to spread H3K9me3 repressive marks across genes and transgenes in an integration site-dependent manner. The deposition of these repressive marks leads to heterochromatin formation and inhibits gene expression, but the underlying mechanism is not fully understood. Here, we show that TASOR silencing or HIV-2 Vpx expression, which induces TASOR degradation, increases the accumulation of transcripts derived from the HIV-1 LTR promoter at a post-transcriptional level. Furthermore, using a yeast 2-hybrid screen, we identify new TASOR partners involved in RNA metabolism including the RNA deadenylase CCR4-NOT complex scaffold CNOT1. TASOR and CNOT1 synergistically repress HIV expression from its LTR. Similar to the RNA-induced transcriptional silencing complex found in fission yeast, we show that TASOR interacts with the RNA exosome and RNA Polymerase II, predominantly under its elongating state. Finally, we show that TASOR facilitates the association of RNA degradation proteins with RNA polymerase II and is detected at transcriptional centers. Altogether, we propose that HUSH operates at the transcriptional and post-transcriptional levels to repress HIV proviral expression.

An introduction to PIWI-interacting RNAs (piRNAs) in the context of metazoan small RNA silencing pathways

PIWI proteins and their associated PIWI-interacting RNAs (piRNAs) constitute a small RNA-based adaptive immune system that restricts the deleterious activity of mobile genetic elements to protect genome integrity. Self/nonself discrimination is at the very core of successful defence and relies on complementary base-pairing in RNA-guided immunity. How the millions of piRNA sequences faithfully discriminate between self and nonself and how they adapt to novel genomic invaders remain key outstanding questions in genome biology. This review aims to introduce principles of piRNA silencing in the context of metazoan small RNA pathways. A distinct feature of piRNAs is their origin from single-stranded instead of double-stranded RNA precursors, and piRNAs require a unique set of processing factors. Novel nucleases, helicases and RNA binding proteins have been identified in piRNA biology, and while we are starting to understand some mechanisms of piRNA biogenesis and function, this diverse and prolific class of small RNAs remains full of surprises.

Pathogenic tau accelerates aging-associated activation of transposable elements in the mouse central nervous system

Transposable elements comprise almost half of the mammalian genome. A growing body of evidence suggests that transposable element dysregulation accompanies brain aging and neurodegenerative disorders, and that transposable element activation is neurotoxic. Recent studies have identified links between pathogenic forms of tau, a protein that accumulates in Alzheimer's disease and related "tauopathies," and transposable element-induced neurotoxicity. Starting with transcriptomic analyses, we find that age- and tau-induced transposable element activation occurs in the mouse brain. Among transposable elements that are activated at the RNA level in the context of brain aging and tauopathy, we find that the endogenous retrovirus (ERV) class of retrotransposons is particularly enriched. We show that protein encoded by Intracisternal A-particle, a highly active mouse ERV, is elevated in brains of tau transgenic mice. Using two complementary approaches, we find that brains of tau transgenic mice contain increased DNA copy number of transposable elements, raising the possibility that these elements actively retrotranspose in the context of tauopathy. Taken together, our study lays the groundwork for future mechanistic studies focused on transposable element regulation in the aging mouse brain and in mouse models of tauopathy and provides support for ongoing therapeutic efforts targeting transposable element activation in patients with Alzheimer's disease.

Spatially Resolved Expression of Transposable Elements in Disease and Somatic Tissue with SpatialTE

Spatial transcriptomics (ST) is transforming the way we can study gene expression and its regulation through position-specific resolution within tissues. However, as in bulk RNA-Seq, transposable elements (TEs) are not being studied due to their highly repetitive nature. In recent years, TEs have been recognized as important regulators of gene expression, and thus, TE expression analysis in a spatially resolved manner could further help to understand their role in gene regulation within tissues. We present SpatialTE, a tool to analyze TE expression from ST datasets and show its application in somatic and diseased tissues. The results indicate that TEs have spatially regulated expression patterns and that their expression profiles are spatially altered in ALS disease, indicating that TEs might perform differential regulatory functions within tissue organs. We have made SpatialTE publicly available as open-source software under an MIT license.

DPPA2 and DPPA4 are dispensable for mouse zygotic genome activation and pre-implantation development

How maternal factors in oocytes initiate zygotic genome activation (ZGA) remains elusive in mammals, partly due to the challenge of de novo identification of key factors using scarce materials. Two-cell (2C)-like cells have been widely used as an in vitro model in order to understand mouse ZGA and totipotency because of their expression of a group of two-cell embryo-specific genes and their simplicity for genetic manipulation. Recent studies indicate that DPPA2 and DPPA4 are required for establishing the 2C-like state in mouse embryonic stem cells in a DUX-dependent manner. These results suggest that DPPA2 and DPPA4 are essential maternal factors that regulate Dux and ZGA in embryos. By analyzing maternal knockout and maternal-zygotic knockout embryos, we unexpectedly found that DPPA2 and DPPA4 are dispensable for Dux activation, ZGA and pre-implantation development. Our study suggests that 2C-like cells do not fully recapitulate two-cell embryos in terms of regulation of two-cell embryo-specific genes, and, therefore, caution should be taken when studying ZGA and totipotency using 2C-like cells as the model system.

Foxd3 controls heterochromatin-mediated repression of repeat elements and 2-cell state transcription

Repeat element transcription plays a vital role in early embryonic development. The expression of repeats such as MERVL characterises mouse embryos at the 2‐cell stage and defines a 2‐cell‐like cell (2CLC) population in a mouse embryonic stem cell culture. Repeat element sequences contain binding sites for numerous transcription factors. We identify the forkhead domain transcription factor FOXD3 as a regulator of major satellite repeats and MERVL transcription in mouse embryonic stem cells. FOXD3 binds to and recruits the histone methyltransferase SUV39H1 to MERVL and major satellite repeats, consequentially repressing the transcription of these repeats by the establishment of the H3K9me3 heterochromatin modification. Notably, depletion of FOXD3 leads to the de‐repression of MERVL and major satellite repeats as well as a subset of genes expressed in the 2‐cell state, shifting the balance between the stem cell and 2‐cell‐like population in culture. Thus, FOXD3 acts as a negative regulator of repeat transcription, ascribing a novel function to this transcription factor.

Transcriptional activity and epigenetic regulation of transposable elements in the symbiotic fungus Rhizophagus irregularis

Arbuscular mycorrhizal (AM) fungi form mutualistic relationships with most land plant species. AM fungi have long been considered as ancient asexuals. Long-term clonal evolution would be remarkable for a eukaryotic lineage and suggests the importance of alternative mechanisms to promote genetic variability facilitating adaptation. Here, we assessed the potential of transposable elements for generating such genomic diversity. The dynamic expression of TEs during Rhizophagus irregularis spore development suggests ongoing TE activity. We find Mutator-like elements located near genes belonging to highly expanded gene families. Whole-genome epigenomic profiling of R. irregularis provides direct evidence of DNA methylation and small RNA production occurring at TE loci. Our results support a model in which TE activity shapes the genome, while DNA methylation and small RNA-mediated silencing keep their overproliferation in check. We propose that a well-controlled TE activity directly contributes to genome evolution in AM fungi.