Widespread plasticity in CTCF occupancy linked to DNA methylation

DNA methylation dynamics: identification and functional annotation

DNA methylation is an epigenetic modification of cytosines that undergoes dynamic changes in a temporal, spatial and cell-type-specific manner. Recent advances in technology have permitted the profiling of high-throughput methylomes in large numbers of biological samples. Various computational tools have been developed to identify and analyze DNA methylation dynamics in a variety of critical biological processes. As DNA methylation is becoming increasingly viewed as a dynamic process, the mechanisms governing DNA methylation dynamics and its roles in the transcriptional regulatory network are of great interest. It has been reported that DNA methylation dynamics plays essential roles in multiple biological processes, including development and cancer. As a functional event, the dynamics of DNA methylation have become increasingly relevant to many researchers. Here, we review state-of-the-art advances at three levels (genome-wide identification, regulatory mechanism investigation and the functional annotation) in the field of DNA methylation dynamics, as well as the future perspective of DNA methylation dynamics.

The Evolution of Epigenetics: From Prokaryotes to Humans and Its Biological Consequences

The evolution process includes genetic alterations that started with prokaryotes and now continues in humans. A distinct difference between prokaryotic chromosomes and eukaryotic chromosomes involves histones. As evolution progressed, genetic alterations accumulated and a mechanism for gene selection developed. It was as if nature was experimenting to optimally utilize the gene pool without changing individual gene sequences. This mechanism is called epigenetics, as it is above the genome. Curiously, the mechanism of epigenetic regulation in prokaryotes is strikingly different from that in eukaryotes, mainly higher eukaryotes, like mammals. In fact, epigenetics plays a significant role in the conserved process of embryogenesis and human development. Malfunction of epigenetic regulation results in many types of undesirable effects, including cardiovascular disease, metabolic disorders, autoimmune diseases, and cancer. This review provides a comparative analysis and new insights into these aspects.

The Many Roles of BAF (mSWI/SNF) and PBAF Complexes in Cancer

During the last decade, a host of epigenetic mechanisms were found to contribute to cancer and other human diseases. Several genomic studies have revealed that ∼20% of malignancies have alterations of the subunits of polymorphic BRG-/BRM-associated factor (BAF) and Polybromo-associated BAF (PBAF) complexes, making them among the most frequently mutated complexes in cancer. Recurrent mutations arise in genes encoding several BAF/PBAF subunits, including ARID1A, ARID2, PBRM1, SMARCA4, and SMARCB1 These subunits share some degree of conservation with subunits from related adenosine triphosphate (ATP)-dependent chromatin remodeling complexes in model organisms, in which a large body of work provides insight into their roles in cancer. Here, we review the roles of BAF- and PBAF-like complexes in these organisms, and relate these findings to recent discoveries in cancer epigenomics. We review several roles of BAF and PBAF complexes in cancer, including transcriptional regulation, DNA repair, and regulation of chromatin architecture and topology. More recent results highlight the need for new techniques to study these complexes.

CRISPR/Cas9 genome editing throws descriptive 3-D genome folding studies for a loop

CRISPR/Cas9 genome editing studies have recently shed new light into the causal link between the linear DNA sequence and 3-D chromatin architecture. Here we describe current models for the folding of genomes into a nested hierarchy of higher-order structures and discuss new insights into the organizing principles governing genome folding at each length scale. WIREs Syst Biol Med 2016, 8:286-299. doi: 10.1002/wsbm.1338 For further resources related to this article, please visit the WIREs website.

Next-generation sequencing identifies major DNA methylation changes during progression of Ph+ chronic myeloid leukemia

Little is known about the impact of DNA methylation on the evolution/progression of Ph+ chronic myeloid leukemia (CML). We investigated the methylome of CML patients in chronic phase (CP-CML), accelerated phase (AP-CML) and blast crisis (BC-CML) as well as in controls by reduced representation bisulfite sequencing. Although only ~600 differentially methylated CpG sites were identified in samples obtained from CP-CML patients compared with controls, ~6500 differentially methylated CpG sites were found in samples from BC-CML patients. In the majority of affected CpG sites, methylation was increased. In CP-CML patients who progressed to AP-CML/BC-CML, we identified up to 897 genes that were methylated at the time of progression but not at the time of diagnosis. Using RNA-sequencing, we observed downregulated expression of many of these genes in BC-CML compared with CP-CML samples. Several of them are well-known tumor-suppressor genes or regulators of cell proliferation, and gene re-expression was observed by the use of epigenetic active drugs. Together, our results demonstrate that CpG site methylation clearly increases during CML progression and that it may provide a useful basis for revealing new targets of therapy in advanced CML.

A Comparative Analysis of Genetic and Epigenetic Events of Breast and Ovarian Cancer Related to Tumorigenesis

Breast cancer persists as the most common cause of cancer death in women worldwide. Ovarian cancer is also a significant source of morbidity and mortality, as the fifth leading cause of cancer death among women. This reflects the continued need for further understanding and innovation in cancer treatment. Though breast and ovarian cancer usually present as distinct clinical entities, the recent explosion of large-scale -omics research has uncovered many overlaps, particularly with respect to genetic and epigenetic alterations. We compared genetic, microenvironmental, stromal, and epigenetic changes common between breast and ovarian cancer cells, as well as the clinical relevance of these changes. Some of the most striking commonalities include genetic alterations of BRCA1 and 2, TP53, RB1, NF1, FAT3, MYC, PTEN, and PIK3CA; down regulation of miRNAs 9, 100, 125a, 125b, and 214; and epigenetic alterations such as H3K27me3, H3K9me2, H3K9me3, H4K20me3, and H3K4me. These parallels suggest shared features of pathogenesis. Furthermore, preliminary evidence suggests a shared epigenetic mechanism of oncogenesis. These similarities, warrant further investigation in order to ultimately inform development of more effective chemotherapeutics, as well as strategies to circumvent drug resistance.

Differential DNA methylation identified in the blood and retina of AMD patients

Age-related macular degeneration (AMD) is a major cause of blindness in the western world. While genetic studies have linked both common and rare variants in genes involved in regulation of the complement system to increased risk of development of AMD, environmental factors, such as smoking and nutrition, can also significantly affect the risk of developing the disease and the rate of disease progression. Since epigenetics has been implicated in mediating, in part, the disease risk associated with some environmental factors, we investigated a possible epigenetic contribution to AMD. We performed genome-wide DNA methylation profiling of blood from AMD patients and controls. No differential methylation site reached genome-wide significance; however, when epigenetic changes in and around known GWAS-defined AMD risk loci were explored, we found small but significant DNA methylation differences in the blood of neovascular AMD patients near age-related maculopathy susceptibility 2 (ARMS2), a top-ranked GWAS locus preferentially associated with neovascular AMD. The methylation level of one of the CpG sites significantly correlated with the genotype of the risk SNP rs10490924, suggesting a possible epigenetic mechanism of risk. Integrating genome-wide DNA methylation analysis of retina samples with and without AMD together with blood samples, we further identified a consistent, replicable change in DNA methylation in the promoter region of protease serine 50 (PRSS50). These methylation changes may identify sites in novel genes that are susceptible to non-genetic factors known to contribute to AMD development and progression.

Reciprocal Regulation of the Cardiac Epigenome by Chromatin Structural Proteins Hmgb and Ctcf: IMPLICATIONS FOR TRANSCRIPTIONAL REGULATION

Transcriptome remodeling in heart disease occurs through the coordinated actions of transcription factors, histone modifications, and other chromatin features at pathology-associated genes. The extent to which genome-wide chromatin reorganization also contributes to the resultant changes in gene expression remains unknown. We examined the roles of two chromatin structural proteins, Ctcf (CCCTC-binding factor) and Hmgb2 (high mobility group protein B2), in regulating pathologic transcription and chromatin remodeling. Our data demonstrate a reciprocal relationship between Hmgb2 and Ctcf in controlling aspects of chromatin structure and gene expression. Both proteins regulate each others' expression as well as transcription in cardiac myocytes; however, only Hmgb2 does so in a manner that involves global reprogramming of chromatin accessibility. We demonstrate that the actions of Hmgb2 on local chromatin accessibility are conserved across genomic loci, whereas the effects on transcription are loci-dependent and emerge in concert with histone modification and other chromatin features. Finally, although both proteins share gene targets, Hmgb2 and Ctcf, neither binds these genes simultaneously nor do they physically colocalize in myocyte nuclei. Our study uncovers a previously unknown relationship between these two ubiquitous chromatin proteins and provides a mechanistic explanation for how Hmgb2 regulates gene expression and cellular phenotype. Furthermore, we provide direct evidence for structural remodeling of chromatin on a genome-wide scale in the setting of cardiac disease.

An ancient yet flexible cis-regulatory architecture allows localized Hedgehog tuning by patched/Ptch1

The Hedgehog signaling pathway is part of the ancient developmental-evolutionary animal toolkit. Frequently co-opted to pattern new structures, the pathway is conserved among eumetazoans yet flexible and pleiotropic in its effects. The Hedgehog receptor, Patched, is transcriptionally activated by Hedgehog, providing essential negative feedback in all tissues. Our locus-wide dissections of the cis-regulatory landscapes of fly patched and mouse Ptch1 reveal abundant, diverse enhancers with stage- and tissue-specific expression patterns. The seemingly simple, constitutive Hedgehog response of patched/Ptch1 is driven by a complex regulatory architecture, with batteries of context-specific enhancers engaged in promoter-specific interactions to tune signaling individually in each tissue, without disturbing patterning elsewhere. This structure—one of the oldest cis-regulatory features discovered in animal genomes—explains how patched/Ptch1 can drive dramatic adaptations in animal morphology while maintaining its essential core function. It may also suggest a general model for the evolutionary flexibility of conserved regulators and pathways.

DOI:http://dx.doi.org/10.7554/eLife.13550.001

Epigenetic silencing of tumor suppressor genes: Paradigms, puzzles, and potential

Cancer constitutes a set of diseases with heterogeneous molecular pathologies. However, there are a number of universal aberrations common to all cancers, one of these being the epigenetic silencing of tumor suppressor genes (TSGs). The silencing of TSGs is thought to be an early, driving event in the oncogenic process. With this in consideration, great efforts have been made to develop small molecules aimed at the restoration of TSGs in order to limit tumor cell proliferation and survival. However, the molecular forces that drive the broad epigenetic reprogramming and transcriptional repression of these genes remain ill-defined. Undoubtedly, understanding the molecular underpinnings of transcriptionally silenced TSGs will aid us in our ability to reactivate these key anti-cancer targets. Here, we describe what we consider to be the five most logical molecular mechanisms that may account for this widely observed phenomenon: 1) ablation of transcription factor binding, 2) overexpression of DNA methyltransferases, 3) disruption of CTCF binding, 4) elevation of EZH2 activity, 5) aberrant expression of long non-coding RNAs. The strengths and weaknesses of each proposed mechanism is highlighted, followed by an overview of clinical efforts to target these processes.

CRISPR Double Cutting through the Labyrinthine Architecture of 3D Genomes

The genomes are organized into ordered and hierarchical topological structures in interphase nuclei. Within discrete territories of each chromosome, topologically associated domains (TADs) play important roles in various nuclear processes such as gene regulation. Inside TADs separated by relatively constitutive boundaries, distal elements regulate their gene targets through specific chromatin-looping contacts such as long-distance enhancer-promoter interactions. High-throughput sequencing studies have revealed millions of potential regulatory DNA elements, which are much more abundant than the mere ∼20,000 genes they control. The recently emerged CRISPR-Cas9 genome editing technologies have enabled efficient and precise genetic and epigenetic manipulations of genomes. The multiplexed and high-throughput CRISPR capabilities facilitate the discovery and dissection of gene regulatory elements. Here, we describe the applications of CRISPR for genome, epigenome, and 3D genome editing, focusing on CRISPR DNA-fragment editing with Cas9 and a pair of sgRNAs to investigate topological folding of chromatin TADs and developmental gene regulation.

DNA Methylation in Osteoarthritis

Osteoarthritis (OA) is a prevalent disease of articular joints and primarily characterized by degradation and calcification of articular cartilage. Presently, no effective treatment other than pain relief exists and patients ultimately need to undergo replacement surgery of the affected joint. During disease progression articular chondrocytes, the single cell type present in articular cartilage, show altered transcriptional profiles and undergo phenotypic changes that resemble the terminal differentiation route apparent in growth plate chondrocytes. Hence, given its prominent function in both regulating gene expression and maintaining cellular phenotypes, DNA methylation of CpG dinucleotides is intensively studied in the context of OA. An increasing number of studies have been published that employed a targeted approach on genes known to play a role in OA pathophysiology. As of such, it has become clear that OA responsive DNA methylation changes seem to mediate disease associated aberrant gene expression. Furthermore, established OA susceptibility alleles such as GDF5 and DIO2 appear to confer OA risk via DNA methylation and respective pathophysiological expression changes. In more recent years, genome wide profiling of DNA methylation in OA affected articular cartilage has emerged as a powerful tool to address the epigenetic changes in their entirety, which has resulted in the identification of putative patient subgroups as well as generic OA associated pathways.

CTCF-mediated chromatin loops enclose inducible gene regulatory domains

Background

The CCTC-binding factor (CTCF) protein is involved in genome organization, including mediating three-dimensional chromatin interactions. Human patient lymphocytes with mutations in a single copy of the CTCF gene have reduced expression of enhancer-associated genes involved in response to stimuli. We hypothesize that CTCF interactions stabilize enhancer-promoter chromatin interaction domains, facilitating increased expression of genes in response to stimuli. Here we systematically investigate this model using computational analyses.

Results

We use CTCF ChIA-PET data from the ENCODE project to show that CTCF-associated chromatin loops have a tendency to enclose regions of enhancer-regulated stimulus responsive genes, insulating them from neighboring regions of constitutively expressed housekeeping genes. To facilitate cell type-specific CTCF loop identification, we develop an algorithm to predict CTCF loops from ChIP-seq data alone by exploiting the CTCF motif directionality in loop anchors. We apply this algorithm to a hundred ENCODE cell line datasets, confirming the universality of our observations as well as identifying a general distinction between primary and immortal cells in loop-enclosed gene content. Finally, we combine the existing evidence to propose a model for the formation of CTCF loops in which partner sites are brought together by chromatin template reeling through stationary RNA polymerases, consistent with the transcription factory hypothesis.

Conclusions

We provide computational evidence that CTCF-mediated chromatin interactions enclose domains of stimulus responsive enhancer-regulated genes, insulating them from nearby housekeeping genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2516-6) contains supplementary material, which is available to authorized users.

Epigenetic silencing of miR-181c by DNA methylation in glioblastoma cell lines

Background

Post-transcriptional regulation by microRNAs is recognized as one of the major pathways for the control of cellular homeostasis. Less well understood is the transcriptional and epigenetic regulation of genes encoding microRNAs. In the present study we addressed the epigenetic regulation of the miR-181c in normal and malignant brain cells.

Methods

To explore the epigenetic regulation of the miR-181c we evaluated its expression using RT-qPCR and the in vivo binding of the CCCTC-binding factor (CTCF) to its regulatory region in different glioblastoma cell lines. DNA methylation survey, chromatin immunoprecipitation and RNA interference assays were used to assess the role of CTCF in the miR-181c epigenetic silencing.

Results

We found that miR-181c is downregulated in glioblastoma cell lines, as compared to normal brain tissues. Loss of expression correlated with a notorious gain of DNA methylation at the miR-181c promoter region and the dissociation of the multifunctional nuclear factor CTCF. Taking advantage of the genomic distribution of CTCF in different cell types we propose that CTCF has a local and cell type specific regulatory role over the miR-181c and not an architectural one through chromatin loop formation. This is supported by the depletion of CTCF in glioblastoma cells affecting the expression levels of NOTCH2 as a target of miR-181c.

Conclusion

Together, our results point to the epigenetic role of CTCF in the regulation of microRNAs implicated in tumorigenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2273-6) contains supplementary material, which is available to authorized users.

Recent advances in understanding nuclear size and shape

Size and shape are important aspects of nuclear structure. While normal cells maintain nuclear size within a defined range, altered nuclear size and shape are associated with a variety of diseases. It is unknown if altered nuclear morphology contributes to pathology, and answering this question requires a better understanding of the mechanisms that control nuclear size and shape. In this review, we discuss recent advances in our understanding of the mechanisms that regulate nuclear morphology, focusing on nucleocytoplasmic transport, nuclear lamins, the endoplasmic reticulum, the cell cycle, and potential links between nuclear size and size regulation of other organelles. We then discuss the functional significance of nuclear morphology in the context of early embryonic development. Looking toward the future, we review new experimental approaches that promise to provide new insights into mechanisms of nuclear size control, in particular microfluidic-based technologies, and discuss how altered nuclear morphology might impact chromatin organization and physiology of diseased cells.

DNA methylation profiling in human Huntington's disease brain

Despite extensive progress in Huntington's disease (HD) research, very little is known about the association of epigenetic variation and HD pathogenesis in human brain tissues. Moreover, its contribution to the tissue-specific transcriptional regulation of the huntingtin gene (HTT), in which HTT expression levels are highest in brain and testes, is currently unknown. To investigate the role of DNA methylation in HD pathogenesis and tissue-specific expression of HTT, we utilized the Illumina HumanMethylation450K BeadChip array to measure DNA methylation in a cohort of age-matched HD and control human cortex and liver tissues. In cortex samples, we found minimal evidence of HD-associated DNA methylation at probed sites after correction for cell heterogeneity but did observe an association with the age of disease onset. In contrast, comparison of matched cortex and liver samples revealed tissue-specific DNA methylation of the HTT gene region at 38 sites (FDR < 0.05). Importantly, we identified a novel differentially methylated binding site in the HTT proximal promoter for the transcription factor CTCF. This CTCF site displayed increased occupancy in cortex, where HTT expression is higher, compared with the liver. Additionally, CTCF silencing reduced the activity of an HTT promoter-reporter construct, suggesting that CTCF plays a role in regulating HTT promoter function. Overall, although we were unable to detect HD-associated DNA methylation alterations at queried sites, we found that DNA methylation may be correlated to the age of disease onset in cortex tissues. Moreover, our data suggest that DNA methylation may, in part, contribute to tissue-specific HTT transcription through differential CTCF occupancy.

Methylation profile of bovine Oct4 gene coding region in relation to three germ layers

Previous studies have shown that octamer-binding transcription factor 4 (Oct4) plays a significant role in early embryonic development of mammalian animals, and different Oct4 expression levels induce multi-lineage differentiation which are regulated by DNA methylation. To explore the relationship between the methylation pattern of Oct4 gene exon 1 and embryonic development, in this work, five different tissues (heart, liver, lung, cerebrum and cerebellum) from three germ layers were chosen from low age (50-60 d) and advanced age (60-70 d) of fetal cattle and the differences between tissues or ages were analyzed, respectively. The result showed that the DNA methylation level of Oct4 gene exon 1 was significant different (P<0.01) between any two of three germ layers in low age (<60 d), but kept steady of advanced age (P>0.05) (>60 d), suggesting that 60-d post coital was an important boundary for embryonic development. In addition, in ectoderm (cerebrum and cerebellum), there was no significant methylation difference of Oct4 gene exon 1 between low age and advanced age (P>0.05), but the result of endoderm (liver and lung) and mesoderm (heart) were on the contrary (P<0.01), which indicated the development of ectoderm was earlier than endoderm and mesoderm. The methylation differences from the 3rd, 5th and 9th CpG-dinucleotide loci of Oct4 gene exon 1 were significantly different between each two of three germ layers (P<0.05), indicating that these three loci may have important influence on bovine embryonic development. This study showed that bovine germ layers differentiation was significantly related to the DNA methylation status of Oct4 gene exon 1. This work firstly identified the DNA methylation profile of bovine Oct4 gene exon 1 and its association with germ layers development in fetus and adult of cattle. Moreover, the work also provided epigenetic information for further studying bovine embryonic development and cellular reprogramming.

3D Chromosome Regulatory Landscape of Human Pluripotent Cells

In this study, we describe the 3D chromosome regulatory landscape of human naive and primed embryonic stem cells. To devise this map, we identified transcriptional enhancers and insulators in these cells and placed them within the context of cohesin-associated CTCF-CTCF loops using cohesin ChIA-PET data. The CTCF-CTCF loops we identified form a chromosomal framework of insulated neighborhoods, which in turn form topologically associating domains (TADs) that are largely preserved during the transition between the naive and primed states. Regulatory changes in enhancer-promoter interactions occur within insulated neighborhoods during cell state transition. The CTCF anchor regions we identified are conserved across species, influence gene expression, and are a frequent site of mutations in cancer cells, underscoring their functional importance in cellular regulation. These 3D regulatory maps of human pluripotent cells therefore provide a foundation for future interrogation of the relationships between chromosome structure and gene control in development and disease.

The three-dimensional cancer genome

The past decade of cancer research has ushered in a comprehensive understanding of the way that the sequence of the genome can be co-opted during the process of tumorigenesis. However, only recently has the epigenome, and in particular the three-dimensional topology of chromatin, been implicated in cancer progression. Here we review recent findings of how the cancer genome is regulated and dysregulated to effect changes in 3D genome topology. We discuss the impact of the spatial organization of the genome on the frequency of tumorigenic chromosomal translocations and the effects of disruption of the proteins responsible for the establishment of chromatin loops. Alteration of the three-dimensional cancer genome is a rapidly emerging hallmark of multiple cancer subtypes.

Insulator dysfunction and oncogene activation in IDH mutant gliomas

Gain-of-function IDH mutations are initiating events that define major clinical and prognostic classes of gliomas. Mutant IDH protein produces a new onco-metabolite, 2-hydroxyglutarate, which interferes with iron-dependent hydroxylases, including the TET family of 5'-methylcytosine hydroxylases. TET enzymes catalyse a key step in the removal of DNA methylation. IDH mutant gliomas thus manifest a CpG island methylator phenotype (G-CIMP), although the functional importance of this altered epigenetic state remains unclear. Here we show that human IDH mutant gliomas exhibit hypermethylation at cohesin and CCCTC-binding factor (CTCF)-binding sites, compromising binding of this methylation-sensitive insulator protein. Reduced CTCF binding is associated with loss of insulation between topological domains and aberrant gene activation. We specifically demonstrate that loss of CTCF at a domain boundary permits a constitutive enhancer to interact aberrantly with the receptor tyrosine kinase gene PDGFRA, a prominent glioma oncogene. Treatment of IDH mutant gliomaspheres with a demethylating agent partially restores insulator function and downregulates PDGFRA. Conversely, CRISPR-mediated disruption of the CTCF motif in IDH wild-type gliomaspheres upregulates PDGFRA and increases proliferation. Our study suggests that IDH mutations promote gliomagenesis by disrupting chromosomal topology and allowing aberrant regulatory interactions that induce oncogene expression.

Stress and glucocorticoid receptor transcriptional programming in time and space: Implications for the brain-gut axis

BACKGROUND

Chronic psychological stress is associated with enhanced abdominal pain and altered intestinal barrier function that may result from a perturbation in the hypothalamic-pituitary-adrenal (HPA) axis. The glucocorticoid receptor (GR) exploits diverse mechanisms to activate or suppress congeneric gene expression, with regulatory variation associated with stress-related disorders in psychiatry and gastroenterology.

PURPOSE

During acute and chronic stress, corticotropin-releasing hormone drives secretion of adrenocorticotropic hormone from the pituitary, ultimately leading to the release of cortisol (human) and corticosterone (rodent) from the adrenal glands. Cortisol binds with the GR in the cytosol, translocates to the nucleus, and activates the NR3C1 (nuclear receptor subfamily 3, group C, member 1 [GR]) gene. This review focuses on the rapidly developing observations that cortisol is responsible for driving circadian and ultradian bursts of transcriptional activity in the CLOCK (clock circadian regulator) and PER (period circadian clock 1) gene families, and this rhythm is disrupted in major depressive disorder, bipolar disorder, and stress-related gastrointestinal and immune disorders. Glucocorticoid receptor regulates different sets of transcripts in a tissue-specific manner, through pulsatile waves of gene expression that includes occupancy of glucocorticoid response elements located within constitutively open spatial domains in chromatin. Emerging evidence supports a potentially pivotal role for epigenetic regulation of how GR interacts with other chromatin regulators to control the expression of its target genes. Dysregulation of the central and peripheral GR regulome has potentially significant consequences for stress-related disorders affecting the brain-gut axis.

Methylation QTLs in the developing brain and their enrichment in schizophrenia risk loci

We characterized DNA methylation quantitative trait loci (mQTLs) in a large collection (n = 166) of human fetal brain samples spanning 56-166 d post-conception, identifying >16,000 fetal brain mQTLs. Fetal brain mQTLs were primarily cis-acting, enriched in regulatory chromatin domains and transcription factor binding sites, and showed substantial overlap with genetic variants that were also associated with gene expression in the brain. Using tissue from three distinct regions of the adult brain (prefrontal cortex, striatum and cerebellum), we found that most fetal brain mQTLs were developmentally stable, although a subset was characterized by fetal-specific effects. Fetal brain mQTLs were enriched amongst risk loci identified in a recent large-scale genome-wide association study (GWAS) of schizophrenia, a severe psychiatric disorder with a hypothesized neurodevelopmental component. Finally, we found that mQTLs can be used to refine GWAS loci through the identification of discrete sites of variable fetal brain methylation associated with schizophrenia risk variants.

Towards a predictive model of chromatin 3D organization

Architectural proteins mediate interactions between distant regions in the genome to bring together different regulatory elements while establishing a specific three-dimensional organization of the genetic material. Depletion of specific architectural proteins leads to miss regulation of gene expression and alterations in nuclear organization. The specificity of interactions mediated by architectural proteins depends on the nature, number, and orientation of their binding site at individual genomic locations. Knowledge of the mechanisms and rules governing interactions among architectural proteins may provide a code to predict the 3D organization of the genome.

Systematic identification and annotation of human methylation marks based on bisulfite sequencing methylomes reveals distinct roles of cell type-specific hypomethylation in the regulation of cell identity genes

DNA methylation is a key epigenetic mark that is critical for gene regulation in multicellular eukaryotes. Although various human cell types may have the same genome, these cells have different methylomes. The systematic identification and characterization of methylation marks across cell types are crucial to understand the complex regulatory network for cell fate determination. In this study, we proposed an entropy-based framework termed SMART to integrate the whole genome bisulfite sequencing methylomes across 42 human tissues/cells and identified 757 887 genome segments. Nearly 75% of the segments showed uniform methylation across all cell types. From the remaining 25% of the segments, we identified cell type-specific hypo/hypermethylation marks that were specifically hypo/hypermethylated in a minority of cell types using a statistical approach and presented an atlas of the human methylation marks. Further analysis revealed that the cell type-specific hypomethylation marks were enriched through H3K27ac and transcription factor binding sites in cell type-specific manner. In particular, we observed that the cell type-specific hypomethylation marks are associated with the cell type-specific super-enhancers that drive the expression of cell identity genes. This framework provides a complementary, functional annotation of the human genome and helps to elucidate the critical features and functions of cell type-specific hypomethylation.

Large-scale identification of sequence variants influencing human transcription factor occupancy in vivo

The function of human regulatory regions depends exquisitely on their local genomic environment and on cellular context, complicating experimental analysis of common disease- and trait-associated variants that localize within regulatory DNA. We use allelically resolved genomic DNase I footprinting data encompassing 166 individuals and 114 cell types to identify >60,000 common variants that directly influence transcription factor occupancy and regulatory DNA accessibility in vivo. The unprecedented scale of these data enables systematic analysis of the impact of sequence variation on transcription factor occupancy in vivo. We leverage this analysis to develop accurate models of variation affecting the recognition sites for diverse transcription factors and apply these models to discriminate nearly 500,000 common regulatory variants likely to affect transcription factor occupancy across the human genome. The approach and results provide a new foundation for the analysis and interpretation of noncoding variation in complete human genomes and for systems-level investigation of disease-associated variants.

Genome-wide DNA methylation map of human neutrophils reveals widespread inter-individual epigenetic variation

The extent of variation in DNA methylation patterns in healthy individuals is not yet well documented. Identification of inter-individual epigenetic variation is important for understanding phenotypic variation and disease susceptibility. Using neutrophils from a cohort of healthy individuals, we generated base-resolution DNA methylation maps to document inter-individual epigenetic variation. We identified 12851 autosomal inter-individual variably methylated fragments (iVMFs). Gene promoters were the least variable, whereas gene body and upstream regions showed higher variation in DNA methylation. The iVMFs were relatively enriched in repetitive elements compared to non-iVMFs, and were associated with genome regulation and chromatin function elements. Further, variably methylated genes were disproportionately associated with regulation of transcription, responsive function and signal transduction pathways. Transcriptome analysis indicates that iVMF methylation at differentially expressed exons has a positive correlation and local effect on the inclusion of that exon in the mRNA transcript.

Trans effects of chromosome aneuploidies on DNA methylation patterns in human Down syndrome and mouse models

Background

Trisomy 21 causes Down syndrome (DS), but the mechanisms by which the extra chromosome leads to deficient intellectual and immune function are not well understood.

Results

Here, we profile CpG methylation in DS and control cerebral and cerebellar cortex of adults and cerebrum of fetuses. We purify neuronal and non-neuronal nuclei and T lymphocytes and find biologically relevant genes with DS-specific methylation (DS-DM) in each of these cell types. Some genes show brain-specific DS-DM, while others show stronger DS-DM in T cells. Both 5-methyl-cytosine and 5-hydroxy-methyl-cytosine contribute to the DS-DM. Thirty percent of genes with DS-DM in adult brain cells also show DS-DM in fetal brains, indicating early onset of these epigenetic changes, and we find early maturation of methylation patterns in DS brain and lymphocytes. Some, but not all, of the DS-DM genes show differential expression. DS-DM preferentially affected CpGs in or near specific transcription factor binding sites (TFBSs), implicating a mechanism involving altered TFBS occupancy. Methyl-seq of brain DNA from mouse models with sub-chromosomal duplications mimicking DS reveals partial but significant overlaps with human DS-DM and shows that multiple chromosome 21 genes contribute to the downstream epigenetic effects.

Conclusions

These data point to novel biological mechanisms in DS and have general implications for trans effects of chromosomal duplications and aneuploidies on epigenetic patterning.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0827-6) contains supplementary material, which is available to authorized users.

A promising hypothesis of c-KIT methylation/ expression paradox in c-KIT (+) squamous cell carcinoma of uterine cervix ----- CTCF transcriptional repressor regulates c-KIT proto-oncogene expression

We recently reported one interesting case showing mutation-free c-KIT proto-oncogene overexpression and paradoxical hypermethylation in 54 cases of primary squamous cell carcinoma of uterine cervix (SCC). However, its molecular mechanisms still remain unknown. We propose the hypothesis that increased methylation at the CpG islands on the promoter near the first exon region might interfere with the binding of CTCF repressor with c-KIT promoter that regulates c-KIT proto-oncogene expression in such case. Further studies focusing on the status of epigenetic modifications of mutation-free c-KIT (+) tumors are encouraged.

The ubiquitous transcription factor CTCF promotes lineage-specific epigenomic remodeling and establishment of transcriptional networks driving cell differentiation

Cell differentiation relies on tissue-specific transcription factors (TFs) that cooperate to establish unique transcriptomes and phenotypes. However, the role of ubiquitous TFs in these processes remains poorly defined. Recently, we have shown that the CCCTC-binding factor (CTCF) is required for adipocyte differentiation through epigenomic remodelling of adipose tissue-specific enhancers and transcriptional activation of Peroxisome proliferator-activated receptor gamma (PPARG), the main driver of the adipogenic program (PPARG), and its target genes. Here, we discuss how these findings, together with the recent literature, illuminate a functional role for ubiquitous TFs in lineage-determining transcriptional networks.

CTCF negatively regulates HOXA10 expression in breast cancer cells

HOX genes not only play important roles in defining body patterning during embryonic development, but also control numerous cellular events in adult cells. Deregulated HOX gene expression in different cancers including breast cancer is now increasingly being reported. Given that human HOXA cluster is marked with several CTCF binding sites, we investigated whether the presence of CTCF is associated directly with expression of HOXA genes in breast cancer cells. Several HOX genes, such as HOXA4, HOXA5 and HOXA10, were deregulated by CTCF overexpression and knockdown in MCF-7 cells. Among these genes, HOXA10 is an emerging tumor suppressor for its role in activation of p53 and in countering tumorigenesis in breast cancer. Here we provided evidences that CTCF functions as a negative regulator of HOXA10 in breast cancer cells. The putative promoter region of HOXA10 lies between 5.3 and 6.1 kb upstream of its start codon and its promoter activity was negatively regulated by CTCF. Together with in-silico analysis and in vitro mutation assay we identified a 20 bp CTCF binding motif flanking with core promoter element of HOXA10. HOXA10 promoter region was kept inactivated by maintaining H3K27me3 inactivation marks in the presence of CTCF. Epigenetic silencing of HOXA10 by CTCF in breast cancer cells may contribute towards tumorigenesis by decreasing apoptosis and promoting metastasis.

Topologically Associating Domains: An invariant framework or a dynamic scaffold?

Metazoan genomes are organized into regions of topologically associating domains (TADs). TADs are demarcated by border elements, which are enriched for active genes and high occupancy architectural protein binding sites. We recently demonstrated that 3D chromatin architecture is dynamic in response to heat shock, a physiological stress that downregulates transcription and causes a global redistribution of architectural proteins. We utilized a quantitative measure of border strength after heat shock, transcriptional inhibition, and architectural protein knockdown to demonstrate that changes in both transcription and architectural protein occupancy contribute to heat shock-induced TAD dynamics. Notably, architectural proteins appear to play a more important role in altering 3D chromatin architecture. Here, we discuss the implications of our findings on previous studies evaluating the dynamics of TAD structure during cellular differentiation. We propose that the subset of variable TADs observed after differentiation are representative of cell-type specific gene expression and are biologically significant.

Functional diversity of CTCFs is encoded in their binding motifs

BACKGROUND

The CCCTC-binding factor (CTCF) has diverse regulatory functions. However, the definitive characteristics of the CTCF binding motif required for its functional diversity still remains elusive.

RESULTS

Here, we describe a new motif discovery workflow by which we have identified three CTCF binding motif variations with highly divergent functionalities. Supported by transcriptomic, epigenomic and chromatin-interactomic data, we show that the functional diversity of the CTCF binding motifs is strongly associated with their GC content, CpG dinucleotide coverage and relative DNA methylation level at the 12th position of the motifs. Further analysis suggested that the co-localization of cohesin, the key factor in cohesion of sister chromatids, is negatively correlated with the CpG coverage and the relative DNA methylation level at the 12th position. Finally, we present evidences for a hypothetical model in which chromatin interactions between promoters and distal regulatory regions are likely mediated by CTCFs binding to sequences with high CpG.

CONCLUSION

These results demonstrate the existence of definitive CTCF binding motifs corresponding to CTCF's diverse functions, and that the functional diversity of the motifs is strongly associated with genetic and epigenetic features at the 12th position of the motifs.

DNA-Methylation Patterns in Trisomy 21 Using Cells from Monozygotic Twins

DNA methylation is essential in mammalian development. We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in Down syndrome (DS). In order to determine the methylation differences in T21 without interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. Reduced Representation Bisulfite Sequencing (RRBS) revealed 35 differentially methylated promoter regions (DMRs) (Absolute methylation differences = 25%, FDR < 0.001) in MZ twins discordant for T21 that have also been observed in comparison between unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis (FDR = 1.60 e -03) and include genes of the HOXB and HOXD clusters. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) observed in the iPSC of the T21 versus normal twin. Altogether, the results of this study highlight the epigenetic effects of the extra chromosome 21 in T21 on loci outside of this chromosome that are relevant to DS associated phenotypes.

Argonaute 2 Binds Directly to tRNA Genes and Promotes Gene Repression in cis

To further our understanding of the RNAi machinery within the human nucleus, we analyzed the chromatin and RNA binding of Argonaute 2 (AGO2) within human cancer cell lines. Our data indicated that AGO2 binds directly to nascent tRNA and 5S rRNA, and to the genomic loci from which these RNAs are transcribed, in a small RNA- and DICER-independent manner. AGO2 chromatin binding was not observed at non-TFIIIC-dependent RNA polymerase III (Pol III) genes or at extra-TFIIIC (ETC) sites, indicating that the interaction is specific for TFIIIC-dependent Pol III genes. A genome-wide analysis indicated that loss of AGO2 caused a global increase in mRNA expression level among genes that flank AGO2-bound tRNA genes. This effect was shown to be distinct from that of the disruption of DICER, DROSHA, or CTCF. We propose that AGO2 binding to tRNA genes has a novel and important regulatory role in human cells.

Comparative analyses of CTCF and BORIS occupancies uncover two distinct classes of CTCF binding genomic regions

Background

CTCF and BORIS (CTCFL), two paralogous mammalian proteins sharing nearly identical DNA binding domains, are thought to function in a mutually exclusive manner in DNA binding and transcriptional regulation.

Results

Here we show that these two proteins co-occupy a specific subset of regulatory elements consisting of clustered CTCF binding motifs (termed 2xCTSes). BORIS occupancy at 2xCTSes is largely invariant in BORIS-positive cancer cells, with the genomic pattern recapitulating the germline-specific BORIS binding to chromatin. In contrast to the single-motif CTCF target sites (1xCTSes), the 2xCTS elements are preferentially found at active promoters and enhancers, both in cancer and germ cells. 2xCTSes are also enriched in genomic regions that escape histone to protamine replacement in human and mouse sperm. Depletion of the BORIS gene leads to altered transcription of a large number of genes and the differentiation of K562 cells, while the ectopic expression of this CTCF paralog leads to specific changes in transcription in MCF7 cells.

Conclusions

We discover two functionally and structurally different classes of CTCF binding regions, 2xCTSes and 1xCTSes, revealed by their predisposition to bind BORIS. We propose that 2xCTSes play key roles in the transcriptional program of cancer and germ cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0736-8) contains supplementary material, which is available to authorized users.

CRISPR Inversion of CTCF Sites Alters Genome Topology and Enhancer/Promoter Function

CTCF and the associated cohesin complex play a central role in insulator function and higher-order chromatin organization of mammalian genomes. Recent studies identified a correlation between the orientation of CTCF-binding sites (CBSs) and chromatin loops. To test the functional significance of this observation, we combined CRISPR/Cas9-based genomic-DNA-fragment editing with chromosome-conformation-capture experiments to show that the location and relative orientations of CBSs determine the specificity of long-range chromatin looping in mammalian genomes, using protocadherin (Pcdh) and β-globin as model genes. Inversion of CBS elements within the Pcdh enhancer reconfigures the topology of chromatin loops between the distal enhancer and target promoters and alters gene-expression patterns. Thus, although enhancers can function in an orientation-independent manner in reporter assays, in the native chromosome context, the orientation of at least some enhancers carrying CBSs can determine both the architecture of topological chromatin domains and enhancer/promoter specificity. These findings reveal how 3D chromosome architecture can be encoded by linear genome sequences.

Role of DNA Methylation in Modulating Transcription Factor Occupancy

Although DNA methylation is commonly invoked as a mechanism for transcriptional repression, the extent to which it actively silences transcription factor (TF) occupancy sites in vivo is unknown. To study the role of DNA methylation in the active modulation of TF binding, we quantified the effect of DNA methylation depletion on the genomic occupancy patterns of CTCF, an abundant TF with known methylation sensitivity that is capable of autonomous binding to its target sites in chromatin. Here, we show that the vast majority (>98.5%) of the tens of thousands of unoccupied, methylated CTCF recognition sequences remain unbound upon abrogation of DNA methylation. The small fraction of sites that show methylation-dependent binding in vivo are in turn characterized by highly variable CTCF occupancy across cell types. Our results suggest that DNA methylation is not a primary groundskeeper of genomic TF landscapes, but rather a specialized mechanism for stabilizing intrinsically labile sites.

De novo DNA Methyltransferases Dnmt3a and Dnmt3b regulate the onset of Igκ light chain rearrangement during early B-cell development

Immunoglobulin genes V(D)J rearrangement during early lymphopoiesis is a critical process involving sequential recombination of the heavy and light chain loci. A number of transcription factors act together with temporally activated recombinases and chromatin accessibility changes to regulate this complex process. Here, we deleted the de novo DNA methyltransferases Dnmt3a and Dnmt3b in early B cells of conditionally targeted mice, and monitored the process of V(D)J recombination. Dnmt3a and Dnmt3b deletion resulted in precocious recombination of the immunoglobulin κ light chain without impairing the differentiation of mature B cells or overall B-cell development. Ex vivo culture of IL-7 restricted early B-cell progenitors lacking Dnmt3a and Dnmt3b showed precocious Vκ-Jκ rearrangements that are limited to the proximal Vκ genes. Furthermore, B-cell progenitors deficient in Dnmt3a and Dnmt3b showed elevated levels of germline transcripts at the proximal Vκ genes, alterations in methylation patterns at Igκ enhancer sites and increased expression of the transcription factor E2A. Our data suggest that Dnmt3a and Dnmt3b are critical to regulate the onset of Igκ light chain rearrangement during early B-cell development.

Targeting CTCF to Control Virus Gene Expression: A Common Theme amongst Diverse DNA Viruses

All viruses target host cell factors for successful life cycle completion. Transcriptional control of DNA viruses by host cell factors is important in the temporal and spatial regulation of virus gene expression. Many of these factors are recruited to enhance virus gene expression and thereby increase virus production, but host cell factors can also restrict virus gene expression and productivity of infection. CCCTC binding factor (CTCF) is a host cell DNA binding protein important for the regulation of genomic chromatin boundaries, transcriptional control and enhancer element usage. CTCF also functions in RNA polymerase II regulation and in doing so can influence co-transcriptional splicing events. Several DNA viruses, including Kaposi's sarcoma-associated herpesvirus (KSHV), Epstein-Barr virus (EBV) and human papillomavirus (HPV) utilize CTCF to control virus gene expression and many studies have highlighted a role for CTCF in the persistence of these diverse oncogenic viruses. CTCF can both enhance and repress virus gene expression and in some cases CTCF increases the complexity of alternatively spliced transcripts. This review article will discuss the function of CTCF in the life cycle of DNA viruses in the context of known host cell CTCF functions.

Quantification of Retinogenesis in 3D Cultures Reveals Epigenetic Memory and Higher Efficiency in iPSCs Derived from Rod Photoreceptors

Cell-based therapies to treat retinal degeneration are now being tested in clinical trials. However, it is not known whether the source of stem cells is important for the production of differentiated cells suitable for transplantation. To test this, we generated induced pluripotent stem cells (iPSCs) from murine rod photoreceptors (r-iPSCs) and scored their ability to make retinae by using a standardized quantitative protocol called STEM-RET. We discovered that r-iPSCs more efficiently produced differentiated retinae than did embryonic stem cells (ESCs) or fibroblast-derived iPSCs (f-iPSCs). Retinae derived from f-iPSCs had fewer amacrine cells and other inner nuclear layer cells. Integrated epigenetic analysis showed that DNA methylation contributes to the defects in f-iPSC retinogenesis and that rod-specific CTCF insulator protein-binding sites may promote r-iPSC retinogenesis. Together, our data suggest that the source of stem cells is important for producing retinal neurons in three-dimensional (3D) organ cultures.

CTCF as a multifunctional protein in genome regulation and gene expression

CCCTC-binding factor (CTCF) is a highly conserved zinc finger protein and is best known as a transcription factor. It can function as a transcriptional activator, a repressor or an insulator protein, blocking the communication between enhancers and promoters. CTCF can also recruit other transcription factors while bound to chromatin domain boundaries. The three-dimensional organization of the eukaryotic genome dictates its function, and CTCF serves as one of the core architectural proteins that help establish this organization. The mapping of CTCF-binding sites in diverse species has revealed that the genome is covered with CTCF-binding sites. Here we briefly describe the diverse roles of CTCF that contribute to genome organization and gene expression.

Circulating DNA of HOTAIR in serum is a novel biomarker for breast cancer

Long non-coding HOX transcript antisense intergenic RNA (HOTAIR) plays an important role in breast cancer. The purpose of this study was to determine whether circulating HOTAIR can be used for breast cancer diagnosis. HOTAIR in serum was measured by PCR-based direct detection. Reverse transcriptase and DNase I treatment were used to distinguish the DNA and RNA forms of HOTAIR. To determine whether circulating HOTAIR is a biomarker for breast cancer, the DNA of HOTAIR from breast cancer patients and healthy controls was measured at both the discovery stage (48 individuals) and an independent validation stage (156 individuals). The diagnostic accuracy was assessed by the receiver operating characteristic curve (ROC) and the area under the curve (AUC). We showed that the major form of HOTAIR-derived fragment in serum is DNA rather than RNA in our study, the same as for MALAT-1, another well-described lincRNA. A higher circulating DNA level of HOTAIR was found in patients at the discovery stage (P = 0.0008). ROC analysis revealed that the circulating HOTAIR DNA distinguished breast cancer patients from healthy individuals (AUC = 0.799). This finding was confirmed at the validation stage. Though circulating MALAT-1 DNA was altered in the discovery stage, it showed no significant difference in the validation stage. In the entire set of 204 samples, the circulating HOTAIR DNA showed a 2.15-fold change in patients compared with healthy controls (P < 0.0001, AUC = 0.786). The optimal cutoff value for diagnosis was 0.30 with sensitivity of 80.0 % and specificity of 68.3 %. Moreover, a correlation between the DNA level of circulating HOTAIR and the progress of breast cancer was established. We have demonstrated that the circulating DNA of HOTAIR is a potential biomarker for breast cancer.

Architectural proteins, transcription, and the three-dimensional organization of the genome

Architectural proteins mediate interactions between distant sequences in the genome. Two well-characterized functions of architectural protein interactions include the tethering of enhancers to promoters and bringing together Polycomb-containing sites to facilitate silencing. The nature of which sequences interact genome-wide appears to be determined by the orientation of the architectural protein binding sites as well as the number and identity of architectural proteins present. Ultimately, long range chromatin interactions result in the formation of loops within the chromatin fiber. In this review, we discuss data suggesting that architectural proteins mediate long range chromatin interactions that both facilitate and hinder neighboring interactions, compartmentalizing the genome into regions of highly interacting chromatin domains.

A chromatin code for alternative splicing involving a putative association between CTCF and HP1α proteins

Background

Alternative splicing is primarily controlled by the activity of splicing factors and by the elongation of the RNA polymerase II (RNAPII). Recent experiments have suggested a new complex network of splicing regulation involving chromatin, transcription and multiple protein factors. In particular, the CCCTC-binding factor (CTCF), the Argonaute protein AGO1, and members of the heterochromatin protein 1 (HP1) family have been implicated in the regulation of splicing associated with chromatin and the elongation of RNAPII. These results raise the question of whether these proteins may associate at the chromatin level to modulate alternative splicing.

Results

Using chromatin immunoprecipitation sequencing (ChIP-Seq) data for CTCF, AGO1, HP1α, H3K27me3, H3K9me2, H3K36me3, RNAPII, total H3 and 5metC and alternative splicing arrays from two cell lines, we have analyzed the combinatorial code of their binding to chromatin in relation to the alternative splicing patterns between two cell lines, MCF7 and MCF10. Using Machine Learning techniques, we identified the changes in chromatin signals that are most significantly associated with splicing regulation between these two cell lines. Moreover, we have built a map of the chromatin signals on the pre-mRNA, that is, a chromatin-based RNA-map, which can explain 606 (68.55%) of the regulated events between MCF7 and MCF10. This chromatin code involves the presence of HP1α, CTCF, AGO1, RNAPII and histone marks around regulated exons and can differentiate patterns of skipping and inclusion. Additionally, we found a significant association of HP1α and CTCF activities around the regulated exons and a putative DNA binding site for HP1α.

Conclusions

Our results show that a considerable number of alternative splicing events could have a chromatin-dependent regulation involving the association of HP1α and CTCF near regulated exons. Additionally, we find further evidence for the involvement of HP1α and AGO1 in chromatin-related splicing regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0141-5) contains supplementary material, which is available to authorized users.

Vanin-1 expression and methylation discriminate pediatric asthma corticosteroid treatment response

BACKGROUND

There is considerable heterogeneity in asthma treatment response.

OBJECTIVE

We sought to identify biomarkers of corticosteroid treatment response in children with asthma and evaluate the utility and mechanistic basis of these biomarkers.

METHODS

Children (5-18 years) presenting to the emergency department with an acute asthma exacerbation were recruited and followed during hospitalization. Nasal epithelial cells were collected on presentation to the emergency department (T0) and 18 to 24 hours later (T1), and T1/T0 gene expression ratios were analyzed to identify genes associated with good and poor corticosteroid treatment response phenotypes. The utility of these genes in discriminating between systemic corticosteroid treatment response groups was then tested prospectively in a new cohort of patients. A gene candidate (vanin-1 [VNN1]) that consistently distinguished good versus poor response phenotypes was further studied in an experimental asthma model, and VNN1 promoter methylation was measured by means of bisulfite pyrosequencing in patients.

RESULTS

VNN1 mRNA expression changes were associated with systemic corticosteroid treatment response in children with acute asthma, and VNN1 was required for optimal response to corticosteroid treatment in an experimental asthma model. A CpG site within the VNN1 promoter was differentially methylated between good versus poor treatment response groups, and methylation at this site correlated with VNN1 mRNA expression.

CONCLUSIONS

We have identified a biological basis for poor corticosteroid treatment response that can be used to distinguish a subgroup of asthmatic children who respond poorly to systemic corticosteroid treatment. VNN1 contributes to corticosteroid responsiveness, and changes in VNN1 nasal epithelial mRNA expression and VNN1 promoter methylation might be clinically useful biomarkers of treatment response in asthmatic children.

Spatial enhancer clustering and regulation of enhancer-proximal genes by cohesin

In addition to mediating sister chromatid cohesion during the cell cycle, the cohesin complex associates with CTCF and with active gene regulatory elements to form long-range interactions between its binding sites. Genome-wide chromosome conformation capture had shown that cohesin's main role in interphase genome organization is in mediating interactions within architectural chromosome compartments, rather than specifying compartments per se. However, it remains unclear how cohesin-mediated interactions contribute to the regulation of gene expression. We have found that the binding of CTCF and cohesin is highly enriched at enhancers and in particular at enhancer arrays or "super-enhancers" in mouse thymocytes. Using local and global chromosome conformation capture, we demonstrate that enhancer elements associate not just in linear sequence, but also in 3D, and that spatial enhancer clustering is facilitated by cohesin. The conditional deletion of cohesin from noncycling thymocytes preserved enhancer position, H3K27ac, H4K4me1, and enhancer transcription, but weakened interactions between enhancers. Interestingly, ∼ 50% of deregulated genes reside in the vicinity of enhancer elements, suggesting that cohesin regulates gene expression through spatial clustering of enhancer elements. We propose a model for cohesin-dependent gene regulation in which spatial clustering of enhancer elements acts as a unified mechanism for both enhancer-promoter "connections" and "insulation."

Using epigenomic studies in monozygotic twins to improve our understanding of cancer

Cancer is a set of diseases that exhibit not only genetic mutations but also a profoundly distorted epigenetic landscape. Over the last two decades, great advances have been made in identifying these alterations and their importance in the initiation and progression of cancer. Epigenetic changes can be seen from the very early stages in tumorigenesis and dysregulation of the epigenome has an increasingly acknowledged pathogenic role. Epigenomic twin studies have great potential to contribute to our understanding of complex diseases, such as cancer. This is because the use of monozygotic twins discordant for cancer enables epigenetic variation analysis without the confounding influence of the constitutive genetic background, age or cohort effects. It therefore allows the identification of susceptibility loci that may be sensitive to modification by the environment. These studies into cancer etiology will potentially lead to robust epigenetic markers for the detection and risk assessment of cancer.

Modeling the relationship of epigenetic modifications to transcription factor binding

Transcription factors (TFs) and epigenetic modifications play crucial roles in the regulation of gene expression, and correlations between the two types of factors have been discovered. However, methods for quantitatively studying the correlations remain limited. Here, we present a computational approach to systematically investigating how epigenetic changes in chromatin architectures or DNA sequences relate to TF binding. We implemented statistical analyses to illustrate that epigenetic modifications are predictive of TF binding affinities, without the need of sequence information. Intriguingly, by considering genome locations relative to transcription start sites (TSSs) or enhancer midpoints, our analyses show that different locations display various relationship patterns. For instance, H3K4me3, H3k9ac and H3k27ac contribute more in the regions near TSSs, whereas H3K4me1 and H3k79me2 dominate in the regions far from TSSs. DNA methylation plays relatively important roles when close to TSSs than in other regions. In addition, the results show that epigenetic modification models for the predictions of TF binding affinities are cell line-specific. Taken together, our study elucidates highly coordinated, but location- and cell type-specific relationships between epigenetic modifications and binding affinities of TFs.