ChIP on SNP-chip for genome-wide analysis of human histone H4 hyperacetylation

Allele-specific transcriptional elongation regulates monoallelic expression of the IGF2BP1 gene

Background

Random monoallelic expression contributes to phenotypic variation of cells and organisms. However, the epigenetic mechanisms by which individual alleles are randomly selected for expression are not known. Taking cues from chromatin signatures at imprinted gene loci such as the insulin-like growth factor 2 gene 2 (IGF2), we evaluated the contribution of CTCF, a zinc finger protein required for parent-of-origin-specific expression of the IGF2 gene, as well as a role for allele-specific association with DNA methylation, histone modification and RNA polymerase II.

Results

Using array-based chromatin immunoprecipitation, we identified 293 genomic loci that are associated with both CTCF and histone H3 trimethylated at lysine 9 (H3K9me3). A comparison of their genomic positions with those of previously published monoallelically expressed genes revealed no significant overlap between allele-specifically expressed genes and colocalized CTCF/H3K9me3. To analyze the contributions of CTCF and H3K9me3 to gene regulation in more detail, we focused on the monoallelically expressed IGF2BP1 gene. In vitro binding assays using the CTCF target motif at the IGF2BP1 gene, as well as allele-specific analysis of cytosine methylation and CTCF binding, revealed that CTCF does not regulate mono- or biallelic IGF2BP1 expression. Surprisingly, we found that RNA polymerase II is detected on both the maternal and paternal alleles in B lymphoblasts that express IGF2BP1 primarily from one allele. Thus, allele-specific control of RNA polymerase II elongation regulates the allelic bias of IGF2BP1 gene expression.

Conclusions

Colocalization of CTCF and H3K9me3 does not represent a reliable chromatin signature indicative of monoallelic expression. Moreover, association of individual alleles with both active (H3K4me3) and silent (H3K27me3) chromatin modifications (allelic bivalent chromatin) or with RNA polymerase II also fails to identify monoallelically expressed gene loci. The selection of individual alleles for expression occurs in part during transcription elongation.

Allele specific analysis of the ADRBK2 gene in lymphoblastoid cells from bipolar disorder patients

G-protein coupled receptor kinase-3 (GRK3), translated from the gene, ADRBK2 has been implicated as a candidate molecule for bipolar disorder through multiple, converging lines of evidence. In some individuals, the ADRBK2 gene harbors the A-haplotype, a collection of single nucleotide polymorphisms (SNPs) previously associated with an increased risk for bipolar disorder. Because the A-haplotype encompasses the ADRBK2 promoter, we hypothesized that it may alter the regulation of gene expression. Using histone H3 acetylation to infer promoter activity in lymphoblastoid cells from patients with bipolar disorder, we examined the A-haplotype within its genomic context and determined that at least four of its SNPs are present in transcriptionally active portions of the promoter. However, using chromatin immunoprecipitation followed by allele-specific PCR in samples heterozygous for the A-haplotype, we found no evidence of altered levels of acetylated histone H3 at the affected allele compared to the common allele. Similarly, using a transcribed SNP to discriminate expressed ADRBK2 mRNA strands by allele of origin; we found that the A-haplotype did not confer an allelic-expression imbalance. Our data suggest that while the A-haplotype is situated in active regulatory sequence, the risk-associated SNPs do not appear to affect ADRBK2 gene regulation at the level of histone H3 acetylation nor do they confer measurable changes in transcription in lymphoblastoid cells. However, tissue-specific mechanisms by which the A-haplotype could affect ADRBK2 in the central nervous system cannot be excluded.

DNaseI hypersensitivity at gene-poor, FSH dystrophy-linked 4q35.2

A subtelomeric region, 4q35.2, is implicated in facioscapulohumeral muscular dystrophy (FSHD), a dominant disease thought to involve local pathogenic changes in chromatin. FSHD patients have too few copies of a tandem 3.3-kb repeat (D4Z4) at 4q35.2. No phenotype is associated with having few copies of an almost identical repeat at 10q26.3. Standard expression analyses have not given definitive answers as to the genes involved. To investigate the pathogenic effects of short D4Z4 arrays on gene expression in the very gene-poor 4q35.2 and to find chromatin landmarks there for transcription control, unannotated genes and chromatin structure, we mapped DNaseI-hypersensitive (DH) sites in FSHD and control myoblasts. Using custom tiling arrays (DNase-chip), we found unexpectedly many DH sites in the two large gene deserts in this 4-Mb region. One site was seen preferentially in FSHD myoblasts. Several others were mapped >0.7 Mb from genes known to be active in the muscle lineage and were also observed in cultured fibroblasts, but not in lymphoid, myeloid or hepatic cells. Their selective occurrence in cells derived from mesoderm suggests functionality. Our findings indicate that the gene desert regions of 4q35.2 may have functional significance, possibly also to FSHD, despite their paucity of known genes.

Use of SNP-arrays for ChIP assays: computational aspects

The simultaneous genotyping of thousands of single nucleotide polymorphisms (SNPs) in a genome using SNP-Arrays is a very important tool that is revolutionizing genetics and molecular biology. We expanded the utility of this technique by using it following chromatin immunoprecipitation (ChIP) to assess the multiple genomic locations protected by a protein complex recognized by an antibody. The power of this technique is illustrated through an analysis of the changes in histone H4 acetylation, a marker of open chromatin and transcriptionally active genomic regions, which occur during differentiation of human myoblasts into myotubes. The findings have been validated by the observation of a significant correlation between the detected histone modifications and the expression of the nearby genes, as measured by DNA expression microarrays. This chapter focuses on the computational analysis of the data.

Overlapping euchromatin/heterochromatin- associated marks are enriched in imprinted gene regions and predict allele-specific modification

Most genome-level analysis treats the two parental alleles equivalently, yet diploid genomes contain two parental genomes that are often epigenetically distinct. While single nucleotide polymorphisms (SNPs) can be used to distinguish these genomes, it would be useful to develop a generalized strategy for identifying candidate genes or regions showing allele-specific differences, independent of SNPs. We have explored this problem by looking for overlapping marks in the genome related to both euchromatin (histone H3 dimethyl lysine-4 [H3K4Me2]) and heterochromatin (DNA methylation [DNAm]). "Double hits" were defined by the intersection of H3K4Me2 and DNAm. For the top 5% of marks, defined by a sliding window, imprinted gene regions were enriched for double hits 5.4-fold. When the location information of CTCF binding sites were integrated, the "triple hits" were enriched 76-fold for known imprinted genes in the regions studied. The double hits in imprinted genes were found to occur usually at the site of alternative or antisense transcripts. In addition, four of four imprinted genes tested showing double hits also showed allele-specific methylation. We suggest that overlapping euchromatin/heterochromatin marks are common and are enriched for epigenetically distinct parental chromosome regions. Furthermore, we developed a novel approach to identifying allele-specific marks that is SNP independent, by fractionating using H3K4Me2 antibodies followed by DNA methylation analysis.

Specific tumour-associated methylation in normal human term placenta and first-trimester cytotrophoblasts

Human placentation displays many similarities with tumourigenesis, including rapid cell division, migration and invasion, overlapping gene expression profiles and escape from immune detection. Recent data have identified promoter methylation in the Ras association factor and adenomatous polyposis coli tumour suppressor genes as part of this process. However, the extent of tumour-associated methylation in the placenta remains unclear. Using whole genome methylation data as a starting point, we have examined this phenomenon in placental tissue. We found no evidence for methylation of the majority of common tumour suppressor genes in term placentas, but identified methylation in several genes previously described in some human tumours. Notably, promoter methylation of four independent negative regulators of Wnt signalling has now been identified in human placental tissue and purified trophoblasts. Methylation is present in baboon, but not in mouse placentas. This supports a role for elevated Wnt signalling in primate trophoblast invasiveness and placentation. Examination of invasive choriocarcinoma cell lines revealed altered methylation patterns consistent with a role of methylation change in gestational trophoblastic disease. This distinct pattern of tumour-associated methylation implicates a coordinated series of epigenetic silencing events, similar to those associated with some tumours, in the distinct features of normal human placental invasion and function.

Genome-wide mapping of allele-specific protein-DNA interactions in human cells

We describe a high-throughput method, named ChIP-SNP, for the identification of allele-specific protein-DNA interactions throughout the human genome. ChIP-SNP combines chromatin immunoprecipitation (ChIP) with whole-genome single nucleotide polymorphism (SNP) genotyping microarray analysis. We demonstrated that it can be used to accurately identify allele-specific binding of RNA polymerase II (RNAP) in the human fibroblast genome, uncovering imprinted genes and other allele-specific binding events. ChIP-SNP will facilitate the study of mechanisms of allele-specific gene expression.