Proteogenomic characterization and mapping of nucleosomes decoded by Brd and HP1 proteins
- Gary LeRoy†1,
- Iouri Chepelev†2,
- Peter A DiMaggio1,
- Mario A Blanco1,
- Barry M Zee1, 3,
- Keji Zhao2 and
- Benjamin A Garcia1, 3, 4, 5Email author
© LeRoy et al.; licensee BioMed Central Ltd. 2012
Received: 7 December 2011
Accepted: 16 August 2012
Published: 16 August 2012
Histone post-translational modifications (PTMs) constitute a branch of epigenetic mechanisms that can control the expression of eukaryotic genes in a heritable manner. Recent studies have identified several PTM-binding proteins containing diverse specialized domains whose recognition of specific PTM sites leads to gene activation or repression. Here, we present a high-throughput proteogenomic platform designed to characterize the nucleosomal make-up of chromatin enriched with a set of histone PTM binding proteins known as histone PTM readers. We support our findings with gene expression data correlating to PTM distribution.
We isolated human mononucleosomes bound by the bromodomain-containing proteins Brd2, Brd3 and Brd4, and by the chromodomain-containing heterochromatin proteins HP1β and HP1α. Histone PTMs were quantified by mass spectrometry (ChIP-qMS), and their associated DNAs were mapped using deep sequencing. Our results reveal that Brd- and HP1-bound nucleosomes are enriched in histone PTMs consistent with actively transcribed euchromatin and silent heterochromatin, respectively. Data collected using RNA-Seq show that Brd-bound sites correlate with highly expressed genes. In particular, Brd3 and Brd4 are most enriched on nucleosomes located within HOX gene clusters, whose expression is reduced upon Brd4 depletion by short hairpin RNA.
Proteogenomic mapping of histone PTM readers, alongside the characterization of their local chromatin environments and transcriptional information, should prove useful for determining how histone PTMs are bound by these readers and how they contribute to distinct transcriptional states.
A cell's transcriptional program is governed not only by cis-acting DNA sequences, but also by chromatin structure, DNA methylation and histone post-translational modifications (PTMs). Chromatin modifications and the gene expression patterns that accompany them are maintained when a cell divides and are thus known as 'epigenetic' [1, 2]. Histone lysine methylation is associated with both transcribed and non-transcribed chromatin. For example, H3K4me3 is found in the 5' region of active genes, whereas H3K36me2 and H3K36me3 are enriched within the body and the 3' end of active genes, respectively [3, 4]. However, H3K9 and H3K27 trimethylation are implicated in both the formation and spreading of constitutive heterochromatin and silencing of euchromatin. In contrast to methylation, histone H3 and H4 acetylation is a general characteristic of transcriptionally active chromatin. Depending on the residue modified, acetylations can be localized to nucleosomes near the 5' promoter region or enriched throughout coding regions of actively transcribed genes. For instance, chromatin immunoprecipitation experiments have shown that H3K27ac and H3K14ac are enriched on only a few nucleosomes at the promoters of active genes, whereas H4 acetylations correlate with large euchromatic regions containing transcribed loci [3, 5].
The 'histone code' hypothesis proposes that specific histone PTMs encode regulatory information that is read by the binding of accessory proteins. These accessory proteins, termed 'readers', bind via specialized histone-PTM-binding domains such as bromodomains, chromodomains and plant homeodomains (PHDs) . The bromodomain binds directly to acetylated lysines in histones and is commonly found in proteins associated with gene activation, such as Brg1, a subunit of the hSWI/SNF remodeling complex, the acetyltransferase hGcn5, and the Brd proteins, which are mammalian homologs of the Drosophila trithorax Fsh1 and Fsh2 proteins . The Brd2, Brd3 and Brd4 proteins are referred to as BET proteins as they contain tandem bromodomains and an extraterminal domain of unknown function . In vitro studies have demonstrated that Brd2 and Brd3 possess nucleosome chaperone activities that allow RNA polymerase II to elongate transcripts through hyperacetylated nucleosomes, implicating their direct role in transcription . Several reports have also implicated Brd regulation of cell cycle progression and inflammatory response [9, 10]. Thus, the unique protein architecture and function of Brd proteins qualify them as promising targets for selective pharmacological inhibitors such as immunosuppressants and anticancer drugs. In fact, recent publications have reported two Brd protein selective small molecule inhibitors (I-BET and JQ-1) that specifically inhibit bromodomain binding to acetylated histones, thereby blocking a lipopolysaccharide-induced cytokine storm and the growth of Brd-dependent tumors [9, 11].
In contrast to bromodomains, chromodomains are methyl-lysine binding domains found in proteins such as the heterochromatic proteins HP1α and HP1β that bind to H3K9me3 . Chromodomains are also found in the polycomb group, CBX proteins, which mediate silencing by packaging specifically methylated nucleosomes into heterochromatin-like clusters . In addition to a chromodomain, HP1α and HP1β contain a chromoshadow domain thought to be involved in interactions with other proteins . The HP1α and HP1β proteins form both homo- and heterodimers, and although these proteins colocalize at many heterochromatic loci, they also bind to distinct loci .
Characterization of PTMs on histone H4 from Brd- and HP1-associated nucleosomes
PTMs observed on the H4 molecules from Brd-bound nucleosomes contained marks consistent with gene activation. Most striking was the enrichment in the overall degree of acetylation. All the potentially acetylated lysines (K5, K8, K12 and K16) occur on the same propionylated tryptic peptide, so we first quantified H4 tail peptides that were unmodified, mono, di-, tri- and tetra-acetylated, irrespective of the residue modified (Figure 2b, table). When compared to whole genome chromatin, all Brd nucleosomes were enriched in H4 acetyl modifications; 1.3- to 2.3-fold for the diacetylated species, 5.4- to 8.3-fold for the triacetylated species and 12.9- to 39.5-fold for the tetracetylated form (Figure 2b, heatmap). Notably, the Brd nucleosomes contained only 15 to 29% unacetylated H4 molecules as opposed to genomic histone H4, which was 50% unacetylated. As nucleosomes contain two copies of each histone, the unacetylated species within these nucleosomes is probably not directly bound by the Brds. On the other hand, HP1α- and HP1β-bound nucleosomes had relatively low levels of H4 acetyls. In fact, approximately 70% of the H4 tails were unmodified on HP1-bound nucleosomes (Figure 2b, table). However, the HP1 bound nucleosomes have moderate levels of H4 acetylation, such as tails with a monoacetyl (approximately 25%). The single most abundant H4 tail acetylation in HEK 293 cells is H4K16, occurring in approximately 41% of all histone H4 molecules (Additional file 2). This particular acetylation mark was only slightly enriched on Brd nucleosomes (1.2- to 1.5-fold) (Figure 2c). We suspect that Brd proteins may not directly read this PTM; rather, its slight enrichment in Brd nucleosomes is a consequence of being associated with other active chromatin marks. Nucleosomes bound by all three Brd proteins were significantly enriched in acetylations at the other sites on histone H4; H4K5ac (4.7- to 6.3-fold), H4K8ac (4.1- to 8.7-fold) and H4K12ac (2.3- to 3.3-fold) (Figure 2c) consistent with X-ray crystallography structural studies [19, 20].
Our in-house software (see Materials and methods) accurately deconvoluted the mixed tandem mass spectra resulting from isobaric acetylation combinations that occur on H4 (that is, K5, K8, K12 and K16). When compared to whole genome chromatin, nucleosomes bound by all three Brd proteins were enriched with the H4 diacetylated forms (K5unK8acK12unK16ac, K5unK8acK12acK16un and K5acK8unK12unK16ac), while only Brd3- and Brd4-bound nucleosomes were enriched in the K5acK8unK12acK16un species (Figure 3b, heatmap; Additional file 3). These results suggest that Brds will tolerate H4 containing K16ac, although it is not required for binding by the Brds. As depicted in Figure 3b, the nucleosomes bound separately by all three Brd family member proteins were also highly enriched with all H4 triacetylated forms when compared to whole-genome chromatin, such as the H4K5unK8acK12acK16ac (4.4- to 7.2-fold), H4K5acK8unK12acK16ac (4.7- to 7.0-fold), H4K5acK8acK12unK16ac (4.3- to 20.3-fold) and especially H4K5acK8acK12acK16un (11.3- to 44.9-fold). t-Tests performed with data from three independent experiments for each ChIP revealed that these fold changes are indeed significant as the P-values were all <1 × 10-2 (Additional file 4). The finding that H4K5acK8acK12acK16un is the most enriched triacetylated H4 code further supports the hypothesis that H4K16 acetylation is not a Brd binding site. The tetracetylated H4 code, H4K5acK8acK12acK16ac, is found on 0.5% of all H4 in whole genome chromatin in HEK 293 cells; however, this form is strikingly enriched in nucleosomes bound by all three Brds - Brd2 (20.8%), Brd3 (6.8%) and Brd4 (9.3%) - and nearly undetectable in HP1-bound nucleosomes. Moreover, HP1 nucleosomes were depleted of most combinatorial forms of H4 acetylations compared to whole genome chromatin (Figure 3b).
H4K20me2 is the most abundant mark on human histone H4, found on approximately 58% of all H4 in HEK 293 cells (Additional file 2), consistent with its incorporation in heterochromatin as ≥70% of mammalian genomes are thought to be heterochromatic. H4K20me2 was enriched in HP1-bound nucleosomes and depleted in Brd nucleosomes (Figure 3b). H4K20me2 was especially depleted in Brd4-bound nucleosomes (approximately 8.5% abundance). Although this mark is extremely high in HP1-bound nucleosomes, approximately 80% for HP1β and approximately 73% for HP1α, we lack direct evidence to indicate that the HP1 chromodomains are reading this mark (Additional file 2). Similarly to H4K16ac, H4K20me2 may be part of a more complex heterochromatic histone code recognized by an HP1 complex. Additionally, as HP1 is known to bind H3K9me3, it is reasonable to speculate that the genomic location of this heterochromatic mark overlaps with H4K20me2.
Characterization of PTMs on histone H3 from Brd- and HP1-associated nucleosomes
As indicated in Figures 2c and 3a, Brd-bound nucleosomes all have elevated levels of H3K4me1, a gene activation mark. HP1-bound nucleosomes had relatively low levels of this mark. Moreover, H3K4me2 and H3K4me3 were also elevated in the Brd-bound nucleosomes and depleted in HP1-bound nucleosomes. When compared to whole genome chromatin, nucleosomes bound by Brd2 were the most enriched in these modifications; 12.4-fold (P-value = 0.01) for H3K4me2 and 23.6-fold (P-value = 0.0009) for H3K4me3 (Figures 2c and 3a; Additional file 5). These particular PTMs occur on the nucleosomes surrounding the promoters of actively transcribed genes and are at very low levels genome-wide (≤1%).
The chromodomains of HP1α/β bind directly to silencing marks H3K9me2 and H3K9me3 [12, 21]. HP1-bound nucleosomes were enriched with both these K9 methylations (Figure 2c). For these modifications, a 1.5- to 2-fold increase is significant, as these modifications are extremely abundant in whole-genome chromatin, 31% and 21%, respectively (Additional file 2). We also found that Brd-bound nucleosomes were low in H3K9me3, but have moderate levels of H3K9me2 (Figure 2c).
Another modification, H3K23me1, is also found at low abundance (0.2%), but was highly enriched on both HP1s (Figure 2c; Additional file 2). This result is not surprising as a recent study employing a histone peptide methyl-lysine array found that HP1β binds H3K23me1 in addition to H3K9me3  and colocalizes with H3K9me3 and HP1β at heterochromatic loci. We believe that a subset of HP1 proteins may bind H3K23me1 via its chromodomain independent of H3K9me3 within these loci. H3K18me1 is another modification that is of very low relative abundance in whole genome chromatin (0.3%). Here we find that Brd-bound nucleosomes were highly enriched with this modification; Brd2 (5.23%), Brd3 (4.17%) and Brd4 (3.10%) (Additional file 2). H3K18 can also be acetylated by the CBP/p300 acetyltransferase, and H3K18ac localizes with H3K4me3 by ChIP at the promoters of some actively transcribed genes [3, 23–25]. As depicted in Figure 2c, nucleosomes bound by Brds were enriched in H3K18ac (4.2- to 6.9-fold; P-values <0.001; Additional file 5) compared to bulk chromatin. This result is consistent with reports that H3K18ac localizes with H3K4me3, as Brd nucleosomes were also enriched with H3K4me3. The nucleosomes bound by the Brds also exhibited high levels of acetylation at H3K9ac, H3K14ac, and H3K23ac, whereas acetylation of these residues was reduced on HP1-associated nucleosomes (Figure 2c).
Two other abundant gene repressive PTMs are H3K27me3 and H3K27me2. Nucleosomes bound by both HP1 proteins were modestly enriched with these modifications (Figure 2c). H3K27me3 and H3K27me2 are thought to demarcate the promoter regions and bodies, respectively, of particular repressed genes. These methylations are added by the polycomb group proteins Ezh1 and Ezh2 and bound by other polycomb group chromodomain-containing CBX proteins [26, 27]. Interestingly, our observation that HP1-bound nucleosomes were enriched with H3K27me2/3 suggests that there may be some overlap amongst heterochromatin bound by the HP1 and polycomb CBX chromodomain proteins. Future ChIP-seq experiments with the polycomb CBX proteins will be required to determine their genomic overlap with the HP1 proteins.
As some H3 PTMs can be observed on the same peptides, a few combinations of PTMs were also characterized (Additional file 6). We believe that there is useful combinatorial information in tandem PTMs; for example, we previously reported that the histone H3K9 euchromatic dimethytransferase G9a preferentially methylates H3 peptides that contain K14 acetylation . This modification, H3K9me2K14ac, shares similar characteristics to the 'bivalent domain' (H3K4me3K27me3) found on silent genes poised for activation in pluripotent stem cells, containing PTMs with opposing epigenetic functions . As mentioned previously, Brd-bound nucleosomes contained some H3K9me2, but the bulk of this modification occurred on peptides that were acetylated on the adjacent lysine (H3K9me2K14ac), the tandem marks effected by G9a (Figure 3a). This result is consistent with the observation that G9a localizes primarily to euchromatin, and that H3K9me2K14ac may also act as a combinatorial switch, similar to the bivalent domain. In contrast, H3K9me3 was depleted in Brd nucleosomes even when H3K14 was acetylated. On the other hand, there was substantial amounts of this combinatorial modification, H3K9me3K14ac, on nucleosomes bound by HP1α (13.7%) and HP1β (16.3%), which occurs on only 9.8% of the nucleosomes in the genome (Additional file 6).
Brd-bound nucleosomes also contain histone H3 enriched with other combinatorial PTMs consistent with gene activation, such as H3K9acK14ac (2.9- to 4.6-fold; P-values <0.038) and H3K18acK23ac (5.0- to 8.2-fold; P-values <0.0008; Additional file 7) and are depleted in marks that correspond with gene repression, such as H3K9me3K14un (Figure 3c). H3K36me3 is another modification associated with the bodies of actively transcribed genes. We found that Brd4-bound nucleosomes were enriched with this modification (2.2-fold), but only when the H3K27 was monomethylated (H3K27me1K36me3) (Figures 2c and 3a).
Chromosomal locations of Brd- and HP1-bound nucleosomes (ChIP-Seq)
At the gene level (Figure 4b), approximately 65% of Brd-binding sites were found on genes (promoter plus gene body). The remaining approximately 35% was found on intergenic regions, very often in regions flanking genes. The promoter regions as defined in this study are the +2 kb regions upstream of the transcription start sites. We observed that 9 to 13% of the Brd-binding sites were found within such promoter regions. On the other hand, HP1α and β binding sites were predominately found in intergenic regions (73% and 66%, respectively) and rarely on promoter regions of genes (1 to 2%). Complete lists of promoters bound by Brd and HP1 ChIPs can be found in Additional file 8. As shown in Figure 4c-e, patterns of Brd binding indicate that these proteins are predominately found on high and moderately expressed genes and rarely on transcriptionally silent genes. Genes were ranked based on their expression levels as determined by RNA-Seq, where high expression is defined as the top third of the genes, the middle third as moderate and the bottom third as silent. As predicted, the vast majority of genes bound by the HP1s were classified as either silent or having moderate expression (Figure 4f,g). A heatmap depicting the DNA consensus motifs enriched in the ChIPs can be found in Additional file 9 with an accompanying list of motifs (Additional file 10). We have also included lists of Gene Ontology terms enriched in the Brd ChIPs (Additional file 11).
Brd-bound nucleosomes are associated with the promoters of HOX genes
HP1 bound nucleosomes are associated with the promoters of ZNF genes on chromosome 19
As stated above, the HP1s were found at the promoter regions of 1 to 2% of the genes they bind and most of these are transcriptionally silent (Figure 4b,f,g). Moreover, neither HP1 protein was found on any HOX gene promoters. We also did not observe an equal distribution amongst the chromosomes with genes whose promoters were bound by the HP1 proteins. For both HP1 proteins, approximately 30% were found on chromosome 19 and 20% were found on the × chromosome. As shown in Figure 5b, when we only considered the top 35 gene promoters bound by HP1 proteins, HP1α and HP1β bound 18 and 15 ZNF gene promoters, respectively, from chromosome 19. ZNF genes encode proteins that belong to a large superfamily of putative transcription factors (approximately 800 members in humans), many of which are thought to function as transcriptional repressors . There are six clusters of ZNF genes found on chromosome 19. These clusters all span from 1 to 3 Mb of DNA each and are located along chromosome 19 as follows: cluster 1, 10 to 12 Mb; cluster 2, 19 to 23 Mb; cluster 3, 40 to 43 Mb; cluster 4, 49 to 50 Mb; cluster 5, 57 to 59 Mb; and cluster 6, 61 to 63 Mb (Figure 5c). When considering the top 100 bound promoters, we observed that HP1α bound to 37 gene promoters (23 known ZNF genes), and HP1β to 55 gene promoters (27 known ZNF genes), all located within these clusters on chromosome 19. However, we did not see complete overlap amongst the promoter binding of both HP1 proteins. While we found extensive binding of HP1β to the promoters of genes in all six clusters, HP1α was only found predominantly in clusters 3, 4 and 6, with very little binding in the other three clusters. Interestingly, the promoters of genes found in cluster 5 are extensively bound by HP1β and almost completely devoid of HP1α binding. Within this cluster, HP1β was bound to the promoters of 59 genes, 16 of which are ZNF genes, as well as a subcluster of 20 microRNA genes in the distal half of this cluster. These data are consistent with a study in 2006 that used a DNA adenine methyltransferase identification (DamID) method coupled with cDNA oligo microarrays to identify the target genes bound by HP1β . In agreement with these findings, ChIP-seq experiments with antibodies directed against H3K9me3 have found distinct islands of this modification map to the ZNF clusters .
Brd4 is a positive regulator of the genes it binds
HP1β is a negative regulator of the genes it binds
In order to determine the direct effects of HP1 binding to nucleosomes, we preformed gene expression arrays on HEK293 cells depleted of HP1β by shRNA. We chose to knockdown HP1β rather than HP1α because it was enriched on more promoters, including most of the HP1α-bound promoters. Our custom shRNA knocked down HP1β mRNA by 3.5-fold in HEK293 cells as determined by RT-PCR. The knockdown was also verified at the protein level by western blot (Additional file 13). In contrast to the effects of Brd knockdown, knockdown of HP1β resulted in at least two-fold increases expression in two-thirds of the top 100 HP1β promoter-bound genes (Figure 6a, left heatmap). These results suggest that the binding of HP1β to histone methylated nucleosomes is required for effective silencing of associated genes. We also looked at the expression of the ZNF genes found on chromosome 19 whose promoters were enriched with HP1β-bound nucleosomes (Figure 6b, left heatmap). Surprisingly, we found that the expression of only one-third of these genes significantly increased, whereas two-thirds remained the same or decreased. This result is in agreement with the aforementioned study that mapped HP1β to the ZNF clusters on chromosome 19 with the DamID method . They also reported that individual depletion of HP1α/β did not lead to an appreciable change in the expression of most of the ZNF clusters. We hypothesize that this is due to the redundancy of the three HP1s (HP1α, β and γ) and attempts to circumvent this by depleting more than one at a time have been unsuccessful. Also, removing the heterochromatic binding proteins alone may not be enough to fully turn on repressed genes, as there are other steps required for gene activation, such as recruitment of transcription factors or gene activating histone-modifying enzymes or other binding proteins. This hypothesis would also be consistent with our finding that the expression of the majority of genes whose promoter regions were bound by HP1β did not change significantly when HP1β was depleted.
Brd-bound nucleosomes are associated throughout the body of transcribed HOX genes
HP1-bound nucleosomes are associated with the bodies of transcriptionally silent ZNF genes
We next analyzed the binding patterns of HP1 proteins on the ZNF560 and ZNF599 genes (Figure 7e,f). Both genes were determined to be silent in control HEK293 cells by RNA-Seq analysis (Figure 7e,f, bottom panels). ZNF560 and ZNF599 are found in the HP1-binding clusters 3 and 6 on chromosome 19, which are significantly bound by both HP1α and HP1β. We chose these ZNF genes for further analysis as their expression was the most upregulated of all HP1-bound ZNF genes upon depletion of HP1β by shRNA (Figure 6b). As depicted in Figure 7e,f, both HP1s bound nucleosomes throughout the bodies and promoter regions of ZNF599 and ZNF560, with HP1β being slightly more enriched within the promoter and 5' region of ZNF599. These results, along with our array analysis, suggest that the HP1 proteins directly repress the expression of these genes. Moreover, nucleosomes found within both of these genes are relatively devoid of binding by all three Brds. Interestingly, we did observe a short cluster of Brd-bound nucleosomes at the end of the ZNF599 gene and continuing in the 3' direction on chromosome 19; however, this result is rather perplexing as there are no actively transcribed genes in this region.
Since the completion of the Human Genome Project in 2000, there have been several calls for a similar project to catalog human epigenomes. In 2008 the Alliance for the Human Epigenome and Disease (AHEAD) outlined the need for a project to provide high-resolution epigenome maps and suggested guidelines for the first pass of a human epigenome reference . At the heart of these initiatives is the use of ChIPs with antibodies that recognize specific PTMs. Such studies are informative; however, ChIPs with modification-specific antibodies will only provide epigenomic maps with little knowledge on how such information is translated. Additionally, site-specific antibodies often suffer from cross-reactivity and epitope occlusion. Our less biased proteomic-genomic approaches complement such ongoing studies and shed light on how histone PTMs are translated into transcriptional outcomes.
Two different models have been proposed for the function of Brd proteins in transcription. In one model Brd4 recruits the transcription elongation factor pTefb to the transcription complex at the transcription start site, which then phosphorylates the CTD of RNA polymerase II, leading to more processive transcript elongation . In the second model, Brd proteins are recruited to acetylated nucleosomes via their bromodomains and facilitate the passage of RNA polymerase II through nucleosomes by their intrinsic chaperone activity . Our ChIP-seq results support both models as we have observed an enrichment of the Brd proteins at both the promoter (approximately 12%) and body (approximately 52%) of active genes. These results suggest an equal distribution of Brd proteins on promoter regions and gene bodies, as we have defined promoters regions as 2 kb upstream of the transcription start site and the average gene body is 10 to 20 kb. Moreover, the Brd-bound nucleosomes were enriched with patterns of histone PTMs found in both the promoter region and body of actively transcribed genes. ChIP-seq experiments with specific antibodies have found that the PTMs H3K9ac and H3K4me3 are associated with only a few nucleosomes at active gene promoters, whereas H3K79me1 and H4K20me1 are exclusively associated with the bodies of active genes [3, 4]. When compared to whole-genome chromatin, the Brd-bound nucleosomes were enriched in all of these modifications, further suggesting that they represent a mixture of promoter and gene body nucleosomes. The Brd-bound nucleosomes were also enriched with the histone PTMs H3K14ac (approximately 55%), H3K23ac (approximately 50%), H4K12ac (approximately 50%) and H4K16ac (approximately 55%). These modifications were previously shown to be enriched on both the promoters and bodies of active genes [3, 4]. Our data suggest that Brd nucleosomes are a mixture of the following: promoter bound with the predominant histone PTM patterns H3(K4me3, K9ac, K14ac, K23ac) and H4(K5ac, K12ac, K16ac), and gene body enriched with H3(K4me1, K14ac, K23ac, K79me1) and H4(K8ac, K12ac, K16ac, K20me1).
HP1-bound nucleosomes were most enriched with histones containing the repressive PTMs, H3K9me2/3 (approximately 85% combined) and H4K20me2/3 (approximately 75% combined). In addition to binding the ZNF clusters of chromosome 19, our ChIP-seq results have revealed that the HP1 proteins reside on centromeric heterochromatin and inactive × chromatin, both of which harbor these marks [3, 4, 17]. Also in line with our results, previous ChIP-seq experiments have colocalized these modifications to the ZNF repeat clusters on chromosome 19 .
Silent Hox genes are a defining hallmark of pluripotent stem cells and their expression is tightly regulated during differentiation . The differential expression patterns of the Hox genes are the master regulators that direct cellular differentiation and maintain particular cell phenotypes . In fact, all of the protein coding Hox genes contained bivalent domains, suggesting that they are poised for activation . Moreover, another ChIP-seq study found that nearly half of all bivalent domains were bound by the polycomb PRC1 complex, which contains a chromodomain protein that binds directly to H3K27me3 . Such bivalent domains were defined as containing overlapping PTMs, H3K4me3 and H3K27me3. It is yet unknown what activating or repressive marks on histone H4 are found within these bivalent domains. Future experiments will be necessary to investigate whether H4 acetylations that are bound by the Brd proteins are also present in bivalent domains. If so, it will be interesting to investigate whether Brd proteins are also bound to bivalent domain chromatin concurrent with the PRC1 binding or whether PRC1 complex is released prior to H4 acetylation and Brd binding upon gene activation. In conclusion, the proteomic and genomic techniques presented here can be used to probe the local chromatin environment bound by any chromatin-associated protein. In particular, we have examined the chromatin from genomic regions where bromodomain- or chromodomain-containing proteins are found, and quantified the histone PTMs present. These data can be used to generate hypotheses on how chromatin-associated proteins (or histone code readers) could influence the transcription of specific genes, hence yielding glimpses into how histone codes could be translated. Such experimental approaches could also be applied to any of the histone code readers in systems of biomedical importance, such as during stem cell differentiation or the metastatic progression of cancers.
Materials and methods
Mononucleosomes were purified after (micrococcal nuclease) Mnase digestion with FLAG-Brd or FLAG-HP1 proteins as described in  from HEK293 cells. FLAG-HP1 vectors were described in . All proteins were generated by transient expression as previously described [8, 44], and expression assessed by western blots as shown in Additional file 1. The quality of the ChIPs was analyzed by 17% SDS-PAGE and Coommassie blue staining. The histone protein ChIP inputs and elutions were propionyl derivatized with either 'light' d 0 - or 'heavy' d 10 -propionic anhydride, trypsinized and prepared for mass spectrometry as previously described in  with minor adjustments. 2-Propanol was substituted for methanol during the derivatization reaction as we have found that methanol can modify peptides at low pH.
Liquid chromatography-mass spectrometry-based proteomics
Histone peptides were separated by reverse phase nanospray liquid chromatography on C18 resin with an Agilent 1200 series HPLC. Mass spectrometry was performed on a LTQ-Orbitrap mass spectrometer (ThermoFisher Scientific Inc. Waltham, MA, USA as previously described .
Mass spectrometry data analysis
The LC-MS/MS data sets were analyzed using in-house software. Briefly, the characterization of the histone PTMs was performed using an optimization-based algorithmic framework that simultaneously utilizes mass spectrometry, MS/MS and chromatographic information to identify and quantify the PTMs on the histone peptides. The method automatically considers multiple charge states of all modified peptides and is able to resolve mixed tandem mass spectra resulting from co-eluting and isobaric modified peptides using a superposition method, as previously described . The statistical significance of the relative abundance of histone modifications in various ChIPs as compared to the relative abundance of genomic histone PTMs was calculated using a two-tailed, unpaired Student's t-test, and the associated P-values were generated using Microsoft Excel. Mass spectrometric data were converted from RAW format into mzXML format using ReAdW. The mzXML files were deposited into the Tranche repository and are available for download from the Proteome Commons .
Preparation of DNA for deep sequencing analysis
Nucleosomal DNAs were treated with alkaline phosphatase to remove the 3' phosphates left after MNase digestion and purified by phenol extraction. The DNA fragments (approximately 1 μg) were further prepared using the Solexa® Illumina library kit directly following the manufacturer's protocol. The purified DNA was used directly for cluster generation and sequencing using a Solexa 1G Genome Analyzer (Illumina).
Preparation of RNA for RNA-Seq analysis
The raw sequence data for the RNA-Seq in HEK293 cells were obtained from  (Gene Expression Omnibus sample ID GSM301568). The RNA preparation and sequencing protocols for this RNA-Seq data can be found in .
Deep sequencing analysis
The raw sequence reads were aligned to the hg18 human genome using the Bowtie short read aligner . To avoid possible PCR artifacts, we kept at most one uniquely aligned read at each genomic location per strand. To determine the number of ChIP-Seq sequencing reads in a genomic interval, we shifted the locations of forward and reverse strand tags by 75 bp, half the size of the nucleosomal DNA.
For the chromosomal binding patterns of Brd and HP1 proteins (Figures 4a and 5c), reads were counted in sliding windows of 1 Mb with a step size of 100 kb. The read counts were normalized to the sequencing depth of 5 million non-redundant reads. The signal level in each 1-Mb window was computed as the difference between the normalized read counts in the window in the sample and the control libraries (Nsample - Ncontrol) with the negative values set to zero. A clustering approach described in  was used to identify Brd and HP1 protein-enriched domains ('islands') in the human genome from the ChIP-Seq data. The percentage overlap (Figure 4b) of the Brd and HP1 islands with the promoter, the gene-body and the inter-genic regions were determined using the gene annotations from the Ensembl database.
The sequence read density (per 100-bp window and 5 million non-redundant reads) was determined along the gene body and in the upstream/downstream 2-kb regions in the sample and control libraries. The signal level in each window was computed as the difference between the read densities in the sample and control libraries, with the negative values set to zero. Genes were ranked based on their expression levels estimated as reads per kilobase per million mapped reads (RPKM  and stratified into three expression groups: the 'silent' group (the bottom 30% genes in the ranked list), the 'high' group (the top 30% genes), and the 'medium' group (the remaining genes). The mean signal levels over all genes in each group were determined and are displayed as gene-body plots in Figure 4c-g.
The patterns of the combinatorial binding of Brd and HP1 proteins to the gene promoters and their correlation to the gene expression levels were determined as follows. A gene promoter, defined as a 2-kb region upstream of the transcription start site, is deemed bound by the protein if the island (determined as described above) overlaps the promoter. Genes with a particular combinatorial promoter-binding pattern of the five proteins (such as 'red-red-red-blue-blue' in Figure 5a) were grouped. The groups were then ranked based on the number of their member genes and displayed in the descending order in Figure 5a (except for the most populous group, 'blue-blue-blue-blue-blue', displayed at the bottom). The expression status of each gene is indicated as expressed ('red') if its RPKM value is at least 1 and not expressed ('blue') otherwise. The UCSC genome browser  bedGraph files displayed in Figure 7a-f were generated as follows. For the ChIP-Seq data, the 75-bp-shifted non-redundant reads were counted in 200-bp tiling windows along the genome and the read counts were normalized to a library size of 5 million non-redundant uniquely mapped reads in each library. For the RNA-Seq data, the centers of the reads were counted in 20-bp tiling windows along the human genome. We deposited raw and analyzed ChIP-Seq and microarray data in the Gene Expression Omnibus .
Brd4 and HP1β knockdown cells
HEK293 cells were infected with lentiviral MISSION® (Sigma-Aldrich Corporation, St. Louis, MO, USA) shRNA vectors (pLKO.1-puro) harboring shRNAs to Brd4 (clone TRCN0000062223), HP1β (clone TRCN0000021427) or control (empty pLKO.1-puro vector). After 72 hours the infected cells were selected with puromycin according to the manufacturer's protocols (Sigma Aldrich). Knockdowns of target RNAs were analyzed by qRT-PCR. For qRT-PCR, RNA from cells was purified with an RNeasy kit (Qiagen) and prepared with a cDNA synthesis kit (Invitrogen Carlsbad, CA, USA). qRT-PCR reactions were performed with SYBR FAST qPCR supermix (KAPA Biosystems Woburn, MA, USA).
RNA was extracted and purified with an RNeasy kit (Qiagen) from Brd4, HP1β and control shRNA knockdown HEK293 cells. For each sample 400 ng RNA was linearly amplified and labeled with Cy5-CTP using Low RNA Linear Amplification reagents (Agilent Technologies Santa Clara, CA, USA). Reference RNA from wild-type HEK293 cells was purified, amplified and labeled with Cy3 in the same manner. Equal quantities of the Cy5- and Cy3-labeled cRNAs were mixed and competitively hybridized to human GE 4x44K microarrays (Agilent, G4112F) for 18 h using an Agilent hybridization kit. Slides were washed according to the manufacturer's protocols and scanned using an Agilent G2565BA scanner. Raw image data were extracted and subjected to standard background subtraction and linear and lowess normalization using Agilent Feature extraction software v9.5. Heatmaps were generated using MATLAB® software (MathWorks Natick, MA, USA). More detailed methods are described in Additional file 14.
chromatin immunoprecipitation-high throughput deep sequencing
DNA adenine methyltransferase identification
nano scale liquid chromatography-tandem mass spectrometry
polymerase chain reaction
quantitative mass spectrometry
RNA-high-throughput deep sequencing
reads per kilobase per million mapped reads
short hairpin RNA.
Support from a National Science Foundation (NSF) Early Faculty CAREER award, an NIH Innovator award (DP2OD007447) from the Office of the Director, National Institutes of Health, an NSF grant (CBET-0941143) and the American Society for Mass Spectrometry Research Award to BAG is gratefully acknowledged. PAD is supported by an NIH NRSA F32 Postdoctoral Fellowship. This research was also supported by an NIH/NIGMS P50GM071508 to the Lewis-Sigler Institute at Princeton University (PI D Botstein). BMZ is supported by an NSF graduate fellowship. The work of IC and KZ was supported by the Division of Intramural Research, National Heart, Lung and Blood Institute, National Institutes of Health, USA.
- Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.PubMedView ArticleGoogle Scholar
- Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, Zhang X, Bernstein BE, Nusbaum C, Jaffe DB, Jaenisch R, Lander ES: Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008, 454: 766-770.PubMedPubMed CentralGoogle Scholar
- Wang Z, Zang C, Rosenfeld JA, Schones DE, Barski A, Cuddapah S, Cui K, Roh TY, Peng W, Zhang MQ, Zhao K: Combinatorial patterns of histone acetylations and methylations in the human genome. Nat Genet. 2008, 40: 897-903. 10.1038/ng.154.PubMedPubMed CentralView ArticleGoogle Scholar
- Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129: 823-837. 10.1016/j.cell.2007.05.009.PubMedView ArticleGoogle Scholar
- Hebbes TR, Thorne AW, Crane-Robinson C: A direct link between core histone acetylation and transcriptionally active chromatin. EMBO J. 1988, 7: 1395-1402.PubMedPubMed CentralGoogle Scholar
- Chang YL, King B, Lin SC, Kennison JA, Huang DH: A double-bromodomain protein, FSH-S, activates the homeotic gene ultrabithorax through a critical promoter-proximal region. Mol Cell Biol. 2007, 27: 5486-5498. 10.1128/MCB.00692-07.PubMedPubMed CentralView ArticleGoogle Scholar
- Florence B, Faller DV: You bet-cha: a novel family of transcriptional regulators. Front Biosci. 2001, 6: D1008-1018. 10.2741/Florence.PubMedView ArticleGoogle Scholar
- LeRoy G, Rickards B, Flint SJ: The double bromodomain proteins Brd2 and Brd3 couple histone acetylation to transcription. Mol Cell. 2008, 30: 51-60. 10.1016/j.molcel.2008.01.018.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicodeme E, Jeffrey KL, Schaefer U, Beinke S, Dewell S, Chung CW, Chandwani R, Marazzi I, Wilson P, Coste H, White J, Kirilovsky J, Rice CM, Lora JM, Prinjha RK, Lee K, Tarakhovsky A: Suppression of inflammation by a synthetic histone mimic. Nature. 2010, 468: 1119-1123. 10.1038/nature09589.PubMedView ArticleGoogle Scholar
- Denis GV, Vaziri C, Guo N, Faller DV: RING3 kinase transactivates promoters of cell cycle regulatory genes through E2F. Cell Growth Differ. 2000, 11: 417-424.PubMedPubMed CentralGoogle Scholar
- Filippakopoulos P, Qi J, Picaud S, Shen Y, Smith WB, Fedorov O, Morse EM, Keates T, Hickman TT, Felletar I, Philpott M, Muro S, McKeown MR, Wang Y, Christie AL, West N, Cameron MJ, Schwartz B, Heightman TD, La Thangue N, French CA, Wiest O, Kung AL, Knapp S, Bradner JE: Selective inhibition of BET bromodomains. Nature. 2010, 468: 1067-1073. 10.1038/nature09504.PubMedPubMed CentralView ArticleGoogle Scholar
- Jacobs SA, Khorasanizadeh S: Structure of HP1 chromodomain bound to a lysine 9-methylated histone H3 tail. Science. 2002, 295: 2080-2083. 10.1126/science.1069473.PubMedView ArticleGoogle Scholar
- Levine SS, King IF, Kingston RE: Division of labor in polycomb group repression. Trends Biochem Sci. 2004, 29: 478-485. 10.1016/j.tibs.2004.07.007.PubMedView ArticleGoogle Scholar
- Nielsen AL, Oulad-Abdelghani M, Ortiz JA, Remboutsika E, Chambon P, Losson R: Heterochromatin formation in mammalian cells: interaction between histones and HP1 proteins. Mol Cell. 2001, 7: 729-739. 10.1016/S1097-2765(01)00218-0.PubMedView ArticleGoogle Scholar
- Dialynas GK, Terjung S, Brown JP, Aucott RL, Baron-Luhr B, Singh PB, Georgatos SD: Plasticity of HP1 proteins in mammalian cells. J Cell Sci. 2007, 120: 3415-3424. 10.1242/jcs.012914.PubMedView ArticleGoogle Scholar
- Cammas F, Janoshazi A, Lerouge T, Losson R: Dynamic and selective interactions of the transcriptional corepressor TIF1 beta with the heterochromatin protein HP1 isotypes during cell differentiation. Differentiation. 2007, 75: 627-637. 10.1111/j.1432-0436.2007.00166.x.PubMedView ArticleGoogle Scholar
- Bartova E, Pachernik J, Kozubik A, Kozubek S: Differentiation-specific association of HP1alpha and HP1beta with chromocentres is correlated with clustering of TIF1beta at these sites. Histochem Cell Biol. 2007, 127: 375-388. 10.1007/s00418-006-0259-1.PubMedView ArticleGoogle Scholar
- Plazas-Mayorca MD, Zee BM, Young NL, Fingerman IM, LeRoy G, Briggs SD, Garcia BA: One-pot shotgun quantitative mass spectrometry characterization of histones. J Proteome Res. 2009, 8: 5367-5374. 10.1021/pr900777e.PubMedPubMed CentralView ArticleGoogle Scholar
- Umehara T, Nakamura Y, Jang MK, Nakano K, Tanaka A, Ozato K, Padmanabhan B, Yokoyama S: Structural basis for acetylated histone H4 recognition by the human BRD2 bromodomain. J Biol Chem. 2010, 285: 7610-7618. 10.1074/jbc.M109.062422.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu Y, Wang X, Zhang J, Huang H, Ding B, Wu J, Shi Y: Structural basis and binding properties of the second bromodomain of Brd4 with acetylated histone tails. Biochemistry. 2008, 47: 6403-6417. 10.1021/bi8001659.PubMedView ArticleGoogle Scholar
- Lachner M, O'Carroll D, Rea S, Mechtler K, Jenuwein T: Methylation of histone H3 lysine 9 creates a binding site for HP1 proteins. Nature. 2001, 410: 116-120. 10.1038/35065132.PubMedView ArticleGoogle Scholar
- Liu H, Galka M, Iberg A, Wang Z, Li L, Voss C, Jiang X, Lajoie G, Huang Z, Bedford MT, Li SS: Systematic identification of methyllysine-driven interactions for histone and nonhistone targets. J Proteome Res. 2010, 9: 5827-5836. 10.1021/pr100597b.PubMedView ArticleGoogle Scholar
- Jin Q, Yu LR, Wang L, Zhang Z, Kasper LH, Lee JE, Wang C, Brindle PK, Dent SY, Ge K: Distinct roles of GCN5/PCAF-mediated H3K9ac and CBP/p300-mediated H3K18/27ac in nuclear receptor transactivation. EMBO J. 2011, 30: 249-262. 10.1038/emboj.2010.318.PubMedPubMed CentralView ArticleGoogle Scholar
- van Oevelen CJ, van Teeffelen HA, van Werven FJ, Timmers HT: Snf1p-dependent Spt-Ada-Gcn5-acetyltransferase (SAGA) recruitment and chromatin remodeling activities on the HXT2 and HXT4 promoters. J Biol Chem. 2006, 281: 4523-4531. 10.1074/jbc.M509330200.PubMedView ArticleGoogle Scholar
- Wang D, Xia X, Weiss RE, Refetoff S, Yen PM: Distinct and histone-specific modifications mediate positive versus negative transcriptional regulation of TSHalpha promoter. PLoS One. 2010, 5: e9853-10.1371/journal.pone.0009853.PubMedPubMed CentralView ArticleGoogle Scholar
- Margueron R, Li G, Sarma K, Blais A, Zavadil J, Woodcock CL, Dynlacht BD, Reinberg D: Ezh1 and Ezh2 maintain repressive chromatin through different mechanisms. Mol Cell. 2008, 32: 503-518. 10.1016/j.molcel.2008.11.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Bernstein E, Duncan EM, Masui O, Gil J, Heard E, Allis CD: Mouse polycomb proteins bind differentially to methylated histone H3 and RNA and are enriched in facultative heterochromatin. Mol Cell Biol. 2006, 26: 2560-2569. 10.1128/MCB.26.7.2560-2569.2006.PubMedPubMed CentralView ArticleGoogle Scholar
- Plazas-Mayorca MD, Bloom JS, Zeissler U, Leroy G, Young NL, DiMaggio PA, Krugylak L, Schneider R, Garcia BA: Quantitative proteomics reveals direct and indirect alterations in the histone code following methyltransferase knockdown. Mol Biosyst. 2010, 6: 1719-1729. 10.1039/c003307c.PubMedPubMed CentralView ArticleGoogle Scholar
- Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, Jaenisch R, Wagschal A, Feil R, Schreiber SL, Lander ES: A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006, 125: 315-326. 10.1016/j.cell.2006.02.041.PubMedView ArticleGoogle Scholar
- Brock HW, Hodgson JW, Petruk S, Mazo A: Regulatory noncoding RNAs at Hox loci. Biochem Cell Biol. 2009, 87: 27-34. 10.1139/O08-108.PubMedView ArticleGoogle Scholar
- Lewis EB: A gene complex controlling segmentation in Drosophila. Nature. 1978, 276: 565-570. 10.1038/276565a0.PubMedView ArticleGoogle Scholar
- Ingham PW: trithorax and the regulation of homeotic gene expression in Drosophila: a historical perspective. Int J Dev Biol. 1998, 42: 423-429.PubMedGoogle Scholar
- Mallo M, Wellik DM, Deschamps J: Hox genes and regional patterning of the vertebrate body plan. Dev Biol. 2010, 344: 7-15. 10.1016/j.ydbio.2010.04.024.PubMedPubMed CentralView ArticleGoogle Scholar
- Petruk S, Sedkov Y, Smith S, Tillib S, Kraevski V, Nakamura T, Canaani E, Croce CM, Mazo A: Trithorax and dCBP acting in a complex to maintain expression of a homeotic gene. Science. 2001, 294: 1331-1334. 10.1126/science.1065683.PubMedView ArticleGoogle Scholar
- Yokoyama A, Wang Z, Wysocka J, Sanyal M, Aufiero DJ, Kitabayashi I, Herr W, Cleary ML: Leukemia proto-oncoprotein MLL forms a SET1-like histone methyltransferase complex with menin to regulate Hox gene expression. Mol Cell Biol. 2004, 24: 5639-5649. 10.1128/MCB.24.13.5639-5649.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Papoulas O, Beek SJ, Moseley SL, McCallum CM, Sarte M, Shearn A, Tamkun JW: The Drosophila trithorax group proteins BRM, ASH1 and ASH2 are subunits of distinct protein complexes. Development. 1998, 125: 3955-3966.PubMedGoogle Scholar
- Shannon M, Ashworth LK, Mucenski ML, Lamerdin JE, Branscomb E, Stubbs L: Comparative analysis of a conserved zinc finger gene cluster on human chromosome 19q and mouse chromosome 7. Genomics. 1996, 33: 112-120. 10.1006/geno.1996.0166.PubMedView ArticleGoogle Scholar
- Vogel MJ, Guelen L, de Wit E, Peric-Hupkes D, Loden M, Talhout W, Feenstra M, Abbas B, Classen AK, van Steensel B: Human heterochromatin proteins form large domains containing KRAB-ZNF genes. Genome Res. 2006, 16: 1493-1504. 10.1101/gr.5391806.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z, Yik JH, Chen R, He N, Jang MK, Ozato K, Zhou Q: Recruitment of P-TEFb for stimulation of transcriptional elongation by the bromodomain protein Brd4. Mol Cell. 2005, 19: 535-545. 10.1016/j.molcel.2005.06.029.PubMedView ArticleGoogle Scholar
- Atkinson SP, Koch CM, Clelland GK, Willcox S, Fowler JC, Stewart R, Lako M, Dunham I, Armstrong L: Epigenetic marking prepares the human HOXA cluster for activation during differentiation of pluripotent cells. Stem Cells. 2008, 26: 1174-1185. 10.1634/stemcells.2007-0497.PubMedView ArticleGoogle Scholar
- Moving AHEAD with an international human epigenome project. Nature. 2008, 454: 711-715. 10.1038/454711a.
- Kashyap V, Gudas LJ, Brenet F, Funk P, Viale A, Scandura JM: Epigenomic reorganization of the clustered Hox genes in embryonic stem cells induced by retinoic acid. J Biol Chem. 2011, 286: 3250-3260. 10.1074/jbc.M110.157545.PubMedPubMed CentralView ArticleGoogle Scholar
- Ku M, Koche RP, Rheinbay E, Mendenhall EM, Endoh M, Mikkelsen TS, Presser A, Nusbaum C, Xie X, Chi AS, Adi M, Kasif S, Ptaszek LM, Cowan CA, Lander ES, Bernstein BE: Genomewide analysis of PRC1 and PRC2 occupancy identifies two classes of bivalent domains. PLoS Genet. 2008, 4: e1000242-10.1371/journal.pgen.1000242.PubMedPubMed CentralView ArticleGoogle Scholar
- LeRoy G, Weston JT, Zee BM, Young NL, Plazas-Mayorca MD, Garcia BA: Heterochromatin protein 1 is extensively decorated with histone code-like post-translational modifications. Mol Cell Proteomics. 2009, 8: 2432-2442. 10.1074/mcp.M900160-MCP200.PubMedPubMed CentralView ArticleGoogle Scholar
- DiMaggio PA, Young NL, Baliban RC, Garcia BA, Floudas CA: A mixed integer linear optimization framework for the identification and quantification of targeted post-translational modifications of highly modified proteins using multiplexed electron transfer dissociation tandem mass spectrometry. Mol Cell Proteomics. 2009, 8: 2527-2543. 10.1074/mcp.M900144-MCP200.PubMedPubMed CentralView ArticleGoogle Scholar
- Proteome Commons using the following Tranche hash: ihSXDmjmtLL5W8cn2t5+NkkENajlU9T0D9hLPkJ3nvLiF9MzW0FIgwyjbQlJXUUMiMiSWB9zK/7hUyqKSt50ViC/lYgAAAAAAAABtg==. [https://proteomecommons.org/]
- Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, Schmidt D, O'Keeffe S, Haas S, Vingron M, Lehrach H, Yaspo ML: A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008, 321: 956-960. 10.1126/science.1160342.PubMedView ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMedPubMed CentralView ArticleGoogle Scholar
- Zang C, Schones DE, Zeng C, Cui K, Zhao K, Peng W: A clustering approach for identification of enriched domains from histone modification ChIP-Seq data. Bioinformatics. 2009, 25: 1952-1958. 10.1093/bioinformatics/btp340.PubMedPubMed CentralView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.PubMedView ArticleGoogle Scholar
- UCSC Genome Browser. [http://genome.ucsc.edu/]
- Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.