Waves of retrotransposon expansion remodel genome organization and CTCF binding in multiple mammalian lineages

Permeable TAD boundaries and their impact on genome-associated functions

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

The characteristics of CTCF binding sequences contribute to enhancer blocking activity

While the elements encoding enhancers and promoters have been relatively well studied, the full spectrum of insulator elements which bind the CCCTC binding factor (CTCF), is relatively poorly characterized. This is partly due to the genomic context of CTCF sites greatly influencing their roles and activity. Here we have developed an experimental system to determine the ability of minimal, consistently sized, individual CTCF elements to interpose between enhancers and promoters and thereby reduce gene expression during differentiation. Importantly, each element is tested in the identical location thereby minimising the effect of genomic context. We found no correlation between the ability of CTCF elements to block enhancer-promoter activity with the degree of evolutionary conservation; their resemblance to the consensus core sequences; or the number of CTCF core motifs harboured in the element. Nevertheless, we have shown that the strongest enhancer-promoter blockers include a previously described bound element lying upstream of the CTCF core motif. In addition, we found other uncharacterised DNaseI footprints located close to the core motif that may affect function. We have developed an assay of CTCF sequences which will enable researchers to sub-classify individual CTCF elements in a uniform and unbiased way.

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.

Loops are geometric catalysts for DNA integration

The insertion of DNA elements within genomes underpins both genetic diversity and disease when unregulated. Most of DNA insertions are not random and the physical mechanisms underlying the integration site selection are poorly understood. Here, we perform Molecular Dynamics simulations to study the insertion of DNA elements, such as viral DNA or transposons, into naked DNA or chromatin substrates. More specifically, we explore the role of loops within the polymeric substrate and discover that they act as 'geometric catalysts' for DNA integration by reducing the energy barrier for substrate deformation. Additionally, we discover that the 1D pattern and 3D conformation of loops have a marked effect on the distribution of integration sites. Finally, we show that loops may compete with nucleosomes to attract DNA integrations. These results may be tested in vitro and they may help to understand patterns of DNA insertions with implications in genome evolution and engineering.

Accurate allocation of multimapped reads enables regulatory element analysis at repeats

Transposable elements (TEs) and other repetitive regions have been shown to contain gene regulatory elements, including transcription factor binding sites. However, regulatory elements harbored by repeats have proven difficult to characterize using short-read sequencing assays such as ChIP-seq or ATAC-seq. Most regulatory genomics analysis pipelines discard "multimapped" reads that align equally well to multiple genomic locations. Because multimapped reads arise predominantly from repeats, current analysis pipelines fail to detect a substantial portion of regulatory events that occur in repetitive regions. To address this shortcoming, we developed Allo, a new approach to allocate multimapped reads in an efficient, accurate, and user-friendly manner. Allo combines probabilistic mapping of multimapped reads with a convolutional neural network that recognizes the read distribution features of potential peaks, offering enhanced accuracy in multimapping read assignment. Allo also provides read-level output in the form of a corrected alignment file, making it compatible with existing regulatory genomics analysis pipelines and downstream peak-finders. In a demonstration application on CTCF ChIP-seq data, we show that Allo results in the discovery of thousands of new CTCF peaks. Many of these peaks contain the expected cognate motif and/or serve as TAD boundaries. We additionally apply Allo to a diverse collection of ENCODE ChIP-seq data sets, resulting in multiple previously unidentified interactions between transcription factors and repetitive element families. Finally, we show that Allo may be particularly beneficial in identifying ChIP-seq peaks at centromeres, near segmentally duplicated genes, and in younger TEs, enabling new regulatory analyses in these regions.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

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.

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.

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.

Transposable elements as scaffold/matrix attachment regions: shaping organization and functions in genomes

The hierarchical structure of eukaryotic genomes has regulatory layers, one of them being epigenetic "indexing" of the genome that leads to cell-type-specific patterns of gene expression. By establishing loops and defining chromatin domains, cells can achieve coordinated control over multi-locus segments of the genome. This is thought to be achieved using scaffold/matrix attachment regions (S/MARs) that establish structural and functional loops and topologically associating domains (TADs) that define a self-interacting region of the genome. Large-scale genome-wide mapping of S/MARs has begun to uncover these aspects of genome organization. A recent genome-wide study showed the association of transposable elements (TEs) with a significant fraction of S/MARs, suggesting that the multitude of TE-derived repeats constitute a class of anchorage sites of chromatin loops to nuclear architecture. In this study, we provide an insight that TE-driven dispersal of S/MARs has the potential to restructure the chromosomes by creating novel loops and domains. The combination of TEs and S/MARs, as elements that can hop through the genome along with regulatory capabilities, may provide an active mechanism of genome evolution leading to the emergence of novel features in biological systems. The significance is that a genome-wide study mapping developmental S/MARs reveals an intriguing link between these elements and TEs. This article highlights the potential of the TE-S/MAR combination to drive evolution by restructuring and shaping the genome.

DNA satellite and chromatin organization at mouse centromeres and pericentromeres

BACKGROUND

Centromeres are essential for faithful chromosome segregation during mitosis and meiosis. However, the organization of satellite DNA and chromatin at mouse centromeres and pericentromeres is poorly understood due to the challenges of assembling repetitive genomic regions.

RESULTS

Using recently available PacBio long-read sequencing data from the C57BL/6 strain, we find that contrary to the previous reports of their homogeneous nature, both centromeric minor satellites and pericentromeric major satellites exhibit a high degree of variation in sequence and organization within and between arrays. While most arrays are continuous, a significant fraction is interspersed with non-satellite sequences, including transposable elements. Using chromatin immunoprecipitation sequencing (ChIP-seq), we find that the occupancy of CENP-A and H3K9me3 chromatin at centromeric and pericentric regions, respectively, is associated with increased sequence enrichment and homogeneity at these regions. The transposable elements at centromeric regions are not part of functional centromeres as they lack significant CENP-A enrichment. Furthermore, both CENP-A and H3K9me3 nucleosomes occupy minor and major satellites spanning centromeric-pericentric junctions and a low yet significant amount of CENP-A spreads locally at centromere junctions on both pericentric and telocentric sides. Finally, while H3K9me3 nucleosomes display a well-phased organization on major satellite arrays, CENP-A nucleosomes on minor satellite arrays are poorly phased. Interestingly, the homogeneous class of major satellites also phase CENP-A and H3K27me3 nucleosomes, indicating that the nucleosome phasing is an inherent property of homogeneous major satellites.

CONCLUSIONS

Our findings reveal that mouse centromeres and pericentromeres display a high diversity in satellite sequence, organization, and chromatin structure.

Mouse retrotransposons: sequence structure, evolutionary age, genomic distribution and function

Retrotransposons are transposable elements that are transposed via transcription and reverse transcription. Their copies have accumulated in the genome of mammals, occupying approximately 40% of mammalian genomic mass. These copies are often involved in numerous phenomena, such as chromatin spatial organization, gene expression, development and disease, and have been recognized as a driving force in evolution. Different organisms have gained specific retrotransposon subfamilies and retrotransposed copies, such as hundreds of Mus-specific subfamilies with diverse sequences and genomic locations. Despite this complexity, basic information is still necessary for present-day genomic and epigenomic studies. Herein, we describe the characteristics of each subfamily of Mus-specific retrotransposons in terms of sequence structure, phylogenetic relationships, evolutionary age, and preference for A or B compartments of chromatin.

Towards targeting transposable elements for cancer therapy

Transposable elements (TEs) represent almost half of the human genome. Historically deemed 'junk DNA', recent technological advancements have stimulated a wave of research into the functional impact of TEs on gene-regulatory networks in evolution and development, as well as in diseases including cancer. The genetic and epigenetic evolution of cancer involves the exploitation of TEs, whereby TEs contribute directly to cancer-specific gene activities. This Review provides a perspective on the role of TEs in cancer as being a 'double-edged sword', both promoting cancer evolution and representing a vulnerability that could be exploited in cancer therapy. We discuss how TEs affect transcriptome regulation and other cellular processes in cancer. We highlight the potential of TEs as therapeutic targets for cancer. We also summarize technical hurdles in the characterization of TEs with genomic assays. Last, we outline open questions and exciting future research avenues.

ZNF143 deletion alters enhancer/promoter looping and CTCF/cohesin geometry

The transcription factor ZNF143 contains a central domain of seven zinc fingers in a tandem array and is involved in 3D genome construction. However, the mechanism by which ZNF143 functions in chromatin looping remains unclear. Here, we show that ZNF143 directionally recognizes a diverse range of genomic sites directly within enhancers and promoters and is required for chromatin looping between these sites. In addition, ZNF143 is located between CTCF and cohesin at numerous CTCF sites, and ZNF143 removal narrows the space between CTCF and cohesin. Moreover, genetic deletion of ZNF143, in conjunction with acute CTCF degradation, reveals that ZNF143 and CTCF collaborate to regulate higher-order topological chromatin organization. Finally, CTCF depletion enlarges direct ZNF143 chromatin looping. Thus, ZNF143 is recruited by CTCF to the CTCF sites to regulate CTCF/cohesin configuration and TAD (topologically associating domain) formation, whereas directional recognition of genomic DNA motifs directly by ZNF143 itself regulates promoter activity via chromatin looping.

A SINE-VNTR-Alu at the LRIG2 locus is associated with proximal and distal gene expression in CRISPR and population models

SINE-VNTR-Alu (SVA) retrotransposons represent mobile regulatory elements that have the potential to influence the surrounding genome when they insert into a locus. Evolutionarily recent mobilisation has resulted in loci in the human genome where a given retrotransposon might be observed to be present or absent, termed a retrotransposon insertion polymorphism (RIP). We previously observed that an SVA RIP ~ 2 kb upstream of LRIG2 on chromosome 1, the 'LRIG2 SVA', was associated with differences in local gene expression and methylation, and that the two were correlated. Here, we have used CRISPR-mediated deletion of the LRIG2 SVA in a cell line model to validate that presence of the retrotransposon is directly affecting local expression and provide evidence that is suggestive of a modest role for the SVA in modulating nearby methylation. Additionally, in leveraging an available Hi-C dataset we observed that the LRIG2 SVA was also involved in long-range chromatin interactions with a cluster of genes ~ 300 kb away, and that expression of these genes was to varying degrees associated with dosage of the SVA in both CRISPR cell line and population models. Altogether, these data support a regulatory role for SVAs in the modulation of gene expression, with the latter potentially involving chromatin looping, consistent with the model that RIPs may contribute to interpersonal differences in transcriptional networks.

Nutrient inputs and social metabolic control of T cell fate

Cells in multicellular organisms experience diverse neighbors, signals, and evolving physical environments that drive functional and metabolic demands. To maintain proper development and homeostasis while avoiding inappropriate cell proliferation or death, individual cells interact with their neighbors via "social" cues to share and partition available nutrients. Metabolic signals also contribute to cell fate by providing biochemical links between cell-extrinsic signals and available resources. In addition to metabolic checkpoints that sense nutrients and directly supply molecular intermediates for biosynthetic pathways, many metabolites directly signal or provide the basis for post-translational modifications of target proteins and chromatin. In this review, we survey the landscape of T cell nutrient sensing and metabolic signaling that supports proper immunity while avoiding immunodeficiency or autoimmunity. The integration of cell-extrinsic microenvironmental cues with cell-intrinsic metabolic signaling provides a social metabolic control model to integrate cell signaling, metabolism, and fate.

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.

Position-specific evolution in transcription factor binding sites, and a fast likelihood calculation for the F81 model

Transcription factor binding sites (TFBS), like other DNA sequence, evolve via mutation and selection relating to their function. Models of nucleotide evolution describe DNA evolution via single-nucleotide mutation. A stationary vector of such a model is the long-term distribution of nucleotides, unchanging under the model. Neutrally evolving sites may have uniform stationary vectors, but one expects that sites within a TFBS instead have stationary vectors reflective of the fitness of various nucleotides at those positions. We introduce 'position-specific stationary vectors' (PSSVs), the collection of stationary vectors at each site in a TFBS locus, analogous to the position weight matrix (PWM) commonly used to describe TFBS. We infer PSSVs for human TFs using two evolutionary models (Felsenstein 1981 and Hasegawa-Kishino-Yano 1985). We find that PSSVs reflect the nucleotide distribution from PWMs, but with reduced specificity. We infer ancestral nucleotide distributions at individual positions and calculate 'conditional PSSVs' conditioned on specific choices of majority ancestral nucleotide. We find that certain ancestral nucleotides exert a strong evolutionary pressure on neighbouring sequence while others have a negligible effect. Finally, we present a fast likelihood calculation for the F81 model on moderate-sized trees that makes this approach feasible for large-scale studies along these lines.

Cohesin and CTCF do not assemble TADs in Xenopus sperm and male pronuclei

Paternal genomes are compacted during spermiogenesis and decompacted following fertilization. These processes are fundamental for inheritance but incompletely understood. We analyzed these processes in the frog Xenopus laevis, whose sperm can be assembled into functional pronuclei in egg extracts in vitro. In such extracts, cohesin extrudes DNA into loops, but in vivo cohesin only assembles topologically associating domains (TADs) at the mid-blastula transition (MBT). Why cohesin assembles TADs only at this stage is unknown. We first analyzed genome architecture in frog sperm and compared it to human and mouse. Our results indicate that sperm genome organization is conserved between frogs and humans and occurs without formation of TADs. TADs can be detected in mouse sperm samples, as reported, but these structures might originate from somatic chromatin contaminations. We therefore discuss the possibility that the absence of TADs might be a general feature of vertebrate sperm. To analyze sperm genome remodeling upon fertilization, we reconstituted male pronuclei in Xenopus egg extracts. In pronuclei, chromatin compartmentalization increases, but cohesin does not accumulate at CTCF sites and assemble TADs. However, if pronuclei are formed in the presence of exogenous CTCF, CTCF binds to its consensus sites, and cohesin accumulates at these and forms short-range chromatin loops, which are preferentially anchored at CTCF's N terminus. These results indicate that TADs are only assembled at MBT because before this stage CTCF sites are not occupied and cohesin only forms short-range chromatin loops.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Transposon wave remodeled the epigenomic landscape in the rapid evolution of X-Chromosome dosage compensation

Sex chromosome dosage compensation is a model to understand the coordinated evolution of transcription; however, the advanced age of the sex chromosomes in model systems makes it difficult to study how the complex regulatory mechanisms underlying chromosome-wide dosage compensation can evolve. The sex chromosomes of Poecilia picta have undergone recent and rapid divergence, resulting in widespread gene loss on the male Y, coupled with complete X Chromosome dosage compensation, the first case reported in a fish. The recent de novo origin of dosage compensation presents a unique opportunity to understand the genetic and evolutionary basis of coordinated chromosomal gene regulation. By combining a new chromosome-level assembly of P. picta with whole-genome bisulfite sequencing and RNA-seq data, we determine that the YY1 transcription factor (YY1) DNA binding motif is associated with male-specific hypomethylated regions on the X, but not the autosomes. These YY1 motifs are the result of a recent and rapid repetitive element expansion on the P. picta X Chromosome, which is absent in closely related species that lack dosage compensation. Taken together, our results present compelling support that a disruptive wave of repetitive element insertions carrying YY1 motifs resulted in the remodeling of the X Chromosome epigenomic landscape and the rapid de novo origin of a dosage compensation system.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

NR5A2 connects zygotic genome activation to the first lineage segregation in totipotent embryos

Zygotic genome activation (ZGA) marks the beginning of the embryonic program for a totipotent embryo, which gives rise to the inner cell mass (ICM) where pluripotent epiblast arises, and extraembryonic trophectoderm. However, how ZGA is connected to the first lineage segregation in mammalian embryos remains elusive. Here, we investigated the role of nuclear receptor (NR) transcription factors (TFs), whose motifs are highly enriched and accessible from the 2-cell (2C) to 8-cell (8C) stages in mouse embryos. We found that NR5A2, an NR TF strongly induced upon ZGA, was required for this connection. Upon Nr5a2 knockdown or knockout, embryos developed beyond 2C normally with the zygotic genome largely activated. However, 4-8C-specific gene activation was substantially impaired and Nr5a2-deficient embryos subsequently arrested at the morula stage. Genome-wide chromatin binding analysis showed that NR5A2-bound cis-regulatory elements in both 2C and 8C embryos are strongly enriched for B1 elements where its binding motif is embedded. NR5A2 was not required for the global opening of its binding sites in 2C embryos but was essential to the opening of its 8C-specific binding sites. These 8C-specific, but not 2C-specific, binding sites are enriched near genes involved in blastocyst and stem cell regulation, and are often bound by master pluripotency TFs in blastocysts and embryonic stem cells (ESCs). Importantly, NR5A2 regulated key pluripotency genes Nanog and Pou5f1/Oct4, and primitive endoderm regulatory genes including Gata6 among many early ICM genes, as well as key trophectoderm regulatory genes including Tead4 and Gata3 at the 8C stage. By contrast, master pluripotency TFs NANOG, SOX2, and OCT4 targeted both early and late ICM genes in mouse ESCs. Taken together, these data identify NR5A2 as a key regulator in totipotent embryos that bridges ZGA to the first lineage segregation during mouse early development.

Out of the Silence: Insights into How Genes Escape X-Chromosome Inactivation

The silencing of all but one X chromosome in mammalian cells is a remarkable epigenetic process leading to near dosage equivalence in X-linked gene products between the sexes. However, equally remarkable is the ability of a subset of genes to continue to be expressed from the otherwise inactive X chromosome-in some cases constitutively, while other genes are variable between individuals, tissues or cells. In this review we discuss the advantages and disadvantages of the approaches that have been used to identify escapees. The identity of escapees provides important clues to mechanisms underlying escape from XCI, an arena of study now moving from correlation to functional studies. As most escapees show greater expression in females, the not-so-inactive X chromosome is a substantial contributor to sex differences in humans, and we highlight some examples of such impact.

The Molecular Impacts of Retrotransposons in Development and Diseases

Retrotransposons are invasive genetic elements that constitute substantial portions of mammalian genomes. They have the potential to influence nearby gene expression through their cis-regulatory sequences, reverse transcription machinery, and the ability to mold higher-order chromatin structures. Due to their multifaceted functions, it is crucial for host fitness to maintain strict regulation of these parasitic sequences to ensure proper growth and development. This review explores how subsets of retrotransposons have undergone evolutionary exaptation to enhance the complexity of mammalian genomes. It also highlights the significance of regulating these elements, drawing on recent studies conducted in human and murine systems.

In silico discovery of repetitive elements as key sequence determinants of 3D genome folding

Natural and experimental genetic variants can modify DNA loops and insulating boundaries to tune transcription, but it is unknown how sequence perturbations affect chromatin organization genome wide. We developed a deep-learning strategy to quantify the effect of any insertion, deletion, or substitution on chromatin contacts and systematically scored millions of synthetic variants. While most genetic manipulations have little impact, regions with CTCF motifs and active transcription are highly sensitive, as expected. Our unbiased screen and subsequent targeted experiments also point to noncoding RNA genes and several families of repetitive elements as CTCF-motif-free DNA sequences with particularly large effects on nearby chromatin interactions, sometimes exceeding the effects of CTCF sites and explaining interactions that lack CTCF. We anticipate that our disruption tracks may be of broad interest and utility as a measure of 3D genome sensitivity, and our computational strategies may serve as a template for biological inquiry with deep learning.

Transposable elements in mammalian chromatin organization

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

CRISPR deletion of a SINE-VNTR-Alu (SVA_67) retrotransposon demonstrates its ability to differentially modulate gene expression at the MAPT locus

Background

SINE-VNTR-Alu (SVA) retrotransposons are hominid-specific elements which have been shown to play important roles in processes such as chromatin structure remodelling and regulation of gene expression demonstrating that these repetitive elements exert regulatory functions. We have previously shown that the presence or absence of a specific SVA element, termed SVA_67, was associated with differential expression of several genes at the MAPT locus, a locus associated with Parkinson's Disease (PD) and frontotemporal dementia. However, we were not able to demonstrate that causation of differential gene expression was directed by the SVA due to lack of functional validation.

Methods

We performed CRISPR to delete SVA_67 in the HEK293 cell line. Quantification of target gene expression was performed using qPCR to assess the effects on expression in response to the deletion of SVA_67. Differences between CRISPR edit and control cell lines were analysed using two-tailed t-test with a minimum 95% confidence interval to determine statistical significance.

Results

In this study, we provide data highlighting the SVA-specific effect on differential gene expression. We demonstrate that the hemizygous deletion of the endogenous SVA_67 in CRISPR edited cell lines was associated with differential expression of several genes at the MAPT locus associated with neurodegenerative diseases including KANSL1, MAPT and LRRC37A.

Discussion

This data is consistent with our previous bioinformatic work of differential gene expression analysis using transcriptomic data from the Parkinson's Progression Markers Initiative (PPMI) cohort. As SVAs have regulatory influences on gene expression, and insertion polymorphisms contribute to interpersonal differences in expression patterns, these results highlight the potential contribution of these elements to complex diseases with potentially many genetic components, such as PD.

Regulation of the three-dimensional chromatin organization by transposable elements in pig spleen

Like other mammalian species, the pig genome is abundant with transposable elements (TEs). The importance of TEs for three-dimensional (3D) chromatin organization has been observed in species like human and mouse, yet current understanding about pig TEs is absent. Here, we investigated the contribution of TEs for the 3D chromatin organization in three pig tissues, focusing on spleen which is crucial for both adaptive and innate immunity. We identified dozens of TE families overrepresented with CTCF binding sites, including LTR22_SS, LTR15_SS and LTR16_SSc which are pig-specific families of endogenous retroviruses (ERVs). Interestingly, LTR22_SS elements harbor a CTCF motif and create hundreds of CTCF binding sites that are associated with adaptive immunity. We further applied Hi-C to profile the 3D chromatin structure in spleen and found that TE-derived CTCF binding sites correlate with chromatin insulation and frequently overlap TAD borders and loop anchors. Notably, one LTR22_SS-derived CTCF binding site demarcate a TAD boundary upstream of XCL1, which is a spleen-enriched chemokine gene important for lymphocyte trafficking and inflammation. Overall, this study represents a first step toward understanding the function of TEs on 3D chromatin organization regulation in pigs and expands our understanding about the functional importance of TEs in mammals.

DNA satellite and chromatin organization at house mouse centromeres and pericentromeres

Centromeres are essential for faithful chromosome segregation during mitosis and meiosis. However, the organization of satellite DNA and chromatin at mouse centromeres and pericentromeres is poorly understood due to the challenges of sequencing and assembling repetitive genomic regions. Using recently available PacBio long-read sequencing data from the C57BL/6 strain and chromatin profiling, we found that contrary to the previous reports of their highly homogeneous nature, centromeric and pericentromeric satellites display varied sequences and organization. We find that both centromeric minor satellites and pericentromeric major satellites exhibited sequence variations within and between arrays. While most arrays are continuous, a significant fraction is interspersed with non-satellite sequences, including transposable elements. Additionally, we investigated CENP-A and H3K9me3 chromatin organization at centromeres and pericentromeres using Chromatin immunoprecipitation sequencing (ChIP-seq). We found that the occupancy of CENP-A and H3K9me3 chromatin at centromeric and pericentric regions, respectively, is associated with increased sequence abundance and homogeneity at these regions. Furthermore, the transposable elements at centromeric regions are not part of functional centromeres as they lack CENP-A enrichment. Finally, we found that while H3K9me3 nucleosomes display a well-phased organization on major satellite arrays, CENP-A nucleosomes on minor satellite arrays lack phased organization. Interestingly, the homogeneous class of major satellites phase CENP-A and H3K27me3 nucleosomes as well, indicating that the nucleosome phasing is an inherent property of homogeneous major satellites. Overall, our findings reveal that house mouse centromeres and pericentromeres, which were previously thought to be highly homogenous, display significant diversity in satellite sequence, organization, and chromatin structure.

Transposable Elements: Epigenetic Silencing Mechanisms or Modulating Tools for Vertebrate Adaptations? Two Sides of the Same Coin

Transposable elements constitute one of the main components of eukaryotic genomes. In vertebrates, they differ in content, typology, and family diversity and played a crucial role in the evolution of this taxon. However, due to their transposition ability, TEs can be responsible for genome instability, and thus silencing mechanisms were evolved to allow the coexistence between TEs and eukaryotic host-coding genes. Several papers are highlighting in TEs the presence of regulatory elements involved in regulating nearby genes in a tissue-specific fashion. This suggests that TEs are not sequences merely to silence; rather, they can be domesticated for the regulation of host-coding gene expression, permitting species adaptation and resilience as well as ensuring human health. This review presents the main silencing mechanisms acting in vertebrates and the importance of exploiting these mechanisms for TE control to rewire gene expression networks, challenging the general view of TEs as threatening elements.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

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.

Refining the genomic determinants underlying escape from X-chromosome inactivation

X-chromosome inactivation (XCI) epigenetically silences one X chromosome in every cell in female mammals. Although the majority of X-linked genes are silenced, in humans 20% or more are able to escape inactivation and continue to be expressed. Such escape genes are important contributors to sex differences in gene expression, and may impact the phenotypes of X aneuploidies; yet the mechanisms regulating escape from XCI are not understood. We have performed an enrichment analysis of transcription factor binding on the X chromosome, providing new evidence for enriched factors at the transcription start sites of escape genes. The top escape-enriched transcription factors were detected at the RPS4X promoter, a well-described human escape gene previously demonstrated to escape from XCI in a transgenic mouse model. Using a cell line model system that allows for targeted integration and inactivation of transgenes on the mouse X chromosome, we further assessed combinations of RPS4X promoter and genic elements for their ability to drive escape from XCI. We identified a small transgenic construct of only 6 kb capable of robust escape from XCI, establishing that gene-proximal elements are sufficient to permit escape, and highlighting the additive effect of multiple elements that work together in a context-specific fashion.

CpG island turnover events predict evolutionary changes in enhancer activity

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

CTCF and Its Multi-Partner Network for Chromatin Regulation

Architectural proteins are essential epigenetic regulators that play a critical role in organizing chromatin and controlling gene expression. CTCF (CCCTC-binding factor) is a key architectural protein responsible for maintaining the intricate 3D structure of chromatin. Because of its multivalent properties and plasticity to bind various sequences, CTCF is similar to a Swiss knife for genome organization. Despite the importance of this protein, its mechanisms of action are not fully elucidated. It has been hypothesized that its versatility is achieved through interaction with multiple partners, forming a complex network that regulates chromatin folding within the nucleus. In this review, we delve into CTCF's interactions with other molecules involved in epigenetic processes, particularly histone and DNA demethylases, as well as several long non-coding RNAs (lncRNAs) that are able to recruit CTCF. Our review highlights the importance of CTCF partners to shed light on chromatin regulation and pave the way for future exploration of the mechanisms that enable the finely-tuned role of CTCF as a master regulator of chromatin.

Mouse B2 SINE elements function as IFN-inducible enhancers

Regulatory networks underlying innate immunity continually face selective pressures to adapt to new and evolving pathogens. Transposable elements (TEs) can affect immune gene expression as a source of inducible regulatory elements, but the significance of these elements in facilitating evolutionary diversification of innate immunity remains largely unexplored. Here, we investigated the mouse epigenomic response to type II interferon (IFN) signaling and discovered that elements from a subfamily of B2 SINE (B2_Mm2) contain STAT1 binding sites and function as IFN-inducible enhancers. CRISPR deletion experiments in mouse cells demonstrated that a B2_Mm2 element has been co-opted as an enhancer driving IFN-inducible expression of Dicer1. The rodent-specific B2 SINE family is highly abundant in the mouse genome and elements have been previously characterized to exhibit promoter, insulator, and non-coding RNA activity. Our work establishes a new role for B2 elements as inducible enhancer elements that influence mouse immunity, and exemplifies how lineage-specific TEs can facilitate evolutionary turnover and divergence of innate immune regulatory networks.

Comparing genomic and epigenomic features across species using the WashU Comparative Epigenome Browser

Genome browsers have become an intuitive and critical tool to visualize and analyze genomic features and data. Conventional genome browsers display data/annotations on a single reference genome/assembly; there are also genomic alignment viewer/browsers that help users visualize alignment, mismatch, and rearrangement between syntenic regions. However, there is a growing need for a comparative epigenome browser that can display genomic and epigenomic data sets across different species and enable users to compare them between syntenic regions. Here, we present the WashU Comparative Epigenome Browser. It allows users to load functional genomic data sets/annotations mapped to different genomes and display them over syntenic regions simultaneously. The browser also displays genetic differences between the genomes from single-nucleotide variants (SNVs) to structural variants (SVs) to visualize the association between epigenomic differences and genetic differences. Instead of anchoring all data sets to the reference genome coordinates, it creates independent coordinates of different genome assemblies to faithfully present features and data mapped to different genomes. It uses a simple, intuitive genome-align track to illustrate the syntenic relationship between different species. It extends the widely used WashU Epigenome Browser infrastructure and can be expanded to support multiple species. This new browser function will greatly facilitate comparative genomic/epigenomic research, as well as support the recent growing needs to directly compare and benchmark the T2T CHM13 assembly and other human genome assemblies.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

Topologically associating domain boundaries are required for normal genome function

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes^1-3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

Comparative single cell epigenomic analysis of gene regulatory programs in the rodent and primate neocortex

Sequence divergence of cis- regulatory elements drives species-specific traits, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains to be elucidated. We investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset, and mouse with single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome, and chromosomal conformation profiles from a total of over 180,000 cells. For each modality, we determined species-specific, divergent, and conserved gene expression and epigenetic features at multiple levels. We find that cell type-specific gene expression evolves more rapidly than broadly expressed genes and that epigenetic status at distal candidate cis -regulatory elements (cCREs) evolves faster than promoters. Strikingly, transposable elements (TEs) contribute to nearly 80% of the human-specific cCREs in cortical cells. Through machine learning, we develop sequence-based predictors of cCREs in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Lastly, we show that epigenetic conservation combined with sequence similarity helps uncover functional cis -regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Conserved microRNAs and Flipons Shape Gene Expression during Development by Altering Promoter Conformations

The classical view of gene regulation draws from prokaryotic models, where responses to environmental changes involve operons regulated by sequence-specific protein interactions with DNA, although it is now known that operons are also modulated by small RNAs. In eukaryotes, pathways based on microRNAs (miR) regulate the readout of genomic information from transcripts, while alternative nucleic acid structures encoded by flipons influence the readout of genetic programs from DNA. Here, we provide evidence that miR- and flipon-based mechanisms are deeply connected. We analyze the connection between flipon conformation and the 211 highly conserved human miR that are shared with other placental and other bilateral species. The direct interaction between conserved miR (c-miR) and flipons is supported by sequence alignments and the engagement of argonaute proteins by experimentally validated flipons as well as their enrichment in promoters of coding transcripts important in multicellular development, cell surface glycosylation and glutamatergic synapse specification with significant enrichments at false discovery rates as low as 10^-116. We also identify a second subset of c-miR that targets flipons essential for retrotransposon replication, exploiting that vulnerability to limit their spread. We propose that miR can act in a combinatorial manner to regulate the readout of genetic information by specifying when and where flipons form non-B DNA (NoB) conformations, providing the interactions of the conserved hsa-miR-324-3p with RELA and the conserved hsa-miR-744 with ARHGAP5 genes as examples.

Murine Gammaherpesvirus 68 ORF45 Stimulates B2 Retrotransposon and Pre-tRNA Activation in a Manner Dependent on Mitogen-Activated Protein Kinase (MAPK) Signaling

RNA polymerase III (RNAPIII) transcribes a variety of noncoding RNAs, including tRNA (tRNA) and the B2 family of short interspersed nuclear elements (SINEs). B2 SINEs are noncoding retrotransposons that possess tRNA-like promoters and are normally silenced in healthy somatic tissue. Infection with the murine gammaherpesvirus MHV68 induces transcription of both SINEs and tRNAs, in part through the activity of the viral protein kinase ORF36. Here, we identify the conserved MHV68 tegument protein ORF45 as an additional activator of these RNAPIII loci. MHV68 ORF45 and ORF36 form a complex, resulting in an additive induction RNAPIII and increased ORF45 expression. ORF45-induced RNAPIII transcription is dependent on its activation of the extracellular signal-regulated kinase (ERK) mitogen-activated protein kinase (MAPK) signaling pathway, which in turn increases the abundance of the RNAPIII transcription factor Brf1. Other viral and nonviral activators of MAPK/ERK signaling also increase the levels of Brf1 protein, B2 SINE RNA, and tRNA, suggesting that this is a common strategy to increase RNAPIII activity. IMPORTANCE Gammaherpesviral infection alters the gene expression landscape of a host cell, including through the induction of noncoding RNAs transcribed by RNA polymerase III (RNAPIII). Among these are a class of repetitive genes known as retrotransposons, which are normally silenced elements and can copy and spread throughout the genome, and transfer RNAs (tRNAs), which are fundamental components of protein translation machinery. How these loci are activated during infection is not well understood. Here, we identify ORF45 from the model murine gammaherpesvirus MHV68 as a novel activator of RNAPIII transcription. To do so, it engages the MAPK/ERK signaling pathway, which is a central regulator of cellular response to environmental stimuli. Activation of this pathway leads to the upregulation of a key factor required for RNAPIII activity, Brf1. These findings expand our understanding of the regulation and dysregulation of RNAPIII transcription and highlight how viral cooption of key signaling pathways can impact host gene expression.