HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

Analysis of DNA Methylation Patterns in Single Blastocysts by Pyrosequencing®

Extensive epigenetic reprogramming occurs during mammalian gametogenesis and preimplantation development. DNA methylation patterns that are laid down during these stages are essential for subsequent normal foetal development. The requirement for more precise assessment of the epigenetic programming of in vitro-derived human preimplantation embryo has become of paramount importance following the identification of epigenetic diseases that are associated with assisted reproduction and/or infertility. Such techniques are also useful and applicable to experimental reproductive biology. In order to expand our knowledge of epigenetic marks, including DNA methylation, during mammalian reproduction and early development, it is necessary to test new and sufficiently sensitive protocols. There are, however, unique challenges to obtain DNA methylation data from the small cell numbers that are present in the preimplantation embryo. In this protocol, we describe the successful application of Pyrosequencing(®) to yield quantitative DNA methylation data over several CpG sites at differentially methylated regions (DMRs) at imprinted loci in single blastocysts, in this case, human blastocysts. Future developments of the protocol will allow DNA methylation analysis of a more extensive panel of genes for each embryo and at the same time, since the protocol allows for the extraction of mRNA from the embryo, the comparison between DNA methylation and gene expression.

Quantitative analysis of DNA methylation of imprinted genes in single human blastocysts by pyrosequencing

We report the first quantitative assessment of DNA methylation for any gene in the human preimplantation embryo to reveal that imprints exist at KvDMR1, RB1, SNRPN, and GRB10 in the human blastocyst. For comparison, in two human embryonic stem cell lines, imprints were also observed at KvDMR1, SNRPN, GRB10, and other imprinted loci, whereas RB1 and MEG3 were hypermethylated.

Quantitative analysis of DNA methylation at all human imprinted regions reveals preservation of epigenetic stability in adult somatic tissue

BACKGROUND

Genes subject to genomic imprinting are mono-allelically expressed in a parent-of-origin dependent manner. Each imprinted locus has at least one differentially methylated region (DMR) which has allele specific DNA methylation and contributes to imprinted gene expression. Once DMRs are established, they are potentially able to withstand normal genome reprogramming events that occur during cell differentiation and germ-line DMRs are stably maintained throughout development. These DMRs, in addition to being either maternally or paternally methylated, have differences in whether methylation was acquired in the germ-line or post fertilization and are present in a variety of genomic locations with different Cytosine-phosphate guanine (CpG) densities and CTCF binding capacities. We therefore examined the stability of maintenance of DNA methylation imprints and determined the normal baseline DNA methylation levels in several adult tissues for all imprinted genes. In order to do this, we first developed and validated 50 highly specific, quantitative DNA methylation pyrosequencing assays for the known DMRs associated with human imprinted genes.

RESULTS

Remarkable stability of the DNA methylation imprint was observed in all germ-line DMRs and paternally methylated somatic DMRs (which maintained average methylation levels of between 35% - 65% in all somatic tissues, independent of gene expression). Maternally methylated somatic DMRs were found to have more variation with tissue specific methylation patterns. Most DMRs, however, showed some intra-individual variability for DNA methylation levels in peripheral blood, suggesting that more than one DMR needs to be examined in order to get an overall impression of the epigenetic stability in a tissue. The plasticity of DNA methylation at imprinted genes was examined in a panel of normal and cancer cell lines. All cell lines showed changes in DNA methylation, especially at the paternal germ-line and the somatic DMRs.

CONCLUSIONS

Our validated pyrosequencing methylation assays can be widely used as a tool to investigate DNA methylation levels of imprinted genes in clinical samples. This first comprehensive analysis of normal methylation levels in adult somatic tissues at human imprinted regions confirm that, despite intra-individual variability and tissue specific expression, imprinted genes faithfully maintain their DNA methylation in healthy adult tissue. DNA methylation levels of a selection of imprinted genes are, therefore, a valuable indicator for epigenetic stability.

Cohesin is required for higher-order chromatin conformation at the imprinted IGF2-H19 locus

Cohesin is a chromatin-associated protein complex that mediates sister chromatid cohesion by connecting replicated DNA molecules. Cohesin also has important roles in gene regulation, but the mechanistic basis of this function is poorly understood. In mammalian genomes, cohesin co-localizes with CCCTC binding factor (CTCF), a zinc finger protein implicated in multiple gene regulatory events. At the imprinted IGF2-H19 locus, CTCF plays an important role in organizing allele-specific higher-order chromatin conformation and functions as an enhancer blocking transcriptional insulator. Here we have used chromosome conformation capture (3C) assays and RNAi-mediated depletion of cohesin to address whether cohesin affects higher order chromatin conformation at the IGF2-H19 locus in human cells. Our data show that cohesin has a critical role in maintaining CTCF-mediated chromatin conformation at the locus and that disruption of this conformation coincides with changes in IGF2 expression. We show that the cohesin-dependent, higher-order chromatin conformation of the locus exists in both G1 and G2 phases of the cell cycle and is therefore independent of cohesin's function in sister chromatid cohesion. We propose that cohesin can mediate interactions between DNA molecules in cis to insulate genes through the formation of chromatin loops, analogous to the cohesin mediated interaction with sister chromatids in trans to establish cohesion.

Generation of an epigenetic signature by chronic hypoxia in prostate cells

Increasing levels of tissue hypoxia have been reported as a natural feature of the aging prostate gland and may be a risk factor for the development of prostate cancer. In this study, we have used PwR-1E benign prostate epithelial cells and an equivalently aged hypoxia-adapted PwR-1E sub-line to identify phenotypic and epigenetic consequences of chronic hypoxia in prostate cells. We have identified a significantly altered cellular phenotype in response to chronic hypoxia as characterized by increased receptor-mediated apoptotic resistance, the induction of cellular senescence, increased invasion and the increased secretion of IL-1 beta, IL6, IL8 and TNFalpha cytokines. In association with these phenotypic changes and the absence of HIF-1 alpha protein expression, we have demonstrated significant increases in global levels of DNA methylation and H3K9 histone acetylation in these cells, concomitant with the increased expression of DNA methyltransferase DMNT3b and gene-specific changes in DNA methylation at key imprinting loci. In conclusion, we have demonstrated a genome-wide adjustment of DNA methylation and histone acetylation under chronic hypoxic conditions in the prostate. These epigenetic signatures may represent an additional mechanism to promote and maintain a hypoxic-adapted cellular phenotype with a potential role in tumour development.

Transcript- and tissue-specific imprinting of a tumour suppressor gene

The Bladder Cancer-Associated Protein gene (BLCAP; previously BC10) is a tumour suppressor that limits cell proliferation and stimulates apoptosis. BLCAP protein or message are downregulated or absent in a variety of human cancers. In mouse and human, the first intron of Blcap/BLCAP contains the distinct Neuronatin (Nnat/NNAT) gene. Nnat is an imprinted gene that is exclusively expressed from the paternally inherited allele. Previous studies found no evidence for imprinting of Blcap in mouse or human. Here we show that Blcap is imprinted in mouse and human brain, but not in other mouse tissues. Moreover, Blcap produces multiple distinct transcripts that exhibit reciprocal allele-specific expression in both mouse and human. We propose that the tissue-specific imprinting of Blcap is due to the particularly high transcriptional activity of Nnat in brain, as has been suggested previously for the similarly organized and imprinted murine Commd1/U2af1-rs1 locus. For Commd1/U2af1-rs1, we show that it too produces distinct transcript variants with reciprocal allele-specific expression. The imprinted expression of BLCAP and its interplay with NNAT at the transcriptional level may be relevant to human carcinogenesis.

Regulation of alternative polyadenylation by genomic imprinting

Maternally and paternally derived alleles can utilize different promoters, but allele-specific differences in cotranscriptional processes have not been reported. We show that alternative polyadenylation sites at a novel murine imprinted gene (H13) are utilized in an allele-specific manner. A differentially methylated CpG island separates polyA sites utilized on maternal and paternal alleles, and contains an internal promoter. Two genetic systems show that alleles lacking methylation generate truncated H13 transcripts that undergo internal polyadenylation. On methylated alleles, the internal promoter is inactive and elongation proceeds to downstream polyadenylation sites. This demonstrates that epigenetic modifications can influence utilization of alternative polyadenylation sites.

Somatically acquired hypomethylation of IGF2 in breast and colorectal cancer

The imprinted insulin-like growth factor 2 (IGF2) gene is expressed predominantly from the paternal allele. Loss of imprinting (LOI) associated with hypomethylation at the promoter proximal sequence (DMR0) of the IGF2 gene was proposed as a predisposing constitutive risk biomarker for colorectal cancer. We used pyrosequencing to assess whether IGF2 DMR0 methylation is either present constitutively prior to cancer or whether it is acquired tissue-specifically after the onset of cancer. DNA samples from tumour tissues and matched non-tumour tissues from 22 breast and 42 colorectal cancer patients as well as peripheral blood samples obtained from colorectal cancer patients [SEARCH (n=case 192, controls 96)], breast cancer patients [ABC (n=case 364, controls 96)] and the European Prospective Investigation of Cancer [EPIC-Norfolk (n=breast 228, colorectal 225, controls 895)] were analysed. The EPIC samples were collected 2–5 years prior to diagnosis of breast or colorectal cancer. IGF2 DMR0 methylation levels in tumours were lower than matched non-tumour tissue. Hypomethylation of DMR0 was detected in breast (33%) and colorectal (80%) tumour tissues with a higher frequency than LOI indicating that methylation levels are a better indicator of cancer than LOI. In the EPIC population, the prevalence of IGF2 DMR0 hypomethylation was 9.5% and this correlated with increased age not cancer risk. Thus, IGF2 DMR0 hypomethylation occurs as an acquired tissue-specific somatic event rather than a constitutive innate epimutation. These results indicate that IGF2 DMR0 hypomethylation has diagnostic potential for colon cancer rather than value as a surrogate biomarker for constitutive LOI.

Distinct methylation changes at the IGF2-H19 locus in congenital growth disorders and cancer

Background

Differentially methylated regions (DMRs) are associated with many imprinted genes. In mice methylation at a DMR upstream of the H19 gene known as the Imprint Control region (IC1) is acquired in the male germline and influences the methylation status of DMRs 100 kb away in the adjacent Insulin-like growth factor 2 (Igf2) gene through long-range interactions. In humans, germline-derived or post-zygotically acquired imprinting defects at IC1 are associated with aberrant activation or repression of IGF2, resulting in the congenital growth disorders Beckwith-Wiedemann (BWS) and Silver-Russell (SRS) syndromes, respectively. In Wilms tumour and colorectal cancer, biallelic expression of IGF2 has been observed in association with loss of methylation at a DMR in IGF2. This DMR, known as DMR0, has been shown to be methylated on the silent maternal IGF2 allele presumably with a role in repression. The effect of IGF2 DMR0 methylation changes in the aetiology of BWS or SRS is unknown.

Methodology/Principal Findings

We analysed the methylation status of the DMR0 in BWS, SRS and Wilms tumour patients by conventional bisulphite sequencing and pyrosequencing. We show here that, contrary to previous reports, the IGF2 DMR0 is actually methylated on the active paternal allele in peripheral blood and kidney. This is similar to the IC1 methylation status and is inconsistent with the proposed silencing function of the maternal IGF2 allele. Beckwith-Wiedemann and Silver-Russell patients with IC1 methylation defects have similar methylation defects at the IGF2 DMR0, consistent with IC1 regulating methylation at IGF2 in cis. In Wilms tumour, however, methylation profiles of IC1 and IGF2 DMR0 are indicative of methylation changes occurring on both parental alleles rather than in cis.

Conclusions/Significance

These results support a model in which DMR0 and IC1 have opposite susceptibilities to global hyper and hypomethylation during tumorigenesis independent of the parent of origin imprint. In contrast, during embryogenesis DMR0 is methylated or demethylated according to the germline methylation imprint at the IC1, indicating different mechanisms of imprinting loss in neoplastic and non-neoplastic cells.

WAMIDEX: a web atlas of murine genomic imprinting and differential expression

The mouse is an established model organism for the study of genomic imprinting. Mice with genetic material originating from only one parent (e.g., mice with uniparental chromosomal duplications) or gene mutations leading to epigenetic deficiencies have proven to be particularly useful tools. In the process of our studies we have accumulated a large set of expression microarray measurements in samples derived from these types of mice. Here, we present the collation of these and third-party microarray data that are relevant to genomic imprinting into a Web Atlas of Murine genomic Imprinting and Differential EXpression (WAMIDEX: https://atlas.genetics.kcl.ac.uk). WAMIDEX integrates the most comprehensive literature-derived catalog of murine imprinted genes to date with a genome browser that makes the microarray data immediately accessible in annotation-rich genomic context. In addition, WAMIDEX exemplifies the use of the self-organizing map method for the discovery of novel imprinted genes from microarray data. The parent-of-origin-specific expression of imprinted genes is frequently limited to specific tissues or developmental stages, a fact that the atlas reflects in its design and data content.

Chromosome-wide identification of novel imprinted genes using microarrays and uniparental disomies

Genomic imprinting refers to a specialized form of epigenetic gene regulation whereby the expression of a given allele is dictated by parental origin. Defining the extent and distribution of imprinting across genomes will be crucial for understanding the roles played by imprinting in normal mammalian growth and development. Using mice carrying uniparental disomies or duplications, microarray screening and stringent bioinformatics, we have developed the first large-scale tissue-specific screen for imprinted gene detection. We quantify the stringency of our methodology and relate it to previous non-tissue-specific large-scale studies. We report the identification in mouse of four brain-specific novel paternally expressed transcripts and an additional three genes that show maternal expression in the placenta. The regions of conserved linkage in the human genome are associated with the Prader-Willi Syndrome (PWS) and Beckwith-Wiedemann Syndrome (BWS) where imprinting is known to be a contributing factor. We conclude that large-scale systematic analyses of this genre are necessary for the full impact of genomic imprinting on mammalian gene expression and phenotype to be elucidated.

The DNA sequence and biological annotation of human chromosome 1

The reference sequence for each human chromosome provides the framework for understanding genome function, variation and evolution. Here we report the finished sequence and biological annotation of human chromosome 1. Chromosome 1 is gene-dense, with 3,141 genes and 991 pseudogenes, and many coding sequences overlap. Rearrangements and mutations of chromosome 1 are prevalent in cancer and many other diseases. Patterns of sequence variation reveal signals of recent selection in specific genes that may contribute to human fitness, and also in regions where no function is evident. Fine-scale recombination occurs in hotspots of varying intensity along the sequence, and is enriched near genes. These and other studies of human biology and disease encoded within chromosome 1 are made possible with the highly accurate annotated sequence, as part of the completed set of chromosome sequences that comprise the reference human genome.

Investigating chromosome organization with genomic microarrays

DNA microarrays are increasingly being used to investigate the functional role of chromatin. These studies are enhanced by the development of high-resolution arrays covering either the whole genome or specific regions of selected chromosomes with large insert clones, PCR products or oligonucleotides of around 100 bp or less. In combination with chromatin immunoprecipitation, this approach allows identification of protein binding for transcription factors, proteins involved in DNA replication and repair as well as sites of chromatin modification. Furthermore, by application of S phase fractions to genomic microarrays, replication timing can be estimated. Thus, microarrays can provide new information about chromosome structure and gene regulation.

Replication timing of human chromosome 6

Genomic microarrays have been used to assess DNA replication timing in a variety of eukaryotic organisms. A replication timing map of the human genome has already been published at a 1Mb resolution. Here we describe how the same method can be used to assess the replication timing of chromosome 6 with a greater resolution using an array of overlapping tile path clones. We report the replication timing map of the whole of chromosome 6 in general, and the MHC region in particular. Positive correlations are observed between replication timing and a number of genomic features including GC content, repeat content and transcriptional activity.

A novel 5q11.2 deletion detected by microarray comparative genomic hybridisation in a child referred as a case of suspected 22q11 deletion syndrome

The 22q11 deletion syndrome (22q11DS) is a developmental syndrome comprising of heart, palate, thymus and parathyroid glands defects. Individuals with 22q11DS usually carry a 1.5- to 3-Mb heterozygous deletion on chromosome 22q11.2. However, there are many patients with features of 22q11DS without a known cause from conventional karyotype and FISH analysis. Six patients with features of 22q11DS, a normal chromosomal and FISH 22q11 analysis, were selected for investigation by microarray genomic comparative hybridisation (array CGH). Array-CGH is a powerful technology enabling detection of submicroscopic chromosome duplications and deletions by comparing a differentially labelled test sample to a control. The samples are co-hybridised to a microarray containing genomic clones and the resulting ratio of fluorescence intensities on each array element is proportional to the DNA copy number difference. No chromosomal changes were detected by hybridisation to a high resolution array representing chromosome 22q. However, one patient was found to have a 6-Mb deletion on 5q11.2 detected by a whole genome 1-Mb array. This deletion was confirmed with fluorescence in-situ hybridisation (FISH) and microsatellite marker analysis. It is the first deletion described in this region. The patient had tetralogy of Fallot, a bifid uvula and velopharyngeal insufficiency, short stature, learning and behavioural difficulties. This case shows the increased sensitivity of array CGH over detailed karyotype analysis for detection of chromosomal changes. It is anticipated that array CGH will improve the clinician's capacity to diagnose congenital syndromes with an unknown aetiology.

Chromatin architecture of the human genome: gene-rich domains are enriched in open chromatin fibers

We present an analysis of chromatin fiber structure across the human genome. Compact and open chromatin fiber structures were separated by sucrose sedimentation and their distributions analyzed by hybridization to metaphase chromosomes and genomic microarrays. We show that compact chromatin fibers originate from some sites of heterochromatin (C-bands), and G-bands (euchromatin). Open chromatin fibers correlate with regions of highest gene density, but not with gene expression since inactive genes can be in domains of open chromatin, and active genes in regions of low gene density can be embedded in compact chromatin fibers. Moreover, we show that chromatin fiber structure impacts on further levels of chromatin condensation. Regions of open chromatin fibers are cytologically decondensed and have a distinctive nuclear organization. We suggest that domains of open chromatin may create an environment that facilitates transcriptional activation and could provide an evolutionary constraint to maintain clusters of genes together along chromosomes.

Replication timing of the human genome

We have developed a directly quantitative method utilizing genomic clone DNA microarrays to assess the replication timing of sequences during the S phase of the cell cycle. The genomic resolution of the replication timing measurements is limited only by the genomic clone size and density. We demonstrate the power of this approach by constructing a genome-wide map of replication timing in human lymphoblastoid cells using an array with clones spaced at 1 Mb intervals and a high-resolution replication timing map of 22q with an array utilizing overlapping sequencing tile path clones. We show a positive correlation, both genome-wide and at a high resolution, between replication timing and a range of genome parameters including GC content, gene density and transcriptional activity.