CTCF: master weaver of the genome

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Enhanced performance of gene expression predictive models with protein-mediated spatial chromatin interactions

There have been multiple attempts to predict the expression of the genes based on the sequence, epigenetics, and various other factors. To improve those predictions, we have decided to investigate adding protein-specific 3D interactions that play a significant role in the condensation of the chromatin structure in the cell nucleus. To achieve this, we have used the architecture of one of the state-of-the-art algorithms, ExPecto, and investigated the changes in the model metrics upon adding the spatially relevant data. We have used ChIA-PET interactions that are mediated by cohesin (24 cell lines), CTCF (4 cell lines), and RNAPOL2 (4 cell lines). As the output of the study, we have developed the Spatial Gene Expression (SpEx) algorithm that shows statistically significant improvements in most cell lines. We have compared ourselves to the baseline ExPecto model, which obtained a 0.82 Spearman's rank correlation coefficient (SCC) score, and 0.85, which is reported by newer Enformer were able to obtain the average correlation score of 0.83. However, in some cases (e.g. RNAPOL2 on GM12878), our improvement reached 0.04, and in some cases (e.g. RNAPOL2 on H1), we reached an SCC of 0.86.

DNA strand asymmetry generated by CpG hemimethylation has opposing effects on CTCF binding

CpG methylation generally occurs on both DNA strands and is essential for mammalian development and differentiation. Until recently, hemimethylation, in which only one strand is methylated, was considered to be simply a transitory state generated during DNA synthesis. The discovery that a subset of CCCTC-binding factor (CTCF) binding sites is heritably hemimethylated suggests that hemimethylation might have an unknown biological function. Here we show that the binding of CTCF is profoundly altered by which DNA strand is methylated and by the specific CTCF binding motif. CpG methylation on the motif strand can inhibit CTCF binding by up to 7-fold, whereas methylation on the opposite strand can stimulate binding by up to 4-fold. Thus, hemimethylation can alter binding by up to 28-fold in a strand-specific manner. The mechanism for sensing methylation on the opposite strand requires two critical residues, V454 and S364, within CTCF zinc fingers 7 and 4. Similar to methylation, CpG hydroxymethylation on the motif strand can inhibit CTCF binding by up to 4-fold. However, hydroxymethylation on the opposite strand removes the stimulatory effect. Strand-specific methylation states may therefore provide a mechanism to explain the transient and dynamic nature of CTCF-mediated chromatin interactions.

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.

Single-cell multi-omics analysis reveals dysfunctional Wnt signaling of spermatogonia in non-obstructive azoospermia

Background

Non-obstructive azoospermia (NOA) is the most severe type that leads to 1% of male infertility. Wnt signaling governs normal sperm maturation. However, the role of Wnt signaling in spermatogonia in NOA has incompletely been uncovered, and upstream molecules regulating Wnt signaling remain unclear.

Methods

Bulk RNA sequencing (RNA-seq) of NOA was used to identify the hub gene module in NOA utilizing weighted gene co-expression network analyses (WGCNAs). Single-cell RNA sequencing (scRNA-seq) of NOA was employed to explore dysfunctional signaling pathways in the specific cell type with gene sets of signaling pathways. Single-cell regulatory network inference and clustering (pySCENIC) for Python analysis was applied to speculate putative transcription factors in spermatogonia. Moreover, single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) determined the regulated genes of these transcription factors. Finally, spatial transcriptomic data were used to analyze cell type and Wnt signaling spatial distribution.

Results

The Wnt signaling pathway was demonstrated to be enriched in the hub gene module of NOA by bulk RNA-seq. Then, scRNA-seq data revealed the downregulated activity and dysfunction of Wnt signaling of spermatogonia in NOA samples. Conjoint analyses of the pySCENIC algorithm and scATAC-seq data indicated that three transcription factors (CTCF, AR, and ARNTL) were related to the activities of Wnt signaling in NOA. Eventually, spatial expression localization of Wnt signaling was identified to be in accordance with the distribution patterns of spermatogonia, Sertoli cells, and Leydig cells.

Conclusion

In conclusion, we identified that downregulated Wnt signaling of spermatogonia in NOA and three transcription factors (CTCF, AR, and ARNTL) may be involved in this dysfunctional Wnt signaling. These findings provide new mechanisms for NOA and new therapeutic targets for NOA patients.

A complex network of transcription factors and epigenetic regulators involved in bovine leukemia virus transcriptional regulation

Bovine Leukemia Virus (BLV) is the etiological agent of enzootic bovine leukosis, a disease characterized by the neoplastic proliferation of B cells in cattle. While most European countries have introduced efficient eradication programs, BLV is still present worldwide and no treatment is available. A major feature of BLV infection is the viral latency, which enables the escape from the host immune system, the maintenance of a persistent infection and ultimately the tumoral development. BLV latency is a multifactorial phenomenon resulting in the silencing of viral genes due to genetic and epigenetic repressions of the viral promoter located in the 5' Long Terminal Repeat (5'LTR). However, viral miRNAs and antisense transcripts are expressed from two different proviral regions, respectively the miRNA cluster and the 3'LTR. These latter transcripts are expressed despite the viral latency affecting the 5'LTR and are increasingly considered to take part in tumoral development. In the present review, we provide a summary of the experimental evidence that has enabled to characterize the molecular mechanisms regulating each of the three BLV transcriptional units, either through cis-regulatory elements or through epigenetic modifications. Additionally, we describe the recently identified BLV miRNAs and antisense transcripts and their implications in BLV-induced tumorigenesis. Finally, we discuss the relevance of BLV as an experimental model for the closely related human T-lymphotropic virus HTLV-1.

Single-nucleus multi-omics of human stem cell-derived islets identifies deficiencies in lineage specification

Insulin-producing β cells created from human pluripotent stem cells have potential as a therapy for insulin-dependent diabetes, but human pluripotent stem cell-derived islets (SC-islets) still differ from their in vivo counterparts. To better understand the state of cell types within SC-islets and identify lineage specification deficiencies, we used single-nucleus multi-omic sequencing to analyse chromatin accessibility and transcriptional profiles of SC-islets and primary human islets. Here we provide an analysis that enabled the derivation of gene lists and activity for identifying each SC-islet cell type compared with primary islets. Within SC-islets, we found that the difference between β cells and awry enterochromaffin-like cells is a gradient of cell states rather than a stark difference in identity. Furthermore, transplantation of SC-islets in vivo improved cellular identities overtime, while long-term in vitro culture did not. Collectively, our results highlight the importance of chromatin and transcriptional landscapes during islet cell specification and maturation.

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.

Association of mutation and expression of the brother of the regulator of imprinted sites (BORIS) gene with breast cancer progression

INTRODUCTION

The BORIS, 11 zinc-finger transcription factors, is a member of the cancer-testis antigen (CTA) family. It is mapped to chromosome number 20q13.2 and this region is genetically linked to the early onset of breast cancer. The current study analyzed the correlation between BORIS mutations and the expression of the protein in breast cancer cases.

MATERIALS AND METHODS

A population-based study including a total of 155 breast cancer tissue samples and an equal number of normal adjacent tissues from Indian female breast cancer patients was carried out. Mutations of the BORIS gene were detected by polymerase chain reaction-single standard confirmation polymorphisms (PCR-SSCP) and automated DNA sequencing and by immunohistochemistry for BORIS protein expression were performed. The observed findings were correlated with several clinicopathological parameters to find out the clinical relevance of associations.

RESULTS

Of all the cases 16.12% (25/155) showed mutations in the BORIS gene. The observed mutations present on codon 329 are missense, leading to Val> Ile (G>A) change on exon 5 of the BORIS gene. A significant association was observed between mutations of the BORIS gene and some clinicopathological features like nodal status (p = 0.013), estrogen receptor (ER) expression (p = 0.008), progesterone receptor (PR) expression (p = 0.039), clinical stage (p = 0.010) and menopausal status (p = 0.023). The protein expression analysis showed 20.64% (32/155) samples showing low or no expression (+), 34.19% (53/155) with moderate expression (++), and 45.17% (70/155) showing high expression (+++) of BORIS protein. A significant association was observed between the expression of BORIS protein and clinicopathological features like clinical stage (p = 0.013), nodal status (p = 0.049), ER expression (p = 0.039), and PR expression (p = 0.027). When mutation and protein expression were correlated in combination with clinicopathological parameters a significant association was observed in the category of high (+++) level of BORIS protein expression (p = 0.017).

CONCLUSION

The BORIS mutations and high protein expression occur frequently in carcinoma of the breast suggesting their association with the onset and progression of breast carcinoma. Further, the BORIS has the potential to be used as a biomarker.

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.

Structure and Functions of HMGB2 Protein

High-Mobility Group (HMG) chromosomal proteins are the most numerous nuclear non-histone proteins. HMGB domain proteins are the most abundant and well-studied HMG proteins. They are involved in variety of biological processes. HMGB1 and HMGB2 were the first members of HMGB-family to be discovered and are found in all studied eukaryotes. Despite the high degree of homology, HMGB1 and HMGB2 proteins differ from each other both in structure and functions. In contrast to HMGB2, there is a large pool of works devoted to the HMGB1 protein whose structure-function properties have been described in detail in our previous review in 2020. In this review, we attempted to bring together diverse data about the structure and functions of the HMGB2 protein. The review also describes post-translational modifications of the HMGB2 protein and its role in the development of a number of diseases. Particular attention is paid to its interaction with various targets, including DNA and protein partners. The influence of the level of HMGB2 expression on various processes associated with cell differentiation and aging and its ability to mediate the differentiation of embryonic and adult stem cells are also discussed.

Loop stacking organizes genome folding from TADs to chromosomes

Although population-level analyses revealed significant roles for CTCF and cohesin in mammalian genome organization, their contributions at the single-cell level remain incompletely understood. Here, we used a super-resolution microscopy approach to measure the effects of removal of CTCF or cohesin in mouse embryonic stem cells. Single-chromosome traces revealed cohesin-dependent loops, frequently stacked at their loop anchors forming multi-way contacts (hubs), bridging across TAD boundaries. Despite these bridging interactions, chromatin in intervening TADs was not intermixed, remaining separated in distinct loops around the hub. At the multi-TAD scale, steric effects from loop stacking insulated local chromatin from ultra-long range (>4 Mb) contacts. Upon cohesin removal, the chromosomes were more disordered and increased cell-cell variability in gene expression. Our data revise the TAD-centric understanding of CTCF and cohesin and provide a multi-scale, structural picture of how they organize the genome on the single-cell level through distinct contributions to loop stacking.

Parallel sequencing of extrachromosomal circular DNAs and transcriptomes in single cancer cells

Extrachromosomal DNAs (ecDNAs) are common in cancer, but many questions about their origin, structural dynamics and impact on intratumor heterogeneity are still unresolved. Here we describe single-cell extrachromosomal circular DNA and transcriptome sequencing (scEC&T-seq), a method for parallel sequencing of circular DNAs and full-length mRNA from single cells. By applying scEC&T-seq to cancer cells, we describe intercellular differences in ecDNA content while investigating their structural heterogeneity and transcriptional impact. Oncogene-containing ecDNAs were clonally present in cancer cells and drove intercellular oncogene expression differences. In contrast, other small circular DNAs were exclusive to individual cells, indicating differences in their selection and propagation. Intercellular differences in ecDNA structure pointed to circular recombination as a mechanism of ecDNA evolution. These results demonstrate scEC&T-seq as an approach to systematically characterize both small and large circular DNA in cancer cells, which will facilitate the analysis of these DNA elements in cancer and beyond.

Landscape of chromatin remodeling gene alterations in endometrial carcinoma

OBJECTIVE

Chromatin remodeling genes (CRGs) encode components of epigenetic regulatory mechanisms and alterations in these genes have been identified in several tumor types, including gynecologic cancers. In this study, we sought to investigate the prevalence and clinicopathological associations of CRG alterations in endometrial carcinoma (EC).

METHODS

We performed a retrospective analysis of 660 ECs sequenced using a clinical massively parallel sequencing assay targeting up to 468 genes, including 25 CRGs, and defined the presence of somatic CRG alterations. Clinicopathologic features were obtained for all cases. Immunohistochemical interrogation of ARID1A and PTEN proteins was performed in a subset of samples.

RESULTS

Of the 660 ECs sequenced, 438 (66.4%) harbored CRG alterations covered by our panel. The most commonly altered CRG was ARID1A (46%), followed by CTCF (21%), KMT2D (18%), KMT2B (17%), BCOR (16%), ARID1B (12%) and SMARCA4 (11%). We found that ARID1A genetic alterations were preferentially bi-allelic and often corresponded to altered ARID1A protein expression in ECs. We further observed that ARID1A alterations were often subclonal when compared to PTEN alterations, which were primarily clonal in ECs harboring both mutations. Finally, CRG alterations were associated with an increased likelihood of myometrial and lymphovascular invasion in endometrioid ECs.

CONCLUSION

CRG alterations are common in EC and are associated with clinicopathologic features and likely play a crucial role in EC.

Three-Dimensional Chromatin Structure of the EBV Genome: A Crucial Factor in Viral Infection

Epstein-Barr Virus (EBV) is a human gamma-herpesvirus that is widespread worldwide. To this day, about 200,000 cancer cases per year are attributed to EBV infection. EBV is capable of infecting both B cells and epithelial cells. Upon entry, viral DNA reaches the nucleus and undergoes a process of circularization and chromatinization and establishes a latent lifelong infection in host cells. There are different types of latency all characterized by different expressions of latent viral genes correlated with a different three-dimensional architecture of the viral genome. There are multiple factors involved in the regulation and maintenance of this three-dimensional organization, such as CTCF, PARP1, MYC and Nuclear Lamina, emphasizing its central role in latency maintenance.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

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.

SATB1 regulates 3D genome architecture in T cells by constraining chromatin interactions surrounding CTCF-binding sites

Special AT-rich sequence binding protein 1 (SATB1) has long been proposed to act as a global chromatin loop organizer in T cells. However, the exact functions of SATB1 in spatial genome organization remain elusive. Here we show that the depletion of SATB1 in human and murine T cells leads to transcriptional dysregulation for genes involved in T cell activation, as well as alterations of 3D genome architecture at multiple levels, including compartments, topologically associating domains, and loops. Importantly, SATB1 extensively colocalizes with CTCF throughout the genome. Depletion of SATB1 leads to increased chromatin contacts among and across the SATB1/CTCF co-occupied sites, thereby affecting the transcription of critical regulators of T cell activation. The loss of SATB1 does not affect CTCF occupancy but significantly reduces the retention of CTCF in the nuclear matrix. Collectively, our data show that SATB1 contributes to 3D genome organization by constraining chromatin topology surrounding CTCF-binding sites.

BORIS/CTCFL-mediated chromatin accessibility alterations promote a pro-invasive transcriptional signature in melanoma cells

Melanoma is the deadliest form of skin cancer, due to its tendency to metastasize early. Brother of regulator of imprinted sites (BORIS), also known as CCCTC binding factor-like (CTCFL), is a transcription regulator that becomes ectopically expressed in melanoma. We recently showed that BORIS contributes to melanoma phenotype switching by altering the gene expression program of melanoma cells from an intermediate melanocytic state toward a more mesenchymal-like state. However, the mechanism underlying this transcriptional switch remains unclear. Here, ATAC-seq was used to study BORIS-mediated chromatin accessibility alterations in melanoma cells harboring an intermediate melanocytic state. The gene set that gained promoter accessibility, following ectopic BORIS expression, showed enrichment for biological processes associated with melanoma invasion, while promoters of genes associated with proliferation showed reduced accessibility. Integration of ATAC-seq and RNA-seq data demonstrated that increased chromatin accessibility was associated with transcriptional upregulation of genes involved in tumor progression processes, and the aberrant activation of oncogenic transcription factors, while reduced chromatin accessibility and downregulated genes were associated with repressed activity of tumor suppressors and proliferation factors. Together, these findings indicate that BORIS mediates transcriptional reprogramming in melanoma cells by altering chromatin accessibility and gene expression, shifting the cellular transcription landscape of melanoma cells toward a mesenchymal-like genetic signature.

The long-range interaction between two GNAS imprinting control regions delineates pseudohypoparathyroidism type 1B pathogenesis

Genetic defects of GNAS, the imprinted gene encoding the stimulatory G protein α-subunit, are responsible for multiple diseases. Abnormal GNAS imprinting causes pseudohypoparathyroidism type 1B (PHP1B), a prototype of mammalian end-organ hormone resistance. Hypomethylation at the maternally methylated GNAS A/B region is the only shared defect in patients with PHP1B. In autosomal dominant (AD) PHP1B kindreds, A/B hypomethylation is associated with maternal microdeletions at either the GNAS NESP55 differentially methylated region or the STX16 gene located approximately 170 kb upstream. Functional evidence is meager regarding the causality of these microdeletions. Moreover, the mechanisms linking A/B methylation and the putative imprinting control regions (ICRs) NESP-ICR and STX16-ICR remain unknown. Here, we generated a human embryonic stem cell model of AD-PHP1B by introducing ICR deletions using CRISPR/Cas9. With this model, we showed that the NESP-ICR is required for methylation and transcriptional silencing of A/B on the maternal allele. We also found that the SXT16-ICR is a long-range enhancer of NESP55 transcription, which originates from the maternal NESP-ICR. Furthermore, we demonstrated that the STX16-ICR is an embryonic stage-specific enhancer enabled by the direct binding of pluripotency factors. Our findings uncover an essential GNAS imprinting control mechanism and advance the molecular understanding of PHP1B pathogenesis.

Actively transcribed rDNA and distal junction (DJ) sequence are involved in association of NORs with nucleoli

Although they are organelles without a limiting membrane, nucleoli have an exclusive structure, built upon the rDNA-rich acrocentric short arms of five human chromosomes (nucleolar organizer regions or NORs). This has raised the question: what are the structural features of a chromosome required for its inclusion in a nucleolus? Previous work has suggested that sequences adjacent to the tandemly repeated rDNA repeat units (DJ, distal junction sequence) may be involved, and we have extended such studies by addressing several issues related to the requirements for the association of NORs with nucleoli. We exploited both a set of somatic cell hybrids containing individual human acrocentric chromosomes and a set of Human Artificial Chromosomes (HACs) carrying different parts of a NOR, including an rDNA unit or DJ or PJ (proximal junction) sequence. Association of NORs with nucleoli was increased when constituent rDNA was transcribed and may be also affected by the status of heterochromatin blocks formed next to the rDNA arrays. Furthermore, our data suggest that a relatively small size DJ region, highly conserved in evolution, is also involved, along with the rDNA repeats, in the localization of p-arms of acrocentric chromosomes in nucleoli. Thus, we infer a cooperative action of rDNA sequence-stimulated by its activity-and sequences distal to rDNA contributing to incorporation into nucleoli. Analysis of NOR sequences also identified LncRNA_038958 in the DJ, a candidate transcript with the region of the suggested promoter that is located close to the DJ/rDNA boundary and contains CTCF binding sites. This LncRNA may affect RNA Polymerase I and/or nucleolar activity. Our findings provide the basis for future studies to determine which RNAs and proteins interact critically with NOR sequences to organize the higher-order structure of nucleoli and their function in normal cells and pathological states.

Enhanced performance of gene expression predictive models with protein-mediated spatial chromatin interactions

There have been multiple attempts to predict the expression of the genes based on the sequence, epigenetics, and various other factors. To improve those predictions, we have decided to investigate adding protein-specific 3D interactions that play a major role in the compensation of the chromatin structure in the cell nucleus. To achieve this, we have used the architecture of one of the state-of-the-art algorithms, ExPecto (J. Zhou et al., 2018), and investigated the changes in the model metrics upon adding the spatially relevant data. We have used ChIA-PET interactions that are mediated by cohesin (24 cell lines), CTCF (4 cell lines), and RNAPOL2 (4 cell lines). As the output of the study, we have developed the Spatial Gene Expression (SpEx) algorithm that shows statistically significant improvements in most cell lines.

Prediction of CTCF loop anchor based on machine learning

Introduction: Various activities in biological cells are affected by three-dimensional genome structure. The insulators play an important role in the organization of higher-order structure. CTCF is a representative of mammalian insulators, which can produce barriers to prevent the continuous extrusion of chromatin loop. As a multifunctional protein, CTCF has tens of thousands of binding sites in the genome, but only a portion of them can be used as anchors of chromatin loops. It is still unclear how cells select the anchor in the process of chromatin looping.

Methods: In this paper, a comparative analysis is performed to investigate the sequence preference and binding strength of anchor and non-anchor CTCF binding sites. Furthermore, a machine learning model based on the CTCF binding intensity and DNA sequence is proposed to predict which CTCF sites can form chromatin loop anchors.

Results: The accuracy of the machine learning model that we constructed for predicting the anchor of the chromatin loop mediated by CTCF reached 0.8646. And we find that the formation of loop anchor is mainly influenced by the CTCF binding strength and binding pattern (which can be interpreted as the binding of different zinc fingers).

Discussion: In conclusion, our results suggest that The CTCF core motif and it’s flanking sequence may be responsible for the binding specificity. This work contributes to understanding the mechanism of loop anchor selection and provides a reference for the prediction of CTCF-mediated chromatin loops.

Structural and functional properties of mSWI/SNF chromatin remodeling complexes revealed through single-cell perturbation screens

The mammalian SWI/SNF (mSWI/SNF or BAF) family of chromatin remodeling complexes play critical roles in regulating DNA accessibility and gene expression. The three final-form subcomplexes-cBAF, PBAF, and ncBAF-are distinct in biochemical componentry, chromatin targeting, and roles in disease; however, the contributions of their constituent subunits to gene expression remain incompletely defined. Here, we performed Perturb-seq-based CRISPR-Cas9 knockout screens targeting mSWI/SNF subunits individually and in select combinations, followed by single-cell RNA-seq and SHARE-seq. We uncovered complex-, module-, and subunit-specific contributions to distinct regulatory networks and defined paralog subunit relationships and shifted subcomplex functions upon perturbations. Synergistic, intra-complex genetic interactions between subunits reveal functional redundancy and modularity. Importantly, single-cell subunit perturbation signatures mapped across bulk primary human tumor expression profiles both mirror and predict cBAF loss-of-function status in cancer. Our findings highlight the utility of Perturb-seq to dissect disease-relevant gene regulatory impacts of heterogeneous, multi-component master regulatory complexes.

Associations of CTCF and FOXA1 with androgen and IGF pathways in men with localized prostate cancer

AIMS

To examine associations between the transcription factors CCCTC-binding factor (CTCF) and forkhead box protein A1 (FOXA1) and the androgen receptor (AR) and their association with components of the insulin-like growth factor (IGF)-pathway in a cohort of men with localized prostate cancer.

METHODS

Using prostate tissue samples collected during the Prostate cancer: Evidence of Exercise and Nutrition Trial (PrEvENT) trial (N = 70 to 92, depending on section availability), we assessed the abundance of CTCF, FOXA1, AR, IGFIR, p-mTOR, PTEN and IGFBP-2 proteins using a modified version of the Allred scoring system. Validation studies were performed using large, publicly available datasets (TCGA) (N = 489).

RESULTS

We identified a strong correlation between CTCF and AR staining with benign prostate tissue. CTCF also strongly associated with the IGFIR, with PTEN and with phospho-mTOR. FOXA1 was also correlated with staining for the IGF-IR, with IGFBP-2 and with staining for activated phosphor-mTOR. The staining for the IGF-IR was strongly correlated with the AR.

CONCLUSION

Our findings emphasise the close and complex links between the endocrine controls, well known to play an important role in prostate cancer, and the transcription factors implicated by the recent genetic evidence.

CDK9 inhibition induces epigenetic reprogramming revealing strategies to circumvent resistance in lymphoma

Diffuse large B-cell lymphoma (DLBCL) exhibits significant genetic heterogeneity which contributes to drug resistance, necessitating development of novel therapeutic approaches. Pharmacological inhibitors of cyclin-dependent kinases (CDK) demonstrated pre-clinical activity in DLBCL, however many stalled in clinical development. Here we show that AZD4573, a selective inhibitor of CDK9, restricted growth of DLBCL cells. CDK9 inhibition (CDK9i) resulted in rapid changes in the transcriptome and proteome, with downmodulation of multiple oncoproteins (eg, MYC, Mcl-1, JunB, PIM3) and deregulation of phosphoinotiside-3 kinase (PI3K) and senescence pathways. Following initial transcriptional repression due to RNAPII pausing, we observed transcriptional recovery of several oncogenes, including MYC and PIM3. ATAC-Seq and ChIP-Seq experiments revealed that CDK9i induced epigenetic remodeling with bi-directional changes in chromatin accessibility, suppressed promoter activation and led to sustained reprograming of the super-enhancer landscape. A CRISPR library screen suggested that SE-associated genes in the Mediator complex, as well as AKT1, confer resistance to CDK9i. Consistent with this, sgRNA-mediated knockout of MED12 sensitized cells to CDK9i. Informed by our mechanistic findings, we combined AZD4573 with either PIM kinase or PI3K inhibitors. Both combinations decreased proliferation and induced apoptosis in DLBCL and primary lymphoma cells in vitro as well as resulted in delayed tumor progression and extended survival of mice xenografted with DLBCL in vivo. Thus, CDK9i induces reprogramming of the epigenetic landscape, and super-enhancer driven recovery of select oncogenes may contribute to resistance to CDK9i. PIM and PI3K represent potential targets to circumvent resistance to CDK9i in the heterogeneous landscape of DLBCL.

High-throughput Pore-C reveals the single-allele topology and cell type-specificity of 3D genome folding

Canonical three-dimensional (3D) genome structures represent the ensemble average of pairwise chromatin interactions but not the single-allele topologies in populations of cells. Recently developed Pore-C can capture multiway chromatin contacts that reflect regional topologies of single chromosomes. By carrying out high-throughput Pore-C, we reveal extensive but regionally restricted clusters of single-allele topologies that aggregate into canonical 3D genome structures in two human cell types. We show that fragments in multi-contact reads generally coexist in the same TAD. In contrast, a concurrent significant proportion of multi-contact reads span multiple compartments of the same chromatin type over megabase distances. Synergistic chromatin looping between multiple sites in multi-contact reads is rare compared to pairwise interactions. Interestingly, the single-allele topology clusters are cell type-specific even inside highly conserved TADs in different types of cells. In summary, HiPore-C enables global characterization of single-allele topologies at an unprecedented depth to reveal elusive genome folding principles.

Myoepithelial carcinoma of the parotid gland with a novel CTCF::NCOA2 fusion

We describe a case of a myoepithelial carcinoma of the superficial parotid gland in a 46-year-old male harboring a novel CTCF::NCOA2 gene fusion. To our knowledge, this novel gene fusion has not been described previously in myoepithelial carcinoma. A 46-year-old male patient presented with a mass involving the superficial left parotid gland with extension into the external auditory canal (EAC) and erosion of the conchal cartilage. Histologically, the neoplasm was composed of uniform spindled, epithelioid/ovoid cells arranged in cords and nests within hyalinized to myxoid stroma. On immunohistochemistry (IHC), the tumor cells demonstrated patchy and variable staining for low molecular weight cytokeratin (CAM5.2), pan-cytokeratin (OSCAR), and S-100. Overall, the morphological and immunohistochemical attributes supported a locally aggressive tumor of myoepithelial differentiation consistent with myoepithelial carcinoma. Molecular analysis using a custom 115-gene gene panel by targeted RNA sequencing, showed an in-frame CTCF::NCOA2 fusion. In addition to reporting this novel fusion in myoepithelial carcinoma, we also discuss relevant differential diagnosis, and provide a brief review of NCOA2 gene function in both normal and neoplastic contexts.

Expansion of the genotypic and phenotypic spectrum of CTCF-related disorder guides clinical management: 43 new subjects and a comprehensive literature review

Monoallelic variants of CTCF cause an autosomal dominant neurodevelopmental disorder with a wide range of features, including impacts on the brain, growth, and craniofacial development. A growing number of subjects with CTCF-related disorder (CRD) have been identified due to the increased application of exome sequencing, and further delineation of the clinical spectrum of CRD is needed. Here, we examined the clinical features, including facial profiles, and genotypic spectrum of 107 subjects with identified CTCF variants, including 43 new and 64 previously described subjects. Among the 43 new subjects, 23 novel variants were reported. The cardinal clinical features in subjects with CRD included intellectual disability/developmental delay (91%) with speech delay (65%), motor delay (53%), feeding difficulties/failure to thrive (66%), ocular abnormalities (56%), musculoskeletal anomalies (53%), and behavioral problems (52%). Other congenital anomalies were also reported, but none of them were common. Our findings expanded the genotypic and phenotypic spectrum of CRD that will guide genetic counseling, management, and surveillance care for patients with CRD. Additionally, a newly built facial gestalt on the Face2Gene tool will facilitate prompt recognition of CRD by physicians and shorten a patient's diagnostic odyssey.

Reactivation of Epstein-Barr Virus from Latency Involves Increased RNA Polymerase Activity at CTCF Binding Sites on the Viral Genome

The ability of Epstein-Barr virus (EBV) to switch between latent and lytic infection is key to its long-term persistence, yet the molecular mechanisms behind this switch remain unclear. To investigate transcriptional events during the latent-to-lytic switch, we utilized Precision nuclear Run On followed by deep Sequencing (PRO-Seq) to map cellular RNA polymerase (Pol) activity to single-nucleotide resolution on the host and EBV genome in three different models of EBV latency and reactivation. In latently infected Mutu-I Burkitt lymphoma (BL) cells, Pol activity was enriched at the Qp promoter, the EBER region, and the BHLF1/LF3 transcripts. Upon reactivation with phorbol ester and sodium butyrate, early-phase Pol activity occurred bidirectionally at CTCF sites within the LMP-2A, EBER-1, and RPMS1 loci. PRO-Seq analysis of Akata cells reactivated from latency with anti-IgG and a lymphoblastoid cell line (LCL) reactivated with small molecule C60 showed a similar pattern of early bidirectional transcription initiating around CTCF binding sites, although the specific CTCF sites and viral genes were different for each latency model. The functional importance of CTCF binding, transcription, and reactivation was confirmed using an EBV mutant lacking the LMP-2A CTCF binding site. This virus was unable to reactivate and had disrupted Pol activity at multiple CTCF binding sites relative to the wild-type (WT) virus. Overall, these data suggest that CTCF regulates the viral early transcripts during reactivation from latency. These activities likely help maintain the accessibility of the viral genome to initiate productive replication. IMPORTANCE The ability of EBV to switch between latent and lytic infection is key to its long-term persistence in memory B cells, and its ability to persist in proliferating cells is strongly linked to oncogenesis. During latency, most viral genes are epigenetically silenced, and the virus must overcome this repression to reactivate lytic replication. Reactivation occurs once the immediate early (IE) EBV lytic genes are expressed. However, the molecular mechanisms behind the switch from the latent transcriptional program to begin transcription of the IE genes remain unknown. In this study, we mapped RNA Pol positioning and activity during latency and reactivation. Unexpectedly, Pol activity accumulated at distinct regions characteristic of transcription initiation on the EBV genome previously shown to be associated with CTCF. We propose that CTCF binding at these regions retains Pol to maintain a stable latent chromosome conformation and a rapid response to various reactivation signals.

Reduced methylation correlates with diabetic nephropathy risk in type 1 diabetes

Diabetic nephropathy (DN) is a polygenic disorder with few risk variants showing robust replication in large-scale genome-wide association studies. To understand the role of DNA methylation, it is important to have the prevailing genomic view to distinguish key sequence elements that influence gene expression. This is particularly challenging for DN because genome-wide methylation patterns are poorly defined. While methylation is known to alter gene expression, the importance of this causal relationship is obscured by array-based technologies since coverage outside promoter regions is low. To overcome these challenges, we performed methylation sequencing using leukocytes derived from participants of the Finnish Diabetic Nephropathy (FinnDiane) type 1 diabetes (T1D) study (n = 39) that was subsequently replicated in a larger validation cohort (n = 296). Gene body-related regions made up more than 60% of the methylation differences and emphasized the importance of methylation sequencing. We observed differentially methylated genes associated with DN in 3 independent T1D registries originating from Denmark (n = 445), Hong Kong (n = 107), and Thailand (n = 130). Reduced DNA methylation at CTCF and Pol2B sites was tightly connected with DN pathways that include insulin signaling, lipid metabolism, and fibrosis. To define the pathophysiological significance of these population findings, methylation indices were assessed in human renal cells such as podocytes and proximal convoluted tubule cells. The expression of core genes was associated with reduced methylation, elevated CTCF and Pol2B binding, and the activation of insulin-signaling phosphoproteins in hyperglycemic cells. These experimental observations also closely parallel methylation-mediated regulation in human macrophages and vascular endothelial cells.

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

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

Profiling human pathogenic repeat expansion regions by synergistic and multi-level impacts on molecular connections

Whilst DNA repeat expansions cause numerous heritable human disorders, their origins and underlying pathological mechanisms are often unclear. We collated a dataset comprising 224 human repeat expansions encompassing 203 different genes, and performed a systematic analysis with respect to key topological features at the DNA, RNA and protein levels. Comparison with controls without known pathogenicity and genomic regions lacking repeats, allowed the construction of the first tool to discriminate repeat regions harboring pathogenic repeat expansions (DPREx). At the DNA level, pathogenic repeat expansions exhibited stronger signals for DNA regulatory factors (e.g. H3K4me3, transcription factor-binding sites) in exons, promoters, 5'UTRs and 5'genes but were not significantly different from controls in introns, 3'UTRs and 3'genes. Additionally, pathogenic repeat expansions were also found to be enriched in non-B DNA structures. At the RNA level, pathogenic repeat expansions were characterized by lower free energy for forming RNA secondary structure and were closer to splice sites in introns, exons, promoters and 5'genes than controls. At the protein level, pathogenic repeat expansions exhibited a preference to form coil rather than other types of secondary structure, and tended to encode surface-located protein domains. Guided by these features, DPREx ( http://biomed.nscc-gz.cn/zhaolab/geneprediction/# ) achieved an Area Under the Curve (AUC) value of 0.88 in a test on an independent dataset. Pathogenic repeat expansions are thus located such that they exert a synergistic influence on the gene expression pathway involving inter-molecular connections at the DNA, RNA and protein levels.

Auxin-inducible degron 2 system deciphers functions of CTCF domains in transcriptional regulation

BACKGROUND

CTCF is a well-established chromatin architectural protein that also plays various roles in transcriptional regulation. While CTCF biology has been extensively studied, how the domains of CTCF function to regulate transcription remains unknown. Additionally, the original auxin-inducible degron 1 (AID1) system has limitations in investigating the function of CTCF.

RESULTS

We employ an improved auxin-inducible degron technology, AID2, to facilitate the study of acute depletion of CTCF while overcoming the limitations of the previous AID system. As previously observed through the AID1 system and steady-state RNA analysis, the new AID2 system combined with SLAM-seq confirms that CTCF depletion leads to modest nascent and steady-state transcript changes. A CTCF domain sgRNA library screening identifies the zinc finger (ZF) domain as the region within CTCF with the most functional relevance, including ZFs 1 and 10. Removal of ZFs 1 and 10 reveals genomic regions that independently require these ZFs for DNA binding and transcriptional regulation. Notably, loci regulated by either ZF1 or ZF10 exhibit unique CTCF binding motifs specific to each ZF.

CONCLUSIONS

By extensively comparing the AID1 and AID2 systems for CTCF degradation in SEM cells, we confirm that AID2 degradation is superior for achieving miniAID-tagged protein degradation without the limitations of the AID1 system. The model we create that combines AID2 depletion of CTCF with exogenous overexpression of CTCF mutants allows us to demonstrate how peripheral ZFs intricately orchestrate transcriptional regulation in a cellular context for the first time.

Mechanisms of chromate carcinogenesis by chromatin alterations

In a dynamic environment, organisms must constantly mount an adaptive response to new environmental conditions in order to survive. Novel patterns of gene expression, driven by attendant changes in chromatin architecture, aid in adaptation and survival. Critical mechanisms in the control of gene transcription govern new spatiotemporal chromatin-chromatin interactions that make regulatory DNA elements accessible to the transcription factors that control the response. Consequently, agents that disrupt chromatin structure are likely to have a direct impact on the transcriptional programs of cells and organisms and to drive alterations in fundamental physiological processes. In this regard, hexavalent chromium (Cr(VI)) is of special interest because it interacts directly with cellular proteins, DNA, and other macromolecules, and is likely to upset cell functions that may cause generalized damage to the organism. Here, we will highlight chromium-mediated mechanisms that disrupt chromatin architecture and discuss how these mechanisms are integral to its carcinogenic properties. Emerging evidence indicates that Cr(VI) targets euchromatin, particularly in genomic locations flanking the binding sites of the essential transcription factors CTCF and AP1, and, in so doing, they disrupt nucleosomal architecture. Ultimately, the ensuing changes, if occurring in critical regulatory domains, may establish a new chromatin state, either toxic or adaptive, that will be governed by the corresponding gene transcription changes in key biological processes associated with that state.

Spatiotemporal Epigenetic Control of the Histone Gene Chromatin Landscape during the Cell Cycle

Higher-order genomic organization supports the activation of histone genes in response to cell cycle regulatory cues that epigenetically mediates stringent control of transcription at the G1/S-phase transition. Histone locus bodies (HLBs) are dynamic, non-membranous, phase-separated nuclear domains where the regulatory machinery for histone gene expression is organized and assembled to support spatiotemporal epigenetic control of histone genes. HLBs provide molecular hubs that support synthesis and processing of DNA replication-dependent histone mRNAs. These regulatory microenvironments support long-range genomic interactions among non-contiguous histone genes within a single topologically associating domain (TAD). HLBs respond to activation of the cyclin E/CDK2/NPAT/HINFP pathway at the G1/S transition. HINFP and its coactivator NPAT form a complex within HLBs that controls histone mRNA transcription to support histone protein synthesis and packaging of newly replicated DNA. Loss of HINFP compromises H4 gene expression and chromatin formation, which may result in DNA damage and impede cell cycle progression. HLBs provide a paradigm for higher-order genomic organization of a subnuclear domain that executes an obligatory cell cycle-controlled function in response to cyclin E/CDK2 signaling. Understanding the coordinately and spatiotemporally organized regulatory programs in focally defined nuclear domains provides insight into molecular infrastructure for responsiveness to cell signaling pathways that mediate biological control of growth, differentiation phenotype, and are compromised in cancer.

Role of primary aging hallmarks in Alzheimer´s disease

Alzheimer's disease (AD) is the most common neurodegenerative disease, which severely threatens the health of the elderly and causes significant economic and social burdens. The causes of AD are complex and include heritable but mostly aging-related factors. The primary aging hallmarks include genomic instability, telomere wear, epigenetic changes, and loss of protein stability, which play a dominant role in the aging process. Although AD is closely associated with the aging process, the underlying mechanisms involved in AD pathogenesis have not been well characterized. This review summarizes the available literature about primary aging hallmarks and their roles in AD pathogenesis. By analyzing published literature, we attempted to uncover the possible mechanisms of aberrant epigenetic markers with related enzymes, transcription factors, and loss of proteostasis in AD. In particular, the importance of oxidative stress-induced DNA methylation and DNA methylation-directed histone modifications and proteostasis are highlighted. A molecular network of gene regulatory elements that undergoes a dynamic change with age may underlie age-dependent AD pathogenesis, and can be used as a new drug target to treat AD.

Identification of genomic binding sites and direct target genes for the transcription factor DDIT3/CHOP

DDIT3 is a tightly regulated basic leucine zipper (bZIP) transcription factor and key regulator in cellular stress responses. It is involved in a variety of pathological conditions and may cause cell cycle block and apoptosis. It is also implicated in differentiation of some specialized cell types and as an oncogene in several types of cancer. DDIT3 was originally believed to act as a dominant-negative inhibitor by forming heterodimers with other bZIP transcription factors, preventing their DNA binding and transactivating functions. DDIT3 has, however, been reported to bind DNA and regulate target genes. Here, we employed ChIP sequencing combined with microarray-based expression analysis to identify direct binding motifs and target genes of DDIT3. The results reveal DDIT3 binding to motifs similar to other bZIP transcription factors, known to form heterodimers with DDIT3. Binding to a class III satellite DNA repeat sequence was also detected. DDIT3 acted as a DNA-binding transcription factor and bound mainly to the promotor region of regulated genes. ChIP sequencing analysis of histone H3K27 methylation and acetylation showed a strong overlap between H3K27-acetylated marks and DDIT3 binding. These results support a role for DDIT3 as a transcriptional regulator of H3K27ac-marked genes in transcriptionally active chromatin.

Diverse silent chromatin states modulate genome compartmentalization and loop extrusion barriers

The relationships between chromosomal compartmentalization, chromatin state and function are poorly understood. Here by profiling long-range contact frequencies in HCT116 colon cancer cells, we distinguish three silent chromatin states, comprising two types of heterochromatin and a state enriched for H3K9me2 and H2A.Z that exhibits neutral three-dimensional interaction preferences and which, to our knowledge, has not previously been characterized. We find that heterochromatin marked by H3K9me3, HP1α and HP1β correlates with strong compartmentalization. We demonstrate that disruption of DNA methyltransferase activity greatly remodels genome compartmentalization whereby domains lose H3K9me3-HP1α/β binding and acquire the neutrally interacting state while retaining late replication timing. Furthermore, we show that H3K9me3-HP1α/β heterochromatin is permissive to loop extrusion by cohesin but refractory to CTCF binding. Together, our work reveals a dynamic structural and organizational diversity of the silent portion of the genome and establishes connections between the regulation of chromatin state and chromosome organization, including an interplay between DNA methylation, compartmentalization and loop extrusion.

Activation of HIV-1 proviruses increases downstream chromatin accessibility

It is unclear how the activation of HIV-1 transcription affects chromatin structure. We interrogated chromatin organization both genome-wide and nearby HIV-1 integration sites using Hi-C and ATAC-seq. In conjunction, we analyzed the transcription of the HIV-1 genome and neighboring genes. We found that long-range chromatin contacts did not differ significantly between uninfected cells and those harboring an integrated HIV-1 genome, whether the HIV-1 genome was actively transcribed or inactive. Instead, the activation of HIV-1 transcription changes chromatin accessibility immediately downstream of the provirus, demonstrating that HIV-1 can alter local cellular chromatin structure. Finally, we examined HIV-1 and neighboring host gene transcripts with long-read sequencing and found populations of chimeric RNAs both virus-to-host and host-to-virus. Thus, multiomics profiling revealed that the activation of HIV-1 transcription led to local changes in chromatin organization and altered the expression of neighboring host genes.

CTCF: an R/bioconductor data package of human and mouse CTCF binding sites

Summary

CTCF (CCCTC-binding factor) is an 11-zinc-finger DNA binding protein which regulates much of the eukaryotic genome's 3D structure and function. The diversity of CTCF binding motifs has led to a fragmented landscape of CTCF binding data. We collected position weight matrices of CTCF binding motifs and defined strand-oriented CTCF binding sites in the human and mouse genomes, including the recent Telomere to Telomere and mm39 assemblies. We included selected experimentally determined and predicted CTCF binding sites, such as CTCF-bound cis-regulatory elements from SCREEN ENCODE. We recommend filtering strategies for CTCF binding motifs and demonstrate that liftOver is a viable alternative to convert CTCF coordinates between assemblies. Our comprehensive data resource and usage recommendations can serve to harmonize and strengthen the reproducibility of genomic studies utilizing CTCF binding data.

Availability and implementation

https://bioconductor.org/packages/CTCF. Companion website: https://dozmorovlab.github.io/CTCF/; Code to reproduce the analyses: https://github.com/dozmorovlab/CTCF.dev.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

The Interplay of Transcription and Genome Topology Programs T Cell Development and Differentiation

T cells are essential for mounting defense against various pathogens and malignantly transformed cells. Thymic development and peripheral T cell differentiation are highly orchestrated biological processes that require precise gene regulation. Higher-order genome organization on multiple scales, in the form of chromatin loops, topologically associating domains and compartments, provides pivotal control of T cell gene expression. CTCF and the cohesin machinery are ubiquitously expressed architectural proteins responsible for establishing chromatin structures. Recent studies indicate that transcription factors, such as T lineage-defining Tcf1 and TCR-induced Batf, may have intrinsic ability and/or engage CTCF to shape chromatin architecture. In this article, we summarize current knowledge on the dynamic changes in genome topology that underlie normal or leukemic T cell development, CD4+ helper T cell differentiation, and CD8+ cytotoxic T cell functions. The knowledge lays a solid foundation for elucidating the causative link of spatial chromatin configuration to transcriptional and functional output in T cells.

Nickel-induced alterations to chromatin structure and function

Nickel (Ni), a heavy metal is prevalent in the atmosphere due to both natural and anthropogenic activities. Ni is a carcinogen implicated in the development of lung and nasal cancers in humans. Furthermore, Ni exposure is associated with a number of chronic lung diseases in humans including asthma, chronic bronchitis, emphysema, pulmonary fibrosis, pulmonary edema and chronic obstructive pulmonary disease (COPD). While Ni compounds are weak mutagens, a number of studies have demonstrated the potential of Ni to alter the epigenome, suggesting epigenomic dysregulation as an important underlying cause for its pathogenicity. In the eukaryotic nucleus, the DNA is organized in a three-dimensional (3D) space through assembly of higher order chromatin structures. Such an organization is critically important for transcription and other biological activities. Accumulating evidence suggests that by negatively affecting various cellular regulatory processes, Ni could potentially affect chromatin organization. In this review, we discuss the role of Ni in altering the chromatin architecture, which potentially plays a major role in Ni pathogenicity.

Chondrocyte De-Differentiation: Biophysical Cues to Nuclear Alterations

Autologous chondrocyte implantation (ACI) is a cell therapy to repair cartilage defects. In ACI a biopsy is taken from a non-load bearing area of the knee and expanded in-vitro. The expansion process provides the benefit of generating a large number of cells required for implantation; however, during the expansion these cells de-differentiate and lose their chondrocyte phenotype. In this review we focus on examining the de-differentiation phenotype from a mechanobiology and biophysical perspective, highlighting some of the nuclear mechanics and chromatin changes in chondrocytes seen during the expansion process and how this relates to the gene expression profile. We propose that manipulating chondrocyte nuclear architecture and chromatin organization will highlight mechanisms that will help to preserve the chondrocyte phenotype.

CW198 acts as a genetic insulator to block enhancer-promoter interaction in plants

Insulators in vertebrates play a role in genome architecture and orchestrate temporo-spatial enhancer-promoter interactions. In plants, insulators and their associated binding factors have not been documented as of yet, largely as a result of a lack of characterized insulators. In this study, we took a comprehensive strategy to identify and validate the enhancer-blocking insulator CW198. We show that a 1.08-kb CW198 fragment from Arabidopsis can, when interposed between an enhancer and a promoter, efficiently abrogate the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and tobacco plants. In plants, both transcriptional crosstalk and spreading of histone modifications were rarely detectable across CW198, which resembles the insulation property observed across the CTCF insulator in the mammalian genome. Taken together, our findings support that CW198 acts as an enhancer-blocking insulator in both Arabidopsis and tobacco. The significance of the present findings and their relevance to the mitigation of mutual interference between enhancers and promoters, as well as multiple promoters in transgenes, is discussed.

From Nucleus to Organs: Insights of Aryl Hydrocarbon Receptor Molecular Mechanisms

The aryl hydrocarbon receptor (AHR) is a markedly established regulator of a plethora of cellular and molecular processes. Its initial role in the detoxification of xenobiotic compounds has been partially overshadowed by its involvement in homeostatic and organ physiology processes. In fact, the discovery of its ability to bind specific target regulatory sequences has allowed for the understanding of how AHR modulates such processes. Thereby, AHR presents functions in transcriptional regulation, chromatin architecture modifications and participation in different key signaling pathways. Interestingly, such fields of influence end up affecting organ and tissue homeostasis, including regenerative response both to endogenous and exogenous stimuli. Therefore, from classical spheres such as canonical transcriptional regulation in embryonic development, cell migration, differentiation or tumor progression to modern approaches in epigenetics, senescence, immune system or microbiome, this review covers all aspects derived from the balance between regulation/deregulation of AHR and its physio-pathological consequences.

CTCF loss induces giant lamellar bodies in Purkinje cell dendrites

CCCTC-binding factor (CTCF) has a key role in higher-order chromatin architecture that is important for establishing and maintaining cell identity by controlling gene expression. In the mature cerebellum, CTCF is highly expressed in Purkinje cells (PCs) as compared with other cerebellar neurons. The cerebellum plays an important role in motor function by regulating PCs, which are the sole output neurons, and defects in PCs cause motor dysfunction. However, the role of CTCF in PCs has not yet been explored. Here we found that the absence of CTCF in mouse PCs led to progressive motor dysfunction and abnormal dendritic morphology in those cells, which included dendritic self-avoidance defects and a proximal shift in the climbing fibre innervation territory on PC dendrites. Furthermore, we found the peculiar lamellar structures known as "giant lamellar bodies" (GLBs), which have been reported in PCs of patients with Werdnig-Hoffman disease, 13q deletion syndrome, and Krabbe disease. GLBs are localized to PC dendrites and are assumed to be associated with neurodegeneration. They have been noted, however, only in case reports following autopsy, and reports of their existence have been very limited. Here we show that GLBs were reproducibly formed in PC dendrites of a mouse model in which CTCF was deleted. GLBs were not noted in PC dendrites at infancy but instead developed over time. In conjunction with GLB development in PC dendrites, the endoplasmic reticulum was almost absent around the nuclei, the mitochondria were markedly swollen and their cristae had decreased drastically, and almost all PCs eventually disappeared as severe motor deficits manifested. Our results revealed the important role of CTCF during normal development and in maintaining PCs and provide new insights into the molecular mechanism of GLB formation during neurodegenerative disease.

Multilevel view on chromatin architecture alterations in cancer

Chromosomes inside the nucleus are not located in the form of linear molecules. Instead, there is a complex multilevel genome folding that includes nucleosomes packaging, formation of chromatin loops, domains, compartments, and finally, chromosomal territories. Proper spatial organization play an essential role for the correct functioning of the genome, and is therefore dynamically changed during development or disease. Here we discuss how the organization of the cancer cell genome differs from the healthy genome at various levels. A better understanding of how malignization affects genome organization and long-range gene regulation will help to reveal the molecular mechanisms underlying cancer development and evolution.

The 3D enhancer network of the developing T cell genome is shaped by SATB1

Mechanisms of tissue-specific gene expression regulation via 3D genome organization are poorly understood. Here we uncover the regulatory chromatin network of developing T cells and identify SATB1, a tissue-specific genome organizer, enriched at the anchors of promoter-enhancer loops. We have generated a T-cell specific Satb1 conditional knockout mouse which allows us to infer the molecular mechanisms responsible for the deregulation of its immune system. H3K27ac HiChIP and Hi-C experiments indicate that SATB1-dependent promoter-enhancer loops regulate expression of master regulator genes (such as Bcl6), the T cell receptor locus and adhesion molecule genes, collectively being critical for cell lineage specification and immune system homeostasis. SATB1-dependent regulatory chromatin loops represent a more refined layer of genome organization built upon a high-order scaffold provided by CTCF and other factors. Overall, our findings unravel the function of a tissue-specific factor that controls transcription programs, via spatial chromatin arrangements complementary to the chromatin structure imposed by ubiquitously expressed genome organizers.