The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening.

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Role of the CTCF Binding Site in Human T-Cell Leukemia Virus-1 Pathogenesis.

During HTLV-1 infection, the virus integrates into the host cell genome as a provirus with a single CCCTC binding protein (CTCF) binding site (vCTCF-BS), which acts as an insulator between transcriptionally active and inactive regions. Previous studies have shown that the vCTCF-BS is important for maintenance of chromatin structure, regulation of viral expression, and DNA and histone methylation. Here, we show that the vCTCF-BS also regulates viral infection and pathogenesis in vivo in a humanized (Hu) mouse model of adult T-cell leukemia/lymphoma. Three cell lines were used to initiate infection of the Hu-mice, i) HTLV-1-WT which carries an intact HTLV-1 provirus genome, ii) HTLV-1-CTCF, which contains a provirus with a mutated vCTCF-BS which abolishes CTCF binding, and a stop codon immediate upstream of the mutated vCTCF-BS which deletes the last 23 amino acids of p12, and iii) HTLV-1-p12stop that contains the intact vCTCF-BS, but retains the same stop codon in p12 as in the HTLV-1-CTCF cell line. Hu-mice were infected with mitomycin treated or irradiated HTLV-1 producing cell lines. There was a delay in pathogenicity when Hu-mice were infected with the CTCF virus compared to mice infected with either p12 stop or WT virus. Proviral load (PVL), spleen weights, and CD4 T cell counts were significantly lower in HTLV-1-CTCF infected mice compared to HTLV-1-p12stop infected mice. Furthermore, we found a direct correlation between the PVL in peripheral blood and death of HTLV-1-CTCF infected mice. In cell lines, we found that the vCTCF-BS regulates Tax expression in a time-dependent manner. The scRNAseq analysis of splenocytes from infected mice suggests that the vCTCF-BS plays an important role in activation and expansion of T lymphocytes in vivo. Overall, these findings indicate that the vCTCF-BS regulates Tax expression, proviral load, and HTLV pathogenicity in vivo.

Hypoxia-induced CTCF promotes EMT in breast cancer.

Cancer cells experiencing hypoxic stress employ epithelial-mesenchymal transition (EMT) to undergo metastasis through rewiring of the chromatin landscape, epigenetics, and importantly, gene expression. Here, we showed that hypoxia modulates the epigenetic landscape on CTCF promoter and upregulates its expression. Hypoxia-driven epigenetic regulation, specifically DNA demethylation mediated by TET2, is a prerequisite for CTCF induction. Mechanistically, in hypoxic conditions, Hypoxia-inducible factor 1-alpha (HIF1α) binds to the unmethylated CTCF promoter, causing transcriptional upregulation. Further, we uncover the pivotal role of CTCF in promoting EMT as loss of CTCF abrogated invasiveness of hypoxic breast cancer cells. These findings highlight the functional contribution of HIF1α-CTCF axis in promoting EMT in hypoxic breast cancer cells. Lastly, CTCF expression is alleviated and the potential for EMT is diminished when the HIF1α binding is particularly disrupted through the dCas9-DNMT3A system-mediated maintenance of DNA methylation on the CTCF promoter. This axis may offer a unique therapeutic target in breast cancer.

Dissecting SOX2 expression and function reveals an association with multiple signaling pathways during embryonic development and in cancer progression.

SRY (Sex Determining Region) box 2 (SOX2) is an essential transcription factor that plays crucial roles in activating genes involved in pre- and post-embryonic development, adult tissue homeostasis, and lineage specifications. SOX2 maintains the self-renewal property of stem cells and is involved in the generation of induced pluripotency stem cells. SOX2 protein contains a particular high-mobility group domain that enables SOX2 to achieve the capacity to participate in a broad variety of functions. The information about the involvement of SOX2 with gene regulatory elements, signaling networks, and microRNA is gradually emerging, and the higher expression of SOX2 is functionally relevant to various cancer types. SOX2 facilitates the oncogenic phenotype via cellular proliferation and enhancement of invasive tumor properties. Evidence are accumulating in favor of three dimensional (higher order) folding of chromatin and epigenetic control of the SOX2 gene by chromatin modifications, which implies that the expression level of SOX2 can be modulated by epigenetic regulatory mechanisms, specifically, via DNA methylation and histone H3 modification. In view of this, and to focus further insights into the roles SOX2 plays in physiological functions, involvement of SOX2 during development, precisely, the advances of our knowledge in pre- and post-embryonic development, and interactions of SOX2 in this scenario with various signaling pathways in tumor development and cancer progression, its potential as a therapeutic target against many cancers are summarized and discussed in this article.

Regulation of 3D genome organization during T cell activation.

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

Overexpression of programmed cell death ligand 1 reduces LPS-induced inflammatory cytokine upregulation and enhances osteo/odontogenic-differentiation of human dental pulp stem cells via upregulation of CCCTC-binding factor.

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

A novel cis-regulatory element regulates αD and αA-globin gene expression in chicken erythroid cells.

Background: Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis. Methods: Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An α-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the α-globin domain in erythroid chicken cells. Results: Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An α-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over αD and αA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that α-globin CRE-2 is required for the induction of αD and αA. Analysis of TF binding motifs at α-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding. Conclusion: Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells.

Differential 3D genome architecture and imprinted gene expression: cause or consequence?

Imprinted genes provide an attractive paradigm to unravel links between transcription and genome architecture. The parental allele-specific expression of these essential genes - which are clustered in chromosomal domains - is mediated by parental methylation imprints at key regulatory DNA sequences. Recent chromatin conformation capture (3C)-based studies show differential organization of topologically associating domains between the parental chromosomes at imprinted domains, in embryonic stem and differentiated cells. At several imprinted domains, differentially methylated regions show allelic binding of the insulator protein CTCF, and linked focal retention of cohesin, at the non-methylated allele only. This generates differential patterns of chromatin looping between the parental chromosomes, already in the early embryo, and thereby facilitates the allelic gene expression. Recent research evokes also the opposite scenario, in which allelic transcription contributes to the differential genome organization, similarly as reported for imprinted X chromosome inactivation. This may occur through epigenetic effects on CTCF binding, through structural effects of RNA Polymerase II, or through imprinted long non-coding RNAs that have chromatin repressive functions. The emerging picture is that epigenetically-controlled differential genome architecture precedes and facilitates imprinted gene expression during development, and that at some domains, conversely, the mono-allelic gene expression also influences genome architecture.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales.

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization.

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

YTHDC1-Mediated lncRNA MSC-AS1 m6A Modification Potentiates Laryngeal Squamous Cell Carcinoma Development via Repressing ATXN7 Transcription.

Laryngeal squamous cell carcinoma (LSCC) has the highest mortality rate among head and neck squamous cell carcinoma. This study was designed to investigate the biological effect of long noncoding RNA (lncRNA) MSC antisense RNA 1 (MSC-AS1) on LSCC development and the underlying mechanism. The expression and prognostic value of lncRNAs in head and neck squamous cell carcinoma were predicted in the bioinformatics tools. The overexpression of MSC-AS1 in LSCC patients predicted a poor prognosis. Depletion of MSC-AS1 using shRNA repressed the malignant phenotype of AMC-HN-8 and TU-177 cells. MSC-AS1, mainly localized in the nucleus, interacted closely with the transcription factor CCCTC-binding factor (CTCF). CTCF played anti-tumor effects in vitro and in vivo. Ataxin-7 (ATXN7) was predicted to be a downstream target of CTCF, whose expression was negatively controlled by MSC-AS1. MSC-AS1 was found to block the expression of CTCF, thereby repressing ATXN7. Finally, MSC-AS1 overexpression in LSCC was governed by YTH domain-containing protein 1 (YTHDC1)-mediated m6A modification. In summary, our research identified the YTHDC1/MSC-AS1/CTCF/ATXN7 axis in LSCC development, which indicated that MSC-AS1 is an attractive biomarker in the LSCC treatment.

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.

LATS1 controls CTCF chromatin occupancy and hormonal response of 3D-grown breast cancer cells.

The cancer epigenome has been studied in cells cultured in two-dimensional (2D) monolayers, but recent studies highlight the impact of the extracellular matrix and the three-dimensional (3D) environment on multiple cellular functions. Here, we report the physical, biochemical, and genomic differences between T47D breast cancer cells cultured in 2D and as 3D spheroids. Cells within 3D spheroids exhibit a rounder nucleus with less accessible, more compacted chromatin, as well as altered expression of ~2000 genes, the majority of which become repressed. Hi-C analysis reveals that cells in 3D are enriched for regions belonging to the B compartment, have decreased chromatin-bound CTCF and increased fusion of topologically associating domains (TADs). Upregulation of the Hippo pathway in 3D spheroids results in the activation of the LATS1 kinase, which promotes phosphorylation and displacement of CTCF from DNA, thereby likely causing the observed TAD fusions. 3D cells show higher chromatin binding of progesterone receptor (PR), leading to an increase in the number of hormone-regulated genes. This effect is in part mediated by LATS1 activation, which favors cytoplasmic retention of YAP and CTCF removal.

ADNP modulates SINE B2-derived CTCF-binding sites during blastocyst formation in mice.

CTCF is crucial for chromatin structure and transcription regulation in early embryonic development. However, the kinetics of CTCF chromatin occupation in preimplantation embryos have remained unclear. In this study, we used CUT&RUN technology to investigate CTCF occupancy in mouse preimplantation development. Our findings revealed that CTCF begins binding to the genome prior to zygotic genome activation (ZGA), with a preference for CTCF-anchored chromatin loops. Although the majority of CTCF occupancy is consistently maintained, we identified a specific set of binding sites enriched in the mouse-specific short interspersed element (SINE) family B2 that are restricted to the cleavage stages. Notably, we discovered that the neuroprotective protein ADNP counteracts the stable association of CTCF at SINE B2-derived CTCF-binding sites. Knockout of Adnp in the zygote led to impaired CTCF binding signal recovery, failed deposition of H3K9me3, and transcriptional derepression of SINE B2 during the morula-to-blastocyst transition, which further led to unfaithful cell differentiation in embryos around implantation. Our analysis highlights an ADNP-dependent restriction of CTCF binding during cell differentiation in preimplantation embryos. Furthermore, our findings shed light on the functional importance of transposable elements (TEs) in promoting genetic innovation and actively shaping the early embryo developmental process specific to mammals.

RPL35A drives ovarian cancer progression by promoting the binding of YY1 to CTCF promoter.

Ovarian cancer is one of the most common gynaecological malignancies with poor prognosis and lack of effective treatment. The improvement of the situation of ovarian cancer urgently requires the exploration of its molecular mechanism to develop more effective molecular targeted drugs. In this study, the role of human ribosomal protein l35a (RPL35A) in ovarian cancer was explored in vitro and in vivo. Our data identified that RPL35A expression was abnormally elevated in ovarian cancer. Clinically, high expression of RPL35A predicted short survival and poor TNM staging in patients with ovarian cancer. Functionally, RPL35A knock down inhibited ovarian cancer cell proliferation and migration, enhanced apoptosis, while overexpression had the opposite effect. Mechanically, RPL35A promoted the direct binding of transcription factor YY1 to CTCF in ovarian cancer cells. Consistently, RPL35A regulated ovarian cancer progression depending on CTCF in vitro and in vivo. Furthermore, RPL35A affected the proliferation and apoptosis of ovarian cancer cells through PPAR signalling pathway. In conclusion, RPL35A drove ovarian cancer progression by promoting the binding of YY1 and CTCF promoter, and inhibiting this process may be an effective strategy for targeted therapy of this disease.

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.

DeepCBS: shedding light on the impact of mutations occurring at CTCF binding sites.

CTCF-mediated chromatin loops create insulated neighborhoods that constrain promoter-enhancer interactions, serving as a unit of gene regulation. Disruption of the CTCF binding sites (CBS) will lead to the destruction of insulated neighborhoods, which in turn can cause dysregulation of the contained genes. In a recent study, it is found that CTCF/cohesin binding sites are a major mutational hotspot in the cancer genome. Mutations can affect CTCF binding, causing the disruption of insulated neighborhoods. And our analysis reveals a significant enrichment of well-known proto-oncogenes in insulated neighborhoods with mutations specifically occurring in anchor regions. It can be assumed that some mutations disrupt CTCF binding, leading to the disruption of insulated neighborhoods and subsequent activation of proto-oncogenes within these insulated neighborhoods. To explore the consequences of such mutations, we develop DeepCBS, a computational tool capable of analyzing mutations at CTCF binding sites, predicting their influence on insulated neighborhoods, and investigating the potential activation of proto-oncogenes. Futhermore, DeepCBS is applied to somatic mutation data of liver cancer. As a result, 87 mutations that disrupt CTCF binding sites are identified, which leads to the identification of 237 disrupted insulated neighborhoods containing a total of 135 genes. Integrative analysis of gene expression differences in liver cancer further highlights three genes: ARHGEF39, UBE2C and DQX1. Among them, ARHGEF39 and UBE2C have been reported in the literature as potential oncogenes involved in the development of liver cancer. The results indicate that DQX1 may be a potential oncogene in liver cancer and may contribute to tumor immune escape. In conclusion, DeepCBS is a promising method to analyze impacts of mutations occurring at CTCF binding sites on the insulator function of CTCF, with potential extensions to shed light on the effects of mutations on other functions of CTCF.

G-quadruplex landscape and its regulation revealed by a new antibody capture method.

Our understanding of DNA G-quadruplexes (G4s) from in vitro studies has been complemented by genome-wide G4 landscapes from cultured cells. Conventionally, the formation of G4s is accepted to depend on G-repeats such that they form tetrads. However, genome-wide G4s characterized through high-throughput sequencing suggest that these structures form at a large number of regions with no such canonical G4-forming signatures. Many G4-binding proteins have been described with no evidence for any protein that binds to and stabilizes G4s. It remains unknown what fraction of G4s formed in human cells are protein-bound. The G4-chromatin immunoprecipitation (G4-ChIP) method hitherto employed to describe G4 landscapes preferentially reports G4s that get crosslinked to proteins in their proximity. Our current understanding of the G4 landscape is biased against representation of G4s which escape crosslinking as they are not stabilized by protein-binding and presumably transient. We report a protocol that captures G4s from the cells efficiently without any bias as well as eliminates the detection of G4s formed artifactually on crosslinked sheared chromatin post-fixation. We discover that G4s form sparingly at SINEs. An application of this method shows that depletion of a repeat-binding protein CGGBP1 enhances net G4 capture at CGGBP1-dependent CTCF-binding sites and regions of sharp interstrand G/C-skew transitions. Thus, we present an improved method for G4 landscape determination and by applying it we show that sequence property-specific constraints of the nuclear environment mitigate G4 formation.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction.

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.