CTCF and cellular heterogeneity

Deep learning models incorporating endogenous factors beyond DNA sequences improve the prediction accuracy of base editing outcomes

Adenine base editors (ABEs) and cytosine base editors (CBEs) enable the single nucleotide editing of targeted DNA sites avoiding generation of double strand breaks, however, the genomic features that influence the outcomes of base editing in vivo still remain to be characterized. High-throughput datasets from lentiviral integrated libraries were used to investigate the sequence features affecting base editing outcomes, but the effects of endogenous factors beyond the DNA sequences are still largely unknown. Here the base editing outcomes of ABE and CBE were evaluated in mammalian cells for 5012 endogenous genomic sites and 11,868 genome-integrated target sequences, with 4654 genomic sites sharing the same target sequences. The comparative analyses revealed that the editing outcomes of ABE and CBE at endogenous sites were substantially different from those obtained using genome-integrated sequences. We found that the base editing efficiency at endogenous target sites of both ABE and CBE was influenced by endogenous factors, including epigenetic modifications and transcriptional activity. A deep-learning algorithm referred as BE_Endo, was developed based on the endogenous factors and sequence information from our genomic datasets, and it yielded unprecedented accuracy in predicting the base editing outcomes. These findings along with the developed computational algorithms may facilitate future application of BEs for scientific research and clinical gene therapy.

An oligodendrocyte silencer element underlies the pathogenic impact of lamin B1 structural variants

The role of non-coding regulatory elements and how they might contribute to tissue type specificity of disease phenotypes is poorly understood. Autosomal Dominant Leukodystrophy (ADLD) is a fatal, adult-onset, neurological disorder that is characterized by extensive CNS demyelination. Most cases of ADLD are caused by tandem genomic duplications involving the lamin B1 gene ( LMNB1 ) while a small subset are caused by genomic deletions upstream of the gene. Utilizing data from recently identified families that carry LMNB1 gene duplications but do not exhibit demyelination, ADLD patient tissues, CRISPR modified cell lines and mouse models, we have identified a novel silencer element that is lost in ADLD patients and that specifically targets overexpression to oligodendrocytes. This element consists of CTCF binding sites that mediate three-dimensional chromatin looping involving the LMNB1 and the recruitment of the PRC2 repressor complex. Loss of the silencer element in ADLD identifies a previously unknown role for silencer elements in tissue specificity and disease causation.

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.

3D chromatin architecture and transcription regulation in cancer

Chromatin has distinct three-dimensional (3D) architectures important in key biological processes, such as cell cycle, replication, differentiation, and transcription regulation. In turn, aberrant 3D structures play a vital role in developing abnormalities and diseases such as cancer. This review discusses key 3D chromatin structures (topologically associating domain, lamina-associated domain, and enhancer-promoter interactions) and corresponding structural protein elements mediating 3D chromatin interactions [CCCTC-binding factor, polycomb group protein, cohesin, and Brother of the Regulator of Imprinted Sites (BORIS) protein] with a highlight of their associations with cancer. We also summarise the recent development of technologies and bioinformatics approaches to study the 3D chromatin interactions in gene expression regulation, including crosslinking and proximity ligation methods in the bulk cell population (ChIA-PET and HiChIP) or single-molecule resolution (ChIA-drop), and methods other than proximity ligation, such as GAM, SPRITE, and super-resolution microscopy techniques.

Long noncoding RNA ELDR promotes cell cycle progression in normal oral keratinocytes through induction of a CTCF-FOXM1-AURKA signaling axis

Long noncoding RNAs (lncRNAs) have gained widespread attention as a new layer of regulation in biological processes during development and disease. The lncRNA ELDR (EGFR long noncoding downstream RNA) was recently shown to be highly expressed in oral cancers as compared to adjacent nontumor tissue, and we previously reported that ELDR may be an oncogene as inhibition of ELDR reduces tumor growth in oral cancer models. Furthermore, overexpression of ELDR induces proliferation and colony formation in normal oral keratinocytes (NOKs). In this study, we examined in further detail how ELDR drives the neoplastic transformation of normal keratinocytes. We performed RNA-seq analysis on NOKs stably expressing ELDR (NOK-ELDR), which revealed that ELDR enhances the expression of cell cycle-related genes. Expression of Aurora kinase A and its downstream targets Polo-like kinase 1, cell division cycle 25C, cyclin-dependent kinase 1, and cyclin B1 (CCNB1) are significantly increased in NOK-ELDR cells, suggesting induction of G2/M progression. We further identified CCCTC-binding factor (CTCF) as a binding partner of ELDR in NOK-ELDR cells. We show that ELDR stabilizes CTCF and increases its expression. Finally, we demonstrate the ELDR-CTCF axis upregulates transcription factor Forkhead box M1, which induces Aurora kinase A expression and downstream G2/M transition. These findings provide mechanistic insights into the role of the lncRNA ELDR as a potential driver of oral cancer during neoplastic transformation of normal keratinocytes.

A novel workflow for the qualitative analysis of DNA methylation data

•

A novel R package (EpiStatProfiler) for the qualitative analysis of DNA methylation data.

•

A novel workflow for the analysis of CG and non-CG epialleles starting from any type of bisulfite sequencing data.

•

EpiStatProfiler can perform strand-specific characterization of epialleles composition.

•

Important loci can be annotated along with their biological role and potential functions.

•

EpiStatProfiler has the ability to identify loci whose epiallelic profile is associated with disease pathogenesis.

DNA methylation is an epigenetic modification that plays a pivotal role in major biological mechanisms, such as gene regulation, genomic imprinting, and genome stability. Different combinations of methylated cytosines for a given DNA locus generate different epialleles and alterations of these latter have been associated with several pathological conditions. Existing computational methods and statistical tests relevant to DNA methylation analysis are mostly based on the comparison of average CpG sites methylation levels and they often neglect non-CG methylation. Here, we present EpiStatProfiler, an R package that allows the analysis of CpG and non-CpG based epialleles starting from bisulfite sequencing data through a collection of dedicated extraction functions and statistical tests. EpiStatProfiler is provided with a set of useful auxiliary features, such as customizable genomic ranges, strand-specific epialleles analysis, locus annotation and gene set enrichment analysis. We showcase the package functionalities on two public datasets by identifying putative relevant loci in mice harboring the Huntington’s disease-causing Htt gene mutation and in Ctcf +/− mice compared to their wild-type counterparts. To our knowledge, EpiStatProfiler is the first package providing functionalities dedicated to the analysis of epialleles composition derived from any kind of bisulfite sequencing experiment.

CGGBP1-dependent CTCF-binding sites restrict ectopic transcription

Binding sites of the chromatin regulator protein CTCF function as important landmarks in the human genome. The recently characterized CTCF-binding sites at LINE-1 repeats depend on another repeat-regulatory protein CGGBP1. These CGGBP1-dependent CTCF-binding sites serve as potential barrier elements for epigenetic marks such as H3K9me3. Such CTCF-binding sites are associated with asymmetric H3K9me3 levels as well as RNA levels in their flanks. The functions of these CGGBP1-dependent CTCF-binding sites remain unknown. By performing targeted studies on candidate CGGBP1-dependent CTCF-binding sites cloned in an SV40 promoter-enhancer episomal system we show that these regions act as inhibitors of ectopic transcription from the SV40 promoter. CGGBP1-dependent CTCF-binding sites that recapitulate their genomic function of loss of CTCF binding upon CGGBP1 depletion and H3K9me3 asymmetry in immediate flanks are also the ones that show the strongest inhibition of ectopic transcription. By performing a series of strand-specific reverse transcription PCRs we demonstrate that this ectopic transcription results in the synthesis of RNA from the SV40 promoter in a direction opposite to the downstream reporter gene in a strand-specific manner. The unleashing of the bidirectionality of the SV40 promoter activity and a breach of the transcription barrier seems to depend on depletion of CGGBP1 and loss of CTCF binding proximal to the SV40 promoter. RNA-sequencing reveals that CGGBP1-regulated CTCF-binding sites act as barriers to transcription at multiple locations genome-wide. These findings suggest a role of CGGBP1-dependent binding sites in restricting ectopic transcription.

The epigenetic basis of cellular heterogeneity

Single-cell sequencing-based methods for profiling gene transcript levels have revealed substantial heterogeneity in expression levels among morphologically indistinguishable cells. This variability has important functional implications for tissue biology and disease states such as cancer. Mapping of epigenomic information such as chromatin accessibility, nucleosome positioning, histone tail modifications and enhancer-promoter interactions in both bulk-cell and single-cell samples has shown that these characteristics of chromatin state contribute to expression or repression of associated genes. Advances in single-cell epigenomic profiling methods are enabling high-resolution mapping of chromatin states in individual cells. Recent studies using these techniques provide evidence that variations in different aspects of chromatin organization collectively define gene expression heterogeneity among otherwise highly similar cells.

Effects of DNA Methylation on TFs in Human Embryonic Stem Cells

DNA methylation is an important epigenetic mechanism for gene regulation. The conventional view of DNA methylation is that DNA methylation could disrupt protein-DNA interactions and repress gene expression. Several recent studies reported that DNA methylation could alter transcription factors (TFs) binding sequence specificity in vitro. Here, we took advantage of the large sets of ChIP-seq data for TFs and whole-genome bisulfite sequencing data in many cell types to perform a systematic analysis of the protein-DNA methylation in vivo. We observed that many TFs could bind methylated DNA regions, especially in H1-hESC cells. By locating binding sites, we confirmed that some TFs could bind to methylated CpGs directly. The different proportion of CpGs at TF binding specificity motifs in different methylation statuses shows that some TFs are sensitive to methylation and some could bind to the methylated DNA with different motifs, such as CEBPB and CTCF. At the same time, TF binding could interactively alter local DNA methylation. The TF hypermethylation binding sites extensively overlap with enhancers. And we also found that some DNase I hypersensitive sites were specifically hypermethylated in H1-hESC cells. At last, compared with TFs' binding regions in multiple cell types, we observed that CTCF binding to high methylated regions in H1-hESC were not conservative. These pieces of evidence indicate that TFs that bind to hypermethylation DNA in H1-hESC cells may associate with enhancers to regulate special biological functions.

CTCF As an Example of DNA-Binding Transcription Factors Containing Clusters of C2H2-Type Zinc Fingers

In mammals, most of the boundaries of topologically associating domains and all well-studied insulators are rich in binding sites for the CTCF protein. According to existing experimental data, CTCF is a key factor in the organization of the architecture of mammalian chromosomes. A characteristic feature of the CTCF is that the central part of the protein contains a cluster consisting of eleven domains of C2H2-type zinc fingers, five of which specifically bind to a long DNA sequence conserved in most animals. The class of transcription factors that carry a cluster of C2H2-type zinc fingers consisting of five or more domains (C2H2 proteins) is widely represented in all groups of animals. The functions of most C2H2 proteins still remain unknown. This review presents data on the structure and possible functions of these proteins, using the example of the vertebrate CTCF protein and several well- characterized C2H2 proteins in Drosophila and mammals.

Dynamic regulation of CTCF stability and sub-nuclear localization in response to stress

The nuclear protein CCCTC-binding factor (CTCF) has diverse roles in chromatin architecture and gene regulation. Functionally, CTCF associates with thousands of genomic sites and interacts with proteins, such as cohesin, or non-coding RNAs to facilitate specific transcriptional programming. In this study, we examined CTCF during the cellular stress response in human primary cells using immune-blotting, quantitative real time-PCR, chromatin immunoprecipitation-sequence (ChIP-seq) analysis, mass spectrometry, RNA immunoprecipitation-sequence analysis (RIP-seq), and Airyscan confocal microscopy. Unexpectedly, we found that CTCF is exquisitely sensitive to diverse forms of stress in normal patient-derived human mammary epithelial cells (HMECs). In HMECs, a subset of CTCF protein forms complexes that localize to Serine/arginine-rich splicing factor (SC-35)-containing nuclear speckles. Upon stress, this species of CTCF protein is rapidly downregulated by changes in protein stability, resulting in loss of CTCF from SC-35 nuclear speckles and changes in CTCF-RNA interactions. Our ChIP-seq analysis indicated that CTCF binding to genomic DNA is largely unchanged. Restoration of the stress-sensitive pool of CTCF protein abundance and re-localization to nuclear speckles can be achieved by inhibition of proteasome-mediated degradation. Surprisingly, we observed the same characteristics of the stress response during neuronal differentiation of human pluripotent stem cells (hPSCs). CTCF forms stress-sensitive complexes that localize to SC-35 nuclear speckles during a specific stage of neuronal commitment/development but not in differentiated neurons. We speculate that these particular CTCF complexes serve a role in RNA processing that may be intimately linked with specific genes in the vicinity of nuclear speckles, potentially to maintain cells in a certain differentiation state, that is dynamically regulated by environmental signals. The stress-regulated activity of CTCF is uncoupled in persistently stressed, epigenetically re-programmed "variant" HMECs and certain cancer cell lines. These results reveal new insights into CTCF function in cell differentiation and the stress-response with implications for oxidative damage-induced cancer initiation and neuro-degenerative diseases.

Cancer evolution: A means by which tumors evade treatment

Although various methods have been tried to study and treat cancer, the cancer remains a major challenge for human medicine today. One important reason for this is the presence of cancer evolution. Cancer evolution is a process in which tumor cells adapt to the external environment, which can suppress the human immune system's ability to recognize and attack tumors, and also reduce the reproducibility of cancer research. Among them, heterogeneity of the tumor provides intrinsic motivation for this process. Recently, with the development of related technologies such as liquid biopsy, more and more knowledge about cancer evolution has been gained and interest in this topic has also increased. Therefore, starting from the causes of tumorigenesis, this paper introduces several tumorigenesis processes and pathways, as well as treatment options for different targets.

Exploring the changing landscape of cell-to-cell variation after CTCF knockdown via single cell RNA-seq

BACKGROUND

CCCTC-Binding Factor (CTCF), also known as 11-zinc finger protein, participates in many cellular processes, including insulator activity, transcriptional regulation and organization of chromatin architecture. Based on single cell flow cytometry and single cell RNA-FISH analyses, our previous study showed that deletion of CTCF binding site led to a significantly increase of cellular variation of its target gene. However, the effect of CTCF on genome-wide landscape of cell-to-cell variation remains unclear.

RESULTS

We knocked down CTCF in EL4 cells using shRNA, and conducted single cell RNA-seq on both wild type (WT) cells and CTCF-Knockdown (CTCF-KD) cells using Fluidigm C1 system. Principal component analysis of single cell RNA-seq data showed that WT and CTCF-KD cells concentrated in two different clusters on PC1, indicating that gene expression profiles of WT and CTCF-KD cells were systematically different. Interestingly, GO terms including regulation of transcription, DNA binding, zinc finger and transcription factor binding were significantly enriched in CTCF-KD-specific highly variable genes, implying tissue-specific genes such as transcription factors were highly sensitive to CTCF level. The dysregulation of transcription factors potentially explains why knockdown of CTCF leads to systematic change of gene expression. In contrast, housekeeping genes such as rRNA processing, DNA repair and tRNA processing were significantly enriched in WT-specific highly variable genes, potentially due to a higher cellular variation of cell activity in WT cells compared to CTCF-KD cells. We further found that cellular variation-increased genes were significantly enriched in down-regulated genes, indicating CTCF knockdown simultaneously reduced the expression levels and increased the expression noise of its regulated genes.

CONCLUSIONS

To our knowledge, this is the first attempt to explore genome-wide landscape of cellular variation after CTCF knockdown. Our study not only advances our understanding of CTCF function in maintaining gene expression and reducing expression noise, but also provides a framework for examining gene function.