Genome-wide analysis of the relationships between DNaseI HS, histone modifications and gene expression reveals distinct modes of chromatin domains

Regulation of chromatin accessibility and gene expression in the developing hippocampal primordium by LIM-HD transcription factor LHX2

In the mammalian cerebral cortex, the hippocampal primordium (Hcp) occupies a discrete position in the dorsal telencephalic neuroepithelium adjacent to the neocortical primordium (Ncp). We examined transcriptomic and chromatin-level features that distinguish the Hcp from the Ncp in the mouse during the early neurogenic period, embryonic day (E)12.5. ATAC-seq revealed that the Hcp was more accessible than the Ncp at this stage. Motif analysis of the differentially accessible loci in these tissues revealed LHX2 as a candidate transcription factor for modulating gene regulatory networks (GRNs). We analyzed LHX2 occupancy profiles and compared these with transcriptomic data from control and Lhx2 mutant Hcp and Ncp at E12.5. Our results revealed that LHX2 directly regulates distinct genes in the Hcp and Ncp within a set of common pathways that control fundamental aspects of development namely pluripotency, axon pathfinding, Wnt, and Hippo signaling. Loss of Lhx2 caused a decrease in accessibility, specifically in hippocampal chromatin, suggesting that this factor may play a unique role in hippocampal development. We identified 14 genes that were preferentially enriched in the Hcp, for which LHX2 regulates both chromatin accessibility and mRNA expression, which have not thus far been examined in hippocampal development. Together, these results provide mechanistic insight into how LHX2 function in the Hcp may contribute to the process by which the hippocampus acquires features distinct from the neocortex.

Chromatin Dynamics in Digestive System Cancer: Commander and Regulator

Digestive system tumors have a poor prognosis due to complex anatomy, insidious onset, challenges in early diagnosis, and chemoresistance. Epidemiological statistics has verified that digestive system tumors rank first in tumor-related death. Although a great number of studies are devoted to the molecular biological mechanism, early diagnostic markers, and application of new targeted drugs in digestive system tumors, the therapeutic effect is still not satisfactory. Epigenomic alterations including histone modification and chromatin remodeling are present in human cancers and are now known to cooperate with genetic changes to drive the cancer phenotype. Chromatin is the carrier of genetic information and consists of DNA, histones, non-histone proteins, and a small amount of RNA. Chromatin and nucleosomes control the stability of the eukaryotic genome and regulate DNA processes such as transcription, replication, and repair. The dynamic structure of chromatin plays a key role in this regulatory function. Structural fluctuations expose internal DNA and thus provide access to the nuclear machinery. The dynamic changes are affected by various complexes and epigenetic modifications. Variation of chromatin dynamics produces early and superior regulation of the expression of related genes and downstream pathways, thereby controlling tumor development. Intervention at the chromatin level can change the process of cancer earlier and is a feasible option for future tumor diagnosis and treatment. In this review, we introduced chromatin dynamics including chromatin remodeling, histone modifications, and chromatin accessibility, and current research on chromatin regulation in digestive system tumors was also summarized.

Bioinformatics strategies for studying the molecular mechanisms of fungal extracellular vesicles with a focus on infection and immune responses

Fungal extracellular vesicles (EVs) are released during pathogenesis and are found to be an opportunistic infection in most cases. EVs are immunocompetent with their host and have paved the way for new biomedical approaches to drug delivery and the treatment of complex diseases including cancer. With computing and processing advancements, the rise of bioinformatics tools for the evaluation of various parameters involved in fungal EVs has blossomed. In this review, we have complied and explored the bioinformatics tools to analyze the host-pathogen interaction, toxicity, omics and pathogenesis with an array of specific tools that have depicted the ability of EVs as vector/carrier for therapeutic agents and as a potential theme for immunotherapy. We have also discussed the generation and pathways involved in the production, transport, pathogenic action and immunological interactions of EVs in the host system. The incorporation of network pharmacology approaches has been discussed regarding fungal pathogens and their significance in drug discovery. To represent the overview, we have presented and demonstrated an in silico study model to portray the human Cryptococcal interactions.

Epitome: predicting epigenetic events in novel cell types with multi-cell deep ensemble learning

The accumulation of large epigenomics data consortiums provides us with the opportunity to extrapolate existing knowledge to new cell types and conditions. We propose Epitome, a deep neural network that learns similarities of chromatin accessibility between well characterized reference cell types and a query cellular context, and copies over signal of transcription factor binding and modification of histones from reference cell types when chromatin profiles are similar to the query. Epitome achieves state-of-the-art accuracy when predicting transcription factor binding sites on novel cellular contexts and can further improve predictions as more epigenetic signals are collected from both reference cell types and the query cellular context of interest.

A Newly Assigned Role of CTCF in Cellular Response to Broken DNAs

Best known as a transcriptional factor, CCCTC-binding factor (CTCF) is a highly conserved multifunctional DNA-binding protein with 11 zinc fingers. It functions in diverse genomic processes, including transcriptional activation/repression, insulation, genome imprinting and three-dimensional genome organization. A big surprise has recently emerged with the identification of CTCF engaging in the repair of DNA double-strand breaks (DSBs) and in the maintenance of genome fidelity. This discovery now adds a new dimension to the multifaceted attributes of this protein. CTCF facilitates the most accurate DSB repair via homologous recombination (HR) that occurs through an elaborate pathway, which entails a chain of timely assembly/disassembly of various HR-repair complexes and chromatin modifications and coordinates multistep HR processes to faithfully restore the original DNA sequences of broken DNA sites. Understanding the functional crosstalks between CTCF and other HR factors will illuminate the molecular basis of various human diseases that range from developmental disorders to cancer and arise from impaired repair. Such knowledge will also help understand the molecular mechanisms underlying the diverse functions of CTCF in genome biology. In this review, we discuss the recent advances regarding this newly assigned versatile role of CTCF and the mechanism whereby CTCF functions in DSB repair.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Integrating distal and proximal information to predict gene expression via a densely connected convolutional neural network

MOTIVATION

Interactions among cis-regulatory elements such as enhancers and promoters are main driving forces shaping context-specific chromatin structure and gene expression. Although there have been computational methods for predicting gene expression from genomic and epigenomic information, most of them neglect long-range enhancer-promoter interactions, due to the difficulty in precisely linking regulatory enhancers to target genes. Recently, HiChIP, a novel high-throughput experimental approach, has generated comprehensive data on high-resolution interactions between promoters and distal enhancers. Moreover, plenty of studies suggest that deep learning achieves state-of-the-art performance in epigenomic signal prediction, and thus promoting the understanding of regulatory elements. In consideration of these two factors, we integrate proximal promoter sequences and HiChIP distal enhancer-promoter interactions to accurately predict gene expression.

RESULTS

We propose DeepExpression, a densely connected convolutional neural network, to predict gene expression using both promoter sequences and enhancer-promoter interactions. We demonstrate that our model consistently outperforms baseline methods, not only in the classification of binary gene expression status but also in regression of continuous gene expression levels, in both cross-validation experiments and cross-cell line predictions. We show that the sequential promoter information is more informative than the experimental enhancer information; meanwhile, the enhancer-promoter interactions within ±100 kbp around the TSS of a gene are most beneficial. We finally visualize motifs in both promoter and enhancer regions and show the match of identified sequence signatures with known motifs. We expect to see a wide spectrum of applications using HiChIP data in deciphering the mechanism of gene regulation.

AVAILABILITY AND IMPLEMENTATION

DeepExpression is freely available at https://github.com/wanwenzeng/DeepExpression.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CORE GREML for estimating covariance between random effects in linear mixed models for complex trait analyses

As a key variance partitioning tool, linear mixed models (LMMs) using genome-based restricted maximum likelihood (GREML) allow both fixed and random effects. Classic LMMs assume independence between random effects, which can be violated, causing bias. Here we introduce a generalized GREML, named CORE GREML, that explicitly estimates the covariance between random effects. Using extensive simulations, we show that CORE GREML outperforms the conventional GREML, providing variance and covariance estimates free from bias due to correlated random effects. Applying CORE GREML to UK Biobank data, we find, for example, that the transcriptome, imputed using genotype data, explains a significant proportion of phenotypic variance for height (0.15, p-value = 1.5e-283), and that these transcriptomic effects correlate with the genomic effects (genome-transcriptome correlation = 0.35, p-value = 1.2e-14). We conclude that the covariance between random effects is a key parameter for estimation, especially when partitioning phenotypic variance by multi-omics layers.

Linear mixed models have bias due to the assumed independence between random effects. Here, the authors describe a genome-based restricted maximum likelihood, CORE GREML, which estimates covariance between random effects. Application to UK Biobank data highlights this as an important parameter for multi-omics analyses of phenotypic variance.

A practical guide for DNase-seq data analysis: from data management to common applications

Deoxyribonuclease I (DNase I)-hypersensitive site sequencing (DNase-seq) has been widely used to determine chromatin accessibility and its underlying regulatory lexicon. However, exploring DNase-seq data requires sophisticated downstream bioinformatics analyses. In this study, we first review computational methods for all of the major steps in DNase-seq data analysis, including experimental design, quality control, read alignment, peak calling, annotation of cis-regulatory elements, genomic footprinting and visualization. The challenges associated with each step are highlighted. Next, we provide a practical guideline and a computational pipeline for DNase-seq data analysis by integrating some of these tools. We also discuss the competing techniques and the potential applications of this pipeline for the analysis of analogous experimental data. Finally, we discuss the integration of DNase-seq with other functional genomics techniques.

Maize E2F transcription factors. Expression, association to promoters of S-phase genes and interaction with the RBR1 protein in chromatin during seed germination

For seed germination, it is necessary to restart the cell cycle, a process regulated at multiple levels including transcriptional control, that is executed by the E2F family of transcription factors. We identified 12 genes of the E2F family in maize that are expressed differentially during the first 28 h post imbibition (HAI). E2Fa/b1;1 and E2Fc proteins were characterized as an activator and a putative repressor respectively, both forming heterodimers with DPb2 that bind differentially to consensus E2F response elements in promoters of E2F target genes. Transcripts of target genes for these transcription factors accumulate during germination; in dry seeds E2Fc protein is enriched in the target promoters and is replaced by E2Fa/b1;1 as germination advances. RBR1 is found in the same promoters in non-imbibed and 28 HAI seeds, when DNA replication has concluded, and transcription of the E2F targets should stop. During germination promoters of these target genes seem to be decorated with histone marks related to relaxed chromatin structure. Therefore, E2Fs appear to occupy their target genes in a context of open chromatin, with RBR1 fine tuning the progression between the phases.

Construction of the gene regulatory network identifies MYC as a transcriptional regulator of SWI/SNF complex

Precise positioning of nucleosomes at the gene regulatory elements mediated by the SWI/SNF family of remodelling complex is important for the transcriptional regulation of genes. A wide set of genes are either positively or negatively regulated by SWI/SNF. In higher eukaryotes, around thirty genes were found to code for SWI/SNF subunits. The construction of a gene regulatory network of SWI/SNF subunits identifies MYC as a common regulator for many of the SWI/SNF subunit genes. A meta-analysis study was conducted to investigate the MYC dependent regulation of SWI/SNF remodelling complex. Subunit information and the promoter sequences of the subunit genes were used to find the canonical E-box motif and its variants. Detailed analysis of mouse and human ChIP-Seq at the SWI/SNF subunit loci indicates the presence of MYC binding peaks overlapping with E-boxes. The co-expression correlation and the differential expression analysis of wt vs. MYC perturbed MEFs indicate the MYC dependent regulation of some of the SWI/SNF subunits. The extension of the analysis was done on MYC proficient and MYC deficient embryonic fibroblast cell lines, TGR1 and HO15, and in one of the MYC amplified cancer types, Medulloblastoma. A transcriptional regulatory feedback loop between MYC and SWI/SNF could be a major factor contributing to the aggressiveness of MYC dependent cancers.

Transcriptional Regulation of INSR, the Insulin Receptor Gene

The insulin receptor gene encodes an evolutionarily conserved signaling protein with a wide spectrum of functions in metazoan development. The insulin signaling pathway plays key roles in processes such as metabolic regulation, growth control, and neuronal function. Misregulation of the pathway features in diabetes, cancer, and neurodegenerative diseases, making it an important target for clinical interventions. While much attention has been focused on differential pathway activation through ligand availability, sensitization of overall signaling may also be mediated by differential expression of the insulin receptor itself. Although first characterized as a “housekeeping” gene with stable expression, comparative studies have shown that expression levels of the human INSR mRNA differ by tissue and in response to environmental signals. Our recent analysis of the transcriptional controls affecting expression of the Drosophila insulin receptor gene indicates that a remarkable amount of DNA is dedicated to encoding sophisticated feedback and feed forward signals. The human INSR gene is likely to contain a similar level of transcriptional complexity; here, we summarize over three decades of molecular biology and genetic research that points to a still incompletely understood regulatory control system. Further elucidation of transcriptional controls of INSR will provide the basis for understanding human genetic variation that underlies population-level physiological differences and disease.

Principles and methods of integrative chromatin analysis in primary tissues and tumors

Recent methodological advances have enabled the genome-wide interrogation of chromatin from primary tumor tissues. Integrative analysis of histone post-translational modifications, transcription factor (TF) binding and open chromatin sites in tumors across cancer stages can elucidate the aberrant epigenetic states accompanying tumor progression. Cancer-associated chromatin alterations can activate or inactivate enhancers at genes involved in cancer while still respecting cell-of-origin constrictions. Accordingly, enhancer analysis in cancer could have uses for biomarker discovery to further refine patient diagnosis and potentially sub-classify patients for tailored therapy. Methodologies used for chromatin analyses of primary tissues need to address issues distinct from cell line studies including the specific sources of variability coming from the heterogeneous cellular composition of tissues and from inter-individual (epi)genetic differences. This leads to requirements for careful histological analysis to select the specific samples and cells of interest. In analyzing tumors somatic changes should be taken into account to distinguish the genuine epigenetic changes across tumor specimens from any genetic alterations such as copy number variations (CNV). In this contribution we review a selection of current results from chromatin profiling, examine experimental methodologies and discuss specific analysis approaches. We also review specific considerations regarding tissue preparation for epigenetic analysis and conclude with our perspectives on emerging approaches that will impact studies of chromatin landscapes of clinical samples in the future.

Identification of gene specific cis-regulatory elements during differentiation of mouse embryonic stem cells: An integrative approach using high-throughput datasets

Gene expression governs cell fate, and is regulated via a complex interplay of transcription factors and molecules that change chromatin structure. Advances in sequencing-based assays have enabled investigation of these processes genome-wide, leading to large datasets that combine information on the dynamics of gene expression, transcription factor binding and chromatin structure as cells differentiate. While numerous studies focus on the effects of these features on broader gene regulation, less work has been done on the mechanisms of gene-specific transcriptional control. In this study, we have focussed on the latter by integrating gene expression data for the in vitro differentiation of murine ES cells to macrophages and cardiomyocytes, with dynamic data on chromatin structure, epigenetics and transcription factor binding. Combining a novel strategy to identify communities of related control elements with a penalized regression approach, we developed individual models to identify the potential control elements predictive of the expression of each gene. Our models were compared to an existing method and evaluated using the existing literature and new experimental data from embryonic stem cell differentiation reporter assays. Our method is able to identify transcriptional control elements in a gene specific manner that reflect known regulatory relationships and to generate useful hypotheses for further testing.

Identification and functional evaluation of accessible chromatin associated with wood formation in Eucalyptus grandis

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue.

Widespread epigenetic changes to the enhancer landscape of mouse liver induced by a specific xenobiotic agonist ligand of the nuclear receptor CAR

CAR (Nr1i3), a liver nuclear receptor and xenobiotic sensor, induces drug, steroid and lipid metabolism and dysregulates genes linked to hepatocellular carcinogenesis, but its impact on the liver epigenome is poorly understood. TCPOBOP, a halogenated xenochemical and highly specific CAR agonist ligand, induces localized chromatin opening or closing at several thousand mouse liver genomic regions, discovered as differential DNase-hypersensitive sites (ΔDHS). Active enhancer and promoter histone marks induced by TCPOBOP were enriched at opening DHS and TCPOBOP-inducible genes. Enrichment of CAR binding and CAR motifs was seen at opening DHS and their inducible drug/lipid metabolism gene targets, and at many constitutively open DHS located nearby. TCPOBOP-responsive cell cycle and DNA replication genes co-dependent on MET/EGFR signaling for induction were also enriched for CAR binding. A subset of opening DHS and many closing DHS mapping to TCPOBOP-responsive target genes did not bind CAR, indicating an indirect mechanism for their changes in chromatin accessibility. TCPOBOP-responsive DHS were also enriched for induced binding of RXRA, CEBPA and CEBPB, and for motifs for liver-enriched factors that may contribute to liver-specific transcriptional responses to TCPOBOP exposure. These studies elucidate the enhancer landscape of TCPOBOP-exposed liver and the widespread epigenetic changes that are induced by both direct and indirect mechanisms linked to CAR activation. The global maps of thousands of environmental chemical-induced epigenetic changes described here constitute a rich resource for further research on xenochemical effects on liver chromatin states and the epigenome.

Complex Relationships between Chromatin Accessibility, Sequence Divergence, and Gene Expression in Arabidopsis thaliana

Variation in regulatory DNA is thought to drive phenotypic variation, evolution, and disease. Prior studies of regulatory DNA and transcription factors across animal species highlighted a fundamental conundrum: Transcription factor binding domains and cognate binding sites are conserved, while regulatory DNA sequences are not. It remains unclear how conserved transcription factors and dynamic regulatory sites produce conserved expression patterns across species. Here, we explore regulatory DNA variation and its functional consequences within Arabidopsis thaliana, using chromatin accessibility to delineate regulatory DNA genome-wide. Unlike in previous cross-species comparisons, the positional homology of regulatory DNA is maintained among A. thaliana ecotypes and less nucleotide divergence has occurred. Of the ∼50,000 regulatory sites in A. thaliana, we found that 15% varied in accessibility among ecotypes. Some of these accessibility differences were associated with extensive, previously unannotated sequence variation, encompassing many deletions and ancient hypervariable alleles. Unexpectedly, for the majority of such regulatory sites, nearby gene expression was unaffected. Nevertheless, regulatory sites with high levels of sequence variation and differential chromatin accessibility were the most likely to be associated with differential gene expression. Finally, and most surprising, we found that the vast majority of differentially accessible sites show no underlying sequence variation. We argue that these surprising results highlight the necessity to consider higher-order regulatory context in evaluating regulatory variation and predicting its phenotypic consequences.

GLANET: genomic loci annotation and enrichment tool

Motivation

Genomic studies identify genomic loci representing genetic variations, transcription factor (TF) occupancy, or histone modification through next generation sequencing (NGS) technologies. Interpreting these loci requires evaluating them with known genomic and epigenomic annotations.

Results

We present GLANET as a comprehensive annotation and enrichment analysis tool which implements a sampling-based enrichment test that accounts for GC content and/or mappability biases, jointly or separately. GLANET annotates and performs enrichment analysis on these loci with a rich library. We introduce and perform novel data-driven computational experiments for assessing the power and Type-I error of its enrichment procedure which show that GLANET has attained high statistical power and well-controlled Type-I error rate. As a key feature, users can easily extend its library with new gene sets and genomic intervals. Other key features include assessment of impact of single nucleotide variants (SNPs) on TF binding sites and regulation based pathway enrichment analysis.

Availability and implementation

GLANET can be run using its GUI or on command line. GLANET's source code is available at https://github.com/burcakotlu/GLANET . Tutorials are provided at https://glanet.readthedocs.org .

Contact

burcak@ceng.metu.edu.tr or oznur.tastan@cs.bilkent.edu.tr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Denoising genome-wide histone ChIP-seq with convolutional neural networks

Motivation

Chromatin immune-precipitation sequencing (ChIP-seq) experiments are commonly used to obtain genome-wide profiles of histone modifications associated with different types of functional genomic elements. However, the quality of histone ChIP-seq data is affected by many experimental parameters such as the amount of input DNA, antibody specificity, ChIP enrichment and sequencing depth. Making accurate inferences from chromatin profiling experiments that involve diverse experimental parameters is challenging.

Results

We introduce a convolutional denoising algorithm, Coda, that uses convolutional neural networks to learn a mapping from suboptimal to high-quality histone ChIP-seq data. This overcomes various sources of noise and variability, substantially enhancing and recovering signal when applied to low-quality chromatin profiling datasets across individuals, cell types and species. Our method has the potential to improve data quality at reduced costs. More broadly, this approach-using a high-dimensional discriminative model to encode a generative noise process-is generally applicable to other biological domains where it is easy to generate noisy data but difficult to analytically characterize the noise or underlying data distribution.

Availability and implementation

https://github.com/kundajelab/coda .

Contact

akundaje@stanford.edu.

NicE-seq: high resolution open chromatin profiling

Open chromatin profiling integrates information across diverse regulatory elements to reveal the transcriptionally active genome. Tn5 transposase and DNase I sequencing-based methods prefer native or high cell numbers. Here, we describe NicE-seq (nicking enzyme assisted sequencing) for high-resolution open chromatin profiling on both native and formaldehyde-fixed cells. NicE-seq captures and reveals open chromatin sites (OCSs) and transcription factor occupancy at single nucleotide resolution, coincident with DNase hypersensitive and ATAC-seq sites at a low sequencing burden. OCSs correlate with RNA polymerase II occupancy and active chromatin marks, while displaying a contrasting pattern to CpG methylation. Decitabine-mediated hypomethylation of HCT116 displays higher numbers of OCSs.

Insights into inner ear-specific gene regulation: Epigenetics and non-coding RNAs in inner ear development and regeneration

The vertebrate inner ear houses highly specialized sensory organs, tuned to detect and encode sound, head motion and gravity. Gene expression programs under the control of transcription factors orchestrate the formation and specialization of the non-sensory inner ear labyrinth and its sensory constituents. More recently, epigenetic factors and non-coding RNAs emerged as an additional layer of gene regulation, both in inner ear development and disease. In this review, we provide an overview on how epigenetic modifications and non-coding RNAs, in particular microRNAs (miRNAs), influence gene expression and summarize recent discoveries that highlight their critical role in the proper formation of the inner ear labyrinth and its sensory organs. Finally, we discuss recent insights into how epigenetic factors and miRNAs may facilitate, or in the case of mammals, restrict inner ear sensory hair cell regeneration.

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.

RGmatch: matching genomic regions to proximal genes in omics data integration

Background

The integrative analysis of multiple genomics data often requires that genome coordinates-based signals have to be associated with proximal genes. The relative location of a genomic region with respect to the gene (gene area) is important for functional data interpretation; hence algorithms that match regions to genes should be able to deliver insight into this information.

Results

In this work we review the tools that are publicly available for making region-to-gene associations. We also present a novel method, RGmatch, a flexible and easy-to-use Python tool that computes associations either at the gene, transcript, or exon level, applying a set of rules to annotate each region-gene association with the region location within the gene. RGmatch can be applied to any organism as long as genome annotation is available. Furthermore, we qualitatively and quantitatively compare RGmatch to other tools.

Conclusions

RGmatch simplifies the association of a genomic region with its closest gene. At the same time, it is a powerful tool because the rules used to annotate these associations are very easy to modify according to the researcher’s specific interests. Some important differences between RGmatch and other similar tools already in existence are RGmatch’s flexibility, its wide range of user options, compatibility with any annotatable organism, and its comprehensive and user-friendly output.

Electronic supplementary material

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

Prenatal exposures and DNA methylation in newborns: a pilot study in Durban, South Africa

The in utero environment has the potential to influence epigenetic programming and subsequently the health of offspring. Even though pregnant women living in urban Africa are exposed to multiple chemicals and infectious agents that may impact their developing children, the neonatal epigenome has not been studied in these regions. We assessed whether prenatal exposures to air pollution and maternal human immunodeficiency virus (HIV) are associated with changes to DNA methylation throughout the epigenome using a pilot sample from the Mother and Child Environmental (MACE) birth cohort, of which 36% of the mothers are HIV positive. Families living in a high air pollution region (south Durban, n = 11) and a low air pollution region (north Durban, n = 11) with comparable socioeconomic characteristics were selected for analysis. DNA methylation was quantified in cord blood plasma DNA at >430 000 CpG sites using the Infinium HumanMethylation450 BeadChip. Sites associated with living in south Durban or maternal HIV infection (p < 0.001) were more likely to be hypomethylated and located in CpG islands. Top differentially methylated sites by region of Durban were enriched in pathways related to xenobiotic metabolism, oxygen and gas transport, and sensory perception of chemical stimuli when performing gene set enrichment testing with LRpath. Differentially methylated sites by maternal HIV status were enriched in cytochrome P450s, pathways involved in detection of chemical stimuli, metabolic processes, and viral regulation and processing. Given the small sample size of the study, future work examining the impact of prenatal exposures to air pollution, maternal infection, and antiviral treatment on the epigenome and downstream health implications is merited in Sub-Saharan African populations.

The genome-wide DNA sequence specificity of the anti-tumour drug bleomycin in human cells

The cancer chemotherapeutic agent, bleomycin, cleaves DNA at specific sites. For the first time, the genome-wide DNA sequence specificity of bleomycin breakage was determined in human cells. Utilising Illumina next-generation DNA sequencing techniques, over 200 million bleomycin cleavage sites were examined to elucidate the bleomycin genome-wide DNA selectivity. The genome-wide bleomycin cleavage data were analysed by four different methods to determine the cellular DNA sequence specificity of bleomycin strand breakage. For the most highly cleaved DNA sequences, the preferred site of bleomycin breakage was at 5'-GT* dinucleotide sequences (where the asterisk indicates the bleomycin cleavage site), with lesser cleavage at 5'-GC* dinucleotides. This investigation also determined longer bleomycin cleavage sequences, with preferred cleavage at 5'-GT*A and 5'- TGT* trinucleotide sequences, and 5'-TGT*A tetranucleotides. For cellular DNA, the hexanucleotide DNA sequence 5'-RTGT*AY (where R is a purine and Y is a pyrimidine) was the most highly cleaved DNA sequence. It was striking that alternating purine-pyrimidine sequences were highly cleaved by bleomycin. The highest intensity cleavage sites in cellular and purified DNA were very similar although there were some minor differences. Statistical nucleotide frequency analysis indicated a G nucleotide was present at the -3 position (relative to the cleavage site) in cellular DNA but was absent in purified DNA.

Massively parallel cis-regulatory analysis in the mammalian central nervous system

Cis-regulatory elements (CREs, e.g., promoters and enhancers) regulate gene expression, and variants within CREs can modulate disease risk. Next-generation sequencing has enabled the rapid generation of genomic data that predict the locations of CREs, but a bottleneck lies in functionally interpreting these data. To address this issue, massively parallel reporter assays (MPRAs) have emerged, in which barcoded reporter libraries are introduced into cells, and the resulting barcoded transcripts are quantified by next-generation sequencing. Thus far, MPRAs have been largely restricted to assaying short CREs in a limited repertoire of cultured cell types. Here, we present two advances that extend the biological relevance and applicability of MPRAs. First, we adapt exome capture technology to instead capture candidate CREs, thereby tiling across the targeted regions and markedly increasing the length of CREs that can be readily assayed. Second, we package the library into adeno-associated virus (AAV), thereby allowing delivery to target organs in vivo. As a proof of concept, we introduce a capture library of about 46,000 constructs, corresponding to roughly 3500 DNase I hypersensitive (DHS) sites, into the mouse retina by ex vivo plasmid electroporation and into the mouse cerebral cortex by in vivo AAV injection. We demonstrate tissue-specific cis-regulatory activity of DHSs and provide examples of high-resolution truncation mutation analysis for multiplex parsing of CREs. Our approach should enable massively parallel functional analysis of a wide range of CREs in any organ or species that can be infected by AAV, such as nonhuman primates and human stem cell–derived organoids.

Chromatin accessibility: a window into the genome

Transcriptional activation throughout the eukaryotic lineage has been tightly linked with disruption of nucleosome organization at promoters, enhancers, silencers, insulators and locus control regions due to transcription factor binding. Regulatory DNA thus coincides with open or accessible genomic sites of remodeled chromatin. Current chromatin accessibility assays are used to separate the genome by enzymatic or chemical means and isolate either the accessible or protected locations. The isolated DNA is then quantified using a next-generation sequencing platform. Wide application of these assays has recently focused on the identification of the instrumental epigenetic changes responsible for differential gene expression, cell proliferation, functional diversification and disease development. Here we discuss the limitations and advantages of current genome-wide chromatin accessibility assays with especial attention on experimental precautions and sequence data analysis. We conclude with our perspective on future improvements necessary for moving the field of chromatin profiling forward.

Non-targeted transcription factors motifs are a systemic component of ChIP-seq datasets

Background

The global effort to annotate the non-coding portion of the human genome relies heavily on chromatin immunoprecipitation data generated with high-throughput DNA sequencing (ChIP-seq). ChIP-seq is generally successful in detailing the segments of the genome bound by the immunoprecipitated transcription factor (TF), however almost all datasets contain genomic regions devoid of the canonical motif for the TF. It remains to be determined if these regions are related to the immunoprecipitated TF or whether, despite the use of controls, there is a portion of peaks that can be attributed to other causes.

Results

Analyses across hundreds of ChIP-seq datasets generated for sequence-specific DNA binding TFs reveal a small set of TF binding profiles for which predicted TF binding site motifs are repeatedly observed to be significantly enriched. Grouping related binding profiles, the set includes: CTCF-like, ETS-like, JUN-like, and THAP11 profiles. These frequently enriched profiles are termed ‘zingers’ to highlight their unanticipated enrichment in datasets for which they were not the targeted TF, and their potential impact on the interpretation and analysis of TF ChIP-seq data. Peaks with zinger motifs and lacking the ChIPped TF’s motif are observed to compose up to 45% of a ChIP-seq dataset. There is substantial overlap of zinger motif containing regions between diverse TF datasets, suggesting a mechanism that is not TF-specific for the recovery of these regions.

Conclusions

Based on the zinger regions proximity to cohesin-bound segments, a loading station model is proposed. Further study of zingers will advance understanding of gene regulation.

Electronic supplementary material

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

Seminoma and embryonal carcinoma footprints identified by analysis of integrated genome-wide epigenetic and expression profiles of germ cell cancer cell lines

Background

Originating from Primordial Germ Cells/gonocytes and developing via a precursor lesion called Carcinoma In Situ (CIS), Germ Cell Cancers (GCC) are the most common cancer in young men, subdivided in seminoma (SE) and non-seminoma (NS). During physiological germ cell formation/maturation, epigenetic processes guard homeostasis by regulating the accessibility of the DNA to facilitate transcription. Epigenetic deregulation through genetic and environmental parameters (i.e. genvironment) could disrupt embryonic germ cell development, resulting in delayed or blocked maturation. This potentially facilitates the formation of CIS and progression to invasive GCC. Therefore, determining the epigenetic and functional genomic landscape in GCC cell lines could provide insight into the pathophysiology and etiology of GCC and provide guidance for targeted functional experiments.

Results

This study aims at identifying epigenetic footprints in SE and EC cell lines in genome-wide profiles by studying the interaction between gene expression, DNA CpG methylation and histone modifications, and their function in the pathophysiology and etiology of GCC. Two well characterized GCC-derived cell lines were compared, one representative for SE (TCam-2) and the other for EC (NCCIT). Data were acquired using the Illumina HumanHT-12-v4 (gene expression) and HumanMethylation450 BeadChip (methylation) microarrays as well as ChIP-sequencing (activating histone modifications (H3K4me3, H3K27ac)). Results indicate known germ cell markers not only to be differentiating between SE and NS at the expression level, but also in the epigenetic landscape.

Conclusion

The overall similarity between TCam-2/NCCIT support an erased embryonic germ cell arrested in early gonadal development as common cell of origin although the exact developmental stage from which the tumor cells are derived might differ. Indeed, subtle difference in the (integrated) epigenetic and expression profiles indicate TCam-2 to exhibit a more germ cell-like profile, whereas NCCIT shows a more pluripotent phenotype. The results provide insight into the functional genome in GCC cell lines.

A systematic description of MLL fusion gene formation

Rearrangements of the MLL gene involve multiple partners and are implicated in both therapy related acute leukemia [tAL] and infant acute leukemia. For these diseases, recently compiled clinical data confirms an elevated frequency of such breakpoints within a 4 kb tract between exon 11 and a region of structural instability adjacent to exon 12. Linked primarily to cases of tAL, interference with topoisomerase II activity may either contribute to the initial DNA lesion directly or indirectly by, for example, providing a physical block to transcription progression. Alternatively, sites of fragmentation may be mis-repaired, guided by intergenic spliced transcripts of the participating genes. Co-transcription of MLL and potential fusion partners may provide the localization that enhances the probability of gene interaction. An indirect role for the leukemogenic activity of topoisomerase II inhibitors would imply that the negative consequences of their use may be separated from their therapeutic effects.

The human imprintome: regulatory mechanisms, methods of ascertainment, and roles in disease susceptibility

Imprinted genes form a special subset of the genome, exhibiting monoallelic expression in a parent-of-origin-dependent fashion. This monoallelic expression is controlled by parental-specific epigenetic marks, which are established in gametogenesis and early embryonic development and are persistent in all somatic cells throughout life. We define this specific set of cis-acting epigenetic regulatory elements as the imprintome, a distinct and specially tasked subset of the epigenome. Imprintome elements contain DNA methylation and histone modifications that regulate monoallelic expression by affecting promoter accessibility, chromatin structure, and chromatin configuration. Understanding their regulation is critical because a significant proportion of human imprinted genes are implicated in complex diseases. Significant species variation in the repertoire of imprinted genes and their epigenetic regulation, however, will not allow model organisms solely to be used for this crucial purpose. Ultimately, only the human will suffice to accurately define the human imprintome.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples.

We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype.

GEO Series accession number: GSE33784, GSE33867.

Comprehensive characterization of erythroid-specific enhancers in the genomic regions of human Krüppel-like factors

Background

Mapping of DNase I hypersensitive sites (DHSs) is a powerful tool to experimentally identify cis-regulatory elements (CREs). Among CREs, enhancers are abundant and predominantly act in driving cell-specific gene expression. Krüppel-like factors (KLFs) are a family of eukaryotic transcription factors. Several KLFs have been demonstrated to play important roles in hematopoiesis. However, transcriptional regulation of KLFs via CREs, particularly enhancers, in erythroid cells has been poorly understood.

Results

In this study, 23 erythroid-specific or putative erythroid-specific DHSs were identified by DNase-seq in the genomic regions of 17 human KLFs, and their enhancer activities were evaluated using dual-luciferase reporter (DLR) assay. Of the 23 erythroid-specific DHSs, the enhancer activities of 15 DHSs were comparable to that of the classical enhancer HS2 in driving minimal promoter (minP). Fifteen DHSs, some overlapping those that increased minP activities, acted as enhancers when driving the corresponding KLF promoters (KLF-Ps) in erythroid cells; of these, 10 DHSs were finally characterized as erythroid-specific KLF enhancers. These 10 erythroid-specific KLF enhancers were further confirmed using chromatin immunoprecipitation coupled to sequencing (ChIP-seq) data-based bioinformatic and biochemical analyses.

Conclusion

Our present findings provide a feasible strategy to extensively identify gene- and cell-specific enhancers from DHSs obtained by high-throughput sequencing, which will help reveal the transcriptional regulation and biological functions of genes in some specific cells.

Targeting H3K4 trimethylation in Huntington disease

Transcriptional dysregulation is an early feature of Huntington disease (HD). We observed gene-specific changes in histone H3 lysine 4 trimethylation (H3K4me3) at transcriptionally repressed promoters in R6/2 mouse and human HD brain. Genome-wide analysis showed a chromatin signature for this mark. Reducing the levels of the H3K4 demethylase SMCX/Jarid1c in primary neurons reversed down-regulation of key neuronal genes caused by mutant Huntingtin expression. Finally, reduction of SMCX/Jarid1c in primary neurons from BACHD mice or the single Jarid1 in a Drosophila HD model was protective. Therefore, targeting this epigenetic signature may be an effective strategy to ameliorate the consequences of HD.

New insights from existing sequence data: generating breakthroughs without a pipette

With the rapidly declining cost of data generation and the accumulation of massive data sets, molecular biology is entering an era in which incisive analysis of existing data will play an increasingly prominent role in the discovery of new biological phenomena and the elucidation of molecular mechanisms. Here, we discuss resources of publicly available sequencing data most useful for interrogating the mechanisms of gene expression. Existing next-generation sequence data sets, however, come with significant challenges in the form of technical and bioinformatic artifacts, which we discuss in detail. We also recount several breakthroughs made largely through the analysis of existing data, primarily in the RNA field.

Cell-Type Specific Determinants of NRAMP1 Expression in Professional Phagocytes

The Natural resistance-associated macrophage protein 1 (Nramp1 or Solute carrier 11 member 1, Slc11a1) transports divalent metals across the membrane of late endosomes and lysosomes in professional phagocytes. Nramp1 represents an ancient eukaryotic cell-autonomous defense whereas the gene duplication that yielded Nramp1 and Nramp2 predated the origin of Sarcopterygians (lobe-finned fishes and tetrapods). SLC11A1 genetic polymorphisms associated with human resistance to tuberculosis consist of potential regulatory variants. Herein, current knowledge of the regulation of SLC11A1 gene expression is reviewed and comprehensive analysis of ENCODE data available for hematopoietic cell-types suggests a hypothesis for the regulation of SLC11A1 expression during myeloid development and phagocyte functional polarization. SLC11A1 is part of a 34.6 kb CTCF-insulated locus scattered with predicted regulatory elements: a 3' enhancer, a large 5' enhancer domain and four elements spread around the transcription start site (TSS), including several C/EBP and PU.1 sites. SLC11A1 locus ends appear mobilized by ETS-related factors early during myelopoiesis; activation of both 5' and 3' enhancers in myelo-monocytic cells correlate with transcription factor binding at the TSS. Characterizing the corresponding cis/trans determinants functionally will establish the mechanisms involved and possibly reveal genetic variation that impacts susceptibility to infectious or immune diseases.

Long distance relationships: enhancer-promoter communication and dynamic gene transcription

The three-dimensional regulation of gene transcription involves loop formation between enhancer and promoter elements, controlling spatiotemporal gene expression in multicellular organisms. Enhancers are usually located in non-coding DNA and can activate gene transcription by recruiting transcription factors, chromatin remodeling factors and RNA Polymerase II. Research over the last few years has revealed that enhancers have tell-tale characteristics that facilitate their detection by several approaches, although the hallmarks of enhancers are not always uniform. Enhancers likely play an important role in the activation of genes by functioning as a primary point of contact for transcriptional activators, and by making physical contact with gene promoters often by means of a chromatin loop. Although numerous transcriptional regulators participate in the formation of chromatin loops that bring enhancers into proximity with promoters, the mechanism(s) of enhancer-promoter connectivity remain enigmatic. Here we discuss enhancer function, review some of the many proteins shown to be involved in establishing enhancer-promoter loops, and describe the dynamics of enhancer-promoter contacts during development, differentiation and in specific cell types.

Comparative genome-wide DNA methylation analysis of colorectal tumor and matched normal tissues

Aberrant DNA methylation often occurs in colorectal cancer (CRC). In our study we applied a genome-wide DNA methylation analysis approach, MethylCap-seq, to map the differentially methylated regions (DMRs) in 24 tumors and matched normal colon samples. In total, 2687 frequently hypermethylated and 468 frequently hypomethylated regions were identified, which include potential biomarkers for CRC diagnosis. Hypermethylation in the tumor samples was enriched at CpG islands and gene promoters, while hypomethylation was distributed throughout the genome. Using epigenetic data from human embryonic stem cells, we show that frequently hypermethylated regions coincide with bivalent loci in human embryonic stem cells. DNA methylation is commonly thought to lead to gene silencing; however, integration of publically available gene expression data indicates that 75% of the frequently hypermethylated genes were most likely already lowly or not expressed in normal tissue. Collectively, our study provides genome-wide DNA methylation maps of CRC, comprehensive lists of DMRs, and gives insights into the role of aberrant DNA methylation in CRC formation.

Comprehensive identification and annotation of cell type-specific and ubiquitous CTCF-binding sites in the human genome

Chromatin insulators are DNA elements that regulate the level of gene expression either by preventing gene silencing through the maintenance of heterochromatin boundaries or by preventing gene activation by blocking interactions between enhancers and promoters. CCCTC-binding factor (CTCF), a ubiquitously expressed 11-zinc-finger DNA-binding protein, is the only protein implicated in the establishment of insulators in vertebrates. While CTCF has been implicated in diverse regulatory functions, CTCF has only been studied in a limited number of cell types across human genome. Thus, it is not clear whether the identified cell type-specific differences in CTCF-binding sites are functionally significant. Here, we identify and characterize cell type-specific and ubiquitous CTCF-binding sites in the human genome across 38 cell types designated by the Encyclopedia of DNA Elements (ENCODE) consortium. These cell type-specific and ubiquitous CTCF-binding sites show uniquely versatile transcriptional functions and characteristic chromatin features. In addition, we confirm the insulator barrier function of CTCF-binding and explore the novel function of CTCF in DNA replication. These results represent a critical step toward the comprehensive and systematic understanding of CTCF-dependent insulators and their versatile roles in the human genome.

Intra-nuclear mobility and target search mechanisms of transcription factors: a single-molecule perspective on gene expression

Precise expression of specific genes in time and space is at the basis of cellular viability as well as correct development of organisms. Understanding the mechanisms of gene regulation is fundamental and still one of the great challenges for biology. Gene expression is regulated also by specific transcription factors that recognize and bind to specific DNA sequences. Transcription factors dynamics, and especially the way they sample the nucleoplasmic space during the search for their specific target in the genome, are a key aspect for regulation and it has been puzzling researchers for forty years. The scope of this review is to give a state-of-the-art perspective over the intra-nuclear mobility and the target search mechanisms of specific transcription factors at the molecular level. Going through the seminal biochemical experiments that have raised the first questions about target localization and the theoretical grounds concerning target search processes, we describe the most recent experimental achievements and current challenges in understanding transcription factors dynamics and interactions with DNA using in vitro assays as well as in live prokaryotic and eukaryotic cells. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Identification of CTCF as a master regulator of the clustered protocadherin genes

The brain is a large and complex network of neurons. Specific neuronal connectivity is thought to be based on the combinatorial expression of the 52 protocadherins (Pcdh) membrane adhesion proteins, whereby each neuron expresses only a specific subset. Pcdh genes are arranged in tandem, in a cluster of three families: Pcdhα, Pcdhβ and Pcdhγ. The expression of each Pcdh gene is regulated by a promoter that has a regulatory conserved sequence element (CSE), common to all 52 genes. The mechanism and factors controlling individual Pcdh gene expression are currently unknown. Here we show that the promoter of each Pcdh gene contains a gene-specific conserved control region, termed specific sequence element (SSE), located adjacent and upstream to the CSE and activates transcription together with the CSE. We purified the complex that specifically binds the SSE-CSE region and identified the CCTC binding-factor (CTCF) as a key molecule that binds and activates Pcdh promoters. Our findings point to CTCF as a factor essential for Pcdh expression and probably governing neuronal connectivity.