Analysis of the vertebrate insulator protein CTCF-binding sites in the human genome

CGLoop: a neural network framework for chromatin loop prediction

BACKGROUND

Chromosomes of species exhibit a variety of high-dimensional organizational features, and chromatin loops, which are fundamental structures in the three-dimensional (3D) structure of the genome. Chromatin loops are visible speckled patterns on Hi-C contact matrix generated by chromosome conformation capture methods. The chromatin loops play an important role in gene expression, and predicting the chromatin loops generated during whole genome interactions is crucial for a deeper understanding of the 3D genome structure and function.

RESULTS

Here, we propose CGLoop, a deep learning based neural network framework that detects chromatin loops in Hi-C contact matrix. CGLoop combines the convolutional neural network (CNN) with Convolutional Block Attention Module (CBAM) and the Bidirectional Gated Recurrent Unit (BiGRU) to capture important features related to chromatin loops by comprehensively analyzing the Hi-C contact matrix, enabling the prediction of candidate chromatin loops. And CGLoop employs a density based clustering method to filter the candidate chromatin loops predicted by the neural network model. Finally, we compared CGloop with other chromatin loops prediction methods on several cell line including GM12878, K562, IMR90, and mESC. The code is available from https://github.com/wllwuliliwll/CGLoop .

CONCLUSIONS

The experimental results show that, loops predicted by CGLoop show high APA scores and there is an enrichment of multiple transcription factors and binding proteins at the predicted loops anchors, which outperforms other methods in terms of accuracy and validity of chromatin loops prediction.

Identifying deleterious noncoding variation through gain and loss of CTCF binding activity

CCCTC binding factor (CTCF) regulates gene expression through DNA binding at thousands of genomic loci. Genetic variation in these CTCF binding sites (CBSs) is an important driver of phenotypic variation, yet extracting those that are likely to have functional consequences in whole-genome sequencing remains challenging. To address this, we develop a hypothesis-driven framework to identify and prioritize CBS variants in gnomAD. We synthesize CTCF's binding patterns at 1,063,878 genomic loci across 214 biological contexts into a summary of binding activity. We find that high binding activity significantly correlates with both conserved nucleotides (Pearson R = 0.35, p < 2.2 × 10-16) and sequences that contain high-quality CTCF binding motifs (Pearson R = 0.63, p = 2.9 × 10-12). We then use binding activity to evaluate high-confidence allelic binding predictions for 1,253,329 single-nucleotide variations (SNVs) in gnomAD that disrupt a CBS. We find a strong, positive relationship between the mutability-adjusted proportion of singletons (MAPS) metric and the loss of CTCF binding at loci with high in vitro activity (Pearson R = 0.74, p < 2.2 × 10-16). To contextualize these findings, we apply MAPS to other functional classes of variation and find that a subset of 339,380 loss of CTCF binding variants is observed as infrequently as missense variants are. This work nominates these thousands of rare, noncoding variants that disrupt CTCF binding for further functional studies while providing a blueprint for prioritizing variation in other transcription factor binding sequences.

CHD7 binds to insulators during neuronal differentiation

Spiral ganglion neurons (SGNs) are crucial for hearing, and the loss of SGNs causes hearing loss. Stem cell-based therapies offer a promising approach for SGN regeneration and require understanding the mechanisms governing SGN differentiation. We investigated the chromatin remodeler CHD7 in neuronal differentiation using immortalized multipotent otic progenitor (iMOP) cells. We demonstrated that CHD7 knockdown impaired neuronal differentiation. Genome-wide analysis revealed CHD7 binding at diverse cis-regulatory elements, with notable enrichment at sites marked by the insulator-binding protein CTCF between topologically associating domains (TADs). Insulators marked by the enrichment of CHD7 and CTCF resided near genes critical for neuronal differentiation, including Mir9-2. Targeting these regulatory regions in iMOPs with CRISPR interference (CRISPRi) and activation (CRISPRa) increased miR-9 transcription, irrespective of the method. Blocking the CHD7 and CTCF marked sites suggested that the elements function as insulators to regulate gene expression. The study highlights CHD7 activity at insulators and underscores an unreported mechanism for promoting neuronal differentiation.

Insulation between adjacent TADs is controlled by the width of their boundaries through distinct mechanisms

Topologically associating domains (TADs) are sub-Megabase regions in vertebrate genomes with enriched intradomain interactions that restrict enhancer-promoter contacts across their boundaries. However, the mechanisms that separate TADs remain incompletely understood. Most boundaries between TADs contain CTCF binding sites (CBSs), which individually contribute to the blocking of Cohesin-mediated loop extrusion. Using genome-wide classification, here we show that the width of TAD boundaries forms a continuum from narrow to highly extended and correlates with CBSs distribution, chromatin features, and gene regulatory elements. To investigate how these boundary widths emerge, we modified the random crosslinker polymer model to incorporate specific boundary configurations, enabling us to evaluate the differential impact of boundary composition on TAD insulation. Our analysis, using three generic boundary categories, identifies differential influence on TAD insulation, with varying local and distal effects on neighboring domains. Notably, we find that increasing boundary width reduces long-range inter-TAD contacts, as confirmed by Hi-C data. While blocking loop extrusion at boundaries indirectly promotes spurious intermingling of neighboring TADs, extended boundaries counteract this effect, emphasizing their role in establishing genome organization. In conclusion, TAD boundary width not only enhances the efficiency of loop extrusion blocking but may also modulate enhancer-promoter contacts over long distances across TAD boundaries, providing a further mechanism for transcriptional regulation.

Systematic evaluation of single-cell multimodal data integration for comprehensive human reference atlas

The integration of multimodal single-cell data enables comprehensive organ reference atlases, yet its impact remains largely unexplored, particularly in complex tissues. We generated a benchmarking dataset for the renal cortex by integrating 3' and 5' scRNA-seq with joint snRNA-seq and snATAC-seq, profiling 119,744 high-quality nuclei/cells from 19 donors. To align cell identities and enable consistent comparisons, we developed the interpretable machine learning tool scOMM (single-cell Omics Multimodal Mapping) and systematically assessed integration strategies. "Horizontal" integration of scRNA and snRNA-seq improved cell-type identification, while "vertical" integration of snRNA-seq and snATAC-seq had an additive effect, enhancing resolution in homogeneous populations and difficult-to-identify states. Global integration was especially effective in identifying adaptive states and rare cell types, including WFDC2-expressing Thick Ascending Limb and Norn cells, previously undetected in kidney atlases. Our work establishes a robust framework for multimodal reference atlas generation, advancing single-cell analysis and extending its applicability to diverse tissues.

Integrative analysis of the 3D genome and epigenome in mouse embryonic tissues

While a rich set of putative cis-regulatory sequences involved in mouse fetal development have been annotated recently on the basis of chromatin accessibility and histone modification patterns, delineating their role in developmentally regulated gene expression continues to be challenging. To fill this gap, here we mapped chromatin contacts between gene promoters and distal sequences across the genome in seven mouse fetal tissues and across six developmental stages of the forebrain. We identified 248,620 long-range chromatin interactions centered at 14,138 protein-coding genes and characterized their tissue-to-tissue variations and developmental dynamics. Integrative analysis of the interactome with previous epigenome and transcriptome datasets from the same tissues revealed a strong correlation between the chromatin contacts and chromatin state at distal enhancers, as well as gene expression patterns at predicted target genes. We predicted target genes of 15,098 candidate enhancers and used them to annotate target genes of homologous candidate enhancers in the human genome that harbor risk variants of human diseases. We present evidence that schizophrenia and other adult disease risk variants are frequently found in fetal enhancers, providing support for the hypothesis of fetal origins of adult diseases.

HDAC4-AS1/CTCF Transcriptionally Represses HDAC4 Under Stress, Whereas HDAC4 Inhibits Stress-Induced Syncytiotrophoblast Cellular Pyroptosis by Deacetylating NLRP3 and GSDMD

Our previous study reported that histone deacetylase 4 (HDAC4) expression is significantly downregulated in placental tissues of pre-eclampsia (PE) pregnancies. Cellular pyroptosis is a key event in the pathogenesis of PE that induces the release of factors into the maternal circulation. The aim of this study is to analyze the role and related molecular mechanisms of HDAC4 in PE trophoblast cell pyroptosis. Hypoxia and lipopolysaccharide (LPS)/ATP-treated immortalized human placental villous trophoblast cells HTR-8/SVneo were utilized to mimic the placental trophoblast cell state in PE. Both hypoxia and LPS/ATP treatments induced significant HTR-8/SVneo cell pyroptosis, whereas HDAC4 overexpression inhibited the induced cell pyroptosis. HDAC4 could bind to NLRP3 and GSDMD proteins, and lead to a decrease in acetylated NLRP3 and GSDMD proteins, thereby inhibiting cell pyroptosis. Hypoxia and LPS/ATP treatment significantly upregulated HDAC4-AS1 levels in HRT-8/SVneo cells. HDAC4-AS1 could bind to HDAC4 gene promoter sequences as well as CTCF protein. HDAC4-AS1 overexpression recruited the enrichment of CTCF on HDAC4 promoter sequences and further repressed HDAC4 transcription and expression. Targeting the transcriptional regulatory mechanism of HDAC4-AS1/HDAC4 may be able to ameliorate the clinical symptoms of PE maternal by inhibiting cellular pyroptosis in syncytiotrophoblast cells under stress.

Multi-locus CRISPRi targeting with a single truncated guide RNA

A critical goal in functional genomics is evaluating which non-coding elements contribute to gene expression, cellular function, and disease. Functional characterization remains a challenge due to the abundance and complexity of candidate elements. Here, we develop a CRISPRi-based approach for multi-locus screening of putative transcription factor binding sites with a single truncated guide. A truncated guide with hundreds of sequence match sites can reliably disrupt enhancer activity, which expands the targeting scope of CRISPRi while maintaining repressive efficacy. We screen over 13,000 possible CTCF binding sites with 24 guides at 10 nucleotides in spacer length. These truncated guides direct CRISPRi-mediated deposition of repressive H3K9me3 marks and disrupt transcription factor binding at most sequence match target sites. This approach can be a valuable screening step for testing transcription factor binding motifs or other repeated genomic sequences and is easily implemented with existing tools.

Predicting CTCF cell type active binding sites in human genome

The CCCTC-binding factor (CTCF) is pivotal in orchestrating diverse biological functions across the human genome, yet the mechanisms driving its cell type-active DNA binding affinity remain underexplored. Here, we collected ChIP-seq data from 67 cell lines in ENCODE, constructed a unique dataset of cell type-active CTCF binding sites (CBS), and trained convolutional neural networks (CNN) to dissect the patterns of CTCF binding activity. Our analysis reveals that transcription factors RAD21/SMC3 and chromatin accessibility are more predictive compared to sequence motifs and histone modifications. Integrating them together achieved AUPRC values consistently above 0.868, highlighting their utility in deciphering CTCF transcription factor binding dynamics. This study provides a deeper understanding of the regulatory functions of CTCF via machine learning framework.

GHT-SELEX demonstrates unexpectedly high intrinsic sequence specificity and complex DNA binding of many human transcription factors

A long-standing challenge in human regulatory genomics is that transcription factor (TF) DNA-binding motifs are short and degenerate, while the genome is large. Motif scans therefore produce many false-positive binding site predictions. By surveying 179 TFs across 25 families using >1,500 cyclic in vitro selection experiments with fragmented, naked, and unmodified genomic DNA - a method we term GHT-SELEX (Genomic HT-SELEX) - we find that many human TFs possess much higher sequence specificity than anticipated. Moreover, genomic binding regions from GHT-SELEX are often surprisingly similar to those obtained in vivo (i.e. ChIP-seq peaks). We find that comparable specificity can also be obtained from motif scans, but performance is highly dependent on derivation and use of the motifs, including accounting for multiple local matches in the scans. We also observe alternative engagement of multiple DNA-binding domains within the same protein: long C2H2 zinc finger proteins often utilize modular DNA recognition, engaging different subsets of their DNA binding domain (DBD) arrays to recognize multiple types of distinct target sites, frequently evolving via internal duplication and divergence of one or more DBDs. Thus, contrary to conventional wisdom, it is common for TFs to possess sufficient intrinsic specificity to independently delineate cellular targets.

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.

Pushing the TAD boundary: Decoding insulator codes of clustered CTCF sites in 3D genomes

Topologically associating domain (TAD) boundaries are the flanking edges of TADs, also known as insulated neighborhoods, within the 3D structure of genomes. A prominent feature of TAD boundaries in mammalian genomes is the enrichment of clustered CTCF sites often with mixed orientations, which can either block or facilitate enhancer-promoter (E-P) interactions within or across distinct TADs, respectively. We will discuss recent progress in the understanding of fundamental organizing principles of the clustered CTCF insulator codes at TAD boundaries. Specifically, both inward- and outward-oriented CTCF sites function as topological chromatin insulators by asymmetrically blocking improper TAD-boundary-crossing cohesin loop extrusion. In addition, boundary stacking and enhancer clustering facilitate long-distance E-P interactions across multiple TADs. Finally, we provide a unified mechanism for RNA-mediated TAD boundary function via R-loop formation for both insulation and facilitation. This mechanism of TAD boundary formation and insulation has interesting implications not only on how the 3D genome folds in the Euclidean nuclear space but also on how the specificity of E-P interactions is developmentally regulated.

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.

A self-complementary AAV proviral plasmid that reduces cross-packaging and ITR promoter activity in AAV vector preparations

Adeno-associated viral vectors (AAVs) are a leading delivery system for gene therapy in animal models and humans. With several Food and Drug Administration-approved AAV gene therapies on the market, issues related to vector manufacturing have become increasingly important. In this study, we focused on potentially toxic DNA contaminants that can arise from AAV proviral plasmids, the raw materials required for manufacturing recombinant AAV in eukaryotic cells. Typical AAV proviral plasmids are circular DNAs containing a therapeutic gene cassette flanked by natural AAV inverted terminal repeat (ITR) sequences, and a plasmid backbone carrying prokaryotic sequences required for plasmid replication and selection in bacteria. While the majority of AAV particles package the intended therapeutic payload, some capsids instead package the bacterial sequences located on the proviral plasmid backbone. Since ITR sequences also have promoter activity, potentially toxic bacterial open reading frames can be produced in vivo, thereby representing a safety risk. In this study, we describe a new AAV proviral plasmid for vector manufacturing that (1) significantly decreases cross-packaged bacterial sequences, (2) increases correctly packaged AAV payloads, and (3) blunts ITR-driven transcription of cross-packaged material to avoid expressing potentially toxic bacterial sequences. This system may help improve the safety of AAV vector products.

Cell-free DNA from germline TP53 mutation carriers reflect cancer-like fragmentation patterns

Germline pathogenic TP53 variants predispose individuals to a high lifetime risk of developing multiple cancers and are the hallmark feature of Li-Fraumeni syndrome (LFS). Our group has previously shown that LFS patients harbor shorter plasma cell-free DNA fragmentation; independent of cancer status. To understand the functional underpinning of cfDNA fragmentation in LFS, we conducted a fragmentomic analysis of 199 cfDNA samples from 82 TP53 mutation carriers and 30 healthy TP53-wildtype controls. We find that LFS individuals exhibit an increased prevalence of A/T nucleotides at fragment ends, dysregulated nucleosome positioning at p53 binding sites, and loci-specific changes in chromatin accessibility at development-associated transcription factor binding sites and at cancer-associated open chromatin regions. Machine learning classification resulted in robust differentiation between TP53 mutant versus wildtype cfDNA samples (AUC-ROC = 0.710-1.000) and intra-patient longitudinal analysis of ctDNA fragmentation signal enabled early cancer detection. These results suggest that cfDNA fragmentation may be a useful diagnostic tool in LFS patients and provides an important baseline for cancer early detection.

guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

MYC Binding Near Transcriptional End Sites Regulates Basal Gene Expression, Read-Through Transcription and Intragenic Contacts

The MYC oncoprotein regulates numerous genes involved in cellular processes such as cell cycle and mitochondrial and ribosomal structure and function. This requires heterodimerization with its partner, MAX, and binding to specific promoter and enhancer elements. Here, we show that MYC and MAX also bind near transcriptional end sites (TESs) of over one-sixth of all annotated genes. These interactions are dose-dependent, evolutionarily conserved, stabilize the normally short-lived MYC protein and regulate expression both in concert with and independent of MYC's binding elsewhere. MYC's TES binding occurs in association with other transcription factors, alters the chromatin landscape, increases nuclease susceptibility and can alter transcriptional read-through, particularly in response to certain stresses. MYC-bound TESs can directly contact promoters and may fine-tune gene expression in response to both physiologic and pathologic stimuli. Collectively, these findings support a previously unrecognized role for MYC in regulating transcription and its read-through via direct intragenic contacts between TESs and promoters.

Endonucleosis mediates internalization of cytoplasm into the nucleus

Setd8 regulates transcription elongation, mitotic DNA condensation, DNA damage response and replication licensing. Here we show that, in mitogen-stimulated liver-specific Setd8-KO mice, most of the hepatocytes are eliminated by necrosis but a significant number of them survive via entering a stage exhibiting several senescence-related features. Setd8-deficient hepatocytes had enlarged nuclei, chromosomal hyperploidy and nuclear engulfments progressing to the formation of intranuclear vesicles surrounded by nuclear lamina. These vesicles contain glycogen, cytoplasmic proteins and even entire organelles. We term this process "endonucleosis". Intranuclear vesicles are absent in hepatocytes of Setd8/Atg5 knockout mice, suggesting that the process requires the function of the canonical autophagy machinery. Endonucleosis and hyperploidization are temporary, early events in the surviving Setd8-deficient cells. Larger vesicles break down into microvesicles over time and are eventually eliminated. The results reveal sequential events in cells with extensive DNA damage, which function as part of survival mechanisms to prevent necrotic death.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Intertwining roles of R-loops and G-quadruplexes in DNA repair, transcription and genome organization

R-loops are three-stranded nucleic acid structures that are abundant and widespread across the genome and that have important physiological roles in many nuclear processes. Their accumulation is observed in cancers and neurodegenerative disorders. Recent studies have implicated a function for R-loops and G-quadruplex (G4) structures, which can form on the displaced single strand of R-loops, in three-dimensional genome organization in both physiological and pathological contexts. Here we discuss the interconnected functions of DNA:RNA hybrids and G4s within R-loops, their impact on DNA repair and gene regulatory networks, and their emerging roles in genome organization during development and disease.

Single-molecule imaging reveals a direct role of CTCF's zinc fingers in SA interaction and cluster-dependent RNA recruitment

CTCF is a zinc finger protein associated with transcription regulation that also acts as a barrier factor for topologically associated domains (TADs) generated by cohesin via loop extrusion. These processes require different properties of CTCF-DNA interaction, and it is still unclear how CTCF's structural features may modulate its diverse roles. Here, we employ single-molecule imaging to study both full-length CTCF and truncation mutants. We show that CTCF enriches at CTCF binding sites (CBSs), displaying a longer lifetime than observed previously. We demonstrate that the zinc finger domains mediate CTCF clustering and that clustering enables RNA recruitment, possibly creating a scaffold for interaction with RNA-binding proteins like cohesin's subunit SA. We further reveal a direct recruitment and an increase of SA residence time by CTCF bound at CBSs, suggesting that CTCF-SA interactions are crucial for cohesin stability on chromatin at TAD borders. Furthermore, we establish a single-molecule T7 transcription assay and show that although a transcribing polymerase can remove CTCF from CBSs, transcription is impaired. Our study shows that context-dependent nucleic acid binding determines the multifaceted CTCF roles in genome organization and transcription regulation.

Massively parallel dissection of RNA in RNA-protein interactions in vivo

Many of the biological functions performed by RNA are mediated by RNA-binding proteins (RBPs), and understanding the molecular basis of these interactions is fundamental to biology. Here, we present massively parallel RNA assay combined with immunoprecipitation (MPRNA-IP) for in vivo high-throughput dissection of RNA-protein interactions and describe statistical models for identifying RNA domains and parsing the structural contributions of RNA. By using custom pools of tens of thousands of RNA sequences containing systematically designed truncations and mutations, MPRNA-IP is able to identify RNA domains, sequences, and secondary structures necessary and sufficient for protein binding in a single experiment. We show that this approach is successful for multiple RNAs of interest, including the long noncoding RNA NORAD, bacteriophage MS2 RNA, and human telomerase RNA, and we use it to interrogate the hitherto unknown sequence or structural RNA-binding preferences of the DNA-looping factor CTCF. By integrating systematic mutation analysis with crosslinking immunoprecipitation, MPRNA-IP provides a novel high-throughput way to elucidate RNA-based mechanisms behind RNA-protein interactions in vivo.

Atypical Modes of CTCF Binding Facilitate Tissue-Specific and Neuronal Activity-Dependent Gene Expression States

Multivalent binding of CTCF to variable DNA sequences is thought to underlie its ability to mediate diverse cellular functions. CTCF typically binds a 20 base-pair consensus DNA sequence, but the full diversity of CTCF binding sites (CBS) within the genome has not been interrogated. We assessed CTCF occupancy in cultured cortical neurons and observed surprisingly that ~ 22% of CBS lack the consensus CTCF motif. We report here that sequence diversity at most of these atypical CBS involves degeneracy at specific nucleotide positions within the consensus CTCF motif, which likely affect the binding of CTCF zinc fingers 6 and 7. This mode of atypical CTCF binding defines most CBS at gene promoters, as well as CBS that are dynamically altered during neural differentiation and following neuronal stimulation, revealing how atypical CTCF binding could influence gene activity. Dynamic CBS are distributed both within and outside loop anchors and TAD boundaries, suggesting both looping-dependent and independent roles for CTCF. Finally, we describe a second mode of atypical CTCF binding to DNA sequences that are completely unrelated to the consensus CTCF motif, which are enriched within the bodies of tissue-specific genes. These tissue-specific atypical CBS are also enriched in H3K27ac, which marks cis-regulatory elements within chromatin, including enhancers. Overall, these results indicate how atypical CBS could dynamically regulate gene activity patterns during differentiation, development, and in response to environmental cues.

Autophagy related gene 5 polymorphism rs17587319 (C/G) in asthmatic patients in North Indian population

Objective: Genetic background and environmental stimuli play an important role in asthma, which is an individual's hyper-responsiveness to these stimuli leading to airway inflammation. Autophagy Related Gene 5 (ATG5) plays a critical role in the autophagy pathway and has been shown to be involved in asthma. The genetic polymorphisms in the ATG5 have been reported to predispose individuals to asthma. The role of single nucleotide polymorphism rs17587319 (C/G) of ATG5 in asthma has not been studied so far. Materials and methods: In this study, we in silico analysed rs17587319 (C/G) using web-based tools Human Splice Finder (HSF) and RegulomeDB and further a case-control study was conducted that included 187 blood samples (94 asthmatic and 93 healthy controls). Results: In silico analysis suggested alteration of splicing signals by this intronic variant. The samples were genotyped by applying the PCR-RFLP method. The MAF obtained was 0.022 and 0.043 in healthy controls and asthmatic individuals, respectively. The statistical analysis revealed no association (allelic model, OR = 2.02, 95%CI = 0.59-6.83, p = 0.25; co-dominant model, OR = 2.06, 95%CI = 0.6-7.12, p = 0.24) of rs17587319 (C/G) with the susceptibility to asthma in the north Indian population. Conclusions: In conclusion, rs17587319 (C/G) of ATG5 does not predispose individuals to asthma in our part of the world. Further studies are needed including more number of samples to ascertain the role of this polymorphism in asthma.

The PTM profiling of CTCF reveals the regulation of 3D chromatin structure by O-GlcNAcylation

CCCTC-binding factor (CTCF), a ubiquitously expressed and highly conserved protein, is known to play a critical role in chromatin structure. Post-translational modifications (PTMs) diversify the functions of protein to regulate numerous cellular processes. However, the effects of PTMs on the genome-wide binding of CTCF and the organization of three-dimensional (3D) chromatin structure have not been fully understood. In this study, we uncovered the PTM profiling of CTCF and demonstrated that CTCF can be O-GlcNAcylated and arginine methylated. Functionally, we demonstrated that O-GlcNAcylation inhibits CTCF binding to chromatin. Meanwhile, deficiency of CTCF O-GlcNAcylation results in the disruption of loop domains and the alteration of chromatin loops associated with cellular development. Furthermore, the deficiency of CTCF O-GlcNAcylation increases the expression of developmental genes and negatively regulates maintenance and establishment of stem cell pluripotency. In conclusion, these results provide key insights into the role of PTMs for the 3D chromatin structure.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

Gastrulation-stage gene expression in Nipbl+/- mouse embryos foreshadows the development of syndromic birth defects

In animal models, Nipbl deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange syndrome, the most common cause of which is Nipbl haploinsufficiency. Previous studies in Nipbl+/- mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA sequencing on wild-type and Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than wild-type and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and Cornelia de Lange syndrome.

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

Logical design of synthetic cis-regulatory DNA for genetic tracing of cell identities and state changes

Descriptive data are rapidly expanding in biomedical research. Instead, functional validation methods with sufficient complexity remain underdeveloped. Transcriptional reporters allow experimental characterization and manipulation of developmental and disease cell states, but their design lacks flexibility. Here, we report logical design of synthetic cis-regulatory DNA (LSD), a computational framework leveraging phenotypic biomarkers and trans-regulatory networks as input to design reporters marking the activity of selected cellular states and pathways. LSD uses bulk or single-cell biomarkers and a reference genome or custom cis-regulatory DNA datasets with user-defined boundary regions. By benchmarking validated reporters, we integrate LSD with a computational ranking of phenotypic specificity of putative cis-regulatory DNA. Experimentally, LSD-designed reporters targeting a wide range of cell states are functional without minimal promoters. Applied to broadly expressed genes from human and mouse tissues, LSD generates functional housekeeper-like sLCRs compatible with size constraints of AAV vectors for gene therapy applications. A mesenchymal glioblastoma reporter designed by LSD outperforms previously validated ones and canonical cell surface markers. In genome-scale CRISPRa screens, LSD facilitates the discovery of known and novel bona fide cell-state drivers. Thus, LSD captures core principles of cis-regulation and is broadly applicable to studying complex cell states and mechanisms of transcriptional regulation.

Dynamic methylation pattern of H19DMR and KvDMR1 in bovine oocytes and preimplantation embryos

PURPOSE

This study aimed to evaluate the epigenetic reprogramming of ICR1 (KvDMR1) and ICR2 (H19DMR) and expression of genes controlled by them as well as those involved in methylation, demethylation, and pluripotency.

METHODS

We collected germinal vesicle (GV) and metaphase II (MII) oocytes, and preimplantation embryos at five stages [zygote, 4-8 cells, 8-16 cells, morula, and expanded blastocysts (ExB)]. DNA methylation was assessed by BiSeq, and the gene expression was evaluated using qPCR.

RESULTS

H19DMR showed an increased DNA methylation from GV to MII oocytes (68.04% and 98.05%, respectively), decreasing in zygotes (85.83%) until morula (61.65%), and ExB (63.63%). H19 and IGF2 showed increased expression in zygotes, which decreased in further stages. KvDMR1 was hypermethylated in both GV (71.82%) and MII (69.43%) and in zygotes (73.70%) up to morula (77.84%), with a loss of methylation at the ExB (36.64%). The zygote had higher expression of most genes, except for CDKN1C and PHLDA2, which were highly expressed in MII and GV oocytes, respectively. DNMTs showed increased expression in oocytes, followed by a reduction in the earliest stages of embryo development. TET1 was downregulated until 4-8-cell and upregulated in 8-16-cell embryos. TET2 and TET3 showed higher expression in oocytes, and a downregulation in MII oocytes and 4-8-cell embryo.

CONCLUSION

We highlighted the heterogeneity in the DNA methylation of H19DMR and KvDMR1 and a dynamic expression pattern of genes controlled by them. The expression of DNMTs and TETs genes was also dynamic owing to epigenetic reprogramming.

Differential DNA methylation and CTCF binding between the ESR1 promoter a of MCF-7 and MDA-MB-231 breast cancer cells

BACKGROUND

ESR1 is expressed by 60-70% of breast tumours. it's a good prognosis factor and the target of hormone therapy. Optimization of ESR1 reactivation therapy is currently ongoing. Here we probe if the transcription factor CTCF plays a role in the differential expression of ESR1 in the breast cancer cell lines MCF-7 (ESR1+) and MDA-MB-231 (ESR1-).

METHODS AND RESULTS

Knockdown of CTCF in MCF-7 resulted in decreased ESR1 gene expression. CTCF binds to the promoter of ESR1 in MCF-7 but not in MDA-MB-231 cells. CTCF ESR1 binding sites are unmethylated in MCF7 but methylated in MDA-MB-231 cells.

CONCLUSION

ESR1 expression in MCF7 cells is dependent on CTCF expression. CTCF can bind to specific regions of the promotor of ESR1 gene in MCF-7 cells but not in MDA-MB-231 cells, this correlates with the methylation status of these regions and could be involved in the transcriptional regulation of ESR1.

The intrinsically disordered tail of ESCO1 binds DNA in a charge-dependent manner

ESCO1 is an acetyltransferase enzyme that regulates chromosome organization and gene expression. It does this by modifying the Smc3 subunit of the Cohesin complex. Although ESCO1 is enriched at the base of chromatin loops in a Cohesin-dependent manner, precisely how it interacts with chromatin is unknown. Here we show that the basic and intrinsically disordered tail of ESCO1 binds DNA with very high affinity, likely through electrostatic interaction. We show that neutralization of positive residues in the N-tail reduces both DNA binding in vitro and association of the enzyme with chromatin in cells. Additionally, disruption of the chromatin state and charge distribution reduces chromatin bound ESCO1. Strikingly, defects in DNA binding do not affect total SMC3 acetylation or sister chromatid cohesion, suggesting that ESCO1-dependent acetylation can occur independently of direct chromatin association. We conclude that the intrinsically disordered tail of ESCO1 binds DNA with both high affinity and turnover, but surprisingly, ESCO1 catalytic activity occurs independently of direct DNA binding by the enzyme.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

CTCF deletion alters the pluripotency and DNA methylation profile of human iPSCs

Pluripotent stem cells are characterized by their differentiation potential toward endoderm, mesoderm, and ectoderm. However, it is still largely unclear how these cell-fate decisions are mediated by epigenetic mechanisms. In this study, we explored the relevance of CCCTC-binding factor (CTCF), a zinc finger-containing DNA-binding protein, which mediates long-range chromatin organization, for directed cell-fate determination. We generated human induced pluripotent stem cell (iPSC) lines with deletions in the protein-coding region in exon 3 of CTCF, resulting in shorter transcripts and overall reduced protein expression. Chromatin immunoprecipitation showed a considerable loss of CTCF binding to target sites. The CTCF deletions resulted in slower growth and modest global changes in gene expression, with downregulation of a subset of pluripotency-associated genes and neuroectodermal genes. CTCF deletion also evoked DNA methylation changes, which were moderately associated with differential gene expression. Notably, CTCF-deletions lead to upregulation of endo-mesodermal associated marker genes and epigenetic signatures, whereas ectodermal differentiation was defective. These results indicate that CTCF plays an important role in the maintenance of pluripotency and differentiation, especially towards ectodermal lineages.

Molecular Mechanisms for the Regulation of Nuclear Membrane Integrity

The nuclear membrane serves a critical role in protecting the contents of the nucleus and facilitating material and signal exchange between the nucleus and cytoplasm. While extensive research has been dedicated to topics such as nuclear membrane assembly and disassembly during cell division, as well as interactions between nuclear transmembrane proteins and both nucleoskeletal and cytoskeletal components, there has been comparatively less emphasis on exploring the regulation of nuclear morphology through nuclear membrane integrity. In particular, the role of type II integral proteins, which also function as transcription factors, within the nuclear membrane remains an area of research that is yet to be fully explored. The integrity of the nuclear membrane is pivotal not only during cell division but also in the regulation of gene expression and the communication between the nucleus and cytoplasm. Importantly, it plays a significant role in the development of various diseases. This review paper seeks to illuminate the biomolecules responsible for maintaining the integrity of the nuclear membrane. It will delve into the mechanisms that influence nuclear membrane integrity and provide insights into the role of type II membrane protein transcription factors in this context. Understanding these aspects is of utmost importance, as it can offer valuable insights into the intricate processes governing nuclear membrane integrity. Such insights have broad-reaching implications for cellular function and our understanding of disease pathogenesis.

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.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Understanding fundamental principles of enhancer biology at a model locus: Analysing the structure and function of an enhancer cluster at the α-globin locus

Despite ever-increasing accumulation of genomic data, the fundamental question of how individual genes are switched on during development, lineage-specification and differentiation is not fully answered. It is widely accepted that this involves the interaction between at least three fundamental regulatory elements: enhancers, promoters and insulators. Enhancers contain transcription factor binding sites which are bound by transcription factors (TFs) and co-factors expressed during cell fate decisions and maintain imposed patterns of activation, at least in part, via their epigenetic modification. This information is transferred from enhancers to their cognate promoters often by coming into close physical proximity to form a 'transcriptional hub' containing a high concentration of TFs and co-factors. The mechanisms underlying these stages of transcriptional activation are not fully explained. This review focuses on how enhancers and promoters are activated during differentiation and how multiple enhancers work together to regulate gene expression. We illustrate the currently understood principles of how mammalian enhancers work and how they may be perturbed in enhanceropathies using expression of the α-globin gene cluster during erythropoiesis, as a model.

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.

G-quadruplexes associated with R-loops promote CTCF binding

CTCF is a critical regulator of genome architecture and gene expression that binds thousands of sites on chromatin. CTCF genomic localization is controlled by the recognition of a DNA sequence motif and regulated by DNA modifications. However, CTCF does not bind to all its potential sites in all cell types, raising the question of whether the underlying chromatin structure can regulate CTCF occupancy. Here, we report that R-loops facilitate CTCF binding through the formation of associated G-quadruplex (G4) structures. R-loops and G4s co-localize with CTCF at many genomic regions in mouse embryonic stem cells and promote CTCF binding to its cognate DNA motif in vitro. R-loop attenuation reduces CTCF binding in vivo. Deletion of a specific G4-forming motif in a gene reduces CTCF binding and alters gene expression. Conversely, chemical stabilization of G4s results in CTCF gains and accompanying alterations in chromatin organization, suggesting a pivotal role for G4 structures in reinforcing long-range genome interactions through CTCF.

The crucial role of CTCF in mitotic progression during early development of sea urchin

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

Molecular mechanisms of multi-omic regulation in breast cancer

Breast cancer is a complex disease that is influenced by the concurrent influence of multiple genetic and environmental factors. Recent advances in genomics and other high throughput biomolecular techniques (-omics) have provided numerous insights into the molecular mechanisms underlying breast cancer development and progression. A number of these mechanisms involve multiple layers of regulation. In this review, we summarize the current knowledge on the role of multiple omics in the regulation of breast cancer, including the effects of DNA methylation, non-coding RNA, and other epigenomic changes. We comment on how integrating such diverse mechanisms is envisioned as key to a more comprehensive understanding of breast carcinogenesis and cancer biology with relevance to prognostics, diagnostics and therapeutics. We also discuss the potential clinical implications of these findings and highlight areas for future research. Overall, our understanding of the molecular mechanisms of multi-omic regulation in breast cancer is rapidly increasing and has the potential to inform the development of novel therapeutic approaches for this disease.

Recent advances in chromosome capture techniques unraveling 3D genome architecture in germ cells, health, and disease

In eukaryotes, the genome does not emerge in a specific shape but rather as a hierarchial bundle within the nucleus. This multifaceted genome organization consists of multiresolution cellular structures, such as chromosome territories, compartments, and topologically associating domains, which are frequently defined by architecture, design proteins including CTCF and cohesin, and chromatin loops. This review briefly discusses the advances in understanding the basic rules of control, chromatin folding, and functional areas in early embryogenesis. With the use of chromosome capture techniques, the latest advancements in technologies for visualizing chromatin interactions come close to revealing 3D genome formation frameworks with incredible detail throughout all genomic levels, including at single-cell resolution. The possibility of detecting variations in chromatin architecture might open up new opportunities for disease diagnosis and prevention, infertility treatments, therapeutic approaches, desired exploration, and many other application scenarios.

Chromatin Organization and Transcriptional Programming of Breast Cancer Cell Identity

The advent of sequencing technologies for assessing chromosome conformations has provided a wealth of information on the organization of the 3-dimensional genome and its role in cancer progression. It is now known that changes in chromatin folding and accessibility can promote aberrant activation or repression of transcriptional programs that can drive tumorigenesis and progression in diverse cancers. This includes breast cancer, which comprises several distinct subtypes defined by their unique transcriptomes that dictate treatment response and patient outcomes. Of these, basal-like breast cancer is an aggressive subtype controlled by a pluripotency-enforcing transcriptome. Meanwhile, the more differentiated luminal subtype of breast cancer is driven by an estrogen receptor-dominated transcriptome that underlies its responsiveness to antihormone therapies and conveys improved patient outcomes. Despite the clear differences in molecular signatures, the genesis of each subtype from normal mammary epithelial cells remains unclear. Recent technical advances have revealed key distinctions in chromatin folding and organization between subtypes that could underlie their transcriptomic and, hence, phenotypic differences. These studies also suggest that proteins controlling particular chromatin states may be useful targets for treating aggressive disease. In this review, we explore the current state of understanding of chromatin architecture in breast cancer subtypes and its potential role in defining their phenotypic characteristics.

On the dependent recognition of some long zinc finger proteins

The human genome contains about 800 C2H2 zinc finger proteins (ZFPs), and most of them are composed of long arrays of zinc fingers. Standard ZFP recognition model asserts longer finger arrays should recognize longer DNA-binding sites. However, recent experimental efforts to identify in vivo ZFP binding sites contradict this assumption, with many exhibiting short motifs. Here we use ZFY, CTCF, ZIM3, and ZNF343 as examples to address three closely related questions: What are the reasons that impede current motif discovery methods? What are the functions of those seemingly unused fingers and how can we improve the motif discovery algorithms based on long ZFPs' biophysical properties? Using ZFY, we employed a variety of methods and find evidence for 'dependent recognition' where downstream fingers can recognize some previously undiscovered motifs only in the presence of an intact core site. For CTCF, high-throughput measurements revealed its upstream specificity profile depends on the strength of its core. Moreover, the binding strength of the upstream site modulates CTCF's sensitivity to different epigenetic modifications within the core, providing new insight into how the previously identified intellectual disability-causing and cancer-related mutant R567W disrupts upstream recognition and deregulates the epigenetic control by CTCF. Our results establish that, because of irregular motif structures, variable spacing and dependent recognition between sub-motifs, the specificities of long ZFPs are significantly underestimated, so we developed an algorithm, ModeMap, to infer the motifs and recognition models of ZIM3 and ZNF343, which facilitates high-confidence identification of specific binding sites, including repeats-derived elements. With revised concept, technique, and algorithm, we can discover the overlooked specificities and functions of those 'extra' fingers, and therefore decipher their broader roles in human biology and diseases.

An updated catalog of CTCF variants associated with neurodevelopmental disorder phenotypes

Introduction

CTCF-related disorder (CRD) is a neurodevelopmental disorder (NDD) caused by monoallelic pathogenic variants in CTCF. The first CTCF variants in CRD cases were documented in 2013. To date, 76 CTCF variants have been further described in the literature. In recent years, due to the increased application of next-generation sequencing (NGS), growing numbers of CTCF variants are being identified, and multiple genotype-phenotype databases cataloging such variants are emerging.

Methods

In this study, we aimed to expand the genotypic spectrum of CRD, by cataloging NDD phenotypes associated with reported CTCF variants. Here, we systematically reviewed all known CTCF variants reported in case studies and large-scale exome sequencing cohorts. We also conducted a meta-analysis using public variant data from genotype-phenotype databases to identify additional CTCF variants, which we then curated and annotated.

Results

From this combined approach, we report an additional 86 CTCF variants associated with NDD phenotypes that have not yet been described in the literature. Furthermore, we describe and explain inconsistencies in the quality of reported variants, which impairs the reuse of data for research of NDDs and other pathologies.

Discussion

From this integrated analysis, we provide a comprehensive and annotated catalog of all currently known CTCF mutations associated with NDD phenotypes, to aid diagnostic applications, as well as translational and basic research.

Benefiting from the intrinsic role of epigenetics to predict patterns of CTCF binding

Motivation

One of the most relevant mechanisms involved in the determination of chromatin structure is the formation of structural loops that are also related with the conservation of chromatin states. Many of these loops are stabilized by CCCTC-binding factor (CTCF) proteins at their base. Despite the relevance of chromatin structure and the key role of CTCF, the role of the epigenetic factors that are involved in the regulation of CTCF binding, and thus, in the formation of structural loops in the chromatin, is not thoroughly understood.

Results

Here we describe a CTCF binding predictor based on Random Forest that employs different epigenetic data and genomic features. Importantly, given the ability of Random Forests to determine the relevance of features for the prediction, our approach also shows how the different types of descriptors impact the binding of CTCF, confirming previous knowledge on the relevance of chromatin accessibility and DNA methylation, but demonstrating the effect of epigenetic modifications on the activity of CTCF. We compared our approach against other predictors and found improved performance in terms of areas under PR and ROC curves (PRAUC-ROCAUC), outperforming current state-of-the-art methods.

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.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

Function and Evolution of the Loop Extrusion Machinery in Animals

Structural maintenance of chromosomes (SMC) complexes are essential proteins found in genomes of all cellular organisms. Essential functions of these proteins, such as mitotic chromosome formation and sister chromatid cohesion, were discovered a long time ago. Recent advances in chromatin biology showed that SMC proteins are involved in many other genomic processes, acting as active motors extruding DNA, which leads to the formation of chromatin loops. Some loops formed by SMC proteins are highly cell type and developmental stage specific, such as SMC-mediated DNA loops required for VDJ recombination in B-cell progenitors, or dosage compensation in Caenorhabditis elegans and X-chromosome inactivation in mice. In this review, we focus on the extrusion-based mechanisms that are common for multiple cell types and species. We will first describe an anatomy of SMC complexes and their accessory proteins. Next, we provide biochemical details of the extrusion process. We follow this by the sections describing the role of SMC complexes in gene regulation, DNA repair, and chromatin topology.