We performed a genome-wide scan of 27,578 CpG loci covering 14,475 genes to identify differentially methylated loci (DML) in colorectal carcinoma (CRC).

Background We examined whether differences in tumor DNA methylation were associated with more aggressive hormone receptor-negative breast cancer in an ethnically diverse group of patients in the Breast Cancer Care in Chicago (BCCC) study and using data from The Cancer Genome Atlas (TCGA). Results DNA was extracted from formalin-fixed, paraffin-embedded samples on 75 patients (21 White, 31 African-American, and 23 Hispanic) (training dataset) enrolled in the BCCC. Hormone receptor status was defined as negative if tumors were negative for both estrogen and progesterone (ER/PR) receptors (N = 22/75). DNA methylation was analyzed at 1505 CpG sites within 807 gene promoters using the Illumina GoldenGate assay. Differential DNA methylation as a predictor of hormone receptor status was tested while controlling for false discovery rate and assigned to the gene closest to the respective CpG site. Next, those genes that predicted ER/PR status were validated using TCGA data with respect to DNA methylation (validation dataset), and correlations between CpG methylation and gene expression were examined. In the training dataset, 5.7 % of promoter mean methylation values (46/807) were associated with receptor status at P < 0.05; for 88 % of these (38/46), hypermethylation was associated with receptor-positive disease. Hypermethylation for FZD9, MME, BCAP31, HDAC9, PAX6, SCGB3A1, PDGFRA, IGFBP3, and PTGS2 genes most strongly predicted receptor-positive disease. Twenty-one of 24 predictor genes from the training dataset were confirmed in the validation dataset. The level of DNA methylation at 19 out 22 genes, for which gene expression data were available, was associated with gene activity. Conclusions Higher levels of promoter methylation strongly correlate with hormone receptor positive status of breast tumors. For most of the genes identified in our training dataset as ER/PR receptor status predictors, DNA methylation correlated with stable gene expression level. The predictors performed well when evaluated on independent set of samples, with different racioethnic distribution, thus providing evidence that this set of DNA methylation biomarkers will likely generalize to prospective patient samples.

Mature microRNAs (miRNAs) are a class of small non-coding RNAs involved in posttranslational gene silencing. Previous studies found that downregulation of miRNAs is a common feature observed in solid tumors, including hepatocellular carcinoma (HCC). We employed a genome-wide approach to test the hypothesis that DNA methylation alterations in miRNA host genes may cause deregulated miRNA expression in HCC. We analyzed tumor and adjacent non-tumor tissues from 62 Taiwanese HCC cases using Infinium HumanMethylation27 DNA Analysis BeadChips that include 254 CpG sites covering 110 miRNAs from 64 host genes. Expression levels of three identified miRNAs (miR-10a, miR-10b and miR-196b) were measured in a subset of 37 HCC tumor and non-tumor tissues. After Bonferroni adjustment, a total of 54 CpG sites from 27 host genes significantly differed in DNA methylation levels between tumor and adjacent non-tumor tissues with 53 sites significantly hypermethylated in tumor tissues. Among the 54 significant CpG sites, 15 sites had more than 2-fold tumor/non-tumor changes, 17 sites had differences > 10%, and 10 sites had both features [including 8 significantly hypermethylated CpG sites in the host genes of miR-10a, miR-10b and miR-196b (HOXB4, HOXD4 and HOXA9, respectively)]. Significant downregulation of miR-10a was observed in tumor compared with non-tumor tissues (0.50 vs. 1.73, p = 0.031). The concordance for HOXB4 methylation alteration and dysregulation of miR-10a was 73.5%. No significant change was observed for miR-10b expression. Unexpectedly, miR-196b was significantly upregulated in tumor compared with non-tumor tissues (p = 0.0001). These data suggest that aberrant DNA methylation may lead to dysregulation of miR-10a in HCC tumor tissues.

Background: Formalin-fixed, paraffin-embedded (FFPE) samples are a highly desirable resource for epigenetic studies, but there is no suitable platform to assay genome-wide methylation in these widely available resources. Recently, Thirlwell et al. (2010) have reported a modified ligation-based DNA repair protocol to prepare FFPE DNA for the Infinium methylation assay. In this study, we have tested the accuracy of methylation data obtained with this modification by comparing paired fresh-frozen (FF) and FFPE colon tissue (normal and tumor) from colorectal cancer patients. We report locus-specific correlation and concordance of tumor-specific differentially methylated loci (DML), both of which were not previously assessed. Methods: We used Illumina's Infinium Methylation 27K chip for 12 pairs of FF and 12 pairs of FFPE tissue from tumor and surrounding healthy tissue from the resected colon of the same individual, after repairing the FFPE DNA using Thirlwell's modified protocol. Results: For both tumor and normal tissue, overall correlation of β values between all loci in paired FF and FFPE was comparable to previous studies. Tissue storage type (FF or FFPE) was found to be the most significant source of variation rather than tissue type (normal or tumor). We found a large number of DML between FF and FFPE DNA. Using ANOVA, we also identified DML in tumor compared to normal tissue in both FF and FFPE samples, and out of the top 50 loci in both groups only 7 were common, indicating poor concordance. Likewise, while looking at the correlation of individual loci between FFPE and FF across the patients, less than 10% of loci showed strong correlation (r ≥ 0.6). Finally, we checked the effect of the ligation-based modification on the Infinium chemistry for SNP genotyping on an independent set of samples, which also showed poor performance. Conclusion: Ligation of FFPE DNA prior to the Infinium genome-wide methylation assay may detect a reasonable number of loci, but the numbers of detected loci are much fewer than in FF samples. More importantly, the concordance of DML detected between FF and FFPE DNA is suboptimal, and DML from FFPE tissues should be interpreted with great caution.

Background: We evaluated (a) the feasibility of whole genome cDNA-mediated Annealing, Selection, extension and Ligation (DASL) assay on formalin-fixed paraffin-embedded (FFPE) tissue and (b) whether similar conclusions can be drawn by examining FFPE samples as proxies for fresh frozen (FF) tissues. We used a whole genome DASL assay (addressing 18,391 genes) on a total of 72 samples from paired breast tumor and surrounding healthy tissues from both FF and FFPE samples. Results: Gene detection was very good with comparable success between the FFPE and FF samples. Reproducibility was also high (r² = 0.98); however, concordance between the two types of samples was low. Only one-third of the differentially expressed genes in tumor tissues (compared to corresponding normal) from FF samples could be detected in FFPE samples and conversely only one-fourth of the differentially expressed genes from FFPE samples could be detected in FF samples. GO-enrichment analysis, gene set enrichment analysis (GSEA) and GO-ANOVA analyses also suggested small overlap between the lead functional groups that were differentially expressed in tumor detectable by examining FFPE and FF samples. In other words, FFPE samples may not be ideal for picking individual target gene(s), but may be used to identify some of the lead functional group(s) of genes that are differentially expressed in tumor. The differentially expressed genes in breast cancer found in our study were biologically meaningful. The "cell cycle" & "cell division" related genes were up-regulated and genes related to "regulation of epithelial cell proliferation" were down-regulated. Conclusions: Gene expression experiments using the DASL assay can efficiently handle fragmentation issues in the FFPE tissues. However, formalin fixation seems to change RNA and consequently significantly alters gene expression in a number of genes which may not be uniform between tumor and normal tissues. Therefore, considerable caution needs to be taken when interpreting gene expression data from FFPE tissues, especially in relation to specific genes.

Exposure to arsenic in drinking water is a global health problem and arsenic-induced skin lesions are hallmark of chronic arsenic toxicity. We and others have reported germline genetic variations as risk factors for such skin lesions. The role of copy number variation (CNV) in the germline DNA in this regard is unknown. From a large prospectively followed-up cohort, exposed to arsenic, we randomly selected 2171 subjects without arsenic-induced skin lesions at enrollment and genotyped their whole blood DNA samples on Illumina Cyto12v2.1 SNP chips to generate DNA copy number. Participants were followed up every 2 years for a total of 8 years, especially for the development of skin lesions. In Cox regression models, each CNV segment was used as a predictor, accounting for other potential covariates, for incidence of skin lesions. The presence of genomic deletion(s) in a number of genes (OR5J2, GOLGA6L7P, APBA2, GALNTL5, VN1R31P, PHKG1P2, SGCZ, ZNF658) and lincRNA genes (RP11-76I14.1, CTC-535 M15.2, RP11-73B2.2) were associated with higher risk [HR between 1.67 (CI 1.3-2.1) and 2.15 (CI 1.5-2.9) for different CNVs] for development of skin lesions independent of gender, age, and arsenic exposure. Some deletions had stronger effect in a specific gender (ZNF658 in males, SGCZ in females) and some had stronger effect in higher arsenic exposure (lincRNA CTD-3179P9.1) suggesting a possible gene-environment interaction. This first genome-wide CNV study in a prospectively followed-up large cohort, exposed to arsenic, suggests that DNA deletion in several genes and lincRNA genes may predispose an individual to a higher risk of development of arsenic-induced skin lesions.