TDP-43 chronic deficiency leads to dysregulation of transposable elements and gene expression by affecting R-loop and 5hmC crosstalk

TDP-43 is an RNA/DNA-binding protein that forms aggregates in various brain disorders. TDP-43 engages in many aspects of RNA metabolism, but its molecular roles in regulating genes and transposable elements (TEs) have not been extensively explored. Chronic TDP-43 knockdown impairs cell proliferation and cellular responses to DNA damage. At the molecular level, TDP-43 chronic deficiency affects gene expression either locally or distally by concomitantly altering the crosstalk between R-loops and 5-hydroxymethylcytosine (5hmC) in gene bodies and long-range enhancer/promoter interactions. Furthermore, TDP-43 knockdown induces substantial disease-relevant TE activation by influencing their R-loop and 5hmC homeostasis in a locus-specific manner. Together, our findings highlight the genomic roles of TDP-43 in modulating R-loop-5hmC coordination in coding genes, distal regulatory elements, and TEs, presenting a general and broad molecular mechanism underlying the contributions of proteinopathies to the etiology of neurodegenerative disorders.

Ectopic transcription due to inappropriately inherited histone methylation may interfere with the ongoing function of terminally differentiated cells

Many human neurodevelopmental disorders are caused by de novo mutations in histone modifying enzymes. These patients have craniofacial defects, developmental delay, intellectual disability and behavioral abnormalities, but it remains unclear how the mutations lead to such developmental defects. Here we take advantage of the invariant C. elegans lineage along with a unique double mutant in the H3K4me1/2 demethylase SPR-5/LSD1/KDM1A and the H3K9 methyltransferase MET-2/SETDB1 to address this question. We demonstrate that spr-5; met-2 double mutant worms have a severe chemotaxis defect that is dependent upon the ectopic expression of germline genes in somatic tissues. In addition, by performing single-cell RNAseq, we find that germline genes begin to be ectopically expression widely in spr-5; met-2 embryos. However, surprisingly we found that spr-5; met-2 mutants have no somatic lineage defects prior to the 200-cell stage of embryogenesis. This suggests that the altered chemotaxis behavior may be due to ongoing defect in terminally differentiated cells rather than a defect in development. To test this directly, we used RNAi to shut off the ectopic expression of germline genes in L2 spr-5; met-2 larvae, which have a fully formed nervous system. Remarkably, we find that shutting off the ectopic germline expression rescues normal chemotaxis behavior in the same adult worms that previously had a chemotaxis defect at the L2 stage. This suggests that ongoing ectopic transcription can block normal behavior in a fully intact nervous system. These data raise the possibility that intellectual disability and altered behavior in neurodevelopmental syndromes, caused by mutations in histone modifying enzymes, could be due to ongoing ectopic transcription and may be reversible.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Structures of CTCF-DNA complexes including all 11 zinc fingers

The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12-15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1-ZF7 and ZF3-ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1-ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9-ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1-ZF7 and ZF8-ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1-ZF7 and bulkier and polar/charged residues for ZF8-ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine-guanine and glutamine-adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9-ZF11. The differences between ZF1-ZF7 and ZF8-ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.

Characterization of a strain-specific CD-1 reference genome reveals potential inter- and intra-strain functional variability

BACKGROUND

CD-1 is an outbred mouse stock that is frequently used in toxicology, pharmacology, and fundamental biomedical research. Although inbred strains are typically better suited for such studies due to minimal genetic variability, outbred stocks confer practical advantages over inbred strains, such as improved breeding performance and low cost. Knowledge of the full genetic variability of CD-1 would make it more useful in toxicology, pharmacology, and fundamental biomedical research.

RESULTS

We performed deep genomic DNA sequencing of CD-1 mice and used the data to identify genome-wide SNPs, indels, and germline transposable elements relative to the mm10 reference genome. We used multiple genome-wide sequencing data types and previously published CD-1 SNPs to validate our called variants. We used the called variants to construct a strain-specific CD-1 reference genome, which we show can improve mappability and reduce experimental biases from genome-wide sequencing data derived from CD-1 mice. Based on previously published ChIP-seq and ATAC-seq data, we find evidence that genetic variation between CD-1 mice can lead to alterations in transcription factor binding. We also identified a number of variants in the coding region of genes which could have effects on translation of genes.

CONCLUSIONS

We have identified millions of previously unidentified CD-1 variants with the potential to confound studies involving CD-1. We used the identified variants to construct a CD-1-specific reference genome, which can improve accuracy and reduce bias when aligning genomics data derived from CD-1 mice.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

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

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

Single nucleus multiome analysis of the prefrontal cortex from C9orf72 ALS/FTD patients illuminates pathways affected during disease progression

Repeat expansions in the C9orf72 gene are the most common genetic cause of amyotrophic lateral sclerosis and familial frontotemporal dementia (ALS/FTD). To identify molecular defects that take place in the dorsolateral frontal cortex of patients with C9orf72 ALS/FTD, we compared healthy controls with C9orf72 ALS/FTD donor samples staged based on the levels of cortical phosphorylated TAR DNA binding protein (pTDP-43), a neuropathological hallmark of disease progression. We identified distinct molecular changes in different cell types that take place during disease progression. These alterations include downregulation of nuclear and mitochondrial ribosomal protein genes in early disease stages that become upregulated as the disease progresses. High ratios of premature oligodendrocytes expressing low levels of genes encoding major myelin protein components are characteristic of late disease stages and may represent a unique signature of C9orf72 ALS/FTD. Microglia with increased reactivity and astrocyte specific transcriptome changes in genes involved in glucose/glycogen metabolism are also associated with disease progression. Late stages of C9orf72 ALS/FTD correlate with sequential changes in the regulatory landscape of several genes in glial cells, namely MBP/MAG/MOG in oligodendrocytes, CD83/IRF8 in microglia, and GLUT1/GYS2/AGL in astrocytes. Only layer 2-3 cortical projection neurons with high expression of CUX2/LAMP5 are significantly reduced in C9orf72 ALS/FTD patients with respect to controls. Our findings reveal previously unknown progressive functional changes in cortical cells of C9orf72 ALS/FTD patients that shed light on the mechanisms underlying the pathology of this disease.

Activation of HIV-1 proviruses increases downstream chromatin accessibility

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

Recruitment of CTCF to an Fto enhancer is responsible for transgenerational inheritance of BPA-induced obesity

The mechanisms by which environmentally-induced epiphenotypes are transmitted transgenerationally in mammals are poorly understood. Here we show that exposure of pregnant mouse females to bisphenol A (BPA) results in obesity in the F2 progeny due to increased food intake. This epiphenotype can be transmitted up to the F6 generation. Analysis of chromatin accessibility in sperm of the F1-F6 generations reveals alterations at sites containing binding motifs for CCCTC-binding factor (CTCF) at two cis-regulatory elements (CREs) of the Fto gene that correlate with transmission of obesity. These CREs show increased interactions in sperm of obese mice with the Irx3 and Irx5 genes, which are involved in the differentiation of appetite-controlling neurons. Deletion of the CTCF site in Fto results in mice that have normal food intake and fail to become obese when ancestrally exposed to BPA. The results suggest that epigenetic alterations of Fto can lead to the same phenotypes as genetic variants.

Beyond Genes: Germline Disruption in the Etiology of Autism Spectrum Disorders

Investigations into the etiology of autism spectrum disorders have been largely confined to two realms: variations in DNA sequence and somatic developmental exposures. Here we suggest a third route-disruption of the germline epigenome induced by exogenous toxicants during a parent's gamete development. Similar to cases of germline mutation, these molecular perturbations may produce dysregulated transcription of brain-related genes during fetal and early development, resulting in abnormal neurobehavioral phenotypes in offspring. Many types of exposures may have these impacts, and here we discuss examples of anesthetic gases, tobacco components, synthetic steroids, and valproic acid. Alterations in parental germline could help explain some unsolved phenomena of autism, including increased prevalence, missing heritability, skewed sex ratio, and heterogeneity of neurobiology and behavior.

The Cupid shuffle: Do enhancers prefer specific promoters?

Two recent reports (Martinez-Ara et al., 2022; Bergman et al., 2022) explore the compatibility between enhancers and promoters and find that enhancers preferentially activate promoters with low intrinsic activity rather than favoring housekeeping or cell-type-specific promoters.

Introduction of Mutant GNAQ into Endothelial Cells Induces a Vascular Malformation Phenotype with Therapeutic Response to Imatinib

Simple Summary

Mutations in GNAQ underlie vascular malformations, including Sturge-Weber disease. In order to develop novel therapies for lesions with mutant GNAQ, we introduced mutant GNAQ into MS1 endothelial cells. Mutant GNAQ conferred a novel phenotype of progressive vascular malformations in mice. Chromatin analysis revealed upregulation of C-Kit in the vascular endothelial cells, and we found C-Kit to be highly expressed in Sturge-Weber disease. Given that imatinib is an FDA approved multikinase inhibitor that blocks C-Kit, we evaluated it in our mouse model, and showed that imatinib had activity against these vascular malformations. Repurposing imatinib should be evaluated in clinical trials, including Sturge-Weber disease.

Abstract

GNAQ is mutated in vascular and melanocytic lesions, including vascular malformations and nevi. No in vivo model of GNAQ activation in endothelial cells has previously been described. We introduce mutant GNAQ into a murine endothelial cell line, MS1. The resultant transduced cells exhibit a novel phenotype in vivo, with extensive vasoformative endothelial cells forming aberrant lumens similar to those seen in vascular malformations. ATAC-seq analysis reveals activation of c-Kit in the novel vascular malformations. We demonstrate that c-Kit is expressed in authentic human Sturge–Weber vascular malformations, indicating a novel druggable target for Sturge–Weber syndrome. Since c-Kit is targeted by the FDA-approved drug imatinib, we tested the ability of imatinib on the phenotype of the vascular malformations in vivo. Imatinib treated vascular malformations are significantly smaller and have decreased supporting stromal cells surrounding the lumen. Imatinib may be useful in the treatment of human vascular malformations that express c-Kit, including Sturge–Weber syndrome.

Is developmental synchrony enabled by CTCF residence time?

Depletion of CTCF in cultured cells has minor effects on transcription whereas its mutation leads to embryonic lethality and developmental defects. In a recent issue of Nature Cell Biology, Soochit et al. (2021) show that the residence time of CTCF on DNA may explain its critical role in cell differentiation.

Principles of 3D compartmentalization of the human genome

Chromatin is organized in the nucleus via CTCF loops and compartmental domains. Here, we compare different cell types to identify distinct paradigms of compartmental domain formation in human tissues. We identify and quantify compartmental forces correlated with histone modifications characteristic of transcriptional activity and previously underappreciated roles for distinct compartmental domains correlated with the presence of H3K27me3 and H3K9me3, respectively. We present a computer simulation model capable of predicting compartmental organization based on the biochemical characteristics of independent chromatin features. Using this model, we show that the underlying forces responsible for compartmental domain formation in human cells are conserved and that the diverse compartmentalization patterns seen across cell types are due to differences in chromatin features. We extend these findings to Drosophila to suggest that the same principles are at work beyond humans. These results offer mechanistic insights into the fundamental forces driving the 3D organization of the genome.

Using high-resolution Hi-C data and computer simulations, Nichols and Corces show that compartments arise as a consequence of interactions among proteins that correlate with the presence of H3K27ac, H3K27me3, and H3K9me3, suggesting that human cells contain at least three distinct compartments. The same principles apply to other organisms.

Exposure to sevoflurane results in changes of transcription factor occupancy in sperm and inheritance of autism

One in 54 children in the U.S. is diagnosed with Autism Spectrum Disorder (ASD). De novo germline and somatic mutations cannot account for all cases of ASD, suggesting that epigenetic alterations triggered by environmental exposures may be responsible for a subset of ASD cases. Human and animal studies have shown that exposure of the developing brain to general anesthetic (GA) agents can trigger neurodegeneration and neurobehavioral abnormalities but the effects of general anesthetics on the germ line have not been explored in detail. We exposed pregnant mice to sevoflurane during the time of embryonic development when the germ cells undergo epigenetic reprogramming and found that more than 38% of the directly exposed F1 animals exhibit impairments in anxiety and social interactions. Strikingly, 44-47% of the F2 and F3 animals, which were not directly exposed to sevoflurane, show the same behavioral problems. We performed ATAC-seq and identified more than 1200 differentially accessible sites in the sperm of F1 animals, 69 of which are also present in the sperm of F2 animals. These sites are located in regulatory regions of genes strongly associated with ASD, including Arid1b, Ntrk2, and Stmn2. These findings suggest that epimutations caused by exposing germ cells to sevoflurane can lead to ASD in the offspring, and this effect can be transmitted through the male germline inter and trans-generationally.

Replication timing maintains the global epigenetic state in human cells

The temporal order of DNA replication [replication timing (RT)] is correlated with chromatin modifications and three-dimensional genome architecture; however, causal links have not been established, largely because of an inability to manipulate the global RT program. We show that loss of RIF1 causes near-complete elimination of the RT program by increasing heterogeneity between individual cells. RT changes are coupled with widespread alterations in chromatin modifications and genome compartmentalization. Conditional depletion of RIF1 causes replication-dependent disruption of histone modifications and alterations in genome architecture. These effects were magnified with successive cycles of altered RT. These results support models in which the timing of chromatin replication and thus assembly plays a key role in maintaining the global epigenetic state.

Regulation of 3D chromatin organization by CTCF

Studies of nuclear architecture using chromosome conformation capture methods have provided a detailed view of how chromatin folds in the 3D nuclear space. New variants of this technology now afford unprecedented resolution and allow the identification of ever smaller folding domains that offer new insights into the mechanisms by which this organization is established and maintained. Here we review recent results in this rapidly evolving field with an emphasis on CTCF function, with the goal of gaining a mechanistic understanding of the principles by which chromatin is folded in the eukaryotic nucleus.

Haspin kinase modulates nuclear architecture and Polycomb-dependent gene silencing

Haspin, a highly conserved kinase in eukaryotes, has been shown to be responsible for phosphorylation of histone H3 at threonine 3 (H3T3ph) during mitosis, in mammals and yeast. Here we report that haspin is the kinase that phosphorylates H3T3 in Drosophila melanogaster and it is involved in sister chromatid cohesion during mitosis. Our data reveal that haspin also phosphorylates H3T3 in interphase. H3T3ph localizes in broad silenced domains at heterochromatin and lamin-enriched euchromatic regions. Loss of haspin compromises insulator activity in enhancer-blocking assays and triggers a decrease in nuclear size that is accompanied by changes in nuclear envelope morphology. We show that haspin is a suppressor of position-effect variegation involved in heterochromatin organization. Our results also demonstrate that haspin is necessary for pairing-sensitive silencing and it is required for robust Polycomb-dependent homeotic gene silencing. Haspin associates with the cohesin complex in interphase, mediates Pds5 binding to chromatin and cooperates with Pds5-cohesin to modify Polycomb-dependent homeotic transformations. Therefore, this study uncovers an unanticipated role for haspin kinase in genome organization of interphase cells and demonstrates that haspin is required for homeotic gene regulation.

Protection from DNA re-methylation by transcription factors in primordial germ cells and pre-implantation embryos can explain trans-generational epigenetic inheritance

Background

A growing body of evidence suggests that certain epiphenotypes can be passed across generations via both the male and female germlines of mammals. These observations have been difficult to explain owing to a global loss of the majority of known epigenetic marks present in parental chromosomes during primordial germ cell development and after fertilization.

Results

By integrating previously published BS-seq, DNase-seq, ATAC-seq, and RNA-seq data collected during multiple stages of primordial germ cell and pre-implantation development, we find that the methylation status of the majority of CpGs genome-wide is restored after global de-methylation, despite the fact that global CpG methylation drops to 10% in primordial germ cells and 20% in the inner cell mass of the blastocyst. We estimate the proportion of such CpGs with preserved methylation status to be 78%. Further, we find that CpGs at sites bound by transcription factors during the global re-methylation phases of germline and embryonic development remain hypomethylated across all developmental stages observed. On the other hand, CpGs at sites not bound by transcription factors during the global re-methylation phase have high methylation levels prior to global de-methylation, become de-methylated during global de-methylation, and then become re-methylated.

Conclusions

The results suggest that transcription factors can act as carriers of epigenetic information during germ cell and pre-implantation development by ensuring that the methylation status of CpGs is maintained. These findings provide the basis for a mechanistic description of trans-generational inheritance of epigenetic information in mammals.

Author Correction: Widespread long-range cis-regulatory elements in the maize genome

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mini-review: Epigenetic mechanisms that promote transgenerational actions of endocrine disrupting chemicals: Applications to behavioral neuroendocrinology

It is our hope this mini-review will stimulate discussion and new research. Here we briefly examine the literature on transgenerational actions of endocrine disrupting chemicals (EDCs) on brain and behavior and their underlying epigenetic mechanisms including: DNA methylation, histone modifications, and non-coding RNAs. We stress that epigenetic modifications need to be examined in a synergistic manner, as they act together in situ on chromatin to change transcription. Next we highlight recent work from one of our laboratories (VGC). The data provide new evidence that the sperm genome is poised for transcription. In developing sperm, gene enhancers and promoters are accessible for transcription and these activating motifs are also found in preimplantation embryos. Thus, DNA modifications associated with transcription factors during fertilization, in primordial germ cells (PGCs), and/or during germ cell maturation may be passed to offspring. We discuss the implications of this model to EDC exposures and speculate on whether natural variation in hormone levels during fertilization and PGC migration may impart transgenerational effects on brain and behavior. Lastly we discuss how this mechanism could apply to neural sexual differentiation.

Analysis of Hi-C data using SIP effectively identifies loops in organisms from C. elegans to mammals

Chromatin loops are a major component of 3D nuclear organization, visually apparent as intense point-to-point interactions in Hi-C maps. Identification of these loops is a critical part of most Hi-C analyses. However, current methods often miss visually evident CTCF loops in Hi-C data sets from mammals, and they completely fail to identify high intensity loops in other organisms. We present SIP, Significant Interaction Peak caller, and SIPMeta, which are platform independent programs to identify and characterize these loops in a time- and memory-efficient manner. We show that SIP is resistant to noise and sequencing depth, and can be used to detect loops that were previously missed in human cells as well as loops in other organisms. SIPMeta corrects for a common visualization artifact by accounting for Manhattan distance to create average plots of Hi-C and HiChIP data. We then demonstrate that the use of SIP and SIPMeta can lead to biological insights by characterizing the contribution of several transcription factors to CTCF loop stability in human cells. We also annotate loops associated with the SMC component of the dosage compensation complex (DCC) in Caenorhabditis elegans and demonstrate that loop anchors represent bidirectional blocks for symmetrical loop extrusion. This is in contrast to the asymmetrical extrusion until unidirectional blockage by CTCF that is presumed to occur in mammals. Using HiChIP and multiway ligation events, we then show that DCC loops form a network of strong interactions that may contribute to X Chromosome-wide condensation in C. elegans hermaphrodites.

Widespread long-range cis-regulatory elements in the maize genome

Genetic mapping studies on crops suggest that agronomic traits can be controlled by gene-distal intergenic loci. Despite the biological importance and the potential agronomic utility of these loci, they remain virtually uncharacterized in all crop species to date. Here, we provide genetic, epigenomic and functional molecular evidence to support the widespread existence of gene-distal (hereafter, distal) loci that act as long-range transcriptional cis-regulatory elements (CREs) in the maize genome. Such loci are enriched for euchromatic features that suggest their regulatory functions. Chromatin loops link together putative CREs with genes and recapitulate genetic interactions. Putative CREs also display elevated transcriptional enhancer activities, as measured by self-transcribing active regulatory region sequencing. These results provide functional support for the widespread existence of CREs that act over large genomic distances to control gene expression.

Chromatin accessibility and transcription dynamics during in vitro astrocyte differentiation of Huntington's Disease Monkey pluripotent stem cells

BACKGROUND

Huntington's Disease (HD) is a fatal neurodegenerative disorder caused by a CAG repeat expansion, resulting in a mutant huntingtin protein. While it is now clear that astrocytes are affected by HD and significantly contribute to neuronal dysfunction and pathogenesis, the alterations in the transcriptional and epigenetic profiles in HD astrocytes have yet to be characterized. Here, we examine global transcription and chromatin accessibility dynamics during in vitro astrocyte differentiation in a transgenic non-human primate model of HD.

RESULTS

We found global changes in accessibility and transcription across different stages of HD pluripotent stem cell differentiation, with distinct trends first observed in neural progenitor cells (NPCs), once cells have committed to a neural lineage. Transcription of p53 signaling and cell cycle pathway genes was highly impacted during differentiation, with depletion in HD NPCs and upregulation in HD astrocytes. E2F target genes also displayed this inverse expression pattern, and strong associations between E2F target gene expression and accessibility at nearby putative enhancers were observed.

CONCLUSIONS

The results suggest that chromatin accessibility and transcription are altered throughout in vitro HD astrocyte differentiation and provide evidence that E2F dysregulation contributes to aberrant cell-cycle re-entry and apoptosis throughout the progression from NPCs to astrocytes.

Immediate and deferred epigenomic signatures of in vivo neuronal activation in mouse hippocampus

Activity-driven transcription plays an important role in many brain processes, including those underlying memory and epilepsy. Here we combine genetic tagging of nuclei and ribosomes with RNA sequencing, chromatin immunoprecipitation with sequencing, assay for transposase-accessible chromatin using sequencing and Hi-C to investigate transcriptional and chromatin changes occurring in mouse hippocampal excitatory neurons at different time points after synchronous activation during seizure and sparse activation by novel context exploration. The transcriptional burst is associated with an increase in chromatin accessibility of activity-regulated genes and enhancers, de novo binding of activity-regulated transcription factors, augmented promoter-enhancer interactions and the formation of gene loops that bring together the transcription start site and transcription termination site of induced genes and may sustain the fast reloading of RNA polymerase complexes. Some chromatin occupancy changes and interactions, particularly those driven by AP1, remain long after neuronal activation and could underlie the changes in neuronal responsiveness and circuit connectivity observed in these neuroplasticity paradigms, perhaps thereby contributing to metaplasticity in the adult brain.

Maintenance of CTCF- and Transcription Factor-Mediated Interactions from the Gametes to the Early Mouse Embryo

The epigenetic information present in mammalian gametes and whether it is transmitted to the progeny are relatively unknown. We find that many promoters in mouse sperm are occupied by RNA polymerase II (Pol II) and Mediator. The same promoters are accessible in GV and MII oocytes and preimplantation embryos. Sperm distal ATAC-seq sites containing motifs for various transcription factors are conserved in monkeys and humans. ChIP-seq analyses confirm that Foxa1, ERα, and AR occupy distal enhancers in sperm. Accessible sperm enhancers containing H3.3 and H2A.Z are also accessible in oocytes and preimplantation embryos. Furthermore, their interactions with promoters in the gametes persist during early development. Sperm- or oocyte-specific interactions mediated by CTCF and cohesin are only present in the paternal or maternal chromosomes, respectively, in the zygote and 2-cell stages. These interactions converge in both chromosomes by the 8-cell stage. Thus, mammalian gametes contain complex patterns of 3D interactions that can be transmitted to the zygote after fertilization.

Rapid Irreversible Transcriptional Reprogramming in Human Stem Cells Accompanied by Discordance between Replication Timing and Chromatin Compartment

The temporal order of DNA replication is regulated during development and is highly correlated with gene expression, histone modifications and 3D genome architecture. We tracked changes in replication timing, gene expression, and chromatin conformation capture (Hi-C) A/B compartments over the first two cell cycles during differentiation of human embryonic stem cells to definitive endoderm. Remarkably, transcriptional programs were irreversibly reprogrammed within the first cell cycle and were largely but not universally coordinated with replication timing changes. Moreover, changes in A/B compartment and several histone modifications that normally correlate strongly with replication timing showed weak correlation during the early cell cycles of differentiation but showed increased alignment in later differentiation stages and in terminally differentiated cell lines. Thus, epigenetic cell fate transitions during early differentiation can occur despite dynamic and discordant changes in otherwise highly correlated genomic properties.

Pach-ing It in: The Peculiar Organization of Mammalian Pachytene Chromosomes

In this issue of Molecular Cell, Wang et al. (2019) use Hi-C to visualize at high resolution the complex reprogramming of chromatin architecture during spermatogenesis in rhesus monkeys and mice. They find that pachytene spermatocytes have a unique chromosome organization that may result from the presence of the synaptonemal complex and transcription-associated proteins.

Resolution of the DNA methylation state of single CpG dyads using in silico strand annealing and WGBS data

Whole-genome bisulfite sequencing (WGBS) has been widely used to quantify cytosine DNA methylation frequency in an expanding array of cell and tissue types. Because of the denaturing conditions used, this method ultimately leads to the measurement of methylation frequencies at single cytosines. Hence, the methylation frequency of CpG dyads (two complementary CG dinucleotides) can be only indirectly inferred by overlaying the methylation frequency of two cytosines measured independently. Furthermore, hemi-methylated CpGs (hemiCpGs) have not been previously analyzed in WGBS studies. We recently developed in silico strand annealing (iSA), a bioinformatics method applicable to WGBS data, to resolve the methylation status of CpG dyads into unmethylated, hemi-methylated, and methylated. HemiCpGs account for 4-20% of the DNA methylome in different cell types, and some can be inherited across cell divisions, suggesting a role as a stable epigenetic mark. Therefore, it is important to resolve hemiCpGs from fully methylated CpGs in WGBS studies. This protocol describes step-by-step commands to accomplish this task, including dividing alignments by strand, pairing alignments between strands, and extracting single-fragment methylation calls. The versatility of iSA enables its application downstream of other WGBS-related methods such as nasBS-seq (nascent DNA bisulfite sequencing), ChIP-BS-seq (ChIP followed by bisulfite sequencing), TAB-seq, oxBS-seq, and fCAB-seq. iSA is also tunable for analyzing the methylation status of cytosines in any sequence context. We exemplify this flexibility by uncovering the single-fragment non-CpG methylome. This protocol provides enough details for users with little experience in bioinformatic analysis and takes 2-7 h.

Condensin II Counteracts Cohesin and RNA Polymerase II in the Establishment of 3D Chromatin Organization

Interaction domains in Drosophila chromosomes form by segregation of active and inactive chromatin in the absence of CTCF loops, but the role of transcription versus other architectural proteins in chromatin organization is unclear. Here, we find that positioning of RNAPII via transcription elongation is essential in the formation of gene loops, which in turn interact to form compartmental domains. Inhibition of transcription elongation or depletion of cohesin decreases gene looping and formation of active compartmental domains. In contrast, depletion of condensin II, which also localizes to active chromatin, causes increased gene looping, formation of compartmental domains, and stronger intra-chromosomal compartmental interactions. Condensin II has a similar role in maintaining inter-chromosomal interactions responsible for pairing between homologous chromosomes, whereas inhibition of transcription elongation or cohesin depletion has little effect on homolog pairing. The results suggest distinct roles for cohesin and condensin II in the establishment of 3D nuclear organization in Drosophila.

Chromatin changes in Anopheles gambiae induced by Plasmodium falciparum infection

BACKGROUND

Infection by the human malaria parasite leads to important changes in mosquito phenotypic traits related to vector competence. However, we still lack a clear understanding of the underlying mechanisms and, in particular, of the epigenetic basis for these changes. We have examined genome-wide distribution maps of H3K27ac, H3K9ac, H3K9me3 and H3K4me3 by ChIP-seq and the transcriptome by RNA-seq, of midguts from Anopheles gambiae mosquitoes blood-fed uninfected and infected with natural isolates of the human malaria parasite Plasmodium falciparum in Burkina Faso.

RESULTS

We report 15,916 regions containing differential histone modification enrichment between infected and uninfected, of which 8339 locate at promoters and/or intersect with genes. The functional annotation of these regions allowed us to identify infection-responsive genes showing differential enrichment in various histone modifications, such as CLIP proteases, antimicrobial peptides-encoding genes, and genes related to melanization responses and the complement system. Further, the motif analysis of regions differentially enriched in various histone modifications predicts binding sites that might be involved in the cis-regulation of these regions, such as Deaf1, Pangolin and Dorsal transcription factors (TFs). Some of these TFs are known to regulate immunity gene expression in Drosophila and are involved in the Notch and JAK/STAT signaling pathways.

CONCLUSIONS

The analysis of malaria infection-induced chromatin changes in mosquitoes is important not only to identify regulatory elements and genes underlying mosquito responses to P. falciparum infection, but also for possible applications to the genetic manipulation of mosquitoes and to other mosquito-borne systems.

Organizational principles of 3D genome architecture

Studies of 3D chromatin organization have suggested that chromosomes are hierarchically organized into large compartments composed of smaller domains called topologically associating domains (TADs). Recent evidence suggests that compartments are smaller than previously thought and that the transcriptional or chromatin state is responsible for interactions leading to the formation of small compartmental domains in all organisms. In vertebrates, CTCF forms loop domains, probably via an extrusion process involving cohesin. CTCF loops cooperate with compartmental domains to establish the 3D organization of the genome. The continuous extrusion of the chromatin fibre by cohesin may also be responsible for the establishment of enhancer-promoter interactions and stochastic aspects of the transcription process. These observations suggest that the 3D organization of the genome is an emergent property of chromatin and its components, and thus may not be only a determinant but also a consequence of its function.

A tethered-inchworm model of SMC DNA translocation

The DNA loop extrusion model is a provocative new concept explaining the formation of chromatin loops that revolutionizes understanding of genome organization. Central to this model is the structural maintenance of chromosomes (SMC) protein family, which is now thought to function as a DNA motor. In this Perspective, we review and reinterpret the current knowledge of SMC structure and function and propose a novel mechanism for SMC motor activity.

Architectural Proteins and Pluripotency Factors Cooperate to Orchestrate the Transcriptional Response of hESCs to Temperature Stress

Cells respond to temperature stress via up- and downregulation of hundreds of genes. This process is thought to be regulated by the heat shock factor HSF1, which controls the release of RNAPII from promoter-proximal pausing. Here, we analyze the events taking place in hESCs upstream of RNAPII release. We find that temperature stress results in the activation or decommissioning of thousands of enhancers. This process involves alterations in the occupancy of transcription factors HSF1, AP-1, NANOG, KLF4, and OCT4 accompanied by nucleosome remodeling by BRG1 and changes in H3K27ac. Furthermore, redistribution of RAD21 and CTCF results in the formation and disassembly of interactions mediated by these two proteins. These alterations tether and untether enhancers to their cognate promoters or refashion insulated neighborhoods, thus transforming the landscape of enhancer-promoter interactions. Details of the 3D interactome remodeling process support loop extrusion initiating at random sites as a mechanism for the establishment of CTCF/cohesin loops.

Developing in 3D: the role of CTCF in cell differentiation

CTCF is a highly conserved zinc-finger DNA-binding protein that mediates interactions between distant sequences in the genome. As a consequence, CTCF regulates enhancer-promoter interactions and contributes to the three-dimensional organization of the genome. Recent studies indicate that CTCF is developmentally regulated, suggesting that it plays a role in cell type-specific genome organization. Here, we review these studies and discuss how CTCF functions during the development of various cell and tissue types, ranging from embryonic stem cells and gametes, to neural, muscle and cardiac cells. We propose that the lineage-specific control of CTCF levels, and its partnership with lineage-specific transcription factors, allows for the control of cell type-specific gene expression via chromatin looping.

Nascent DNA methylome mapping reveals inheritance of hemimethylation at CTCF/cohesin sites

The faithful inheritance of the epigenome is critical for cells to maintain gene expression programs and cellular identity across cell divisions. We mapped strand-specific DNA methylation after replication forks and show maintenance of the vast majority of the DNA methylome within 20 minutes of replication and inheritance of some hemimethylated CpG dinucleotides (hemiCpGs). Mapping the nascent DNA methylome targeted by each of the three DNA methyltransferases (DNMTs) reveals interactions between DNMTs and substrate daughter cytosines en route to maintenance methylation or hemimethylation. Finally, we show the inheritance of hemiCpGs at short regions flanking CCCTC-binding factor (CTCF)/cohesin binding sites in pluripotent cells. Elimination of hemimethylation causes reduced frequency of chromatin interactions emanating from these sites, suggesting a role for hemimethylation as a stable epigenetic mark regulating CTCF-mediated chromatin interactions.

Not just heads and tails: The complexity of the sperm epigenome

Transgenerational inheritance requires mechanisms by which epigenetic information is transferred via gametes. Canonical thought holds that mammalian sperm chromatin would be incapable of carrying epigenetic information as post-translational modifications of histones because of their replacement with protamine proteins. Furthermore, compaction of the sperm genome would hinder DNA accessibility of proteins involved in transcriptional regulation and genome architecture. In this Minireview, we delineate the paternal chromatin remodeling events during spermatogenesis and fertilization. Sperm chromatin is epigenetically modified at various time points throughout its development. This allows for the addition of environment-specific modifications that can be passed from parents to offspring.

Analysis of Chromatin Interactions Mediated by Specific Architectural Proteins in Drosophila Cells

Chromosome conformation capture assays have been established, modified, and enhanced for over a decade with the purpose of studying nuclear organization. A recently published method uses in situ Hi-C followed by chromatin immunoprecipitation (HiChIP) to enrich the overall yield of significant genome-wide interactions mediated by a specific protein. Here we applied a modified version of the HiChIP protocol to retrieve the significant contacts mediated by architectural protein CP190 in D. melanogaster cells.

Genome-Wide Mapping of Protein-DNA Interactions on Nascent Chromatin

Chromatin immunoprecipitation (ChIP) is the most widely used method to analyze protein-DNA interactions in vivo. Coupled with next generation sequencing, ChIP-seq experiments map protein-DNA interactions in a genome-wide fashion. Here we describe a novel method called nasChIP-seq for mapping genome-wide occupancy of posttranslationally modified histones or transcription factors on newly replicated DNA.

Different enhancer classes in Drosophila bind distinct architectural proteins and mediate unique chromatin interactions and 3D architecture

Eukaryotic gene expression is regulated by enhancer–promoter interactions but the molecular mechanisms that govern specificity have remained elusive. Genome-wide studies utilizing STARR-seq identified two enhancer classes in Drosophila that interact with different core promoters: housekeeping enhancers (hkCP) and developmental enhancers (dCP). We hypothesized that the two enhancer classes are occupied by distinct architectural proteins, affecting their enhancer–promoter contacts. By evaluating ChIP-seq occupancy of architectural proteins, typical enhancer-associated proteins, and histone modifications, we determine that both enhancer classes are enriched for RNA Polymerase II, CBP, and architectural proteins but there are also distinctions. hkCP enhancers contain H3K4me3 and exclusively bind Cap-H2, Chromator, DREF and Z4, whereas dCP enhancers contain H3K4me1 and are more enriched for Rad21 and Fs(1)h-L. Additionally, we map the interactions of each enhancer class utilizing a Hi-C dataset with <1 kb resolution. Results suggest that hkCP enhancers are more likely to form multi-TSS interaction networks and be associated with topologically associating domain (TAD) borders, while dCP enhancers are more often bound to one or two TSSs and are enriched at chromatin loop anchors. The data support a model suggesting that the unique architectural protein occupancy within enhancers is one contributor to enhancer–promoter interaction specificity.

Minute-Made Data Analysis: Tools for Rapid Interrogation of Hi-C Contacts

Juicer and Juicebox, described by Durand et al. (2016a, 2016b), are two new tools for fast and reliable processing of Hi-C data, providing approaches for read processing, multiple normalization schemes, feature annotation, and dynamic browsing of chromatin contacts, thus reducing arduous Hi-C analysis into an easy yet flexible pipeline.

Evolutionarily Conserved Principles Predict 3D Chromatin Organization

Topologically associating domains (TADs), CTCF loop domains, and A/B compartments have been identified as important structural and functional components of 3D chromatin organization, yet the relationship between these features is not well understood. Using high-resolution Hi-C and HiChIP, we show that Drosophila chromatin is organized into domains we term compartmental domains that correspond precisely with A/B compartments at high resolution. We find that transcriptional state is a major predictor of Hi-C contact maps in several eukaryotes tested, including C. elegans and A. thaliana. Architectural proteins insulate compartmental domains by reducing interaction frequencies between neighboring regions in Drosophila, but CTCF loops do not play a distinct role in this organism. In mammals, compartmental domains exist alongside CTCF loop domains to form topological domains. The results suggest that compartmental domains are responsible for domain structure in all eukaryotes, with CTCF playing an important role in domain formation in mammals.

Local Genome Topology Can Exhibit an Incompletely Rewired 3D-Folding State during Somatic Cell Reprogramming

Pluripotent genomes are folded in a topological hierarchy that reorganizes during differentiation. The extent to which chromatin architecture is reconfigured during somatic cell reprogramming is poorly understood. Here we integrate fine-resolution architecture maps with epigenetic marks and gene expression in embryonic stem cells (ESCs), neural progenitor cells (NPCs), and NPC-derived induced pluripotent stem cells (iPSCs). We find that most pluripotency genes reconnect to target enhancers during reprogramming. Unexpectedly, some NPC interactions around pluripotency genes persist in our iPSC clone. Pluripotency genes engaged in both "fully-reprogrammed" and "persistent-NPC" interactions exhibit over/undershooting of target expression levels in iPSCs. Additionally, we identify a subset of "poorly reprogrammed" interactions that do not reconnect in iPSCs and display only partially recovered, ESC-specific CTCF occupancy. 2i/LIF can abrogate persistent-NPC interactions, recover poorly reprogrammed interactions, reinstate CTCF occupancy, and restore expression levels. Our results demonstrate that iPSC genomes can exhibit imperfectly rewired 3D-folding linked to inaccurately reprogrammed gene expression.

Structural Basis for the Versatile and Methylation-Dependent Binding of CTCF to DNA

The multidomain CCCTC-binding factor (CTCF), containing a tandem array of 11 zinc fingers (ZFs), modulates the three-dimensional organization of chromatin. We crystallized the human CTCF DNA-binding domain in complex with a known CTCF-binding site. While ZF2 does not make sequence-specific contacts, each finger of ZF3-7 contacts three bases of the 15-bp consensus sequence. Each conserved nucleotide makes base-specific hydrogen bonds with a particular residue. Most of the variable base pairs within the core sequence also engage in interactions with the protein. These interactions compensate for deviations from the consensus sequence, allowing CTCF to adapt to sequence variations. CTCF is sensitive to cytosine methylation at position 2, but insensitive at position 12 of the 15-bp core sequence. These differences can be rationalized structurally. Although included in crystallizations, ZF10 and ZF11 are not visible, while ZF8 and ZF9 span the backbone of the DNA duplex, conferring no sequence specificity but adding to overall binding stability.

Chromatin States in Mouse Sperm Correlate with Embryonic and Adult Regulatory Landscapes

The mammalian sperm genome is thought to lack substantial information for the regulation of future expression after fertilization. Here, we show that most promoters in mouse sperm are flanked by well-positioned nucleosomes marked by active histone modifications. Analysis of these modifications suggests that many enhancers and super-enhancers functional in embryonic and adult tissues are already specified in sperm. The sperm genome is bound by CTCF and cohesin at sites that are also present in round spermatids and embryonic stem cells (ESCs). These sites mediate interactions that organize the sperm genome into domains and compartments that overlap extensively with those found in mESCs. These results suggest that sperm carry a rich source of regulatory information, encoded in part by its three-dimensional folding specified by CTCF and cohesin. This information may contribute to future expression during embryonic and adult life, suggesting mechanisms by which environmental effects on the paternal germline are transmitted transgenerationally.

Chromatin States in Mouse Sperm Correlate with Embryonic and Adult Regulatory Landscapes

The mammalian sperm genome is thought to lack substantial information for the regulation of future expression after fertilization. Here, we show that most promoters in mouse sperm are flanked by well-positioned nucleosomes marked by active histone modifications. Analysis of these modifications suggests that many enhancers and super-enhancers functional in embryonic and adult tissues are already specified in sperm. The sperm genome is bound by CTCF and cohesin at sites that are also present in round spermatids and embryonic stem cells (ESCs). These sites mediate interactions that organize the sperm genome into domains and compartments that overlap extensively with those found in mESCs. These results suggest that sperm carry a rich source of regulatory information, encoded in part by its three-dimensional folding specified by CTCF and cohesin. This information may contribute to future expression during embryonic and adult life, suggesting mechanisms by which environmental effects on the paternal germline are transmitted transgenerationally.

Epigenetic regulation of Plasmodium falciparum clonally variant gene expression during development in Anopheles gambiae

P. falciparum phenotypic plasticity is linked to the variant expression of clonal multigene families such as the var genes. We have examined changes in transcription and histone modifications that occur during sporogonic development of P. falciparum in the mosquito host. All var genes are silenced or transcribed at low levels in blood stages (gametocyte/ring) of the parasite in the human host. After infection of mosquitoes, a single var gene is selected for expression in the oocyst, and transcription of this gene increases dramatically in the sporozoite. The same PF3D7_1255200 var gene was activated in 4 different experimental infections. Transcription of this var gene during parasite development in the mosquito correlates with the presence of low levels of H3K9me3 at the binding site for the PF3D7_1466400 AP2 transcription factor. This chromatin state in the sporozoite also correlates with the expression of an antisense long non-coding RNA (lncRNA) that has previously been shown to promote var gene transcription during the intraerythrocytic cycle in vitro. Expression of both the sense protein-coding transcript and the antisense lncRNA increase dramatically in sporozoites. The findings suggest a complex process for the activation of a single particular var gene that involves AP2 transcription factors and lncRNAs.

Capturing native interactions: intrinsic methods to study chromatin conformation

The 3D organization of chromatin controls gene expression through spatial interactions between genomic loci. FISH and 3C‐based methods that are commonly used to study chromatin organization utilize chemical crosslinking, a step that may introduce biases in detectable chromatin interactions. In their recent study, Papantonis and colleagues (Brant et al, 2016) developed alternative new methods of detecting chromatin contacts without the use of chemical crosslinking agents. These tools increase the resolution and confidence at which interactions can be identified, and may be informative for chromatin interaction dynamics.