CTCF mediates chromatin looping via N-terminal domain-dependent cohesin retention

Zinc finger homeobox-3 (ZFHX3) orchestrates genome-wide daily gene expression in the suprachiasmatic nucleus

The mammalian suprachiasmatic nucleus (SCN), situated in the ventral hypothalamus, directs daily cellular and physiological rhythms across the body. The SCN clockwork is a self-sustaining transcriptional-translational feedback loop (TTFL) that in turn coordinates the expression of clock-controlled genes (CCGs) directing circadian programmes of SCN cellular activity. In the mouse, the transcription factor, ZFHX3 (zinc finger homeobox-3), is necessary for the development of the SCN and influences circadian behaviour in the adult. The molecular mechanisms by which ZFHX3 affects the SCN at transcriptomic and genomic levels are, however, poorly defined. Here, we used chromatin immunoprecipitation sequencing to map the genomic localization of ZFHX3-binding sites in SCN chromatin. To test for function, we then conducted comprehensive RNA sequencing at six distinct times-of-day to compare the SCN transcriptional profiles of control and ZFHX3-conditional null mutants. We show that the genome-wide occupancy of ZFHX3 occurs predominantly around gene transcription start sites, co-localizing with known histone modifications, and preferentially partnering with clock transcription factors (CLOCK, BMAL1) to regulate clock gene(s) transcription. Correspondingly, we show that the conditional loss of ZFHX3 in the adult has a dramatic effect on the SCN transcriptome, including changes in the levels of transcripts encoding elements of numerous neuropeptide neurotransmitter systems while attenuating the daily oscillation of the clock TF Bmal1. Furthermore, various TTFL genes and CCGs exhibited altered circadian expression profiles, consistent with an advanced in daily behavioural rhythms under 12 h light-12 h dark conditions. Together, these findings reveal the extensive genome-wide regulation mediated by ZFHX3 in the central clock that orchestrates daily timekeeping in mammals.

CTCF-anchored chromatin loop dynamics during human meiosis

BACKGROUND

During meiosis, the mammalian genome is organised within chromatin loops, which facilitate synapsis, crossing over and chromosome segregation, setting the stage for recombination events and the generation of genetic diversity. Chromatin looping is thought to play a major role in the establishment of cross overs during prophase I of meiosis, in diploid early primary spermatocytes. However, chromatin conformation dynamics during human meiosis are difficult to study experimentally, due to the transience of each cell division and the difficulty of obtaining stage-resolved cell populations. Here, we employed a machine learning framework trained on single cell ATAC-seq and RNA-seq data to predict CTCF-anchored looping during spermatogenesis, including cell types at different stages of meiosis.

RESULTS

We find dramatic changes in genome-wide looping patterns throughout meiosis: compared to pre-and-post meiotic germline cell types, loops in meiotic early primary spermatocytes are more abundant, more variable between individual cells, and more evenly spread throughout the genome. In preparation for the first meiotic division, loops also include longer stretches of DNA, encompassing more than half of the total genome. These loop structures then influence the rate of recombination initiation and resolution as cross overs. In contrast, in later mature sperm stages, we find evidence of genome compaction, with loops being confined to the telomeric ends of the chromosomes.

CONCLUSION

Overall, we find that chromatin loops do not orchestrate the gene expression dynamics seen during spermatogenesis, but loops do play important roles in recombination, influencing the positions of DNA breakage and cross over events.

CCCTC-binding factor N-terminal domain regulates clustered protocadherin gene expression by enhancing cohesin processivity

CTCF (CCCTC-binding factor) instructs 3D genome folding by anchoring or forestalling cohesin loop extrusion, but the exact mechanism remains obscure. Here, using clustered protocadherins (cPcdh) as model genes, we report that CTCF assists or facilitates cohesin loop extrusion by enhancing its processivity. Specifically, we show that, compared with the Pcdh α and γ gene clusters, the Pcdhβ cluster is greatly affected upon CTCFY226A/F228A mutation in the N-terminal domain. Given the long-range distance of the Pcdhβ cluster from the distal enhancer, this finding has interesting implications in CTCF regulation of cohesin processivity along the linear chromatin during DNA loop extrusion. In particular, the effect on cohesin processivity upon CTCFY226A/F228A mutation is conspicuously similar to that of WAPL overexpression, suggesting that, in contrast to the general view of blocking or forestalling cohesin, CTCF may actually enhance or facilitate cohesin loop extrusion during 3D genome folding. We conclude that CTCF enhances cohesin enrichments via the N-terminal YDF motif in clustered protocadherin genes in a genomic-distance biased manner.

Increasingly efficient chromatin binding of cohesin and CTCF supports chromatin architecture formation during zebrafish embryogenesis

The three-dimensional folding of chromosomes is essential for nuclear functions such as DNA replication and gene regulation. The emergence of chromatin architecture is thus an important process during embryogenesis. To shed light on the molecular and kinetic underpinnings of chromatin architecture formation, we characterized biophysical properties of cohesin and CTCF binding to chromatin and their changes upon cofactor depletion using single-molecule imaging in live developing zebrafish embryos. We found that chromatin-bound fractions of both cohesin and CTCF increased significantly between the 1000-cell and shield stages, which we could explain through changes in both their association and dissociation rates. Moreover, increasing binding of cohesin restricted chromatin motion, potentially via loop extrusion, and showed distinct stage-dependent nuclear distribution. Polymer simulations with experimentally derived parameters recapitulated the experimentally observed gradual emergence of chromatin architecture. Our findings reveal molecular kinetics underlying chromatin architecture formation during zebrafish embryogenesis.

CTCF/RAD21 organize the ground state of chromatin-nuclear speckle association

Recent findings indicate that nuclear speckles, a distinct type of nuclear body, interact with certain chromatin regions in a ground state. Here, we report that the chromatin structural factors CTCF and cohesin are required for full ground-state association between DNA and nuclear speckles. We identified a putative speckle-targeting motif (STM) within cohesin subunit RAD21 and demonstrated that the STM is required for chromatin-nuclear speckle association, disruption of which also impaired induction of speckle-associated genes. Depletion of the cohesin-releasing factor WAPL, which stabilizes cohesin on chromatin, resulted in reinforcement of DNA-speckle contacts and enhanced inducibility of speckle-associated genes. Additionally, we observed disruption of chromatin-nuclear speckle association in patient-derived cells with Cornelia de Lange syndrome, a congenital neurodevelopmental disorder involving defective cohesin pathways. In summary, our findings reveal a mechanism for establishing the ground state of chromatin-speckle association and promoting gene inducibility, with relevance to human disease.

Binding domain mutations provide insight into CTCF's relationship with chromatin and its contribution to gene regulation

Here we used a series of CTCF mutations to explore CTCF's relationship with chromatin and its contribution to gene regulation. CTCF's impact depends on the genomic context of bound sites and the unique binding properties of WT and mutant CTCF proteins. Specifically, CTCF's signal strength is linked to changes in accessibility, and the ability to block cohesin is linked to its binding stability. Multivariate modelling reveals that both CTCF and accessibility contribute independently to cohesin binding and insulation, however CTCF signal strength has a stronger effect. CTCF and chromatin have a bidirectional relationship such that at CTCF sites, accessibility is reduced in a cohesin-dependent, mutant specific fashion. In addition, each mutant alters TF binding and accessibility in an indirect manner, changes which impart the most influence on rewiring transcriptional networks and the cell's ability to differentiate. Collectively, the mutant perturbations provide a rich resource for determining CTCF's site-specific effects.

A negatively charged region within carboxy-terminal domain maintains proper CTCF DNA binding

As an essential regulator of higher-order chromatin structures, CCCTC-binding factor (CTCF) is a highly conserved protein with a central DNA-binding domain of 11 tandem zinc fingers (ZFs), which are flanked by amino (N-) and carboxy (C-) terminal domains of intrinsically disordered regions. Here we report that CRISPR deletion of the entire C-terminal domain of alternating charge blocks decreases CTCF DNA binding but deletion of the C-terminal fragment of 116 amino acids results in increased CTCF DNA binding and aberrant gene regulation. Through a series of genetic targeting experiments, in conjunction with electrophoretic mobility shift assay (EMSA), circularized chromosome conformation capture (4C), qPCR, chromatin immunoprecipitation with sequencing (ChIP-seq), and assay for transposase-accessible chromatin with sequencing (ATAC-seq), we uncovered a negatively charged region (NCR) responsible for weakening CTCF DNA binding and chromatin accessibility. AlphaFold prediction suggests an autoinhibitory mechanism of CTCF via NCR as a flexible DNA mimic domain, possibly competing with DNA binding for the positively charged ZF surface area. Thus, the unstructured C-terminal domain plays an intricate role in maintaining proper CTCF-DNA interactions and 3D genome organization.

A genome wide code to define cell-type specific CTCF binding and chromatin organization

CTCF-mediated chromatin folding plays a key role in gene regulation, however the mechanisms underlying cell type-specific control are not fully elucidated. Comprehensive analyses reveal that CTCF binding stability and cohesin overlap in mice and humans is regulated by species specific differences in accessibility, the presence of CTCF extended Up and Downstream binding site sequences and motifs corresponding to expressed TFs enriched at most bound sites. CTCF repositions nucleosomes asymmetrically, with linker spacing altered by accessibility, while cohesin-mediated nucleosome phasing is affected by surrounding motifs that similarly impact chromatin insulation. Importantly, cell type-specific transcriptional programs determine motif enrichment at CTCF bound sites, reflecting the stabilizing/destabilizing effect of individual TFs. These studies identify mechanisms underlying cell type-specific CTCF profiles, linked to local and long-range chromatin organization.

The role of high mobility group proteins in cellular senescence mechanisms

Aging is a universal physiological phenomenon, and chronic age-related diseases have become one of the leading causes of human mortality, accounting for nearly half of all deaths. Studies have shown that reducing the incidence of these diseases can not only extend lifespan but also promote healthy aging. In recent years, the potential role of non-histone high-mobility group proteins (HMGs) in the regulation of aging and lifespan has attracted widespread attention. HMGs play critical roles in cellular senescence and associated diseases through various pathways, encompassing multi-layered mechanisms involving protein interactions, molecular regulation, and chromatin dynamics. This review provides a comprehensive analysis of the interactions between HMG family proteins and senescence-associated secretory phenotype (SASP), chromatin structure, and histone modifications, offering a deeper exploration of the pivotal functions and impacts of HMGs in the aging process. Furthermore, we summarize recent findings on the contributions of HMG proteins to aging and age-related diseases. HMG proteins not only regulate senescence-associated inflammation through modulating the SASP but also influence genomic stability and cell fate decisions via interactions with chromatin and histones. Targeting HMG proteins holds great potential in delaying the progression of aging and its associated diseases. This review aims to provide a systematic overview of HMG proteins' roles in aging and to lay a solid foundation for future anti-aging drug development and therapeutic strategies. With the advancing understanding of the mechanisms by which HMGs regulate aging, developing therapeutic interventions targeting HMGs may emerge as a promising approach to extending lifespan and enhancing healthspan.

Interaction of CTCF and CTCFL in genome regulation through chromatin architecture during the spermatogenesis and carcinogenesis

The zinc finger protein CTCF is ubiquitously expressed and is integral to the regulation of chromatin architecture through its interaction with cohesin. Conversely, CTCFL expression is predominantly restricted to the adult male testis but is aberrantly expressed in certain cancers. Despite their distinct expression patterns, the cooperative and competitive mechanisms by which CTCF and CTCFL regulate target gene expression in spermatocytes and cancer cells remain inadequately understood. In this review, we comprehensively examine the literature on the divergent amino acid sequences, target sites, expression profiles and functions of CTCF and CTCFL in normal tissues and cancers. We further elucidate the mechanisms by which CTCFL competitively or cooperatively binds to CTCF target sites during spermatogenesis and carcinogenesis to modulate chromatin architecture. We mainly focus on the role of CTCFL in testicular and cancer development, highlighting its interaction with CTCF at CTCF binding sites to regulate target genes. In the testis, CTCF and CTCFL cooperate to regulate the expression of testis-specific genes, essential for proper germ cell progression. In cancers, CTCFL overexpression competes with CTCF for DNA binding, leading to aberrant gene expression, a more relaxed chromatin state, and altered chromatin loops. By uncovering the roles of CTCF and CTCFL in spermatogenesis and carcinogenesis, we can better understand the implications of aberrant CTCFL expression in altering chromatin loops and its contribution to disease pathogenesis.

Dynamic barriers modulate cohesin positioning and genome folding at fixed occupancy

In mammalian interphase cells, genomes are folded by cohesin loop extrusion limited by directional CTCF barriers. This interplay leads to the enrichment of cohesin at barriers, isolation between neighboring topologically associating domains, and elevated contact frequency between convergent CTCF barriers across the genome. However, recent in vivo measurements present a puzzle: reported residence times for CTCF on chromatin are in the range of a few minutes, while lifetimes for cohesin are much longer. Can the observed features of genome folding result from the action of relatively transient barriers? To address this question, we developed a dynamic barrier model, where CTCF sites switch between bound and unbound states with rates that can be directly compared with biophysical measurements. Using this model, we investigated how barrier dynamics would impact observables for a range of experimental genomic and imaging datasets, including ChIP-seq, Hi-C, and microscopy. We found the interplay of CTCF and cohesin binding timescales influence the strength of each of these features, leaving a signature of barrier dynamics even in the population-averaged snapshots offered by genomic datasets. First, in addition to barrier occupancy, barrier bound times are crucial for instructing features of genome folding. Second, the ratio of boundary to extruder lifetime greatly alters simulated ChIP-seq and simulated Hi-C. Third, large-scale changes in chromosome morphology observed experimentally after increasing extruder lifetime require dynamic barriers. By integrating multiple sources of experimental data, our biophysical model argues that CTCF barrier bound times effectively approach those of cohesin extruder lifetimes. Together, we demonstrate how models that are informed by biophysically measured protein dynamics broaden our understanding of genome folding.

Multifaceted role of CTCF in X-chromosome inactivation

Therian female mammals compensate for the dosage of X-linked gene expression by inactivating one of the X-chromosomes. X-inactivation is facilitated by the master regulator Xist long non-coding RNA, which coats the inactive-X and facilitates heterochromatinization through recruiting different chromatin modifiers and changing the X-chromosome 3D conformation. However, many mechanistic aspects behind the X-inactivation process remain poorly understood. Among the many contributing players, CTCF has emerged as one of the key players in orchestrating various aspects related to X-chromosome inactivation by interacting with several other protein and RNA partners. In general, CTCF is a well-known architectural protein, which plays an important role in chromatin organization and transcriptional regulation. Here, we provide significant insight into the role of CTCF in orchestrating X-chromosome inactivation and highlight future perspectives.

Factors that determine cell type-specific CTCF binding in health and disease

A number of factors contribute to cell type-specific CTCF chromatin binding, but how they act in concert to determine binding stability and functionality has not been fully elucidated. In this review, we tie together different layers of regulation to provide a holistic view of what is known. What emerges from these studies is a multifaceted system in which DNA sequence, DNA and chromatin accessibility, and cell type-specific transcription factors together contribute to CTCF binding profile and function. We discuss these findings in the light of disease settings in which changes in the chromatin landscape and transcriptional programming can disrupt CTCF's binding profile and involvement in looping.

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.

Enhancing Vaxign-DL for Vaccine Candidate Prediction with added ESM-Generated Features

Many vaccine design programs have been developed, including our own machine learning approaches Vaxign-ML and Vaxign-DL. Using deep learning techniques, Vaxign-DL predicts bacterial protective antigens by calculating 509 biological and biomedical features from protein sequences. In this study, we first used the protein folding ESM program to calculate a set of 1,280 features from individual protein sequences, and then utilized the new set of features separately or in combination with the traditional set of 509 features to predict protective antigens. Our result showed that the usage of ESM-derived features alone was able to accurately predict vaccine antigens with a performance similar to the orginal Vaxign-DL prediction method, and the usage of the combined ESM-derived and orginal Vaxign-DL features significantly improved the prediction performance according to a set of seven scores including specificity, sensitivity, and AUROC. To further evaluate the updated methods, we conducted a Leave-One-Pathogen-Out Validation (LOPOV) study, and found that the usage of ESM-derived features significantly improved the the prediction of vaccine antigens from 10 bacterial pathogens. This research is the first reported study demonstrating the added value of protein folding features for vaccine antigen prediction.

Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

Enhancers on the edge - how the nuclear envelope controls gene regulatory elements

Precise temporal and sequential control of gene expression during development and in response to environmental stimuli requires tight regulation of the physical contact between gene regulatory elements and promoters. Current models describing how the genome folds in 3D space to establish these interactions often ignore the role of the most stable structural nuclear feature - the nuclear envelope. While contributions of 3D folding within/between topologically associated domains (TADs) have been extensively described, mechanical contributions from the nuclear envelope can impact enhancer-promoter interactions both directly and indirectly through influencing intra/inter-TAD interactions. Importantly, these nuclear envelope contributions clearly link this mechanism to development and, when defective, to human disease. Here, we discuss evidence for nuclear envelope regulation of tissue-specific enhancer-promoter pairings, potential mechanisms for this regulation, exciting recent findings that other regulatory elements such as microRNAs and long noncoding RNAs are under nuclear envelope regulation, the possible involvement of condensates, and how disruption of this regulation can lead to disease.

Brain and cancer associated binding domain mutations provide insight into CTCF's relationship with chromatin and its ability to act as a chromatin organizer

Although only a fraction of CTCF motifs are bound in any cell type, and approximately half of the occupied sites overlap cohesin, the mechanisms underlying cell-type specific attachment and ability to function as a chromatin organizer remain unknown. To investigate the relationship between CTCF and chromatin we applied a combination of imaging, structural and molecular approaches, using a series of brain and cancer associated CTCF mutations that act as CTCF perturbations. We demonstrate that binding and the functional impact of WT and mutant CTCF depend not only on the unique properties of each protein, but also on the genomic context of bound sites. Our studies also highlight the reciprocal relationship between CTCF and chromatin, demonstrating that the unique binding properties of WT and mutant proteins have a distinct impact on accessibility, TF binding, cohesin overlap, chromatin interactivity and gene expression programs, providing insight into their cancer and brain related effects.

Connecting Chromatin Structures to Gene Regulation Using Dynamic Polymer Simulations

The transfer of regulatory information between distal loci on chromatin is thought to involve physical proximity, but key biophysical features of these contacts remain unclear. For instance, it is unknown how close and for how long two loci need to be in order to productively interact. The main challenge is that it is currently impossible to measure chromatin dynamics with high spatiotemporal resolution at scale. Polymer simulations provide an accessible and rigorous way to test biophysical models of chromatin regulation, yet there is a lack of simple and general methods for extracting the values of model parameters. Here we adapt the Nelder-Mead simplex optimization algorithm to select the best polymer model matching a given Hi-C dataset, using the MYC locus as an example. The model's biophysical parameters predict a compartmental rearrangement of the MYC locus in leukemia, which we validate with single-cell measurements. Leveraging trajectories predicted by the model, we find that loci with similar Hi-C contact frequencies can exhibit widely different contact dynamics. Interestingly, the frequency of productive interactions between loci exhibits a non-linear relationship with their Hi-C contact frequency when we enforce a specific capture radius and contact duration. These observations are consistent with recent experimental observations and suggest that the dynamic ensemble of chromatin configurations, rather than average contact matrices, is required to fully predict productive long-range chromatin interactions.

CTCF is essential for proper mitotic spindle structure and anaphase segregation

Mitosis is an essential process in which the duplicated genome is segregated equally into two daughter cells. CTCF has been reported to be present in mitosis and has a role in localizing CENP-E, but its importance for mitotic fidelity remains to be determined. To evaluate the importance of CTCF in mitosis, we tracked mitotic behaviors in wild-type and two different CTCF CRISPR-based genetic knockdowns. We find that knockdown of CTCF results in prolonged mitoses and failed anaphase segregation via time-lapse imaging of SiR-DNA. CTCF knockdown did not alter cell cycling or the mitotic checkpoint, which was activated upon nocodazole treatment. Immunofluorescence imaging of the mitotic spindle in CTCF knockdowns revealed disorganization via tri/tetrapolar spindles and chromosomes behind the spindle pole. Imaging of interphase nuclei showed that nuclear size increased drastically, consistent with failure to divide the duplicated genome in anaphase. Long-term inhibition of CNEP-E via GSK923295 recapitulates CTCF knockdown abnormal mitotic spindles with polar chromosomes and increased nuclear sizes. Population measurements of nuclear shape in CTCF knockdowns do not display decreased circularity or increased nuclear blebbing relative to wild-type. However, failed mitoses do display abnormal nuclear morphologies relative to successful mitoses, suggesting that population images do not capture individual behaviors. Thus, CTCF is important for both proper metaphase organization and anaphase segregation which impacts the size and shape of the interphase nucleus likely through its known role in recruiting CENP-E.

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.

Strand-resolved mutagenicity of DNA damage and repair

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

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.

Attraction and disruption: how loop extrusion and compartmentalisation shape the nuclear genome

Chromatin loops, which bring two distal loci of the same chromosome into close physical proximity, are the ubiquitous units of the three-dimensional genome. Recent advances in understanding the spatial organisation of chromatin suggest that several distinct mechanisms control chromatin interactions, such as loop extrusion by cohesin complexes, compartmentalisation by phase separation, direct protein-protein interactions and others. Here, we review different types of chromatin loops and highlight the factors and processes involved in their regulation. We discuss how loop extrusion and compartmentalisation shape chromatin interactions and how these two processes can either positively or negatively influence each other.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

Nanoscale Interaction of Endonuclease APE1 with DNA

Apurinic/apyrimidinic endonuclease 1 (APE1) is involved in DNA repair and transcriptional regulation mechanisms. This multifunctional activity of APE1 should be supported by specific structural properties of APE1 that have not yet been elucidated. Herein, we applied atomic force microscopy (AFM) to characterize the interactions of APE1 with DNA containing two well-separated G-rich segments. Complexes of APE1 with DNA containing G-rich segments were visualized, and analysis of the complexes revealed the affinity of APE1 to G-rich DNA sequences, and their yield was as high as 53%. Furthermore, APE1 is capable of binding two DNA segments leading to the formation of loops in the DNA-APE1 complexes. The analysis of looped APE1-DNA complexes revealed that APE1 can bridge G-rich segments of DNA. The yield of loops bridging two G-rich DNA segments was 41%. Analysis of protein size in various complexes was performed, and these data showed that loops are formed by APE1 monomer, suggesting that APE1 has two DNA binding sites. The data led us to a model for the interaction of APE1 with DNA and the search for the specific sites. The implication of these new APE1 properties in organizing DNA, by bringing two distant sites together, for facilitating the scanning for damage and coordinating repair and transcription is discussed.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Dissection of a CTCF topological boundary uncovers principles of enhancer-oncogene regulation

Enhancer-gene communication is dependent on topologically associating domains (TADs) and boundaries enforced by the CCCTC-binding factor (CTCF) insulator, but the underlying structures and mechanisms remain controversial. Here, we investigate a boundary that typically insulates fibroblast growth factor (FGF) oncogenes but is disrupted by DNA hypermethylation in gastrointestinal stromal tumors (GISTs). The boundary contains an array of CTCF sites that enforce adjacent TADs, one containing FGF genes and the other containing ANO1 and its putative enhancers, which are specifically active in GIST and its likely cell of origin. We show that coordinate disruption of four CTCF motifs in the boundary fuses the adjacent TADs, allows the ANO1 enhancer to contact FGF3, and causes its robust induction. High-resolution micro-C maps reveal specific contact between transcription initiation sites in the ANO1 enhancer and FGF3 promoter that quantitatively scales with FGF3 induction such that modest changes in contact frequency result in strong changes in expression, consistent with a causal relationship.

Cohesin-Dependent Loop Extrusion: Molecular Mechanics and Role in Cell Physiology

The most prominent representatives of multisubunit SMC complexes, cohesin and condensin, are best known as structural components of mitotic chromosomes. It turned out that these complexes, as well as their bacterial homologues, are molecular motors, the ATP-dependent movement of these complexes along DNA threads leads to the formation of DNA loops. In recent years, we have witnessed an avalanche-like accumulation of data on the process of SMC dependent DNA looping, also known as loop extrusion. This review briefly summarizes the current understanding of the place and role of cohesin-dependent extrusion in cell physiology and presents a number of models describing the potential molecular mechanism of extrusion in a most compelling way. We conclude the review with a discussion of how the capacity of cohesin to extrude DNA loops may be mechanistically linked to its involvement in sister chromatid cohesion.

The N-terminal dimerization domains of human and Drosophila CTCF have similar functionality

BACKGROUND

CTCF is highly likely to be the ancestor of proteins that contain large clusters of C2H2 zinc finger domains, and its conservation is observed across most bilaterian organisms. In mammals, CTCF is the primary architectural protein involved in organizing chromosome topology and mediating enhancer-promoter interactions over long distances. In Drosophila, CTCF (dCTCF) cooperates with other architectural proteins to establish long-range interactions and chromatin boundaries. CTCFs of various organisms contain an unstructured N-terminal dimerization domain (DD) and clusters comprising eleven zinc-finger domains of the C2H2 type. The Drosophila (dCTCF) and human (hCTCF) CTCFs share sequence homology in only five C2H2 domains that specifically bind to a conserved 15 bp motif.

RESULTS

Previously, we demonstrated that CTCFs from different organisms carry unstructured N-terminal dimerization domains (DDs) that lack sequence homology. Here we used the CTCFattP(mCh) platform to introduce desired changes in the Drosophila CTCF gene and generated a series of transgenic lines expressing dCTCF with different variants of the N-terminal domain. Our findings revealed that the functionality of dCTCF is significantly affected by the deletion of the N-terminal DD. Additionally, we observed a strong impact on the binding of the dCTCF mutant to chromatin upon deletion of the DD. However, chromatin binding was restored in transgenic flies expressing a chimeric CTCF protein with the DD of hCTCF. Although the chimeric protein exhibited lower expression levels than those of the dCTCF variants, it efficiently bound to chromatin similarly to the wild type (wt) protein.

CONCLUSIONS

Our findings suggest that one of the evolutionarily conserved functions of the unstructured N-terminal dimerization domain is to recruit dCTCF to its genomic sites in vivo.

Identification of DNA methylation episignature for the intellectual developmental disorder, autosomal dominant 21 syndrome, caused by variants in the CTCF gene

PURPOSE

The main objective of this study was to assess clinical features and genome-wide DNA methylation profiles in individuals affected by intellectual developmental disorder, autosomal dominant 21 (IDD21) syndrome, caused by variants in the CCCTC-binding factor (CTCF) gene.

METHODS

DNA samples were extracted from peripheral blood of 16 individuals with clinical features and genetic findings consistent with IDD21. DNA methylation analysis was performed using the Illumina Infinium Methylation EPIC Bead Chip microarrays. The methylation levels were fitted in a multivariate linear regression model to identify the differentially methylated probes. A binary support vector machine classification model was constructed to differentiate IDD21 samples from controls.

RESULTS

We identified a highly specific, reproducible, and sensitive episignature associated with CTCF variants. Six variants of uncertain significance were tested, of which 2 mapped to the IDD21 episignature and clustered alongside IDD21 cases in both heatmap and multidimensional scaling plots. Comparison of the genomic DNA methylation profile of IDD21 with that of 56 other neurodevelopmental disorders provided insights into the underlying molecular pathophysiology of this disorder.

CONCLUSION

The robust and specific CTCF/IDD21 episignature expands the growing list of neurodevelopmental disorders with distinct DNA methylation profiles, which can be applied as supporting evidence in variant classification.

SpecLoop predicts cell type-specific chromatin loop via transcription factor cooperation

Cell-type-Specific Chromatin Loops (CSCLs) are crucial for gene regulation and cell fate determination. However, the mechanisms governing their establishment remain elusive. Here, we present SpecLoop, a network regularization-based machine learning framework, to investigate the role of transcription factors (TFs) cooperation in CSCL formation. SpecLoop integrates multi-omics data, including gene expression, chromatin accessibility, sequence, protein-protein interaction, and TF binding motif data, to predict CSCLs and identify TF cooperations. Using high resolution Hi-C data as the gold standard, SpecLoop accurately predicts CSCL in GM12878, IMR90, HeLa-S3, K562, HUVEC, HMEC, and NHEK seven cell types, with the AUROC values ranging from 0.8645 to 0.9852 and AUPR values ranging from 0.8654 to 0.9734. Notably SpecLoop demonstrates improved accuracy in predicting long-distance CSCLs and identifies TF complexes with strong predictive ability. Our study systematically explores the TFs and TF pairs associated with CSCL through effective integration of diverse omics data. SpecLoop is freely available at https://github.com/AMSSwanglab/SpecLoop.

BORIS/CTCFL epigenetically reprograms clustered CTCF binding sites into alternative transcriptional start sites

BACKGROUND

Pervasive usage of alternative promoters leads to the deregulation of gene expression in carcinogenesis and may drive the emergence of new genes in spermatogenesis. However, little is known regarding the mechanisms underpinning the activation of alternative promoters.

RESULTS

Here we describe how alternative cancer-testis-specific transcription is activated. We show that intergenic and intronic CTCF binding sites, which are transcriptionally inert in normal somatic cells, could be epigenetically reprogrammed into active de novo promoters in germ and cancer cells. BORIS/CTCFL, the testis-specific paralog of the ubiquitously expressed CTCF, triggers the epigenetic reprogramming of CTCF sites into units of active transcription. BORIS binding initiates the recruitment of the chromatin remodeling factor, SRCAP, followed by the replacement of H2A histone with H2A.Z, resulting in a more relaxed chromatin state in the nucleosomes flanking the CTCF binding sites. The relaxation of chromatin around CTCF binding sites facilitates the recruitment of multiple additional transcription factors, thereby activating transcription from a given binding site. We demonstrate that the epigenetically reprogrammed CTCF binding sites can drive the expression of cancer-testis genes, long noncoding RNAs, retro-pseudogenes, and dormant transposable elements.

CONCLUSIONS

Thus, BORIS functions as a transcription factor that epigenetically reprograms clustered CTCF binding sites into transcriptional start sites, promoting transcription from alternative promoters in both germ cells and cancer cells.

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 Comprehensive Evaluation of Generalizability of Deep Learning-Based Hi-C Resolution Improvement Methods

Hi-C is a widely used technique to study the 3D organization of the genome. Due to its high sequencing cost, most of the generated datasets are of a coarse resolution, which makes it impractical to study finer chromatin features such as Topologically Associating Domains (TADs) and chromatin loops. Multiple deep learning-based methods have recently been proposed to increase the resolution of these datasets by imputing Hi-C reads (typically called upscaling). However, the existing works evaluate these methods on either synthetically downsampled datasets, or a small subset of experimentally generated sparse Hi-C datasets, making it hard to establish their generalizability in the real-world use case. We present our framework-Hi-CY-that compares existing Hi-C resolution upscaling methods on seven experimentally generated low-resolution Hi-C datasets belonging to various levels of read sparsities originating from three cell lines on a comprehensive set of evaluation metrics. Hi-CY also includes four downstream analysis tasks, such as TAD and chromatin loops recall, to provide a thorough report on the generalizability of these methods. We observe that existing deep learning methods fail to generalize to experimentally generated sparse Hi-C datasets, showing a performance reduction of up to 57%. As a potential solution, we find that retraining deep learning-based methods with experimentally generated Hi-C datasets improves performance by up to 31%. More importantly, Hi-CY shows that even with retraining, the existing deep learning-based methods struggle to recover biological features such as chromatin loops and TADs when provided with sparse Hi-C datasets. Our study, through the Hi-CY framework, highlights the need for rigorous evaluation in the future. We identify specific avenues for improvements in the current deep learning-based Hi-C upscaling methods, including but not limited to using experimentally generated datasets for training.

TET1 regulates responses to house dust mite by altering chromatin accessibility, DNA methylation, and gene expression in airway epithelial cells

Background

Previous studies have identified TET1 as a potential key regulator of genes linked to asthma. TET1 has been shown to transcriptionally respond to house dust mite extract, an allergen known to directly cause allergic asthma development, and regulate the expression of genes involved in asthma. How TET1 regulates expression of these genes, however, is unknown. TET1 is a DNA demethylase; therefore, most prior research on TET1-based gene regulation has focused on how TET1 affects methylation. However, TET1 can also interact directly with transcription factors and histone modifiers to regulate gene expression. Understanding how TET1 regulates expression to contribute to allergic responses and asthma development thus requires a comprehensive approach. To this end, we measured mRNA expression, DNA methylation, chromatin accessibility and histone modifications in control and TET1 knockdown human bronchial epithelial cells treated or untreated with house dust mite extract.

Results

Throughout our analyses, we detected strong similarities between the effects of TET1 knockdown alone and the effects of HDM treatment alone. One especially striking pattern was that both TET1 knockdown and HDM treatment generally led to decreased chromatin accessibility at largely the same genomic loci. Transcription factor enrichment analyses indicated that altered chromatin accessibility following the loss of TET1 may affect, or be affected by, CTCF and CEBP binding. TET1 loss also led to changes in DNA methylation, but these changes were generally in regions where accessibility was not changing.

Conclusions

TET1 regulates gene expression through different mechanisms (DNA methylation and chromatin accessibility) in different parts of the genome in the airway epithelial cells, which mediates inflammatory responses to allergen. Collectively, our data suggest novel molecular mechanisms through which TET1 regulates critical pathways following allergen challenges and contributes to the development of asthma.

The PENGUIN approach to reconstruct protein interactions at enhancer-promoter regions and its application to prostate cancer

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

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.

CTCF/cohesin organize the ground state of chromatin-nuclear speckle association

The interchromatin space in the cell nucleus contains various membrane-less nuclear bodies. Recent findings indicate that nuclear speckles, comprising a distinct nuclear body, exhibit interactions with certain chromatin regions in a ground state. Key questions are how this ground state of chromatin-nuclear speckle association is established and what are the gene regulatory roles of this layer of nuclear organization. We report here that chromatin structural factors CTCF and cohesin are required for full ground state association between DNA and nuclear speckles. Disruption of ground state DNA-speckle contacts via either CTCF depletion or cohesin depletion had minor effects on basal level expression of speckle-associated genes, however we show strong negative effects on stimulus-dependent induction of speckle-associated genes. We identified a putative speckle targeting motif (STM) within cohesin subunit RAD21 and demonstrated that the STM is required for chromatin-nuclear speckle association. In contrast to reduction of CTCF or RAD21, depletion of the cohesin releasing factor WAPL stabilized cohesin on chromatin and DNA-speckle contacts, resulting in enhanced inducibility of speckle-associated genes. In addition, we observed disruption of chromatin-nuclear speckle association in patient derived cells with Cornelia de Lange syndrome (CdLS), a congenital neurodevelopmental diagnosis involving defective cohesin pathways, thus revealing nuclear speckles as an avenue for therapeutic inquiry. In summary, our findings reveal a mechanism to establish the ground organizational state of chromatin-speckle association, to promote gene inducibility, and with relevance to human disease.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

Genome control by SMC complexes

Many cellular processes require large-scale rearrangements of chromatin structure. Structural maintenance of chromosomes (SMC) protein complexes are molecular machines that can provide structure to chromatin. These complexes can connect DNA elements in cis, walk along DNA, build and processively enlarge DNA loops and connect DNA molecules in trans to hold together the sister chromatids. These DNA-shaping abilities place SMC complexes at the heart of many DNA-based processes, including chromosome segregation in mitosis, transcription control and DNA replication, repair and recombination. In this Review, we discuss the latest insights into how SMC complexes such as cohesin, condensin and the SMC5-SMC6 complex shape DNA to direct these fundamental chromosomal processes. We also consider how SMC complexes, by building chromatin loops, can counteract the natural tendency of alike chromatin regions to cluster. SMC complexes thus control nuclear organization by participating in a molecular tug of war that determines the architecture of our genome.

CTCF and R-loops are boundaries of cohesin-mediated DNA looping

Cohesin and CCCTC-binding factor (CTCF) are key regulatory proteins of three-dimensional (3D) genome organization. Cohesin extrudes DNA loops that are anchored by CTCF in a polar orientation. Here, we present direct evidence that CTCF binding polarity controls cohesin-mediated DNA looping. Using single-molecule imaging, we demonstrate that a critical N-terminal motif of CTCF blocks cohesin translocation and DNA looping. The cryo-EM structure of the cohesin-CTCF complex reveals that this CTCF motif ahead of zinc fingers can only reach its binding site on the STAG1 cohesin subunit when the N terminus of CTCF faces cohesin. Remarkably, a C-terminally oriented CTCF accelerates DNA compaction by cohesin. DNA-bound Cas9 and Cas12a ribonucleoproteins are also polar cohesin barriers, indicating that stalling may be intrinsic to cohesin itself. Finally, we show that RNA-DNA hybrids (R-loops) block cohesin-mediated DNA compaction in vitro and are enriched with cohesin subunits in vivo, likely forming TAD boundaries.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Cobind: quantitative analysis of the genomic overlaps

Motivation

Analyzing the overlap between two sets of genomic intervals is a frequent task in the field of bioinformatics. Typically, this is accomplished by counting the number (or proportion) of overlapped regions, which applies an arbitrary threshold to determine if two genomic intervals are overlapped. By making binary calls but disregarding the magnitude of the overlap, such an approach often leads to biased, non-reproducible, and incomparable results.

Results

We developed the cobind package, which incorporates six statistical measures: the Jaccard coefficient, Sørensen-Dice coefficient, Szymkiewicz-Simpson coefficient, collocation coefficient, pointwise mutual information (PMI), and normalized PMI. These measures allow for a quantitative assessment of the collocation strength between two sets of genomic intervals. To demonstrate the effectiveness of these methods, we applied them to analyze CTCF's binding sites identified from ChIP-seq, cancer-specific open-chromatin regions (OCRs) identified from ATAC-seq of 17 cancer types, and oligodendrocytes-specific OCRs identified from scATAC-seq. Our results indicated that these new approaches effectively re-discover CTCF's cofactors, as well as cancer-specific and oligodendrocytes-specific master regulators implicated in disease and cell type development.

Availability and implementation

The cobind package is implemented in Python and freely available at https://cobind.readthedocs.io/en/latest/.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

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.

PyMEGABASE: Predicting cell-type-specific structural annotations of chromosomes using the epigenome

The folding patterns of interphase genomes in higher eukaryotes, as obtained from DNA-proximity-ligation or Hi-C experiments, are used to classify loci into structural classes called compartments and subcompartments. These structurally annotated (sub)compartments are known to exhibit specific epigenomic characteristics and cell-type-specific variations. To explore the relationship between genome structure and the epigenome, we present PyMEGABASE (PYMB), a maximum-entropy-based neural network model that predicts (sub)compartment annotations of a locus based solely on the local epigenome, such as ChIP-Seq of histone post-translational modifications. PYMB builds upon our previous model while improving robustness, capability to handle diverse inputs and user-friendly implementation. We employed PYMB to predict subcompartments for over a hundred human cell types available in ENCODE, shedding light on the links between subcompartments, cell identity, and epigenomic signals. The fact that PYMB, trained on data for human cells, can accurately predict compartments in mice suggests that the model is learning underlying physicochemical principles transferable across cell types and species. Reliable at higher resolutions (up to 5 kbp), PYMB is used to investigate compartment-specific gene expression. Not only can PYMB generate (sub)compartment information without Hi-C experiments, but its predictions are also interpretable. Analyzing PYMB's trained parameters, we explore the importance of various epigenomic marks in each subcompartment prediction. Furthermore, the predictions of the model can be used as input for OpenMiChroM software, which has been calibrated to generate three-dimensional structures of the genome. Detailed documentation of PYMB is available at https://pymegabase.readthedocs.io, including an installation guide using pip or conda, and Jupyter/Colab notebook tutorials.

Single-nucleus multiregion transcriptomic analysis of brain vasculature in Alzheimer's disease

Cerebrovascular dysregulation is a hallmark of Alzheimer's disease (AD), but the changes that occur in specific cell types have not been fully characterized. Here, we profile single-nucleus transcriptomes in the human cerebrovasculature in six brain regions from 220 individuals with AD and 208 age-matched controls. We annotate 22,514 cerebrovascular cells, including 11 subtypes of endothelial, pericyte, smooth muscle, perivascular fibroblast and ependymal cells. We identify 2,676 differentially expressed genes in AD, including downregulation of PDGFRB in pericytes, and of ABCB1 and ATP10A in endothelial cells, and validate the downregulation of SLC6A1 and upregulation of APOD, INSR and COL4A1 in postmortem AD brain tissues. We detect vasculature, glial and neuronal coexpressed gene modules, suggesting coordinated neurovascular unit dysregulation in AD. Integration with AD genetics reveals 125 AD differentially expressed genes directly linked to AD-associated genetic variants. Lastly, we show that APOE4 genotype-associated differences are significantly enriched among AD-associated genes in capillary and venule endothelial cells, as well as subsets of pericytes and fibroblasts.

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.