Widespread plasticity in CTCF occupancy linked to DNA methylation

Permeable TAD boundaries and their impact on genome-associated functions

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

Progress in multifactorial single-cell chromatin profiling methods

Chromatin states play a key role in shaping overall cellular states and fates. Building a complete picture of the functional state of chromatin in cells requires the co-detection of several distinct biochemical aspects. These span DNA methylation, chromatin accessibility, chromosomal conformation, histone posttranslational modifications, and more. While this certainly presents a challenging task, over the past few years many new and creative methods have been developed that now enable co-assay of these different aspects of chromatin at single cell resolution. This field is entering an exciting phase, where a confluence of technological improvements, decreased sequencing costs, and computational innovation are presenting new opportunities to dissect the diversity of chromatin states present in tissues, and how these states may influence gene regulation. In this review, I discuss the spectrum of current experimental approaches for multifactorial chromatin profiling, highlight some of the experimental and analytical challenges, as well as some areas for further innovation.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

SETDB1 regulates short interspersed nuclear elements and chromatin loop organization in mouse neural precursor cells

BACKGROUND

Transposable elements play a critical role in maintaining genome architecture during neurodevelopment. Short Interspersed Nuclear Elements (SINEs), a major subtype of transposable elements, are known to harbor binding sites for the CCCTC-binding factor (CTCF) and pivotal in orchestrating chromatin organization. However, the regulatory mechanisms controlling the activity of SINEs in the developing brain remains elusive.

RESULTS

In our study, we conduct a comprehensive genome-wide epigenetic analysis in mouse neural precursor cells using ATAC-seq, ChIP-seq, whole genome bisulfite sequencing, in situ Hi-C, and RNA-seq. Our findings reveal that the SET domain bifurcated histone lysine methyltransferase 1 (SETDB1)-mediated H3K9me3, in conjunction with DNA methylation, restricts chromatin accessibility on a selective subset of SINEs in neural precursor cells. Mechanistically, loss of Setdb1 increases CTCF access to these SINE elements and contributes to chromatin loop reorganization. Moreover, de novo loop formation contributes to differential gene expression, including the dysregulation of genes enriched in mitotic pathways. This leads to the disruptions of cell proliferation in the embryonic brain after genetic ablation of Setdb1 both in vitro and in vivo.

CONCLUSIONS

In summary, our study sheds light on the epigenetic regulation of SINEs in mouse neural precursor cells, suggesting their role in maintaining chromatin organization and cell proliferation during neurodevelopment.

cBAF generates subnucleosomes that expand OCT4 binding and function beyond DNA motifs at enhancers

The canonical BRG/BRM-associated factor (cBAF) complex is essential for chromatin opening at enhancers in mammalian cells. However, the nature of the open chromatin remains unclear. Here, we show that, in addition to producing histone-free DNA, cBAF generates stable hemisome-like subnucleosomal particles containing the four core histones associated with 50-80 bp of DNA. Our genome-wide analysis indicates that cBAF makes these particles by targeting and splitting fragile nucleosomes. In mouse embryonic stem cells, these subnucleosomes become an in vivo binding substrate for the master transcription factor OCT4 independently of the presence of OCT4 DNA motifs. At enhancers, the OCT4-subnucleosome interaction increases OCT4 occupancy and amplifies the genomic interval bound by OCT4 by up to one order of magnitude compared to the region occupied on histone-free DNA. We propose that cBAF-dependent subnucleosomes orchestrate a molecular mechanism that projects OCT4 function in chromatin opening beyond its DNA motifs.

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

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

The long non-coding RNA Meg3 mediates imprinted gene expression during stem cell differentiation

The imprinted Dlk1-Dio3 domain comprises the developmental genes Dlk1 and Rtl1, which are silenced on the maternal chromosome in different cell types. On this parental chromosome, the domain's imprinting control region activates a polycistron that produces the lncRNA Meg3 and many miRNAs (Mirg) and C/D-box snoRNAs (Rian). Although Meg3 lncRNA is nuclear and associates with the maternal chromosome, it is unknown whether it controls gene repression in cis. We created mouse embryonic stem cells (mESCs) that carry an ectopic poly(A) signal, reducing RNA levels along the polycistron, and generated Rian-/- mESCs as well. Upon ESC differentiation, we found that Meg3 lncRNA (but not Rian) is required for Dlk1 repression on the maternal chromosome. Biallelic Meg3 expression acquired through CRISPR-mediated demethylation of the paternal Meg3 promoter led to biallelic Dlk1 repression, and to loss of Rtl1 expression. lncRNA expression also correlated with DNA hypomethylation and CTCF binding at the 5'-side of Meg3. Using Capture Hi-C, we found that this creates a Topologically Associating Domain (TAD) organization that brings Meg3 close to Dlk1 on the maternal chromosome. The requirement of Meg3 for gene repression and TAD structure may explain how aberrant MEG3 expression at the human DLK1-DIO3 locus associates with imprinting disorders.

Regeneration-specific promoter switching facilitates Mest expression in the mouse digit tip to modulate neutrophil response

The mouse digit tip regenerates following amputation, a process mediated by a cellularly heterogeneous blastema. We previously found the gene Mest to be highly expressed in mesenchymal cells of the blastema and a strong candidate pro-regenerative gene. We now show Mest digit expression is regeneration-specific and not upregulated in post-amputation fibrosing proximal digits. Mest homozygous knockout mice exhibit delayed bone regeneration though no phenotype is found in paternal knockout mice, inconsistent with the defined maternal genomic imprinting of Mest. We demonstrate that promoter switching, not loss of imprinting, regulates biallelic Mest expression in the blastema and does not occur during embryogenesis, indicating a regeneration-specific mechanism. Requirement for Mest expression is tied to modulating neutrophil response, as revealed by scRNAseq and FACS comparing wildtype and knockout blastemas. Collectively, the imprinted gene Mest is required for proper digit tip regeneration and its blastema expression is facilitated by promoter switching for biallelic expression.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Job-related exhaustion risk variant in UST is associated with dementia and DNA methylation

Previous genome-wide association and replication study for job-related exhaustion indicated a risk variant, rs13219957 in the UST gene. Epidemiological studies suggest connection of stress-related conditions and dementia risk. Therefore, we first studied association of rs13219957 and register-based incident dementia using survival models in the Finnish National FINRISK study surveys (N = 26,693). The AA genotype of rs13219957 was significantly associated with 40% increased risk of all-cause dementia. Then we analysed the UST locus association with brain pathology in the Vantaa 85+ cohort and found association with tau pathology (Braak stage) but not with amyloid pathology. Finally, in the functional analyses, rs13219957 showed a highly significant association with two DNA methylation sites of UST, and UST expression. Thus, the results suggest a common risk variant for a stress-related condition and dementia. Mechanisms to mediate the connection may include differential DNA methylation and transcriptional regulation of UST.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Molecular basis for differential Igk versus Igh V(D)J joining mechanisms

In developing B cells, V(D)J recombination assembles exons encoding IgH and Igκ variable regions from hundreds of gene segments clustered across Igh and Igk loci. V, D and J gene segments are flanked by conserved recombination signal sequences (RSSs) that target RAG endonuclease1. RAG orchestrates Igh V(D)J recombination upon capturing a JH-RSS within the JH-RSS-based recombination centre1-3 (RC). JH-RSS orientation programmes RAG to scan upstream D- and VH-containing chromatin that is presented in a linear manner by cohesin-mediated loop extrusion4-7. During Igh scanning, RAG robustly utilizes only D-RSSs or VH-RSSs in convergent (deletional) orientation with JH-RSSs4-7. However, for Vκ-to-Jκ joining, RAG utilizes Vκ-RSSs from deletional- and inversional-oriented clusters8, inconsistent with linear scanning2. Here we characterize the Vκ-to-Jκ joining mechanism. Igk undergoes robust primary and secondary rearrangements9,10, which confounds scanning assays. We therefore engineered cells to undergo only primary Vκ-to-Jκ rearrangements and found that RAG scanning from the primary Jκ-RC terminates just 8 kb upstream within the CTCF-site-based Sis element11. Whereas Sis and the Jκ-RC barely interacted with the Vκ locus, the CTCF-site-based Cer element12 4 kb upstream of Sis interacted with various loop extrusion impediments across the locus. Similar to VH locus inversion7, DJH inversion abrogated VH-to-DJH joining; yet Vκ locus or Jκ inversion allowed robust Vκ-to-Jκ joining. Together, these experiments implicated loop extrusion in bringing Vκ segments near Cer for short-range diffusion-mediated capture by RC-based RAG. To identify key mechanistic elements for diffusional V(D)J recombination in Igk versus Igh, we assayed Vκ-to-JH and D-to-Jκ rearrangements in hybrid Igh-Igk loci generated by targeted chromosomal translocations, and pinpointed remarkably strong Vκ and Jκ RSSs. Indeed, RSS replacements in hybrid or normal Igk and Igh loci confirmed the ability of Igk-RSSs to promote robust diffusional joining compared with Igh-RSSs. We propose that Igk evolved strong RSSs to mediate diffusional Vκ-to-Jκ joining, whereas Igh evolved weaker RSSs requisite for modulating VH joining by RAG-scanning impediments.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

The PTM profiling of CTCF reveals the regulation of 3D chromatin structure by O-GlcNAcylation

CCCTC-binding factor (CTCF), a ubiquitously expressed and highly conserved protein, is known to play a critical role in chromatin structure. Post-translational modifications (PTMs) diversify the functions of protein to regulate numerous cellular processes. However, the effects of PTMs on the genome-wide binding of CTCF and the organization of three-dimensional (3D) chromatin structure have not been fully understood. In this study, we uncovered the PTM profiling of CTCF and demonstrated that CTCF can be O-GlcNAcylated and arginine methylated. Functionally, we demonstrated that O-GlcNAcylation inhibits CTCF binding to chromatin. Meanwhile, deficiency of CTCF O-GlcNAcylation results in the disruption of loop domains and the alteration of chromatin loops associated with cellular development. Furthermore, the deficiency of CTCF O-GlcNAcylation increases the expression of developmental genes and negatively regulates maintenance and establishment of stem cell pluripotency. In conclusion, these results provide key insights into the role of PTMs for the 3D chromatin structure.

DNA methylation-dependent and -independent binding of CDX2 directs activation of distinct developmental and homeostatic genes

Precise spatiotemporal and cell type-specific gene expression is essential for proper tissue development and function. Transcription factors (TFs) guide this process by binding to developmental stage-specific targets and establishing an appropriate enhancer landscape. In turn, DNA and chromatin modifications direct the genomic binding of TFs. However, how TFs navigate various chromatin features and selectively bind a small portion of the millions of possible genomic target loci is still not well understood. Here we show that Cdx2 - a pioneer TF that binds distinct targets in developing versus adult intestinal epithelial cells - has a preferential affinity for a non-canonical CpG-containing motif in vivo. A higher frequency of this motif at embryonic and fetal Cdx2 target loci and the specifically methylated state of the CpG during development allows selective Cdx2 binding and activation of developmental enhancers and linked genes. Conversely, demethylation at these enhancers prohibits ectopic Cdx2 binding in adult cells, where Cdx2 binds its canonical motif without a CpG. This differential Cdx2 binding allows for corecruitment of Ctcf and Hnf4, facilitating the establishment of intestinal superenhancers during development and enhancers mediating adult homeostatic functions, respectively. Induced gain of DNA methylation in the adult mouse epithelium or cultured cells causes ectopic recruitment of Cdx2 to the developmental target loci and facilitates cobinding of the partner TFs. Together, our results demonstrate that the differential CpG motif requirements for Cdx2 binding to developmental versus adult target sites allow it to navigate different DNA methylation profiles and activate cell type-specific genes at appropriate times.

Identification of DNA methylation episignature for the intellectual developmental disorder, autosomal dominant 21 syndrome, caused by variants in the CTCF gene

PURPOSE

The main objective of this study was to assess clinical features and genome-wide DNA methylation profiles in individuals affected by intellectual developmental disorder, autosomal dominant 21 (IDD21) syndrome, caused by variants in the CCCTC-binding factor (CTCF) gene.

METHODS

DNA samples were extracted from peripheral blood of 16 individuals with clinical features and genetic findings consistent with IDD21. DNA methylation analysis was performed using the Illumina Infinium Methylation EPIC Bead Chip microarrays. The methylation levels were fitted in a multivariate linear regression model to identify the differentially methylated probes. A binary support vector machine classification model was constructed to differentiate IDD21 samples from controls.

RESULTS

We identified a highly specific, reproducible, and sensitive episignature associated with CTCF variants. Six variants of uncertain significance were tested, of which 2 mapped to the IDD21 episignature and clustered alongside IDD21 cases in both heatmap and multidimensional scaling plots. Comparison of the genomic DNA methylation profile of IDD21 with that of 56 other neurodevelopmental disorders provided insights into the underlying molecular pathophysiology of this disorder.

CONCLUSION

The robust and specific CTCF/IDD21 episignature expands the growing list of neurodevelopmental disorders with distinct DNA methylation profiles, which can be applied as supporting evidence in variant classification.

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

Noncanonical functions of microRNAs in the nucleus

MicroRNAs (miRNAs) are small noncoding RNAs (ncRNAs) that play their roles in the regulation of physiological and pathological processes. Originally, it was assumed that miRNAs only modulate gene expression posttranscriptionally in the cytoplasm by inducing target mRNA degradation. However, with further research, evidence shows that mature miRNAs also exist in the cell nucleus, where they can impact gene transcription and ncRNA maturation in several ways. This review provides an overview of novel models of nuclear miRNA functions. Some of the models remain to be verified by experimental evidence, and more details of the miRNA regulation network remain to be discovered in the future.

Epiallelic variation of non-coding RNA genes and their phenotypic consequences

Epigenetic variations contribute greatly to the phenotypic plasticity and diversity. Current functional studies on epialleles have predominantly focused on protein-coding genes, leaving the epialleles of non-coding RNA (ncRNA) genes largely understudied. Here, we uncover abundant DNA methylation variations of ncRNA genes and their significant correlations with plant adaptation among 1001 natural Arabidopsis accessions. Through genome-wide association study (GWAS), we identify large numbers of methylation QTL (methylQTL) that are independent of known DNA methyltransferases and enriched in specific chromatin states. Proximal methylQTL closely located to ncRNA genes have a larger effect on DNA methylation than distal methylQTL. We ectopically tether a DNA methyltransferase MQ1v to miR157a by CRISPR-dCas9 and show de novo establishment of DNA methylation accompanied with decreased miR157a abundance and early flowering. These findings provide important insights into the genetic basis of epigenetic variations and highlight the contribution of epigenetic variations of ncRNA genes to plant phenotypes and diversity.

Cell-type differential targeting of SETDB1 prevents aberrant CTCF binding, chromatin looping, and cis-regulatory interactions

SETDB1 is an essential histone methyltransferase that deposits histone H3 lysine 9 trimethylation (H3K9me3) to transcriptionally repress genes and repetitive elements. The function of differential H3K9me3 enrichment between cell-types remains unclear. Here, we demonstrate mutual exclusivity of H3K9me3 and CTCF across mouse tissues from different developmental timepoints. We analyze SETDB1 depleted cells and discover that H3K9me3 prevents aberrant CTCF binding independently of DNA methylation and H3K9me2. Such sites are enriched with SINE B2 retrotransposons. Moreover, analysis of higher-order genome architecture reveals that large chromatin structures including topologically associated domains and subnuclear compartments, remain intact in SETDB1 depleted cells. However, chromatin loops and local 3D interactions are disrupted, leading to transcriptional changes by modifying pre-existing chromatin landscapes. Specific genes with altered expression show differential interactions with dysregulated cis-regulatory elements. Collectively, we find that cell-type specific targets of SETDB1 maintain cellular identities by modulating CTCF binding, which shape nuclear architecture and transcriptomic networks.

Transcriptional regulation of FACT involves Coordination of chromatin accessibility and CTCF binding

Histone chaperone FACT (facilitates chromatin transcription) is well known to promote chromatin recovery during transcription. However, the mechanism how FACT regulates genome-wide chromatin accessibility and transcription factor binding has not been fully elucidated. Through loss-of-function studies, we show here that FACT component Ssrp1 is required for DNA replication and DNA damage repair and is also essential for progression of cell phase transition and cell proliferation in mouse embryonic fibroblast cells. On the molecular level, absence of the Ssrp1 leads to increased chromatin accessibility, enhanced CTCF binding, and a remarkable change in dynamic range of gene expression. Our study thus unequivocally uncovers a unique mechanism by which FACT complex regulates transcription by coordinating genome-wide chromatin accessibility and CTCF binding.

TET1 regulates responses to house dust mite by altering chromatin accessibility, DNA methylation, and gene expression in airway epithelial cells

Background

Previous studies have identified TET1 as a potential key regulator of genes linked to asthma. TET1 has been shown to transcriptionally respond to house dust mite extract, an allergen known to directly cause allergic asthma development, and regulate the expression of genes involved in asthma. How TET1 regulates expression of these genes, however, is unknown. TET1 is a DNA demethylase; therefore, most prior research on TET1-based gene regulation has focused on how TET1 affects methylation. However, TET1 can also interact directly with transcription factors and histone modifiers to regulate gene expression. Understanding how TET1 regulates expression to contribute to allergic responses and asthma development thus requires a comprehensive approach. To this end, we measured mRNA expression, DNA methylation, chromatin accessibility and histone modifications in control and TET1 knockdown human bronchial epithelial cells treated or untreated with house dust mite extract.

Results

Throughout our analyses, we detected strong similarities between the effects of TET1 knockdown alone and the effects of HDM treatment alone. One especially striking pattern was that both TET1 knockdown and HDM treatment generally led to decreased chromatin accessibility at largely the same genomic loci. Transcription factor enrichment analyses indicated that altered chromatin accessibility following the loss of TET1 may affect, or be affected by, CTCF and CEBP binding. TET1 loss also led to changes in DNA methylation, but these changes were generally in regions where accessibility was not changing.

Conclusions

TET1 regulates gene expression through different mechanisms (DNA methylation and chromatin accessibility) in different parts of the genome in the airway epithelial cells, which mediates inflammatory responses to allergen. Collectively, our data suggest novel molecular mechanisms through which TET1 regulates critical pathways following allergen challenges and contributes to the development of asthma.

DNA methylation modulated genetic variant effect on gene transcriptional regulation

BACKGROUND

Expression quantitative trait locus (eQTL) analysis has emerged as an important tool in elucidating the link between genetic variants and gene expression, thereby bridging the gap between risk SNPs and associated diseases. We recently identified and validated a specific case where the methylation of a CpG site influences the relationship between the genetic variant and gene expression.

RESULTS

Here, to systematically evaluate this regulatory mechanism, we develop an extended eQTL mapping method, termed DNA methylation modulated eQTL (memo-eQTL). Applying this memo-eQTL mapping method to 128 normal prostate samples enables identification of 1063 memo-eQTLs, the majority of which are not recognized as conventional eQTLs in the same cohort. We observe that the methylation of the memo-eQTL CpG sites can either enhance or insulate the interaction between SNP and gene expression by altering CTCF-based chromatin 3D structure.

CONCLUSIONS

This study demonstrates the prevalence of memo-eQTLs paving the way to identify novel causal genes for traits or diseases associated with genetic variations.

Single-cell bisulfite-free 5mC and 5hmC sequencing with high sensitivity and scalability

Existing single-cell bisulfite-based DNA methylation analysis is limited by low DNA recovery, and the measurement of 5hmC at single-base resolution remains challenging. Here, we present a bisulfite-free single-cell whole-genome 5mC and 5hmC profiling technique, named Cabernet, which can characterize 5mC and 5hmC at single-base resolution with high genomic coverage. Cabernet utilizes Tn5 transposome for DNA fragmentation, which enables the discrimination between different alleles for measuring hemi-methylation status. Using Cabernet, we revealed the 5mC, hemi-5mC and 5hmC dynamics during early mouse embryo development, uncovering genomic regions exclusively governed by active or passive demethylation. We show that hemi-methylation status can be used to distinguish between pre- and post-replication cells, enabling more efficient cell grouping when integrated with 5mC profiles. The property of Tn5 naturally enables Cabernet to achieve high-throughput single-cell methylome profiling, where we probed mouse cortical neurons and embryonic day 7.5 (E7.5) embryos, and constructed the library for thousands of single cells at high efficiency, demonstrating its potential for analyzing complex tissues at substantially low cost. Together, we present a way of high-throughput methylome and hydroxymethylome detection at single-cell resolution, enabling efficient analysis of the epigenetic status of biological systems with complicated nature such as neurons and cancer cells.

Defining a core configuration for human centromeres during mitosis

The centromere components cohesin, CENP-A, and centromeric DNA are essential for biorientation of sister chromatids on the mitotic spindle and accurate sister chromatid segregation. Insight into the 3D organization of centromere components would help resolve how centromeres function on the mitotic spindle. We use ChIP-seq and super-resolution microscopy with single particle averaging to examine the geometry of essential centromeric components on human chromosomes. Both modalities suggest cohesin is enriched at pericentromeric DNA. CENP-A localizes to a subset of the α-satellite DNA, with clusters separated by ~562 nm and a perpendicular intervening ~190 nM wide axis of cohesin in metaphase chromosomes. Differently sized α-satellite arrays achieve a similar core structure. Here we present a working model for a common core configuration of essential centromeric components that includes CENP-A nucleosomes, α-satellite DNA and pericentromeric cohesion. This configuration helps reconcile how centromeres function and serves as a foundation to add components of the chromosome segregation machinery.

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

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

Hypomethylation and Overexpression of Th17-Associated Genes is a Hallmark of Intestinal CD4+ Lymphocytes in Crohn's Disease

BACKGROUND

The development of Crohn's disease [CD] involves immune cell signalling pathways regulated by epigenetic modifications. Aberrant DNA methylation has been identified in peripheral blood and bulk intestinal tissue from CD patients. However, the DNA methylome of disease-associated intestinal CD4+ lymphocytes has not been evaluated.

MATERIALS AND METHODS

Genome-wide DNA methylation sequencing was performed from terminal ileum CD4+ cells from 21 CD patients and 12 age- and sex-matched controls. Data were analysed for differentially methylated CpGs [DMCs] and methylated regions [DMRs]. Integration was performed with RNA-sequencing data to evaluate the functional impact of DNA methylation changes on gene expression. DMRs were overlapped with regions of differentially open chromatin [by ATAC-seq] and CCCTC-binding factor [CTCF] binding sites [by ChIP-seq] between peripherally derived Th17 and Treg cells.

RESULTS

CD4+ cells in CD patients had significantly increased DNA methylation compared to those from the controls. A total of 119 051 DMCs and 8113 DMRs were detected. While hypermethylated genes were mostly related to cell metabolism and homeostasis, hypomethylated genes were significantly enriched within the Th17 signalling pathway. The differentially enriched ATAC regions in Th17 cells [compared to Tregs] were hypomethylated in CD patients, suggesting heightened Th17 activity. There was significant overlap between hypomethylated DNA regions and CTCF-associated binding sites.

CONCLUSIONS

The methylome of CD patients shows an overall dominant hypermethylation yet hypomethylation is more concentrated in proinflammatory pathways, including Th17 differentiation. Hypomethylation of Th17-related genes associated with areas of open chromatin and CTCF binding sites constitutes a hallmark of CD-associated intestinal CD4+ cells.

The Molecular Impacts of Retrotransposons in Development and Diseases

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

eRNA co-expression network uncovers TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced by transcriptional enhancers that are highly correlated with their activity. Using a capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across 67 individuals, we identified inter-individual variation in the expression of over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers are associated with each other and with promoters. Mid- to long-range co-expression showed a distance-dependent decay that was modified by TF occupancy. In particular, we found a class of "bivalent" TFs, including Cohesin, that both facilitate and isolate the interaction between enhancers and/or promoters, depending on their topology. At short distances, we observed strand-specific correlations between nearby eRNAs in both convergent and divergent orientations. Our results support a cooperative model of convergent eRNAs, consistent with eRNAs facilitating adjacent enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions as a function of distance, orientation, and binding landscapes of TFs.

CTCF coordinates cell fate specification via orchestrating regulatory hubs with pioneer transcription factors

CCCTC-binding factor (CTCF), a ubiquitously expressed architectural protein, has emerged as a key regulator of cell identity gene transcription. However, the precise molecular mechanism underlying specialized functions of CTCF remains elusive. Here, we investigate the mechanism through integrative analyses of primary hepatocytes, myocytes, and B cells from mouse and human. We demonstrate that CTCF cooperates with lineage-specific pioneer transcription factors (TFs), including MyoD, FOXA, and PU.1, to control cell identity at 1D and 3D levels. At the 1D level, pioneer TFs facilitate lineage-specific CTCF occupancy via opening chromatin. At the 3D level, CTCF and pioneer TFs form regulatory hubs to govern the expression of cell identity genes. This mechanism is validated using MyoD-null mice, CTCF knockout mice, and CRISPR editing during myogenic differentiation. Collectively, these findings uncover a general mechanism whereby CTCF acts as a cell identity cofactor to control cell identity genes via orchestrating regulatory hubs with pioneer TFs.

Decitabine disrupts EBV genomic epiallele DNA methylation patterns around CTCF binding sites to increase chromatin accessibility and lytic transcription in gastric cancer

IMPORTANCE

Epstein-Barr virus (EBV) latency is controlled by epigenetic silencing by DNA methylation [5-methyl cytosine (5mC)], histone modifications, and chromatin looping. However, how they dictate the transcriptional program in EBV-associated gastric cancers remains incompletely understood. EBV-associated gastric cancer displays a 5mC hypermethylated phenotype. A potential treatment for this cancer subtype is the DNA hypomethylating agent, which induces EBV lytic reactivation and targets hypermethylation of the cellular DNA. In this study, we identified a heterogeneous pool of EBV epialleles within two tumor-derived gastric cancer cell lines that are disrupted with a hypomethylating agent. Stochastic DNA methylation patterning at critical regulatory regions may be an underlying mechanism for spontaneous reactivation. Our results highlight the critical role of epigenetic modulation on EBV latency and life cycle, which is maintained through the interaction between 5mC and the host protein CCCTC-binding factor.

Molecular Mechanisms for the Regulation of Nuclear Membrane Integrity

The nuclear membrane serves a critical role in protecting the contents of the nucleus and facilitating material and signal exchange between the nucleus and cytoplasm. While extensive research has been dedicated to topics such as nuclear membrane assembly and disassembly during cell division, as well as interactions between nuclear transmembrane proteins and both nucleoskeletal and cytoskeletal components, there has been comparatively less emphasis on exploring the regulation of nuclear morphology through nuclear membrane integrity. In particular, the role of type II integral proteins, which also function as transcription factors, within the nuclear membrane remains an area of research that is yet to be fully explored. The integrity of the nuclear membrane is pivotal not only during cell division but also in the regulation of gene expression and the communication between the nucleus and cytoplasm. Importantly, it plays a significant role in the development of various diseases. This review paper seeks to illuminate the biomolecules responsible for maintaining the integrity of the nuclear membrane. It will delve into the mechanisms that influence nuclear membrane integrity and provide insights into the role of type II membrane protein transcription factors in this context. Understanding these aspects is of utmost importance, as it can offer valuable insights into the intricate processes governing nuclear membrane integrity. Such insights have broad-reaching implications for cellular function and our understanding of disease pathogenesis.

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.

G-quadruplexes associated with R-loops promote CTCF binding

CTCF is a critical regulator of genome architecture and gene expression that binds thousands of sites on chromatin. CTCF genomic localization is controlled by the recognition of a DNA sequence motif and regulated by DNA modifications. However, CTCF does not bind to all its potential sites in all cell types, raising the question of whether the underlying chromatin structure can regulate CTCF occupancy. Here, we report that R-loops facilitate CTCF binding through the formation of associated G-quadruplex (G4) structures. R-loops and G4s co-localize with CTCF at many genomic regions in mouse embryonic stem cells and promote CTCF binding to its cognate DNA motif in vitro. R-loop attenuation reduces CTCF binding in vivo. Deletion of a specific G4-forming motif in a gene reduces CTCF binding and alters gene expression. Conversely, chemical stabilization of G4s results in CTCF gains and accompanying alterations in chromatin organization, suggesting a pivotal role for G4 structures in reinforcing long-range genome interactions through CTCF.

Allelic chromatin structure precedes imprinted expression of Kcnk9 during neurogenesis

Differences in chromatin state inherited from the parental gametes influence the regulation of maternal and paternal alleles in offspring. This phenomenon, known as genomic imprinting, results in genes preferentially transcribed from one parental allele. While local epigenetic factors such as DNA methylation are known to be important for the establishment of imprinted gene expression, less is known about the mechanisms by which differentially methylated regions (DMRs) lead to differences in allelic expression across broad stretches of chromatin. Allele-specific higher-order chromatin structure has been observed at multiple imprinted loci, consistent with the observation of allelic binding of the chromatin-organizing factor CTCF at multiple DMRs. However, whether allelic chromatin structure impacts allelic gene expression is not known for most imprinted loci. Here we characterize the mechanisms underlying brain-specific imprinted expression of the Peg13-Kcnk9 locus, an imprinted region associated with intellectual disability. We performed region capture Hi-C on mouse brains from reciprocal hybrid crosses and found imprinted higher-order chromatin structure caused by the allelic binding of CTCF to the Peg13 DMR. Using an in vitro neuron differentiation system, we showed that imprinted chromatin structure precedes imprinted expression at the locus. Additionally, activation of a distal enhancer induced imprinted expression of Kcnk9 in an allelic chromatin structure-dependent manner. This work provides a high-resolution map of imprinted chromatin structure and demonstrates that chromatin state established in early development can promote imprinted expression upon differentiation.

Extensive evaluation of DNA methylation of functional elements in the murine Fkbp5 locus using high-accuracy DNA methylation measurement via targeted bisulfite sequencing

FKBP5 is an important stress-regulatory gene implicated in stress-related psychiatric diseases. Single nucleotide polymorphisms of the FKBP5 gene were shown to interact with early life stress to alter the glucocorticoid-related stress response and moderate disease risk. Demethylation of cytosine-phosphate-guanine-dinucleotides (CpGs) in regulatory glucocorticoid-responsive elements was suggested to be the mediating epigenetic mechanism for long-term stress effects, but studies on Fkbp5 DNA methylation (DNAm) in rodents are so far limited. We evaluated the applicability of high-accuracy DNA methylation measurement via targeted bisulfite sequencing (HAM-TBS), a next-generation sequencing-based technology, to allow a more in-depth characterisation of the DNA methylation of the murine Fkbp5 locus in three different tissues (blood, frontal cortex and hippocampus). In this study, we not only increased the number of evaluated sites in previously described regulatory regions (in introns 1 and 5), but also extended the evaluation to novel, possibly relevant regulatory regions of the gene (in intron 8, the transcriptional start site, the proximal enhancer and CTCF-binding sites within the 5'UTR). We here describe the assessment of HAM-TBS assays for a panel of 157 CpGs with possible functional relevance in the murine Fkbp5 gene. DNAm profiles were tissue-specific, with lesser differences between the two brain regions than between the brain and blood. Moreover, we identified DNAm changes in the Fkbp5 locus after early life stress exposure in the frontal cortex and blood. Our findings indicate that HAM-TBS is a valuable tool for broader exploration of the DNAm of the murine Fkbp5 locus and its involvement in the stress response.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes

BACKGROUND

Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk.

RESULTS

We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites).

CONCLUSIONS

Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

The Role of FAS Receptor Methylation in Osteosarcoma Metastasis

Osteosarcoma is the most frequent primary malignant bone tumor with an annual incidence of about 400 cases in the United States. Osteosarcoma primarily metastasizes to the lungs, where FAS ligand (FASL) is constitutively expressed. The interaction of FASL and its cell surface receptor, FAS, triggers apoptosis in normal cells; however, this function is altered in cancer cells. DNA methylation has previously been explored as a mechanism for altering FAS expression, but no variability was identified in the CpG island (CGI) overlapping the promoter. Analysis of an expanded region, including CGI shores and shelves, revealed high variability in the methylation of certain CpG sites that correlated significantly with FAS mRNA expression in a negative manner. Bisulfite sequencing revealed additional CpG sites, which were highly methylated in the metastatic LM7 cell line but unmethylated in its parental non-metastatic SaOS-2 cell line. Treatment with the demethylating agent, 5-azacytidine, resulted in a loss of methylation in CpG sites located within the FAS promoter and restored FAS protein expression in LM7 cells, resulting in reduced migration. Orthotopic implantation of 5-azacytidine treated LM7 cells into severe combined immunodeficient mice led to decreased lung metastases. These results suggest that DNA methylation of CGI shore sites may regulate FAS expression and constitute a potential target for osteosarcoma therapy, utilizing demethylating agents currently approved for the treatment of other cancers.

Genome-wide DNA methylation and transcriptome analyses reveal the key gene for wool type variation in sheep

BACKGROUND

Wool fibers are valuable materials for textile industry. Typical wool fibers are divided into medullated and non-medullated types, with the former generated from primary wool follicles and the latter by either primary or secondary wool follicles. The medullated wool is a common wool type in the ancestors of fine wool sheep before breeding. The fine wool sheep have a non-medullated coat. However, the critical period determining the type of wool follicles is the embryonic stage, which limits the phenotypic observation and variant contrast, making both selection and studies of wool type variation fairly difficult.

RESULTS

During the breeding of a modern fine (MF) wool sheep population with multiple-ovulation and embryo transfer technique, we serendipitously discovered lambs with ancestral-like coarse (ALC) wool. Whole-genome resequencing confirmed ALC wool lambs as a variant type from the MF wool population. We mapped the significantly associated methylation locus on chromosome 4 by using whole genome bisulfite sequencing signals, and in turn identified the SOSTDC1 gene as exons hypermethylated in ALC wool lambs compare to their half/full sibling MF wool lambs. Transcriptome sequencing found that SOSTDC1 was expressed dozens of times more in ALC wool lamb skin than that of MF and was at the top of all differentially expressed genes. An analogy with the transcriptome of coarse/fine wool breeds revealed that differentially expressed genes and enriched pathways at postnatal lamb stage in ALC/MF were highly similar to those at the embryonic stage in the former. Further experiments validated that the SOSTDC1 gene was specifically highly expressed in the nucleus of the dermal papilla of primary wool follicles.

CONCLUSION

In this study, we conducted genome-wide differential methylation site association analysis on differential wool type trait, and located the only CpG locus that strongly associated with primary wool follicle development. Combined with transcriptome analysis, SOSTDC1 was identified as the only gene at this locus that was specifically overexpressed in the primary wool follicle stem cells of ALC wool lamb skin. The discovery of this key gene and its epigenetic regulation contributes to understanding the domestication and breeding of fine wool sheep.

DNA strand asymmetry generated by CpG hemimethylation has opposing effects on CTCF binding

CpG methylation generally occurs on both DNA strands and is essential for mammalian development and differentiation. Until recently, hemimethylation, in which only one strand is methylated, was considered to be simply a transitory state generated during DNA synthesis. The discovery that a subset of CCCTC-binding factor (CTCF) binding sites is heritably hemimethylated suggests that hemimethylation might have an unknown biological function. Here we show that the binding of CTCF is profoundly altered by which DNA strand is methylated and by the specific CTCF binding motif. CpG methylation on the motif strand can inhibit CTCF binding by up to 7-fold, whereas methylation on the opposite strand can stimulate binding by up to 4-fold. Thus, hemimethylation can alter binding by up to 28-fold in a strand-specific manner. The mechanism for sensing methylation on the opposite strand requires two critical residues, V454 and S364, within CTCF zinc fingers 7 and 4. Similar to methylation, CpG hydroxymethylation on the motif strand can inhibit CTCF binding by up to 4-fold. However, hydroxymethylation on the opposite strand removes the stimulatory effect. Strand-specific methylation states may therefore provide a mechanism to explain the transient and dynamic nature of CTCF-mediated chromatin interactions.

Loss of PBAF promotes expansion and effector differentiation of CD8+ T cells during chronic viral infection and cancer

During chronic viral infection and cancer, it has been established that a subset of progenitor CD8+ T cells continuously gives rise to terminally exhausted cells and cytotoxic effector cells. Although multiple transcriptional programs governing the bifurcated differentiation trajectories have been previously studied, little is known about the chromatin structure changes regulating CD8+ T cell-fate decision. In this study, we demonstrate that the chromatin remodeling complex PBAF restrains expansion and promotes exhaustion of CD8+ T cells during chronic viral infection and cancer. Mechanistically, transcriptomic and epigenomic analyses reveal the role of PBAF in maintaining chromatin accessibility of multiple genetic pathways and transcriptional programs to restrain proliferation and promote T cell exhaustion. Harnessing this knowledge, we demonstrate that perturbation of PBAF complex constrained exhaustion and promoted expansion of tumor-specific CD8+ T cells resulting in antitumor immunity in a preclinical melanoma model, implicating PBAF as an attractive target for cancer immunotherapeutic.

An epigenetic signature of advanced colorectal cancer metastasis

Colorectal cancer (CRC) is a leading cause of morbidity and mortality worldwide. The majority of CRC deaths are caused by tumor metastasis, even following treatment. There is strong evidence for epigenetic changes, such as DNA methylation, accompanying CRC metastasis and poorer patient survival. Earlier detection and a better understanding of molecular drivers for CRC metastasis are of critical clinical importance. Here, we identify a signature of advanced CRC metastasis by performing whole genome-scale DNA methylation and full transcriptome analyses of paired primary cancers and liver metastases from CRC patients. We observed striking methylation differences between primary and metastatic pairs. A subset of loci showed coordinated methylation-expression changes, suggesting these are potentially epigenetic drivers that control the expression of critical genes in the metastatic cascade. The identification of CRC epigenomic markers of metastasis has the potential to enable better outcome prediction and lead to the discovery of new therapeutic targets.

Allelic chromatin structure primes imprinted expression of Kcnk9 during neurogenesis

Differences in chromatin state inherited from the parental gametes influence the regulation of maternal and paternal alleles in offspring. This phenomenon, known as genomic imprinting, results in genes preferentially transcribed from one parental allele. While local epigenetic factors such as DNA methylation are known to be important for the establishment of imprinted gene expression, less is known about the mechanisms by which differentially methylated regions (DMRs) lead to differences in allelic expression across broad stretches of chromatin. Allele-specific higher-order chromatin structure has been observed at multiple imprinted loci, consistent with the observation of allelic binding of the chromatin-organizing factor CTCF at multiple DMRs. However, whether allelic chromatin structure impacts allelic gene expression is not known for most imprinted loci. Here we characterize the mechanisms underlying brain-specific imprinted expression of the Peg13-Kcnk9 locus, an imprinted region associated with intellectual disability. We performed region capture Hi-C on mouse brain from reciprocal hybrid crosses and found imprinted higher-order chromatin structure caused by the allelic binding of CTCF to the Peg13 DMR. Using an in vitro neuron differentiation system, we show that on the maternal allele enhancer-promoter contacts formed early in development prime the brain-specific potassium leak channel Kcnk9 for maternal expression prior to neurogenesis. In contrast, these enhancer-promoter contacts are blocked by CTCF on the paternal allele, preventing paternal Kcnk9 activation. This work provides a high-resolution map of imprinted chromatin structure and demonstrates that chromatin state established in early development can promote imprinted expression upon differentiation.

Benefiting from the intrinsic role of epigenetics to predict patterns of CTCF binding

Motivation

One of the most relevant mechanisms involved in the determination of chromatin structure is the formation of structural loops that are also related with the conservation of chromatin states. Many of these loops are stabilized by CCCTC-binding factor (CTCF) proteins at their base. Despite the relevance of chromatin structure and the key role of CTCF, the role of the epigenetic factors that are involved in the regulation of CTCF binding, and thus, in the formation of structural loops in the chromatin, is not thoroughly understood.

Results

Here we describe a CTCF binding predictor based on Random Forest that employs different epigenetic data and genomic features. Importantly, given the ability of Random Forests to determine the relevance of features for the prediction, our approach also shows how the different types of descriptors impact the binding of CTCF, confirming previous knowledge on the relevance of chromatin accessibility and DNA methylation, but demonstrating the effect of epigenetic modifications on the activity of CTCF. We compared our approach against other predictors and found improved performance in terms of areas under PR and ROC curves (PRAUC-ROCAUC), outperforming current state-of-the-art methods.

Defining a core configuration for human centromeres during mitosis

The biorientation of sister chromatids on the mitotic spindle, essential for accurate sister chromatid segregation, relies on critical centromere components including cohesin, the centromere-specific H3 variant CENP-A, and centromeric DNA. Centromeric DNA is highly variable between chromosomes yet must accomplish a similar function. Moreover, how the 50 nm cohesin ring, proposed to encircle sister chromatids, accommodates inter-sister centromeric distances of hundreds of nanometers on the metaphase spindle is a conundrum. Insight into the 3D organization of centromere components would help resolve how centromeres function on the mitotic spindle. We used ChIP-seq and super-resolution microscopy to examine the geometry of essential centromeric components on human chromosomes. ChIP-seq demonstrates that cohesin subunits are depleted in α-satellite arrays where CENP-A nucleosomes and kinetochores assemble. Cohesin is instead enriched at pericentromeric DNA. Structured illumination microscopy of sister centromeres is consistent, revealing a non-overlapping pattern of CENP-A and cohesin. We used single particle averaging of hundreds of mitotic sister chromatids to develop an average centromere model. CENP-A clusters on sister chromatids, connected by α-satellite, are separated by ~562 nm with a perpendicular intervening ~190 nM wide axis of cohesin. Two differently sized α-satellite arrays on chromosome 7 display similar inter-sister CENP-A cluster distance, demonstrating different sized arrays can achieve a common spacing. Our data suggest a working model for a common core configuration of essential centromeric components that includes CENP-A nucleosomes at the outer edge of extensible α-satellite DNA and pericentromeric cohesion. This configuration helps reconcile how centromeres function and serves as a foundation for future studies of additional components required for centromere function.

Targeted DNA Demethylation: Vectors, Effectors and Perspectives

Aberrant DNA hypermethylation at regulatory cis-elements of particular genes is seen in a plethora of pathological conditions including cardiovascular, neurological, immunological, gastrointestinal and renal diseases, as well as in cancer, diabetes and others. Thus, approaches for experimental and therapeutic DNA demethylation have a great potential to demonstrate mechanistic importance, and even causality of epigenetic alterations, and may open novel avenues to epigenetic cures. However, existing methods based on DNA methyltransferase inhibitors that elicit genome-wide demethylation are not suitable for treatment of diseases with specific epimutations and provide a limited experimental value. Therefore, gene-specific epigenetic editing is a critical approach for epigenetic re-activation of silenced genes. Site-specific demethylation can be achieved by utilizing sequence-dependent DNA-binding molecules such as zinc finger protein array (ZFA), transcription activator-like effector (TALE) and clustered regularly interspaced short palindromic repeat-associated dead Cas9 (CRISPR/dCas9). Synthetic proteins, where these DNA-binding domains are fused with the DNA demethylases such as ten-eleven translocation (Tet) and thymine DNA glycosylase (TDG) enzymes, successfully induced or enhanced transcriptional responsiveness at targeted loci. However, a number of challenges, including the dependence on transgenesis for delivery of the fusion constructs, remain issues to be solved. In this review, we detail current and potential approaches to gene-specific DNA demethylation as a novel epigenetic editing-based therapeutic strategy.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.

DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation

BACKGROUND

Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions.

RESULTS

The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression.

CONCLUSIONS

In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.