Henriette O'Geen

Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.

Discovering Hematopoietic Mechanisms Through Genome-Wide Analysis of GATA Factor Chromatin Occupancy

GATA factors interact with simple DNA motifs (WGATAR) to regulate critical processes, including hematopoiesis, but very few WGATAR motifs are occupied in genomes. Given the rudimentary knowledge of mechanisms underlying this restriction and how GATA factors establish genetic networks, we used ChIP-seq to define GATA-1 and GATA-2 occupancy genome-wide in erythroid cells. Coupled with genetic complementation analysis and transcriptional profiling, these studies revealed a rich collection of targets containing a characteristic binding motif of greater complexity than WGATAR. GATA factors occupied loci encoding multiple components of the Scl/TAL1 complex, a master regulator of hematopoiesis and leukemogenic target. Mechanistic analyses provided evidence for crossregulatory and autoregulatory interactions among components of this complex, including GATA-2 induction of the hematopoietic corepressor ETO-2 and an ETO-2-negative autoregulatory loop. These results establish fundamental principles underlying GATA factor mechanisms in chromatin and illustrate a complex network of considerable importance for the control of hematopoiesis.

Comparison of sample preparation methods for ChIP-chip assays

A single chromatin immunoprecipitation (ChIP) sample does not provide enough DNA for hybridization to a genomic tiling array. A commonly used technique for amplifying the DNA obtained from ChIP assays is ligation-mediated PCR (LM-PCR). However; using this amplification method, we could not identify Oct4 binding sites on genomic tiling arrays representing 1% of the human genome (ENCODE arrays). In contrast, hybridization of a pool of 10 ChIP samples to the arrays produced reproducible binding patterns and low background signals. However the pooling method would greatly increase the number of ChIP reactions needed to analyze the entire human genome. Therefore, we have adapted the GenomePlex whole genome amplification (WGA) method for use in ChIP-chip assays; detailed ChIP and amplification protocols used for these analyses are provided as supplementary material. When applied to ENCODE arrays, the products prepared using this new method resulted in an Oct4 binding pattern similar to that from the pooled Oct4 ChIP samples. Importantly, the signal-to-noise ratio using the GenomePlex WGA method is superior to the LM-PCR amplification method.

Sole-Search: an integrated analysis program for peak detection and functional annotation using ChIP-seq data

Next-generation sequencing is revolutionizing the identification of transcription factor binding sites throughout the human genome. However, the bioinformatics analysis of large datasets collected using chromatin immunoprecipitation and high-throughput sequencing is often a roadblock that impedes researchers in their attempts to gain biological insights from their experiments. We have developed integrated peak-calling and analysis software (Sole-Search) which is available through a user-friendly interface and (i) converts raw data into a format for visualization on a genome browser, (ii) outputs ranked peak locations using a statistically based method that overcomes the significant problem of false positives, (iii) identifies the gene nearest to each peak, (iv) classifies the location of each peak relative to gene structure, (v) provides information such as the number of binding sites per chromosome and per gene and (vi) allows the user to determine overlap between two different experiments. In addition, the program performs an analysis of amplified and deleted regions of the input genome. This software is web-based and automated, allowing easy and immediate access to all investigators. We demonstrate the utility of our software by collecting, analyzing and comparing ChIP-seq data for six different human transcription factors/cell line combinations.

DNA hypomethylation within specific transposable element families associates with tissue-specific enhancer landscape

Transposable element (TE)-derived sequences comprise half of the human genome and DNA methylome and are presumed to be densely methylated and inactive. Examination of genome-wide DNA methylation status within 928 TE subfamilies in human embryonic and adult tissues identified unexpected tissue-specific and subfamily-specific hypomethylation signatures. Genes proximal to tissue-specific hypomethylated TE sequences were enriched for functions important for the relevant tissue type, and their expression correlated strongly with hypomethylation within the TEs. When hypomethylated, these TE sequences gained tissue-specific enhancer marks, including monomethylation of histone H3 at lysine 4 (H3K4me1) and occupancy by p300, and a majority exhibited enhancer activity in reporter gene assays. Many such TEs also harbored binding sites for transcription factors that are important for tissue-specific functions and showed evidence of evolutionary selection. These data suggest that sequences derived from TEs may be responsible for wiring tissue type–specific regulatory networks and may have acquired tissue-specific epigenetic regulation.

ZNF274 recruits the histone methyltransferase SETDB1 to the 3' ends of ZNF genes.

Only a small percentage of human transcription factors (e.g. those associated with a specific differentiation program) are expressed in a given cell type. Thus, cell fate is mainly determined by cell type-specific silencing of transcription factors that drive different cellular lineages. Several histone modifications have been associated with gene silencing, including H3K27me3 and H3K9me3. We have previously shown that genes for the two largest classes of mammalian transcription factors are marked by distinct histone modifications; homeobox genes are marked by H3K27me3 and zinc finger genes are marked by H3K9me3. Several histone methyltransferases (e.g. G9a and SETDB1) may be involved in mediating the H3K9me3 silencing mark. We have used ChIP-chip and ChIP-seq to demonstrate that SETDB1, but not G9a, is associated with regions of the genome enriched for H3K9me3. One current model is that SETDB1 is recruited to specific genomic locations via interaction with the corepressor TRIM28 (KAP1), which is in turn recruited to the genome via interaction with zinc finger transcription factors that contain a Kruppel-associated box (KRAB) domain. However, specific KRAB-ZNFs that recruit TRIM28 (KAP1) and SETDB1 to the genome have not been identified. We now show that ZNF274 (a KRAB-ZNF that contains 5 C2H2 zinc finger domains), can interact with KAP1 both in vivo and in vitro and, using ChIP-seq, we show that ZNF274 binding sites co-localize with SETDB1, KAP1, and H3K9me3 at the 3′ ends of zinc finger genes. Knockdown of ZNF274 with siRNAs reduced the levels of KAP1 and SETDB1 recruitment to the binding sites. These studies provide the first identification of a KRAB domain-containing ZNF that is involved in recruitment of the KAP1 and SETDB1 to specific regions of the human genome.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

A genome-wide analysis of Cas9 binding specificity using ChIP-seq and targeted sequence capture

Clustered regularly interspaced short palindromic repeat (CRISPR) RNA-guided nucleases have gathered considerable excitement as a tool for genome engineering. However, questions remain about the specificity of target site recognition. Cleavage specificity is typically evaluated by low throughput assays (T7 endonuclease I assay, target amplification followed by high-throughput sequencing), which are limited to a subset of potential off-target sites. Here, we used ChIP-seq to examine genome-wide CRISPR binding specificity at gRNA-specific and gRNA-independent sites for two guide RNAs. RNA-guided Cas9 binding was highly specific to the target site while off-target binding occurred at much lower intensities. Cas9-bound regions were highly enriched in NGG sites, a sequence required for target site recognition by Streptococcus pyogenes Cas9. To determine the relationship between Cas9 binding and endonuclease activity, we applied targeted sequence capture, which allowed us to survey 1200 genomic loci simultaneously including potential off-target sites identified by ChIP-seq and by computational prediction. A high frequency of indels was observed at both target sites and one off-target site, while no cleavage activity could be detected at other ChIP-bound regions. Our results confirm the high-specificity of CRISPR endonucleases and demonstrate that sequence capture can be used as a high-throughput genome-wide approach to identify off-target activity.

Autophagy driven by a master regulator of hematopoiesis.

ABSTRACT Developmental and homeostatic remodeling of cellular organelles is mediated by a complex process termed autophagy. The cohort of proteins that constitute the autophagy machinery functions in a multistep biochemical pathway. Though components of the autophagy machinery are broadly expressed, autophagy can occur in specialized cellular contexts, and mechanisms underlying cell-type-specific autophagy are poorly understood. We demonstrate that the master regulator of hematopoiesis, GATA-1, directly activates transcription of genes encoding the essential autophagy component microtubule-associated protein 1 light chain 3B (LC3B) and its homologs (MAP1LC3A, GABARAP, GABARAPL1, and GATE-16). In addition, GATA-1 directly activates genes involved in the biogenesis/function of lysosomes, which mediate autophagic protein turnover. We demonstrate that GATA-1 utilizes the forkhead protein FoxO3 to activate select autophagy genes. GATA-1-dependent LC3B induction is tightly coupled to accumulation of the active form of LC3B and autophagosomes, which mediate mitochondrial clearance as a critical step in erythropoiesis. These results illustrate a novel mechanism by which a master regulator of development establishes a genetic network to instigate cell-type-specific autophagy.

Requirement of the C3HC4 zinc RING finger of the Arabidopsis PEX10 for photorespiration and leaf peroxisome contact with chloroplasts

Plant peroxisomes perform multiple vital metabolic processes including lipid mobilization in oil-storing seeds, photorespiration, and hormone biosynthesis. Peroxisome biogenesis requires the function of peroxin (PEX) proteins, including PEX10, a C3HC4 Zn RING finger peroxisomal membrane protein. Loss of function of PEX10 causes embryo lethality at the heart stage. We investigated the function of PEX10 with conditional sublethal mutants. Four T-DNA insertion lines expressing pex10 with a dysfunctional RING finger were created in an Arabidopsis WT background (ΔZn plants). They could be normalized by growth in an atmosphere of high CO2 partial pressure, indicating a defect in photorespiration. β-Oxidation in mutant glyoxysomes was not affected. However, an abnormal accumulation of the photorespiratory metabolite glyoxylate, a lowered content of carotenoids and chlorophyll a and b, and a decreased quantum yield of photosystem II were detected under normal atmosphere, suggesting impaired leaf peroxisomes. Light and transmission electron microscopy demonstrated leaf peroxisomes of the ΔZn plants to be more numerous, multilobed, clustered, and not appressed to the chloroplast envelope as in WT. We suggest that inactivation of the RING finger domain in PEX10 has eliminated protein interaction required for attachment of peroxisomes to chloroplasts and movement of metabolites between peroxisomes and chloroplasts.