CTCF: master weaver of the genome

Motorized chain models of the ideal chromosome

An array of motor proteins consumes chemical energy in setting up the architectures of chromosomes. Here, we explore how the structure of ideal polymer chains is influenced by two classes of motors. The first class which we call "swimming motors" acts to propel the chromatin fiber through three-dimensional space. They represent a caricature of motors such as RNA polymerases. Previously, they have often been described by adding a persistent flow onto Brownian diffusion of the chain. The second class of motors, which we call "grappling motors" caricatures the loop extrusion processes in which segments of chromatin fibers some distance apart are brought together. We analyze these models using a self-consistent variational phonon approximation to a many-body Master equation incorporating motor activities. We show that whether the swimming motors lead to contraction or expansion depends on the susceptibility of the motors, that is, how their activity depends on the forces they must exert. Grappling motors in contrast to swimming motors lead to long-ranged correlations that resemble those first suggested for fractal globules and that are consistent with the effective interactions inferred by energy landscape analyses of Hi-C data on the interphase chromosome.

CTCF mutation at R567 causes developmental disorders via 3D genome rearrangement and abnormal neurodevelopment

The three-dimensional genome structure organized by CTCF is required for development. Clinically identified mutations in CTCF have been linked to adverse developmental outcomes. Nevertheless, the underlying mechanism remains elusive. In this investigation, we explore the regulatory roles of a clinically relevant R567W point mutation, located within the 11th zinc finger of CTCF, by introducing this mutation into both murine models and human embryonic stem cell-derived cortical organoid models. Mice with homozygous CTCFR567W mutation exhibit growth impediments, resulting in postnatal mortality, and deviations in brain, heart, and lung development at the pathological and single-cell transcriptome levels. This mutation induces premature stem-like cell exhaustion, accelerates the maturation of GABAergic neurons, and disrupts neurodevelopmental and synaptic pathways. Additionally, it specifically hinders CTCF binding to peripheral motifs upstream to the core consensus site, causing alterations in local chromatin structure and gene expression, particularly at the clustered protocadherin locus. Comparative analysis using human cortical organoids mirrors the consequences induced by this mutation. In summary, this study elucidates the influence of the CTCFR567W mutation on human neurodevelopmental disorders, paving the way for potential therapeutic interventions.

Yin Yang 1 regulates cohesin complex protein SMC3 in mouse hematopoietic stem cells

ABSTRACT

Yin Yang 1 (YY1) and structural maintenance of chromosomes 3 (SMC3) are 2 critical chromatin structural factors that mediate long-distance enhancer-promoter interactions and promote developmentally regulated changes in chromatin architecture in hematopoietic stem/progenitor cells (HSPCs). Although YY1 has critical functions in promoting hematopoietic stem cell (HSC) self-renewal and maintaining HSC quiescence, SMC3 is required for proper myeloid lineage differentiation. However, many questions remain unanswered regarding how YY1 and SMC3 interact with each other and affect hematopoiesis. We found that YY1 physically interacts with SMC3 and cooccupies with SMC3 at a large cohort of promoters genome wide, and YY1 deficiency deregulates the genetic network governing cell metabolism. YY1 occupies the Smc3 promoter and represses SMC3 expression in HSPCs. Although deletion of 1 Smc3 allele partially restores HSC numbers and quiescence in YY1 knockout mice, Yy1-/-Smc3+/- HSCs fail to reconstitute blood after bone marrow transplant. YY1 regulates HSC metabolic pathways and maintains proper intracellular reactive oxygen species levels in HSCs, and this regulation is independent of the YY1-SMC3 axis. Our results establish a distinct YY1-SMC3 axis and its impact on HSC quiescence and metabolism.

An atlas of chromatin landscape in KSHV-infected cells during de novo infection and reactivation

Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic γ-herpesvirus with a double-stranded DNA capable of establishing latent infection in the host cell. During latency, only a limited number of viral genes are expressed in infected host cells, and that helps the virus to evade host immune cell response. During primary infection, the KSHV genome is chromatinized and maintained as an episome, which is tethered to the host chromosome via Latency Associated Nuclear Antigen (LANA). The KSHV episome undergoes the same chromatin modification with the host cell chromosome and, therefore, is regulated by various epigenetic modifications, such as DNA methylation, histone methylation, and histone acetylation. The KSHV genome is also organized in a spatiotemporal manner by forming genomic loops, which enable simultaneous and coordinated control of dynamic gene transcription, particularly during the lytic replication phase. The genome-wide approaches and advancing bioinformatic tools have increased the resolution of studies on the dynamic transcriptional control and our understanding of KSHV latency-lytic switch regulation. We will summarize our current understanding of the epigenetic gene regulation on the KSHV chromatin.

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.

Effects of blood-processing protocols on cell-free DNA fragmentomics in plasma: Comparisons of one- and two-step centrifugations

BACKGROUND

Cell-free DNA (cfDNA) fragmentomic characteristics are promising analytes with abundant physiological signals for non-invasive disease diagnosis and monitoring. Previous studies on plasma cfDNA fragmentomics commonly employed a two-step centrifugation process for removing cell debris, involving a low-speed centrifugation followed by a high-speed centrifugation. However, the effects of centrifugation conditions on the analysis of cfDNA fragmentome remain uncertain.

METHODS

We collected blood samples from 10 healthy individuals and divided each sample into two aliquots for plasma preparation with one- and two-step centrifugation processes. We performed whole genome sequencing (WGS) of the plasma cfDNA in the two groups and comprehensively compared the cfDNA fragmentomic features. Additionally, we reanalyzed the fragmentomic features of cfDNA from 16 healthy individuals and 16 COVID-19 patients, processed through one- and two-step centrifugation in our previous study, to investigate the impact of centrifugation on disease signals.

RESULTS

Our results showed that there were no significant differences observed in the characteristics of nuclear cfDNA, including size, motif diversity score (MDS) of end motifs, and genome distribution, between plasma samples treated with one- and two-step centrifugation. The cfDNA size shortening in COVID-19 patients was observed in plasma samples with one- and two-step centrifugation methods. However, we observed a significantly higher relative abundance and longer size of cell-free mitochondrial DNA (mtDNA) in the one-step samples compared to the two-step samples. This difference in mtDNA caused by the one- and two-step centrifugation methods surpasses the pathological difference between COVID-19 patients and healthy individuals.

CONCLUSIONS

Our findings indicate that one-step low-speed centrifugation is a simple and potentially suitable method for analyzing nuclear cfDNA fragmentation characteristics. These results offer valuable guidance for cfDNA research in various clinical scenarios.

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

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

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.

RUNX1::ETO and CBFβ::MYH11 converge on aberrant activation of BCAT1 to confer a therapeutic vulnerability in core-binding factor-acute myeloid leukaemia

Effectively targeting transcription factors in therapeutic interventions remains challenging, especially in core-binding factor-acute myeloid leukaemia (CBF-AML) characterized by RUNX1::ETO and CBFβ::MYH11 fusions. However, recent studies have drawn attention towards aberrant amino acid metabolisms as actionable therapeutic targets. Here, by integrating the expression profile and genetic makeup in AML cohort, we found higher BCAT1 expression in CBF-AML patients compared with other subtypes. Metabolic profiling revealed that high BCAT1 expression led to reprogrammed branch amino acid metabolism in CBF-AML and was associated with sphingolipid pathway relating to the fitness of leukaemia cells, supported by transcriptomic profiling. Mechanistically, we demonstrated in cell lines and primary patient samples that BCAT1 was directly activated by RUNX1::ETO and CBFβ::MYH11 fusion proteins similarly in a RUNX1-dependent manner through rewiring chromatin conformation at the BCAT1 gene locus. Furthermore, BCAT1 inhibition resulted in blunted cell cycle, enhanced apoptosis and myeloid differentiation of CBF-AML cells in vitro, and alleviated leukaemia burden and prolonged survival in vivo. Importantly, pharmacological inhibition of BCAT1 using the specific inhibitor Gabapentin demonstrated therapeutic effects, as evidenced by delayed leukaemia progression and improved survival in vivo. In conclusion, our study uncovers BCAT1 as a genetic vulnerability and a promising targeted therapeutic opportunity for CBF-AML.

SATB1, senescence and senescence-related diseases

Aging leads to an accumulation of cellular mutations and damage, increasing the risk of senescence, apoptosis, and malignant transformation. Cellular senescence, which is pivotal in aging, acts as both a guard against cellular transformation and as a check against cancer progression. It is marked by stable cell cycle arrest, widespread macromolecular changes, a pro-inflammatory profile, and altered gene expression. However, it remains to be determined whether these differing subsets of senescent cells result from unique intrinsic programs or are influenced by their environmental contexts. Multiple transcription regulators and chromatin modifiers contribute to these alterations. Special AT-rich sequence-binding protein 1 (SATB1) stands out as a crucial regulator in this process, orchestrating gene expression by structuring chromatin into loop domains and anchoring DNA elements. This review provides an overview of cellular senescence and delves into the role of SATB1 in senescence-related diseases. It highlights SATB1's potential in developing antiaging and anticancer strategies, potentially contributing to improved quality of life and addressing aging-related diseases.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

Jump-starting life: balancing transposable element co-option and genome integrity in the developing mammalian embryo

Remnants of transposable elements (TEs) are widely expressed throughout mammalian embryo development. Originally infesting our genomes as selfish elements and acting as a source of genome instability, several of these elements have been co-opted as part of a complex system of genome regulation. Many TEs have lost transposition ability and their transcriptional potential has been tampered as a result of interactions with the host throughout evolutionary time. It has been proposed that TEs have been ultimately repurposed to function as gene regulatory hubs scattered throughout our genomes. In the early embryo in particular, TEs find a perfect environment of naïve chromatin to escape transcriptional repression by the host. As a consequence, it is thought that hosts found ways to co-opt TE sequences to regulate large-scale changes in chromatin and transcription state of their genomes. In this review, we discuss several examples of TEs expressed during embryo development, their potential for co-option in genome regulation and the evolutionary pressures on TEs and on our genomes.

Towards Development of the 4C-Based Method Detecting Interactions of Plasmid DNA with Host Genome

Chromosome conformation capture techniques have revolutionized our understanding of chromatin architecture and dynamics at the genome-wide scale. In recent years, these methods have been applied to a diverse array of species, revealing fundamental principles of chromosomal organization. However, structural organization of the extrachromosomal entities, like viral genomes or plasmids, and their interactions with the host genome, remain relatively underexplored. In this work, we introduce an enhanced 4C-protocol tailored for probing plasmid DNA interactions. We design specific plasmid vector and optimize protocol to allow high detection rate of contacts between the plasmid and host DNA.

Viral remodeling of the 4D nucleome

The dynamic spatial organization of genomes across time, referred to as the four-dimensional nucleome (4DN), is a key component of gene regulation and biological fate. Viral infections can lead to a reconfiguration of viral and host genomes, impacting gene expression, replication, latency, and oncogenic transformation. This review provides a summary of recent research employing three-dimensional genomic methods such as Hi-C, 4C, ChIA-PET, and HiChIP in virology. We review how viruses induce changes in gene loop formation between regulatory elements, modify chromatin accessibility, and trigger shifts between A and B compartments in the host genome. We highlight the central role of cellular chromatin organizing factors, such as CTCF and cohesin, that reshape the 3D structure of both viral and cellular genomes. We consider how viral episomes, viral proteins, and viral integration sites can alter the host epigenome and how host cell type and conditions determine viral epigenomes. This review consolidates current knowledge of the diverse host-viral interactions that impact the 4DN.

C2H2 Zinc Finger Transcription Factors Associated with Hemoglobinopathies

In humans, specific aberrations in β-globin results in sickle cell disease and β-thalassemia, symptoms of which can be ameliorated by increased expression of fetal globin (HbF). Two recent CRISPR-Cas9 screens, centered on ∼1500 annotated sequence-specific DNA binding proteins and performed in a human erythroid cell line that expresses adult hemoglobin, uncovered four groups of candidate regulators of HbF gene expression. They are (1) members of the nucleosome remodeling and deacetylase (NuRD) complex proteins that are already known for HbF control; (2) seven C2H2 zinc finger (ZF) proteins, including some (ZBTB7A and BCL11A) already known for directly silencing the fetal γ-globin genes in adult human erythroid cells; (3) a few other transcription factors of different structural classes that might indirectly influence HbF gene expression; and (4) DNA methyltransferase 1 (DNMT1) that maintains the DNA methylation marks that attract the MBD2-associated NuRD complex to DNA as well as associated histone H3 lysine 9 methylation. Here we briefly discuss the effects of these regulators, particularly C2H2 ZFs, in inducing HbF expression for treating β-hemoglobin disorders, together with recent advances in developing safe and effective small-molecule therapeutics for the regulation of this well-conserved hemoglobin switch.

A CTCF-dependent mechanism underlies the Hox timer: relation to a segmented body plan

During gastrulation, Hox genes are activated in a time-sequence that follows the order of the genes along their clusters. This property, which is observed in all animals that develop following a progressive rostral-to-caudal morphogenesis, is associated with changes in the chromatin structure and epigenetic profiles of Hox clusters, suggesting a process at least partly based on sequential gene accessibility. Here, we discuss recent work on this issue, as well as a possible mechanism based on the surprising conservation in both the distribution and orientation of CTCF sites inside vertebrate Hox clusters.

Single-cell RNA-seq data analysis reveals functionally relevant biomarkers of early brain development and their regulatory footprints in human embryonic stem cells (hESCs)

The complicated process of neuronal development is initiated early in life, with the genetic mechanisms governing this process yet to be fully elucidated. Single-cell RNA sequencing (scRNA-seq) is a potent instrument for pinpointing biomarkers that exhibit differential expression across various cell types and developmental stages. By employing scRNA-seq on human embryonic stem cells, we aim to identify differentially expressed genes (DEGs) crucial for early-stage neuronal development. Our focus extends beyond simply identifying DEGs. We strive to investigate the functional roles of these genes through enrichment analysis and construct gene regulatory networks to understand their interactions. Ultimately, this comprehensive approach aspires to illuminate the molecular mechanisms and transcriptional dynamics governing early human brain development. By uncovering potential links between these DEGs and intelligence, mental disorders, and neurodevelopmental disorders, we hope to shed light on human neurological health and disease. In this study, we have used scRNA-seq to identify DEGs involved in early-stage neuronal development in hESCs. The scRNA-seq data, collected on days 26 (D26) and 54 (D54), of the in vitro differentiation of hESCs to neurons were analyzed. Our analysis identified 539 DEGs between D26 and D54. Functional enrichment of those DEG biomarkers indicated that the up-regulated DEGs participated in neurogenesis, while the down-regulated DEGs were linked to synapse regulation. The Reactome pathway analysis revealed that down-regulated DEGs were involved in the interactions between proteins located in synapse pathways. We also discovered interactions between DEGs and miRNA, transcriptional factors (TFs) and DEGs, and between TF and miRNA. Our study identified 20 significant transcription factors, shedding light on early brain development genetics. The identified DEGs and gene regulatory networks are valuable resources for future research into human brain development and neurodevelopmental disorders.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

The impact of HTLV-1 expression on the 3D structure and expression of host chromatin

A typical HTLV-1-infected individual carries >104 different HTLV-1-infected T cell clones, each with a single-copy provirus integrated in a unique genomic site. We previously showed that the HTLV-1 provirus causes aberrant transcription in the flanking host genome and, by binding the chromatin architectural protein CTCF, forms abnormal chromatin loops with the host genome. However, it remained unknown whether these effects were exerted simply by the presence of the provirus or were induced by its transcription. To answer this question, we sorted HTLV-1-infected T-cell clones into cells positive or negative for proviral plus-strand expression, and then quantified host and provirus transcription using RNA-seq, and chromatin looping using quantitative chromosome conformation capture (q4C), in each cell population. We found that proviral plus-strand transcription induces aberrant transcription and splicing in the flanking genome but suppresses aberrant chromatin loop formation with the nearby host chromatin. Reducing provirus-induced host transcription with an inhibitor of transcriptional elongation allows recovery of chromatin loops in the plus-strand-expressing population. We conclude that aberrant host transcription induced by proviral expression causes temporary, reversible disruption of chromatin looping in the vicinity of the provirus.

Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

PARP1 as an Epigenetic Modulator: Implications for the Regulation of Host-Viral Dynamics

The principal understanding of the Poly(ADP-ribose) polymerase (PARP) regulation of genomes has been focused on its role in DNA repair; however, in the past few years, an additional role for PARPs and PARylation has emerged in regulating viral-host interactions. In particular, in the context of DNA virus infection, PARP1-mediated mechanisms of gene regulations, such as the involvement with cellular protein complexes responsible for the folding of the genome into the nucleus, the formation of chromatin loops connecting distant regulatory genomic regions, and other methods of transcriptional regulation, provide additional ways through which PARPs can modulate the function of both the host and the viral genomes during viral infection. In addition, potential viral amplification of the activity of PARPs on the host genome can contribute to the pathogenic effect of viral infection, such as viral-driven oncogenesis, opening the possibility that PARP inhibition may represent a potential therapeutic approach to target viral infection. This review will focus on the role of PARPs, particularly PARP1, in regulating the infection of DNA viruses.

Complementing Hi-C information for 3D chromatin reconstruction by ChromStruct

A multiscale method proposed elsewhere for reconstructing plausible 3D configurations of the chromatin in cell nuclei is recalled, based on the integration of contact data from Hi-C experiments and additional information coming from ChIP-seq, RNA-seq and ChIA-PET experiments. Provided that the additional data come from independent experiments, this kind of approach is supposed to leverage them to complement possibly noisy, biased or missing Hi-C records. When the different data sources are mutually concurrent, the resulting solutions are corroborated; otherwise, their validity would be weakened. Here, a problem of reliability arises, entailing an appropriate choice of the relative weights to be assigned to the different informational contributions. A series of experiments is presented that help to quantify the advantages and the limitations offered by this strategy. Whereas the advantages in accuracy are not always significant, the case of missing Hi-C data demonstrates the effectiveness of additional information in reconstructing the highly packed segments of the structure.

A cis-regulatory element regulates ERAP2 expression through autoimmune disease risk SNPs

Single-nucleotide polymorphisms (SNPs) near the ERAP2 gene are associated with various autoimmune conditions, as well as protection against lethal infections. Due to high linkage disequilibrium, numerous trait-associated SNPs are correlated with ERAP2 expression; however, their functional mechanisms remain unidentified. We show by reciprocal allelic replacement that ERAP2 expression is directly controlled by the splice region variant rs2248374. However, disease-associated variants in the downstream LNPEP gene promoter are independently associated with ERAP2 expression. Allele-specific conformation capture assays revealed long-range chromatin contacts between the gene promoters of LNPEP and ERAP2 and showed that interactions were stronger in patients carrying the alleles that increase susceptibility to autoimmune diseases. Replacing the SNPs in the LNPEP promoter by reference sequences lowered ERAP2 expression. These findings show that multiple SNPs act in concert to regulate ERAP2 expression and that disease-associated variants can convert a gene promoter region into a potent enhancer of a distal gene.

Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation

Background

The Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during CGN differentiation.

Results

We first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally-regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs.

Conclusion

Our data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

Position-specific evolution in transcription factor binding sites, and a fast likelihood calculation for the F81 model

Transcription factor binding sites (TFBS), like other DNA sequence, evolve via mutation and selection relating to their function. Models of nucleotide evolution describe DNA evolution via single-nucleotide mutation. A stationary vector of such a model is the long-term distribution of nucleotides, unchanging under the model. Neutrally evolving sites may have uniform stationary vectors, but one expects that sites within a TFBS instead have stationary vectors reflective of the fitness of various nucleotides at those positions. We introduce 'position-specific stationary vectors' (PSSVs), the collection of stationary vectors at each site in a TFBS locus, analogous to the position weight matrix (PWM) commonly used to describe TFBS. We infer PSSVs for human TFs using two evolutionary models (Felsenstein 1981 and Hasegawa-Kishino-Yano 1985). We find that PSSVs reflect the nucleotide distribution from PWMs, but with reduced specificity. We infer ancestral nucleotide distributions at individual positions and calculate 'conditional PSSVs' conditioned on specific choices of majority ancestral nucleotide. We find that certain ancestral nucleotides exert a strong evolutionary pressure on neighbouring sequence while others have a negligible effect. Finally, we present a fast likelihood calculation for the F81 model on moderate-sized trees that makes this approach feasible for large-scale studies along these lines.

Transcriptional Control of Subcutaneous Adipose Tissue by the Transcription Factor CTCF Modulates Heterogeneity in Fat Distribution in Women

Determining the mechanism driving body fat distribution will provide insights into obesity-related health risks. We used functional genomics tools to profile the epigenomic landscape to help infer the differential transcriptional potential of apple- and pear-shaped women's subcutaneous adipose-derived stem cells (ADSCs). We found that CCCTC-binding factor (CTCF) expression and its chromatin binding were increased in ADSCs from pear donors compared to those from apple donors. Interestingly, the pear enriched CTCF binding sites were located predominantly at the active transcription start sites (TSSs) of genes with active histone marks and YY1 motifs and were also associated with pear enriched RNAPII binding. In contrast, apple enriched CTCF binding sites were mainly found at intergenic regions and when identified at TSS, they were enriched with the bivalent chromatin signatures. Altogether, we provide evidence that CTCF plays an important role in differential regulation of subcutaneous ADSCs gene expression and may influence the development of apple vs. pear body shape.

Outward-oriented sites within clustered CTCF boundaries are key for intra-TAD chromatin interactions and gene regulation

CTCF plays an important role in 3D genome organization by adjusting the strength of chromatin insulation at TAD boundaries, where clustered CBS (CTCF-binding site) elements are often arranged in a tandem array with a complex divergent or convergent orientation. Here, using Pcdh and HOXD loci as a paradigm, we look into the clustered CTCF TAD boundaries and find that, counterintuitively, outward-oriented CBS elements are crucial for inward enhancer-promoter interactions as well as for gene regulation. Specifically, by combinatorial deletions of a series of putative enhancer elements in mice in vivo or CBS elements in cultured cells in vitro, in conjunction with chromosome conformation capture and RNA-seq analyses, we show that deletions of outward-oriented CBS elements weaken the strength of long-distance intra-TAD promoter-enhancer interactions and enhancer activation of target genes. Our data highlight the crucial role of outward-oriented CBS elements within the clustered CTCF TAD boundaries in developmental gene regulation and have interesting implications on the organization principles of clustered CTCF sites within TAD boundaries.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer

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

Association of DDX5/p68 protein with the upstream erythroid enhancer element (EHS1) of the gene encoding the KLF1 transcription factor

EKLF/KLF1 is an essential transcription factor that plays a global role in erythroid transcriptional activation. Regulation of KLF1 is of interest, as it displays a highly restricted expression pattern, limited to erythroid cells and its progenitors. Here we use biochemical affinity purification to identify the DDX5/p68 protein as an activator of KLF1 by virtue of its interaction with the erythroid-specific DNAse hypersensitive site upstream enhancer element (EHS1). We further show that this protein associates with DEK and CTCF. We postulate that the range of interactions of DDX5/p68 with these and other proteins known to interact with this element render it part of the enhanseosome complex critical for optimal expression of KLF1 and enables the formation of a proper chromatin configuration at the Klf1 locus. These individual interactions provide quantitative contributions that, in sum, establish the high-level activity of the Klf1 promoter and suggest they can be selectively manipulated for clinical benefit.

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

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

The Molecular Impacts of Retrotransposons in Development and Diseases

Retrotransposons are invasive genetic elements that constitute substantial portions of mammalian genomes. They have the potential to influence nearby gene expression through their cis-regulatory sequences, reverse transcription machinery, and the ability to mold higher-order chromatin structures. Due to their multifaceted functions, it is crucial for host fitness to maintain strict regulation of these parasitic sequences to ensure proper growth and development. This review explores how subsets of retrotransposons have undergone evolutionary exaptation to enhance the complexity of mammalian genomes. It also highlights the significance of regulating these elements, drawing on recent studies conducted in human and murine systems.

Abdominal surgery induces long-lasting changes in expression and binding of CTCF with impact on Major Histocompatibility Complex II transcription in circulating human monocytes

BACKGROUND

Postoperative immunosuppression has been recognized as an important driver of surgery-related morbidity and mortality. It is characterized by lymphocyte depression and impaired monocyte capability to present foreign antigens to T-cells via Major Histocompatibility Complex, Class II (MHC-II) molecules. In patients with postoperative abdominal sepsis, we previously detected a persisting differential binding of the CCCTC-Binding Factor (CTCF), a superordinate regulator of transcription, inside the MHC-II region with specific impact on human leucocyte antigen (HLA) gene expression. In this prospective exploratory study, we investigated to which extent major surgery affects the MHC-II region of circulating CD14+-monocytes.

RESULTS

In non-immunocompromised patients undergoing elective major abdominal surgery, a postoperative loss of monocyte HLA-DR surface receptor density was accompanied by a decline in the transcription levels of the classical MHC-II genes HLA-DRA, HLA-DRB1, HLA-DPA1 and HLA-DPB1. The surgical event decreased the expression of the transcriptional MHC-II regulators CIITA and CTCF and led to a lower CTCF enrichment at an intergenic sequence within the HLA-DR subregion. During the observation period, we found a slow and only incomplete restoration of monocyte HLA-DR surface receptor density as well as a partial recovery of CIITA, HLA-DRA and HLA-DRB1 expression. In contrast, transcription of HLA-DPA1, HLA-DPB1, CTCF and binding of CTCF within the MHC-II remained altered.

CONCLUSION

In circulating monocytes, major surgery does not globally affect MHC-II transcription but rather induces specific changes in the expression of selected HLA genes, followed by differential recovery patterns and accompanied by a prolonged reduction of CTCF expression and binding within the MHC-II region. Our results hint toward a long-lasting impact of a major surgical intervention on monocyte functionality, possibly mediated by epigenetic changes that endure the life span of the individual cell.

Cross-tissue patterns of DNA hypomethylation reveal genetically distinct histories of cell development

BACKGROUND

Establishment of DNA methylation (DNAme) patterns is essential for balanced multi-lineage cellular differentiation, but exactly how these patterns drive cellular phenotypes is unclear. While > 80% of CpG sites are stably methylated, tens of thousands of discrete CpG loci form hypomethylated regions (HMRs). Because they lack DNAme, HMRs are considered transcriptionally permissive, but not all HMRs actively regulate genes. Unlike promoter HMRs, a subset of non-coding HMRs is cell type-specific and enriched for tissue-specific gene regulatory functions. Our data further argues not only that HMR establishment is an important step in enforcing cell identity, but also that cross-cell type and spatial HMR patterns are functionally informative of gene regulation.

RESULTS

To understand the significance of non-coding HMRs, we systematically dissected HMR patterns across diverse human cell types and developmental timepoints, including embryonic, fetal, and adult tissues. Unsupervised clustering of 126,104 distinct HMRs revealed that levels of HMR specificity reflects a developmental hierarchy supported by enrichment of stage-specific transcription factors and gene ontologies. Using a pseudo-time course of development from embryonic stem cells to adult stem and mature hematopoietic cells, we find that most HMRs observed in differentiated cells (~ 60%) are established at early developmental stages and accumulate as development progresses. HMRs that arise during differentiation frequently (~ 35%) establish near existing HMRs (≤ 6 kb away), leading to the formation of HMR clusters associated with stronger enhancer activity. Using SNP-based partitioned heritability from GWAS summary statistics across diverse traits and clinical lab values, we discovered that genetic contribution to trait heritability is enriched within HMRs. Moreover, the contribution of heritability to cell-relevant traits increases with both increasing HMR specificity and HMR clustering, supporting the role of distinct HMR subsets in regulating normal cell function.

CONCLUSIONS

Our results demonstrate that the entire HMR repertoire within a cell-type, rather than just the cell type-specific HMRs, stores information that is key to understanding and predicting cellular phenotypes. Ultimately, these data provide novel insights into how DNA hypo-methylation provides genetically distinct historical records of a cell's journey through development, highlighting HMRs as functionally distinct from other epigenomic annotations.

Reprogramming of cis-regulatory networks during skeletal muscle atrophy in male mice

A comprehensive atlas of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of cellular structure maintenance, metabolism, and responses to the environment. Here we show, using matched single-nucleus chromatin accessibility and RNA-sequencing from juvenile male C57BL6 mice, an atlas of accessible chromatin regions in both normal and denervated skeletal muscles. We identified cell-type-specific cis-regulatory networks, highlighting the dynamic regulatory circuits mediating transitions between myonuclear types. Through comparison of normal and perturbed muscle, we delineated the reprogramming of cis-regulatory networks in response to denervation, described the interplay of promoters/enhancers and target genes. We further unveil a hierarchical structure of transcription factors that delineate a regulatory network in atrophic muscle, identifying ELK4 as a key atrophy-related transcription factor that instigates muscle atrophy through TGF-β1 regulation. This study furnishes a rich genomic resource, essential for decoding the regulatory dynamics of skeletal muscle in both physiological and pathological states.

Epigenetic balance ensures mechanistic control of MLL amplification and rearrangement

MLL/KMT2A amplifications and translocations are prevalent in infant, adult, and therapy-induced leukemia. However, the molecular contributor(s) to these alterations are unclear. Here, we demonstrate that histone H3 lysine 9 mono- and di-methylation (H3K9me1/2) balance at the MLL/KMT2A locus regulates these amplifications and rearrangements. This balance is controlled by the crosstalk between lysine demethylase KDM3B and methyltransferase G9a/EHMT2. KDM3B depletion increases H3K9me1/2 levels and reduces CTCF occupancy at the MLL/KMT2A locus, in turn promoting amplification and rearrangements. Depleting CTCF is also sufficient to generate these focal alterations. Furthermore, the chemotherapy doxorubicin (Dox), which associates with therapy-induced leukemia and promotes MLL/KMT2A amplifications and rearrangements, suppresses KDM3B and CTCF protein levels. KDM3B and CTCF overexpression rescues Dox-induced MLL/KMT2A alterations. G9a inhibition in human cells or mice also suppresses MLL/KMT2A events accompanying Dox treatment. Therefore, MLL/KMT2A amplifications and rearrangements are controlled by epigenetic regulators that are tractable drug targets, which has clinical implications.

LncRNA EILA promotes CDK4/6 inhibitor resistance in breast cancer by stabilizing cyclin E1 protein

CDK4/6 inhibitors (CDK4/6i) plus endocrine therapy are now standard first-line therapy for advanced HR+/HER2- breast cancer, but developing resistance is just a matter of time in these patients. Here, we report that a cyclin E1-interacting lncRNA (EILA) is up-regulated in CDK4/6i-resistant breast cancer cells and contributes to CDK4/6i resistance by stabilizing cyclin E1 protein. EILA overexpression correlates with accelerated cell cycle progression and poor prognosis in breast cancer. Silencing EILA reduces cyclin E1 protein and restores CDK4/6i sensitivity both in vitro and in vivo. Mechanistically, hairpin A of EILA binds to the carboxyl terminus of cyclin E1 protein and hinders its binding to FBXW7, thereby blocking its ubiquitination and degradation. EILA is transcriptionally regulated by CTCF/CDK8/TFII-I complexes and can be inhibited by CDK8 inhibitors. This study unveils the role of EILA in regulating cyclin E1 stability and CDK4/6i resistance, which may serve as a biomarker to predict therapy response and a potential therapeutic target to overcome resistance.

Deciphering a shared transcriptomic regulation and the relative contribution of each regulator type through endometrial gene expression signatures

BACKGORUND

While various endometrial biomarkers have been characterized at the transcriptomic and functional level, there is generally a poor overlap among studies, making it unclear to what extent their upstream regulators (e.g., ovarian hormones, transcription factors (TFs) and microRNAs (miRNAs)) realistically contribute to menstrual cycle progression and function. Unmasking the intricacies of the molecular interactions in the endometrium from a novel systemic point of view will help gain a more accurate perspective of endometrial regulation and a better explanation the molecular etiology of endometrial-factor infertility.

METHODS

An in-silico analysis was carried out to identify which regulators consistently target the gene biomarkers proposed in studies related to endometrial progression and implantation failure (19 gene lists/signatures were included). The roles of these regulators, and of genes related to progesterone and estrogens, were then analysed in transcriptomic datasets compiled from samples collected throughout the menstrual cycle (n = 129), and the expression of selected TFs were prospectively validated in an independent cohort of healthy participants (n = 19).

RESULTS

A total of 3,608 distinct genes from the 19 gene lists were associated with endometrial progression and implantation failure. The lists' regulation was significantly favoured by TFs (89% (17/19) of gene lists) and progesterone (47% (8 /19) of gene lists), rather than miRNAs (5% (1/19) of gene lists) or estrogen (0% (0/19) of gene lists), respectively (FDR < 0.05). Exceptionally, two gene lists that were previously associated with implantation failure and unexplained infertility were less hormone-dependent, but primarily regulated by estrogen. Although endometrial progression genes were mainly targeted by hormones rather than non-hormonal contributors (odds ratio = 91.94, FDR < 0.05), we identified 311 TFs and 595 miRNAs not previously associated with ovarian hormones. We highlight CTCF, GATA6, hsa-miR-15a-5p, hsa-miR-218-5p, hsa-miR-107, hsa-miR-103a-3p, and hsa-miR-128-3p, as overlapping novel master regulators of endometrial function. The gene expression changes of selected regulators throughout the menstrual cycle (FDR < 0.05), dually validated in-silico and through endometrial biopsies, corroborated their potential regulatory roles in the endometrium.

CONCLUSIONS

This study revealed novel hormonal and non-hormonal regulators and their relative contributions to endometrial progression and pathology, providing new leads for the potential causes of endometrial-factor infertility.

LncRNA ARGI Contributes to Virus-Induced Pancreatic β Cell Inflammation Through Transcriptional Activation of IFN-Stimulated Genes

Type 1 diabetes (T1D) is a complex autoimmune disease that develops in genetically susceptible individuals. Most T1D-associated single nucleotide polymorphisms (SNPs) are located in non-coding regions of the human genome. Interestingly, SNPs in long non-coding RNAs (lncRNAs) may result in the disruption of their secondary structure, affecting their function, and in turn, the expression of potentially pathogenic pathways. In the present work, the function of a virus-induced T1D-associated lncRNA named ARGI (Antiviral Response Gene Inducer) is characterized. Upon a viral insult, ARGI is upregulated in the nuclei of pancreatic β cells and binds to CTCF to interact with the promoter and enhancer regions of IFNβ and interferon-stimulated genes, promoting their transcriptional activation in an allele-specific manner. The presence of the T1D risk allele in ARGI induces a change in its secondary structure. Interestingly, the T1D risk genotype induces hyperactivation of type I IFN response in pancreatic β cells, an expression signature that is present in the pancreas of T1D patients. These data shed light on the molecular mechanisms by which T1D-related SNPs in lncRNAs influence pathogenesis at the pancreatic β cell level and opens the door for the development of therapeutic strategies based on lncRNA modulation to delay or avoid pancreatic β cell inflammation in T1D.

Methylation across the central dogma in health and diseases: new therapeutic strategies

The proper transfer of genetic information from DNA to RNA to protein is essential for cell-fate control, development, and health. Methylation of DNA, RNAs, histones, and non-histone proteins is a reversible post-synthesis modification that finetunes gene expression and function in diverse physiological processes. Aberrant methylation caused by genetic mutations or environmental stimuli promotes various diseases and accelerates aging, necessitating the development of therapies to correct the disease-driver methylation imbalance. In this Review, we summarize the operating system of methylation across the central dogma, which includes writers, erasers, readers, and reader-independent outputs. We then discuss how dysregulation of the system contributes to neurological disorders, cancer, and aging. Current small-molecule compounds that target the modifiers show modest success in certain cancers. The methylome-wide action and lack of specificity lead to undesirable biological effects and cytotoxicity, limiting their therapeutic application, especially for diseases with a monogenic cause or different directions of methylation changes. Emerging tools capable of site-specific methylation manipulation hold great promise to solve this dilemma. With the refinement of delivery vehicles, these new tools are well positioned to advance the basic research and clinical translation of the methylation field.

A novel SATB1 protein isoform with different biophysical properties

Intra-thymic T cell development is coordinated by the regulatory actions of SATB1 genome organizer. In this report, we show that SATB1 is involved in the regulation of transcription and splicing, both of which displayed deregulation in Satb1 knockout murine thymocytes. More importantly, we characterized a novel SATB1 protein isoform and described its distinct biophysical behavior, implicating potential functional differences compared to the commonly studied isoform. SATB1 utilized its prion-like domains to transition through liquid-like states to aggregated structures. This behavior was dependent on protein concentration as well as phosphorylation and interaction with nuclear RNA. Notably, the long SATB1 isoform was more prone to aggregate following phase separation. Thus, the tight regulation of SATB1 isoforms expression levels alongside with protein post-translational modifications, are imperative for SATB1's mode of action in T cell development. Our data indicate that deregulation of these processes may also be linked to disorders such as cancer.

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

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

Gene Interactions in Human Disease Studies-Evidence Is Mounting

Despite monumental advances in molecular technology to generate genome sequence data at scale, there is still a considerable proportion of heritability in most complex diseases that remains unexplained. Because many of the discoveries have been single-nucleotide variants with small to moderate effects on disease, the functional implication of many of the variants is still unknown and, thus, we have limited new drug targets and therapeutics. We, and many others, posit that one primary factor that has limited our ability to identify novel drug targets from genome-wide association studies may be due to gene interactions (epistasis), gene-environment interactions, network/pathway effects, or multiomic relationships. We propose that many of these complex models explain much of the underlying genetic architecture of complex disease. In this review, we discuss the evidence from multiple research avenues, ranging from pairs of alleles to multiomic integration studies and pharmacogenomics, that supports the need for further investigation of gene interactions (or epistasis) in genetic and genomic studies of human disease. Our goal is to catalog the mounting evidence for epistasis in genetic studies and the connections between genetic interactions and human health and disease that could enable precision medicine of the future.

Three-dimensional genome structure and function

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

Epigenomic analysis of Alzheimer's disease brains reveals diminished CTCF binding on genes involved in synaptic organization

Epigenetic aberrations are suggested to play an important role in transcriptional alterations in Alzheimer's disease (AD). One of the key mechanisms of epigenetic regulation of gene expression is through the dynamic organization of chromatin structure via the master genome architecture protein, CCCTC-binding factor (CTCF). By forming chromatin loops, CTCF can influence gene transcription in a complex manner. To find out whether genome-wide DNA binding sites for CTCF are altered in AD, we compared CTCF chromatin immunoprecipitation sequencing (ChIP-Seq) data from frontal cortex of human AD patients and normal controls (n = 9 pairs, all females). We have revealed that CTCF-binding affinity on many genes is significantly reduced in AD patients, and these genes are enriched in synaptic organization, cell adhesion, and actin cytoskeleton, including synaptic scaffolding molecules and receptors, such as SHANK2, HOMER1, NRXN1, CNTNAP2 and GRIN2A, and protocadherin (PCDH) and cadherin (CDH) family members. By comparing transcriptomic data from AD patients, we have discovered that many of the synaptic and adhesion genes with reduced CTCF binding in AD are significantly reduced in their mRNA expression. Moreover, a significant overlap of genes with the diminished CTCF binding and the reduced H3K27ac is identified in AD, with the common genes enriched in synaptic organization. These data suggest that the CTCF-controlled 3D chromatin organization is perturbed in AD, which may be linked to the diminished expression of target genes, probably through changes in histone modification.

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.