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

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

Natural killer cells (NK) are the \"professional killer\" of tumors and play a crucial role in anti-tumor immunotherapy. NK cell desensitization is a key mechanism of tumor immune escape. Dysregulated NKG2D-NKG2DL signaling is a primary driver of this desensitization process. However, the factors that regulate NK cell desensitization remain largely uncharacterized. Here, we present the first report that circular RNA circARAP2 (hsa_circ_0069396) is involved in the soluble MICA (sMICA)-induced NKG2D endocytosis in the NK cell desensitization model. CircARAP2 was upregulated during NK cell desensitization and the loss of circARAP2 alleviated NKG2D endocytosis and NK cell desensitization. Using Chromatin isolation by RNA purification (ChIRP) and RNA pull-down approaches, we identified that RAB5A, a molecular marker of early endosomes, was its downstream target. Notably, transcription factor CTCF was an intermediate functional partner of circARAP2. Mechanistically, we discovered that circARAP2 interacted with CTCF and inhibited the recruitment of CTCF-Polycomb Repressive Complex 2 (PRC2) to the promoter region of RAB5A, thereby erasing histone H3K27 and H3K9 methylation suppression to enhance RAB5A transcription. These data demonstrate that inhibition of circARAP2 effectively alleviates sMICA-induced NKG2D endocytosis and NK cell desensitization, providing a novel target for therapeutic intervention in tumor immune evasion.

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.

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

BACKGROUND: Transposable elements play a critical role in maintaining genome architecture during neurodevelopment. Short Interspersed Nuclear Elements (SINEs), a major subtype of transposable elements, are known to harbor binding sites for the CCCTC-binding factor (CTCF) and pivotal in orchestrating chromatin organization. However, the regulatory mechanisms controlling the activity of SINEs in the developing brain remains elusive.
RESULTS: In our study, we conduct a comprehensive genome-wide epigenetic analysis in mouse neural precursor cells using ATAC-seq, ChIP-seq, whole genome bisulfite sequencing, in situ Hi-C, and RNA-seq. Our findings reveal that the SET domain bifurcated histone lysine methyltransferase 1 (SETDB1)-mediated H3K9me3, in conjunction with DNA methylation, restricts chromatin accessibility on a selective subset of SINEs in neural precursor cells. Mechanistically, loss of Setdb1 increases CTCF access to these SINE elements and contributes to chromatin loop reorganization. Moreover, de novo loop formation contributes to differential gene expression, including the dysregulation of genes enriched in mitotic pathways. This leads to the disruptions of cell proliferation in the embryonic brain after genetic ablation of Setdb1 both in vitro and in vivo.
CONCLUSIONS: In summary, our study sheds light on the epigenetic regulation of SINEs in mouse neural precursor cells, suggesting their role in maintaining chromatin organization and cell proliferation during neurodevelopment.

The three-dimensional genome structure organized by CTCF is required for development. Clinically identified mutations in CTCF have been linked to adverse developmental outcomes. Nevertheless, the underlying mechanism remains elusive. In this investigation, we explore the regulatory roles of a clinically relevant R567W point mutation, located within the 11th zinc finger of CTCF, by introducing this mutation into both murine models and human embryonic stem cell-derived cortical organoid models. Mice with homozygous CTCFR567W mutation exhibit growth impediments, resulting in postnatal mortality, and deviations in brain, heart, and lung development at the pathological and single-cell transcriptome levels. This mutation induces premature stem-like cell exhaustion, accelerates the maturation of GABAergic neurons, and disrupts neurodevelopmental and synaptic pathways. Additionally, it specifically hinders CTCF binding to peripheral motifs upstream to the core consensus site, causing alterations in local chromatin structure and gene expression, particularly at the clustered protocadherin locus. Comparative analysis using human cortical organoids mirrors the consequences induced by this mutation. In summary, this study elucidates the influence of the CTCFR567W mutation on human neurodevelopmental disorders, paving the way for potential therapeutic interventions.

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

Transgenerational gene expression depends on both underlying DNA sequences and epigenetic modifications. The latter, which can result in transmission of variegated gene expression patterns across multiple generations without DNA alterations, has been termed epigenetic inheritance and has been documented in plants, worms, flies and mammals. Whereas transcription factors binding to cognate DNA sequence elements regulate gene expression, the molecular basis for epigenetic inheritance has been linked to histone and DNA modifications and non-coding RNA. Here we report that mutation of the CCAAT box promoter element abrogates NF-Y binding and disrupts the stable transgenerational expression of an MHC class I transgene. Transgenic mice with a mutated CCAAT box in the MHC class I transgene display variegated expression of the transgene among littermates and progeny in multiple independently derived transgenic lines. After 4 generations, CCAAT mutant transgenic lines derived from a single founder stably displayed distinct patterns of expression. Histone modifications and RNA polymerase II binding correlate with expression of CCAAT mutant transgenic lines, whereas DNA methylation and nucleosome occupancy do not. Mutation of the CCAAT box also results in changes to CTCF binding and DNA looping patterns across the transgene that correlate with expression status. These studies identify the CCAAT promoter element as a regulator of stable transgenerational gene expression such that mutation of the CCAAT box results in variegated transgenerational inheritance. Considering that the CCAAT box is present in 30% of eukaryotic promoters, this study provides insights into how fidelity of gene expression patterns is maintained through multiple generations.

OBJECTIVE: The aim of this study was to explore the effect and mechanism of programmed cell death ligand 1 (PD-L1) in promoting the proliferation and osteo/odontogenic-differentiation of human dental pulp stem cells (hDPSCs) by mediating CCCTC-binding factor (CTCF) expression.
METHODS: The interaction between PD-L1 and CTCF was verified through co-immunoprecipitation. hDPSCs transfected with PD-L1 overexpression and CTCF knockdown vectors were treated with lipopolysaccharide or an osteogenic-inducing medium. Inflammatory cytokines and osteo/odontogenic-differentiation related genes were measured. Osteo/odontogenic-differentiation of hDPSCs was assessed using alkaline phosphatase (ALP) and alizarin red S staining.
RESULTS: Overexpression of PD-L1 inhibited LPS-induced pro-inflammatory cytokine upregulation, cell proliferation, ALP activity, and calcium deposition in hDPSCs and elevated the expression of osteo/odontogenic-differentiation related genes; however, such expression patterns could be reversed by CTCF knockdown. Co-immunoprecipitation results confirmed the binding of PD-L1 to CTCF, indicating that PD-L1 overexpression in hDPSCs increases CTCF expression, thus inhibiting the inflammatory response and increasing osteo/odontogenic-differentiation of hDPSCs.
CONCLUSIONS: PD-L1 overexpression in hDPSCs enhances the proliferation and osteo/odontogenic-differentiation of hDPSCs and inhibit the inflammatory response by upregulating CTCF expression.