The majority of primate-specific regulatory sequences are derived from transposable elements

The latent cis-regulatory potential of mobile DNA in Escherichia coli

Transposable elements can alter gene regulation in their host genome, either when they integrate into a genome, or when they accrue mutations after integration. However, the extent to which transposons can alter gene expression, as well as the necessary mutational steps, are not well characterized. Here we study the gene regulatory potential of the prominent IS3 family of transposable elements in E.coli. We started with 10 sequences from the ends of 5 IS3 sequences, created 18,537 random mutations in them, and measured their promoter activity using a massively parallel reporter assay. All 10 sequences could evolve de-novo promoter activity from single point mutations. De-novo promoters mostly emerge from existing proto-promoter sequences when mutations create new -10 boxes downstream of preexisting -35 boxes. The ends of IS3s harbor ~1.5 times as many such proto-promoter sequences than the E. coli genome. We also estimate that at least 26% of the 706 characterized IS3s already encode promoters. Our study shows that transposable elements can have a high latent cis-regulatory potential. Our observations can help to explain why mobile DNA may persist in prokaryotic genomes. They also underline the potential use of transposable elements as a substrate for evolving new gene expression.

Binding of NF-Y to transposable elements in mouse and human cells

BACKGROUND

Transposable Elements (TEs) represent a sizeable amount of mammalian genomes, providing regulatory sequences involved in shaping gene expression patterns. NF-Y is a Transcription factor -TF- trimer that binds to the CCAAT box, belonging to a selected group implicated in determining initiation of coding and noncoding RNAs.

RESULTS

We focus on NF-Y TE locations in 8 human and 8 mouse cells. Binding is exclusive for retroviral LTR12, MLT1 and MER in human and RLTR10 and IAPLTR in mouse cells. Cobinding and analysis of the DNA matrices signal enrichment of distinct TFs neighboring CCAAT in the three TE classes: MAFK/F/G in LTR12 and USF1/2 in MLT1 with precise alignment of sites, PKNOX1, MEIS2, PBX2/3 TALE TFs in MER57. The presence of "epigenetic" marks in human cells indicate prevalent co-association with open chromatin in MER, closed in LTR12 and mixed in MLT1. Based on chromatin features, these locations are mostly marked as enhancers, as confirmed by analysis of loci predicted to generate eRNAs.

CONCLUSIONS

These results are discussed in the context of functional data, suggesting a complex -positive and potentially-negative role of NF-Y on distinct classes of repetitive sequences.

Transposable elements as instructors of the immune system

Transposable elements (TEs) are mobile repetitive nucleic acid sequences that have been incorporated into the genome through spontaneous integration, accounting for almost 50% of human DNA. Even though most TEs are no longer mobile today, studies have demonstrated that they have important roles in different biological processes, such as ageing, embryonic development, and cancer. TEs influence these processes through various mechanisms, including active transposition of TEs contributing to ongoing evolution, transposon transcription generating RNA or protein, and by influencing gene regulation as enhancers. However, how TEs interact with the immune system remains a largely unexplored field. In this Perspective, we describe how TEs might influence different aspects of the immune system, such as innate immune responses, T cell activation and differentiation, and tissue adaptation. Furthermore, TEs can serve as a source of neoantigens for T cells in antitumour immunity. We suggest that TE biology is an important emerging field of immunology and discuss the potential to harness the TE network therapeutically, for example, to improve immunotherapies for cancer and autoimmune and inflammatory diseases.

Expression spectrum of TE-derived transcripts in human adult tissues

Transposable elements (TEs) are vital components of eukaryotic genomes and have played a critical role in genome evolution. Although most TEs are silenced in the mammalian genome, increasing evidence suggests that certain TEs are actively involved in gene regulation during early developmental stages. However, the extent to which human TEs drive gene transcription in adult tissues remains largely unexplored. In this study, we systematically analyzed 17,329 human transcriptomes to investigate how TEs influence gene transcription across 47 adult tissues. Our findings reveal that TE-derived transcripts are broadly expressed in human tissues, contributing to both housekeeping functions and tissue-specific gene regulation. We identified sex-specific expression of TE-derived transcripts regulated by sex hormones in breast tissue between females and males. Our results demonstrated that TE-derived alternative transcription initiation significantly enhances the variety of translated protein products, e.g., changes in the N-terminal peptide length of WNT2B caused by TE-derived transcription result in isoform-specific subcellular localization. Additionally, we identified 68 human-specific TE-derived transcripts associated with metabolic processes and environmental adaptation. Together, these findings highlight the pivotal evolutionary role of TEs in shaping the human transcriptome, demonstrating how conserved and human-specific TEs contribute to transcriptional and translational innovation in human genome evolution.

Derepression of transposable elements in mouse prefrontal cortex disrupts social behavior

Here, we present a synthetic biology approach to assess the social behavioral consequences of altered function of the Krüppel-associated box zinc finger protein (KZFP) interacting protein TRIM28 within the prefrontal cortex (PFC) of male and female mice. We reprogrammed natural TRIM28 WT by replacing the transcriptionally repressive domain with an enhanced transcriptional activation domain VP64-p65-Rta (TRIM28 VPR ), or by excising the transcriptional regulatory domain (TRIM28 NFD ). In vitro validation confirmed that TRIM28 WT represses, and TRIM28 VPR activates, the expression of a KZFP-regulated luciferase reporter gene. Upon intra-PFC viral-mediated delivery of TRIM28 variants, we observed that inversion of TRIM28 transcriptional control via HSV-TRIM28 VPR reduced the salience of novel social interaction for male and female mice while not affecting non-social behaviors. RNA-sequencing revealed HSV-TRIM28 VPR promoted transcriptional escape of all classes of TEs, particularly those located within intronic and distal enhancer regions of downregulated immune genes. HSV-TRIM28 VPR -driven social deficits were reversible by intra-PFC repletion of interferon cytokines. These novel data point to PFC KZFP-TRIM28 interactions as necessary to stabilize TEs to enable cis-regulation of key immune gene expression and enhance organismal capacity for complex, pro-social behaviors.

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

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

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

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

Single chromatin fiber profiling and nucleosome position mapping in the human brain

We apply a single-molecule chromatin fiber sequencing (Fiber-seq) protocol designed for amplification-free cell-type-specific mapping of the regulatory architecture at nucleosome resolution along extended ∼10-kb chromatin fibers to neuronal and non-neuronal nuclei sorted from human brain tissue. Specifically, application of this method enables the resolution of cell-selective promoter and enhancer architectures on single fibers, including transcription factor footprinting and position mapping, with sequence-specific fixation of nucleosome arrays flanking transcription start sites and regulatory motifs. We uncover haplotype-specific chromatin patterns, multiple regulatory elements cis-aligned on individual fibers, and accessible chromatin at 20,000 unique sites encompassing retrotransposons and other repeat sequences hitherto "unmappable" by short-read epigenomic sequencing. Overall, we show that Fiber-seq is applicable to human brain tissue, offering sharp demarcation of nucleosome-depleted regions at sites of open chromatin in conjunction with multi-kilobase nucleosomal positioning at single-fiber resolution on a genome-wide scale.

Construction of the porcine genome mobile element variations and investigation of its role in population diversity and gene expression

BACKGROUND

Mobile element variants (MEVs) have a significant and complex impact on genomic diversity and phenotypic traits. However, the quantity, distribution, and relationship with gene expression and complex traits of MEVs in the pig genome remain poorly understood.

RESULTS

We constructed the most comprehensive porcine MEV library based on high-depth whole genome sequencing (WGS) data from 747 pigs across 59 breeds worldwide. This database identified a total of 147,993 polymorphic MEVs, including 121,099 short interspersed nuclear elements (SINEs), 26,053 long interspersed nuclear elements (LINEs), 802 long terminal repeats (LTRs), and 39 other transposons, among which 54% are newly discovered. We found that MEVs are unevenly distributed across the genome and are strongly influenced by negative selection effects. Importantly, we identified 514, 530, and 584 candidate MEVs associated with population differentiation, domestication, and breed formation, respectively. For example, a significantly differentiated MEV is located in the ATRX intron between Asian and European pigs, whereas ATRX is also differentially expressed between Asian and European pigs in muscle tissue. In addition, we identified 4,169 expressed MEVs (eMEVs) significantly associated with gene expression and 6,914 splicing MEVs (sMEVs) associated with gene splicing based on RNA-seq data from 266 porcine liver tissues. These eMEVs and sMEVs explain 6.24% and 9.47%, respectively, of the observed cis-heritability and highlight the important role of MEVs in the regulation of gene expression. Finally, we provide a high-quality SNP-MEV reference haplotype panel to impute MEV genotypes from genome-wide SNPs. Notably, we identified a candidate MEV significantly associated with total teat number, demonstrating the functionality of this reference panel.

CONCLUSIONS

The present investigation demonstrated the importance of MEVs in pigs in terms of population diversity, gene expression and phenotypic traits, which may provide useful resources and theoretical support for pig genetics and breeding.

Positive selection in cilia-related genes may facilitate deep-sea adaptation of Thermocollonia jamsteci

Deep-sea hydrothermal vents are characterized by high hydrostatic pressure, hypoxia, darkness and toxic substances. However, how organisms adapt to such extreme marine ecosystems remain poorly understood. We hypothesize that adaptive evolution plays an essential role in generating novelty for evolutionary adaptation to the deep-sea environment because adaptive evolution has been found to be critical for species origin and evolution. In this project, the chromosome-level genome of the deep-sea hydrothermal vent gastropod T. jamsteci was constructed for the first time to examine molecular mechanisms of its adaptation to the deep-sea environment. The genome size was large (2.54 Gb), ranking at the top of all species in the Vetigastropoda subclass, driven primarily by the bursts of transposable elements (TEs). The transposition of TEs may also trigger chromosomal changes including both inter-chromosomal fusions and intra-chromosomal activities involving chromosome inversions, rearrangements and fusions, as revealed by comparing the genomes of T. jamsteci and its closely related shallow-sea species Gibbula magus. Innovative changes including the expansion of the ABC transporter gene family that may facilitate detoxification, duplication of genes related to endocytosis, immunity, apoptosis, and anti-apoptotic domains that may help T. jamsteci fight against microbial pathogens, were identified. Furthermore, comparative analysis identified positive selection signals in a large number of genes including the hypoxia up-regulated protein 1, which is a chaperone that may promote adaptation of the T. jamsteci to hypoxic deepsea environments, hox2, Rx2, Pax6 and cilia-related genes BBS1, BBS2, BBS9 and RFX4. Notably, because of the critical importance of cilia and IFT in development, positive selection in cilia-related genes may play a critical role in facilitating T. jamsteci to adapt to the high-pressure deep-sea ecosystem. Results from this study thus revealed important molecular clues that may facilitate further research on the adaptation of molluscs to deep-sea hydrothermal vents.

Shared and unique patterns of autonomous human endogenous retrovirus loci transcriptomes in CD14 + monocytes from individuals with physical trauma or infection with COVID-19

Since previous studies have suggested that the RNAs of human endogenous retrovirus (HERV) might be involved in regulating innate immunity, it is important to investigate the HERV transcriptome patterns in innate immune cell types such as CD14 + monocytes. Using single cell RNA-seq datasets from resting or stimulated PBMCs mapped to 3,220 known discrete autonomous proviral HERV loci, we found individual-specific variation in HERV transcriptomes between HERV loci in CD14 + monocytes. Analysis of paired datasets from the same individual that were cultured in vitro with LPS or without (i.e. control) revealed 36 HERV loci in CD14 + monocytes that were detected only after activation. To extend our analysis to in vivo activated CD14 + monocytes, we used two scRNA-seq datasets from studies that had demonstrated activation of circulating CD14 + monocytes in patients with physical trauma or patients hospitalized with COVID-19 infections. For direct comparison between the trauma and COVID-19 datasets, we first analyzed 1.625 billion sequence reads from a composite pangenome control of 21 normal individuals. Comparison of the sequence read depth of HERV loci in the trauma or COVID-19 samples to the pangenome control revealed that 39 loci in the COVID-19 and 11 HERV loci in the trauma samples were significantly different (Mann-Whitney U test), with 9 HERV loci shared between the COVID-19 and trauma datasets. The capacity to compare HERV loci transcriptome patterns in innate immune cells, like CD14 + monocytes, across different pathological conditions will lead to greater understanding of the physiological role of HERV expression in health and disease.

Transposable elements contribute to tissue-specific gene regulation in humans

BACKGROUND

Transposable elements (TEs) contribute to approximately half of the human genome, and along with many other functions, they have been known to play a role in gene regulation in the genome. With TEs' active/repressed states varying across tissue and cell types, they have the potential to regulate gene expression in a tissue-specific manner.

OBJECTIVE AND METHODS

To provide a systematic analysis of TEs' contribution in tissue-specific gene regulation, we examined the regulatory elements and genes in association with TE-derived regulatory sequences in 14 human cell lines belonging to 10 different tissue types using the functional genomics data from the ENCODE project. Specifically, we separately analyzed regulatory regions identified by three different approaches (DNase hypersensitive sites (DHS), histone active sites (HA), and histone repressive sites (HR)).

RESULTS

These regulatory regions showed to be distinct from each other by sharing less than 2.5% among all three types and more than 95% showed to be cell line-specific. Despite a lower total TE content overall than the genome average, each regulatory sequence type showed enrichment for one or two specific TE type(s): DHS for long terminal repeats (LTRs) and DNA transposons, HA for short interspersed nucleotide elements (SINEs), and HR for LTRs. In contrast, SINE was shown to be overrepresented in all three types of regulatory sequences located in gene-neighboring regions. TE-regulated genes were mostly shown to have cell line specific pattern, and tissue-specific genes (TSGs) showed higher usage of TE regulatory sequences in the tissue of their expression. While TEs in the regulatory sequences showed to be older than their genome-wide counterparts, younger TEs were shown to be more likely used in cell line specific regulatory sequences.

CONCLUSIONS

Collectively, our study provided further evidence enforcing an important contribution of TEs to tissue-specific gene regulation in humans.

Development and evolution of Drosophila chromatin landscape in a 3D genome context

Little is known about how the epigenomic states change during development and evolution in a 3D genome context. Here we use Drosophila pseudoobscura with complex turnover of sex chromosomes as a model to address this, by collecting massive epigenomic and Hi-C data from five developmental stages and three adult tissues. We reveal that over 60% of the genes and transposable elements (TE) exhibit at least one developmental transition of chromatin state. Transitions on specific but not housekeeping enhancers are associated with specific chromatin loops and topologically associated domain borders (TABs). While evolutionarily young TEs are generally silenced, old TEs more often have been domesticated as interacting TABs or specific enhancers. But on the recently evolved X chromosome, young TEs are instead often active and recruited as TABs, due to acquisition of dosage compensation. Overall we characterize how Drosophila epigenomic landscapes change during development and in response to chromosome evolution, and highlight the important roles of TEs in genome organization and regulation.

Transposable elements as drivers of dedifferentiation: Connections between enhancers in embryonic stem cells, placenta, and cancer

Transposable elements (TEs) have emerged as important factors in establishing the cell type-specific gene regulatory networks and evolutionary novelty of embryonic and placental development. Recently, studies on the role of TEs and their dysregulation in cancers have shed light on the transcriptional, transpositional, and regulatory activity of TEs, revealing that the activation of developmental transcriptional programs by TEs may have a role in the dedifferentiation of cancer cells to the progenitor-like cell states. This essay reviews the recent evidence of the cis-regulatory TEs (henceforth crTE) in normal development and malignancy as well as the key transcription factors and regulatory pathways that are implicated in both cell states, and presents existing gaps remaining to be studied, limitations of current technologies, and therapeutic possibilities.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Endogenous retroviruses mediate transcriptional rewiring in response to oncogenic signaling in colorectal cancer

Cancer cells exhibit rewired transcriptional regulatory networks that promote tumor growth and survival. However, the mechanisms underlying the formation of these pathological networks remain poorly understood. Through a pan-cancer epigenomic analysis, we found that primate-specific endogenous retroviruses (ERVs) are a rich source of enhancers displaying cancer-specific activity. In colorectal cancer and other epithelial tumors, oncogenic MAPK/AP1 signaling drives the activation of enhancers derived from the primate-specific ERV family LTR10. Functional studies in colorectal cancer cells revealed that LTR10 elements regulate tumor-specific expression of multiple genes associated with tumorigenesis, such as ATG12 and XRCC4. Within the human population, individual LTR10 elements exhibit germline and somatic structural variation resulting from a highly mutable internal tandem repeat region, which affects AP1 binding activity. Our findings reveal that ERV-derived enhancers contribute to transcriptional dysregulation in response to oncogenic signaling and shape the evolution of cancer-specific regulatory networks.

Comprehensive Identification and Characterization of HML-9 Group in Chimpanzee Genome

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

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.

Regulation and function of transposable elements in cancer genomes

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

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.

The role of transposable elements in human evolution and methods for their functional analysis: current status and future perspectives

Transposable elements (TEs) are mobile DNA sequences that can insert themselves into various locations within the genome, causing mutations that may provide advantages or disadvantages to individuals and species. The insertion of TEs can result in genetic variation that may affect a wide range of human traits including genetic disorders. Understanding the role of TEs in human biology is crucial for both evolutionary and medical research. This review discusses the involvement of TEs in human traits and disease susceptibility, as well as methods for functional analysis of TEs.

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

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

The Association between Long-Term DDT or DDE Exposures and an Altered Sperm Epigenome-a Cross-Sectional Study of Greenlandic Inuit and South African VhaVenda Men

BACKGROUND

The organochlorine dichlorodiphenyltrichloroethane (DDT) is banned worldwide owing to its negative health effects. It is exceptionally used as an insecticide for malaria control. Exposure occurs in regions where DDT is applied, as well as in the Arctic, where its endocrine disrupting metabolite, p , p ' -dichlorodiphenyldichloroethylene (p , p ' -DDE) accumulates in marine mammals and fish. DDT and p , p ' -DDE exposures are linked to birth defects, infertility, cancer, and neurodevelopmental delays. Of particular concern is the potential of DDT use to impact the health of generations to come via the heritable sperm epigenome.

OBJECTIVES

The objective of this study was to assess the sperm epigenome in relation to p , p ' -DDE serum levels between geographically diverse populations.

METHODS

In the Limpopo Province of South Africa, we recruited 247 VhaVenda South African men and selected 50 paired blood serum and semen samples, and 47 Greenlandic Inuit blood and semen paired samples were selected from a total of 193 samples from the biobank of the INUENDO cohort, an EU Fifth Framework Programme Research and Development project. Sample selection was based on obtaining a range of p , p ' -DDE serum levels (mean = 870.734 ± 134.030  ng / mL ). We assessed the sperm epigenome in relation to serum p , p ' -DDE levels using MethylC-Capture-sequencing (MCC-seq) and chromatin immunoprecipitation followed by sequencing (ChIP-seq). We identified genomic regions with altered DNA methylation (DNAme) and differential enrichment of histone H3 lysine 4 trimethylation (H3K4me3) in sperm.

RESULTS

Differences in DNAme and H3K4me3 enrichment were identified at transposable elements and regulatory regions involved in fertility, disease, development, and neurofunction. A subset of regions with sperm DNAme and H3K4me3 that differed between exposure groups was predicted to persist in the preimplantation embryo and to be associated with embryonic gene expression.

DISCUSSION

These findings suggest that DDT and p , p ' -DDE exposure impacts the sperm epigenome in a dose-response-like manner and may negatively impact the health of future generations through epigenetic mechanisms. Confounding factors, such as other environmental exposures, genetic diversity, and selection bias, cannot be ruled out. https://doi.org/10.1289/EHP12013.

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

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

Statistical learning quantifies transposable element-mediated cis-regulation

BACKGROUND

Transposable elements (TEs) have colonized the genomes of most metazoans, and many TE-embedded sequences function as cis-regulatory elements (CREs) for genes involved in a wide range of biological processes from early embryogenesis to innate immune responses. Because of their repetitive nature, TEs have the potential to form CRE platforms enabling the coordinated and genome-wide regulation of protein-coding genes by only a handful of trans-acting transcription factors (TFs).

RESULTS

Here, we directly test this hypothesis through mathematical modeling and demonstrate that differences in expression at protein-coding genes alone are sufficient to estimate the magnitude and significance of TE-contributed cis-regulatory activities, even in contexts where TE-derived transcription fails to do so. We leverage hundreds of overexpression experiments and estimate that, overall, gene expression is influenced by TE-embedded CREs situated within approximately 500 kb of promoters. Focusing on the cis-regulatory potential of TEs within the gene regulatory network of human embryonic stem cells, we find that pluripotency-specific and evolutionarily young TE subfamilies can be reactivated by TFs involved in post-implantation embryogenesis. Finally, we show that TE subfamilies can be split into truly regulatorily active versus inactive fractions based on additional information such as matched epigenomic data, observing that TF binding may better predict TE cis-regulatory activity than differences in histone marks.

CONCLUSION

Our results suggest that TE-embedded CREs contribute to gene regulation during and beyond gastrulation. On a methodological level, we provide a statistical tool that infers TE-dependent cis-regulation from RNA-seq data alone, thus facilitating the study of TEs in the next-generation sequencing era.

Transposable Elements Are Co-opted as Oncogenic Regulatory Elements by Lineage-Specific Transcription Factors in Prostate Cancer

Transposable elements hold regulatory functions that impact cell fate determination by controlling gene expression. However, little is known about the transcriptional machinery engaged at transposable elements in pluripotent and mature versus oncogenic cell states. Through positional analysis over repetitive DNA sequences of H3K27ac chromatin immunoprecipitation sequencing data from 32 normal cell states, we report pluripotent/stem and mature cell state-specific "regulatory transposable elements." Pluripotent/stem elements are binding sites for pluripotency factors (e.g., NANOG, SOX2, OCT4). Mature cell elements are docking sites for lineage-specific transcription factors, including AR and FOXA1 in prostate epithelium. Expanding the analysis to prostate tumors, we identify a subset of regulatory transposable elements shared with pluripotent/stem cells, including Tigger3a. Using chromatin editing technology, we show how such elements promote prostate cancer growth by regulating AR transcriptional activity. Collectively, our results suggest that oncogenesis arises from lineage-specific transcription factors hijacking pluripotent/stem cell regulatory transposable elements.

SIGNIFICANCE

We show that oncogenesis relies on co-opting transposable elements from pluripotent stem cells as regulatory elements altering the recruitment of lineage-specific transcription factors. We further discover how co-option is dependent on active chromatin states with important implications for developing treatment options against drivers of oncogenesis across the repetitive DNA. This article is featured in Selected Articles from This Issue, p. 2293.

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

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

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators in various biological processes. However, due to their low expression, their systematic characterization is difficult to determine. Here, we performed transcript annotation by a newly developed computational pipeline, termed RNA-seq and small RNA-seq combined strategy (RSCS), in a wide variety of cellular contexts. Thousands of high-confidence potential novel transcripts were identified by the RSCS, and the reliability of the transcriptome was verified by analysis of transcript structure, base composition, and sequence complexity. Evidenced by the length comparison, the frequency of the core promoter and the polyadenylation signal motifs, and the locations of transcription start and end sites, the transcripts appear to be full length. Furthermore, taking advantage of our strategy, we identified a large number of endogenous retrovirus-associated lncRNAs, and a novel endogenous retrovirus-lncRNA that was functionally involved in control of Yap1 expression and essential for early embryogenesis was identified. In summary, the RSCS can generate a more complete and precise transcriptome, and our findings greatly expanded the transcriptome annotation for the mammalian community.

Transposable elements as essential elements in the control of gene expression

Interspersed repetitions called transposable elements (TEs), commonly referred to as mobile elements, make up a significant portion of the genomes of higher animals. TEs contribute in controlling the expression of genes locally and even far away at the transcriptional and post-transcriptional levels, which is one of their significant functional effects on gene function and genome evolution. There are different mechanisms through which TEs control the expression of genes. First, TEs offer cis-regulatory regions in the genome with their inherent regulatory features for their own expression, making them potential factors for controlling the expression of the host genes. Promoter and enhancer elements contain cis-regulatory sites generated from TE, which function as binding sites for a variety of trans-acting factors. Second, a significant portion of miRNAs and long non-coding RNAs (lncRNAs) have been shown to have TEs that encode for regulatory RNAs, revealing the TE origin of these RNAs. Furthermore, it was shown that TE sequences are essential for these RNAs' regulatory actions, which include binding to the target mRNA. By being a member of cis-regulatory and regulatory RNA sequences, TEs therefore play essential regulatory roles. Additionally, it has been suggested that TE-derived regulatory RNAs and cis-regulatory regions both contribute to the evolutionary novelty of gene regulation. Additionally, these regulatory systems arising from TE frequently have tissue-specific functions. The objective of this review is to discuss TE-mediated gene regulation, with a particular emphasis on the processes, contributions of various TE types, differential roles of various tissue types, based mostly on recent studies on humans.

Dynamic DNA 5-hydroxylmethylcytosine and RNA 5-methycytosine Reprogramming During Early Human Development

After implantation, complex and highly specialized molecular events render functionally distinct organ formation, whereas how the epigenome shapes organ-specific development remains to be fully elucidated. Here, nano-hmC-Seal, RNA bisulfite sequencing (RNA-BisSeq), and RNA sequencing (RNA-Seq) were performed, and the first multilayer landscapes of DNA 5-hydroxymethylcytosine (5hmC) and RNA 5-methylcytosine (m5C) epigenomes were obtained in the heart, kidney, liver, and lung of the human foetuses at 13-28 weeks with 123 samples in total. We identified 70,091 and 503 organ- and stage-specific differentially hydroxymethylated regions (DhMRs) and m5C-modified mRNAs, respectively. The key transcription factors (TFs), T-box transcription factor 20 (TBX20), paired box 8 (PAX8), krueppel-like factor 1 (KLF1), transcription factor 21 (TCF21), and CCAAT enhancer binding protein beta (CEBPB), specifically contribute to the formation of distinct organs at different stages. Additionally, 5hmC-enriched Alu elements may participate in the regulation of expression of TF-targeted genes. Our integrated studies reveal a putative essential link between DNA modification and RNA methylation, and illustrate the epigenetic maps during human foetal organogenesis, which provide a foundation for for an in-depth understanding of the epigenetic mechanisms underlying early development and birth defects.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Locus specific endogenous retroviral expression associated with Alzheimer's disease

Introduction

Human endogenous retroviruses (HERVs) are transcriptionally-active remnants of ancient retroviral infections that may play a role in Alzheimer's disease.

Methods

We combined two, publicly available RNA-Seq datasets with a third, novel dataset for a total cohort of 103 patients with Alzheimer's disease and 45 healthy controls. We use telescope to perform HERV quantification for these samples and simultaneously perform gene expression analysis.

Results

We identify differentially expressed genes and differentially expressed HERVs in Alzheimer's disease patients. Differentially expressed HERVs are scattered throughout the genome; many of them are members of the HERV-K superfamily. A number of HERVs are correlated with the expression of dysregulated genes in Alzheimer's and are physically proximal to genes which drive disease pathways.

Discussion

Dysregulated expression of ancient retroviral insertions in the human genome are present in Alzheimer's disease and show localization patterns that may explain how these elements drive pathogenic gene expression.

H4K16ac activates the transcription of transposable elements and contributes to their cis-regulatory function

Mammalian genomes harbor abundant transposable elements (TEs) and their remnants, with numerous epigenetic repression mechanisms enacted to silence TE transcription. However, TEs are upregulated during early development, neuronal lineage, and cancers, although the epigenetic factors contributing to the transcription of TEs have yet to be fully elucidated. Here, we demonstrate that the male-specific lethal (MSL)-complex-mediated histone H4 acetylation at lysine 16 (H4K16ac) is enriched at TEs in human embryonic stem cells (hESCs) and cancer cells. This in turn activates transcription of subsets of full-length long interspersed nuclear elements (LINE1s, L1s) and endogenous retrovirus (ERV) long terminal repeats (LTRs). Furthermore, we show that the H4K16ac-marked L1 and LTR subfamilies display enhancer-like functions and are enriched in genomic locations with chromatin features associated with active enhancers. Importantly, such regions often reside at boundaries of topologically associated domains and loop with genes. CRISPR-based epigenetic perturbation and genetic deletion of L1s reveal that H4K16ac-marked L1s and LTRs regulate the expression of genes in cis. Overall, TEs enriched with H4K16ac contribute to the cis-regulatory landscape at specific genomic locations by maintaining an active chromatin landscape at TEs.

Phylogenomic analyses provide insights into primate evolution

Comparative analysis of primate genomes within a phylogenetic context is essential for understanding the evolution of human genetic architecture and primate diversity. We present such a study of 50 primate species spanning 38 genera and 14 families, including 27 genomes first reported here, with many from previously less well represented groups, the New World monkeys and the Strepsirrhini. Our analyses reveal heterogeneous rates of genomic rearrangement and gene evolution across primate lineages. Thousands of genes under positive selection in different lineages play roles in the nervous, skeletal, and digestive systems and may have contributed to primate innovations and adaptations. Our study reveals that many key genomic innovations occurred in the Simiiformes ancestral node and may have had an impact on the adaptive radiation of the Simiiformes and human evolution.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

Mobile element variation contributes to population-specific genome diversification, gene regulation and disease risk

Mobile genetic elements (MEs) are heritable mutagens that recursively generate structural variants (SVs). ME variants (MEVs) are difficult to genotype and integrate in statistical genetics, obscuring their impact on genome diversification and traits. We developed a tool that accurately genotypes MEVs using short-read whole-genome sequencing (WGS) and applied it to global human populations. We find unexpected population-specific MEV differences, including an Alu insertion distribution distinguishing Japanese from other populations. Integrating MEVs with expression quantitative trait loci (eQTL) maps shows that MEV classes regulate tissue-specific gene expression by shared mechanisms, including creating or attenuating enhancers and recruiting post-transcriptional regulators, supporting class-wide interpretability. MEVs more often associate with gene expression changes than SNVs, thus plausibly impacting traits. Performing genome-wide association study (GWAS) with MEVs pinpoints potential causes of disease risk, including a LINE-1 insertion associated with keloid and fasciitis. This work implicates MEVs as drivers of human divergence and disease risk.

Genome graphs detect human polymorphisms in active epigenomic state during influenza infection

Genetic variants, including mobile element insertions (MEIs), are known to impact the epigenome. We hypothesized that genome graphs, which encapsulate genetic diversity, could reveal missing epigenomic signals. To test this, we sequenced the epigenome of monocyte-derived macrophages from 35 ancestrally diverse individuals before and after influenza infection, allowing us to investigate the role of MEIs in immunity. We characterized genetic variants and MEIs using linked reads and built a genome graph. Mapping epigenetic data revealed 2.3%-3% novel peaks for H3K4me1, H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq), and ATAC-seq. Additionally, the use of a genome graph modified some quantitative trait loci estimates and revealed 375 polymorphic MEIs in an active epigenomic state. Among these is an AluYh3 polymorphism whose chromatin state changed after infection and was associated with the expression of TRIM25, a gene that restricts influenza RNA synthesis. Our results demonstrate that graph genomes can reveal regulatory regions that would have been overlooked by other approaches.

Chiropterans Are a Hotspot for Horizontal Transfer of DNA Transposons in Mammalia

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.

HERVs and Cancer-A Comprehensive Review of the Relationship of Human Endogenous Retroviruses and Human Cancers

Genomic instability and genetic mutations can lead to exhibition of several cancer hallmarks in affected cells such as sustained proliferative signaling, evasion of growth suppression, activated invasion, deregulation of cellular energetics, and avoidance of immune destruction. Similar biological changes have been observed to be a result of pathogenic viruses and, in some cases, have been linked to virus-induced cancers. Human endogenous retroviruses (HERVs), once external pathogens, now occupy more than 8% of the human genome, representing the merge of genomic and external factors. In this review, we outline all reported effects of HERVs on cancer development and discuss the HERV targets most suitable for cancer treatments as well as ongoing clinical trials for HERV-targeting drugs. We reviewed all currently available reports of the effects of HERVs on human cancers including solid tumors, lymphomas, and leukemias. Our review highlights the central roles of HERV genes, such as gag, env, pol, np9, and rec in immune regulation, checkpoint blockade, cell differentiation, cell fusion, proliferation, metastasis, and cell transformation. In addition, we summarize the involvement of HERV long terminal repeat (LTR) regions in transcriptional regulation, creation of fusion proteins, expression of long non-coding RNAs (lncRNAs), and promotion of genome instability through recombination.

Widespread contribution of transposable elements to the rewiring of mammalian 3D genomes

Transposable elements (TEs) are major contributors of genetic material in mammalian genomes. These often include binding sites for architectural proteins, including the multifarious master protein, CTCF, which shapes the 3D genome by creating loops, domains, compartment borders, and RNA-DNA interactions. These play a role in the compact packaging of DNA and have the potential to facilitate regulatory function. In this study, we explore the widespread contribution of TEs to mammalian 3D genomes by quantifying the extent to which they give rise to loops and domain border differences across various cell types and species using several 3D genome mapping technologies. We show that specific families and subfamilies of TEs have contributed to lineage-specific 3D chromatin structures across mammalian species. In many cases, these loops may facilitate sustained interaction between distant cis-regulatory elements and target genes, and domains may segregate chromatin state to impact gene expression in a lineage-specific manner. An experimental validation of our analytical findings using CRISPR-Cas9 to delete a candidate TE resulted in disruption of species-specific 3D chromatin structure. Taken together, we comprehensively quantify and selectively validate our finding that TEs contribute to shaping 3D genome organization and may, in some cases, impact gene regulation during the course of mammalian evolution.

Transposable Elements Co-Option in Genome Evolution and Gene Regulation

The genome is no longer deemed as a fixed and inert item but rather as a moldable matter that is continuously evolving and adapting. Within this frame, Transposable Elements (TEs), ubiquitous, mobile, repetitive elements, are considered an alive portion of the genomes to date, whose functions, although long considered "dark", are now coming to light. Here we will review that, besides the detrimental effects that TE mobilization can induce, TEs have shaped genomes in their current form, promoting genome sizing, genomic rearrangements and shuffling of DNA sequences. Although TEs are mostly represented in the genomes by evolutionarily old, short, degenerated, and sedentary fossils, they have been thoroughly co-opted by the hosts as a prolific and original source of regulatory instruments for the control of gene transcription and genome organization in the nuclear space. For these reasons, the deregulation of TE expression and/or activity is implicated in the onset and progression of several diseases. It is likely that we have just revealed the outermost layers of TE functions. Further studies on this portion of the genome are required to unlock novel regulatory functions that could also be exploited for diagnostic and therapeutic approaches.

Long Noncoding RNA Lx8-SINE B2 Interacts with Eno1 to Regulate Self-Renewal and Metabolism of Embryonic Stem Cells

Long noncoding RNAs (lncRNAs) emerge as important orchestrators of biological processes in embryonic stem cells (ESCs). LncRNA Lx8-SINE B2 was recently identified as an ESC-specific lncRNA that marks pluripotency. Here, we studied the function of lncRNA Lx8-SINE B2 in ESCs. Depletion of Lx8-SINE B2 disrupted ESC proliferation, repressed the expression of pluripotency genes, activated differentiation genes, and inhibited reprogramming to induced pluripotent stem cells. The reduction of the colony formation ability of ESCs upon Lx8-SINE B2 knockdown was accompanied by the elongation of the G1 phase and the shortening of the S phase. Transcriptome analysis revealed that Lx8-SINE B2 deficiency affected multiple metabolic pathways, particularly glycolysis. Mechanistically, Lx8-SINE B2 functions as a cytoplasmic lncRNA and interacts with the glycolytic enzyme Eno1 as shown by RNA pull-down and RNA localization analysis. Lx8-SINE B2 and Eno1 interact with and regulate each other's expression, hence promoting the expression of metabolic genes and influencing glycolysis. In conclusion, we have identified lncRNA Lx8-SINE B2 as a novel regulator of ESC proliferation, cell cycle, and metabolism through working with Eno1.

Endogenous Retrovirus Elements Are Co-Expressed with IFN Stimulation Genes in the JAK-STAT Pathway

Background: Endogenous retrovirus (ERV) elements can act as proximal regulatory elements in promoting interferon (IFN) responses. Previous relevant studies have mainly focused on IFN-stimulated genes (ISGs). However, the role of ERV elements as cis-regulatory motifs in regulating genes of the JAK-STAT pathway remains poorly understood. In our study, we analyzed the changes in ERV elements and genes under both IFN stimulation and blockade of the signaling pathway. Methods: The effects of interferon on cells under normal conditions and knockout of the receptor were compared based on the THP1_IFNAR1_KO and THP1_IFNAR2_mutant cell lines. The correlation between differentially expressed ERVs (DHERVs) and differentially expressed genes (DEGs) as DEHERV-G pairs was explored with construction of gene regulatory networks related to ERV and induced by proinflammatory cytokines. Results: A total of 430 DEHERV loci and 190 DEGs were identified in 842 DEHERV-G pairs that are common to the three groups. More than 87% of DEHERV-G pairs demonstrated a consistent expression pattern. ISGs such as AIM2, IFIT1, IFIT2, IFIT3, STAT1, and IRF were activated via the JAK-STAT pathway in response to interferon stimulation. Thus, STAT1, STAT2, and IRF1 appear to play core roles in regulatory networks and are closely associated with ERVs. Conclusions: The RNA expression of ISGs and ERV elements is correlated, indicating that ERV elements are closely linked to host innate immune responses.

Metagenomic analysis of viral genes integrated in whole genome sequencing data of Thai patients with Brugada syndrome

Brugada syndrome (BS) is an autosomal dominant inheritance cardiac arrhythmia disorder associated with sudden death in young adults. Thailand has the highest prevalence of BS worldwide, and over 60% of patients with BS still have unclear disease etiology. Here, we performeda new viral metagenome analysis pipeline called VIRIN and validated it with whole genome sequencing (WGS) data of HeLa cell lines and hepatocellular carcinoma. Then the VIRIN pipelinewas applied to identify viral integration positions from unmapped WGS data of Thai males, including 100 BS patients (case) and 100 controls. Even though the sample preparation had noviral enrichment step, we can identify several virus genes from our analysis pipeline. The predominance of human endogenous retrovirus K (HERV-K) viruses was found in both cases andcontrols by blastn and blastx analysis. This study is the first report on the full-length HERV-Kassembled genomes in the Thai population. Furthermore, the HERV-K integration breakpointpositions were validated and compared between the case and control datasets. Interestingly,Brugada cases contained HERV-K integration breakpoints at promoters five times more oftenthan controls. Overall, the highlight of this study is the BS-specific HERV-K breakpoint positionsthat were found at the gene coding region "NBPF11" (n = 9), "NBPF12" (n = 8) and longnon-coding RNA (lncRNA) "PCAT14" (n = 4) region. The genes and the lncRNA have been reported to be associated with congenital heart and arterial diseases. These findings provide another aspect of the BS etiology associated with viral genome integrations within the humangenome.

Primate-specific transposable elements shape transcriptional networks during human development

The human genome contains more than 4.5 million inserts derived from transposable elements (TEs), the result of recurrent waves of invasion and internal propagation throughout evolution. For new TE copies to be inherited, they must become integrated in the genome of the germline or pre-implantation embryo, which requires that their source TE be expressed at these stages. Accordingly, many TEs harbor DNA binding sites for the pluripotency factors OCT4, NANOG, SOX2, and KLFs and are transiently expressed during embryonic genome activation. Here, we describe how many primate-restricted TEs have additional binding sites for lineage-specific transcription factors driving their expression during human gastrulation and later steps of fetal development. These TE integrants serve as lineage-specific enhancers fostering the transcription, amongst other targets, of KRAB-zinc finger proteins (KZFPs) of comparable evolutionary age, which in turn corral the activity of TE-embedded regulatory sequences in a similarly lineage-restricted fashion. Thus, TEs and their KZFP controllers play broad roles in shaping transcriptional networks during early human development.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

Background

A comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.

Results

Here we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.

Conclusions

Our data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00767-3.

Young transposable elements rewired gene regulatory networks in human and chimpanzee hippocampal intermediate progenitors

The hippocampus is associated with essential brain functions, such as learning and memory. Human hippocampal volume is significantly greater than expected compared with that of non-human apes, suggesting a recent expansion. Intermediate progenitors, which are able to undergo multiple rounds of proliferative division before a final neurogenic division, may have played a role in evolutionary hippocampal expansion. To investigate the evolution of gene regulatory networks underpinning hippocampal neurogenesis in apes, we leveraged the differentiation of human and chimpanzee induced pluripotent stem cells into TBR2 (or EOMES)-positive hippocampal intermediate progenitor cells (hpIPCs). We found that the gene networks active in hpIPCs are significantly different between humans and chimpanzees, with ∼2500 genes being differentially expressed. We demonstrate that species-specific transposon-derived enhancers contribute to these transcriptomic differences. Young transposons, predominantly endogenous retroviruses and SINE-Vntr-Alus (SVAs), were co-opted as enhancers in a species-specific manner. Human-specific SVAs provided substrates for thousands of novel TBR2-binding sites, and CRISPR-mediated repression of these SVAs attenuated the expression of ∼25% of the genes that are upregulated in human intermediate progenitors relative to the same cell population in the chimpanzee.

Functional characterization of enhancer activity during a long terminal repeat's evolution

Many transposable elements (TEs) contain transcription factor binding sites and are implicated as potential regulatory elements. However, TEs are rarely functionally tested for regulatory activity, which in turn limits our understanding of how TE regulatory activity has evolved. We systematically tested the human LTR18A subfamily for regulatory activity using massively parallel reporter assay (MPRA) and found AP-1- and CEBP-related binding motifs as drivers of enhancer activity. Functional analysis of evolutionarily reconstructed ancestral sequences revealed that LTR18A elements have generally lost regulatory activity over time through sequence changes, with the largest effects occurring owing to mutations in the AP-1 and CEBP motifs. We observed that the two motifs are conserved at higher rates than expected based on neutral evolution. Finally, we identified LTR18A elements as potential enhancers in the human genome, primarily in epithelial cells. Together, our results provide a model for the origin, evolution, and co-option of TE-derived regulatory elements.