Quantification of DNaseI-sensitivity by real-time PCR: quantitative analysis of DNaseI-hypersensitivity of the mouse beta-globin LCR

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.

SIRT6 promotes transcription of a subset of NRF2 targets by mono-ADP-ribosylating BAF170

SIRT6 is critical for activating transcription of Nuclear factor (erythroid-derived 2)-like 2 (NRF2) responsive genes during oxidative stress. However, while the mechanism of SIRT6-mediated silencing is well understood, the mechanism of SIRT6-mediated transcriptional activation is unknown. Here, we employed SIRT6 separation of function mutants to reveal that SIRT6 mono-ADP-ribosylation activity is required for transcriptional activation. We demonstrate that SIRT6 mono-ADP-ribosylation of BAF170, a subunit of BAF chromatin remodeling complex, is critical for activation of a subset of NRF2 responsive genes upon oxidative stress. We show that SIRT6 recruits BAF170 to enhancer region of the Heme oxygenase-1 locus and promotes recruitment of RNA polymerase II. Furthermore, SIRT6 mediates the formation of the active chromatin 10-kb loop at the HO-1 locus, which is absent in SIRT6 deficient tissue. These results provide a novel mechanism for SIRT6-mediated transcriptional activation, where SIRT6 mono-ADP-ribosylates and recruits chromatin remodeling proteins to mediate the formation of active chromatin loop.

Heterochromatin delays CRISPR-Cas9 mutagenesis but does not influence the outcome of mutagenic DNA repair

Genome editing occurs in the context of chromatin, which is heterogeneous in structure and function across the genome. Chromatin heterogeneity is thought to affect genome editing efficiency, but this has been challenging to quantify due to the presence of confounding variables. Here, we develop a method that exploits the allele-specific chromatin status of imprinted genes in order to address this problem in cycling mouse embryonic stem cells (mESCs). Because maternal and paternal alleles of imprinted genes have identical DNA sequence and are situated in the same nucleus, allele-specific differences in the frequency and spectrum of mutations induced by CRISPR-Cas9 can be unequivocally attributed to epigenetic mechanisms. We found that heterochromatin can impede mutagenesis, but to a degree that depends on other key experimental parameters. Mutagenesis was impeded by up to 7-fold when Cas9 exposure was brief and when intracellular Cas9 expression was low. In contrast, the outcome of mutagenic DNA repair was unaffected by chromatin state, with similar efficiencies of homology-directed repair (HDR) and deletion spectra on maternal and paternal chromosomes. Combined, our data show that heterochromatin imposes a permeable barrier that influences the kinetics, but not the endpoint, of CRISPR-Cas9 genome editing and suggest that therapeutic applications involving low-level Cas9 exposure will be particularly affected by chromatin status.

Chromatin status and transcription factor binding to gonadotropin promoters in gonadotrope cell lines

BACKGROUND

Proper expression of key reproductive hormones from gonadotrope cells of the pituitary is required for pubertal onset and reproduction. To further our understanding of the molecular events taking place during embryonic development, leading to expression of the glycoproteins luteinizing hormone (LH) and follicle-stimulating hormone (FSH), we characterized chromatin structure changes, imparted mainly by histone modifications, in model gonadotrope cell lines.

METHODS

We evaluated chromatin status and gene expression profiles by chromatin immunoprecipitation assays, DNase sensitivity assay, and RNA sequencing in three developmentally staged gonadotrope cell lines, αT1-1 (progenitor, expressing Cga), αT3-1 (immature, expressing Cga and Gnrhr), and LβT2 (mature, expressing Cga, Gnrhr, Lhb, and Fshb), to assess changes in chromatin status and transcription factor access of gonadotrope-specific genes.

RESULTS

We found the common mRNA α-subunit of LH and FSH, called Cga, to have an open chromatin conformation in all three cell lines. In contrast, chromatin status of Gnrhr is open only in αT3-1 and LβT2 cells. Lhb begins to open in LβT2 cells and was further opened by activin treatment. Histone H3 modifications associated with active chromatin were high on Gnrhr in αT3-1 and LβT2, and Lhb in LβT2 cells, while H3 modifications associated with repressed chromatin were low on Gnrhr, Lhb, and Fshb in LβT2 cells. Finally, chromatin status correlates with the progressive access of LHX3 to Cga and Gnrhr, followed by PITX1 binding to the Lhb promoter.

CONCLUSION

Our data show the gonadotrope-specific genes Cga, Gnrhr, Lhb, and Fshb are not only controlled by developmental transcription factors, but also by epigenetic mechanisms that include the modulation of chromatin structure, and histone modifications.

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

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

Elevated HMGN4 expression potentiates thyroid tumorigenesis

Thyroid cancer originates from genetic and epigenetic changes that alter gene expression and cellular signaling pathways. Here, we report that altered expression of the nucleosome-binding protein HMGN4 potentiates thyroid tumorigenesis. Bioinformatics analyses reveal increased HMGN4 expression in thyroid cancer. We find that upregulation of HMGN4 expression in mouse and human cells, and in the thyroid of transgenic mice, alters the cellular transcription profile, downregulates the expression of the tumor suppressors Atm, Atrx and Brca2, and elevates the levels of the DNA damage marker γH2AX. Mouse and human cells overexpressing HMGN4 show increased tumorigenicity as measured by colony formation, by tumor generation in nude mice, and by the formation of preneoplastic lesions in the thyroid of transgenic mice. Our study identifies a novel epigenetic factor that potentiates thyroid oncogenesis and raises the possibility that HMGN4 may serve as an additional diagnostic marker, or therapeutic target in certain thyroid cancers.

Remodeling of retrotransposon elements during epigenetic induction of adult visual cortical plasticity by HDAC inhibitors

BACKGROUND

The capacity for plasticity in the adult brain is limited by the anatomical traces laid down during early postnatal life. Removing certain molecular brakes, such as histone deacetylases (HDACs), has proven to be effective in recapitulating juvenile plasticity in the mature visual cortex (V1). We investigated the chromatin structure and transcriptional control by genome-wide sequencing of DNase I hypersensitive sites (DHSS) and cap analysis of gene expression (CAGE) libraries after HDAC inhibition by valproic acid (VPA) in adult V1.

RESULTS

We found that VPA reliably reactivates the critical period plasticity and induces a dramatic change of chromatin organization in V1 yielding significantly greater accessibility distant from promoters, including at enhancer regions. VPA also induces nucleosome eviction specifically from retrotransposon (in particular SINE) elements. The transiently accessible SINE elements overlap with transcription factor-binding sites of the Fox family. Mapping of transcription start site activity using CAGE revealed transcription of epigenetic and neural plasticity-regulating genes following VPA treatment, which may help to re-program the genomic landscape and reactivate plasticity in the adult cortex.

CONCLUSIONS

Treatment with HDAC inhibitors increases accessibility to enhancers and repetitive elements underlying brain-specific gene expression and reactivation of visual cortical plasticity.

Survey of protein-DNA interactions in Aspergillus oryzae on a genomic scale

The genome-scale delineation of in vivo protein–DNA interactions is key to understanding genome function. Only ∼5% of transcription factors (TFs) in the Aspergillus genus have been identified using traditional methods. Although the Aspergillus oryzae genome contains >600 TFs, knowledge of the in vivo genome-wide TF-binding sites (TFBSs) in aspergilli remains limited because of the lack of high-quality antibodies. We investigated the landscape of in vivo protein–DNA interactions across the A. oryzae genome through coupling the DNase I digestion of intact nuclei with massively parallel sequencing and the analysis of cleavage patterns in protein–DNA interactions at single-nucleotide resolution. The resulting map identified overrepresented de novo TF-binding motifs from genomic footprints, and provided the detailed chromatin remodeling patterns and the distribution of digital footprints near transcription start sites. The TFBSs of 19 known Aspergillus TFs were also identified based on DNase I digestion data surrounding potential binding sites in conjunction with TF binding specificity information. We observed that the cleavage patterns of TFBSs were dependent on the orientation of TF motifs and independent of strand orientation, consistent with the DNA shape features of binding motifs with flanking sequences.

Genome-wide mapping of DNase I hypersensitive sites in plants

Genomic regions associated with regulatory proteins are known to be highly sensitive to DNase I digestion and are termed DNase I hypersensitive sites (DHSs). DHSs can be identified by DNase I digestion followed by high-throughput DNA sequencing (DNase-seq). DNase-seq has become a powerful technique for genome-wide mapping of chromatin accessibility in eukaryotes with a sequenced genome. We have developed a DNase-seq procedure in plants. This procedure was adapted from the protocol originally developed for mammalian cell lines. It includes plant nuclei isolation, digestion of purified nuclei with DNase I, recovery of DNase-trimmed DNA fragments, DNase-seq library development, Illumina sequencing and data analysis. We also introduce a barcoding system for library preparation. We have conducted DNase-seq in both Arabidopsis thaliana and rice, and developed genome-wide open chromatin maps in both species. These DHS datasets have been used to detect footprints from regulatory protein binding and to reveal genome-wide nucleosome positioning patterns.

KLF1 stabilizes GATA-1 and TAL1 occupancy in the human β-globin locus

KLF1 is an erythroid specific transcription factor that binds to regulatory regions of erythroid genes. Binding sites of KLF1 are often found near binding sites of GATA-1 and TAL1. In the β-globin locus, KLF1 is required for forming active chromatin structure, although its role is unclear. To explore the role of KLF1 in transcribing the human γ-globin genes, we stably reduced the expression of KLF1 in erythroid K562 cells, compromising its association in the β-globin locus. The γ-globin transcription was reduced with disappearance of active chromatin structure of the locus in the KLF1 knockdown cells. Interestingly, GATA-1 and TAL1 binding was reduced in the β-globin locus, even though their expressions were not affected by KLF1 knockdown. The KLF1-dependent GATA-1 and TAL1 binding was observed in the adult locus transcribing the β-globin gene and in several erythroid genes, where GATA-1 occupancy is independent from TAL1. These results indicate that KLF1 plays a role in facilitating and/or stabilizing GATA-1 and TAL1 occupancy in the erythroid genes, contributing to the generation of active chromatin structure such as histone acetylation and chromatin looping.

Transcription factor SOX9 plays a key role in the regulation of visual cycle gene expression in the retinal pigment epithelium

The retinal pigment epithelium (RPE) performs specialized functions to support retinal photoreceptors, including regeneration of the visual chromophore. Enzymes and carrier proteins in the visual cycle function sequentially to regenerate and continuously supply 11-cis-retinal to retinal photoreceptor cells. However, it is unknown how the expression of the visual cycle genes is coordinated at the transcriptional level. Here, we show that the proximal upstream regions of six visual cycle genes contain chromatin-accessible sex-determining region Y box (SOX) binding sites, that SOX9 and LIM homeobox 2 (LHX2) are coexpressed in the nuclei of mature RPE cells, and that SOX9 acts synergistically with orthodenticle homeobox 2 (OTX2) to activate the RPE65 and retinaldehyde binding protein 1 (RLBP1) promoters and acts synergistically with LHX2 to activate the retinal G protein-coupled receptor (RGR) promoter. ChIP reveals that SOX9 and OTX2 bind to the promoter regions of RPE65, RLBP1, and RGR and that LHX2 binds to those of RPE65 and RGR in bovine RPE. ChIP with human fetal RPE cells shows that SOX9 and OTX2 also bind to the human RPE65, RLBP1, and RGR promoters. Conditional inactivation of Sox9 in mouse RPE results in reduced expression of several visual cycle genes, most dramatically Rpe65 and Rgr. Furthermore, bioinformatic analysis predicts that multiple common microRNAs (miRNAs) regulate visual cycle genes, and cotransfection of miRNA mimics with luciferase reporter constructs validated some of the predicted miRNAs. These results implicate SOX9 as a key regulator of visual cycle genes, reveal for the first time the functional role of LHX2 in the RPE, and suggest the possible regulation of visual cycle genes by common miRNAs.

The hematopoietic regulator TAL1 is required for chromatin looping between the β-globin LCR and human γ-globin genes to activate transcription

TAL1 is a key hematopoietic transcription factor that binds to regulatory regions of a large cohort of erythroid genes as part of a complex with GATA-1, LMO2 and Ldb1. The complex mediates long-range interaction between the β-globin locus control region (LCR) and active globin genes, and although TAL1 is one of the two DNA-binding complex members, its role is unclear. To explore the role of TAL1 in transcription activation of the human γ-globin genes, we reduced the expression of TAL1 in erythroid K562 cells using lentiviral short hairpin RNA, compromising its association in the β-globin locus. In the TAL1 knockdown cells, the γ-globin transcription was reduced to 35% and chromatin looping of the ^Gγ-globin gene with the LCR was disrupted with decreased occupancy of the complex member Ldb1 and LMO2 in the locus. However, GATA-1 binding, DNase I hypersensitive site formation and several histone modifications were largely maintained across the β-globin locus. In addition, overexpression of TAL1 increased the γ-globin transcription and increased interaction frequency between the ^Gγ-globin gene and LCR. These results indicate that TAL1 plays a critical role in chromatin loop formation between the γ-globin genes and LCR, which is a critical step for the transcription of the γ-globin genes.

Measuring Arabidopsis chromatin accessibility using DNase I-polymerase chain reaction and DNase I-chip assays

DNA accessibility is an important layer of regulation of DNA-dependent processes. Methods that measure DNA accessibility at local and genome-wide scales have facilitated a rapid increase in the knowledge of chromatin architecture in animal and yeast systems. In contrast, much less is known about chromatin organization in plants. We developed a robust DNase I-polymerase chain reaction (PCR) protocol for the model plant Arabidopsis (Arabidopsis thaliana). DNA accessibility is probed by digesting nuclei with a gradient of DNase I followed by locus-specific PCR. The reduction in PCR product formation along the gradient of increasing DNase I concentrations is used to determine the accessibility of the chromatin DNA. We explain a strategy to calculate the decay constant of such signal reduction as a function of increasing DNase I concentration. This allows describing DNA accessibility using a single variable: the decay constant. We also used the protocol together with AGRONOMICS1 DNA tiling microarrays to establish genome-wide DNase I sensitivity landscapes.

Heterochromatinization induced by GAA-repeat hyperexpansion in Friedreich's ataxia can be reduced upon HDAC inhibition by vitamin B3

Large intronic expansions of the triplet-repeat sequence (GAA.TTC) cause transcriptional repression of the Frataxin gene (FXN) leading to Friedreich's ataxia (FRDA). We previously found that GAA-triplet expansions stimulate heterochromatinization in vivo in transgenic mice. We report here using chromosome conformation capture (3C) coupled with high-throughput sequencing that the GAA-repeat expansion in FRDA cells stimulates a higher-order structure as a fragment containing the GAA-repeat expansion showed an increased interaction frequency with genomic regions along the FXN locus. This is consistent with a more compacted chromatin and coincided with an increase in both constitutive H3K9me3 and facultative H3K27me3 heterochromatic marks in FRDA. Consistent with this, DNase I accessibility in regions flanking the GAA repeats in patients was decreased compared with healthy controls. Strikingly, this effect could be antagonized with the class III histone deactylase (HDAC) inhibitor vitamin B3 (nicotinamide) which activated the silenced FXN gene in several FRDA models. Examination of the FXN locus revealed a reduction of H3K9me3 and H3K27me3, an increased accessibility to DNase I and an induction of euchromatic H3 and H4 histone acetylations upon nicotinamide treatment. In addition, transcriptomic analysis of nicotinamide treated and untreated FRDA primary lymphocytes revealed that the expression of 67% of genes known to be dysregulated in FRDA was ameliorated by the treatment. These findings show that nictotinamide can up-regulate the FXN gene and reveal a potential mechanism of action for nicotinamide in reactivating the epigenetically silenced FXN gene and therefore support the further assessment of HDAC inhibitors (HDACi's) in FRDA and diseases caused by a similar mechanism.

Causes and consequences of chromatin variation between inbred mice

Variation at regulatory elements, identified through hypersensitivity to digestion by DNase I, is believed to contribute to variation in complex traits, but the extent and consequences of this variation are poorly characterized. Analysis of terminally differentiated erythroblasts in eight inbred strains of mice identified reproducible variation at approximately 6% of DNase I hypersensitive sites (DHS). Only 30% of such variable DHS contain a sequence variant predictive of site variation. Nevertheless, sequence variants within variable DHS are more likely to be associated with complex traits than those in non-variant DHS, and variants associated with complex traits preferentially occur in variable DHS. Changes at a small proportion (less than 10%) of variable DHS are associated with changes in nearby transcriptional activity. Our results show that whilst DNA sequence variation is not the major determinant of variation in open chromatin, where such variants exist they are likely to be causal for complex traits.

Regulatory sites of the genome affect gene expression and complex traits, including disease susceptibility. Variable regulatory sites are potentially interesting because they are a likely cause of phenotypic variation, providing a bridge between sequence and transcriptional variation. In this paper we identify regions of the genome where DNA is not wrapped up in chromatin (hence potentially regulatory) in eight inbred strains of mice. We compare sites that vary among strains and compare them to non-variable sites. We show that more than half of variable sites cannot be attributed to local sequence variation. Functional consequences (in terms of readily detectable changes in gene expression) are associated with less than 10% of variable DNase I hypersensitive sites. We show that variable sites are enriched for sequence variants contributing to complex traits in mice.

Epigenetic modifications unlock the milk protein gene loci during mouse mammary gland development and differentiation

BACKGROUND

Unlike other tissues, development and differentiation of the mammary gland occur mostly after birth. The roles of systemic hormones and local growth factors important for this development and functional differentiation are well-studied. In other tissues, it has been shown that chromatin organization plays a key role in transcriptional regulation and underlies epigenetic regulation during development and differentiation. However, the role of chromatin organization in mammary gland development and differentiation is less well-defined. Here, we have studied the changes in chromatin organization at the milk protein gene loci (casein, whey acidic protein, and others) in the mouse mammary gland before and after functional differentiation.

METHODOLOGY/PRINCIPAL FINDINGS

Distal regulatory elements within the casein gene cluster and whey acidic protein gene region have an open chromatin organization after pubertal development, while proximal promoters only gain open-chromatin marks during pregnancy in conjunction with the major induction of their expression. In contrast, other milk protein genes, such as alpha-lactalbumin, already have an open chromatin organization in the mature virgin gland. Changes in chromatin organization in the casein gene cluster region that are present after puberty persisted after lactation has ceased, while the changes which occurred during pregnancy at the gene promoters were not maintained. In general, mammary gland expressed genes and their regulatory elements exhibit developmental stage- and tissue-specific chromatin organization.

CONCLUSIONS/SIGNIFICANCE

A progressive gain of epigenetic marks indicative of open/active chromatin on genes marking functional differentiation accompanies the development of the mammary gland. These results support a model in which a chromatin organization is established during pubertal development that is then poised to respond to the systemic hormonal signals of pregnancy and lactation to achieve the full functional capacity of the mammary gland.

DNase I digestion of isolated nulcei for genome-wide mapping of DNase hypersensitivity sites in chromatin

DNase I hypersensitivity (DHS) analysis is a powerful method to analyze chromatin structure and identify genomic regulatory elements. Integration of a high-throughput detection method into DHS analysis makes genome-wide mapping of DHS sites possible at a reasonable cost. Here we describe methods for DHS analysis carried out with mouse liver nuclei, involving DNase I digestion followed by isolation of DNase I-released DNA fragments suitable for high-throughput, next generation DNA sequencing (DNase-seq). A real-time PCR-based assay used to optimize DNase I digestion conditions is also described.

Epigenetic expansion of VHL-HIF signal output drives multiorgan metastasis in renal cancer

Inactivation of the von Hippel-Lindau tumor suppressor gene, VHL, is an archetypical tumor-initiating event in clear cell renal carcinoma (ccRCC) that leads to the activation of hypoxia-inducible transcription factors (HIFs). However, VHL mutation status in ccRCC is not correlated with clinical outcome. Here we show that during ccRCC progression, cancer cells exploit diverse epigenetic alterations to empower a branch of the VHL-HIF pathway for metastasis, and the strength of this activation is associated with poor clinical outcome. By analyzing metastatic subpopulations of VHL-deficient ccRCC cells, we discovered an epigenetically altered VHL-HIF response that is specific to metastatic ccRCC. Focusing on the two most prominent pro-metastatic VHL-HIF target genes, we show that loss of Polycomb repressive complex 2 (PRC2)-dependent histone H3 Lys27 trimethylation (H3K27me3) activates HIF-driven chemokine (C-X-C motif) receptor 4 (CXCR4) expression in support of chemotactic cell invasion, whereas loss of DNA methylation enables HIF-driven cytohesin 1 interacting protein (CYTIP) expression to protect cancer cells from death cytokine signals. Thus, metastasis in ccRCC is based on an epigenetically expanded output of the tumor-initiating pathway.

Novel evolutionary-conserved role for the activity-dependent neuroprotective protein (ADNP) family that is important for erythropoiesis

BACKGROUND

ADNP is vital for embryonic development. Is this function conserved for the homologous protein ADNP2?

RESULTS

Down-regulation/silencing of ADNP or ADNP2 in zebrafish embryos or mouse erythroleukemia cells inhibited erythroid maturation, with ADNP directly associating with the β-globin locus control region.

CONCLUSION

ADNPs are novel molecular regulators of erythropoiesis.

SIGNIFICANCE

New regulators of globin synthesis are suggested. Activity-dependent neuroprotective protein (ADNP) and its homologue ADNP2 belong to a homeodomain, the zinc finger-containing protein family. ADNP is essential for mouse embryonic brain formation. ADNP2 is associated with cell survival, but its role in embryogenesis has not been evaluated. Here, we describe the use of the zebrafish model to elucidate the developmental roles of ADNP and ADNP2. Although we expected brain defects, we were astonished to discover that the knockdown zebrafish embryos were actually lacking blood and suffered from defective hemoglobin production. Evolutionary conservation was established using mouse erythroleukemia (MEL) cells, a well studied erythropoiesis model, in which silencing of ADNP or ADNP2 produced similar results as in zebrafish. Exogenous RNA encoding ADNP/ADNP2 rescued the MEL cell undifferentiated state, demonstrating phenotype specificity. Brg1, an ADNP-interacting chromatin-remodeling protein involved in erythropoiesis through regulation of the globin locus, was shown here to interact also with ADNP2. Furthermore, chromatin immunoprecipitation revealed recruitment of ADNP, similar to Brg1, to the mouse β-globin locus control region in MEL cells. This recruitment was apparently diminished upon dimethyl sulfoxide (DMSO)-induced erythrocyte differentiation compared with the nondifferentiated state. Importantly, exogenous RNA encoding ADNP/ADNP2 significantly increased β-globin expression in MEL cells in the absence of any other differentiation factors. Taken together, our results reveal an ancestral role for the ADNP protein family in maturation and differentiation of the erythroid lineage, associated with direct regulation of β-globin expression.

The DNA demethylating agent decitabine activates the TRAIL pathway and induces apoptosis in acute myeloid leukemia

Although epigenetic drugs have been approved for use in selected malignancies, there is significant need for a better understanding of their mechanism of action. Here, we study the action of a clinically approved DNA-methyltransferase inhibitor - decitabine (DAC) - in acute myeloid leukemia (AML) cells. At low doses, DAC treatment induced apoptosis of NB4 Acute Promyelocytic Leukemia (APL) cells, which was associated with the activation of the extrinsic apoptotic pathway. Expression studies of the members of the Death Receptor family demonstrated that DAC induces the expression of TNF-related apoptosis-inducing ligand (TRAIL). Upregulation of TRAIL, upon DAC treatment, was associated with specific epigenetic modifications induced by DAC in the proximity of the TRAIL promoter, as demonstrated by DNA demethylation, increased DNaseI sensitivity and histone acetylation of a non-CpG island, CpG-rich region located 2kb upstream to the transcription start site. Luciferase assay experiments showed that this region behave as a DNA methylation sensitive transcriptional regulatory element. The CpG regulatory element was also found methylated in samples derived from APL patients. These findings have been confirmed in the non-APL, AML Kasumi cell line, suggesting that this regulatory mechanism may be extended to other AMLs. Our study suggests that DNA methylation is a regulatory mechanism relevant for silencing of the TRAIL apoptotic pathway in leukemic cells, and further elucidates the mechanism by which epigenetic drugs mediate their anti-leukemic effects.

Regulation of delta-aminolevulinic acid dehydratase by krüppel-like factor 1

Krüppel-like factor 1(KLF1) is a hematopoietic-specific zinc finger transcription factor essential for erythroid gene expression. In concert with the transacting factor GATA1, KLF1 modulates the coordinate expression of the genes encoding the multi-enzyme heme biosynthetic pathway during erythroid differentiation. To explore the mechanisms underpinning KLF1 action at the gene loci regulating the first 3 steps in this process, we have exploited the K1-ERp erythroid cell line, in which KLF1 translocates rapidly to the nucleus in response to treatment with 4-OH-Tamoxifen (4-OHT). KLF1 acts as a differentiation-independent transcriptional co-regulator of delta-aminolevulinic acid dehydratase (Alad), but not 5-aminolevulinate synthase gene (Alas2) or porphobilinogen deaminase (Pbgd). Similar to its role at the β-globin promoter, KLF1 induces factor recruitment and chromatin changes at the Alad1b promoter in a temporally-specific manner. In contrast to these changes, we observed a distinct mechanism of histone eviction at the Alad1b promoter. Furthermore, KLF1-dependent events were not modulated by GATA1 factor promoter co-occupancy alone. These results not only enhance our understanding of erythroid-specific modulation of heme biosynthetic regulation by KLF1, but provide a model that will facilitate the elucidation of novel KLF1-dependent events at erythroid gene loci that are independent of GATA1 activity.

Chromatin structure of the LCR in the human β-globin locus transcribing the adult δ- and β-globin genes

The β-like globin genes are transcribed in a developmental stage specific fashion in erythroid cells. The specific transcription of globin genes is conferred by the locus control region (LCR), but the chromatin structure of the LCR in the human adult β-globin locus transcribing the δ- and β-globin genes is not clear. Here, we employed hybrid MEL cells that contain a human chromosome 11. The δ- and β-globin genes were highly transcribed in hybrid MEL/ch11 cells after transcriptional induction. LCR HS3 and HS2 were strongly occupied by erythroid specific transcriptional activators and co-factors in the induced locus. These HSs, but not HS4 and HS1, were in close proximity with the active globin genes as revealed by high resolution 3C experiments. The active features at HS3 were markedly established after transcriptional induction, while HS2 was in a relatively active conformation before the induction. Unexpectedly, HS1 did not show notable active features except histone hyperacetylation. Taken together, the LCR of the human β-globin locus transcribing the adult δ- and β-globin genes has HS specific chromatin structure. The structure at each HS, which is different from the locus transcribing the fetal globin genes, might relate to its role in transcribing the adult genes.

High-resolution mapping of open chromatin in the rice genome

Gene expression is controlled by the complex interaction of transcription factors binding to promoters and other regulatory DNA elements. One common characteristic of the genomic regions associated with regulatory proteins is a pronounced sensitivity to DNase I digestion. We generated genome-wide high-resolution maps of DNase I hypersensitive (DH) sites from both seedling and callus tissues of rice (Oryza sativa). Approximately 25% of the DH sites from both tissues were found in putative promoters, indicating that the vast majority of the gene regulatory elements in rice are not located in promoter regions. We found 58% more DH sites in the callus than in the seedling. For DH sites detected in both the seedling and callus, 31% displayed significantly different levels of DNase I sensitivity within the two tissues. Genes that are differentially expressed in the seedling and callus were frequently associated with DH sites in both tissues. The DNA sequences contained within the DH sites were hypomethylated, consistent with what is known about active gene regulatory elements. Interestingly, tissue-specific DH sites located in the promoters showed a higher level of DNA methylation than the average DNA methylation level of all the DH sites located in the promoters. A distinct elevation of H3K27me3 was associated with intergenic DH sites. These results suggest that epigenetic modifications play a role in the dynamic changes of the numbers and DNase I sensitivity of DH sites during development.

Epigenetic priors for identifying active transcription factor binding sites

MOTIVATION

Accurate knowledge of the genome-wide binding of transcription factors in a particular cell type or under a particular condition is necessary for understanding transcriptional regulation. Using epigenetic data such as histone modification and DNase I, accessibility data has been shown to improve motif-based in silico methods for predicting such binding, but this approach has not yet been fully explored.

RESULTS

We describe a probabilistic method for combining one or more tracks of epigenetic data with a standard DNA sequence motif model to improve our ability to identify active transcription factor binding sites (TFBSs). We convert each data type into a position-specific probabilistic prior and combine these priors with a traditional probabilistic motif model to compute a log-posterior odds score. Our experiments, using histone modifications H3K4me1, H3K4me3, H3K9ac and H3K27ac, as well as DNase I sensitivity, show conclusively that the log-posterior odds score consistently outperforms a simple binary filter based on the same data. We also show that our approach performs competitively with a more complex method, CENTIPEDE, and suggest that the relative simplicity of the log-posterior odds scoring method makes it an appealing and very general method for identifying functional TFBSs on the basis of DNA and epigenetic evidence.

AVAILABILITY AND IMPLEMENTATION

FIMO, part of the MEME Suite software toolkit, now supports log-posterior odds scoring using position-specific priors for motif search. A web server and source code are available at http://meme.nbcr.net. Utilities for creating priors are at http://research.imb.uq.edu.au/t.bailey/SD/Cuellar2011.

CONTACT

t.bailey@uq.edu.au

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

NF-κB and BRG1 bind a distal regulatory element in the IL-3/GM-CSF locus

We investigated gene regulation at the IL-3/GM-CSF gene cluster. We found BRG1, a SWI/SNF remodeling ATPase, bound a distal element, CNSa. BRG1 binding was strongest in differentiated, stimulated T helper cells, paralleling IL-3 and GM-CSF expression. Depletion of BRG1 reduced IL-3 and GM-CSF transcription. BAF-specific SWI/SNF subunits bound to this locus and regulated IL-3 expression. CNSa was in closed chromatin in fibroblasts, open chromatin in differentiated T helper cells, and moderately open chromatin in naïve (undifferentiated) T helper cells; BRG1 was required for the most open state. CNSa increased transcription of a reporter in an episomal expression system, in a BRG1-dependent manner. The NF-κB subunit RelA/p65 bound CNSa in activated T helper cells. Inhibition of NF-κB blocked BRG1 binding to CNSa, chromatin opening at CNSa, and activation of IL-3 and GM-CSF. Together, these findings suggest CNSa is a distal enhancer that binds BRG1 and NF-κB.

The distinctive roles of erythroid specific activator GATA-1 and NF-E2 in transcription of the human fetal γ-globin genes

GATA-1 and NF-E2 are erythroid specific activators that bind to the β-globin locus. To explore the roles of these activators in transcription of the human fetal stage specific γ-globin genes, we reduced GATA-1 and p45/NF-E2 using shRNA in erythroid K562 cells. GATA-1 or p45/NF-E2 knockdown inhibited the transcription of the γ-globin genes, hypersensitive site (HS) formation in the LCR and chromatin loop formation of the β-globin locus, but histone acetylation across the locus was decreased only in the case of GATA-1 knockdown. In p45/NF-E2 knockdown cells, GATA-1 binding was maintained at the LCR HSs and γ-globin promoter, but NF-E2 binding at the LCR HSs was reduced by GATA-1 knockdown regardless of the amount of p45/NF-E2 in K562 cells. These results indicate that histone acetylation is dependent on GATA-1 binding, but the binding of GATA-1 is not sufficient for the γ-globin transcription, HS formation and chromatin loop formation and NF-E2 is required. This idea is supported by the distinctive binding pattern of CBP and Brg1 in the β-globin locus. Furthermore GATA-1-dependent loop formation between HS5 and 3′HS1 suggests correlation between histone modifications and chromatin looping.

An embryonic stage-specific enhancer within the murine β-globin locus mediates domain-wide histone hyperacetylation

In mammalian nuclei, a select number of tissue-specific gene loci exhibit broadly distributed patterns of histone modifications, such as histone hyperacetylation, that are normally associated with active gene promoters. Previously, we characterized such hyperacetylated domains within mammalian β-globin gene loci, and determined that within the murine locus, neither the β-globin locus control region nor the gene promoters were required for domain formation. Here, we identify a developmentally specific erythroid enhancer, hypersensitive site-embryonic 1 (HS-E1), located within the embryonic β-globin domain in mouse, which is homologous to a region located downstream of the human embryonic ε-globin gene. This sequence exhibits nuclease hypersensitivity in primitive erythroid cells and acts as an enhancer in gain-of-function assays. Deletion of HS-E1 from the endogenous murine β-globin locus results in significant decrease in the expression of the embryonic β-globin genes and loss of the domain-wide pattern of histone hyperacetylation. The data suggest that HS-E1 is an enhancer that is uniquely required for β-like globin expression in primitive erythroid cells, and that it defines a novel class of enhancer that works in part by domain-wide modulation of chromatin structure.

Dynamic chromatin modifications control GnRH gene expression during neuronal differentiation and protein kinase C signal transduction

GnRH, a neuropeptide produced by rare, specialized hypothalamic secretory neurons, is critical for reproduction. During development, GnRH gene expression increases as neurons migrate from the olfactory placode to the hypothalamus, with highest levels in the mature, postmitotic state. While neuronal differentiation is known to be controlled by chromatin modulations, the role of chromatin dynamics in GnRH gene regulation has not been studied. Here, we use mature and immature GnRH neuronal cell models to show that both neuron-specific and protein kinase C regulation of GnRH expression are mediated by chromatin structure and histone modifications. Only in GT1-7 mature GnRH neuronal cells did GnRH regulatory elements display high sensitivity to DNase and enrichment of active histone markers histone-H3 acetylation and H3 lysine 4 trimethylation (H3K4-Me3), as well as RNA polymerase II (RNAPII) binding and enhancer RNA transcription. In contrast, H3K9-Me2, a marker of inactive chromatin, was highest in nonneuronal cells, low in GT1-7 cells, and intermediate in immature GnRH neuronal cells. The chromatin of the GnRH gene was therefore active in mature GnRH neuronal cells, inactive in nonneuronal cells, but not fully inactive in immature GnRH neuronal cells. Activation of protein kinase C (PKC) potently represses GnRH expression. PKC activation caused closing of the chromatin and decreased RNAPII occupancy at the GnRH minimal promoter (-278/-97). At GnRH-Enhancer-1 (-2404/-2100), PKC activation decreased phosphorylated-RNAPII binding, enhancer RNA transcription, and H3 acetylation, and reciprocally increased H3K9-Me2. Chromatin modifications therefore participate in the dynamic regulation and specification of GnRH expression to differentiated hypothalamic neurons.

Nuclear role of WASp in the pathogenesis of dysregulated TH1 immunity in human Wiskott-Aldrich syndrome

The clinical symptomatology in the X-linked Wiskott-Aldrich syndrome (WAS), a combined immunodeficiency and autoimmune disease resulting from WAS protein (WASp) deficiency, reflects the underlying coexistence of an impaired T helper 1 (TH1) immunity alongside intact TH2 immunity. This suggests a role for WASp in patterning T(H) subtype immunity, yet the molecular basis for the TH1-TH2 imbalance in human WAS is unknown. We have discovered a nuclear role for WASp in the transcriptional regulation of the TH1 regulator gene TBX21 at the chromatin level. In primary TH1-differentiating cells, a fraction of WASp is found in the nucleus, where it is recruited to the proximal promoter locus of the TBX21 gene, but not to the core promoter of GATA3 (a TH2 regulator gene) or RORc (a TH17 regulator gene). Genome-wide mapping demonstrates association of WASp in vivo with the gene-regulatory network that orchestrates TH1 cell fate choice in the human TH cell genome. Functionally, nuclear WASp associates with H3K4 trimethyltransferase [RBBP5 (retinoblastoma-binding protein 5)] and H3K9/H3K36 tridemethylase [JMJD2A (Jumonji domain-containing protein 2A)] proteins, and their enzymatic activity in vitro and in vivo is required for achieving transcription-permissive chromatin dynamics at the TBX21 proximal promoter in primary differentiating TH1 cells. During TH1 differentiation, the loss of WASp accompanies decreased enrichment of RBBP5 and, in a subset of WAS patients, also of filamentous actin at the TBX21 proximal promoter locus. Accordingly, human WASp-deficient TH cells, from natural mutation or RNA interference-mediated depletion, demonstrate repressed TBX21 promoter dynamics when driven under TH1-differentiating conditions. These chromatin derangements accompany deficient T-BET messenger RNA and protein expression and impaired TH1 function, defects that are ameliorated by reintroducing WASp. Our findings reveal a previously unappreciated role of WASp in the epigenetic control of T-BET transcription and provide a new mechanism for the pathogenesis of WAS by linking aberrant histone methylation at the TBX21 promoter to dysregulated adaptive immunity.

A key commitment step in erythropoiesis is synchronized with the cell cycle clock through mutual inhibition between PU.1 and S-phase progression

Hematopoietic progenitors undergo differentiation while navigating several cell division cycles, but it is unknown whether these two processes are coupled. We addressed this question by studying erythropoiesis in mouse fetal liver in vivo. We found that the initial upregulation of cell surface CD71 identifies developmentally matched erythroblasts that are tightly synchronized in S-phase. We show that DNA replication within this but not subsequent cycles is required for a differentiation switch comprising rapid and simultaneous committal transitions whose precise timing was previously unknown. These include the onset of erythropoietin dependence, activation of the erythroid master transcriptional regulator GATA-1, and a switch to an active chromatin conformation at the β-globin locus. Specifically, S-phase progression is required for the formation of DNase I hypersensitive sites and for DNA demethylation at this locus. Mechanistically, we show that S-phase progression during this key committal step is dependent on downregulation of the cyclin-dependent kinase p57(KIP2) and in turn causes the downregulation of PU.1, an antagonist of GATA-1 function. These findings therefore highlight a novel role for a cyclin-dependent kinase inhibitor in differentiation, distinct to their known function in cell cycle exit. Furthermore, we show that a novel, mutual inhibition between PU.1 expression and S-phase progression provides a "synchromesh" mechanism that "locks" the erythroid differentiation program to the cell cycle clock, ensuring precise coordination of critical differentiation events.

A highly conserved SOX6 double binding site mediates SOX6 gene downregulation in erythroid cells

The Sox6 transcription factor plays critical roles in various cell types, including erythroid cells. Sox6-deficient mice are anemic due to impaired red cell maturation and show inappropriate globin gene expression in definitive erythrocytes. To identify new Sox6 target genes in erythroid cells, we used the known repressive double Sox6 consensus within the εy-globin promoter to perform a bioinformatic genome-wide search for similar, evolutionarily conserved motifs located within genes whose expression changes during erythropoiesis. We found a highly conserved Sox6 consensus within the Sox6 human gene promoter itself. This sequence is bound by Sox6 in vitro and in vivo, and mediates transcriptional repression in transient transfections in human erythroleukemic K562 cells and in primary erythroblasts. The binding of a lentiviral transduced Sox6FLAG protein to the endogenous Sox6 promoter is accompanied, in erythroid cells, by strong downregulation of the endogenous Sox6 transcript and by decreased in vivo chromatin accessibility of this region to the PstI restriction enzyme. These observations suggest that the negative Sox6 autoregulation, mediated by the double Sox6 binding site within its own promoter, may be relevant to control the Sox6 transcriptional downregulation that we observe in human erythroid cultures and in mouse bone marrow cells in late erythroid maturation.

An epigenetic chromatin remodeling role for NFATc1 in transcriptional regulation of growth and survival genes in diffuse large B-cell lymphomas

The nuclear factor of activated T cells (NFAT) family of transcription factors functions as integrators of multiple signaling pathways by binding to chromatin in combination with other transcription factors and coactivators to regulate genes central for cell growth and survival in hematopoietic cells. Recent experimental evidence has implicated the calcineurin/NFAT signaling pathway in the pathogenesis of various malignancies, including diffuse large B-cell lymphoma (DLBCL). However, the molecular mechanism(s) underlying NFATc1 regulation of genes controlling lymphoma cell growth and survival is still unclear. In this study, we demonstrate that the transcription factor NFATc1 regulates gene expression in DLBCL cells through a chromatin remodeling mechanism that involves recruitment of the SWItch/Sucrose NonFermentable chromatin remodeling complex ATPase enzyme SMARCA4 (also known as Brahma-related gene 1) to NFATc1 targeted gene promoters. The NFATc1/Brahma-related gene 1 complex induces promoter DNase I hypersensitive sites and recruits other transcription factors to the active chromatin site to regulate gene transcription. Targeting NFATc1 with specific small hairpin RNA inhibits DNase I hypersensitive site formation and down-regulates target gene expression. Our data support a novel epigenetic control mechanism for the transcriptional regulation of growth and survival genes by NFATc1 in the pathophysiology of DLBCL and suggests that targeting NFATc1 could potentially have therapeutic value.

Integration of distinct intracellular signaling pathways at distal regulatory elements directs T-bet expression in human CD4+ T cells

T-bet is a key regulator controlling Th1 cell development. This factor is not expressed in naive CD4(+) T cells, and the mechanisms controlling expression of T-bet are incompletely understood. In this study, we defined regulatory elements at the human T-bet locus and determined how signals originating at the TCR and at cytokine receptors are integrated to induce chromatin modifications and expression of this gene during human Th1 cell differentiation. We found that T cell activation induced two strong DNase I-hypersensitive sites (HS) and rapid histone acetylation at these elements in CD4(+) T cells. Histone acetylation and T-bet expression were strongly inhibited by cyclosporine A, and we detected binding of NF-AT to a HS in vivo. IL-12 and IFN-gamma signaling alone were not sufficient to induce T-bet expression in naive CD4(+) T cells, but enhanced T-bet expression in TCR/CD28-stimulated cells. We detected a third HS 12 kb upstream of the mRNA start site only in developing Th1 cells, which was bound by IL-12-induced STAT4. Our data suggest that T-bet locus remodeling and gene expression are initiated by TCR-induced NF-AT recruitment and amplified by IL-12-mediated STAT4 binding to distinct distal regulatory elements during human Th1 cell differentiation.

BRG1 directly regulates nucleosome structure and chromatin looping of the alpha globin locus to activate transcription

α globin expression must be regulated properly to prevent the occurrence of α-thalassemias, yet many questions remain unanswered regarding the mechanism of transcriptional activation. Identifying factors that regulate chromatin structure of the endogenous α globin locus in developing erythroblasts will provide important mechanistic insight. Here, we demonstrate that the BRG1 catalytic subunit of SWI/SNF-related complexes co-immunoprecipitates with GATA-1 and EKLF in murine fetal liver cells in vivo and is recruited to the far-upstream major-regulatory element (MRE) and α2 promoter. Furthermore, based on our analysis of Brg1^null/ENU1 mutant mice, BRG1 regulates DNase I sensitivity, H3ac, and H3K4me2 but not CpG methylation at both sites. Most importantly, BRG1 is required for chromatin loop formation between the MRE and α2 promoter and for maximal RNA Polymerase II occupancy at the α2 promoter. Consequently, Brg1 mutants express α globin mRNA at only 5–10% of wild-type levels and die at mid-gestation. These data identify BRG1 as a chromatin-modifying factor required for nucleosome remodeling and transcriptional activation of the α globin locus. These data also demonstrate that chromatin looping between the MRE and α2 promoter is required as part of the transcriptional activation mechanism.

The HBS1L-MYB intergenic interval associated with elevated HbF levels shows characteristics of a distal regulatory region in erythroid cells

HBS1L-MYB intergenic polymorphism (HMIP) on chromosome 6q23 is associated with elevated fetal hemoglobin levels and has pleiotropic effects on several hematologic parameters. To investigate potential regulatory activity in the region, we have measured sensitivity of the sequences to DNase I cleavage that identified 3 tissue-specific DNase I hypersensitive sites in the core intergenic interval. Chromatin immunoprecipitation with microarray (ChIP-chip) analysis showed strong histone acetylation in a defined interval of 65 kb corresponding to the core HBS1L-MYB intergenic region in primary human erythroid cells but not in non–MYB-expressing HeLa cells. ChIP-chip analysis also identified several potential cis-regulatory elements as strong GATA-1 signals that coincided with the DNase I hypersensitive sites present in MYB-expressing erythroid cells. We suggest that HMIP contains regulatory sequences that could be important in hematopoiesis by controlling MYB expression. This study provides the functional link between genetic association of HMIP with control of fetal hemoglobin and other hematologic parameters. We also present a large-scale analysis of histone acetylation as well as RNA polymerase II and GATA-1 interactions on chromosome 6q, and α and β globin gene loci. The data suggest that GATA-1 regulates numerous genes of various functions on chromosome 6q.

Genetic and epigenetic networks controlling T helper 1 cell differentiation

Significant progress has been made during the past years in our understanding of the mechanisms that control the differentiation of naïve CD4(+) T cells into effector T-cell subsets with distinct functional properties. Previous work allowed the identification of key molecules involved in regulating this highly complex process, such as cytokines and their receptors, signal transducers and transcription factors. More recently, the emphasis of research in this field has been to elucidate how the multiplicity of signals is integrated to shape a T helper subset-specific gene-expression program controlling differentiation and effector functions. In this review we will highlight advances that have been made in unravelling the genetic and epigenetic networks controlling differentiation of naïve CD4(+) T cells into interferon-gamma(IFN-gamma)-secreting T helper type 1 (Th1) cells.

A 220-nucleotide deletion of the intronic enhancer reveals an epigenetic hierarchy in immunoglobulin heavy chain locus activation

A tissue-specific transcriptional enhancer, Eμ, has been implicated in developmentally regulated recombination and transcription of the immunoglobulin heavy chain (IgH) gene locus. We demonstrate that deleting 220 nucleotides that constitute the core Eμ results in partially active locus, characterized by reduced histone acetylation, chromatin remodeling, transcription, and recombination, whereas other hallmarks of tissue-specific locus activation, such as loss of H3K9 dimethylation or gain of H3K4 dimethylation, are less affected. These observations define Eμ-independent and Eμ-dependent phases of locus activation that reveal an unappreciated epigenetic hierarchy in tissue-specific gene expression.

Chromatin and DNA modifications in the Opaque2-mediated regulation of gene transcription during maize endosperm development

The maize (Zea mays) Opaque2 (O2) gene encodes an endosperm-specific bZIP-type transcription activator. In this study, we analyzed O2 targets for chromatin and DNA modifications and transcription factors binding during endosperm development and in leaves. In leaves, O2 targets exhibit high cytosine methylation levels and transcriptionally silent chromatin, enriched with histones H3 dimethylated at Lys-9 (H3K9me2) and Lys-27 (H3K27me2). Transcriptional activation in the endosperm occurs through a two-step process, with an early potentiated state and a later activated state. The potentiated state has cytosine demethylation at symmetric sites, substitution of H3K9me2 and H3K27me2 with histones H3 acetylated at Lys-14 (H3K14ac) and dimethylated at Lys-4 (H3K4me2), and increased DNaseI sensitivity. During the activated state, the mRNA of O2 targets accumulates in correspondence to RNPII, O2, and Ada2/Gcn5 coactivator binding. The active state also exhibits further increases of H3K14ac/H3K4me2 and DNaseI accessibility levels and deposition of histone H3 acetylated at Lys-9 and trimethylated at Lys-4. Analysis of o2 mutants revealed that O2 targets differ in their dependence on O2 activity for coactivator recruitment and for formation of specific chromatin modification profiles. These results indicate gene-specific involvement of mechanisms that modify chromatin states in the O2-mediated regulation of transcription.

NA-Seq: a discovery tool for the analysis of chromatin structure and dynamics during differentiation

It is well established that epigenetic modulation of genome accessibility in chromatin occurs during biological processes. Here we describe a method based on restriction enzymes and next-generation sequencing for identifying accessible DNA elements using a small amount of starting material, and use it to examine myeloid differentiation of primary human CD34+ cells. The accessibility of several classes of cis-regulatory elements was a predictive marker of in vivo DNA binding by transcription factors, and was associated with distinct patterns of histone posttranslational modifications. We also mapped large chromosomal domains with differential accessibility in progenitors and maturing cells. Accessibility became restricted during differentiation, correlating with a decreased number of expressed genes and loss of regulatory potential. Our data suggest that a permissive chromatin structure in multipotent cells is progressively and selectively closed during differentiation, and illustrate the use of our method for the identification of functional cis-regulatory elements.

Genome-scale identification of Caenorhabditis elegans regulatory elements by tiling-array mapping of DNase I hypersensitive sites

Background

A major goal of post-genomics research is the integrated analysis of genes, regulatory elements and the chromatin architecture on a genome-wide scale. Mapping DNase I hypersensitive sites within the nuclear chromatin is a powerful and well-established method of identifying regulatory element candidates.

Results

Here, we report the first genome-wide analysis of DNase I hypersensitive sites (DHSs) in Caenorhabditis elegans. The data was obtained by hybridizing DNase I-treated and end-captured material from young adult worms to a high-resolution tiling microarray. The data show that C. elegans DHSs were significantly enriched within intergenic regions located 2 kb upstream and downstream of coding genes, and also that a considerable fraction of all DHSs mapped to intergenic positions distant to annotated coding genes. Annotated transcribed loci were generally depleted in DHSs relative to intergenic regions, but DHSs were nonetheless enriched in coding exons and UTRs, whereas introns were significantly depleted in DHSs. Many DHSs appeared to be associated with annotated non-coding RNAs and recently detected transcripts of unknown function. It has been reported that nematode highly conserved non-coding elements were associated with cis-regulatory elements, and we also found that DHSs, particularly distal intergenic DHSs, were significantly enriched in regions that were conserved between the C. elegans and C. briggsae genomes.

Conclusion

We describe the first genome-wide analysis of C. elegans DHSs, and show that the distribution of DHSs is strongly associated with functional elements in the genome.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells

The fraction of experimentally active conserved non-coding sequences within any given cell type is low, so classical assays are unlikely to expose their potential.

Background

Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian.

Results

We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types.

Conclusions

The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Identification of long range regulatory elements of mouse alpha-globin gene cluster by quantitative associated chromatin trap (QACT)

Chromatin from different regions of the genome frequently forms steady associations that play important roles in regulating gene expression. The widely used chromatin conformation capture (3C) assay allows determination of the in vivo structural organization of an active endogenous locus. However, unpredicted chromatin associations within a given genomic locus can not be identified by 3C. Here, we describe a new strategy, quantitative associated chromatin trap (QACT), which incorporates a modified 3C method and a quantitative assay tool, to capture and quantitatively analyzes all possible associated chromatin partners (ACPs) of a given chromatin fragment. Using QACT, we have analyzed the chromatin conformation of the mouse alpha-globin gene cluster and proved the extensive interaction between HS26 and alpha-globin genes. In addition, we have identified a candidate alpha1-globin gene specific silencer 475A8 which shows the differentiation-stage specific DNase I hypersensitivity. Functional analysis suggests that 475A8 may regulate the alpha1-globin gene during terminal differentiation of committed erythroid progenitor cells. ChIP (chromatin immunoprecipitation) and cotransfection assays demonstrate that GATA-1, a hemopoietic specific transcriptional factor, may increase alpha1-globin gene expression by suppressing the function of 475A8 in terminally differentiated erythroid cells.

Long-range enhancer associated with chromatin looping allows AP-1 regulation of the peptidylarginine deiminase 3 gene in differentiated keratinocyte

Transcription control at a distance is a critical mechanism, particularly for contiguous genes. The peptidylarginine deiminases (PADs) catalyse the conversion of protein-bound arginine into citrulline (deimination), a critical reaction in the pathophysiology of multiple sclerosis, Alzheimer's disease and rheumatoid arthritis, and in the metabolism of the major epidermal barrier protein filaggrin, a strong predisposing factor for atopic dermatitis. PADs are encoded by 5 clustered PADI genes (1p35-6). Unclear are the mechanisms controlling the expression of the gene PADI3 encoding the PAD3 isoform, a strong candidate for the deimination of filaggrin in the terminally differentiating epidermal keratinocyte. We describe the first PAD Intergenic Enhancer (PIE), an evolutionary conserved non coding segment located 86-kb from the PADI3 promoter. PIE is a strong enhancer of the PADI3 promoter in Ca2+-differentiated epidermal keratinocytes, and requires bound AP-1 factors, namely c-Jun and c-Fos. As compared to proliferative keratinocytes, calcium stimulation specifically associates with increased local DNase I hypersensitivity around PIE, and increased physical proximity of PIE and PADI3 as assessed by Chromosome Conformation Capture. The specific AP-1 inhibitor nordihydroguaiaretic acid suppresses the calcium-induced increase of PADI3 mRNA levels in keratinocytes. Our findings pave the way to the exploration of deimination control during tumorigenesis and wound healing, two conditions for which AP-1 factors are critical, and disclose that long-range transcription control has a role in the regulation of the gene PADI3. Since invalidation of distant regulators causes a variety of human diseases, PIE results to be a plausible candidate in association studies on deimination-related disorders or atopic disease.

Long-range enhancer differentially regulated by c-Jun and JunD controls peptidylarginine deiminase-3 gene in keratinocytes

Long-range cis elements are critical regulators of transcription, particularly for clustered paralogous genes. Such are the five PADI genes in 1p35-36 encoding peptidylarginine deiminases, which catalyze deimination, a Ca2+-dependent post-translational modification. Deimination has been implicated in the pathophysiology of severe human diseases such as multiple sclerosis and rheumatoid arthritis. The PADI genes present different expression patterns. PADI1-3 are expressed in the epidermis, with increased expression levels in the most differentiated keratinocytes. Previous studies on PADI proximal promoters failed to explain such specificity of expression. We identified a conserved intergenic sequence in the PADI locus (IG1), which may play a role in PADI transcriptional regulation. In this work, we identified two DNase I.hypersensitive sites located in IG1, PAD intergenic enhancer segment 1 (PIE-S1) and PIE-S2, which act in synergy as a bipartite enhancer of the PADI3 and probably PADI1 promoters in normal human epidermal keratinocytes differentiated by a high-calcium-containing medium (1.5 mM). PIE-S1 and PIE-S2 present all the hallmarks of transcriptional enhancers: orientation-independence, copy-number dependence and cell-type specificity. PIE-S1 and PIE-S2 comprise conserved putative binding sites for MIBP1/RFX1 and activator protein 1, respectively. Deletion mutant screening revealed that these sites are crucial for the enhancer activity. Furthermore, chromatin immunoprecipitation assays evidenced differential binding of JunD or c-Jun on the activator protein 1 site depending on the cell differentiation state. Our results reveal the molecular bases of the expression specificity of PADI1 and PADI3 during keratinocyte differentiation through a long-range enhancer and support a model of PADI gene regulation depending on c-Jun-JunD competition.

BRG1-mediated chromatin remodeling regulates differentiation and gene expression of T helper cells

During T helper cell differentiation, distinct programs of gene expression play a key role in defining the immune response to an environmental challenge. How chromatin remodeling events at the associated cytokine loci control differentiation is not known. We found that the ATP-dependent remodeling enzyme subunit BRG1 was required for T helper 2 (Th2) differentiation and Th2 cytokine transcription. BRG1 binding to cytokine genes was regulated by the extent of differentiation, the extent of activation, and cell fate. BRG1 was required for some features of the chromatin structure in target genes (DNase I hypersensitivity and histone acetylation), suggesting that BRG1 remodeling activity was directly responsible for changes in gene expression. NFAT and STAT6 activity were required for BRG1 recruitment to the Th2 locus control region, and STAT6 associated with BRG1 in a differentiation-inducible manner, suggesting direct recruitment of BRG1 to the bound loci. Together, these findings suggest BRG1 interprets differentiation signals and plays a causal role in gene regulation, chromatin structure, and cell fate.

BEST1 expression in the retinal pigment epithelium is modulated by OTX family members

A number of genes preferentially expressed in the retinal pigment epithelium (RPE) are associated with retinal degenerative disease. One of these, BEST1, encodes bestrophin-1, a protein that when mutated causes Best macular dystrophy. As a model for RPE gene regulation, we have been studying the mechanisms that control BEST1 expression, and recently demonstrated that members of the MITF-TFE family modulate BEST1 transcription. The human BEST1 upstream region from -154 to +38 bp is sufficient to direct expression in the RPE, and positive-regulatory elements exist between -154 and -104 bp. Here, we show that the -154 to -104 bp region is necessary for RPE expression in transgenic mice and contains a predicted OTX-binding site (Site 1). Since another non-canonical OTX site (Site 2) is located nearby, we tested the function of these sites using BEST1 promoter/luciferase constructs by in vivo electroporation and found that mutation of both sites reduces promoter activity. Three OTX family proteins - OTX1, OTX2 and CRX - bound to both Sites 1 and 2 in vitro, and all of them increased BEST1 promoter activity. Surprisingly, we found that human and bovine RPE expressed not only OTX2 but also CRX, the CRX genomic region in bovine RPE was hypersensitive to DNase I, consistent with active transcription, and that both OTX2 and CRX bound to the BEST1 proximal promoter in vivo. These results demonstrate for the first time CRX expression in the RPE, and suggest that OTX2 and CRX may act as positive modulators of the BEST1 promoter in the RPE.

Acetylation of EKLF is essential for epigenetic modification and transcriptional activation of the beta-globin locus

Posttranslational modifications of transcription factors provide alternate protein interaction platforms that lead to varied downstream effects. We have investigated how the acetylation of EKLF plays a role in its ability to alter the beta-like globin locus chromatin structure and activate transcription of the adult beta-globin gene. By establishing an EKLF-null erythroid line whose closed beta-locus chromatin structure and silent beta-globin gene status can be rescued by retroviral infection of EKLF, we demonstrate the importance of EKLF acetylation at lysine 288 in the recruitment of CBP to the locus, modification of histone H3, occupancy by EKLF, opening of the chromatin structure, and transcription of adult beta-globin. We also find that EKLF helps to coordinate this process by the specific association of its zinc finger domain with the histone H3 amino terminus. Although EKLF interacts equally well with H3.1 and H3.3, we find that only H3.3 is enriched at the adult beta-globin promoter. These data emphasize the critical nature of lysine acetylation in transcription factor activity and enable us to propose a model of how modified EKLF integrates coactivators, chromatin remodelers, and nucleosomal components to alter epigenetic chromatin structure and stimulate transcription.

The role of transcriptional activator GATA-1 at human beta-globin HS2

GATA-1 is an erythroid activator that binds β-globin gene promoters and DNase I hypersensitive sites (HSs) of the β-globin locus control region (LCR). We investigated the direct role of GATA-1 interaction at the LCR HS2 enhancer by mutating its binding sites within minichromosomes in erythroid cells. Loss of GATA-1 in HS2 did not compromise interaction of NF-E2, a second activator that binds to HS2, nor was DNase I hypersensitivity at HS2 or the promoter of a linked ε-globin gene altered. Reduction of NF-E2 using RNAi confirmed the overall importance of this activator in establishing LCR HSs. However, recruitment of the histone acetyltransferase CBP and RNA pol II to HS2 was diminished by GATA-1 loss. Transcription of ε-globin was severely compromised with loss of RNA pol II from the transcription start site and reduction of H3 acetylation and H3K4 di- and tri-methylation in coding sequences. In contrast, widespread detection of H3K4 mono-methylation was unaffected by loss of GATA-1 in HS2. These results support the idea that GATA-1 interaction in HS2 has a prominent and direct role in co-activator and pol II recruitment conferring active histone tail modifications and transcription activation to a target gene but that it does not, by itself, play a major role in establishing DNase I hypersensitivity.