Shalima S. Nair

Comparison of methyl-DNA immunoprecipitation (MeDIP) and methyl-CpG binding domain (MBD) protein capture for genome-wide DNA methylation analysis reveal CpG sequence coverage bias

DNA methylation primarily occurs at CpG dinucleotides in mammals and is a common epigenetic mark that plays a critical role in the regulation of gene expression. Profiling DNA methylation patterns across the genome is vital to understand DNA methylation changes that occur during development and in disease phenotype. In this study, we compared two commonly used approaches to enrich for methylated DNA regions of the genome, namely methyl-DNA immunoprecipitation (MeDIP) that is based on enrichment with antibodies specific for 5′-methylcytosine (5MeC), and capture of methylated DNA using a methyl-CpG binding domain-based (MBD) protein to discover differentially methylated regions (DMRs) in cancer. The enriched methylated DNA fractions were interrogated on Affymetrix promoter tiling arrays and differentially methylated regions were identified. A detailed validation study of 42 regions was performed using Sequenom MassCLEAVE technique. This detailed analysis revealed that both enrichment techniques are sensitive for detecting DMRs and preferentially identified different CpG rich regions of the prostate cancer genome, with MeDIP commonly enriching for methylated regions with a low CpG density, while MBD capture favors regions of higher CpG density and identifies the greatest proportion of CpG islands. This is the first detailed validation report comparing different methylated DNA enrichment techniques for identifying regions of differential DNA methylation. Our study highlights the importance of understanding the nuances of the methods used for DNA genome-wide methylation analyses so that accurate interpretation of the biology is not overlooked.

Evaluation of affinity-based genome-wide DNA methylation data: Effects of CpG density, amplification bias, and copy number variation

DNA methylation is an essential epigenetic modification that plays a key role associated with the regulation of gene expression during differentiation, but in disease states such as cancer, the DNA methylation landscape is often deregulated. There are now numerous technologies available to interrogate the DNA methylation status of CpG sites in a targeted or genome-wide fashion, but each method, due to intrinsic biases, potentially interrogates different fractions of the genome. In this study, we compare the affinity-purification of methylated DNA between two popular genome-wide techniques, methylated DNA immunoprecipitation (MeDIP) and methyl-CpG binding domain-based capture (MBDCap), and show that each technique operates in a different domain of the CpG density landscape. We explored the effect of whole-genome amplification and illustrate that it can reduce sensitivity for detecting DNA methylation in GC-rich regions of the genome. By using MBDCap, we compare and contrast microarray- and sequencing-based readouts and highlight the impact that copy number variation (CNV) can make in differential comparisons of methylomes. These studies reveal that the analysis of DNA methylation data and genome coverage is highly dependent on the method employed, and consideration must be made in light of the GC content, the extent of DNA amplification, and the copy number.

Acetylation of H2A.Z is a key epigenetic modification associated with gene deregulation and epigenetic remodeling in cancer

Histone H2A.Z (H2A.Z) is an evolutionarily conserved H2A variant implicated in the regulation of gene expression; however, its role in transcriptional deregulation in cancer remains poorly understood. Using genome-wide studies, we investigated the role of promoter-associated H2A.Z and acetylated H2A.Z (acH2A.Z) in gene deregulation and its relationship with DNA methylation and H3K27me3 in prostate cancer. Our results reconcile the conflicting reports of positive and negative roles for histone H2A.Z and gene expression states. We find that H2A.Z is enriched in a bimodal distribution at nucleosomes, surrounding the transcription start sites (TSSs) of both active and poised gene promoters. In addition, H2A.Z spreads across the entire promoter of inactive genes in a deacetylated state. In contrast, acH2A.Z is only localized at the TSSs of active genes. Gene deregulation in cancer is also associated with a reorganization of acH2A.Z and H2A.Z nucleosome occupancy across the promoter region and TSS of genes. Notably, in cancer cells we find that a gain of acH2A.Z at the TSS occurs with an overall decrease of H2A.Z levels, in concert with oncogene activation. Furthermore, deacetylation of H2A.Z at TSSs is increased with silencing of tumor suppressor genes. We also demonstrate that acH2A.Z anti-correlates with promoter H3K27me3 and DNA methylation. We show for the first time, that acetylation of H2A.Z is a key modification associated with gene activity in normal cells and epigenetic gene deregulation in tumorigenesis.

MicroRNA Alterations and Associated Aberrant DNA Methylation Patterns across Multiple Sample Types in Oral Squamous Cell Carcinoma

Background MicroRNA (miRNA) expression is broadly altered in cancer, but few studies have investigated miRNA deregulation in oral squamous cell carcinoma (OSCC). Epigenetic mechanisms are involved in the regulation of >30 miRNA genes in a range of tissues, and we aimed to investigate this further in OSCC. Methods TaqMan® qRT-PCR arrays and individual assays were used to profile miRNA expression in a panel of 25 tumors with matched adjacent tissues from patients with OSCC, and 8 control paired oral stroma and epithelium from healthy volunteers. Associated DNA methylation changes of candidate epigenetically deregulated miRNA genes were measured in the same samples using the MassArray® mass spectrometry platform. MiRNA expression and DNA methylation changes were also investigated in FACS sorted CD44high oral cancer stem cells from primary tumor samples (CSCs), and in oral rinse and saliva from 15 OSCC patients and 7 healthy volunteers. Results MiRNA expression patterns were consistent in healthy oral epithelium and stroma, but broadly altered in both tumor and adjacent tissue from OSCC patients. MiR-375 is repressed and miR-127 activated in OSCC, and we confirm previous reports of miR-137 hypermethylation in oral cancer. The miR-200 s/miR-205 were epigenetically activated in tumors vs normal tissues, but repressed in the absence of DNA hypermethylation specifically in CD44high oral CSCs. Aberrant miR-375 and miR-200a expression and miR-200c-141 methylation could be detected in and distinguish OSCC patient oral rinse and saliva from healthy volunteers, suggesting a potential clinical application for OSCC specific miRNA signatures in oral fluids. Conclusions MiRNA expression and DNA methylation changes are a common event in OSCC, and we suggest miR-375, miR-127, miR-137, the miR-200 family and miR-205 as promising candidates for future investigations. Although overall activated in OSCC, miR-200/miR-205 suppression in oral CSCs indicate that cell specific silencing of these miRNAs may drive tumor expansion and progression.

Epigenetic-induced repression of microRNA-205 is associated with MED1 activation and a poorer prognosis in localized prostate cancer

Deregulation of microRNA (miRNA) expression can have a critical role in carcinogenesis. Here we show in prostate cancer that miRNA-205 (miR-205) transcription is commonly repressed and the MIR-205 locus is hypermethylated. LOC642587, the MIR-205 host gene of unknown function, is also concordantly inactivated. We show that miR-205 targets mediator 1 (MED1, also called TRAP220 and PPARBP) for transcriptional silencing in normal prostate cells, leading to reduction in MED1 mRNA levels, and in total and active phospho-MED1 protein. Overexpression of miR-205 in prostate cancer cells negatively affects cell viability, consistent with a tumor suppressor function. We found that hypermethylation of the MIR-205 locus was strongly related with a decrease in miR-205 expression and an increase in MED1 expression in primary tumor samples (n=14), when compared with matched normal prostate (n=7). An expanded patient cohort (tumor n=149, matched normal n=30) also showed significant MIR-205 DNA methylation in tumors compared with normal, and MIR-205 hypermethylation is significantly associated with biochemical recurrence (hazard ratio=2.005, 95% confidence interval (1.109, 3.625), P=0.02), in patients with low preoperative prostate specific antigen. In summary, these results suggest that miR-205 is an epigenetically regulated tumor suppressor that targets MED1 and may provide a potential biomarker in prostate cancer management.

Methylome sequencing in triple-negative breast cancer reveals distinct methylation clusters with prognostic value

Epigenetic alterations in the cancer methylome are common in breast cancer and provide novel options for tumour stratification. Here, we perform whole-genome methylation capture sequencing on small amounts of DNA isolated from formalin-fixed, paraffin-embedded tissue from triple-negative breast cancer (TNBC) and matched normal samples. We identify differentially methylated regions (DMRs) enriched with promoters associated with transcription factor binding sites and DNA hypersensitive sites. Importantly, we stratify TNBCs into three distinct methylation clusters associated with better or worse prognosis and identify 17 DMRs that show a strong association with overall survival, including DMRs located in the Wilms tumour 1 (WT1) gene, bi-directional-promoter and antisense WT1-AS. Our data reveal that coordinated hypermethylation can occur in oestrogen receptor-negative disease, and that characterizing the epigenetic framework provides a potential signature to stratify TNBCs. Together, our findings demonstrate the feasibility of profiling the cancer methylome with limited archival tissue to identify regulatory regions associated with cancer.

Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer

MicroRNAs (miRNAs) are small non-coding RNAs of ∼20 nt in length that are capable of modulating gene expression post-transcriptionally. Although miRNAs have been implicated in cancer, including breast cancer, the regulation of miRNA transcription and the role of defects in this process in cancer is not well understood. In this study we have mapped the promoters of 93 breast cancer-associated miRNAs, and then looked for associations between DNA methylation of 15 of these promoters and miRNA expression in breast cancer cells. The miRNA promoters with clearest association between DNA methylation and expression included a previously described and a novel promoter of the Hsa-mir-200b cluster. The novel promoter of the Hsa-mir-200b cluster, denoted P2, is located ∼2 kb upstream of the 5′ stemloop and maps within a CpG island. P2 has comparable promoter activity to the previously reported promoter (P1), and is able to drive the expression of miR-200b in its endogenous genomic context. DNA methylation of both P1 and P2 was inversely associated with miR-200b expression in eight out of nine breast cancer cell lines, and in vitro methylation of both promoters repressed their activity in reporter assays. In clinical samples, P1 and P2 were differentially methylated with methylation inversely associated with miR-200b expression. P1 was hypermethylated in metastatic lymph nodes compared with matched primary breast tumours whereas P2 hypermethylation was associated with loss of either oestrogen receptor or progesterone receptor. Hypomethylation of P2 was associated with gain of HER2 and androgen receptor expression. These data suggest an association between miR-200b regulation and breast cancer subtype and a potential use of DNA methylation of miRNA promoters as a component of a suite of breast cancer biomarkers.

Methylation-capture and Next-Generation Sequencing of free circulating DNA from human plasma

BackgroundFree circulating DNA (fcDNA) has many potential clinical applications, due to the non-invasive way in which it is collected. However, because of the low concentration of fcDNA in blood, genome-wide analysis carries many technical challenges that must be overcome before fcDNA studies can reach their full potential. There are currently no definitive standards for fcDNA collection, processing and whole-genome sequencing. We report novel detailed methodology for the capture of high-quality methylated fcDNA, library preparation and downstream genome-wide Next-Generation Sequencing. We also describe the effects of sample storage, processing and scaling on fcDNA recovery and quality.ResultsUse of serum versus plasma, and storage of blood prior to separation resulted in genomic DNA contamination, likely due to leukocyte lysis. Methylated fcDNA fragments were isolated from 5 donors using a methyl-binding protein-based protocol and appear as a discrete band of ~180 bases. This discrete band allows minimal sample loss at the size restriction step in library preparation for Next-Generation Sequencing, allowing for high-quality sequencing from minimal amounts of fcDNA. Following sequencing, we obtained 37×106-86×106 unique mappable reads, representing more than 50% of total mappable reads. The methylation status of 9 genomic regions as determined by DNA capture and sequencing was independently validated by clonal bisulphite sequencing.ConclusionsOur optimized methods provide high-quality methylated fcDNA suitable for whole-genome sequencing, and allow good library complexity and accurate sequencing, despite using less than half of the recommended minimum input DNA.

Comprehensive evaluation of genome-wide 5-hydroxymethylcytosine profiling approaches in human DNA

BackgroundThe discovery that 5-methylcytosine (5mC) can be oxidized to 5-hydroxymethylcytosine (5hmC) by the ten-eleven translocation (TET) proteins has prompted wide interest in the potential role of 5hmC in reshaping the mammalian DNA methylation landscape. The gold-standard bisulphite conversion technologies to study DNA methylation do not distinguish between 5mC and 5hmC. However, new approaches to mapping 5hmC genome-wide have advanced rapidly, although it is unclear how the different methods compare in accurately calling 5hmC. In this study, we provide a comparative analysis on brain DNA using three 5hmC genome-wide approaches, namely whole-genome bisulphite/oxidative bisulphite sequencing (WG Bis/OxBis-seq), Infinium HumanMethylation450 BeadChip arrays coupled with oxidative bisulphite (HM450K Bis/OxBis) and antibody-based immunoprecipitation and sequencing of hydroxymethylated DNA (hMeDIP-seq). We also perform loci-specific TET-assisted bisulphite sequencing (TAB-seq) for validation of candidate regions.ResultsWe show that whole-genome single-base resolution approaches are advantaged in providing precise 5hmC values but require high sequencing depth to accurately measure 5hmC, as this modification is commonly in low abundance in mammalian cells. HM450K arrays coupled with oxidative bisulphite provide a cost-effective representation of 5hmC distribution, at CpG sites with 5hmC levels >~10%. However, 5hmC analysis is restricted to the genomic location of the probes, which is an important consideration as 5hmC modification is commonly enriched at enhancer elements. Finally, we show that the widely used hMeDIP-seq method provides an efficient genome-wide profile of 5hmC and shows high correlation with WG Bis/OxBis-seq 5hmC distribution in brain DNA. However, in cell line DNA with low levels of 5hmC, hMeDIP-seq-enriched regions are not detected by WG Bis/OxBis or HM450K, either suggesting misinterpretation of 5hmC calls by hMeDIP or lack of sensitivity of the latter methods.ConclusionsWe highlight both the advantages and caveats of three commonly used genome-wide 5hmC profiling technologies and show that interpretation of 5hmC data can be significantly influenced by the sensitivity of methods used, especially as the levels of 5hmC are low and vary in different cell types and different genomic locations.

Multiplex bisulfite PCR resequencing of clinical FFPE DNA.

BackgroundThe clinical utility of DNA methylation as a predictive or prognostic biomarker requires scalable resequencing protocols for bisulfite-converted DNA. Key features of any validation method should be adaptability for low- or high-throughput needs and high reproducibility, and should only require minimal amounts of precious clinical sample as input material. Critically, this method should also deliver robust results when working with bisulfite-converted DNA extracted from formalin-fixed, paraffin-embedded (FFPE) blocks.ResultsWe report here for the first time on comparison studies between the Fluidigm Access Array system and multiplex assays for multiplex bisulfite PCR resequencing. The requirement of the Fluidigm Access Array system for high template amounts and its sensitivity to variations in template quality rendered it unsuitable for bisulfite PCR applications utilizing FFPE DNA. In response to this limitation, we established a multiplex bisulfite PCR assay capable of delivering robust methylation data using minimal amounts of FFPE clinical DNA. To evaluate the parameters and reproducibility of this assay, 57 amplicons were used to prepare sequencing libraries in triplicate for 13 FFPE tumour samples, as well as a series of 5 methylated controls (0%, 25%, 50%, 75%, and 100%). Analysis of this data demonstrated that this multiplex assay had high reproducibility (mean standard deviation of 1.4% for methylation values), was low cost, required low sample input (50 ng of DNA or less), and could be scaled for both low- and high-throughput needs. Notably, ExoSAP-IT (exonuclease I) treatment to remove residual primers in bisulfite resequencing libraries appeared to degrade the library and generate a high-molecular weight smear which may impact on the degree of methylation assessed.ConclusionsMultiplex bisulfite PCR assays represent a convenient and scalable method for validation and screening of methylated DNA regions from archival FFPE DNA. Moreover, the library construction process detailed here can be rapidly optimized and implemented with a minimal amount of work, can be performed using the standard equipment found in any molecular biology laboratory, and can be easily adapted for use on both genomic DNA and bisulfite DNA applications. However, in preparing bisulfite libraries for sequencing, the use of ExoSAP-IT is not recommended due to potential off-target nuclease effects which may impact downstream methylation analysis.