Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity

Changes in H3K27ac following lipopolysaccharide stimulation of nasopharyngeal epithelial cells

BACKGROUND

The epithelium is the first line of defense against pathogens. Notably the epithelial cells lining the respiratory track are crucial in sensing airborne microbes and mounting an effective immune response via the expression of target genes such as cytokines and chemokines. Gene expression regulation following microbial recognition is partly regulated by chromatin re-organization and has been described in immune cells but data from epithelial cells is not as detailed. Here, we report genome-wide changes of the H3K27ac mark, characteristic of activated enhancers and promoters, after stimulation of nasopharyngeal epithelial cells with the bacterial endotoxin Lipopolysaccharide (LPS).

RESULTS

In this study, we have identified 626 regions where the H3K27ac mark showed reproducible increase following LPS induction in epithelial cells. This indicated that sensing of LPS led to opening of the chromatin in our system. Moreover, this phenomenon seemed to happen extensively at enhancers regions and we could observe instances of Super-enhancer formation. As expected, LPS-increased H3K27ac regions were found in the vicinity of genes relevant for LPS response and these changes correlated with up-regulation of their expression. In addition, we found the induction of H3K27ac mark to overlap with the binding of one of the NF-kB members and key regulator of the innate immune response, RELA, following LPS sensing. Indeed, inhibiting the NF-kB pathway abolished the deposition of H3K27ac at the TNF locus, a target of RELA, suggesting that these two phenomena are associated.

CONCLUSIONS

Enhancers' selection and activation following microbial or inflammatory stimuli has been described previously and shown to be mediated via the NF-kB pathway. Here, we demonstrate that this is also likely to occur in the case of LPS-sensing by nasopharyngeal epithelial cells as well. In addition to validating previous findings, we generated a valuable data set relevant to the host immune response to epithelial cell colonizing or infecting pathogens.

Gonadal supporting cells acquire sex-specific chromatin landscapes during mammalian sex determination

Cis-regulatory elements are critical for the precise spatiotemporal regulation of genes during development. However, identifying functional regulatory sites that drive cell differentiation in vivo has been complicated by the high numbers of cells required for whole-genome epigenetic assays. Here, we identified putative regulatory elements during sex determination by performing ATAC-seq and ChIP-seq for H3K27ac in purified XX and XY gonadal supporting cells before and after sex determination in mice. We show that XX and XY supporting cells initiate sex determination with similar chromatin landscapes and acquire sex-specific regulatory elements as they commit to the male or female fate. To validate our approach, we identified a functional gonad-specific enhancer downstream of Bmp2, an ovary-promoting gene. This work increases our understanding of the complex regulatory network underlying mammalian sex determination and provides a powerful resource for identifying non-coding regulatory elements that could harbor mutations that lead to Disorders of Sexual Development.

Polycomb repressive complex 1 shapes the nucleosome landscape but not accessibility at target genes

Polycomb group (PcG) proteins are transcriptional repressors that play important roles in regulating gene expression during animal development. In vitro experiments have shown that PcG protein complexes can compact chromatin to limit the activity of chromatin remodeling enzymes and access of the transcriptional machinery to DNA. In fitting with these ideas, gene promoters associated with PcG proteins have been reported to be less accessible than other gene promoters. However, it remains largely untested in vivo whether PcG proteins define chromatin accessibility or other chromatin features. To address this important question, we examine the chromatin accessibility and nucleosome landscape at PcG protein-bound promoters in mouse embryonic stem cells using the assay for transposase accessible chromatin (ATAC)-seq. Combined with genetic ablation strategies, we unexpectedly discover that although PcG protein-occupied gene promoters exhibit reduced accessibility, this does not rely on PcG proteins. Instead, the Polycomb repressive complex 1 (PRC1) appears to play a unique role in driving elevated nucleosome occupancy and decreased nucleosomal spacing in Polycomb chromatin domains. Our new genome-scale observations argue, in contrast to the prevailing view, that PcG proteins do not significantly affect chromatin accessibility and highlight an underappreciated complexity in the relationship between chromatin accessibility, the nucleosome landscape, and PcG-mediated transcriptional repression.

coTRaCTE predicts co-occurring transcription factors within cell-type specific enhancers

Cell-type specific gene expression is regulated by the combinatorial action of transcription factors (TFs). In this study, we predict transcription factor (TF) combinations that cooperatively bind in a cell-type specific manner. We first divide DNase hypersensitive sites into cell-type specifically open vs. ubiquitously open sites in 64 cell types to describe possible cell-type specific enhancers. Based on the pattern contrast between these two groups of sequences we develop "co-occurring TF predictor on Cell-Type specific Enhancers" (coTRaCTE) - a novel statistical method to determine regulatory TF co-occurrences. Contrasting the co-binding of TF pairs between cell-type specific and ubiquitously open chromatin guarantees the high cell-type specificity of the predictions. coTRaCTE predicts more than 2000 co-occurring TF pairs in 64 cell types. The large majority (70%) of these TF pairs is highly cell-type specific and overlaps in TF pair co-occurrence are highly consistent among related cell types. Furthermore, independently validated co-occurring and directly interacting TFs are significantly enriched in our predictions. Focusing on the regulatory network derived from the predicted co-occurring TF pairs in embryonic stem cells (ESCs) we find that it consists of three subnetworks with distinct functions: maintenance of pluripotency governed by OCT4, SOX2 and NANOG, regulation of early development governed by KLF4, STAT3, ZIC3 and ZNF148 and general functions governed by MYC, TCF3 and YY1. In summary, coTRaCTE predicts highly cell-type specific co-occurring TFs which reveal new insights into transcriptional regulatory mechanisms.

Human cardiac cis-regulatory elements, their cognate transcription factors, and regulatory DNA sequence variants

Cis-regulatory elements (CRE), short DNA sequences through which transcription factors (TFs) exert regulatory control on gene expression, are postulated to be the major sites of causal sequence variation underlying the genetics of complex traits and diseases. We present integrative analyses, combining high-throughput genomic and epigenomic data with sequence-based computations, to identify the causal transcriptional components in a given tissue. We use data on adult human hearts to demonstrate that (1) sequence-based predictions detect numerous, active, tissue-specific CREs missed by experimental observations, (2) learned sequence features identify the cognate TFs, (3) CRE variants are specifically associated with cardiac gene expression, and (4) a significant fraction of the heritability of exemplar cardiac traits (QT interval, blood pressure, pulse rate) is attributable to these variants. This general systems approach can thus identify candidate causal variants and the components of gene regulatory networks (GRN) to enable understanding of the mechanisms of complex disorders on a tissue- or cell-type basis.

Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia

Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants.

Apolipoprotein E region molecular signatures of Alzheimer's disease

Although the APOE region is the strongest genetic risk factor for Alzheimer's diseases (ADs), its pathogenic role remains poorly understood. Elucidating genetic predisposition to ADs, a subset of age-related diseases characteristic for postreproductive period, is hampered by the undefined role of evolution in establishing molecular mechanisms of such diseases. This uncertainty is inevitable source of natural-selection-free genetic heterogeneity in predisposition to ADs. We performed first large-scale analysis of linkage disequilibrium (LD) structures characterized by 30 polymorphisms from five genes in the APOE 19q13.3 region (BCAM, NECTIN2, TOMM40, APOE, and APOC1) in 2,673 AD-affected and 16,246 unaffected individuals from five cohorts. Consistent with the undefined role of evolution in age-related diseases, we found that these structures, being highly heterogeneous, are significantly different in subjects with and without ADs. The pattern of the difference represents molecular signature of AD comprised of single nucleotide polymorphisms (SNPs) from all five genes in the APOE region. Significant differences in LD in subjects with and without ADs indicate SNPs from different genes likely involved in AD pathogenesis. Significant and highly heterogeneous molecular signatures of ADs provide unprecedented insight into complex polygenetic predisposition to ADs in the APOE region. These findings are more consistent with a complex haplotype than with a single genetic variant origin of ADs in this region.

Chromatin reprogramming in breast cancer

Reprogramming of the chromatin landscape is a critical component to the transcriptional response in breast cancer. Effects of sex hormones such as estrogens and progesterone have been well described to have a critical impact on breast cancer proliferation. However, the complex network of the chromatin landscape, enhancer regions and mode of function of steroid receptors (SRs) and other transcription factors (TFs), is an intricate web of signaling and functional processes that is still largely misunderstood at the mechanistic level. In this review, we describe what is currently known about the dynamic interplay between TFs with chromatin and the reprogramming of enhancer elements. Emphasis has been placed on characterizing the different modes of action of TFs in regulating enhancer activity, specifically, how different SRs target enhancer regions to reprogram chromatin in breast cancer cells. In addition, we discuss current techniques employed to study enhancer function at a genome-wide level. Further, we have noted recent advances in live cell imaging technology. These single-cell approaches enable the coupling of population-based assays with real-time studies to address many unsolved questions about SRs and chromatin dynamics in breast cancer.

Towards a map of cis-regulatory sequences in the human genome

Accumulating evidence indicates that transcription factor (TF) binding sites, or cis-regulatory elements (CREs), and their clusters termed cis-regulatory modules (CRMs) play a more important role than do gene-coding sequences in specifying complex traits in humans, including the susceptibility to common complex diseases. To fully characterize their roles in deriving the complex traits/diseases, it is necessary to annotate all CREs and CRMs encoded in the human genome. However, the current annotations of CREs and CRMs in the human genome are still very limited and mostly coarse-grained, as they often lack the detailed information of CREs in CRMs. Here, we integrated 620 TF ChIP-seq datasets produced by the ENCODE project for 168 TFs in 79 different cell/tissue types and predicted an unprecedentedly completely map of CREs in CRMs in the human genome at single nucleotide resolution. The map includes 305 912 CRMs containing a total of 1 178 913 CREs belonging to 736 unique TF binding motifs. The predicted CREs and CRMs tend to be subject to either purifying selection or positive selection, thus are likely to be functional. Based on the results, we also examined the status of available ChIP-seq datasets for predicting the entire regulatory genome of humans.

Parent-of-origin-dependent nucleosome organization correlates with genomic imprinting in maize

Genomic imprinting refers to allele-specific expression of genes depending on their parental origin. Nucleosomes, the fundamental units of chromatin, play a critical role in gene transcriptional regulation. However, it remains unknown whether differential nucleosome organization is related to the allele-specific expression of imprinted genes. Here, we generated a genome-wide map of allele-specific nucleosome occupancy in maize endosperm and presented an integrated analysis of its relationship with parent-of-origin-dependent gene expression and DNA methylation. We found that ∼2.3% of nucleosomes showed significant parental bias in maize endosperm. The parent-of-origin-dependent nucleosomes mostly exist as single isolated nucleosomes. Parent-of-origin-dependent nucleosomes were significantly associated with the allele-specific expression of imprinted genes, with nucleosomes positioned preferentially in the promoter of nonexpressed alleles of imprinted genes. Furthermore, we found that most of the paternal specifically positioned nucleosomes (pat-nucleosomes) were associated with parent-of-origin-dependent differential methylated regions, suggesting a functional link between the maternal demethylation and the occurrence of pat-nucleosome. Maternal specifically positioned nucleosomes (mat-nucleosomes) were independent of allele-specific DNA methylation but seem to be associated with allele-specific histone modification. Our study provides the first genome-wide map of allele-specific nucleosome occupancy in plants and suggests a mechanistic connection between chromatin organization and genomic imprinting.

TFAP2C regulates transcription in human naive pluripotency by opening enhancers

Naïve and primed pluripotent hESCs bear transcriptional similarity to pre- and post-implantation epiblast and thus constitute a developmental model for understanding the earliest pluripotent stages in human embryo development. To identify new transcription factors that differentially regulate the unique pluripotent stages, we mapped open chromatin using ATAC-Seq and found enrichment of the AP2 transcription factor binding motif at naïve-specific open chromatin. We determined that the AP2 family member TFAP2C is upregulated during primed to naïve reversion and becomes widespread at naïve-specific enhancers. TFAP2C functions to maintain pluripotency and repress neuroectodermal differentiation during the transition from primed to naïve by facilitating the opening of enhancers proximal to pluripotency factors. Additionally, we identify a previously undiscovered naïve-specific POU5F1 (OCT4) enhancer enriched for TFAP2C binding. Taken together, TFAP2C establishes and maintains naïve human pluripotency and regulates OCT4 expression by mechanisms that are distinct from mouse.

OCEAN-C: mapping hubs of open chromatin interactions across the genome reveals gene regulatory networks

We develop a method called open chromatin enrichment and network Hi-C (OCEAN-C) for antibody-independent mapping of global open chromatin interactions. By integrating FAIRE-seq and Hi-C, OCEAN-C detects open chromatin interactions enriched by active cis-regulatory elements. We identify more than 10,000 hubs of open chromatin interactions (HOCIs) in human cells, which are mainly active promoters and enhancers bound by many DNA-binding proteins and form interaction networks crucial for gene transcription. In addition to identifying large-scale topological structures, including topologically associated domains and A/B compartments, OCEAN-C can detect HOCI-mediated chromatin interactions that are strongly associated with gene expression, super-enhancers, and broad H3K4me3 domains.

Enhancer identification and activity evaluation in the red flour beetle, Tribolium castaneum

Evolution of cis-regulatory elements (such as enhancers) plays an important role in the production of diverse morphology. However, a mechanistic understanding is often limited by the absence of methods for studying enhancers in species other than established model systems. Here, we sought to establish methods to identify and test enhancer activity in the red flour beetle, Tribolium castaneum To identify possible enhancer regions, we first obtained genome-wide chromatin profiles from various tissues and stages of Tribolium using FAIRE (formaldehyde-assisted isolation of regulatory elements)-sequencing. Comparison of these profiles revealed a distinct set of open chromatin regions in each tissue and at each stage. In addition, comparison of the FAIRE data with sets of computationally predicted (i.e. supervised cis-regulatory module-predicted) enhancers revealed a very high overlap between the two datasets. Second, using nubbin in the wing and hunchback in the embryo as case studies, we established the first universal reporter assay system that works in various contexts in Tribolium, and in a cross-species context. Together, these advances will facilitate investigation of cis-evolution and morphological diversity in Tribolium and other insects.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

Androgen Receptor Deregulation Drives Bromodomain-Mediated Chromatin Alterations in Prostate Cancer

Global changes in chromatin accessibility may drive cancer progression by reprogramming transcription factor (TF) binding. In addition, histone acetylation readers such as bromodomain-containing protein 4 (BRD4) have been shown to associate with these TFs and contribute to aggressive cancers including prostate cancer (PC). Here, we show that chromatin accessibility defines castration-resistant prostate cancer (CRPC). We show that the deregulation of androgen receptor (AR) expression is a driver of chromatin relaxation and that AR/androgen-regulated bromodomain-containing proteins (BRDs) mediate this effect. We also report that BRDs are overexpressed in CRPCs and that ATAD2 and BRD2 have prognostic value. Finally, we developed gene stratification signature (BROMO-10) for bromodomain response and PC prognostication, to inform current and future trials with drugs targeting these processes. Our findings provide a compelling rational for combination therapy targeting bromodomains in selected patients in which BRD-mediated TF binding is enhanced or modified as cancer progresses.

Chromatin Regulation of Neuronal Maturation and Plasticity

Neurons are dynamic cells that respond and adapt to stimuli throughout their long postmitotic lives. The structural and functional plasticity of neurons requires the regulated transcription of new gene products, and dysregulation of transcription in either the developing or adult brain impairs cognition. We discuss how mechanisms of chromatin regulation help to orchestrate the transcriptional programs that underlie the maturation of developing neurons and the plasticity of adult neurons. We review how chromatin regulation acts locally to modulate the expression of specific genes and more broadly to coordinate gene expression programs during transitions between cellular states. These data highlight the importance of epigenetic transcriptional mechanisms in postmitotic neurons. We suggest areas where emerging methods may advance understanding in the future.

Pioneer transcription factors shape the epigenetic landscape

Pioneer transcription factors have the unique and important role of unmasking chromatin domains during development to allow the implementation of new cellular programs. Compared with those of other transcription factors, this activity implies that pioneer factors can recognize their target DNA sequences in so-called compacted or "closed" heterochromatin and can trigger remodeling of the adjoining chromatin landscape to provide accessibility to nonpioneer transcription factors. Recent studies identified several steps of pioneer action, namely rapid but weak initial binding to heterochromatin and stabilization of binding followed by chromatin opening and loss of cytosine-phosphate-guanine (CpG) methylation that provides epigenetic memory. Whereas CpG demethylation depends on replication, chromatin opening does not. In this Minireview, we highlight the unique properties of this transcription factor class and the challenges of understanding their mechanism of action.

Relationship between histone modifications and transcription factor binding is protein family specific

The very small fraction of putative binding sites (BSs) that are occupied by transcription factors (TFs) in vivo can be highly variable across different cell types. This observation has been partly attributed to changes in chromatin accessibility and histone modification (HM) patterns surrounding BSs. Previous studies focusing on BSs within DNA regulatory regions found correlations between HM patterns and TF binding specificities. However, a mechanistic understanding of TF-DNA binding specificity determinants is still not available. The ability to predict in vivo TF binding on a genome-wide scale requires the identification of features that determine TF binding based on evolutionary relationships of DNA binding proteins. To reveal protein family-dependent mechanisms of TF binding, we conducted comprehensive comparisons of HM patterns surrounding BSs and non-BSs with exactly matched core motifs for TFs in three cell lines: 33 TFs in GM12878, 37 TFs in K562, and 18 TFs in H1-hESC. These TFs displayed protein family-specific preferences for HM patterns surrounding BSs, with high agreement among cell lines. Moreover, compared to models based on DNA sequence and shape at flanking regions of BSs, HM-augmented quantitative machine-learning methods resulted in increased performance in a TF family-specific manner. Analysis of the relative importance of features in these models indicated that TFs, displaying larger HM pattern differences between BSs and non-BSs, bound DNA in an HM-specific manner on a protein family-specific basis. We propose that TF family-specific HM preferences reveal distinct mechanisms that assist in guiding TFs to their cognate BSs by altering chromatin structure and accessibility.

Genetic risk for Alzheimer's disease is concentrated in specific macrophage and microglial transcriptional networks

Background

Genome-wide association studies of Alzheimer’s disease (AD) have identified a number of significant risk loci, the majority of which lie in non-coding regions of the genome. The lack of causal alleles and considerable polygenicity remains a significant barrier to translation into mechanistic understanding. This includes identifying causal variants and the cell/tissue types in which they operate. A fuller understanding of the cell types and transcriptional networks involved in AD genetic risk mechanisms will provide important insights into pathogenesis.

Methods

We assessed the significance of the overlap between genome-wide significant AD risk variants and sites of open chromatin from data sets representing diverse tissue types. We then focussed on macrophages and microglia to investigate the role of open chromatin sites containing motifs for specific transcription factors. Partitioned heritability using LDscore regression was used to investigate the contribution of specific macrophage and microglia transcription factor motif-containing open chromatin sites to the heritability of AD.

Results

AD risk single nucleotide polymorphisms (SNPs) are preferentially located at sites of open chromatin in immune cells, particularly monocytes (z score = 4.43; corrected P = 5.88 × 10^{− 3}). Similar enrichments are observed for macrophages (z score = 4.10; corrected P < 2.40 × 10^{− 3}) and microglia (z score = 4.34, corrected P = 0.011). In both macrophages and microglia, AD risk variants are enriched at a subset of open chromatin sites that contain DNA binding motifs for specific transcription factors, e.g. SPI1 and MEF2. Genetic variation at many of these motif-containing sites also mediate a substantial proportion of AD heritability, with SPI1-containing sites capturing the majority of the common variant SNP-chip heritability (microglia enrichment = 16.28, corrected enrichment P = 0.0044).

Conclusions

AD risk alleles plausibly operate in immune cells, including microglia, and are concentrated in specific transcriptional networks. Combined with primary genetic association results, the SPI1 and MEF2 transcriptional networks appear central to AD risk mechanisms. Investigation of transcription factors targeting AD risk SNP associated regulatory elements could provide powerful insights into the molecular processes affected by AD polygenic risk. More broadly, our findings support a model of polygenic disease risk that arises from variants located in specific transcriptional networks.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0523-8) contains supplementary material, which is available to authorized users.

Molecular characterization of gene regulatory networks in primary human tracheal and bronchial epithelial cells

BACKGROUND

Robust methods to culture primary airway epithelial cells were developed several decades ago and these cells provide the model of choice to investigate many diseases of the human lung. However, the molecular signature of cells from different regions of the airway epithelium has not been well characterized.

METHODS

We utilize DNase-seq and RNA-seq to examine the molecular signatures of primary cells derived from human tracheal and bronchial tissues, as well as healthy and diseased (cystic fibrosis (CF)) donor lung tissue.

RESULTS

Our data reveal an airway cell signature that is divergent from other epithelial cell types and from common airway epithelial cell lines. The differences between tracheal and bronchial cells are clearly evident as are common regulatory features. Only minor variation is seen between bronchial cells from healthy or CF donors.

CONCLUSIONS

These data are a valuable resource for functional genomics analysis of airway epithelial tissues in human disease.

Computational Methods for Assessing Chromatin Hierarchy

The hierarchical organization of chromatin is known to associate with diverse cellular functions; however, the precise mechanisms and the 3D structure remain to be determined. With recent advances in high-throughput next generation sequencing (NGS) techniques, genome-wide profiling of chromatin structures is made possible. Here, we provide a comprehensive overview of NGS-based methods for profiling "higher-order" and "primary-order" chromatin structures from both experimental and computational aspects. Experimental requirements and considerations specific for each method were highlighted. For computational analysis, we summarized a common analysis strategy for both levels of chromatin assessment, focusing on the characteristic computing steps and the tools. The recently developed single-cell level techniques based on Hi-C and ATAC-seq present great potential to reveal cell-to-cell variability in chromosome architecture. A brief discussion on these methods in terms of experimental and data analysis features is included. We also touch upon the biological relevance of chromatin organization and how the combination with other techniques uncovers the underlying mechanisms. We conclude with a summary and our prospects on necessary improvements of currently available methods in order to advance understanding of chromatin hierarchy. Our review brings together the analyses of both higher- and primary-order chromatin structures, and serves as a roadmap when choosing appropriate experimental and computational methods for assessing chromatin hierarchy.

Estimating the mutational load for cardiovascular diseases in Pakistani population

The deleterious genetic variants contributing to certain diseases may differ in terms of number and allele frequency from population to population depending on their evolutionary background. Here, we prioritize the deleterious variants from Pakistani population in manually curated gene list already reported to be associated with common, Mendelian, and congenital cardiovascular diseases (CVDs) using the genome/exome sequencing data of Pakistani individuals publically available in 1000 Genomes Project (PJL), and Exome Aggregation Consortium (ExAC) South Asia. By applying a set of tools such as Combined Annotation Dependent Depletion (CADD), ANNOVAR, and Variant Effect Predictor (VEP), we highlighted 561 potentially detrimental variants from PJL data, and 7374 variants from ExAC South Asian data. Likewise, filtration from ClinVar for CVDs revealed 03 pathogenic and 02 likely pathogenic variants from PJL and 112 pathogenic and 42 likely pathogenic variants from ExAC South Asians. The comparison of derived allele frequencies (DAF) revealed many of these prioritized variants having two fold and higher DAF in Pakistani individuals than in other populations. The highest number of deleterious variants contributing to common CVDs in descending order includes hypertension, atherosclerosis, heart failure, aneurysm, and coronary heart disease, and for Mendelian and congenital CVDs cardiomyopathies, cardiac arrhythmias, and atrioventricular septal defects.

Chromatin proteomics reveals novel combinatorial histone modification signatures that mark distinct subpopulations of macrophage enhancers

The integrated activity of cis-regulatory elements fine-tunes transcriptional programs of mammalian cells by recruiting cell type–specific as well as ubiquitous transcription factors (TFs). Despite their key role in modulating transcription, enhancers are still poorly characterized at the molecular level, and their limited DNA sequence conservation in evolution and variable distance from target genes make their unbiased identification challenging. The coexistence of high mono-methylation and low tri-methylation levels of lysine 4 of histone H3 is considered a signature of enhancers, but a comprehensive view of histone modifications associated to enhancers is still lacking. By combining chromatin immunoprecipitation (ChIP) with mass spectrometry, we investigated cis-regulatory regions in macrophages to comprehensively identify histone marks specifically associated with enhancers, and to profile their dynamics after transcriptional activation elicited by an inflammatory stimulation. The intersection of the proteomics data with ChIP-seq and RNA-seq analyses revealed the existence of novel subpopulations of enhancers, marked by specific histone modification signatures: specifically, H3K4me1/K36me2 marks transcribed enhancers, while H3K4me1/K36me3 and H3K4me1/K79me2 combinations mark distinct classes of intronic enhancers. Thus, our MS analysis of functionally distinct genomic regions revealed the combinatorial code of histone modifications, highlighting the potential of proteomics in addressing fundamental questions in epigenetics.

Genome-wide mapping of DNase I hypersensitive sites in rare cell populations using single-cell DNase sequencing

Increased chromatin accessibility is a feature of cell-type-specific cis-regulatory elements; therefore, mapping of DNase I hypersensitive sites (DHSs) enables the detection of active regulatory elements of transcription, including promoters, enhancers, insulators and locus-control regions. Single-cell DNase sequencing (scDNase-seq) is a method of detecting genome-wide DHSs when starting with either single cells or <1,000 cells from primary cell sources. This technique enables genome-wide mapping of hypersensitive sites in a wide range of cell populations that cannot be analyzed using conventional DNase I sequencing because of the requirement for millions of starting cells. Fresh cells, formaldehyde-cross-linked cells or cells recovered from formalin-fixed paraffin-embedded (FFPE) tissue slides are suitable for scDNase-seq assays. To generate scDNase-seq libraries, cells are lysed and then digested with DNase I. Circular carrier plasmid DNA is included during subsequent DNA purification and library preparation steps to prevent loss of the small quantity of DHS DNA. Libraries are generated for high-throughput sequencing on the Illumina platform using standard methods. Preparation of scDNase-seq libraries requires only 2 d. The materials and molecular biology techniques described in this protocol should be accessible to any general molecular biology laboratory. Processing of high-throughput sequencing data requires basic bioinformatics skills and uses publicly available bioinformatics software.

Redundant CArG Box Cis-motif Activity Mediates SHATTERPROOF2 Transcriptional Regulation during Arabidopsis thaliana Gynoecium Development

In the Arabidopsis thaliana seed pod, pod shatter and seed dispersal properties are in part determined by the development of a longitudinally orientated dehiscence zone (DZ) that derives from cells of the gynoecial valve margin (VM). Transcriptional regulation of the MADS protein encoding transcription factors genes SHATTERPROOF1 (SHP1) and SHATTERPROOF2 (SHP2) are critical for proper VM identity specification and later on for DZ development. Current models of SHP1 and SHP2 regulation indicate that the transcription factors FRUITFULL (FUL) and REPLUMLESS (RPL) repress these SHP genes in the developing valve and replum domains, respectively. Thus the expression of the SHP genes is restricted to the VM. FUL encodes a MADS-box containing transcription factor that is predicted to act through CArG-box containing cis-regulatory motifs. Here we delimit functional modules within the SHP2 cis-regulatory region and examine the functional importance of CArG box motifs within these regulatory regions. We have characterized a 2.2kb region upstream of the SHP2 translation start site that drives early and late medial domain expression in the gynoecium, as well as expression within the VM and DZ. We identified two separable, independent cis-regulatory modules, a 1kb promoter region and a 700bp enhancer region, that are capable of giving VM and DZ expression. Our results argue for multiple independent cis-regulatory modules that support SHP2 expression during VM development and may contribute to the robustness of SHP2 expression in this tissue. Additionally, three closely positioned CArG box motifs located in the SHP2 upstream regulatory region were mutated in the context of the 2.2kb reporter construct. Mutating simultaneously all three CArG boxes caused a moderate de-repression of the SHP2 reporter that was detected within the valve domain, suggesting that these CArG boxes are involved in SHP2 repression in the valve.

Identification of cis-regulatory sequences reveals potential participation of lola and Deaf1 transcription factors in Anopheles gambiae innate immune response

The innate immune response of Anopheles gambiae involves the transcriptional upregulation of effector genes. Therefore, the cis-regulatory sequences and their cognate binding factors play essential roles in the mosquito’s immune response. However, the genetic control of the mosquito’s innate immune response is not yet fully understood. To gain further insight on the elements, the factors and the potential mechanisms involved, an open chromatin profiling was carried out on A. gambiae-derived immune-responsive cells. Here, we report the identification of cis-regulatory sites, immunity-related transcription factor binding sites, and cis-regulatory modules. A de novo motif discovery carried out on this set of cis-regulatory sequences identified immunity-related motifs and cis-regulatory modules. These modules contain motifs that are similar to binding sites for REL-, STAT-, lola- and Deaf1-type transcription factors. Sequence motifs similar to the binding sites for GAGA were found within a cis-regulatory module, together with immunity-related transcription factor binding sites. The presence of Deaf1- and lola-type binding sites, along with REL- and STAT-type binding sites, suggests that the immunity function of these two factors could have been conserved both in Drosophila and Anopheles gambiae.

Genetic and epigenetic determinants of inter-individual variability in responses to toxicants

It is well established that genetic variability has a major impact on susceptibility to common diseases, responses to drugs and toxicants, and influences disease-related outcomes. The appreciation that epigenetic marks also vary across the population is growing with more data becoming available from studies in humans and model organisms. In addition, the links between genetic variability, toxicity outcomes and epigenetics are being actively explored. Recent studies demonstrate that gene-by-environment interactions involve both chromatin states and transcriptional regulation, and that epigenetics provides important mechanistic clues to connect expression-related quantitative trait loci (QTL) and disease outcomes. However, studies of Gene×Environment×Epigenetics further extend the complexity of the experimental designs and create a challenge for selecting the most informative epigenetic readouts that can be feasibly performed to interrogate multiple individuals, exposures, tissue types and toxicity phenotypes. We propose that among the many possible epigenetic experimental methodologies, assessment of chromatin accessibility coupled with total RNA levels provides a cost-effective and comprehensive option to sufficiently characterize the complexity of epigenetic and regulatory activity in the context of understanding the inter-individual variability in responses to toxicants.

Analysis of chromatin accessibility in human epidermis identifies putative barrier dysfunction-sensing enhancers

To identify putative gene regulatory regions that respond to epidermal injury, we mapped chromatin dynamics in a stratified human epidermis during barrier maturation and disruption. Engineered skin substitutes (ESS) cultured at the air-liquid interface were used as a model of developing human epidermis with incomplete barrier formation. The epidermal barrier stabilized following engraftment onto immunocompromised mice, and was compromised again upon injury. Modified formaldehyde-assisted isolation of regulatory elements (FAIRE) was used to identify accessible genomic regions characteristic of monolayer keratinocytes, ESS in vitro, grafted ESS, and tape-stripped ESS graft. We mapped differentiation- and maturation-associated changes in transcription factor binding sites enriched at each stage and observed overrepresentation of AP-1 gene family motifs in barrier-deficient samples. Transcription of TSLP, an important effector of immunological memory in response to allergen exposure, was dramatically elevated in our barrier-deficient samples. We identified dynamic DNA elements that correlated with TSLP induction and may contain enhancers that regulate TSLP. Two dynamic regions were located near the TSLP promoter and overlapped with allergy-associated SNPs rs17551370 and rs2289877, strongly implicating these loci in the regulation of TSLP expression in allergic disease. Additional dynamic chromatin regions ~250kb upstream of the TSLP promoter were found to be in high linkage disequilibrium with allergic disease SNPs. Taken together, these results define dynamic chromatin accessibility changes during epidermal development and dysfunction.

Topological organization and dynamic regulation of human tRNA genes during macrophage differentiation

Background

The human genome is hierarchically organized into local and long-range structures that help shape cell-type-specific transcription patterns. Transfer RNA (tRNA) genes (tDNAs), which are transcribed by RNA polymerase III (RNAPIII) and encode RNA molecules responsible for translation, are dispersed throughout the genome and, in many cases, linearly organized into genomic clusters with other tDNAs. Whether the location and three-dimensional organization of tDNAs contribute to the activity of these genes has remained difficult to address, due in part to unique challenges related to tRNA sequencing. We therefore devised integrated tDNA expression profiling, a method that combines RNAPIII mapping with biotin-capture of nascent tRNAs. We apply this method to the study of dynamic tRNA gene regulation during macrophage development and further integrate these data with high-resolution maps of 3D chromatin structure.

Results

Integrated tDNA expression profiling reveals domain-level and loop-based organization of tRNA gene transcription during cellular differentiation. tRNA genes connected by DNA loops, which are proximal to CTCF binding sites and expressed at elevated levels compared to non-loop tDNAs, change coordinately with tDNAs and protein-coding genes at distal ends of interactions mapped by in situ Hi-C. We find that downregulated tRNA genes are specifically marked by enhanced promoter-proximal binding of MAF1, a transcriptional repressor of RNAPIII activity, altogether revealing multiple levels of tDNA regulation during cellular differentiation.

Conclusions

We present evidence of both local and coordinated long-range regulation of human tDNA expression, suggesting the location and organization of tRNA genes contribute to dynamic tDNA activity during macrophage development.

Electronic supplementary material

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

The Hox proteins Ubx and AbdA collaborate with the transcription pausing factor M1BP to regulate gene transcription

In metazoans, the pausing of RNA polymerase II at the promoter (paused Pol II) has emerged as a widespread and conserved mechanism in the regulation of gene transcription. While critical in recruiting Pol II to the promoter, the role transcription factors play in transitioning paused Pol II into productive Pol II is, however, little known. By studying how Drosophila Hox transcription factors control transcription, we uncovered a molecular mechanism that increases productive transcription. We found that the Hox proteins AbdA and Ubx target gene promoters previously bound by the transcription pausing factor M1BP, containing paused Pol II and enriched with promoter-proximal Polycomb Group (PcG) proteins, yet lacking the classical H3K27me3 PcG signature. We found that AbdA binding to M1BP-regulated genes results in reduction in PcG binding, the release of paused Pol II, increases in promoter H3K4me3 histone marks and increased gene transcription. Linking transcription factors, PcG proteins and paused Pol II states, these data identify a two-step mechanism of Hox-driven transcription, with M1BP binding leading to Pol II recruitment followed by AbdA targeting, which results in a change in the chromatin landscape and enhanced transcription.

Identification of breast cancer associated variants that modulate transcription factor binding

Genome-wide association studies (GWAS) have discovered thousands loci associated with disease risk and quantitative traits, yet most of the variants responsible for risk remain uncharacterized. The majority of GWAS-identified loci are enriched for non-coding single-nucleotide polymorphisms (SNPs) and defining the molecular mechanism of risk is challenging. Many non-coding causal SNPs are hypothesized to alter transcription factor (TF) binding sites as the mechanism by which they affect organismal phenotypes. We employed an integrative genomics approach to identify candidate TF binding motifs that confer breast cancer-specific phenotypes identified by GWAS. We performed de novo motif analysis of regulatory elements, analyzed evolutionary conservation of identified motifs, and assayed TF footprinting data to identify sequence elements that recruit TFs and maintain chromatin landscape in breast cancer-relevant tissue and cell lines. We identified candidate causal SNPs that are predicted to alter TF binding within breast cancer-relevant regulatory regions that are in strong linkage disequilibrium with significantly associated GWAS SNPs. We confirm that the TFs bind with predicted allele-specific preferences using CTCF ChIP-seq data. We used The Cancer Genome Atlas breast cancer patient data to identify ANKLE1 and ZNF404 as the target genes of candidate TF binding site SNPs in the 19p13.11 and 19q13.31 GWAS-identified loci. These SNPs are associated with the expression of ZNF404 and ANKLE1 in breast tissue. This integrative analysis pipeline is a general framework to identify candidate causal variants within regulatory regions and TF binding sites that confer phenotypic variation and disease risk.

TMPRSS2-ERG fusion co-opts master transcription factors and activates NOTCH signaling in primary prostate cancer

TMPRSS2-ERG (T2E) structural rearrangements typify ∼50% of prostate tumors and result in overexpression of the ERG transcription factor. Using chromatin, genomic and expression data, we show distinct cis-regulatory landscapes between T2E-positive and non-T2E primary prostate tumors, which include clusters of regulatory elements (COREs). This difference is mediated by ERG co-option of HOXB13 and FOXA1, implementing a T2E-specific transcriptional profile. We also report a T2E-specific CORE on the structurally rearranged ERG locus arising from spreading of the TMPRSS2 locus pre-existing CORE, assisting in its overexpression. Finally, we show that the T2E-specific cis-regulatory landscape underlies a vulnerability against the NOTCH pathway. Indeed, NOTCH pathway inhibition antagonizes the growth and invasion of T2E-positive prostate cancer cells. Taken together, our work shows that overexpressed ERG co-opts master transcription factors to deploy a unique cis-regulatory landscape, inducing a druggable dependency on NOTCH signaling in T2E-positive prostate tumors.

Prediction of Chromatin Accessibility in Gene-Regulatory Regions from Transcriptomics Data

The epigenetics landscape of cells plays a key role in the establishment of cell-type specific gene expression programs characteristic of different cellular phenotypes. Different experimental procedures have been developed to obtain insights into the accessible chromatin landscape including DNase-seq, FAIRE-seq and ATAC-seq. However, current downstream computational tools fail to reliably determine regulatory region accessibility from the analysis of these experimental data. In particular, currently available peak calling algorithms are very sensitive to their parameter settings and show highly heterogeneous results, which hampers a trustworthy identification of accessible chromatin regions. Here, we present a novel method that predicts accessible and, more importantly, inaccessible gene-regulatory chromatin regions solely relying on transcriptomics data, which complements and improves the results of currently available computational methods for chromatin accessibility assays. We trained a hierarchical classification tree model on publicly available transcriptomics and DNase-seq data and assessed the predictive power of the model in six gold standard datasets. Our method increases precision and recall compared to traditional peak calling algorithms, while its usage is not limited to the prediction of accessible and inaccessible gene-regulatory chromatin regions, but constitutes a helpful tool for optimizing the parameter settings of peak calling methods in a cell type specific manner.

Genome-wide mapping of DNase I hypersensitive sites reveals chromatin accessibility changes in Arabidopsis euchromatin and heterochromatin regions under extended darkness

Light, as the energy source in photosynthesis, is essential for plant growth and development. Extended darkness causes dramatic gene expression changes. In this study, we applied DNase-seq (DNase I hypersensitive site sequencing) to study changes of chromatin accessibility in euchromatic and heterochromatic regions under extended darkness in Arabidopsis. We generated 27 Gb DNase-seq and 67.6 Gb RNA-seq data to investigate chromatin accessibility changes and global gene expression under extended darkness and control condition in Arabidopsis. We found that ~40% DHSs (DNaseI hypersensitive sites) were diminished under darkness. In non-TE regions, the majority of DHS-changed genes were DHS-diminished under darkness. A total of 519 down-regulated genes were associated with diminished DHSs under darkness, mainly involved in photosynthesis process and retrograde signaling, and were regulated by chloroplast maintenance master regulators such as GLK1. In TE regions, approximately half of the DHS-changed TEs were DHS-increased under darkness and were primarily associated with the LTR/Gypsy retrotransposons in the heterochromatin flanking the centromeres. In contrast, DHS-diminished TEs under darkness were enriched in Copia, LINE, and MuDR dispersed across chromosomes. Together, our results indicated that extended darkness resulted in more increased chromatin compaction in euchromatin and decompaction in heterochromatin, thus further leading to gene expression changes in Arabidopsis.

Dosage compensation and sex-specific epigenetic landscape of the X chromosome in the pea aphid

Background

Heterogametic species display a differential number of sex chromosomes resulting in imbalanced transcription levels for these chromosomes between males and females. To correct this disequilibrium, dosage compensation mechanisms involving gene expression and chromatin accessibility regulations have emerged throughout evolution. In insects, these mechanisms have been extensively characterized only in Drosophila but not in insects of agronomical importance. Aphids are indeed major pests of a wide range of crops. Their remarkable ability to switch from asexual to sexual reproduction during their life cycle largely explains the economic losses they can cause. As heterogametic insects, male aphids are X0, while females (asexual and sexual) are XX.

Results

Here, we analyzed transcriptomic and open chromatin data obtained from whole male and female individuals to evaluate the putative existence of a dosage compensation mechanism involving differential chromatin accessibility of the pea aphid’s X chromosome. Transcriptomic analyses first showed X/AA and XX/AA expression ratios for expressed genes close to 1 in males and females, respectively, suggesting dosage compensation in the pea aphid. Analyses of open chromatin data obtained by Formaldehyde-Assisted Isolation of Regulatory Elements (FAIRE-seq) revealed a X chromosome chromatin accessibility globally and significantly higher in males than in females, while autosomes’ chromatin accessibility is similar between sexes. Moreover, chromatin environment of X-linked genes displaying similar expression levels in males and females—and thus likely to be compensated—is significantly more accessible in males.

Conclusions

Our results suggest the existence of an underlying epigenetic mechanism enhancing the X chromosome chromatin accessibility in males to allow X-linked gene dose correction between sexes in the pea aphid, similar to Drosophila. Our study gives new evidence into the comprehension of dosage compensation in link with chromatin biology in insects and newly in a major crop pest, taking benefits from both transcriptomic and open chromatin data.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0137-1) contains supplementary material, which is available to authorized users.

WSMD: weakly-supervised motif discovery in transcription factor ChIP-seq data

Although discriminative motif discovery (DMD) methods are promising for eliciting motifs from high-throughput experimental data, due to consideration of computational expense, most of existing DMD methods have to choose approximate schemes that greatly restrict the search space, leading to significant loss of predictive accuracy. In this paper, we propose Weakly-Supervised Motif Discovery (WSMD) to discover motifs from ChIP-seq datasets. In contrast to the learning strategies adopted by previous DMD methods, WSMD allows a "global" optimization scheme of the motif parameters in continuous space, thereby reducing the information loss of model representation and improving the quality of resultant motifs. Meanwhile, by exploiting the connection between DMD framework and existing weakly supervised learning (WSL) technologies, we also present highly scalable learning strategies for the proposed method. The experimental results on both real ChIP-seq datasets and synthetic datasets show that WSMD substantially outperforms former DMD methods (including DREME, HOMER, XXmotif, motifRG and DECOD) in terms of predictive accuracy, while also achieving a competitive computational speed.

DeFCoM: analysis and modeling of transcription factor binding sites using a motif-centric genomic footprinter

Motivation

Identifying the locations of transcription factor binding sites is critical for understanding how gene transcription is regulated across different cell types and conditions. Chromatin accessibility experiments such as DNaseI sequencing (DNase-seq) and Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) produce genome-wide data that include distinct 'footprint' patterns at binding sites. Nearly all existing computational methods to detect footprints from these data assume that footprint signals are highly homogeneous across footprint sites. Additionally, a comprehensive and systematic comparison of footprinting methods for specifically identifying which motif sites for a specific factor are bound has not been performed.

Results

Using DNase-seq data from the ENCODE project, we show that a large degree of previously uncharacterized site-to-site variability exists in footprint signal across motif sites for a transcription factor. To model this heterogeneity in the data, we introduce a novel, supervised learning footprinter called Detecting Footprints Containing Motifs (DeFCoM). We compare DeFCoM to nine existing methods using evaluation sets from four human cell-lines and eighteen transcription factors and show that DeFCoM outperforms current methods in determining bound and unbound motif sites. We also analyze the impact of several biological and technical factors on the quality of footprint predictions to highlight important considerations when conducting footprint analyses and assessing the performance of footprint prediction methods. Finally, we show that DeFCoM can detect footprints using ATAC-seq data with similar accuracy as when using DNase-seq data.

Availability and Implementation

Python code available at https://bitbucket.org/bryancquach/defcom.

Contact

bquach@email.unc.edu or tsfurey@email.unc.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

A computational approach for the functional classification of the epigenome

BACKGROUND

In the last decade, advanced functional genomics approaches and deep sequencing have allowed large-scale mapping of histone modifications and other epigenetic marks, highlighting functional relationships between chromatin organization and genome function. Here, we propose a novel approach to explore functional interactions between different epigenetic modifications and extract combinatorial profiles that can be used to annotate the chromatin in a finite number of functional classes. Our method is based on non-negative matrix factorization (NMF), an unsupervised learning technique originally employed to decompose high-dimensional data in a reduced number of meaningful patterns. We applied the NMF algorithm to a set of different epigenetic marks, consisting of ChIP-seq assays for multiple histone modifications, Pol II binding and chromatin accessibility assays from human H1 cells.

RESULTS

We identified a number of chromatin profiles that contain functional information and are biologically interpretable. We also observe that epigenetic profiles are characterized by specific genomic contexts and show significant association with distinct genomic features. Moreover, analysis of RNA-seq data reveals that distinct chromatin signatures correlate with the level of gene expression.

CONCLUSIONS

Overall, our study highlights the utility of NMF in studying functional relationships between different epigenetic modifications and may provide new biological insights for the interpretation of the chromatin dynamics.

PRC2 represses transcribed genes on the imprinted inactive X chromosome in mice

Background

Polycomb repressive complex 2 (PRC2) catalyzes histone H3K27me3, which marks many transcriptionally silent genes throughout the mammalian genome. Although H3K27me3 is associated with silenced gene expression broadly, it remains unclear why some but not other PRC2 target genes require PRC2 and H3K27me3 for silencing.

Results

Here we define the transcriptional and chromatin features that predict which PRC2 target genes require PRC2/H3K27me3 for silencing by interrogating imprinted mouse X-chromosome inactivation. H3K27me3 is enriched at promoters of silenced genes across the inactive X chromosome. To abrogate PRC2 function, we delete the core PRC2 protein EED in F1 hybrid trophoblast stem cells (TSCs), which undergo imprinted inactivation of the paternally inherited X chromosome. Eed^–/– TSCs lack H3K27me3 and Xist lncRNA enrichment on the inactive X chromosome. Despite the absence of H3K27me3 and Xist RNA, only a subset of the inactivated X-linked genes is derepressed in Eed^–/– TSCs. Unexpectedly, in wild-type (WT) TSCs these genes are transcribed and are enriched for active chromatin hallmarks on the inactive-X, including RNA PolII, H3K27ac, and H3K36me3, but not the bivalent mark H3K4me2. By contrast, PRC2 targets that remain repressed in Eed^–/– TSCs are depleted for active chromatin characteristics in WT TSCs.

Conclusions

A comparative analysis of transcriptional and chromatin features of inactive X-linked genes in WT and Eed^–/– TSCs suggests that PRC2 acts as a brake to prevent induction of transcribed genes on the inactive X chromosome, a mode of PRC2 function that may apply broadly.

Electronic supplementary material

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

Dynamics of chromatin accessibility during TGF-β-induced EMT of Ras-transformed mammary gland epithelial cells

Epithelial-mesenchymal transition (EMT) is induced by transforming growth factor (TGF)-β and facilitates tumor progression. We here performed global mapping of accessible chromatin in the mouse mammary gland epithelial EpH4 cell line and its Ras-transformed derivative (EpRas) using formaldehyde-assisted isolation of regulatory element (FAIRE)-sequencing. TGF-β and Ras altered chromatin accessibility either cooperatively or independently, and AP1, ETS, and RUNX binding motifs were enriched in the accessible chromatin regions of EpH4 and EpRas cells. Etv4, an ETS family oncogenic transcription factor, was strongly expressed and bound to more than one-third of the accessible chromatin regions in EpRas cells treated with TGF-β. While knockdown of Etv4 and another ETS family member Etv5 showed limited effects on the decrease in the E-cadherin abundance and stress fiber formation by TGF-β, gene ontology analysis showed that genes encoding extracellular proteins were most strongly down-regulated by Etv4 and Etv5 siRNAs. Accordingly, TGF-β-induced expression of Mmp13 and cell invasiveness were suppressed by Etv4 and Etv5 siRNAs, which were accompanied by the reduced chromatin accessibility at an enhancer region of Mmp13 gene. These findings suggest a mechanism of transcriptional regulation during Ras- and TGF-β-induced EMT that involves alterations of accessible chromatin, which are partly regulated by Etv4 and Etv5.

A time-series analysis of altered histone H3 acetylation and gene expression during the course of MMAIII-induced malignant transformation of urinary bladder cells

Our previous studies have shown that chronic exposure to low doses of monomethylarsonous acid (MMAIII) causes global histone acetylation dysregulation in urothelial cells (UROtsa cells) during the course of malignant transformation. To reveal the relationship between altered histone acetylation patterns and aberrant gene expression, more specifically, the carcinogenic relevance of these alterations, we performed a time-course analysis of the binding patterns of histone 3 lysine 18 acetylation (H3K18ac) across the genome and generated global gene-expression profiles from this UROtsa cell malignant transformation model. We showed that H3K18ac, one of the most significantly upregulated histone acetylation sites following MMAIII exposure, was enriched at gene promoter-specific regions across the genome and that MMAIII-induced upregulation of H3K18ac led to an altered binding pattern in a large number of genes that was most significant during the critical window for MMAIII-induced UROtsa cells' malignant transformation. Some genes identified as having a differential binding pattern with H3K18ac, acted as upstream regulators of critical gene networks with known functions in tumor development and progression. The altered H3K18ac binding patterns not only led to changes in expression of these directly affected upstream regulators but also resulted in gene-expression changes in their regulated networks. Collectively, our data suggest that MMAIII-induced alteration of histone acetylation patterns in UROtsa cells led to a time- and malignant stage-dependent aberrant gene-expression pattern, and that some gene regulatory networks were altered in accordance with their roles in carcinogenesis, probably contributing to MMAIII-induced urothelial cell malignant transformation and carcinogenesis.

Inference of cell type specific regulatory networks on mammalian lineages

Transcriptional regulatory networks are at the core of establishing cell type specific gene expression programs. In mammalian systems, such regulatory networks are determined by multiple levels of regulation, including by transcription factors, chromatin environment, and three-dimensional organization of the genome. Recent efforts to measure diverse regulatory genomic datasets across multiple cell types and tissues offer unprecedented opportunities to examine the context-specificity and dynamics of regulatory networks at a greater resolution and scale than before. In parallel, numerous computational approaches to analyze these data have emerged that serve as important tools for understanding mammalian cell type specific regulation. In this article, we review recent computational approaches to predict the expression and sequence-based regulators of a gene's expression level and examine long-range gene regulation. We highlight promising approaches, insights gained, and open challenges that need to be overcome to build a comprehensive picture of cell type specific transcriptional regulatory networks.

CRISPR-Cas9 epigenome editing enables high-throughput screening for functional regulatory elements in the human genome

Large genome-mapping consortia and thousands of genome-wide association studies have identified non-protein coding elements in the genome as a having a central role in various biological processes. However, decoding the function of the millions of putative regulatory elements discovered in these studies remains challenging. CRISPR–Cas9-based epigenome editing technologies have enabled precise perturbation of the activity of specific regulatory elements. Here we describe CRISPR–Cas9-based epigenomic regulatory element screening (CERES) for improved high-throughput screening of regulatory element activity within the native genomic context. Using dCas9^KRAB repressor and dCas9^p300 activator constructs and lentiviral sgRNA libraries targeting DNase I hypersensitive sites surrounding a gene of interest, we perform both loss- and gain-of-function screens to identify regulatory elements for the β-globin and the HER2 loci in human cells. CERES readily identified known and novel regulatory elements, some of which were dependent on cell type or direction of perturbation. This technology allows the high-throughput functional annotation of putative regulatory elements in their native chromosomal context.

Transcription factor-DNA binding: beyond binding site motifs

Sequence-specific transcription factors (TFs) regulate gene expression by binding to cis-regulatory elements in promoter and enhancer DNA. While studies of TF-DNA binding have focused on TFs' intrinsic preferences for primary nucleotide sequence motifs, recent studies have elucidated additional layers of complexity that modulate TF-DNA binding. In this review, we discuss technological developments for identifying TF binding preferences and highlight recent discoveries that elaborate how TF interactions, local DNA structure, and genomic features influence TF-DNA binding. We highlight novel approaches for characterizing functional binding site motifs that promise to inform our understanding of how TF binding controls gene expression and ultimately contributes to phenotype.

Heterarchy of transcription factors driving basal and luminal cell phenotypes in human urothelium

Cell differentiation is affected by complex networks of transcription factors that co-ordinate re-organisation of the chromatin landscape. The hierarchies of these relationships can be difficult to dissect. During in vitro differentiation of normal human uro-epithelial cells, formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) and RNA-seq was used to identify alterations in chromatin accessibility and gene expression changes following activation of the nuclear receptor peroxisome proliferator-activated receptor gamma (PPARγ) as a differentiation-initiating event. Regions of chromatin identified by FAIRE-seq, as having altered accessibility during differentiation, were found to be enriched with sequence-specific binding motifs for transcription factors predicted to be involved in driving basal and differentiated urothelial cell phenotypes, including forkhead box A1 (FOXA1), P63, GRHL2, CTCF and GATA-binding protein 3 (GATA3). In addition, co-occurrence of GATA3 motifs was observed within subsets of differentiation-specific peaks containing P63 or FOXA1. Changes in abundance of GRHL2, GATA3 and P63 were observed in immunoblots of chromatin-enriched extracts. Transient siRNA knockdown of P63 revealed that P63 favoured a basal-like phenotype by inhibiting differentiation and promoting expression of basal marker genes. GATA3 siRNA prevented differentiation-associated downregulation of P63 protein and transcript, and demonstrated positive feedback of GATA3 on PPARG transcript, but showed no effect on FOXA1 transcript or protein expression. This approach indicates that as a transcriptionally regulated programme, urothelial differentiation operates as a heterarchy, wherein GATA3 is able to co-operate with FOXA1 to drive expression of luminal marker genes, but that P63 has potential to transrepress expression of the same genes.

Cell Type-Specific Epigenomic Analysis Reveals a Uniquely Closed Chromatin Architecture in Mouse Rod Photoreceptors

Rod photoreceptors are specialized neurons that mediate vision in dim light and are the predominant photoreceptor type in nocturnal mammals. The rods of nocturnal mammals are unique among vertebrate cell types in having an ‘inverted’ nuclear architecture, with a dense mass of heterochromatin in the center of the nucleus rather than dispersed clumps at the periphery. To test if this unique nuclear architecture is correlated with a unique epigenomic landscape, we performed ATAC-seq on mouse rods and their most closely related cell type, cone photoreceptors. We find that thousands of loci are selectively closed in rods relative to cones as well as >60 additional cell types. Furthermore, we find that the open chromatin profile of photoreceptors lacking the rod master regulator Nrl is nearly indistinguishable from that of native cones, indicating that Nrl is required for selective chromatin closure in rods. Finally, we identified distinct enrichments of transcription factor binding sites in rods and cones, revealing key differences in the cis-regulatory grammar of these cell types. Taken together, these data provide insight into the development and maintenance of photoreceptor identity, and highlight rods as an attractive system for studying the relationship between nuclear organization and local changes in gene regulation.

Role of the chromatin landscape and sequence in determining cell type-specific genomic glucocorticoid receptor binding and gene regulation

The genomic loci bound by the glucocorticoid receptor (GR), a hormone-activated transcription factor, show little overlap between cell types. To study the role of chromatin and sequence in specifying where GR binds, we used Bayesian modeling within the universe of accessible chromatin. Taken together, our results uncovered that although GR preferentially binds accessible chromatin, its binding is biased against accessible chromatin located at promoter regions. This bias can only be explained partially by the presence of fewer GR recognition sequences, arguing for the existence of additional mechanisms that interfere with GR binding at promoters. Therefore, we tested the role of H3K9ac, the chromatin feature with the strongest negative association with GR binding, but found that this correlation does not reflect a causative link. Finally, we find a higher percentage of promoter-proximal GR binding for genes regulated by GR across cell types than for cell type-specific target genes. Given that GR almost exclusively binds accessible chromatin, we propose that cell type-specific regulation by GR preferentially occurs via distal enhancers, whose chromatin accessibility is typically cell type-specific, whereas ubiquitous target gene regulation is more likely to result from binding to promoter regions, which are often accessible regardless of cell type examined.

Association of 5-hydroxymethylation and 5-methylation of DNA cytosine with tissue-specific gene expression

Differentially methylated or hydroxymethylated regions (DMRs) in mammalian DNA are often associated with tissue-specific gene expression but the functional relationships are still being unraveled. To elucidate these relationships, we studied 16 human genes containing myogenic DMRs by analyzing profiles of their epigenetics and transcription and quantitatively assaying 5-hydroxymethylcytosine (5hmC) and 5-methylcytosine (5mC) at specific sites in these genes in skeletal muscle (SkM), myoblasts, heart, brain, and diverse other samples. Although most human promoters have little or no methylation regardless of expression, more than half of the genes that we chose to study-owing to their myogenic DMRs-overlapped tissue-specific alternative or cryptic promoters displaying corresponding tissue-specific differences in histone modifications. The 5mC levels in myoblast DMRs were significantly associated with 5hmC levels in SkM at the same site. Hypermethylated myogenic DMRs within CDH15, a muscle- and cerebellum-specific cell adhesion gene, and PITX3, a homeobox gene, were used for transfection in reporter gene constructs. These intragenic DMRs had bidirectional tissue-specific promoter activity that was silenced by in vivo-like methylation. The CDH15 DMR, which was previously associated with an imprinted maternal germline DMR in mice, had especially strong promoter activity in myogenic host cells. These findings are consistent with the controversial hypothesis that intragenic DNA methylation can facilitate transcription and is not just a passive consequence of it. Our results support varied roles for tissue-specific 5mC- or 5hmC-enrichment in suppressing inappropriate gene expression from cryptic or alternative promoters and in increasing the plasticity of gene expression required for development and rapid responses to tissue stress or damage.