Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays

A data-adaptive methods in detecting exogenous methyltransferase accessible chromatin in human genome using nanopore sequencing

MOTIVATION

Identifying chromatin accessibility is one of the key steps in studying the regulation of eukaryotic genomes. The combination of exogenous methyltransferase and nanopore sequencing provides an strategy to identify open chromatin over long genomic ranges at the single-molecule scale. However, endogenous methylation, non-open-chromatin-specific exogenous methylation and base-calling errors limit the accuracy and hinders its application to complex genomes.

RESULTS

We systematically evaluated the impact of these three influence factors, and developed a model-based computational method, methyltransferase accessible genome region finder (MAGNIFIER), to address the issues. By incorporating control data, MAGNIFIER attenuates the three influence factors with data-adaptive comparison strategy. We demonstrate that MAGNIFIER is not only sensitive to identify the open chromatin with much improved accuracy, but also able to detect the chromatin accessibility of repetitive regions that are missed by NGS-based methods. By incorporating long-read RNA-seq data, we revealed the association between the accessible Alu elements and non-classic gene isoforms.

AVAILABILITY AND IMPLEMENTATION

Freely available on web at https://github.com/Goatofmountain/MAGNIFIER.

SEAMoD: A fully interpretable neural network for cis-regulatory analysis of differentially expressed genes

A common way to investigate gene regulatory mechanisms is to identify differentially expressed genes using transcriptomics, find their candidate enhancers using epigenomics, and search for over-represented transcription factor (TF) motifs in these enhancers using bioinformatics tools. A related follow-up task is to model gene expression as a function of enhancer sequences and rank TF motifs by their contribution to such models, thus prioritizing among regulators. We present a new computational tool called SEAMoD that performs the above tasks of motif finding and sequence-to-expression modeling simultaneously. It trains a convolutional neural network model to relate enhancer sequences to differential expression in one or more biological conditions. The model uses TF motifs to interpret the sequences, learning these motifs and their relative importance to each biological condition from data. It also utilizes epigenomic information in the form of activity scores of putative enhancers and automatically searches for the most promising enhancer for each gene. Compared to existing neural network models of non-coding sequences, SEAMoD uses far fewer parameters, requires far less training data, and emphasizes biological interpretability. We used SEAMoD to understand regulatory mechanisms underlying the differentiation of neural stem cell (NSC) derived from mouse forebrain. We profiled gene expression and histone modifications in NSC and three differentiated cell types and used SEAMoD to model differential expression of nearly 12,000 genes with an accuracy of 81%, in the process identifying the Olig2, E2f family TFs, Foxo3, and Tcf4 as key transcriptional regulators of the differentiation process.

Genome-wide mapping of DNase I hypersensitive sites in pineapple leaves

Pineapple [Ananas comosus (L.) Merr.] is the most economically important crop possessing crassulacean acid metabolism (CAM) photosynthesis which has a higher water use efficiency by control of nocturnal opening and diurnal closure of stomata. To provide novel insights into the diel regulatory landscape in pineapple leaves, we performed genome-wide mapping of DNase I hypersensitive sites (DHSs) in pineapple leaves at day (2a.m.) and night (10a.m.) using a simplified DNase-seq method. As a result, totally 33340 and 28753 DHSs were found in green-tip tissue, and 29597 and 40068 were identified in white-base tissue at 2a.m. and 10a.m., respectively. We observed that majority of the pineapple genes occupied less than two DHSs with length shorter than 1 kb, and the promotor DHSs showed a proximal trend to the transcription start site (>77% promotor DHSs within 1 kb). In addition, more intergenic DHSs were identified around transcription factors or transcription co-regulators (TFs/TCs) than other functional genes, indicating complex regulatory contexts around TFs/TCs. Through combined analysis of tissue preferential DHSs and genes, we respectively found 839 and 888 coordinately changed genes in green-tip at 2a.m. and 10a.m. (AcG2 and AcG10). Furthermore, AcG2-specific, AcG10-specific and common accessible DHSs were dissected from the total photosynthetic preferential DHSs, and the regulatory networks indicated dynamic regulations with multiple cis-regulatory elements occurred to genes preferentially expressed in photosynthetic tissues. Interestingly, binding motifs of several cycling TFs were identified in the DHSs of key CAM genes, revealing a circadian regulation to CAM coordinately diurnal expression. Our results provide a chromatin regulatory landscape in pineapple leaves during the day and night. This will provide important information to assist with deciphering the circadian regulation of CAM photosynthesis.

3D genome organization and epigenetic regulation in autoimmune diseases

Three-dimensional (3D) genomics is an emerging field of research that investigates the relationship between gene regulatory function and the spatial structure of chromatin. Chromatin folding can be studied using chromosome conformation capture (3C) technology and 3C-based derivative sequencing technologies, including chromosome conformation capture-on-chip (4C), chromosome conformation capture carbon copy (5C), and high-throughput chromosome conformation capture (Hi-C), which allow scientists to capture 3D conformations from a single site to the entire genome. A comprehensive analysis of the relationships between various regulatory components and gene function also requires the integration of multi-omics data such as genomics, transcriptomics, and epigenomics. 3D genome folding is involved in immune cell differentiation, activation, and dysfunction and participates in a wide range of diseases, including autoimmune diseases. We describe hierarchical 3D chromatin organization in this review and conclude with characteristics of C-techniques and multi-omics applications of the 3D genome. In addition, we describe the relationship between 3D genome structure and the differentiation and maturation of immune cells and address how changes in chromosome folding contribute to autoimmune diseases.

Defining Candidate Imprinted loci in Bos taurus

Using a whole-genome assembly of Bos taurus, I applied my bioinformatics strategy to locate candidate imprinting control regions (ICRs) genome-wide. In mammals, genomic imprinting plays essential roles in embryogenesis. In my strategy, peaks in plots mark the locations of known, inferred, and candidate ICRs. Genes in the vicinity of candidate ICRs correspond to potential imprinted genes. By displaying my datasets on the UCSC genome browser, one could view peak positions with respect to genomic landmarks. I give two examples of candidate ICRs in loci that influence spermatogenesis in bulls: CNNM1 and CNR1. I also give examples of candidate ICRs in loci that influence muscle development: SIX1 and BCL6. By examining the ENCODE data reported for mice, I deduced regulatory clues about cattle. I focused on DNase I hypersensitive sites (DHSs). Such sites reveal accessibility of chromatin to regulators of gene expression. For inspection, I chose DHSs in chromatin from mouse embryonic stem cells (ESCs) ES-E14, mesoderm, brain, heart, and skeletal muscle. The ENCODE data revealed that the SIX1 promoter was accessible to the transcription initiation apparatus in mouse ESCs, mesoderm, and skeletal muscles. The data also revealed accessibility of BCL6 locus to regulatory proteins in mouse ESCs and examined tissues.

Genome-Wide Mapping of Active Regulatory Elements Using ATAC-seq

Active cis-regulatory elements (cREs) in eukaryotes are characterized by nucleosomal depletion and, accordingly, higher accessibility. This property has turned out to be immensely useful for identifying cREs genome-wide and tracking their dynamics across different cellular states and is the basis of numerous methods taking advantage of the preferential enzymatic cleavage/labeling of accessible DNA. ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) has emerged as the most versatile and widely adaptable method and has been widely adopted as the standard tool for mapping open chromatin regions. Here, we discuss the current optimal practices and important considerations for carrying out ATAC-seq experiments, primarily in the context of mammalian systems.

Multi-modular structure of the gene regulatory network for specification and commitment of murine T cells

T cells develop from multipotent progenitors by a gradual process dependent on intrathymic Notch signaling and coupled with extensive proliferation. The stages leading them to T-cell lineage commitment are well characterized by single-cell and bulk RNA analyses of sorted populations and by direct measurements of precursor-product relationships. This process depends not only on Notch signaling but also on multiple transcription factors, some associated with stemness and multipotency, some with alternative lineages, and others associated with T-cell fate. These factors interact in opposing or semi-independent T cell gene regulatory network (GRN) subcircuits that are increasingly well defined. A newly comprehensive picture of this network has emerged. Importantly, because key factors in the GRN can bind to markedly different genomic sites at one stage than they do at other stages, the genes they significantly regulate are also stage-specific. Global transcriptome analyses of perturbations have revealed an underlying modular structure to the T-cell commitment GRN, separating decisions to lose "stem-ness" from decisions to block alternative fates. Finally, the updated network sheds light on the intimate relationship between the T-cell program, which depends on the thymus, and the innate lymphoid cell (ILC) program, which does not.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Advances in Single-Cell Multi-Omics and Application in Cardiovascular Research

With the development of ever more powerful and versatile high-throughput sequencing techniques and innovative ways to capture single cells, mapping the multicellular tissues at the single-cell level is becoming routine practice. However, it is still challenging to depict the epigenetic landscape of a single cell, especially the genome-wide chromatin accessibility, histone modifications, and DNA methylation. We summarize the most recent methodologies to profile these epigenetic marks at the single-cell level. We also discuss the development and advancement of several multi-omics sequencing technologies from individual cells. Advantages and limitations of various methods to compare and integrate datasets obtained from different sources are also included with specific practical notes. Understanding the heart tissue at single-cell resolution and multi-modal levels will help to elucidate the cell types and states involved in physiological and pathological events during heart development and disease. The rich information produced from single-cell multi-omics studies will also promote the research of heart regeneration and precision medicine on heart diseases.

GFI1-Dependent Repression of SGPP1 Increases Multiple Myeloma Cell Survival

Simple Summary

New therapies have greatly improved the progression-free and overall survival for patients with “standard risk” multiple myeloma (MM). However, patients with “high risk” MM, in particular patients whose MM cells harbor non-functional p53, have very short survival times because of the early relapse and rapid development of highly therapy-resistant MM. In this report, we identify a novel mechanism responsible for Growth Factor Independence-1 (GFI1) regulation of the growth and survival of MM cells through its modulation of sphingolipid metabolism, regardless of their p53 status. We identify the Sphingosine-1-Phosphate Phosphatase (SGPP1) gene as a novel direct target of GFI1 transcriptional repression in MM cells, thus increasing intracellular sphingosine-1-phosphate levels, which stabilizes c-Myc. Our results support GFI1 as an attractive therapeutic target for all types of MM, including the “high risk” patient population with non-functional p53, as well as a possible therapeutic approach for other types of cancers expressing high levels of c-Myc.

Abstract

Multiple myeloma (MM) remains incurable for most patients due to the emergence of drug resistant clones. Here we report a p53-independent mechanism responsible for Growth Factor Independence-1 (GFI1) support of MM cell survival by its modulation of sphingolipid metabolism to increase the sphingosine-1-phosphate (S1P) level regardless of the p53 status. We found that expression of enzymes that control S1P biosynthesis, SphK1, dephosphorylation, and SGPP1 were differentially correlated with GFI1 levels in MM cells. We detected GFI1 occupancy on the SGGP1 gene in MM cells in a predicted enhancer region at the 5’ end of intron 1, which correlated with decreased SGGP1 expression and increased S1P levels in GFI1 overexpressing cells, regardless of their p53 status. The high S1P:Ceramide intracellular ratio in MM cells protected c-Myc protein stability in a PP2A-dependent manner. The decreased MM viability by SphK1 inhibition was dependent on the induction of autophagy in both p53WT and p53mut MM. An autophagic blockade prevented GFI1 support for viability only in p53mut MM, demonstrating that GFI1 increases MM cell survival via both p53WT inhibition and upregulation of S1P independently. Therefore, GFI1 may be a key therapeutic target for all types of MM that may significantly benefit patients that are highly resistant to current therapies.

Single-Molecule Multikilobase-Scale Profiling of Chromatin Accessibility Using m6A-SMAC-Seq and m6A-CpG-GpC-SMAC-Seq

A hallmark feature of active cis-regulatory elements (CREs) in eukaryotes is their nucleosomal depletion and, accordingly, higher accessibility to enzymatic treatment. This property has been the basis of a number of sequencing-based assays for genome-wide identification and tracking the activity of CREs across different biological conditions, such as DNAse-seq, ATAC-seq , NOMeseq, and others. However, the fragmentation of DNA inherent to many of these assays and the limited read length of short-read sequencing platforms have so far not allowed the simultaneous measurement of the chromatin accessibility state of CREs located distally from each other. The combination of labeling accessible DNA with DNA modifications and nanopore sequencing has made it possible to develop such assays. Here, we provide a detailed protocol for carrying out the SMAC-seq assay (Single-Molecule long-read Accessible Chromatin mapping sequencing), in its m6A-SMAC-seq and m6A-CpG-GpC-SMAC-seq variants, together with methods for data processing and analysis, and discuss key experimental and analytical considerations for working with SMAC-seq datasets.

Redefining the PTEN promoter: Identification of novel upstream transcription start regions

Germline mutation of PTEN is causally observed in Cowden syndrome (CS) and is one of the most common, penetrant risk genes for autism spectrum disorder (ASD). However, the majority of individuals who present with CS-like clinical features are PTEN-mutation negative. Reassessment of PTEN promoter regulation may help explain abnormal PTEN dosage, as only the minimal promoter and coding regions are currently included in diagnostic PTEN mutation analysis. Therefore, we reanalyzed the architecture of the PTEN promoter using next-generation sequencing datasets. Specifically, run-on sequencing assays identified two additional transcription start regions (TSRs) at -2053 and - 1906 basepairs from the canonical start of PTEN, thus extending the PTEN 5'UTR and redefining the PTEN promoter. We show that these novel upstream TSRs are active in cancer cell lines, human cancer, and normal tissue. Further, these TSRs can produce novel PTEN transcripts due to the introduction of new splice donors at -2041, -1826, and - 1355, which may allow for splicing out of the PTEN 5'UTR or the first and second exon in upstream-initiated transcripts. Combining ENCODE ChIP-seq and pertinent literature, we also compile and analyze all transcription factors (TFs) binding at the redefined PTEN locus. Enrichment analyses suggest that TFs bind specifically to the upstream TSRs may be implicated in inflammatory processes. Together, these data redefine the architecture of the PTEN promoter, an important step toward a comprehensive model of PTEN transcription regulation, a basis for future investigations into the new promoters' role in disease pathogenesis.

Genome-Wide Identification of DNase I Hypersensitive Sites in Plants

DNase I hypersensitive site (DHS) mapping combined with high-throughput sequencing (DNase-seq) enables the identification of cis-regulatory DNA elements (CREs) genome wide. However, despite the wide applications of DNase-seq in plants, its application to the highly repetitive genomes of plants has lagged. Here, we describe a modified DNase-seq method, making it more practical for application to plants with genomes enriched with repetitive DNA. This approach adopts a double-hit-based strategy, in which small (<250-bp) DNA fragments digested by DNase I are selected and used for sequencing library construction. Using these protocols, we have conducted DNase-seq in plants with high content of repetitive DNA, including maize, sugarcane, and tetraploid cotton. Genome-wide maps of DHS and CREs have been created using these DNase-seq datasets. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Nuclei isolation Basic Protocol 2: DNase I digestion Basic Protocol 3: Target DNA isolation Basic Protocol 4: Library construction and validation.

Genomic Space of MGMT in Human Glioma Revisited: Novel Motifs, Regulatory RNAs, NRF1, 2, and CTCF Involvement in Gene Expression

Background: The molecular regulation of increased MGMT expression in human brain tumors, the associated regulatory elements, and linkages of these to its epigenetic silencing are not understood. Because the heightened expression or non-expression of MGMT plays a pivotal role in glioma therapeutics, we applied bioinformatics and experimental tools to identify the regulatory elements in the MGMT and neighboring EBF3 gene loci. Results: Extensive genome database analyses showed that the MGMT genomic space was rich in and harbored many undescribed RNA regulatory sequences and recognition motifs. We extended the MGMT’s exon-1 promoter to 2019 bp to include five overlapping alternate promoters. Consensus sequences in the revised promoter for (a) the transcriptional factors CTCF, NRF1/NRF2, GAF, (b) the genetic switch MYC/MAX/MAD, and (c) two well-defined p53 response elements in MGMT intron-1, were identified. A putative protein-coding or non-coding RNA sequence was located in the extended 3′ UTR of the MGMT transcript. Eleven non-coding RNA loci coding for miRNAs, antisense RNA, and lncRNAs were identified in the MGMT-EBF3 region and six of these showed validated potential for curtailing the expression of both MGMT and EBF3 genes. ChIP analysis verified the binding site in MGMT promoter for CTCF which regulates the genomic methylation and chromatin looping. CTCF depletion by a pool of specific siRNA and shRNAs led to a significant attenuation of MGMT expression in human GBM cell lines. Computational analysis of the ChIP sequence data in ENCODE showed the presence of NRF1 in the MGMT promoter and this occurred only in MGMT-proficient cell lines. Further, an enforced NRF2 expression markedly augmented the MGMT mRNA and protein levels in glioma cells. Conclusions: We provide the first evidence for several new regulatory components in the MGMT gene locus which predict complex transcriptional and posttranscriptional controls with potential for new therapeutic avenues.

Chromatin Accessibility Is Associated with Artemisinin Biosynthesis Regulation in Artemisia annua

Glandular trichome (GT) is the dominant site for artemisinin production in Artemisia annua. Several critical genes involved in artemisinin biosynthesis are specifically expressed in GT. However, the molecular mechanism of differential gene expression between GT and other tissue types remains elusive. Chromatin accessibility, defined as the degree to which nuclear molecules are able to interact with chromatin DNA, reflects gene expression capacity to a certain extent. Here, we investigated and compared the landscape of chromatin accessibility in Artemisia annua leaf and GT using the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) technique. We identified 5413 GT high accessible and 4045 GT low accessible regions, and these GT high accessible regions may contribute to GT-specific biological functions. Several GT-specific artemisinin biosynthetic genes, such as DBR2 and CYP71AV1, showed higher accessible regions in GT compared to that in leaf, implying that they might be regulated by chromatin accessibility. In addition, transcription factor binding motifs for MYB, bZIP, C2H2, and AP2 were overrepresented in the highly accessible chromatin regions associated with artemisinin biosynthetic genes in glandular trichomes. Finally, we proposed a working model illustrating the chromatin accessibility dynamics in regulating artemisinin biosynthetic gene expression. This work provided new insights into epigenetic regulation of gene expression in GT.

Interrogating the Accessible Chromatin Landscape of Eukaryote Genomes Using ATAC-seq

The ATAC-seq assay has emerged as the most useful, versatile, and widely adaptable method for profiling accessible chromatin regions and tracking the activity of cis-regulatory elements (cREs) in eukaryotes. Thanks to its great utility, it is now being applied to map active chromatin in the context of a very wide diversity of biological systems and questions. In the course of these studies, considerable experience working with ATAC-seq data has accumulated and a standard set of computational tasks that need to be carried for most ATAC-seq analyses has emerged. Here, we review and provide examples of common such analytical procedures (including data processing, quality control, peak calling, identifying differentially accessible open chromatin regions, and variable transcription factor (TF) motif accessibility) and discuss recommended optimal practices.

Lineage-Resolved Enhancer and Promoter Usage during a Time Course of Embryogenesis

Enhancers are essential drivers of cell states, yet the relationship between accessibility, regulatory activity, and in vivo lineage commitment during embryogenesis remains poorly understood. Here, we measure chromatin accessibility in isolated neural and mesodermal lineages across a time course of Drosophila embryogenesis. Promoters, including tissue-specific genes, are often constitutively open, even in contexts where the gene is not expressed. In contrast, the majority of distal elements have dynamic, tissue-specific accessibility. Enhancer priming appears rarely within a lineage, perhaps reflecting the speed of Drosophila embryogenesis. However, many tissue-specific enhancers are accessible in other lineages early on and become progressively closed as embryogenesis proceeds. We demonstrate the usefulness of this tissue- and time-resolved resource to definitively identify single-cell clusters, to uncover predictive motifs, and to identify many regulators of tissue development. For one such predicted neural regulator, l(3)neo38, we generate a loss-of-function mutant and uncover an essential role for neuromuscular junction and brain development.

Enhancer Predictions and Genome-Wide Regulatory Circuits

Spatiotemporal control of gene expression during development requires orchestrated activities of numerous enhancers, which are cis-regulatory DNA sequences that, when bound by transcription factors, support selective activation or repression of associated genes. Proper activation of enhancers is critical during embryonic development, adult tissue homeostasis, and regeneration, and inappropriate enhancer activity is often associated with pathological conditions such as cancer. Multiple consortia [e.g., the Encyclopedia of DNA Elements (ENCODE) Consortium and National Institutes of Health Roadmap Epigenomics Mapping Consortium] and independent investigators have mapped putative regulatory regions in a large number of cell types and tissues, but the sequence determinants of cell-specific enhancers are not yet fully understood. Machine learning approaches trained on large sets of these regulatory regions can identify core transcription factor binding sites and generate quantitative predictions of enhancer activity and the impact of sequence variants on activity. Here, we review these computational methods in the context of enhancer prediction and gene regulatory network models specifying cell fate.

A novel NGS library preparation method to characterize native termini of fragmented DNA

Biological and chemical DNA fragmentation generates DNA molecules with a variety of termini, including blunt ends and single-stranded overhangs. We have developed a Next Generation Sequencing (NGS) assay, XACTLY, to interrogate the termini of fragmented DNA, information traditionally lost in standard NGS library preparation methods. Here we describe the XACTLY method, showcase its sensitivity and specificity, and demonstrate its utility in in vitro experiments. The XACTLY assay is able to report relative abundances of all lengths and types (5′ and 3′) of single-stranded overhangs, if present, on each DNA fragment with an overall accuracy between 80–90%. In addition, XACTLY retains the sequence of each native DNA molecule after fragmentation and can capture the genomic landscape of cleavage events at single nucleotide resolution. The XACTLY assay can be applied as a novel research and discovery tool for fragmentation analyses and in cell-free DNA.

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.

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

Background

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

Results

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

Conclusions

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

Towards a comprehensive catalogue of validated and target-linked human enhancers

The human gene catalogue is essentially complete, but we lack an equivalently vetted inventory of bona fide human enhancers. Hundreds of thousands of candidate enhancers have been nominated via biochemical annotations; however, only a handful of these have been validated and confidently linked to their target genes. Here we review emerging technologies for discovering, characterizing and validating human enhancers at scale. We furthermore propose a new framework for operationally defining enhancers that accommodates the heterogeneous and complementary results that are emerging from reporter assays, biochemical measurements and CRISPR screens.

Long-range single-molecule mapping of chromatin accessibility in eukaryotes

Mapping open chromatin regions has emerged as a widely used tool for identifying active regulatory elements in eukaryotes. However, existing approaches, limited by reliance on DNA fragmentation and short-read sequencing, cannot provide information about large-scale chromatin states or reveal coordination between the states of distal regulatory elements. We have developed a method for profiling the accessibility of individual chromatin fibers, a single-molecule long-read accessible chromatin mapping sequencing assay (SMAC-seq), enabling the simultaneous, high-resolution, single-molecule assessment of chromatin states at multikilobase length scales. Our strategy is based on combining the preferential methylation of open chromatin regions by DNA methyltransferases with low sequence specificity, in this case EcoGII, an N6-methyladenosine (m6A) methyltransferase, and the ability of nanopore sequencing to directly read DNA modifications. We demonstrate that aggregate SMAC-seq signals match bulk-level accessibility measurements, observe single-molecule nucleosome and transcription factor protection footprints, and quantify the correlation between chromatin states of distal genomic elements.

A framework for transcriptome-wide association studies in breast cancer in diverse study populations

BACKGROUND

The relationship between germline genetic variation and breast cancer survival is largely unknown, especially in understudied minority populations who often have poorer survival. Genome-wide association studies (GWAS) have interrogated breast cancer survival but often are underpowered due to subtype heterogeneity and clinical covariates and detect loci in non-coding regions that are difficult to interpret. Transcriptome-wide association studies (TWAS) show increased power in detecting functionally relevant loci by leveraging expression quantitative trait loci (eQTLs) from external reference panels in relevant tissues. However, ancestry- or race-specific reference panels may be needed to draw correct inference in ancestrally diverse cohorts. Such panels for breast cancer are lacking.

RESULTS

We provide a framework for TWAS for breast cancer in diverse populations, using data from the Carolina Breast Cancer Study (CBCS), a population-based cohort that oversampled black women. We perform eQTL analysis for 406 breast cancer-related genes to train race-stratified predictive models of tumor expression from germline genotypes. Using these models, we impute expression in independent data from CBCS and TCGA, accounting for sampling variability in assessing performance. These models are not applicable across race, and their predictive performance varies across tumor subtype. Within CBCS (N = 3,828), at a false discovery-adjusted significance of 0.10 and stratifying for race, we identify associations in black women near AURKA, CAPN13, PIK3CA, and SERPINB5 via TWAS that are underpowered in GWAS.

CONCLUSIONS

We show that carefully implemented and thoroughly validated TWAS is an efficient approach for understanding the genetics underpinning breast cancer outcomes in diverse populations.

CATA: a comprehensive chromatin accessibility database for cancer

Accessible chromatin refers to the active regions of a chromosome that are bound by many transcription factors (TFs). Changes in chromatin accessibility play a critical role in tumorigenesis. With the emergence of novel methods like Assay for Transposase-accessible Chromatin Sequencing, a sequencing method that maps chromatin-accessible regions (CARs) and enables the computational analysis of TF binding at chromatin-accessible sites, the regulatory landscape in cancer can be dissected. Herein, we developed a comprehensive cancer chromatin accessibility database named CATA, which aims to provide available resources of cancer CARs and to annotate their potential roles in the regulation of genes in a cancer type-specific manner. In this version, CATA stores 2 991 163 CARs from 23 cancer types, binding information of 1398 TFs within the CARs, and provides multiple annotations about these regions, including common single nucleotide polymorphisms (SNPs), risk SNPs, copy number variation, somatic mutations, motif changes, expression quantitative trait loci, methylation and CRISPR/Cas9 target loci. Moreover, CATA supports cancer survival analysis of the CAR-associated genes and provides detailed clinical information of the tumor samples. Database URL: CATA is available at http://www.xiejjlab.bio/cata/.

Genome-Wide Identification of Regulatory DNA Elements in Crop Plants

cis-regulatory DNA elements (CREs) are noncoding but functional DNA sequences. The binding of regulatory proteins into CRE regions leads to chromatin high sensitive to DNase I digestion, which are termed as DNase I hypersensitive sites (DHSs). These DHSs can be efficiently detected through DNase I digestion followed by high-throughput DNA sequencing (DNase-seq). Thus, DNase-seq has become a powerful technique for DHSs mapping at whole-genome level in both plants and animals. Here we describe a DNase-seq procedure modified and developed for crop plants. These plants usually contain large amounts of repetitive sequences and complex organic constituents. With the main improvement in nuclei isolation, this method has been successfully used in mapping DHSs in cotton and sugarcane.

Systematic integration of GATA transcription factors and epigenomes via IDEAS paints the regulatory landscape of hematopoietic cells

Members of the GATA family of transcription factors play key roles in the differentiation of specific cell lineages by regulating the expression of target genes. Three GATA factors play distinct roles in hematopoietic differentiation. In order to better understand how these GATA factors function to regulate genes throughout the genome, we are studying the epigenomic and transcriptional landscapes of hematopoietic cells in a model-driven, integrative fashion. We have formed the collaborative multi-lab VISION project to conduct ValIdated Systematic IntegratiON of epigenomic data in mouse and human hematopoiesis. The epigenomic data included nuclease accessibility in chromatin, CTCF occupancy, and histone H3 modifications for 20 cell types covering hematopoietic stem cells, multilineage progenitor cells, and mature cells across the blood cell lineages of mouse. The analysis used the Integrative and Discriminative Epigenome Annotation System (IDEAS), which learns all common combinations of features (epigenetic states) simultaneously in two dimensions-along chromosomes and across cell types. The result is a segmentation that effectively paints the regulatory landscape in readily interpretable views, revealing constitutively active or silent loci as well as the loci specifically induced or repressed in each stage and lineage. Nuclease accessible DNA segments in active chromatin states were designated candidate cis-regulatory elements in each cell type, providing one of the most comprehensive registries of candidate hematopoietic regulatory elements to date. Applications of VISION resources are illustrated for the regulation of genes encoding GATA1, GATA2, GATA3, and Ikaros. VISION resources are freely available from our website http://usevision.org.

Local Changes in Chromatin Accessibility and Transcriptional Networks Underlying the Nitrate Response in Arabidopsis Roots

Transcriptional regulation, determined by the chromatin structure and regulatory elements interacting at promoter regions, is a key step in plant responses to environmental cues. Nitrate (NO₃^-) is a nutrient signal that regulates the expression of hundreds of genes in Arabidopsis thaliana. Here, we integrate mRNA sequencing, genome-wide RNA polymerase II (RNPII), chromatin immunoprecipitation sequencing, and DNase sequencing datasets to establish the relationship between RNPII occupancy and chromatin accessibility in response to NO₃^- treatments in Arabidopsis roots. Genomic footprinting allowed us to identify in vivo regulatory elements controlling gene expression in response to NO₃^- treatments. NO₃^--modulated transcription factor (TF) footprints are important for a rapid increase in RNPII occupancy and transcript accumulation over time. We mapped key TF regulatory interactions and functionally validated the role of NAP, an NAC-domain containing TF, as a new regulatory factor in NO₃^- transport. Taken together, our study provides a comprehensive view of transcriptional networks in response to a nutrient signal in Arabidopsis roots.

Genomic methods in profiling DNA accessibility and factor localization

Recent advancements in next-generation sequencing technologies and accompanying reductions in cost have led to an explosion of techniques to examine DNA accessibility and protein localization on chromatin genome-wide. Generally, accessible regions of chromatin are permissive for factor binding and are therefore hotspots for regulation of gene expression; conversely, genomic regions that are highly occupied by histone proteins are not permissive for factor binding and are less likely to be active regulatory regions. Identifying regions of differential accessibility can be useful to uncover putative gene regulatory regions, such as enhancers, promoters, and insulators. In addition, DNA-binding proteins, such as transcription factors that preferentially bind certain DNA sequences and histone proteins that form the core of the nucleosome, play essential roles in all DNA-templated processes. Determining the genomic localization of chromatin-bound proteins is therefore essential in determining functional roles, sequence motifs important for factor binding, and regulatory networks controlling gene expression. In this review, we discuss techniques for determining DNA accessibility and nucleosome positioning (DNase-seq, FAIRE-seq, MNase-seq, and ATAC-seq) and techniques for detecting and functionally characterizing chromatin-bound proteins (ChIP-seq, DamID, and CUT&RUN). These methods have been optimized to varying degrees of resolution, specificity, and ease of use. Here, we outline some advantages and disadvantages of these techniques, their general protocols, and a brief discussion of their development. Together, these complimentary approaches have provided an unparalleled view of chromatin architecture and functional gene regulation.

Nucleic Acid Platform Technologies

Every microarray experiment is based on a common format. First, a large number of nucleotide "spots" are arrayed onto a substrate, typically a glass slide, a silicon chip, or microbeads. Second, a complex population of nucleic acids (isolated from cells, selected from in vitro-synthesized libraries, or obtained from another source) is labeled, typically with fluorescent dyes. Third, the labeled nucleic acids are allowed to hybridize to their complementary spot(s) on the microarray. Fourth, the hybridized microarray is washed, allowing the amount of hybridized label to then be quantified. Analysis of the raw data generates a readout of the levels of each species of RNA in the original complex population. This introduction includes several examples of microarray applications and provides a discussion of the basic steps of most microarray experiments.

Chromatin accessibility and the regulatory epigenome

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.

HNF4α and CDX2 Regulate Intestinal YAP1 Promoter Activity

The Hippo pathway is important for tissue homeostasis, regulation of organ size and growth in most tissues. The co-transcription factor yes-associated protein 1 (YAP1) serves as a main downstream effector of the Hippo pathway and its dysregulation increases cancer development and blocks colonic tissue repair. Nevertheless, little is known about the transcriptional regulation of YAP1 in intestinal cells. The aim of this study to identify gene control regions in the YAP1 gene and transcription factors important for intestinal expression. Bioinformatic analysis of caudal type homeobox 2 (CDX2) and hepatocyte nuclear factor 4 alpha (HNF4α) chromatin immunoprecipitated DNA from differentiated Caco-2 cells revealed potential intragenic enhancers in the YAP1 gene. Transfection of luciferase-expressing YAP1 promoter-reporter constructs containing the potential enhancer regions validated one potent enhancer of the YAP1 promoter activity in Caco-2 and T84 cells. Two potential CDX2 and one HNF4α binding sites were identified in the enhancer by in silico transcription factor binding site analysis and protein-DNA binding was confirmed in vitro using electrophoretic mobility shift assay. It was found by chromatin immunoprecipitation experiments that CDX2 and HNF4α bind to the YAP1 enhancer in Caco-2 cells. These results reveal a previously unknown enhancer of the YAP1 promoter activity in the YAP1 gene, with importance for high expression levels in intestinal epithelial cells. Additionally, CDX2 and HNF4α binding are important for the YAP1 enhancer activity in intestinal epithelial cells.

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.

Evolution and Comprehensive Analysis of DNaseI Hypersensitive Sites in Regulatory Regions of Primate Brain-Related Genes

How the human brain differs from those of non-human primates is largely unknown and the complex drivers underlying such differences at the genomic level remain unclear. In this study, we selected 243 brain-related genes, based on Gene Ontology, and identified 184,113 DNaseI hypersensitive sites (DHSs) within their regulatory regions. To performed comprehensive evolutionary analyses, we set strict filtering criteria for alignment quality and filtered 39,132 DHSs for inclusion in the investigation and found that 2,397 (~6%) exhibited evidence of accelerated evolution (aceDHSs), which was a much higher proportion that DHSs genome-wide. Target genes predicted to be regulated by brain-aceDHSs were functionally enriched for brain development and exhibited differential expression between human and chimpanzee. Alignments indicated 61 potential human-specific transcription factor binding sites in brain-aceDHSs, including for CTCF, FOXH1, and FOXQ1. Furthermore, based on GWAS, Hi-C, and eQTL data, 16 GWAS SNPs, and 82 eQTL SNPs were in brain-aceDHSs that regulate genes related to brain development or disease. Among these brain-aceDHSs, we confirmed that one enhanced the expression of GPR133, using CRISPR-Cas9 and western blotting. The GPR133 gene is associated with glioblastoma, indicating that SNPs within DHSs could be related to brain disorders. These findings suggest that brain-related gene regulatory regions are under adaptive evolution and contribute to the differential expression profiles among primates, providing new insights into the genetic basis of brain phenotypes or disorders between humans and other primates.

Integrated Epigenetic Mapping of Human and Mouse Salivary Gene Regulation

Significant effort has been applied to identify the genome-wide gene expression profiles associated with salivary gland development and pathophysiology. However, relatively little is known about the regulators that control salivary gland gene expression. We integrated data from DNase1 digital genomic footprinting, RNA-seq, and gene expression microarrays to comprehensively characterize the cis- and trans-regulatory components controlling gene expression of the healthy submandibular salivary gland. Analysis of 32 human tissues and 87 mouse tissues was performed to identify the highly expressed and tissue-enriched transcription factors driving salivary gland gene expression. Following RNA analysis, protein expression levels and subcellular localization of 39 salivary transcription factors were confirmed by immunohistochemistry. These expression analyses revealed that the salivary gland highly expresses transcription factors associated with endoplasmic reticulum stress, human T-cell lymphotrophic virus 1 expression, and Epstein-Barr virus reactivation. DNase1 digital genomic footprinting to a depth of 333,426,353 reads was performed and utilized to generate a salivary gland gene regulatory network describing the genome-wide chromatin accessibility and transcription factor binding of the salivary gland at a single-nucleotide resolution. Analysis of the DNase1 gene regulatory network identified dense interconnectivity among PLAG1, MYB, and 13 other transcription factors associated with balanced chromosomal translocations and salivary gland tumors. Collectively, these analyses provide a comprehensive atlas of the cis- and trans-regulators of the salivary gland and highlight known aberrantly regulated pathways of diseases affecting the salivary glands.

Chromatin accessibility dynamics across C. elegans development and ageing

An essential step for understanding the transcriptional circuits that control development and physiology is the global identification and characterization of regulatory elements. Here, we present the first map of regulatory elements across the development and ageing of an animal, identifying 42,245 elements accessible in at least one Caenorhabditis elegans stage. Based on nuclear transcription profiles, we define 15,714 protein-coding promoters and 19,231 putative enhancers, and find that both types of element can drive orientation-independent transcription. Additionally, more than 1000 promoters produce transcripts antisense to protein coding genes, suggesting involvement in a widespread regulatory mechanism. We find that the accessibility of most elements changes during development and/or ageing and that patterns of accessibility change are linked to specific developmental or physiological processes. The map and characterization of regulatory elements across C. elegans life provides a platform for understanding how transcription controls development and ageing.

A distinct class of antioxidant response elements is consistently activated in tumors with NRF2 mutations

NRF2 is a redox-responsive transcription factor that regulates expression of cytoprotective genes via its interaction with DNA sequences known as antioxidant response elements (AREs). NRF2 activity is induced by oxidative stress, but oxidative stress is not the only context in which NRF2 can be activated. Mutations that disrupt the interaction between NRF2 and KEAP1, an inhibitor of NRF2, lead to NRF2 hyperactivation and promote oncogenesis. The mechanisms underlying NRF2's oncogenic properties remain unclear, but likely involve aberrant expression of select NRF2 target genes. We tested this possibility using an integrative genomics approach to get a precise view of the direct NRF2 target genes dysregulated in tumors with NRF2 hyperactivating mutations. This approach revealed a core set of 32 direct NRF2 targets that are consistently upregulated in NRF2 hyperactivated tumors. This set of NRF2 "cancer target genes" includes canonical redox-related NRF2 targets, as well as target genes that have not been previously linked to NRF2 activation. Importantly, NRF2-driven upregulation of this gene set is largely independent of the organ system where the tumor developed. One key distinguishing feature of these NRF2 cancer target genes is that they are regulated by high affinity AREs that fall within genomic regions possessing a ubiquitously permissive chromatin signature. This implies that these NRF2 cancer target genes are responsive to oncogenic NRF2 in most tissues because they lack the regulatory constraints that restrict expression of most other NRF2 target genes. This NRF2 cancer target gene set also serves as a reliable proxy for NRF2 activity, and high NRF2 activity is associated with significant decreases in survival in multiple cancer types. Overall, the pervasive upregulation of these NRF2 cancer targets across multiple cancers, and their association with negative outcomes, suggests that these will be central to dissecting the functional implications of NRF2 hyperactivation in several cancer contexts.

Evolution of DNAase I Hypersensitive Sites in MHC Regulatory Regions of Primates

It has been challenging to determine the disease-causing variant(s) for most major histocompatibility complex (MHC)-associated diseases. However, it is becoming increasingly clear that regulatory variation is pervasive and a fundamentally important mechanism governing phenotypic diversity and disease susceptibility. We gathered DNase I data from 136 human cells to characterize the regulatory landscape of the MHC region, including 4867 DNase I hypersensitive sites (DHSs). We identified thousands of regulatory elements that have been gained or lost in the human or chimpanzee genomes since their evolutionary divergence. We compared alignments of the DHS across six primates and found 149 DHSs with convincing evidence of positive and/or purifying selection. Of these DHSs, compared to neutral sequences, 24 evolved rapidly in the human lineage. We identified 15 instances of transcription-factor-binding motif gains, such as USF, MYC, MAX, MAFK, STAT1, PBX3, etc, and observed 16 GWAS (genome-wide association study) SNPs associated with diseases within these 24 DHSs using FIMO (Find Individual Motif Occurrences) and UCSC (University of California, Santa Cruz) ChIP-seq data. Combining eQTL and Hi-C data, our results indicated that there were five SNPs located in human gains motifs affecting the corresponding gene's expression, two of which closely matched DHS target genes. In addition, a significant SNP, rs7756521, at genome-wide significant level likely affects DDR expression and represents a causal genetic variant for HIV-1 control. These results indicated that species-specific motif gains or losses of rapidly evolving DHSs in the primate genomes might play a role during adaptation evolution and provided some new evidence for a potentially causal role for these GWAS SNPs.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Hypomethylated domain-enriched DNA motifs prepattern the accessible nucleosome organization in teleosts

Background

Gene promoters in vertebrate genomes show distinct chromatin features such as stably positioned nucleosome array and DNA hypomethylation. The nucleosomes are known to have certain sequence preferences, and the prediction of nucleosome positioning from DNA sequence has been successful in some organisms such as yeast. However, at gene promoters where nucleosomes are much more stably positioned than in other regions, the sequence-based model has failed to work well, and sequence-independent mechanisms have been proposed.

Results

Using DNase I-seq in medaka embryos, we demonstrated that hypomethylated domains (HMDs) specifically possess accessible nucleosome organization with longer linkers, and we reassessed the DNA sequence preference for nucleosome positioning in these specific regions. Remarkably, we found with a supervised machine learning algorithm, k-mer SVM, that nucleosome positioning in HMDs is accurately predictable from DNA sequence alone. Specific short sequences (6-mers) that contribute to the prediction are specifically enriched in HMDs and distribute periodically with approximately 200-bp intervals which prepattern the position of accessible linkers. Surprisingly, the sequence preference of the nucleosome and linker in HMDs is opposite from that reported previously. Furthermore, the periodicity of specific motifs at hypomethylated promoters was conserved in zebrafish.

Conclusion

This study reveals strong link between nucleosome positioning and DNA sequence at vertebrate promoters, and we propose hypomethylated DNA-specific regulation of nucleosome positioning.

Electronic supplementary material

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

Imprinted control regions include composite DNA elements consisting of the ZFP57 binding site overlapping MLL1 morphemes

Mammalian genomes include DNA segments that are imprinted (CpG-methylated) only on one of the two parental chromosomes, leading to parent-of-origin-specific gene expression. The process is regulated by Imprinting Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs). Previously, ZFP57 was shown to recognize a methylated hexanucleotide in ICRs to maintain allele-specific gene repression. In Bioinformatics analyses, I found that the hexamer occurred frequently in mouse chromosomal DNA, suggesting that beside the ZFP57 binding site (ZFBS), ICRs contained sequence features with unknown characteristics. To identify such features, I examined chromosomal abundance of motifs in which the length of the hexamer was extended by one or several nucleotides. Results led to the discovery of a group of functionally significant composite DNA elements (ZFBS-Morph overlaps) that may play dual roles in the regulation of allele-specific gene expression. Importantly, results of genome-wide evaluations revealed that nearly 90% of the gDMRs included closely-spaced ZFBS-Morph overlaps.

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.

Genome-wide mapping of histone H3K9me2 in acute myeloid leukemia reveals large chromosomal domains associated with massive gene silencing and sites of genome instability

A facultative heterochromatin mark, histone H3 lysine 9 dimethylation (H3K9me2), which is mediated by histone methyltransferases G9a/GLP (EHMT2/1), undergoes dramatic rearrangements during myeloid cell differentiation as observed by chromatin imaging. To determine whether these structural transitions also involve genomic repositioning of H3K9me2, we used ChIP-sequencing to map genome-wide topography of H3K9me2 in normal human granulocytes, normal CD34+ hematopoietic progenitors, primary myeloblasts from acute myeloid leukemia (AML) patients, and a model leukemia cell line K562. We observe that H3K9me2 naturally repositions from the previously designated “repressed” chromatin state in hematopoietic progenitors to predominant association with heterochromatin regions in granulocytes. In contrast, AML cells accumulate H3K9me2 on previously undefined large (> 100 Kb) genomic blocks that are enriched with AML-specific single nucleotide variants, sites of chromosomal translocations, and genes downregulated in AML. Specifically, the AML-specific H3K9me2 blocks are enriched with genes regulated by the proto-oncogene ERG that promotes stem cell characteristics. The AML-enriched H3K9me2 blocks (in contrast to the heterochromatin-associated H3K9me2 blocks enriched in granulocytes) are reduced by pharmacological inhibition of the histone methyltransferase G9a/GLP in K562 cells concomitantly with transcriptional activation of ERG and ETS1 oncogenes. Our data suggest that G9a/GLP mediate formation of transient H3K9me2 blocks that are preserved in AML myeloblasts and may lead to an increased rate of AML-specific mutagenesis and chromosomal translocations.

Genome-wide analysis of p53-regulated transcription in Myc-driven lymphomas

The tumour suppressor p53 is a transcription factor that controls cellular stress responses. Here, we dissected the transcriptional programmes triggered upon restoration of p53 in Myc-driven lymphomas, based on the integrated analysis of p53 genomic occupancy and gene regulation. p53 binding sites were identified at promoters and enhancers, both characterized by the pre-existence of active chromatin marks. Only a small fraction of these sites showed the 20 base-pair p53 consensus motif, suggesting that p53 recruitment to genomic DNA was primarily mediated through protein-protein interactions in a chromatin context. p53 also targeted distal sites devoid of activation marks, at which binding was prevalently driven by sequence recognition. In all instances, the relevant motif was the canonical unsplit consensus element, with no clear evidence for p53 recruitment by split motifs. At promoters, p53 binding to the consensus motif was associated with gene induction, but not repression, indicating that the latter was most likely indirect. Altogether, our data highlight key features of genome recognition by p53 and provide unprecedented insight into the pathways associated with p53 reactivation and tumour regression, paving the way for their therapeutic application.

Distal regulation of c-myb expression during IL-6-induced differentiation in murine myeloid progenitor M1 cells

The c-Myb transcription factor is a major regulator that controls differentiation and proliferation of hematopoietic progenitor cells, which is frequently deregulated in hematological diseases, such as lymphoma and leukemia. Understanding of the mechanisms regulating the transcription of c-myb gene is challenging as it lacks a typical promoter and multiple factors are involved. Our previous studies identified some distal regulatory elements in the upstream regions of c-myb gene in murine myeloid progenitor M1 cells, but the detailed mechanisms still remain unclear. In the present study, we found that a cell differentiation-related DNase1 hypersensitive site is located at a -28k region upstream of c-myb gene and that transcription factors Hoxa9, Meis1 and PU.1 bind to the -28k region. Circular chromosome conformation capture (4C) assay confirmed the interaction between the -28k region and the c-myb promoter, which is supported by the enrichment of CTCF and Cohesin. Our analysis also points to a critical role for Hoxa9 and PU.1 in distal regulation of c-myb expression in murine myeloid cells and cell differentiation. Overexpression of Hoxa9 disrupted the IL-6-induced differentiation of M1 cells and upregulated c-myb expression through binding of the -28k region. Taken together, our results provide an evidence for critical role of the -28k region in distal regulatory mechanism for c-myb gene expression during differentiation of myeloid progenitor M1 cells.

Transcriptional regulation of the proto-oncogene Zfp521 by SPI1 (PU.1) and HOXC13

The mouse zinc‐finger gene Zfp521 (also known as ecotropic viral insertion site 3; Evi3; and ZNF521 in humans) has been identified as a B‐cell proto‐oncogene, causing leukemia in mice following retroviral insertions in its promoter region that drive Zfp521 over‐expression. Furthermore, ZNF521 is expressed in human hematopoietic cells, and translocations between ZNF521 and PAX5 are associated with pediatric acute lymphoblastic leukemia. However, the regulatory factors that control Zfp521 expression directly have not been characterized. Here we demonstrate that the transcription factors SPI1 (PU.1) and HOXC13 synergistically regulate Zfp521 expression, and identify the regions of the Zfp521 promoter required for this transcriptional activity. We also show that SPI1 and HOXC13 activate Zfp521 in a dose‐dependent manner. Our data support a role for this regulatory mechanism in vivo, as transgenic mice over‐expressing Hoxc13 in the fetal liver show a strong correlation between Hoxc13 expression levels and Zfp521 expression. Overall these experiments provide insights into the regulation of Zfp521 expression in a nononcogenic context. The identification of transcription factors capable of activating Zfp521 provides a foundation for further investigation of the regulatory mechanisms involved in ZFP521‐driven cell differentiation processes and diseases linked to Zfp521 mis‐expression.

Chromatin remodeling effects on enhancer activity

During organism development, a diversity of cell types emerges with disparate, yet stable profiles of gene expression with distinctive cellular functions. In addition to gene promoters, the genome contains enhancer regulatory sequences, which are implicated in cellular specialization by facilitating cell-type and tissue-specific gene expression. Enhancers are DNA binding elements characterized by highly sophisticated and various mechanisms of action allowing for the specific interaction of general and tissue-specific transcription factors (TFs). However, eukaryotic organisms package their genetic material into chromatin, generating a physical barrier for TFs to interact with their cognate sequences. The ability of TFs to bind DNA regulatory elements is also modulated by changes in the chromatin structure, including histone modifications, histone variants, ATP-dependent chromatin remodeling, and the methylation status of DNA. Furthermore, it has recently been revealed that enhancer sequences are also transcribed into a set of enhancer RNAs with regulatory potential. These interdependent processes act in the context of a complex network of chromatin interactions, which together contributes to a renewed vision of how gene activation is coordinated in a cell-type-dependent manner. In this review, we describe the interplay between genetic and epigenetic aspects associated with enhancers and discuss their possible roles on enhancer function.

Evidence of selection for an accessible nucleosomal array in human

BACKGROUND

Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.

RESULTS

We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.

CONCLUSION

The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.

Genome-wide analysis of p53 transcriptional programs in B cells upon exposure to genotoxic stress in vivo

The tumor suppressor p53 is a transcription factor that coordinates the cellular response to DNA damage. Here we provide an integrated analysis of p53 genomic occupancy and p53-dependent gene regulation in the splenic B and non-B cell compartments of mice exposed to whole-body ionizing radiation, providing insight into general principles of p53 activity in vivo. In unstressed conditions, p53 bound few genomic targets; induction of p53 by ionizing radiation increased the number of p53 bound sites, leading to highly overlapping profiles in the different cell types. Comparison of these profiles with chromatin features in unstressed B cells revealed that, upon activation, p53 localized at active promoters, distal enhancers, and a smaller set of unmarked distal regions. At promoters, recognition of the canonical p53 motif as well as binding strength were associated with p53-dependent transcriptional activation, but not repression, indicating that the latter was most likely indirect. p53-activated targets constituted the core of a cell type-independent response, superimposed onto a cell type-specific program. Core response genes included most of the known p53-regulated genes, as well as many new ones. Our data represent a unique characterization of the p53-regulated response to ionizing radiation in vivo.

Identifying peaks in *-seq data using shape information

Background

Peak calling is a fundamental step in the analysis of data generated by ChIP-seq or similar techniques to acquire epigenetics information. Current peak callers are often hard to parameterise and may therefore be difficult to use for non-bioinformaticians. In this paper, we present the ChIP-seq analysis tool available in CLC Genomics Workbench and CLC Genomics Server (version 7.5 and up), a user-friendly peak-caller designed to be not specific to a particular *-seq protocol.

Results

We illustrate the advantages of a shape-based approach and describe the algorithmic principles underlying the implementation. Thanks to the generality of the idea and the fact the algorithm is able to learn the peak shape from the data, the implementation requires only minimal user input, while still being applicable to a range of *-seq protocols. Using independently validated benchmark datasets, we compare our implementation to other state-of-the-art algorithms explicitly designed to analyse ChIP-seq data and provide an evaluation in terms of receiver-operator characteristic (ROC) plots. In order to show the applicability of the method to similar *-seq protocols, we also investigate algorithmic performances on DNase-seq data.

Conclusions

The results show that CLC shape-based peak caller ranks well among popular state-of-the-art peak callers while providing flexibility and ease-of-use.