Roles of transposable elements in the regulation of mammalian transcription

From the cauldron of conflict: Endogenous gene regulation by piRNA and other modes of adaptation enabled by selfish transposable elements

Transposable elements (TEs) provide a prime example of genetic conflict because they can proliferate in genomes and populations even if they harm the host. However, numerous studies have shown that TEs, though typically harmful, can also provide fuel for adaptation. This is because they code functional sequences that can be useful for the host in which they reside. In this review, I summarize the "how" and "why" of adaptation enabled by the genetic conflict between TEs and hosts. In addition, focusing on mechanisms of TE control by small piwi-interacting RNAs (piRNAs), I highlight an indirect form of adaptation enabled by conflict. In this case, mechanisms of host defense that regulate TEs have been redeployed for endogenous gene regulation. I propose that the genetic conflict released by meiosis in early eukaryotes may have been important because, among other reasons, it spurred evolutionary innovation on multiple interwoven trajectories - on the part of hosts and also embedded genetic parasites. This form of evolution may function as a complexity generating engine that was a critical player in eukaryotic evolution.

Human lncRNAs harbor conserved modules embedded in different sequence contexts

We analyzed the structure of human long non-coding RNA (lncRNAs) genes to investigate whether the non-coding transcriptome is organized in modular domains, as is the case for protein-coding genes. To this aim, we compared all known human lncRNA exons and identified 340 pairs of exons with high sequence and/or secondary structure similarity but embedded in a dissimilar sequence context. We grouped these pairs in 106 clusters based on their reciprocal similarities. These shared modules are highly conserved between humans and the four great ape species, display evidence of purifying selection and likely arose as a result of recent segmental duplications. Our analysis contributes to the understanding of the mechanisms driving the evolution of the non-coding genome and suggests additional strategies towards deciphering the functional complexity of this class of molecules.

Rise and SINE: roles of transcription factors and retrotransposons in zygotic genome activation

In sexually reproducing organisms, life begins with the fusion of transcriptionally silent gametes, the oocyte and sperm. Although initiation of transcription in the embryo, known as zygotic genome activation (ZGA), is universally required for development, the transcription factors regulating this process are poorly conserved. In this Perspective, we discuss recent insights into the mechanisms of ZGA in totipotent mammalian embryos, namely ZGA regulation by several transcription factors, including by orphan nuclear receptors (OrphNRs) such as the pioneer transcription factor NR5A2, and by factors of the DUX, TPRX and OBOX families. We performed a meta-analysis and compiled a list of pan-ZGA genes, and found that most of these genes are indeed targets of the above transcription factors. Remarkably, more than a third of these ZGA genes appear to be regulated both by OrphNRs such as NR5A2 and by OBOX proteins, whose motifs co-occur in SINE B1 retrotransposable elements, which are enriched near ZGA genes. We propose that ZGA in mice is activated by recruitment of multiple transcription factors to SINE B1 elements that function as enhancers, and discuss a potential relevance of this mechanism to Alu retrotransposable elements in human ZGA.

ERVcancer: a web resource designed for querying activation of human endogenous retroviruses across major cancer types

Human endogenous retroviruses (HERVs) comprise approximately 8% of the human genome, co-opted into the dynamic regulatory network of cellular potency in early embryonic development. In recent studies, resurgent HERVs' transcriptional activity has been frequently observed in many types of human cancers, suggesting their potential functions in the occurrence and progression of malignancy. However, a dedicated web resource for querying the relationship between activation of HERVs and cancer development is lacking. Here, we have constructed a database to explore the sequence information, expression profiles, survival prognosis, and genetic interactions of HERVs in diverse cancer types. Our database currently contains RNA sequencing data of 580 HERVs across 16246 samples, including that of 6478 tumoral and 634 normal tissues, 932 cancer cell lines, as well as 151 early embryonic and 8051 human adult tissues. The primary goal is to provide an easily accessible and user-friendly database for professionals in the fields of bioinformatics, pathology, pharmacology, and related areas, enabling them to efficiently screen the activity of HERVs of interest in normal and cancerous tissues and evaluate the clinical relevance. The ERVcancer database is available at http://kyuanlab.com/ervcancer/.

Adenovirus small E1A directs activation of Alu transcription at YAP/TEAD- and AP-1-bound enhancers through interactions with the EP400 chromatin remodeler

Alu retrotransposons, which form the largest family of mobile DNA elements in the human genome, have recently come to attention as a potential source of regulatory novelties, most notably by participating in enhancer function. Even though Alu transcription by RNA polymerase III is subjected to tight epigenetic silencing, their expression has long been known to increase in response to various types of stress, including viral infection. Here we show that, in primary human fibroblasts, adenovirus small e1a triggered derepression of hundreds of individual Alus by promoting TFIIIB recruitment by Alu-bound TFIIIC. Epigenome profiling revealed an e1a-induced decrease of H3K27 acetylation and increase of H3K4 monomethylation at derepressed Alus, making them resemble poised enhancers. The enhancer nature of e1a-targeted Alus was confirmed by the enrichment, in their upstream regions, of the EP300/CBP acetyltransferase, EP400 chromatin remodeler and YAP1 and FOS transcription factors. The physical interaction of e1a with EP400 was critical for Alu derepression, which was abrogated upon EP400 ablation. Our data suggest that e1a targets a subset of enhancer Alus whose transcriptional activation, which requires EP400 and is mediated by the e1a-EP400 interaction, may participate in the manipulation of enhancer activity by adenoviruses.

Retrotransposons and Diabetes Mellitus

Retrotransposons are invasive genetic elements, which replicate by copying and pasting themselves throughout the genome in a process called retrotransposition. The most abundant retrotransposons by number in the human genome are Alu and LINE-1 elements, which comprise approximately 40% of the human genome. The ability of retrotransposons to expand and colonize eukaryotic genomes has rendered them evolutionarily successful and is responsible for creating genetic alterations leading to significant impacts on their hosts. Previous research suggested that hypomethylation of Alu and LINE-1 elements is associated with global hypomethylation and genomic instability in several types of cancer and diseases, such as neurodegenerative diseases, obesity, osteoporosis, and diabetes mellitus (DM). With the advancement of sequencing technologies and computational tools, the study of the retrotransposon's association with physiology and diseases is becoming a hot topic among researchers. Quantifying Alu and LINE-1 methylation is thought to serve as a surrogate measurement of global DNA methylation level. Although Alu and LINE-1 hypomethylation appears to serve as a cellular senescence biomarker promoting genomic instability, there is sparse information available regarding their potential functional and biological significance in DM. This review article summarizes the current knowledge on the involvement of the main epigenetic alterations in the methylation status of Alu and LINE-1 retrotransposons and their potential role as epigenetic markers of global DNA methylation in the pathogenesis of DM.

A primate-specific endogenous retroviral envelope protein sequesters SFRP2 to regulate human cardiomyocyte development

Endogenous retroviruses (ERVs) occupy a significant part of the human genome, with some encoding proteins that influence the immune system or regulate cell-cell fusion in early extra-embryonic development. However, whether ERV-derived proteins regulate somatic development is unknown. Here, we report a somatic developmental function for the primate-specific ERVH48-1 (SUPYN/Suppressyn). ERVH48-1 encodes a fragment of a viral envelope that is expressed during early embryonic development. Loss of ERVH48-1 led to impaired mesoderm and cardiomyocyte commitment and diverted cells to an ectoderm-like fate. Mechanistically, ERVH48-1 is localized to sub-cellular membrane compartments through a functional N-terminal signal peptide and binds to the WNT antagonist SFRP2 to promote its polyubiquitination and degradation, thus limiting SFRP2 secretion and blocking repression of WNT/β-catenin signaling. Knockdown of SFRP2 or expression of a chimeric SFRP2 with the ERVH48-1 signal peptide rescued cardiomyocyte differentiation. This study demonstrates how ERVH48-1 modulates WNT/β-catenin signaling and cell type commitment in somatic development.

A telomere-to-telomere cotton genome assembly reveals centromere evolution and a Mutator transposon-linked module regulating embryo development

Assembly of complete genomes can reveal functional genetic elements missing from draft sequences. Here we present the near-complete telomere-to-telomere and contiguous genome of the cotton species Gossypium raimondii. Our assembly identified gaps and misoriented or misassembled regions in previous assemblies and produced 13 centromeres, with 25 chromosomal ends having telomeres. In contrast to satellite-rich Arabidopsis and rice centromeres, cotton centromeres lack phased CENH3 nucleosome positioning patterns and probably evolved by invasion from long terminal repeat retrotransposons. In-depth expression profiling of transposable elements revealed a previously unannotated DNA transposon (MuTC01) that interacts with miR2947 to produce trans-acting small interfering RNAs (siRNAs), one of which targets the newly evolved LEC2 (LEC2b) to produce phased siRNAs. Systematic genome editing experiments revealed that this tripartite module, miR2947-MuTC01-LEC2b, controls the morphogenesis of complex folded embryos characteristic of Gossypium and its close relatives in the cotton tribe. Our study reveals a trans-acting siRNA-based tripartite regulatory pathway for embryo development in higher plants.

The neoantigens derived from transposable elements - A hidden treasure for cancer immunotherapy

Neoantigen-based therapy is a promising approach that selectively activates the immune system of the host to recognize and eradicate cancer cells. Preliminary clinical trials have validated the feasibility, safety, and immunogenicity of personalized neoantigen-directed vaccines, enhancing their effectiveness and broad applicability in immunotherapy. While many ongoing oncological trials concentrate on neoantigens derived from mutations, these targets do not consistently provoke an immune response in all patients harboring the mutations. Additionally, tumors like ovarian cancer, which have a low tumor mutational burden (TMB), may be less amenable to mutation-based neoantigen therapies. Recent advancements in next-generation sequencing and bioinformatics have uncovered a rich source of neoantigens from non-canonical RNAs associated with transposable elements (TEs). Considering the substantial presence of TEs in the human genome and the proven immunogenicity of TE-derived neoantigens in various tumor types, this review investigates the latest findings on TE-derived neoantigens, examining their clinical implications, challenges, and unique advantages in enhancing tumor immunotherapy.

Targeting transposable elements in cancer: developments and opportunities

Transposable elements (TEs), comprising nearly 50% of the human genome, have transitioned from being perceived as "genomic junk" to key players in cancer progression. Contemporary research links TE regulatory disruptions with cancer development, underscoring their therapeutic potential. Advances in long-read sequencing, computational analytics, single-cell sequencing, proteomics, and CRISPR-Cas9 technologies have enriched our understanding of TEs' clinical implications, notably their impact on genome architecture, gene regulation, and evolutionary processes. In cancer, TEs, including long interspersed element-1 (LINE-1), Alus, and long terminal repeat (LTR) elements, demonstrate altered patterns, influencing both tumorigenic and tumor-suppressive mechanisms. TE-derived nucleic acids and tumor antigens play critical roles in tumor immunity, bridging innate and adaptive responses. Given their central role in oncology, TE-targeted therapies, particularly through reverse transcriptase inhibitors and epigenetic modulators, represent a novel avenue in cancer treatment. Combining these TE-focused strategies with existing chemotherapy or immunotherapy regimens could enhance efficacy and offer a new dimension in cancer treatment. This review delves into recent TE detection advancements, explores their multifaceted roles in tumorigenesis and immune regulation, discusses emerging diagnostic and therapeutic approaches centered on TEs, and anticipates future directions in cancer research.

Genomic analyses of intricate interaction of TE-lncRNA overlapping genes with miRNAs in human diseases

BACKGROUND

Transposable elements (TEs) are known to be inserted into genome to create transcript isoforms or to generate long non-coding RNA (lncRNA) sequences. The insertion of TEs generates a gene protein sequence within the genome, but also provides a microRNA (miRNA) regulatory region.

OBJECTIVE

To determine the effect of gene sequence changes caused by TE insertion on miRNA binding and to investigate the formation of an overlapping lncRNA that represses it.

METHODS

The distribution of overlapping regions between exons and TE regions with lncRNA was examined using the Bedtools. miRNAs that can bind to those overlapping regions were identified through the miRDB web program. For TE-lncRNA overlapping genes, bioinformatic analysis was conducted using DAVID web database. Differential expression analysis was conducted using data from the GEO dataset and TCGA.

RESULTS

Most TEs were distributed more frequently in untranslated regions than open reading frames. There were 30 annotated TE-lncRNA overlapping genes with same strand that could bind to the same miRNA. As a result of identifying the association between these 30 genes and diseases, TGFB2, FCGR2A, DCTN5, and IFI6 were associated with breast cancer, and HMGCS1, FRMD4A, EDNRB, and SNCA were associated with Alzheimer's disease. Analysis of the GEO and TCGA data showed that the relevant expression of miR-891a and miR-28, which bind to the TE overlapping region of DCTN5 and HMGCS1, decreased.

CONCLUSION

This study indicates that the interaction between TE-lncRNA overlapping genes and miRNAs can affect disease progression.

Regulatory transposable elements in the encyclopedia of DNA elements

Transposable elements (TEs) comprise ~50% of our genome, but knowledge of how TEs affect genome evolution remains incomplete. Leveraging ENCODE4 data, we provide the most comprehensive study to date of TE contributions to the regulatory genome. We find 236,181 (~25%) human candidate cis-regulatory elements (cCREs) are TE-derived, with over 90% lineage-specific since the human-mouse split, accounting for 8-36% of lineage-specific cCREs. Except for SINEs, cCRE-associated transcription factor (TF) motifs in TEs are derived from ancestral TE sequence more than expected by chance. We show that TEs may adopt similar regulatory activities of elements near their integration site. Since human-mouse divergence, TEs have contributed 3-56% of TF binding site turnover events across 30 examined TFs. Finally, TE-derived cCREs are similar to non-TE cCREs in terms of MPRA activity and GWAS variant enrichment. Overall, our results substantiate the notion that TEs have played an important role in shaping the human regulatory genome.

Histone H3.3 lysine 9 and 27 control repressive chromatin at cryptic enhancers and bivalent promoters

Histone modifications are associated with distinct transcriptional states, but it is unclear whether they instruct gene expression. To investigate this, we mutate histone H3.3 K9 and K27 residues in mouse embryonic stem cells (mESCs). Here, we find that H3.3K9 is essential for controlling specific distal intergenic regions and for proper H3K27me3 deposition at promoters. The H3.3K9A mutation resulted in decreased H3K9me3 at regions encompassing endogenous retroviruses and induced a gain of H3K27ac and nascent transcription. These changes in the chromatin environment unleash cryptic enhancers, resulting in the activation of distinctive transcriptional programs and culminating in protein expression normally restricted to specialized immune cell types. The H3.3K27A mutant disrupts the deposition and spreading of the repressive H3K27me3 mark, particularly impacting bivalent genes with higher basal levels of H3.3 at promoters. Therefore, H3.3K9 and K27 crucially orchestrate repressive chromatin states at cis-regulatory elements and bivalent promoters, respectively, and instruct proper transcription in mESCs.

DUSP6 protein action and related hub genes prevention of sepsis-induced lung injury were screened by WGCNA and Venn

During a sepsis infection, the lung is extremely susceptible to damage. A condition known as acute respiratory distress syndrome (ARDS) may develop in extreme circumstances. The primary objective of this research is to identify important genes that are related with both sepsis and lung injury. These genes have the potential to act as novel biomarkers in the investigation of sepsis-induced lung injury prevention strategies. It was possible to download from GEO data both the sepsis-related dataset (GSE64457) and the lung injury-related dataset (GSE40839). In the GSE64457 dataset, using the "limma" package in R revealed 429 differentially expressed genes (DEGs) with logFC values more than or equal to -1 and p values <0.05. There were 266 genes that were up-regulated and 163 genes that were down-regulated. Through the use of Gene Ontology (GO), it was discovered that the majority of the DEGs were associated with the inflammatory response (BP terms), a particular granule lumen (CC terms), and protein binding (MF terms). By doing a pathway enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG), researchers were able to identify DEGs that were mostly associated with the NOD-like receptor signalling pathway, the TNF signalling pathway, and Epstein-Barr virus infection. Within the GSE40839 dataset, Weighted Gene Co-Expression Network Analysis (WGCNA) yielded a total of 7 modules, from which it was possible to screen out 2 critical modules and 693 key genes. The important genes and DEGs were both subjected to a Venn analysis. Finally, 14 genes that overlapped (ARL4A, LAIR1, MTHFD2, TSPAN13, DUSP6, PECR, CBS, TES, ASNS, SYNE1, FGF13, LCN2, KLF10, BCAT1) were closely associated to the incidence and development of sepsis-induced lung injury. This indicates that these genes are the essential genes to avoid the occurrence of sepsis-induced lung injury. This study provides novel strategies for preventing lung harm brought on by sepsis.

Single-cell chromatin accessibility and transposable element landscapes reveal shared features of tissue-residing immune cells

Tissue adaptation is required for regulatory T (Treg) cell function within organs. Whether this program shares aspects with other tissue-localized immune populations is unclear. Here, we analyzed single-cell chromatin accessibility data, including the transposable element (TE) landscape of CD45+ immune cells from colon, skin, adipose tissue, and spleen. We identified features of organ-specific tissue adaptation across different immune cells. Focusing on tissue Treg cells, we found conservation of the Treg tissue adaptation program in other tissue-localized immune cells, such as amphiregulin-producing T helper (Th)17 cells. Accessible TEs can act as regulatory elements, but their contribution to tissue adaptation is not understood. TE landscape analysis revealed an enrichment of specific transcription factor binding motifs in TE regions within accessible chromatin peaks. TEs, specifically from the LTR family, were located in enhancer regions and associated with tissue adaptation. These findings broaden our understanding of immune tissue residency and provide an important step toward organ-specific immune interventions.

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

Retro-age: A unique epigenetic biomarker of aging captured by DNA methylation states of retroelements

Reactivation of retroelements in the human genome has been linked to aging. However, whether the epigenetic state of specific retroelements can predict chronological age remains unknown. We provide evidence that locus-specific retroelement DNA methylation can be used to create retroelement-based epigenetic clocks that accurately measure chronological age in the immune system, across human tissues, and pan-mammalian species. We also developed a highly accurate retroelement epigenetic clock compatible with EPICv.2.0 data that was constructed from CpGs that did not overlap with existing first- and second-generation epigenetic clocks, suggesting a unique signal for epigenetic clocks not previously captured. We found retroelement-based epigenetic clocks were reversed during transient epigenetic reprogramming, accelerated in people living with HIV-1, and responsive to antiretroviral therapy. Our findings highlight the utility of retroelement-based biomarkers of aging and support a renewed emphasis on the role of retroelements in geroscience.

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

Gapped-kmer sequence modeling robustly identifies regulatory vocabularies and distal enhancers conserved between evolutionarily distant mammals

Gene regulatory elements drive complex biological phenomena and their mutations are associated with common human diseases. The impacts of human regulatory variants are often tested using model organisms such as mice. However, mapping human enhancers to conserved elements in mice remains a challenge, due to both rapid enhancer evolution and limitations of current computational methods. We analyze distal enhancers across 45 matched human/mouse cell/tissue pairs from a comprehensive dataset of DNase-seq experiments, and show that while cell-specific regulatory vocabulary is conserved, enhancers evolve more rapidly than promoters and CTCF binding sites. Enhancer conservation rates vary across cell types, in part explainable by tissue specific transposable element activity. We present an improved genome alignment algorithm using gapped-kmer features, called gkm-align, and make genome wide predictions for 1,401,803 orthologous regulatory elements. We show that gkm-align discovers 23,660 novel human/mouse conserved enhancers missed by previous algorithms, with strong evidence of conserved functional activity.

Evolutionary patterns and functional effects of 3D chromatin structures in butterflies with extensive genome rearrangements

Chromosome rearrangements may distort 3D chromatin architectures and thus change gene regulation, yet how 3D chromatin structures evolve in insects is largely unknown. Here, we obtain chromosome-level genomes for four butterfly species, Graphium cloanthus, Graphium sarpedon, Graphium eurypylus with 2n = 30, 40, and 60, respectively, and Papilio bianor with 2n = 60. Together with large-scale Hi-C data, we find that inter-chromosome rearrangements very rarely disrupted the pre-existing 3D chromatin structure of ancestral chromosomes. However, some intra-chromosome rearrangements changed 3D chromatin structures compared to the ancestral configuration. We find that new TADs and subTADs have emerged across the rearrangement sites where their adjacent compartments exhibit uniform types. Two intra-chromosome rearrangements altered Rel and lft regulation, potentially contributing to wing patterning differentiation and host plant choice. Notably, butterflies exhibited chromatin loops between Hox gene cluster ANT-C and BX-C, unlike Drosophila. Our CRISPR-Cas9 experiments in butterflies confirm that knocking out the CTCF binding site of the loops in BX-C affected the phenotypes regulated by Antp in ANT-C, resulting in legless larva. Our results reveal evolutionary patterns of insect 3D chromatin structures and provide evidence that 3D chromatin structure changes can play important roles in the evolution of traits.

RNaseH-based ribodepletion of total planarian RNA improves detection of longer and non-polyadenylated transcripts

The overwhelming majority of RNA species isolated from cells or tissues using organic extraction are ribosomal RNAs (rRNA), whereas a relatively small percentage are messenger RNAs (mRNA). For studies that seek to detect mRNA transcripts and measure changes in their expression, this lopsided ratio of desired transcripts to undesired transcripts creates a significant challenge to obtaining sensitive and reproducible results. One method for improving mRNA detection is to selectively amplify polyadenylated (polyA) mRNA molecules when generating RNA-seq libraries, a strategy that is generally very successful in many species. However, this strategy is less effective when starting with total RNA from some species e.g., the planarian species Schmidtea mediterranea (S.med), as it generates libraries that still contain significant and variable amounts of rRNA reads. Further, commercially available ribodepletion kits do not efficiently deplete rRNAs from these samples because their sequences are divergent from mammalian rRNAs. Here we report a customized, optimized, and economical ribodepletion strategy than allows the generation of comprehensive RNA-seq libraries with less than one percent rRNA contamination. We show that this method improves transcript detection, particularly for those without polyA tails (e.g., core histones) and those that are relatively long (e.g., microtubule motor proteins). Using this custom ribodepletion approach, we also detected many transcripts that are not represented in the most recent set of S.med gene annotations, including a subset that are likely expressed transposable elements (TEs). To facilitate future differential expression analyses of these newly identified loci, we created both an annotation file of the new loci we identified and a bioinformatic pipeline for generating additional annotations from future libraries. As significant recent research shows that TE activation is regulated and functionally important, the resources provided here will provide a starting point for investigating such mechanisms in planarians and other species with less conserved rRNA sequences.

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

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

TEENA: an integrated web server for transposable element enrichment analysis in various model and non-model organisms

Transposable elements (TEs) are abundant in the genomes of various eukaryote organisms. Increasing evidence suggests that TEs can play crucial regulatory roles-usually by creating cis-elements (e.g. enhancers and promoters) bound by distinct transcription factors (TFs). TE-derived cis-elements have gained unprecedented attentions recently, and one key step toward their understanding is to identify the enriched TEs in distinct genomic intervals (e.g. a set of enhancers or TF binding sites) as candidates for further study. Nevertheless, such analysis remains challenging for researchers unfamiliar with TEs or lack strong bioinformatic skills. Here, we present TEENA (Transposable Element ENrichment Analyzer) to streamline TE enrichment analysis in various organisms. It implements an optimized pipeline, hosts the genome/gene/TE annotations of almost one hundred species, and provides multiple parameters to enable its flexibility. Taking genomic interval data as the only user-supplied file, it can automatically retrieve the corresponding annotations and finish a routine analysis in a couple minutes. Multiple case studies demonstrate that it can produce highly reliable results matching previous knowledge. TEENA can be freely accessed at: https://sun-lab.yzu.edu.cn/TEENA. Due to its easy-to-use design, we expect it to facilitate the studies of the regulatory function of TEs in various model and non-model organisms.

Transposable elements-mediated recruitment of KDM1A epigenetically silences HNF4A expression to promote hepatocellular carcinoma

Transposable elements (TEs) contribute to gene expression regulation by acting as cis-regulatory elements that attract transcription factors and epigenetic regulators. This research aims to explore the functional and clinical implications of transposable element-related molecular events in hepatocellular carcinoma, focusing on the mechanism through which liver-specific accessible TEs (liver-TEs) regulate adjacent gene expression. Our findings reveal that the expression of HNF4A is inversely regulated by proximate liver-TEs, which facilitates liver cancer cell proliferation. Mechanistically, liver-TEs are predominantly occupied by the histone demethylase, KDM1A. KDM1A negatively influences the methylation of histone H3 Lys4 (H3K4) of liver-TEs, resulting in the epigenetic silencing of HNF4A expression. The suppression of HNF4A mediated by KDM1A promotes liver cancer cell proliferation. In conclusion, this study uncovers a liver-TE/KDM1A/HNF4A regulatory axis that promotes liver cancer growth and highlights KDM1A as a promising therapeutic target. Our findings provide insight into the transposable element-related molecular mechanisms underlying liver cancer progression.

A spatially resolved multiomic single-cell atlas of soybean development

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DOF11 motif that were co-up-regulated in late peripheral endosperm and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

LINE-1 transcription activates long-range gene expression

Long interspersed nuclear element-1 (LINE-1 or L1) is a retrotransposon group that constitutes 17% of the human genome and shows variable expression across cell types. However, the control of L1 expression and its function in gene regulation are incompletely understood. Here we show that L1 transcription activates long-range gene expression. Genome-wide CRISPR-Cas9 screening using a reporter driven by the L1 5' UTR in human cells identifies functionally diverse genes affecting L1 expression. Unexpectedly, altering L1 expression by knockout of regulatory genes impacts distant gene expression. L1s can physically contact their distal target genes, with these interactions becoming stronger upon L1 activation and weaker when L1 is silenced. Remarkably, L1s contact and activate genes essential for zygotic genome activation (ZGA), and L1 knockdown impairs ZGA, leading to developmental arrest in mouse embryos. These results characterize the regulation and function of L1 in long-range gene activation and reveal its importance in mammalian ZGA.

Precise Gene Knock-In Tools with Minimized Risk of DSBs: A Trend for Gene Manipulation

Gene knock-in refers to the insertion of exogenous functional genes into a target genome to achieve continuous expression. Currently, most knock-in tools are based on site-directed nucleases, which can induce double-strand breaks (DSBs) at the target, following which the designed donors carrying functional genes can be inserted via the endogenous gene repair pathway. The size of donor genes is limited by the characteristics of gene repair, and the DSBs induce risks like genotoxicity. New generation tools, such as prime editing, transposase, and integrase, can insert larger gene fragments while minimizing or eliminating the risk of DSBs, opening new avenues in the development of animal models and gene therapy. However, the elimination of off-target events and the production of delivery carriers with precise requirements remain challenging, restricting the application of the current knock-in treatments to mainly in vitro settings. Here, a comprehensive review of the knock-in tools that do not/minimally rely on DSBs and use other mechanisms is provided. Moreover, the challenges and recent advances of in vivo knock-in treatments in terms of the therapeutic process is discussed. Collectively, the new generation of DSBs-minimizing and large-fragment knock-in tools has revolutionized the field of gene editing, from basic research to clinical treatment.

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

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

Pioneer Transcription Factors: The First Domino in Zygotic Genome Activation

Zygotic genome activation (ZGA) is a pivotal event in mammalian embryogenesis, marking the transition from maternal to zygotic control of development. During the ZGA process that is characterized by the intricate cascade of gene expression, who tipped the first domino in a meticulously arranged sequence is a subject of paramount interest. Recently, Dux, Obox and Nr5a2 were identified as pioneer transcription factors that reside at the top of transcriptional hierarchy. Through co-option of retrotransposon elements as hubs for transcriptional activation, these pioneer transcription factors rewire the gene regulatory network, thus initiating ZGA. In this review, we provide a snapshot of the mechanisms underlying the functions of these pioneer transcription factors. We propose that ZGA is the starting point where the embryo's own genome begins to influence development trajectory, therefore in-depth dissecting the functions of pioneer transcription factors during ZGA will form a cornerstone of our understanding for early embryonic development, which will pave the way for advancing our grasp of mammalian developmental biology and optimizing in vitro production (IVP) techniques.

Structural variants contribute to phenotypic variation in maize

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants, especially transposable elements are likely to affect phenotypic variation but we need better methods in maize for detecting polymorphic structural variants and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic structural variants and then genotyped a large maize diversity panel for these variants using short-read sequencing data. We characterized variation of SVs within the panel and identified SV polymorphisms that are associated with life history traits and genotype-by-environment interactions. While most of the SVs associated with traits contained TEs, only one of the SV's boundaries clearly matched TE breakpoints indicative of a TE insertion, whereas the other polymorphisms were likely caused by deletions. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that this method did not identify variants that would have been missed in a SNP association study.

PR-SET7 epigenetically restrains uterine interferon response and cell death governing proper postnatal stromal development

The differentiation of the stroma is a hallmark event during postnatal uterine development. However, the spatiotemporal changes that occur during this process and the underlying regulatory mechanisms remain elusive. Here, we comprehensively delineated the dynamic development of the neonatal uterus at single-cell resolution and characterized two distinct stromal subpopulations, inner and outer stroma. Furthermore, single-cell RNA sequencing revealed that uterine ablation of Pr-set7, the sole methyltransferase catalyzing H4K20me1, led to a reduced proportion of the inner stroma due to massive cell death, thus impeding uterine development. By combining RNA sequencing and epigenetic profiling of H4K20me1, we demonstrated that PR-SET7-H4K20me1 either directly repressed the transcription of interferon stimulated genes or indirectly restricted the interferon response via silencing endogenous retroviruses. Declined H4K20me1 level caused viral mimicry responses and ZBP1-mediated apoptosis and necroptosis in stromal cells. Collectively, our study provides insight into the epigenetic machinery governing postnatal uterine stromal development mediated by PR-SET7.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

Integrated enhancer regulatory network by enhancer-promoter looping in gastric cancer

Enhancer cis-regulatory elements play critical roles in gene regulation at many stages of cell growth. Enhancers in cancer cells also regulate the transcription of oncogenes. In this study, we performed a comprehensive analysis of long-range chromatin interactions, histone modifications, chromatin accessibility and expression in two gastric cancer (GC) cell lines compared to normal gastric epithelial cells. We found that GC-specific enhancers marked by histone modifications can activate a population of genes, including some oncogenes, by interacting with their proximal promoters. In addition, motif analysis of enhancer-promoter interacting enhancers showed that GC-specific transcription factors are enriched. Among them, we found that MYB is crucial for GC cell growth and activated by the enhancer with an enhancer-promoter loop and TCF7 upregulation. Clinical GC samples showed epigenetic activation of enhancers at the MYB locus and significant upregulation of TCF7 and MYB, regardless of molecular GC subtype and clinicopathological factors. Single-cell RNA sequencing of gastric mucosa with intestinal metaplasia showed high expression of TCF7 and MYB in intestinal stem cells. When we inactivated the loop-forming enhancer at the MYB locus using CRISPR interference (dCas9-KRAB), GC cell growth was significantly inhibited. In conclusion, we identified MYB as an oncogene activated by a loop-forming enhancer and contributing to GC cell growth.

Insertion of short L1 sequences generates inter-strain histone acetylation differences in the mouse

BACKGROUND

Gene expression divergence between populations and between individuals can emerge from genetic variations within the genes and/or in the cis regulatory elements. Since epigenetic modifications regulate gene expression, it is conceivable that epigenetic variations in cis regulatory elements can also be a source of gene expression divergence.

RESULTS

In this study, we compared histone acetylation (namely, H3K9ac) profiles in two mouse strains of different subspecies origin, C57BL/6 J (B6) and MSM/Ms (MSM), as well as their F1 hybrids. This identified 319 regions of strain-specific acetylation, about half of which were observed between the alleles of F1 hybrids. While the allele-specific presence of the interferon regulatory factor 3 (IRF3) binding sequence was associated with allele-specific histone acetylation, we also revealed that B6-specific insertions of a short 3' fragment of LINE-1 (L1) retrotransposon occur within or proximal to MSM-specific acetylated regions. Furthermore, even in hyperacetylated domains, flanking regions of non-polymorphic 3' L1 fragments were hypoacetylated, suggesting a general activity of the 3' L1 fragment to induce hypoacetylation. Indeed, we confirmed the binding of the 3' region of L1 by three Krüppel-associated box domain-containing zinc finger proteins (KZFPs), which interact with histone deacetylases. These results suggest that even a short insertion of L1 would be excluded from gene- and acetylation-rich regions by natural selection. Finally, mRNA-seq analysis for F1 hybrids was carried out, which disclosed a link between allele-specific promoter/enhancer acetylation and gene expression.

CONCLUSIONS

This study disclosed a number of genetic changes that have changed the histone acetylation levels during the evolution of mouse subspecies, a part of which is associated with gene expression changes. Insertions of even a very short L1 fragment can decrease the acetylation level in their neighboring regions and thereby have been counter-selected in gene-rich regions, which may explain a long-standing mystery of discrete genomic distribution of LINEs and SINEs.

Single-cell transcriptome profiling reveals the spatiotemporal distribution of triterpenoid saponin biosynthesis and transposable element activity in Gynostemma pentaphyllum shoot apexes and leaves

Gynostemma pentaphyllum (Thunb.) Makino is an important producer of dammarene-type triterpenoid saponins. These saponins (gypenosides) exhibit diverse pharmacological benefits such as anticancer, antidiabetic, and immunomodulatory effects, and have major potential in the pharmaceutical and health care industries. Here, we employed single-cell RNA sequencing (scRNA-seq) to profile the transcriptomes of more than 50,000 cells derived from G. pentaphyllum shoot apexes and leaves. Following cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Each cell type displayed substantial transcriptomic heterogeneity both within and between tissues. Examining gene expression patterns across various cell types revealed that gypenoside biosynthesis predominantly occurred in mesophyll cells, with heightened activity observed in shoot apexes compared to leaves. Furthermore, we explored the impact of transposable elements (TEs) on G. pentaphyllum transcriptomic landscapes. Our findings the highlighted the unbalanced expression of certain TE families across different cell types in shoot apexes and leaves, marking the first investigation of TE expression at the single-cell level in plants. Additionally, we observed dynamic expression of genes involved in gypenoside biosynthesis and specific TE families during epidermal and vascular cell development. The involvement of TE expression in regulating cell differentiation and gypenoside biosynthesis warrant further exploration. Overall, this study not only provides new insights into the spatiotemporal organization of gypenoside biosynthesis and TE activity in G. pentaphyllum shoot apexes and leaves but also offers valuable cellular and genetic resources for a deeper understanding of developmental and physiological processes at single-cell resolution in this species.

Endogenous retroelement expression in the gut microenvironment of people living with HIV-1

BACKGROUND

Endogenous retroelements (EREs), including human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs), comprise almost half of the human genome. Our previous studies of the interferome in the gut suggest potential mechanisms regarding how IFNb may drive HIV-1 gut pathogenesis. As ERE activity is suggested to partake in type 1 immune responses and is incredibly sensitive to viral infections, we sought to elucidate underlying interactions between ERE expression and gut dynamics in people living with HIV-1 (PLWH).

METHODS

ERE expression profiles from bulk RNA sequencing of colon biopsies and PBMC were compared between a cohort of PLWH not on antiretroviral therapy (ART) and uninfected controls.

FINDINGS

59 EREs were differentially expressed in the colon of PLWH when compared to uninfected controls (padj <0.05 and FC ≤ -1 or ≥ 1) [Wald's Test]. Of these 59, 12 EREs were downregulated in PLWH and 47 were upregulated. Colon expression of the ERE loci LTR19_12p13.31 and L1FLnI_1q23.1s showed significant correlations with certain gut immune cell subset frequencies in the colon. Furthermore L1FLnI_1q23.1s showed a significant upregulation in peripheral blood mononuclear cells (PBMCs) of PLWH when compared to uninfected controls suggesting a common mechanism of differential ERE expression in the colon and PBMC.

INTERPRETATION

ERE activity has been largely understudied in genomic characterizations of human pathologies. We show that the activity of certain EREs in the colon of PLWH is deregulated, supporting our hypotheses that their underlying activity could function as (bio)markers and potential mediators of pathogenesis in HIV-1 reservoirs.

FUNDING

US NIH grants NCI CA260691 (DFN) and NIAID UM1AI164559 (DFN).

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BACKGROUND

Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.

RESULTS

Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.

CONCLUSIONS

Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

The Diverse Evolutionary Histories of Domesticated Metaviral Capsid Genes in Mammals

Selfish genetic elements comprise significant fractions of mammalian genomes. In rare instances, host genomes domesticate segments of these elements for function. Using a complete human genome assembly and 25 additional vertebrate genomes, we re-analyzed the evolutionary trajectories and functional potential of capsid (CA) genes domesticated from Metaviridae, a lineage of retrovirus-like retrotransposons. Our study expands on previous analyses to unearth several new insights about the evolutionary histories of these ancient genes. We find that at least five independent domestication events occurred from diverse Metaviridae, giving rise to three universally retained single-copy genes evolving under purifying selection and two gene families unique to placental mammals, with multiple members showing evidence of rapid evolution. In the SIRH/RTL family, we find diverse amino-terminal domains, widespread loss of protein-coding capacity in RTL10 despite its retention in several mammalian lineages, and differential utilization of an ancient programmed ribosomal frameshift in RTL3 between the domesticated CA and protease domains. Our analyses also reveal that most members of the PNMA family in mammalian genomes encode a conserved putative amino-terminal RNA-binding domain (RBD) both adjoining and independent from domesticated CA domains. Our analyses lead to a significant correction of previous annotations of the essential CCDC8 gene. We show that this putative RBD is also present in several extant Metaviridae, revealing a novel protein domain configuration in retrotransposons. Collectively, our study reveals the divergent outcomes of multiple domestication events from diverse Metaviridae in the common ancestor of placental mammals.

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

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

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

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

Population genetic analysis based on the polymorphisms mediated by transposons in the genomes of pig

Transposable elements (TEs) mobility is capable of generating a large number of structural variants (SVs), which can have considerable potential as molecular markers for genetic analysis and molecular breeding in livestock. Our results showed that the pig genome contains mainly TE-SVs generated by short interspersed nuclear elements (51,873/76.49%), followed by long interspersed nuclear elements (11,131/16.41%), and more than 84% of the common TE-SVs (Minor allele frequency, MAF > 0.10) were validated to be polymorphic. Subsequently, we utilized the identified TE-SVs to gain insights into the population structure, resulting in clear differentiation among the three pig groups and facilitating the identification of relationships within Chinese local pig breeds. In addition, we investigated the frequencies of TEs in the gene coding regions of different pig groups and annotated the respective TE types, related genes, and functional pathways. Through genome-wide comparisons of Large White pigs and Chinese local pigs utilizing the Beijing Black pigs, we identified TE-mediated SVs associated with quantitative trait loci and observed that they were mainly involved in carcass traits and meat quality traits. Lastly, we present the first documented evidence of TE transduction in the pig genome.

Endogenous retroviruses in development and health

Endogenous retroviruses (ERVs) are evolutionary remnants of retroviral infections in which the viral genome became embedded as a dormant regulatory element within the host germline. When ERVs become activated, they comprehensively rewire genomic regulatory networks of the host and facilitate critical developmental events, such as preimplantation development and placentation, in a manner specific to species, developmental stage, and tissues. However, accumulating evidence suggests that aberrant ERV transcription compromises genome stability and has been implicated in cellular senescence and various pathogenic processes, underscoring the significance of host genomic surveillance mechanisms. Here, we revisit the prominent functions of ERVs in early development and highlight their emerging roles in mammalian post-implantation development and organogenesis. We also discuss their implications for aging and pathological processes such as microbial infection, immune response. Furthermore, we discuss recent advances in stem-cell-based models, single-cell omics, and genome editing technologies, which serve as beacons illuminating the versatile nature of ERVs in mammalian development and health.

Condensin-mediated restriction of retrotransposable elements facilitates brain development in Drosophila melanogaster

Neural stem and progenitor cell (NSPC) maintenance is essential for ensuring that organisms are born with proper brain volumes and head sizes. Microcephaly is a disorder in which babies are born with significantly smaller head sizes and cortical volumes. Mutations in subunits of the DNA organizing complex condensin have been identified in microcephaly patients. However, the molecular mechanisms by which condensin insufficiency causes microcephaly remain elusive. We previously identified conserved roles for condensins in repression of retrotransposable elements (RTEs). Here, we show that condensin subunit knockdown in NSPCs of the Drosophila larval central brain increases RTE expression and mobility which causes cell death, and significantly decreases adult head sizes and brain volumes. These findings suggest that unrestricted RTE expression and activity may lead to improper brain development in condensin insufficient organisms, and lay the foundation for future exploration of causative roles for RTEs in other microcephaly models.

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

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

ADNP modulates SINE B2-derived CTCF-binding sites during blastocyst formation in mice

CTCF is crucial for chromatin structure and transcription regulation in early embryonic development. However, the kinetics of CTCF chromatin occupation in preimplantation embryos have remained unclear. In this study, we used CUT&RUN technology to investigate CTCF occupancy in mouse preimplantation development. Our findings revealed that CTCF begins binding to the genome prior to zygotic genome activation (ZGA), with a preference for CTCF-anchored chromatin loops. Although the majority of CTCF occupancy is consistently maintained, we identified a specific set of binding sites enriched in the mouse-specific short interspersed element (SINE) family B2 that are restricted to the cleavage stages. Notably, we discovered that the neuroprotective protein ADNP counteracts the stable association of CTCF at SINE B2-derived CTCF-binding sites. Knockout of Adnp in the zygote led to impaired CTCF binding signal recovery, failed deposition of H3K9me3, and transcriptional derepression of SINE B2 during the morula-to-blastocyst transition, which further led to unfaithful cell differentiation in embryos around implantation. Our analysis highlights an ADNP-dependent restriction of CTCF binding during cell differentiation in preimplantation embryos. Furthermore, our findings shed light on the functional importance of transposable elements (TEs) in promoting genetic innovation and actively shaping the early embryo developmental process specific to mammals.

Activation of human endogenous retroviruses and its physiological consequences

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

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

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

Navigating the brain and aging: exploring the impact of transposable elements from health to disease

Transposable elements (TEs) are mobile genetic elements that constitute on average 45% of mammalian genomes. Their presence and activity in genomes represent a major source of genetic variability. While this is an important driver of genome evolution, TEs can also have deleterious effects on their hosts. A growing number of studies have focused on the role of TEs in the brain, both in physiological and pathological contexts. In the brain, their activity is believed to be important for neuronal plasticity. In neurological and age-related disorders, aberrant activity of TEs may contribute to disease etiology, although this remains unclear. After providing a comprehensive overview of transposable elements and their interactions with the host, this review summarizes the current understanding of TE activity within the brain, during the aging process, and in the context of neurological and age-related conditions.