CTCF: master weaver of the genome

Three-dimensional loop extrusion

Loop extrusion convincingly describes how certain structural maintenance of chromosome (SMC) proteins mediate the formation of large DNA loops. Yet most of the existing computational models cannot reconcile recent in vitro observations showing that condensins can traverse each other, bypass large roadblocks, and perform steps longer than their own size. To fill this gap, we propose a three-dimensional (3D) "trans-grabbing" model for loop extrusion, which not only reproduces the experimental features of loop extrusion by one SMC complex but also predicts the formation of so-called Z-loops via the interaction of two or more SMCs extruding along the same DNA substrate. By performing molecular dynamics simulations of this model, we discover that the experimentally observed asymmetry in the different types of Z-loops is a natural consequence of the DNA tethering in vitro. Intriguingly, our model predicts this bias to disappear in the absence of tethering and a third type of Z-loop, which has not yet been identified in experiments, to appear. Our model naturally explains roadblock bypassing and the appearance of steps larger than the SMC size as a consequence of non-contiguous DNA grabbing. Finally, this study is the first, to our knowledge, to address how Z-loops and bypassing might occur in a way that is broadly consistent with existing cis-only 1D loop extrusion models.

Quantitative-enhancer-FACS-seq (QeFS) reveals epistatic interactions among motifs within transcriptional enhancers in developing Drosophila tissue

Understanding the contributions of transcription factor DNA binding sites to transcriptional enhancers is a significant challenge. We developed Quantitative enhancer-FACS-Seq for highly parallel quantification of enhancer activities from a genomically integrated reporter in Drosophila melanogaster embryos. We investigate the contributions of the DNA binding motifs of four poorly characterized TFs to the activities of twelve embryonic mesodermal enhancers. We measure quantitative changes in enhancer activity and discover a range of epistatic interactions among the motifs, both synergistic and alleviating. We find that understanding the regulatory consequences of TF binding motifs requires that they be investigated in combination across enhancer contexts.

A deletion containing a CTCF-element in intron 8 of the Bbs7 gene is partially responsible for juvenile obesity in the Berlin Fat Mouse

The Berlin Fat Mouse Inbred (BFMI) line is a model for juvenile obesity. Previous studies on crosses between BFMI and C57Bl/6N (B6N) have identified a recessive defect causing juvenile obesity on chromosome 3 (jObes1). Bbs7 was identified as the most likely candidate gene for the observed effect. Comparative sequence analysis showed a 1578 bp deletion in intron 8 of Bbs7 in BFMI mice. A CTCF-element is located inside this deletion. To investigate the functional effect of this deletion, it was introduced into B6N mice using CRISPR/Cas9. Two mice containing the target deletion were obtained (B6N Bbs7^emI8∆1 and Bbs7^emI8∆2) and were subsequently mated to BFMI and B6N to generate two families suitable for complementation. Inherited alleles were determined and body composition was measured by quantitative magnetic resonance. Evidence for a partial complementation (13.1-15.1%) of the jObes1 allele by the CRISPR/Cas9 modified B6N Bbs7^emI8∆1 and Bbs7^emI8∆2 alleles was found. Mice carrying the complementation alleles had a 23-27% higher fat-to-lean ratio compared to animals which have a B6N allele (P_{(Bbs7emI8∆1)} = 4.25 × 10^-7; P_{(Bbs7emI8∆2)} = 3.17 × 10^-5). Consistent with previous findings, the recessive effect of the BFMI allele was also seen for the B6N Bbs7^emI8∆1 and Bbs7^emI8∆2 alleles. However, the effect size of the B6N Bbs7^emI8∆1 and Bbs7^emI8∆2 alleles was smaller than the BFMI allele, and thus showed only a partial complementation. Findings suggest additional variants near Bbs7 in addition to or interacting with the deletion in intron 8.

CCCTC-binding factor is an upstream regulator of the pluripotency factor Oct4 and functions in active transcription of linc1253 and linc1356 genes in pluripotent cells

The embryonic stem cell- (ESC) specific transcription factor Oct4 is a well-known master regulator of pluripotency. The aim of this study was to identify upstream regulators of Oct4 and explore their functional link in regulating lincRNA expression in ESCs. By quantitative real-time PCR (RT-qPCR) analysis upon CCCTC-binding factor (CTCF) or Oct4 knockdown, here, we found that the chromatin insulator CTCF transcriptionally controls Oct4 gene expression in mouse ES cells. Furthermore, co-immunoprecipitation assays showed that CTCF physically interacts with Oct4. By analyzing CTCF and Oct4 ChIP-seq datasets in mouse ES cells and investigating their genomic occupancies, we demonstrated that CTCF and Oct4 share overlapping regulatory functions and are required for active transcription of long intergenic non-coding RNAs (lincRNAs) linc1253 and linc1356, which were reported to repress cellular lineage programs and maintain a pluripotent state. In summary, we propose an integrated model of transcriptional control mediated by CTCF, the master weaver of the genome, for the upstream regulation of Oct4-and ESC-associated genes. These results connect the chromatin insulator CTCF and the pluripotency factor Oct4 in the regulation of lincRNAs in pluripotent cells.

KDM6B promotes activation of the oncogenic CDK4/6-pRB-E2F pathway by maintaining enhancer activity in MYCN-amplified neuroblastoma

The H3K27me2/me3 histone demethylase KDM6B is essential to neuroblastoma cell survival. However, the mechanism of KDM6B action remains poorly defined. We demonstrate that inhibition of KDM6B activity 1) reduces the chromatin accessibility of E2F target genes and MYCN, 2) selectively leads to an increase of H3K27me3 but a decrease of the enhancer mark H3K4me1 at the CTCF and BORIS binding sites, which may, consequently, disrupt the long-range chromatin interaction of MYCN and E2F target genes, and 3) phenocopies the transcriptome induced by the specific CDK4/6 inhibitor palbociclib. Overexpression of CDK4/6 or Rb1 knockout confers neuroblastoma cell resistance to both palbociclib and the KDM6 inhibitor GSK-J4. These data indicate that KDM6B promotes an oncogenic CDK4/6-pRB-E2F pathway in neuroblastoma cells via H3K27me3-dependent enhancer-promoter interactions, providing a rationale to target KDM6B for high-risk neuroblastoma.

The histone demethylase KDM6B is reported to be essential for neuroblastoma cell survival. Here the authors show that KDM6B regulates CDK4/6-pRB-E2F pathway through H3K27me3-dependent enhancer-promoter interactions in neuroblastoma.

A translational exploration of the effects of WNT2 variants on altered cortical structures in autism spectrum disorder

BACKGROUND

Evidence suggests that cortical anatomy may be aytpical in autism spectrum disorder. The wingless-type MMTV integration site family, member 2 (WNT2), a candidate gene for autism spectrum disorder, may regulate cortical development. However, it is unclear whether WNT2 variants are associated with altered cortical thickness in autism spectrum disorder.

METHODS

In a sample of 118 people with autism spectrum disorder and 122 typically developing controls, we investigated cortical thickness using FreeSurfer software. We then examined the main effects of the WNT2 variants and the interactions of group × SNP and age × SNP for each hemisphere and brain region that was altered in people with autism spectrum disorder.

RESULTS

Compared to neurotypical controls, people with autism spectrum disorder showed reduced mean cortical thickness in both hemispheres and 9 cortical regions after false discovery rate correction, including the right cingulate gyrus, the orbital gyrus, the insula, the inferior frontal gyrus (orbital part and triangular part), the lateral occipitotemporal gyrus, the posterior transverse collateral sulcus, the lateral sulcus and the superior temporal sulcus. In the full sample, 2 SNPs of WNT2 (rs6950765 and rs2896218) showed age × SNP interactions for the mean cortical thickness of both hemispheres, the middle-posterior cingulate cortex and the superior temporal cortex.

LIMITATIONS

We examined the genetic effect for each hemisphere and the 9 regions that were altered in autism spectrum disorder. The age effect we found in this cross-sectional study needs to be examined in longitudinal studies.

CONCLUSION

Based on neuroimaging and genetic data, our findings suggest that WNT2 variants might be associated with altered cortical thickness in autism spectrum disorder. Whether and how these WNT2 variants might involve cortical thinning requires further investigation.

TRIAL REGISTRATION

ClinicalTrials.gov no. NCT01582256.

PROTOCOL REGISTRATION

National Institutes of Health no. NCT00494754.

Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.

Downregulation of BORIS/CTCFL leads to ROS-dependent cellular senescence and drug sensitivity in MYCN-amplified neuroblastoma

Brother of Regulator of Imprinted Sites (BORIS) or CCCTC-binding factor like (CTCFL) is a nucleotide-binding protein, aberrantly expressed in various malignancies. Expression of BORIS has been found to be associated with the expression of oncogenes which regulate the reactive oxygen species (ROS) biogenesis, DNA double-strand break repair, regulation of stemness, and induction of cellular senescence. In the present study, we have analyzed the effects of knockdown of BORIS, a potential oncogene, on the induction of senescence and tumor suppression. Loss of BORIS downregulated the expression of critical oncogenes such as BMI1, Akt, MYCN, and STAT3, whereas overexpression increased their respective expression levels in MYCN-amplified neuroblastoma cells. BORIS knockdown exhibited high levels of ROS biogenesis, indicating an upregulated mitochondrial superoxide production and thereby induction of senescence. Our study also showed that the loss of BORIS facilitated cellular senescence through the disruption of telomere integrity via altering the expression of various proteins required for telomere capping (POT1, TRF2, and TIN1). In addition to affecting ROS production and DNA damage, BORIS knockdown sensitized the cells toward chemotherapeutic drugs and induced apoptosis. Tumor induction studies on in vivo xenograft mouse models showed that cells with loss of BORIS/CTCFL failed to induce tumors. From our study, we conclude that silencing BORIS/CTCFL influences tumor growth and proliferation by regulating key oncogenes. The results also indicated that the BORIS knockdown can cause cellular senescence and upon a combinatorial treatment with chemotherapeutic drugs can induce enhanced drug sensitivity in MYCN-amplified neuroblastoma cells.

CircESRP1 inhibits clear cell renal cell carcinoma progression through the CTCF-mediated positive feedback loop

Circular RNA (circRNA), a closed continuous loop formed by back-splicing, has been confirmed to be implicated in a variety of human diseases including cancers. However, the underlying molecular mechanism of circRNA regulating the progression of renal cell carcinoma (RCC) remains largely unclear. In the present study, we identified a novel circular RNA, circESRP1, that derived from the ESRP1 gene locus at 8q22.1 exons. Lower expression of circESRP1 was found in clear cell RCC (ccRCC) tissues and cell lines. Besides, circESRP1 expression level showed inversely correlated with the advanced tumor size, TNM stage and distant metastasis of ccRCC. The expression level of circESRP1 exhibited a positive correlation with CTCF protein but negatively correlated with miR-3942 in 79 ccRCC tissues. In vivo experiments, we found that overexpression of circESRP1 effectively repressed xenograft tumor growth and inhibited c-Myc-mediated EMT progression. CircESRP1 acted as a sponge to competitively bind with miR-3942 as confirmed through RNA pull-down, RIP and dual-luciferase reporter assays. Moreover, CTCF, a downstream target of miR-3942, was validated to specifically promote the circESRP1 transcript expression and regulated by circESRP1/miR-3942 pathway to form a positive feedback loop. We also revealed that the circESRP1/miR-3942/CTCF feedback loop regulated the ccRCC cell functions via c-Myc mediated EMT process. This study provides a novel regulatory model of circRNA via forming a positive-feedback loop that perpetuates the circESRP1/miR-3942/CTCF axis, suggesting that this signaling may serve as a novel target for the treatment of ccRCC.

Transcription Factor CTCFL Promotes Cell Proliferation, Migration, and Invasion in Gastric Cancer via Activating DPPA2

Objective

To explore the relationship between CTCFL and DPPA2 and validate the positive role of CTCFL/DPPA2 in cell malignant behaviors in gastric cancer.

Methods

We predicted gastric cancer-related transcription factors and corresponding target mRNAs through bioinformatics. Levels of CTCFL and DPPA2 were assessed via qRT-PCR and western blot. In vitro experiments were utilized to assay the cell biological behaviors. CHIP was utilized for the assessment of the targeted relationship between CTCFL and DPPA2.

Results

CTCFL and DPPA2 were both highly expressed in gastric cancer cells, and high CTCFLL and DPPA2 could promote cell malignant behaviors. CHIP validated that DPPA2 was a target of CTCFL. In addition, high DPPA2 rescued the repressive impact of CTCFL silencing on the cell proliferation, migration, and invasion in gastric cancer.

Conclusion

The transcription factor CTCFL fosters cell proliferative, migratory, and invasive properties via activating DPPA2 in gastric cancer.

Hexavalent chromium disrupts chromatin architecture

Accessibility of DNA elements and the orchestration of spatiotemporal chromatin-chromatin interactions are critical mechanisms in the regulation of gene transcription. Thus, in an ever-changing milieu, cells mount an adaptive response to environmental stimuli by modulating gene expression that is orchestrated by coordinated changes in chromatin architecture. Correspondingly, agents that alter chromatin structure directly impact transcriptional programs in cells. Heavy metals, including hexavalent chromium (Cr(VI)), are of special interest because of their ability to interact directly with cellular protein, DNA and other macromolecules, resulting in general damage or altered function. In this review we highlight the chromium-mediated mechanisms that promote disruption of chromatin architecture and how these processes are integral to its carcinogenic properties. Emerging evidence shows that Cr(VI) targets nucleosomal architecture in euchromatin, particularly in genomic locations flanking binding sites of the essential transcription factors CTCF and AP1. Ultimately, these changes contribute to an altered chromatin state in critical gene regulatory regions, which disrupts gene transcription in functionally relevant biological processes.

CGGBP1-dependent CTCF-binding sites restrict ectopic transcription

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

Transposable Element Dynamics and Regulation during Zygotic Genome Activation in Mammalian Embryos and Embryonic Stem Cell Model Systems

Transposable elements (TEs) are mobile genetic sequences capable of duplicating and reintegrating at new regions within the genome. A growing body of evidence has demonstrated that these elements play important roles in host genome evolution, despite being traditionally viewed as parasitic elements. To prevent ectopic activation of TE transposition and transcription, they are epigenetically silenced in most somatic tissues. Intriguingly, a specific class of TEs-retrotransposons-is transiently expressed at discrete phases during mammalian development and has been linked to the establishment of totipotency during zygotic genome activation (ZGA). While mechanisms controlling TE regulation in somatic tissues have been extensively studied, the significance underlying the unique transcriptional reactivation of retrotransposons during ZGA is only beginning to be uncovered. In this review, we summarize the expression dynamics of key retrotransposons during ZGA, focusing on findings from in vivo totipotent embryos and in vitro totipotent-like embryonic stem cells (ESCs). We then dissect the functions of retrotransposons and discuss how their transcriptional activities are finetuned during early stages of mammalian development.

Resilience in Long-Term Viral Infection: Genetic Determinants and Interactions

Virus-induced neurological sequelae resulting from infection by Theiler's murine encephalomyelitis virus (TMEV) are used for studying human conditions ranging from epileptic seizures to demyelinating disease. Mouse strains are typically considered susceptible or resistant to TMEV infection based on viral persistence and extreme phenotypes, such as demyelination. We have identified a broader spectrum of phenotypic outcomes by infecting strains of the genetically diverse Collaborative Cross (CC) mouse resource. We evaluated the chronic-infection gene expression profiles of hippocampi and thoracic spinal cords for 19 CC strains in relation to phenotypic severity and TMEV persistence. Strains were clustered based on similar phenotypic profiles and TMEV levels at 90 days post-infection, and we categorized distinct TMEV response profiles. The three most common profiles included "resistant" and "susceptible," as before, as well as a "resilient" TMEV response group which experienced both TMEV persistence and mild neurological phenotypes even at 90 days post-infection. Each profile had a distinct gene expression signature, allowing the identification of pathways and networks specific to each TMEV response group. CC founder haplotypes for genes involved in these pathways/networks revealed candidate response-specific alleles. These alleles demonstrated pleiotropy and epigenetic (miRNA) regulation in long-term TMEV infection, with particular relevance for resilient mouse strains.

Recent Issues in Varicella-Zoster Virus Latency

Varicella-zoster virus (VZV) is a human herpes virus which causes varicella (chicken pox) as a primary infection, and, following a variable period of latency in neurons in the peripheral ganglia, may reactivate to cause herpes zoster (shingles) as well as a variety of neurological syndromes. In this overview we consider some recent issues in alphaherpesvirus latency with special focus on VZV ganglionic latency. A key question is the nature and extent of viral gene transcription during viral latency. While it is known that this is highly restricted, it is only recently that the very high degree of that restriction has been clarified, with both VZV gene 63-encoded transcripts and discovery of a novel VZV transcript (VLT) that maps antisense to the viral transactivator gene 61. It has also emerged in recent years that there is significant epigenetic regulation of VZV gene transcription, and the mechanisms underlying this are complex and being unraveled. The last few years has also seen an increased interest in the immunological aspects of VZV latency and reactivation, in particular from the perspective of inborn errors of host immunity that predispose to different VZV reactivation syndromes.

An Epigenetic Perspective on Intra-Tumour Heterogeneity: Novel Insights and New Challenges from Multiple Fields

Simple Summary

Although research on cancer biology in recent decades has unveiled the main genetic perturbations driving the onset of tumorigenesis, we are still far from properly treating this disease without the occurrence of drug resistance and metastatic burden. This achievement is hampered by the onset of intra-tumour heterogeneity (ITH), which increases cancer cell fitness and plasticity, thereby fostering cell adaptation to foreign environments and stimuli. In this review, we discuss the contribution of the epigenetic factors in sustaining ITH and their interplay with the tumour microenvironment. We also highlight the recent technological advancements that are contributing to defining the epigenetic mechanisms governing tumour heterogeneity at the single-cell level.

Abstract

Cancer is a group of heterogeneous diseases that results from the occurrence of genetic alterations combined with epigenetic changes and environmental stimuli that increase cancer cell plasticity. Indeed, multiple cancer cell populations coexist within the same tumour, favouring cancer progression and metastatic dissemination as well as drug resistance, thereby representing a major obstacle for treatment. Epigenetic changes contribute to the onset of intra-tumour heterogeneity (ITH) as they facilitate cell adaptation to perturbation of the tumour microenvironment. Despite being its central role, the intrinsic multi-layered and reversible epigenetic pattern limits the possibility to uniquely determine its contribution to ITH. In this review, we first describe the major epigenetic mechanisms involved in tumourigenesis and then discuss how single-cell-based approaches contribute to dissecting the key role of epigenetic changes in tumour heterogeneity. Furthermore, we highlight the importance of dissecting the interplay between genetics, epigenetics, and tumour microenvironments to decipher the molecular mechanisms governing tumour progression and drug resistance.

SATB1-dependent mitochondrial ROS production controls TCR signaling in CD4 T cells

SATB1 regulates mitochondrial function and reactive oxygen species (ROS) production through the expression of mitochondrial transcription factor A. SATB1-mediated ROS production is necessary for TCR stimulation and T-cell function.

Special AT-rich sequence binding protein-1 (SATB1) is localized to the nucleus and remodels chromatin structure in T cells. SATB1-deficient CD4 T cells cannot respond to TCR stimulation; however, the cause of this unresponsiveness is to be clarified. Here, we demonstrate that SATB1 is indispensable to proper mitochondrial functioning and necessary for the activation of signal cascades via the TCR in CD4 T cells. Naïve SATB1-deficient CD4 T cells contain fewer mitochondria than WT T cells, as the former do not express mitochondrial transcription factor A (TFAM). Impaired mitochondrial function in SATB1-deficient T cells subverts mitochondrial ROS production and SHP-1 inactivation by constitutive oxidization. Ectopic TFAM expression increases mitochondrial mass and mitochondrial ROS production and rescues defects in the antigen-specific response in the SATB1-deficient T cells. Thus, SATB1 is vital for maintaining mitochondrial mass and function by regulating TFAM expression, which is necessary for TCR signaling.

Stress-Induced Epstein-Barr Virus Reactivation

Epstein-Barr virus (EBV) is typically found in a latent, asymptomatic state in immunocompetent individuals. Perturbations of the host immune system can stimulate viral reactivation. Furthermore, there are a myriad of EBV-associated illnesses including various cancers, post-transplant lymphoproliferative disease, and autoimmune conditions. A thorough understanding of this virus, and the interplay between stress and the immune system, is essential to establish effective treatment. This review will provide a summary of the interaction between both psychological and cellular stressors resulting in EBV reactivation. It will examine mechanisms by which EBV establishes and maintains latency and will conclude with a brief overview of treatments targeting EBV.

Biallelic null variants in ZNF142 cause global developmental delay with familial epilepsy and dysmorphic features

Biallelic variants in ZNF142 at 2q35, which encodes zinc-finger protein 142, cause neurodevelopmental disorder with seizures or dystonia. We identified compound heterozygous null variants in ZNF142, NM_001105537.4:c.[1252C>T];[1274-2A>G],p.[Arg418*];[Glu426*], in Malaysian siblings suffering from global developmental delay with epilepsy and dysmorphism. cDNA analysis showed the marked reduction of ZNF142 transcript level through nonsense-mediated mRNA decay by these novel biallelic variants. The affected siblings present with global developmental delay and epilepsy in common, which were previously described, as well as dysmorphism, which was not recognized. It is important to collect patients with ZNF142 abnormality to define its phenotypic spectrum.

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes.

Long noncoding RNAs: Emerging regulators of normal and malignant hematopoiesis

Genome wide analyses have revealed that long-noncoding RNAs (lncRNAs) are not only passive transcription products, but also major regulators of genome structure and transcription. In particular, lncRNAs exert profound effects on various biological processes, such as chromatin structure, transcription, RNA stability and translation, and protein degradation and localization, which depend on their localization and interacting partners. Recent studies have revealed that thousands of lncRNAs are aberrantly expressed in various cancer types and some of them are associated with malignant transformation. Despite extensive efforts, the diverse functions of lncRNAs and molecular mechanisms in which they act remain elusive. Many hematological disorders and malignancies are primarily resulted from genetic alterations that lead to the dysregulation of gene regulatory networks required for cellular proliferation and differentiation. Consequently, a growing list of lncRNAs has been reported for their involvement in the modulation of hematopoietic gene expression networks and hematopoietic stem and progenitor cell (HS/PC) function. Dysregulation of some of these lncRNAs has been attributed to pathogenesis of hematological malignancies. In this review, we will summarize current advances and knowledge of lncRNAs in gene regulation, focusing on the recent progresses on the role of lncRNAs in CTCF/cohesin mediated three-dimensional (3D) genome organization, and how such genome folding signals in turn regulate transcription, HS/PC function and transformation. The knowledge will provide mechanistic and translational insights into HS/PC biology and myeloid malignancy pathophysiology.

A Novel Regulatory Player in the Innate Immune System: Long Non-Coding RNAs

Long non-coding RNAs (lncRNAs) represent crucial transcriptional and post-transcriptional gene regulators during antimicrobial responses in the host innate immune system. Studies have shown that lncRNAs are expressed in a highly tissue- and cell-specific- manner and are involved in the differentiation and function of innate immune cells, as well as inflammatory and antiviral processes, through versatile molecular mechanisms. These lncRNAs function via the interactions with DNA, RNA, or protein in either cis or trans pattern, relying on their specific sequences or their transcriptions and processing. The dysregulation of lncRNA function is associated with various human non-infectious diseases, such as inflammatory bowel disease, cardiovascular diseases, and diabetes mellitus. Here, we provide an overview of the regulation and mechanisms of lncRNA function in the development and differentiation of innate immune cells, and during the activation or repression of innate immune responses. These elucidations might be beneficial for the development of therapeutic strategies targeting inflammatory and innate immune-mediated diseases.

LAG3 (LAG-3, CD223) DNA methylation correlates with LAG3 expression by tumor and immune cells, immune cell infiltration, and overall survival in clear cell renal cell carcinoma

Background

Lymphocyte activating 3 (LAG3, LAG-3, CD223) is a promising target for immune checkpoint inhibition in clear cell renal cell carcinoma (KIRC). The aim of this study was to investigate the epigenetic regulation of LAG3 in KIRC by methylation.

Methods

We correlated quantitative LAG3 methylation levels with transcriptional activity, immune cell infiltration, and overall survival in a cohort of n=533 patients with KIRC and n=160 normal adjacent tissue (NAT) samples obtained from The Cancer Genome Atlas (TCGA). Furthermore, we analyzed LAG3 methylation in peripheral blood mononuclear cells (PBMCs) and KIRC cell lines. We validated correlations between LAG3 expression, immune cell infiltrates, survival, and methylation in an independent KIRC cohort (University Hospital Bonn (UHB) cohort, n=118) by means of immunohistochemistry and quantitative methylation-specific PCR.

Results

We found differential methylation profiles among PBMCs, NAT, KIRC cell lines, and KIRC tumor tissue. Methylation strongly correlated with LAG3 mRNA expression in KIRCs (TCGA cohort) and KIRC cell lines. In the UHB cohort, methylation correlated with LAG3-positive immune cells and tumor-intrinsic LAG3 protein expression. Furthermore, LAG3 methylation strongly correlated with signatures of distinct immune cell infiltrates, an interferon-y signature (TCGA cohort), and immunohistochemically quantified CD45⁺, CD8⁺, and CD4⁺ immune cell infiltrates (UHB cohort). LAG3 mRNA expression (TCGA cohort), methylation (both cohorts), and tumor cell-intrinsic protein expression (UHB cohort) was significantly associated with overall survival.

Conclusion

Our data suggest an epigenetic regulation of LAG3 expression in tumor and immune cells via DNA methylation. LAG3 expression and methylation is associated with a subset of KIRCs showing a distinct clinical course and immunogenicity. Our study provides rationale for further testing LAG3 DNA methylation as a predictive biomarker for response to LAG3 immune checkpoint inhibitors.

Acute depletion of CTCF rewires genome-wide chromatin accessibility

BACKGROUND

The transcription factor CTCF appears indispensable in defining topologically associated domain boundaries and maintaining chromatin loop structures within these domains, supported by numerous functional studies. However, acute depletion of CTCF globally reduces chromatin interactions but does not significantly alter transcription.

RESULTS

Here, we systematically integrate multi-omics data including ATAC-seq, RNA-seq, WGBS, Hi-C, Cut&Run, and CRISPR-Cas9 survival dropout screens, and time-solved deep proteomic and phosphoproteomic analyses in cells carrying auxin-induced degron at endogenous CTCF locus. Acute CTCF protein degradation markedly rewires genome-wide chromatin accessibility. Increased accessible chromatin regions are frequently located adjacent to CTCF-binding sites at promoter regions and insulator sites associated with enhanced transcription of nearby genes. In addition, we use CTCF-associated multi-omics data to establish a combinatorial data analysis pipeline to discover CTCF co-regulatory partners. We successfully identify 40 candidates, including multiple established partners. Interestingly, many CTCF co-regulators that have alterations of their respective downstream gene expression do not show changes of their own expression levels across the multi-omics measurements upon acute CTCF loss, highlighting the strength of our system to discover hidden co-regulatory partners associated with CTCF-mediated transcription.

CONCLUSIONS

This study highlights that CTCF loss rewires genome-wide chromatin accessibility, which plays a critical role in transcriptional regulation.

Dynamic landscape of chromatin accessibility and transcriptomic changes during differentiation of human embryonic stem cells into dopaminergic neurons

Chromatin architecture influences transcription by modulating the physical access of regulatory factors to DNA, playing fundamental roles in cell identity. Studies on dopaminergic differentiation have identified coding genes, but the relationship with non-coding genes or chromatin accessibility remains elusive. Using RNA-Seq and ATAC-Seq we profiled differentially expressed transcripts and open chromatin regions during early dopaminergic neuron differentiation. Hierarchical clustering of differentially expressed genes, resulted in 6 groups with unique characteristics. Surprisingly, the abundance of long non-coding RNAs (lncRNAs) was high in the most downregulated transcripts, and depicted positive correlations with target mRNAs. We observed that open chromatin regions decrease upon differentiation. Enrichment analyses of accessibility depict an association between open chromatin regions and specific functional pathways and gene-sets. A bioinformatic search for motifs allowed us to identify transcription factors and structural nuclear proteins that potentially regulate dopaminergic differentiation. Interestingly, we also found changes in protein and mRNA abundance of the CCCTC-binding factor, CTCF, which participates in genome organization and gene expression. Furthermore, assays demonstrated co-localization of CTCF with Polycomb-repressed chromatin marked by H3K27me3 in pluripotent cells, progressively decreasing in neural precursor cells and differentiated neurons. Our work provides a unique resource of transcription factors and regulatory elements, potentially involved in the acquisition of human dopaminergic neuron cell identity.

CTCF as a regulator of alternative splicing: new tricks for an old player

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the 'master weaver of the genome'. Given its multiple connections, further characterization of CTCF's emerging role in splicing regulation might extend its functional repertoire towards a 'conductor of the splicing orchestra'.

The Mutational Landscape of Myeloid Leukaemia in Down Syndrome

Children with Down syndrome (DS) are particularly prone to haematopoietic disorders. Paediatric myeloid malignancies in DS occur at an unusually high frequency and generally follow a well-defined stepwise clinical evolution. First, the acquisition of mutations in the GATA1 transcription factor gives rise to a transient myeloproliferative disorder (TMD) in DS newborns. While this condition spontaneously resolves in most cases, some clones can acquire additional mutations, which trigger myeloid leukaemia of Down syndrome (ML-DS). These secondary mutations are predominantly found in chromatin and epigenetic regulators-such as cohesin, CTCF or EZH2-and in signalling mediators of the JAK/STAT and RAS pathways. Most of them are also found in non-DS myeloid malignancies, albeit at extremely different frequencies. Intriguingly, mutations in proteins involved in the three-dimensional organization of the genome are found in nearly 50% of cases. How the resulting mutant proteins cooperate with trisomy 21 and mutant GATA1 to promote ML-DS is not fully understood. In this review, we summarize and discuss current knowledge about the sequential acquisition of genomic alterations in ML-DS.

DNA methylation QTL analysis identifies new regulators of human longevity

Human longevity is a complex trait influenced by both genetic and environmental factors, whose interaction is mediated by epigenetic mechanisms like DNA methylation. Here, we generated genome-wide whole-blood methylome data from 267 individuals, of which 71 were long-lived (90–104 years), by applying reduced representation bisulfite sequencing. We followed a stringent two-stage analysis procedure using discovery and replication samples to detect differentially methylated sites (DMSs) between young and long-lived study participants. Additionally, we performed a DNA methylation quantitative trait loci analysis to identify DMSs that underlie the longevity phenotype. We combined the DMSs results with gene expression data as an indicator of functional relevance. This approach yielded 21 new candidate genes, the majority of which are involved in neurophysiological processes or cancer. Notably, two candidates (PVRL2, ERCC1) are located on chromosome 19q, in close proximity to the well-known longevity- and Alzheimer’s disease-associated loci APOE and TOMM40. We propose this region as a longevity hub, operating on both a genetic (APOE, TOMM40) and an epigenetic (PVRL2, ERCC1) level. We hypothesize that the heritable methylation and associated gene expression changes reported here are overall advantageous for the LLI and may prevent/postpone age-related diseases and facilitate survival into very old age.

Multiple CTCF sites cooperate with each other to maintain a TAD for enhancer-promoter interaction in the β-globin locus

Insulators are cis-regulatory elements that block enhancer activity and prevent heterochromatin spreading. The binding of CCCTC-binding factor (CTCF) protein is essential for insulators to play the roles in a chromatin context. The β-globin locus, consisting of multiple genes and enhancers, is flanked by two insulators 3'HS1 and HS5. However, it has been reported that the absence of these insulators did not affect the β-globin transcription. To explain the unexpected finding, we have deleted a CTCF motif at 3'HS1 or HS5 in the human β-globin locus and analyzed chromatin interactions around the locus. It was found that a topologically associating domain (TAD) containing the β-globin locus is maintained by neighboring CTCF sites in the CTCF motif-deleted loci. The additional deletions of neighboring CTCF motifs disrupted the β-globin TAD, resulting in decrease of the β-globin transcription. Chromatin interactions of the β-globin enhancers with gene promoter were weakened in the multiple CTCF motifs-deleted loci, even though the enhancers have still active chromatin features such as histone H3K27ac and histone H3 depletion. Genome-wide analysis using public CTCF ChIA-PET and ChIP-seq data showed that chromatin domains possessing multiple CTCF binding sites tend to contain super-enhancers like the β-globin enhancers. Taken together, our results show that multiple CTCF sites surrounding the β-globin locus cooperate with each other to maintain a TAD. The β-globin TAD appears to provide a compact spatial environment that enables enhancers to interact with promoter.

Three-dimensional chromosome organization in flowering plants

Research on plant three-dimensional (3D) genome architecture made rapid progress over the past 5 years. Numerous Hi-C interaction data sets were generated in a wide range of plant species, allowing for a comprehensive overview on 3D chromosome folding principles in the plant kingdom. Plants lack important genes reported to be vital for chromosome folding in animals. However, similar 3D structures such as topologically associating domains and chromatin loops were identified. Recent studies in Arabidopsis thaliana revealed how chromosomal regions are positioned within the nucleus by determining their association with both, the nuclear periphery and the nucleolus. Additionally, many plant species exhibit high-frequency interactions among KNOT entangled elements, which are associated with safeguarding the genome from invasive DNA elements. Many of the recently published Hi-C data sets were generated to aid de novo genome assembly and remain to date little explored. These data sets represent a valuable resource for future comparative studies, which may lead to a more profound understanding of the evolution of 3D chromosome organization in plants.

Order and disorder: abnormal 3D chromatin organization in human disease

A precise three-dimensional (3D) organization of chromatin is central to achieve the intricate transcriptional patterns that are required to form complex organisms. Growing evidence supports an important role of 3D chromatin architecture in development and delineates its alterations as prominent causes of disease. In this review, we discuss emerging concepts on the fundamental forces shaping genomes in space and on how their disruption can lead to pathogenic phenotypes. We describe the molecular mechanisms underlying a wide range of diseases, from the systemic effects of coding mutations on 3D architectural factors, to the more tissue-specific phenotypes resulting from genetic and epigenetic modifications at specific loci. Understanding the connection between the 3D organization of the genome and its underlying biological function will allow a better interpretation of human pathogenesis.

A universal framework for detecting cis-regulatory diversity in DNA regulatory regions

High-throughput sequencing-based assays measure different biochemical activities pertaining to gene regulation, genome-wide. These activities include transcription factor (TF)–DNA binding, enhancer activity, open chromatin, and more. A major goal is to understand underlying sequence components, or motifs, that can explain the measured activity. It is usually not one motif but a combination of motifs bound by cooperatively acting proteins that confers activity to such regions. Furthermore, regions can be diverse, governed by different combinations of TFs/motifs. Current approaches do not take into account this issue of combinatorial diversity. We present a new statistical framework, cisDIVERSITY, which models regions as diverse modules characterized by combinations of motifs while simultaneously learning the motifs themselves. Because cisDIVERSITY does not rely on knowledge of motifs, modules, cell type, or organism, it is general enough to be applied to regions reported by most high-throughput assays. For example, in enhancer predictions resulting from different assays—GRO-cap, STARR-seq, and those measuring chromatin structure—cisDIVERSITY discovers distinct modules and combinations of TF binding sites, some specific to the assay. From protein–DNA binding data, cisDIVERSITY identifies potential cofactors of the profiled TF, whereas from ATAC-seq data, it identifies tissue-specific regulatory modules. Finally, analysis of single-cell ATAC-seq data suggests that regions open in one cell-state encode information about future states, with certain modules staying open and others closing down in the next time point.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Molecular Complexes at Euchromatin, Heterochromatin and Centromeric Chromatin

Chromatin consists of a complex of DNA and histone proteins as its core components and plays an important role in both packaging DNA and regulating DNA metabolic pathways such as DNA replication, transcription, recombination, and chromosome segregation. Proper functioning of chromatin further involves a network of interactions among molecular complexes that modify chromatin structure and organization to affect the accessibility of DNA to transcription factors leading to the activation or repression of the transcription of target DNA loci. Based on its structure and compaction state, chromatin is categorized into euchromatin, heterochromatin, and centromeric chromatin. In this review, we discuss distinct chromatin factors and molecular complexes that constitute euchromatin-open chromatin structure associated with active transcription; heterochromatin-less accessible chromatin associated with silencing; centromeric chromatin-the site of spindle binding in chromosome segregation.

The combined action of CTCF and its testis-specific paralog BORIS is essential for spermatogenesis

CTCF is a key organizer of the 3D genome. Its specialized paralog, BORIS, heterodimerizes with CTCF but is expressed only in male germ cells and in cancer states. Unexpectedly, BORIS-null mice have only minimal germ cell defects. To understand the CTCF-BORIS relationship, mouse models with varied CTCF and BORIS levels were generated. Whereas Ctcf+/+Boris+/+, Ctcf+/-Boris+/+, and Ctcf+/+Boris-/- males are fertile, Ctcf+/-Boris-/- (Compound Mutant; CM) males are sterile. Testes with combined depletion of both CTCF and BORIS show reduced size, defective meiotic recombination, increased apoptosis, and malformed spermatozoa. Although CM germ cells exhibit only 25% of CTCF WT expression, chromatin binding of CTCF is preferentially lost from CTCF-BORIS heterodimeric sites. Furthermore, CM testes lose the expression of a large number of spermatogenesis genes and gain the expression of developmentally inappropriate genes that are "toxic" to fertility. Thus, a combined action of CTCF and BORIS is required to both repress pre-meiotic genes and activate post-meiotic genes for a complete spermatogenesis program.

SATB1-mediated chromatin landscape in T cells

The regulatory circuits that define developmental decisions of thymocytes are still incompletely resolved. SATB1 protein is predominantly expressed at the CD4⁺CD8⁺cell stage exerting its broad transcription regulation potential with both activatory and repressive roles. A series of post-translational modifications and the presence of potential SATB1 protein isoforms indicate the complexity of its regulatory potential. The most apparent mechanism of its involvement in gene expression regulation is via the orchestration of long-range chromatin loops between genes and their regulatory elements. Multiple SATB1 perturbations in mice uncovered a link to autoimmune diseases while clinical investigations on cancer research uncovered that SATB1 has a promoting role in several types of cancer and can be used as a prognostic biomarker. SATB1 is a multivalent tissue-specific factor with a broad and yet undetermined regulatory potential. Future investigations on this protein could further uncover T cell-specific regulatory pathways and link them to (patho)physiology.

GAD1 expression and its methylation as indicators of malignant behavior in thymic epithelial tumors

Thymic epithelial tumors (TETs) comprise thymomas and thymic carcinoma (TC). TC has more aggressive features and a poorer prognosis than thymomas. Genetic and epigenetic alterations in thymomas and TC have been investigated in an attempt to identify novel target molecules for TC. In the present study, genome-wide screening was performed on aberrantly methylated CpG islands in thymomas and TC, and the glutamate decarboxylase 1 gene (GAD1) was identified as the 4th significantly hypermethylated CpG island in TC compared with thymomas. GAD1 catalyzes the production of γ-aminobutyric acid from L-glutamic acid. GAD1 expression is abundant in the brain but rare in other tissues, including the thymus. A total of 73 thymomas and 17 TC tissues were obtained from 90 patients who underwent surgery or biopsy at Tokushima University Hospital between 1990 and 2017. DNA methylation was examined by bisulfite pyrosequencing, and the mRNA and protein expression levels of GAD1 were analyzed using reverse transcription-quantitative PCR and immunohistochemistry, respectively. The DNA methylation levels of GAD1 were significantly higher in TC tissues than in the normal thymus and thymoma tissues, and GAD1 methylation exhibited high sensitivity and specificity for discriminating between TC and thymoma. The mRNA and protein expression levels of GAD1 were significantly higher in TC tissues than in thymomas. Patients with TET with high GAD1 DNA hypermethylation and high mRNA and protein expression levels had significantly shorter relapse-free survival rates than those with low levels. In conclusion, significantly more epigenetic alterations were observed in TC tissues compared with in thymomas, which may contribute to the clinical features and prognosis of patients.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Regulatory Noncoding and Predicted Pathogenic Coding Variants of CCR5 Predispose to Severe COVID-19

Genome-wide association studies (GWAS) found locus 3p21.31 associated with severe COVID-19. CCR5 resides at the same locus and, given its known biological role in other infection diseases, we investigated if common noncoding and rare coding variants, affecting CCR5, can predispose to severe COVID-19. We combined single nucleotide polymorphisms (SNPs) that met the suggestive significance level (P ≤ 1 × 10-5) at the 3p21.31 locus in public GWAS datasets (6406 COVID-19 hospitalized patients and 902,088 controls) with gene expression data from 208 lung tissues, Hi-C, and Chip-seq data. Through whole exome sequencing (WES), we explored rare coding variants in 147 severe COVID-19 patients. We identified three SNPs (rs9845542, rs12639314, and rs35951367) associated with severe COVID-19 whose risk alleles correlated with low CCR5 expression in lung tissues. The rs35951367 resided in a CTFC binding site that interacts with CCR5 gene in lung tissues and was confirmed to be associated with severe COVID-19 in two independent datasets. We also identified a rare coding variant (rs34418657) associated with the risk of developing severe COVID-19. Our results suggest a biological role of CCR5 in the progression of COVID-19 as common and rare genetic variants can increase the risk of developing severe COVID-19 by affecting the functions of CCR5.

Mechanism for DPY30 and ASH2L intrinsically disordered regions to modulate the MLL/SET1 activity on chromatin

Recent cryo-EM structures show the highly dynamic nature of the MLL1-NCP (nucleosome core particle) interaction. Functional implication and regulation of such dynamics remain unclear. Here we show that DPY30 and the intrinsically disordered regions (IDRs) of ASH2L work together in restricting the rotational dynamics of the MLL1 complex on the NCP. We show that DPY30 binding to ASH2L leads to stabilization and integration of ASH2L IDRs into the MLL1 complex and establishes new ASH2L-NCP contacts. The significance of ASH2L-DPY30 interactions is demonstrated by requirement of both ASH2L IDRs and DPY30 for dramatic increase of processivity and activity of the MLL1 complex. This DPY30 and ASH2L-IDR dependent regulation is NCP-specific and applies to all members of the MLL/SET1 family of enzymes. We further show that DPY30 is causal for de novo establishment of H3K4me3 in ESCs. Our study provides a paradigm of how H3K4me3 is regulated on chromatin and how H3K4me3 heterogeneity can be modulated by ASH2L IDR interacting proteins.

Regulation of the MLL family of histone H3K4 methyltransferases on the nucleosome core particle (NCP) remains largely unknown. Here the authors show that intrinsically disordered regions of ASH2L and DPY30 restrict the rotational dynamics of MLL1 on the NCP, allowing more efficient enzyme-substrate engagement and higher H3K4 trimethylation activity.

Discovering unknown human and mouse transcription factor binding sites and their characteristics from ChIP-seq data

Transcription factor binding sites (TFBSs) are essential for gene regulation, but the number of known TFBSs remains limited. We aimed to discover and characterize unknown TFBSs by developing a computational pipeline for analyzing ChIP-seq (chromatin immunoprecipitation followed by sequencing) data. Applying it to the latest ENCODE ChIP-seq data for human and mouse, we found that using the irreproducible discovery rate as a quality-control criterion resulted in many experiments being unnecessarily discarded. By contrast, the number of motif occurrences in ChIP-seq peak regions provides a highly effective criterion, which is reliable even if supported by only one experimental replicate. In total, we obtained 2,058 motifs from 1,089 experiments for 354 human TFs and 163 motifs from 101 experiments for 34 mouse TFs. Among these motifs, 487 have not previously been reported. Mapping the canonical motifs to the human genome reveals a high TFBS density ±2 kb around transcription start sites (TSSs) with a peak at -50 bp. On average, a promoter contains 5.7 TFBSs. However, 70% of TFBSs are in introns (41%) and intergenic regions (29%), whereas only 12% are in promoters (-1 kb to +100 bp from TSSs). Notably, some TFs (e.g., CTCF, JUN, JUNB, and NFE2) have motifs enriched in intergenic regions, including enhancers. We inferred 142 cobinding TF pairs and 186 (including 115 completely) tethered binding TF pairs, indicating frequent interactions between TFs and a higher frequency of tethered binding than cobinding. This study provides a large number of previously undocumented motifs and insights into the biological and genomic features of TFBSs.

Cohesin subunit Rad21 binds to the HSV-1 genome near CTCF insulator sites during latency in vivo

Herpes Simplex Virus 1 (HSV-1) is a human pathogen that has the ability to establish a lifelong infection in the host. During latency, HSV-1 genomes are chromatinized and are abundantly associated with histones in sensory neurons, yet the mechanisms that govern the latent-lytic transition remain unclear. We hypothesize that the latent-lytic switch is controlled by CTCF insulators, positioned within the HSV-1 latent genome. CTCF insulators, together with the cohesin complex, have the ability to establish and maintain chromtin loops that allow distance separated gene regions to be spatially oriented for transcriptional control. In this current study, we demonstrated that the cohesin subunit Rad21 was recruited to latent HSV-1 genomes near four of the CTCF insulators during latency. We showed that the CTCF insulator known as CTRS1/2, positioned downstream from the essential transactivating IE region of ICP4 was only enriched in Rad21 prior to but not during latency, suggesting that the CTRS1/2 insulator is not required for the maintenance of latency. Further, deletion of the CTRL2 insulator, positioned downstream from the LAT enhancer, resulted in a loss of Rad21 enrichment at insulators flanking the ICP4 region at early times post-infection in mice ganglia, suggesting that these insulators are interdependent. Finally, deletion of the CTRL2 insulator resulted in a loss of Rad21 enrichment at the CTRL2 insulator in a cell-type specific manner, and this loss of Rad21 enrichment was correlated to decreased LAT expression, suggesting that Rad21 recruitment to viral genomes is important for efficient gene expression.ImportanceCTCF insulators are important for transcriptional control and increasing evidence suggests that that CTCF insulators, together with the cohesin complex, regulate viral transcription in DNA viruses. The CTCF-cohesin interaction is important for the formation of chromatin loops, structures that orient distance separated elements in close spatial proximity for transcriptional control. Herpes Simplex Virus 1 (HSV-1) has seven putative CTCF insulators that flank the LAT and the IE, indicating that CTCF insulators play a role in the transition from latency to reactivation. Contributions from the work presented here include the finding that CTCF insulators in HSV-1 genomes are differentially enriched in the cohesin subunit Rad21, suggesting that CTCF-cohesin interactions could be establishing and anchoring chromatin loop structures to control viral transcription.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Postoperative abdominal sepsis induces selective and persistent changes in CTCF binding within the MHC-II region of human monocytes

Background

Postoperative abdominal infections belong to the most common triggers of sepsis and septic shock in intensive care units worldwide. While monocytes play a central role in mediating the initial host response to infections, sepsis-induced immune dysregulation is characterized by a defective antigen presentation to T-cells via loss of Major Histocompatibility Complex Class II DR (HLA-DR) surface expression. Here, we hypothesized a sepsis-induced differential occupancy of the CCCTC-Binding Factor (CTCF), an architectural protein and superordinate regulator of transcription, inside the Major Histocompatibility Complex Class II (MHC-II) region in patients with postoperative sepsis, contributing to an altered monocytic transcriptional response during critical illness.

Results

Compared to a matched surgical control cohort, postoperative sepsis was associated with selective and enduring increase in CTCF binding within the MHC-II. In detail, increased CTCF binding was detected at four sites adjacent to classical HLA class II genes coding for proteins expressed on monocyte surface. Gene expression analysis revealed a sepsis-associated decreased transcription of (i) the classical HLA genes HLA-DRA, HLA-DRB1, HLA-DPA1 and HLA-DPB1 and (ii) the gene of the MHC-II master regulator, CIITA (Class II Major Histocompatibility Complex Transactivator). Increased CTCF binding persisted in all sepsis patients, while transcriptional recovery CIITA was exclusively found in long-term survivors.

Conclusion

Our experiments demonstrate differential and persisting alterations of CTCF occupancy within the MHC-II, accompanied by selective changes in the expression of spatially related HLA class II genes, indicating an important role of CTCF in modulating the transcriptional response of immunocompromised human monocytes during critical illness.

Regulation of host and viral promoters during human cytomegalovirus latency via US28 and CTCF

Viral latency is an active process during which the host cell environment is optimized for latent carriage and reactivation. This requires control of both viral and host gene promoters and enhancers often at the level of chromatin, and several viruses co-opt the chromatin organiser CTCF to control gene expression during latency. While CTCF has a role in the latencies of alpha- and gamma-herpesviruses, it was not known whether CTCF played a role in the latency of the beta-herpesvirus human cytomegalovirus (HCMV). Here, we show that HCMV latency is associated with increased CTCF expression and CTCF binding to the viral major lytic promoter, the major immediate early promoter (MIEP). This increase in CTCF binding is dependent on the virally encoded G protein coupled receptor, US28, and contributes to suppression of MIEP-driven transcription, a hallmark of latency. Furthermore, we show that latency-associated upregulation of CTCF represses expression of the neutrophil chemoattractants S100A8 and S100A9 which we have previously shown are downregulated during HCMV latency. As with downregulation of the MIEP, CTCF binding to the enhancer region of S100A8/A9 drives their suppression, again in a US28-dependent manner. Taken together, we identify CTCF upregulation as an important mechanism for optimizing latent carriage of HCMV at both the levels of viral and cellular gene expression.

Association between structural brain features and gene expression by weighted gene co-expression network analysis in conversion from MCI to AD

Alzheimer's disease (AD) is a neurodegenerative disease. Mild cognitive impairment (MCI) represents a state of cognitive function between normal cognition and dementia. Longitudinal studies showed that some MCI patients remained in a state of MCI, and some developed AD. The reason for these different conversions from MCI remains to be investigated. 180 MCI participants were followed for eight years. 143 MCI patients maintained the MCI state (MCI_S), and the remaining thirty-seven MCI patients were re-evaluated as having AD (MCI_AD). We obtained 1,036 structural brain characteristics and 15,481 gene expression values from the 180 MCI participants and applied weighted gene co-expression network analysis (WGCNA) to explore the relationship between structural brain features and gene expression. Regulating mediator effect analysis was employed to explore the relationships among gene expression, brain region measurements and clinical phenotypes. We found that 60 genes from the MCI_S group and 18 genes from the MCI_AD group respectively had the most significant correlations with left paracentral lobule and sulcus (L.PTS) and right subparietal sulcus (R.SubPS) thickness; CTCF, UQCR11 and WDR5B were the mutual genes between the two groups. The expression of CTCF gene and clinical score are completely mediated by L.PTS thickness, and the UQCR11 and WDR5B gene expression levels significantly regulate the mediating effect pathway. In conclusion, the factors affecting the different conversions from MCI are closely related to L.PTS thickness and the CTCF, UQCR11 and WDR5B gene expression levels. Our results add a theoretical foundation of imaging genetics for conversion from MCI to AD.

Deletion of CTCF sites in the SHH locus alters enhancer-promoter interactions and leads to acheiropodia

Acheiropodia, congenital limb truncation, is associated with homozygous deletions in the LMBR1 gene around ZRS, an enhancer regulating SHH during limb development. How these deletions lead to this phenotype is unknown. Using whole-genome sequencing, we fine-mapped the acheiropodia-associated region to 12 kb and show that it does not function as an enhancer. CTCF and RAD21 ChIP-seq together with 4C-seq and DNA FISH identify three CTCF sites within the acheiropodia-deleted region that mediate the interaction between the ZRS and the SHH promoter. This interaction is substituted with other CTCF sites centromeric to the ZRS in the disease state. Mouse knockouts of the orthologous 12 kb sequence have no apparent abnormalities, showcasing the challenges in modelling CTCF alterations in animal models due to inherent motif differences between species. Our results show that alterations in CTCF motifs can lead to a Mendelian condition due to altered enhancer-promoter interactions.