Global DNA methylation patterns in Barrett’s esophagus, dysplastic Barrett’s, and esophageal adenocarcinoma are associated with BMI, gender, and tobacco use
Clinical Epigeneticsvolume 8, Article number: 111 (2016)
The Erratum to this article has been published in Clinical Epigenetics 2017 9:23
The risk of developing Barrett’s esophagus (BE) and/or esophageal adenocarcinoma (EAC) is associated with specific demographic and behavioral factors, including gender, obesity/elevated body mass index (BMI), and tobacco use. Alterations in DNA methylation, an epigenetic modification that can affect gene expression and that can be influenced by environmental factors, is frequently present in both BE and EAC and is believed to play a role in the formation of BE and its progression to EAC. It is currently unknown whether obesity or tobacco smoking influences the risk of developing BE/EAC via the induction of alterations in DNA methylation. To investigate this possibility, we assessed the genome-wide methylation status of 81 esophageal tissues, including BE, dysplastic BE, and EAC epithelia using HumanMethylation450 BeadChips (Illumina).
We found numerous differentially methylated loci in the esophagus tissues when comparing males to females, obese to lean individuals, and smokers to nonsmokers. Differences in DNA methylation between these groups were seen in a variety of functional genomic regions and both within and outside of CpG islands. Several cancer-related pathways were found to have differentially methylated genes between these comparison groups.
Our findings suggest obesity and tobacco smoking may influence DNA methylation in the esophagus and raise the possibility that these risk factors affect the development of BE, dysplastic BE, and EAC through influencing the epigenetic status of specific loci that have a biologically plausible role in cancer formation.
The incidence of esophageal adenocarcinoma (EAC) has been increasing in the USA for several decades for reasons that are not entirely clear but may be related to the increasing prevalence of risk factors such as obesity . The precursor lesion for EAC is Barrett’s esophagus (BE), a metaplastic condition where the squamous-lined esophageal mucosa is replaced by specialized intestinal mucosa. A minority of individuals with BE will develop EAC through a progression sequence in which BE transitions to BE with low-grade dysplasia (LGD), BE with high-grade dysplasia (HGD), and ultimately to EAC .
It is recognized that both genetic and epigenetic alterations arise in the esophagus during the development and progression of BE and EAC [3–5]. Epigenetic alterations, primarily in the form of hypermethylated or hypomethylated CpG dinucleotides in the DNA, have been described in BE and EAC using both candidate gene approaches and microarray-based strategies. Hypermethylation of CpGs in CpG islands in promoter regions has been associated with the repression of transcription of some genes, and hypermethylation of CpGs in gene bodies is associated with increased gene expression [6, 7]. The effects of DNA methylation on the regulation of gene expression have supported the plausibility that alterations in DNA methylation can affect disease processes in people.
Aberrant DNA methylation has been shown to occur early in the BE → dysplastic BE → EAC progression sequence . The aberrant methylation of numerous cancer-related genes, such as CDKN2A, as well as global alterations in DNA methylation has been observed in BE, and many of these epigenetic alterations are also found in dysplastic BE and EAC [8–13]. However, despite the near universal observation of altered DNA methylation in BE and EAC, the mechanisms driving aberrant DNA methylation in the esophagus, as in most other pre-neoplastic and neoplastic tissues, remain elusive.
The risk of developing BE and/or EAC is associated with specific demographic and behavioral factors, including obesity/elevated body mass index (BMI) and tobacco use [14, 15]. Numerous mechanisms through which these factors may affect BE and/or EAC formation have been proposed [16, 17]; however, no assessment of effects on the epigenome in the esophagus has been made to date. There is evidence that certain environmental, behavioral, and demographic factors can influence the epigenetic state, which suggests that the behavioral factors associated with BE and EAC may act by inducing alterations in the methylation status of DNA . For example, alterations in the methylation status of CpG islands in the promoter regions of genes implicated in obesity, appetite control, and metabolism have been shown to occur in DNA isolated from blood and breast tissue of obese compared to lean individuals [19–22]. Tobacco smoking, meanwhile, has been associated with alterations in DNA methylation of multiple cancer-related genes in studies focused on single candidate genes as well as in genome-wide methylation studies of prostate cancer, the bronchial epithelium, and peripheral blood mononuclear cells [23–26].
These observations led us to use HumanMethylation450 (HM450) BeadChips to evaluate epigenome-wide patterns of DNA methylation in a collection of human esophageal tissue samples, including cases of Barrett’s esophagus (BE), Barrett’s with low- or high-grade dysplasia, and esophageal adenocarcinoma (EAC). We were interested in determining whether BMI, tobacco smoking, and/or gender were associated with increased or decreased DNA methylation at specific CpG dinucleotides or in particular genomic regions, which would support a possible functional role in the pathogenesis of EAC. We also focused on whether epigenetic alterations linked with these demographic features associated with particular molecular or cancer-related pathways in order to assess for possible mechanisms through which alterations in the DNA methylation status may be involved in the formation of BE and/or EAC.
Differences in the methylation status of genes in obesity-related pathways are associated with BMI status
Obesity has been consistently associated with an increased risk for developing both BE and EAC, yet little is known about the mechanisms involved in this elevated risk [27, 28]. While it is likely that both somatic genetic and epigenetic alterations play a role in the pathogenesis of BE and EAC, there is currently very little information about the relationship between Barrett’s esophagus and esophageal adenocarcinoma, obesity, and aberrant DNA methylation. From the 81 samples we analyzed on the HM450 array, body mass index (BMI) data were available for 46 cases, including 15 BE, 14 LGD, nine HGD, and eight EAC cases. We classified each of these samples as arising in the setting of either high BMI (BMI > 30) or low BMI (BMI ≤ 30). For female patients (N = 7), there were three in the low BMI and four in the high BMI groups. First, we determined whether the BE samples from the high BMI group (N = 11) had global DNA methylation alterations that were more closely related to HGD and/or EAC cases compared to the BE samples from study subjects with low BMI (N = 4). We found that high and low BMI BE cases tend to cluster together and that the high BMI BE cases did not appear to be more related to HGD/EAC than the low BMI BE cases (data not shown).
Next, we assessed for differentially methylated loci (DML) that varied between the combined esophageal tissue samples (BE, LGD, HGD, EAC) from individuals with high vs. low BMI. Using criteria for DML of a p value <0.001 and Δβ between high BMI and low BMI > 0.10, we found a total of 974 DML between the high and low BMI groups, including 226, 471, and 277 DML located in promoter, intragenic, and intergenic regions, respectively. A dendrogram depicting the DML between high and low BMI patients is shown in Fig. 1. One hundred and eighty-two (182) DML were located in CpG islands and 376 were located in CpG island shores (within 2 kb of a transcription start site ). We also found 352 DML (36.1 % of the total 974 DML) that were cancer associated, which we defined as loci that were differentially methylated between the normal squamous (SQ; N = 12) and EAC (N = 24) cases on the HM450 array. In general, the high BMI cases showed increased methylation at the DML, with 872 out of 974 DML (89.5 %) demonstrating elevated methylation in high vs. low BMI cases. The DML with the greatest statistical significance (p < 5 × 10−6) associated with BMI are shown in Table 1.
We also evaluated the association of BMI with tissue DNA methylation in the separate histologic types of esophageal tissues (e.g., BE, BE with LGD, BE with HGD, and EAC). We compared methylation in the high BMI (N = 4) vs. low BMI (N = 11) BE cases, the high BMI (N = 7) vs. low BMI (N = 7) LGD cases, and the high BMI (N = 9) vs. low BMI (N = 8) HGD/EAC cases. Table 2 summarizes the DML found when comparing these groups. The methylation status of the high compared to low BMI BE cases with respect to genomic regions and CpG island location is shown in Fig. 2. In general, in the BE cases, DML located in promoters and CpG islands were hypermethylated in high BMI vs. low BMI cases, whereas DML located elsewhere were hypomethylated in high BMI vs. low BMI cases. In contrast to this, DML in the HGD/EAC cases were hypermethylated at all functional regions as well as CpG island shores, shelves , and open seas in the high BMI vs. low BMI cases but not at CpG islands (Fig. 2).
We also looked to see if any of the DML between the high and low BMI BE cases overlapped with any of the DML when comparing BE to EAC, in order to determine if methylation alterations in obese individuals with BE might be associated with progression to HGD/EAC. We did find nine probes that overlapped between these groups, including those targeting the genes HLA-DPA1, TBR1, OSR2, TMEM63A, CD300E, and UBD/FAT10. UBD/FAT10, which we found to be hypomethylated in high BMI BE patients, is of interest as this gene has been shown to be overexpressed in hepatocellular carcinoma (HCC) and is thought to modulate the β-catenin/TCF4 pathway and drive HCC invasion and metastasis .
Because of the potential for DNA methylation alterations to modify gene expression, we next assessed the methylation status of CpGs located in genes associated with signaling pathways and biological mediators implicated in obesity-associated cancers [17, 32, 33] in the esophageal tissues from the subjects with low vs. high BMI. With regard to the insulin and insulin growth factor 1 (IGF-1) related pathways, we observed increased methylation of IGFBP1 (average beta = 0.11 in low BMI cases and 0.27 in high BMI cases) and IRS2 (average beta = 0.11 in low BMI cases and 0.36 in high BMI cases) in the high BMI compared to low BMI BE cases. Both genes were hypermethylated in the high BMI cases in a CpG island located within exon 1. Unlike with BE cases, genes of the insulin or IGF-1 pathways did not show altered methylation in high vs. low BMI cases in the LGD, HGD, or EAC tissue sets. We also examined molecular pathways associated with adipose inflammation, which has been shown to mediate obesity-related cancer  and found the proinflammatory gene IL-1β (IL1B) to be hypermethylated in high vs. low BMI cases when we assessed the combined esophageal tissue sets. We also found hypermethylation of IL1B in the HGD/EAC cases from high BMI subjects. For the combined cases, the average beta was 0.25 (SD = 0.10, 95 % CI = 0.21–0.30) in low BMI cases and 0.35 (SD = 0.12, 95 % CI = 0.30–0.41) in high BMI cases and for the HGD/EAC cases, average beta was 0.20 (SD = 0.08, 95 % CI = 0.12–0.27) in low BMI cases and 0.38 (SD = 0.11, 95 % CI = 0.30–0.47) in high BMI cases. Of interest, adiponectin and leptin have also been implicated in obesity-associated cancer [34, 35]; however, we did not observe any differences in the DNA methylation status of genes involved in leptin or adiponectin pathways in any of the esophageal tissue sets in the high vs. low BMI subjects.
There are numerous differentially methylated regions (DMR) between individuals with high and low BMI in esophageal tissues
The analysis described above was focused on the methylation status of individual CpG dinucleotides located in promoters, gene bodies, and intergenic regions. In light of recent studies suggesting gene expression changes are highly correlated with the aberrant methylation of large regions of DNA, called differentially methylated regions (genomic areas where numerous contiguous CpGs demonstrate significant concordant methylation alterations) [36, 37], we next assessed for differentially methylated regions (DMR) in the esophageal tissue samples from the low vs. high BMI subjects. Among the BE cases, there were DMR in 10 genes that differed between the high and low BMI groups (FWER < 0.10, Δβ > 0.10, and at least two contiguous CpG dinucleotides differentially methylated). Examples of two of these genes, TFAP2C and DIP2C, are shown in Fig. 3. Among the HGD/EAC cases, 31 DMR were identified using the same criteria, including regions in the genes ZNF790 and SIM2 (Fig. 3). We did not find any DMR within prominent genes in the insulin, IGF-1, TNF-α, or leptin pathways.
A comparison of genes showing differential methylation between high vs. low BMI cases demonstrates the involvement of cancer-related pathways and gene sets
We used the NCI Pathway Interaction Database (NCI-PID), Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and the list of Gene Ontology (GO) terms to identify biological processes or pathways that were over- or under-represented based on genes containing DML between the esophageal tissue sets in the subjects with either high or low BMI status. As mentioned previously, we defined “cancer-associated” probes as those that were differentially methylated between EAC and SQ cases on the array.
Among the BE cases, we found one NCI-PID pathway, “direct p53 effectors”, which includes the differentially methylated gene RDX from our dataset, associated with methylation differences between high and low BMI groups. There were 13 KEGG pathways (including “cell adhesion molecules”) and 77 GO terms (including “response to growth hormone” and “biological adhesion”) that were represented in the differentially methylated genes in the BE samples from the high vs. low BMI subjects. The list of GO terms is shown in Additional file 1: Table S2.
With respect to the HGD/EAC cases, there were no NCI-PID pathways that were significantly associated with methylation differences between high and low BMI status after restricting our analysis to only cancer-related genes. There was one KEGG pathway (“Wnt signaling”) and 87 GO terms (such as “tissue morphogenesis” and “response to TGF-beta”) differentially methylated between HGD/EAC cases from subjects with high BMI vs. low BMI (p value <0.05) (Additional file 2: Table S3).
Gender-related differences in DNA methylation in esophageal tissues
Little is known about gender-specific variations in DNA methylation in most tissues, including the esophagus. Previous studies have shown that repetitive elements and specific CpG dinucleotides isolated from blood samples demonstrate modestly increased methylation in males compared to females [38, 39]. Another study of four candidate genes in colorectal adenocarcinoma cells demonstrated that males had increased methylation of MTHFR, CALCA, and MGMT compared to females .
To the best of our knowledge, a genome-wide analysis of gender differences in DNA methylation in the esophagus has not been previously reported. Using HM450 array analysis of BE, HGD, and EAC esophageal samples from 118 males and 23 females, we found numerous CpG sites that were differentially methylated between the genders after excluding probes on the X and Y chromosomes and after accounting for differences in the age between the men and women in our study. When we combined the BE, HGD, and EAC cases, there were 1092 DML, including 369, 421, and 402 DML located in promoter, intragenic, and intergenic regions, respectively. From this list, there were 402 DML where the mean beta value difference between males and females was >0.10 and p value was <0.001. These DML were associated with CpGs in genes such as DUSP22, a regulator of estrogen receptor alpha mediated signaling, FRG1B, which is involved in pre-mRNA splicing, and CGREF1, which mediates cell-cell adhesion in a calcium-dependent manner. Of these 402 DML, 327 (81.3 %) were more highly methylated in females. The DML with the greatest statistical significance (p < 5 × 10−6) between males and females are listed in Table 3. Of interest, half of the top DML were located in CpG islands.
Tobacco use is associated with DNA hypermethylation in the esophagus
Tobacco smoking, which is a well-known risk factor for Barrett’s esophagus and EAC, has been associated with alterations in DNA methylation in peripheral blood lymphocytes [26, 41]. However, little is known about the relationship between smoking and DNA methylation alterations in esophageal tissues, including BE and EAC. To investigate this further, we assessed the relationship between smoking and aberrant DNA methylation in samples from subjects for which we had data on tobacco use. We divided cases into “smokers” (which included both current and former smokers) and “nonsmokers;” we did not further segregate smokers by current smoking status, pack-years, etc. due to the relatively small number of cases available. We first compared BE nonsmokers (N = 7) to BE smokers (N = 9) using principal component analysis (PCA) to determine whether methylation patterns of BE smokers more closely resembled the patterns we observed in LGD, HGD, and/or EAC cases compared to BE nonsmokers. When we examined the 1000 loci with the most variable methylation between groups, we did not find the BE smokers were grouped more closely with LGD/HGD/EAC than BE nonsmokers (data not shown).
Next, we evaluated 54 esophageal samples of various histologic types (BE, LGD, HGD, and EAC) for global alterations in DNA methylation associated with tobacco smoking. After controlling for differences associated with the histological diagnosis (BE, LGD, or HGD/EAC), we found a total of 256 DML between the smokers (N = 40) and nonsmokers (N = 14) (Δβ > 0.10, p < 0.001). Heatmaps depicting the DML between smokers and nonsmokers are shown in Fig. 4.
These DML included 98, 40, and 118 loci located in promoter, intragenic, and intergenic regions, respectively. Two hundred forty-two (242) of the 256 DML (94.5 %) were more highly methylated in smokers compared to nonsmokers, and 105 of the 256 DML (41.0 %) affected cancer-associated genes, as based on the criteria described above. The DML with the greatest statistical significance (p < 1 × 10−4) associated with smoking are shown in Table 4.
We also evaluated the association of tobacco use with DNA methylation in the separate esophageal tissue types (i.e., BE, LGD, HGD, and EAC). We assessed for DML in the BE smokers (N = 9) vs. BE nonsmokers (N = 7) and in the HGD/EAC smokers (N = 19) vs. HGD/EAC nonsmokers (N = 7) while controlling for age differences. We were not able to compare the LGD cases as all samples were from smokers. Table 5 summarizes the DML we identified for these comparisons and shows the functional genomic locations of the loci when comparing these groups. The methylation status of the BE and HGD/EAC tissues from smokers compared to nonsmokers with respect to genomic regions and CpG island location is shown in Fig. 5. In both BE and HGD/EAC cases, the DML from smokers showed much higher methylation in all genomic regions analyzed (Fig. 5).
There are numerous differentially methylated regions (DMR) in esophageal tissues based on tobacco use status
As with the BMI-based comparison described above, we were interested in extending our analysis of differential DNA methylation between smokers and nonsmokers to include differentially methylated regions in addition to DML, which are single CpG sites. Among the BE cases, there were DMR found involving 13 genes when comparing smokers to nonsmokers (FWER < 0.10, Δβ > 0.10, and at least two contiguous CpG dinucleotides differentially methylated). These DMR were located within the genes TNXB and HOXA4, which are notable because TNXB is a member of the tenascin family and regulates cell-extracellular matrix interactions [42, 43] and HOXA4 is a transcription factor previously shown to inhibit cell motility and to be aberrantly methylated in acute myeloid leukemia [44, 45] (Fig. 6). TNXB is normally more highly expressed in BE tissues compared to normal squamous esophagus (fold change = 3.39) whereas HOXA4 has not been shown to be differentially expressed in BE vs. normal esophagus (expression data obtained from www.oncomine.org). Among the HGD/EAC cases, we identified 29 DMR, including areas with altered methylation in the genes GFI1, which is a transcriptional repressor implicated in the regulation of p53 activity and Notch signaling [46, 47] and CLDN11, a cell adhesion protein involved in cell migration that is commonly altered in cancer  (Fig. 6). Normally, both GFI1 and CLDN11 have been shown to be more highly expressed in EAC tissues vs. normal esophagus (fold changes = 1.30 and 1.11–3.39, respectively; www.oncomine.org).
Differences in esophageal DNA methylation between smokers and nonsmokers are associated with several cancer-related pathways and gene sets
We were interested to see which molecular and cancer-related pathways were associated with the epigenetic differences in the BE and EAC tissues from smokers as compared to nonsmokers. As with the BMI cases, we restricted our NCI Pathway Interaction Database (NCI-PID) analysis to only those DML that we considered to be “cancer related” to improve the likelihood these pathways would contain biologically plausible mechanisms involved in smoking-related BE and/or EAC formation.
Analysis of BE cases alone did not identify any NCI-PID pathways that were differentially methylated in BE smokers vs. nonsmokers. However, there was 1 KEGG pathway (“type 1 diabetes mellitus”) and 20 GO terms (including “positive regulation of mismatch repair” and “enteric smooth muscle cell differentiation”) that were differentially represented between the BE samples from smokers vs. nonsmokers (p < 0.05) (Additional file 3: Table S4).
When we compared DNA methylation in the HGD/EAC tissues of smokers and nonsmokers, we found two NCI-PID pathways associated with alterations in DNA methylation and smoking (FDR ≤ 0.05), including the “neurotrophic factor-mediated Trk receptor signaling” and “SHP2 signaling” pathways. The differentially methylated genes NTRK2 and NTRK3 were notable affected members of both of these pathways. There were no KEGG pathways but there were 217 GO terms (such as “localization of cell” and “regulation of cell migration”) that were differentially represented (Additional file 4: Table S5).
Genetic and epigenetic alterations are commonly found in BE and EAC and likely play a prominent role in driving the initiation and progression of BE to EAC. It is also well known that a variety of environmental factors associate with the risk of developing BE and/or EAC. Thus, we assessed the relationship between DNA methylation in the esophagus and known risk factors for BE and EAC using a genome-wide methylation platform. We also sought to describe the epigenetic differences between males and females in esophageal tissues in light of the known differences in BE and EAC incidence in men vs. women. With respect to demographic and behavioral variables, we were particularly interested in the correlation of BMI and tobacco use with DNA methylation since both are well-established risk factors for BE and EAC.
We assessed the methylation status of more than 485,000 CpG sites located in 99 % of the RefSeq genes in 81 esophageal tissues representative of the stages of esophageal adenocarcinoma development (BE, BE + LGD, BE + HGD, EAC). The annotation of array probes permitted us to determine whether differentially methylated loci were located in specific types of genomic regions (promoter, gene body, or intergenic) and to determine the relationship of differentially methylated loci (DML) to CpG islands (CpG island, shore, shelf, or open sea). Our analysis of the regions outside of promoter-related CpG islands is notable because an understanding of methylation alterations in areas with relatively low CpG density is becoming increasingly recognized to be important in diseases such as cancer [49, 50]. It has been shown that CpG-rich regions (i.e., CpG islands) demonstrate more stable DNA methylation across tissues and cell populations whereas methylation is more dynamic in CpG shores (within 2 kb of a CpG islands) and CpG shelves (within 4 kb of a CpG island). Furthermore, the methylation status of CpG shores and shelves appears to regulate gene expression [29, 51].
We initially investigated the relationship between DNA methylation and BMI in esophageal tissues. Elevated body mass index (BMI) is an established risk factor for BE and EAC, and we demonstrated that DNA isolated from individuals with BMI > 30 was differentially methylated at nearly 1000 CpG sites in combined BE, BE with low- and high-grade dysplasia, and EAC tissues when compared to samples from individuals with a low BMI status. Interestingly, nearly 90 % of the DML showed elevated methylation in the high BMI cases, and over 36 % of the total DML were cancer related. There were 20 % more cancer-related DML in the high BMI group than we would expect by chance alone since just 16 % of the total probes on the array are “cancer related” by our criteria as previously described. In the BE cases, DML located in promoters and CpG islands tended to be hypermethylated in those with high BMI which suggests a possible association between methylation and altered gene expression in those with elevated BMI as promoter hypermethylation has been associated with gene silencing ; this remains speculative given we did not have associated gene expression data. There was also evidence of altered methylation in BE and EAC samples from obese patients when we looked at differentially methylated regions (DMR), which are genomic regions that have multiple adjacent CpG sites showing concordant methylation changes. DMR are potentially more biologically important than differentially methylated individual CpG dinucleotides because they are indicative of larger scale epigenetic alterations that might be more relevant functionally [36, 53].
We were also interested in whether the high BMI BE cases displayed methylation alterations resembling the EAC cases, our rationale being these epigenetic alterations in the obese with BE might be markers for progression to dysplasia or cancer and provide some evidence of a biological role for the genes subjected to aberrant methylation. This was not the case, however, as the high BMI cases clustered more closely with the low BMI cases, not the EAC cases.
We subjected the DML to KEGG, Gene Ontology (GO), and NCI-PID analyses to determine whether particular molecular groups or pathways were associated with the methylation changes in obese individuals with BE, dysplastic BE, or EAC. Among the BE cases, we found epigenetic alterations in the direct p53 effectors pathway in individuals with elevated BMI. This included differentially methylated loci within the RDX gene, which encode a cytoskeletal component that has been shown to inhibit metastasis in gastric cancer .
TP53, the gene for p53, is a well-known tumor suppressor gene that is frequently lost early in BE through mutation or loss of heterozygosity (LOH) . TP53 LOH has been shown to identify a subset of BE patients who are at risk for progression to EAC [56, 57]. The finding of differential methylation involving the p53 pathway in BE from subjects with high vs. low BMI suggests a relationship between obesity and DNA methylation of cancer-related genes in the esophagus. Similar results have been found in other studies comparing methylation in obese to lean individuals. In a recent study of 345 breast cancer cases, the majority (87 %) of CpG sites analyzed showed elevated methylation in obese patients, particularly in estrogen receptor-positive tumors. Obesity was associated with the aberrant methylation of cancer-related genes involved with the immune response, cell growth, and DNA repair . Several prior studies have compared DNA methylation in whole blood or peripheral blood leukocytes among obese and non-obese individuals [58–60]. In two of these studies, the gene HIF3A was found to be hypermethylated in the blood cells and adipose tissue of obese adults, suggesting perturbation of the hypoxia inducible transcription factor pathway in those with elevated BMI.
We were also interested in determining if there were genome-wide differences in esophageal DNA methylation between males and females. Previously, when we used a genome-wide approach to compare methylation in the normal colon between males and females, we found 82 DML between the groups, with females showing increased methylation at 69.5 % of the differentially methylated CpGs . In the present study, we found 402 DML after controlling for age and histology, with 81.3 % showing higher methylation in females. Other studies have shown differences in autosomal DNA methylation by gender in the brain, saliva, and blood [9, 62, 63]. These results suggest that DNA methylation might function in the differentiation or maintenance of sexual dimorphism. An understanding of tissue-specific gender differences is also important in terms of understanding the role of environmental, behavioral, and demographic factors on alterations in DNA methylation in order to appropriately account for potentially confounding effects of gender .
Tobacco smoking is another known risk factor for developing Barrett’s esophagus and esophageal adenocarcinoma . The mechanisms accounting for this risk are only partly understood and are believed to involve carcinogen-mediated mutations. Cigarette smoke contains multiple carcinogens which likely exert their effects via the induction of DNA adducts, aberrant DNA methylation and mutation, and chromosomal translocation [65, 66]. In order to define the association between tobacco use and aberrant DNA methylation in BE/EAC, we analyzed 54 esophageal samples of various histological types for global alterations in DNA methylation associated with tobacco smoking. We found 256 DML in these tissues between smokers and nonsmokers. Ninety-five percent (95 %) of these DML showed elevated methylation in the smoker group and 41.0 % were cancer related, which is 25 % more cancer-related DML than would be expected by chance alone.
The finding of widespread and frequent hypermethylation in BE, dysplastic BE, and EAC tissues of tobacco smokers suggests that tobacco-related epigenetic alterations may be a mechanism through which tobacco affects the development of BE and EAC. After enriching the DML (smokers vs. nonsmokers) for cancer-related genes, we found the Trk and Shp2 pathways to be differentially activated between these groups; these differences were driven by hypermethylation of the NTRK2 and NTRK3 genes in smokers. The differentially methylated NTRK2 locus, located in a promoter CpG island, demonstrated an average methylation level of 36 % in the HGD/EAC samples from smokers vs. 9 % in nonsmokers. The differentially methylated NTRK3 locus, located in the gene body, showed an average methylation of 85 % in the HGD/EAC samples of smokers compared to 62 % in nonsmokers. We previously found the aberrant methylation of NTRK3 in 60 % of colon adenomas and 67 % of colon adenocarcinomas, suggesting NTRK3 is a novel conditional tumor suppressor gene that is commonly inactivated in colorectal cancer by both epigenetic and genetic mechanisms . NTRK2 has also been shown to be hypermethylated in colon cancers as well as prostate cancer cell lines and cancers [68, 69]. With respect to esophageal cancer, NTRK2 was found to have an altered allele frequency in a group of mainly esophageal squamous cell cancers, suggesting a role in esophageal cancer susceptibility and/or development . The effect of DNA methylation on NTRK2 in BE and HGD/EAC is not clear at this time as its expression in BE or HGD is similar to normal esophagus based on publically available gene expression data, whereas NTRK3 is normally overexpressed in EAC (but not BE) vs. normal esophagus (expression level 1.03–1.79; www.oncomine.org).
We did not have mRNA expression data available for these samples to allow us to determine whether methylation alterations were associated with concordant changes in expression, which is a limitation of this study. In order to increase the likelihood that differences in methylation between the groups we studied were biologically relevant, we focused upon cancer-related pathways and pathways known to be involved in obesity and inflammation. Another potential limitation of this study in the EAC cases presumably contained a mix of cell types, including cancer cells, stromal cells, and inflammatory cells. We aimed to reduce the effects of cell heterogeneity by including only samples with >75 % cancer cells and focusing on genes with relatively large differences in methylation.
In summary, we used a microarray-based approach to determine genome-wide methylation profiles of a collection of 81 esophageal specimens, including samples of BE, dysplastic BE, and EAC DNA. With respect to gender, BMI, and tobacco use we found numerous alterations in DNA methylation involving various regions of the genome. These results suggest that obesity and tobacco smoking influence DNA methylation in the esophagus and provide novel insights into the pathways linking these risk factors to the development of BE, dysplastic BE, and EAC.
Primary tissue samples and sample preparation
DNA was extracted from formalin-fixed, paraffin-embedded (FFPE) tissue cores obtained from the Department of Pathology at University Hospitals Case Medical Center using the DNAeasy blood and tissue kit (Qiagen #69504) according to the manufacturer’s instructions with minor modifications . Protocols were approved by the institutional review board. All samples were reviewed by an expert gastrointestinal pathologist (JEW) prior to processing. The total number of samples prepared was: 21 Barrett’s esophagus (BE), 18 Barrett’s with low-grade dysplasia (BE + LGD), 18 Barrett’s with high-grade dysplasia (BE + HGD), and 24 with esophageal adenocarcinoma (EAC) (Additional file 5: Table S1). We also analyzed 12 cases of esophageal squamous epithelia (SQ) and compared methylation of this sample group to the EAC group to generate a list of “cancer-associated” loci.
Epithelial cell layers were identified and subsequently microdissected from glass slides. For the EAC cases, at least 75 % of each sample contained cancer in order to minimize methylation differences that might be due to cellular heterogeneity. After extraction, the DNA concentration was determined using the Quant-iT PicoGreen dsDNA assay kit (Invitrogen/Life Technologies, #P7589), and DNA quality was confirmed using the Illumina FFPE QC kit (Illumina, #WG-321-1001). Next, a total of 250 ng of each sample was sodium bisulfite converted using the EZ DNA methylation kit following the manufacturer’s protocol (ZymoResearch, #D5002), and then the DNA samples were treated with the Infinium HD FFPE DNA restore kit to repair any degraded DNA (Illumina, #WG-321-1002). Bisulfite-converted, restored DNA was submitted to the Genomics Core at the Fred Hutchinson Cancer Research Center (FHCRC) for processing, application, and scanning on the HumanMethylation450 (HM450) BeadChip following the manufacturer’s instructions (Illumina #WG-314-1003; http://www.Illumina.com).
Genome-wide methylation arrays
HM450 BeadChips were used to analyze patterns of DNA methylation in 81 of the esophageal samples listed above. We followed our previously validated protocols for data filtering, normalization, and differential methylation analysis [61, 72] with the following modifications or clarifications: probes with detection p value >0.05, probes on the X chromosome, and probes containing at least one SNP with low minor allele frequency (MAF = 0) in the probe body were filtered out. After filtering, a total of 453,444 probes were available for downstream analysis. The ComBat algorithm was used to correct known batch effects across the three different microarray experiments while retaining the expected variation between the different histological tissue types [73, 74]. Data was analyzed using both “β values,” where 0.0 is equivalent to 0 % methylation and 1.0 is equivalent to 100 % and “M values” which are logarithmic scores similar to those used in gene expression microarrays. We performed clustering analysis using the 3000 most highly variable loci when considering all BE, LGD, HGD, and EAC cases assessed using the HM450 array. We used the limma and minfi Bioconductor packages to compute a refined F-statistic to quantify the difference in DNA methylation based on a probe’s M-value between sample types. We used a false discovery rate (FDR) q value to determine the significance of differentially methylated loci (DML) and considered loci to be differentially methylated if q < 1 × 10−5 . Cancer-associated loci were those that showed differential methylation when comparing EAC and squamous (SQ) samples (q < 0.001).
The relationships of differentially methylated CpG dinucleotides to CpG islands were determined using the HM450 array annotation along with published definitions [29, 30]. Thus, a CpG island shore is the region located 0–2 kb from a CpG island, a CpG island shelf is located 2–4 kb from a CpG island, and an “open sea” is located >4 kb from a CpG island.
Differentially methylated regions (DMR) were calculated using the Bumphunter method as described by Jaffe et al. . To be considered a DMR, regions had to contain at least two contiguous CpG dinucleotides that were differentially methylated with family-wise error rate (FWER) < 0.10 and Δβ > 0.10.
The actual vs. expected distributions of BMI and smoking-associated DML with respect to genomic location and CpG island were calculated using Pearson chi-square tests. In other words, we compared the distribution of differentially hypo- or hypermethylated loci and the distribution of all of the probes on the HM450 array with respect to genomic location/CpG island location to determine whether particular regions were enriched in hypo/hypermethylated loci.
Gene pathway analysis
Pathway enrichment analysis of significantly differentially methylated genes between any two sample groups was performed using pathway definitions derived from the NCI Pathway Interaction Database (NCI-PID), a curated collection of known biomolecular interactions and key signaling pathways associated with cancer . The enrichment analysis was performed using the hypergeometric test to evaluate if genes belonging to a given pathway were enriched among the significantly differentially methylated loci. We elected to increase the possibility that altered molecular pathways would be biologically relevant by restricting our NCI-PID analysis to only cancer-associated DML. Hochberg FDR methodology and pathways with FDR ≤ 0.05 were considered significantly methylated genes. This was followed by followed by assessment of false discovery rate using the Benjamini Hochberg correction . Genes with multiple differentially methylated probes were excluded if the methylation state of any probe was inconsistent between comparison groups. Gene set enrichment analysis was performed using genes from differentially methylated groups to identify affected Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways across different comparison groups using hypergeometric testing provided by the GOstats software . Probes whose target genes were not annotated to at least one GO term in the biological process ontology were filtered out. A gene set was considered altered if its number of differentially methylated CpG sites was higher or lower than expected using a p value <0.05.
Body mass index
Differentially methylated loci/locus
Differentially methylated region
Kyoto Encyclopedia of Genes and Genomes
National Cancer Institute Pathway Interaction Database
Long E, Beales IL. The role of obesity in oesophageal cancer development. Ther Adv Gastroenterol. 2014;7:247–68.
Spechler SJ. Clinical practice. Barrett’s esophagus. N Engl J Med. 2002;346:836–42.
Werner M, Mueller J, Walch A, Hofler H. The molecular pathology of Barrett’s esophagus. Histol Histopathol. 1999;14:553–9.
Flejou JF. Barrett’s oesophagus: from metaplasia to dysplasia and cancer. Gut. 2005;54 Suppl 1:i6–i12.
Reid BJ, Levine DS, Longton G, Blount PL, Rabinovitch PS. Predictors of progression to cancer in Barrett’s esophagus: baseline histology and flow cytometry identify low- and high-risk patient subsets. Am J Gastroenterol. 2000;95:1669–76.
Baylin SB, Jones PA. A decade of exploring the cancer epigenome—biological and translational implications. Nat Rev Cancer. 2011;11:726–34.
Yang X, Han H, De Carvalho DD, Lay FD, Jones PA, Liang G. Gene body methylation can alter gene expression and is a therapeutic target in cancer. Cancer Cell. 2014;26:577–90.
Kaz AM, Wong CJ, Luo Y, Virgin JB, Washington MK, Willis JE, Leidner RS, Chak A, Grady WM. DNA methylation profiling in Barrett’s esophagus and esophageal adenocarcinoma reveals unique methylation signatures and molecular subclasses. Epigenetics. 2011;6:1403–12.
Xu E, Gu J, Hawk ET, Wang KK, Lai M, Huang M, Ajani J, Wu X. Genome-wide methylation analysis shows similar patterns in Barrett’s esophagus and esophageal adenocarcinoma. Carcinogenesis. 2013;34:2750–6.
Alvi MA, Liu X, O’Donovan M, Newton R, Wernisch L, Shannon NB, Shariff K, di Pietro M, Bergman JJ, Ragunath K, et al. DNA methylation as an adjunct to histopathology to detect prevalent, inconspicuous dysplasia and early-stage neoplasia in Barrett’s esophagus. Clin Cancer Res. 2013;19:878–88.
Wu W, Bhagat TD, Yang X, Song JH, Cheng Y, Agarwal R, Abraham JM, Ibrahim S, Bartenstein M, Hussain Z, et al. Hypomethylation of noncoding DNA regions and overexpression of the long noncoding RNA, AFAP1-AS1, in Barrett’s esophagus and esophageal adenocarcinoma. Gastroenterology. 2013;144(956–966):e954.
Kaz AM, Grady WM. Epigenetic biomarkers in esophageal cancer. Cancer Lett. 2012;in press.
Eads C, Lord R, Wickramasinghe K, Long T, Kurumboor S, Bernstein L, Peters J, DeMeester S, DeMeester T, Skinner K. Epigenetic patterns in the progression of esophageal adenocarcinoma. Cancer Res. 2001;61:3410–8.
Spechler SJ, Sharma P, Souza RF, Inadomi JM, Shaheen NJ. American Gastroenterological Association medical position statement on the management of Barrett’s esophagus. Gastroenterology. 2011;140:1084–91.
Rubenstein JH, Thrift AP. Risk factors and populations at risk: selection of patients for screening for Barrett’s oesophagus. Best Pract Res Clin Gastroenterol. 2015;29:41–50.
Kendall BJ, Thrift AP. Unravelling the riddle of gastroesophageal reflux disease, obesity, and Barrett’s esophagus. Clin Gastroenterol Hepatol. 2015; in press.
Hardikar S, Song X, Risques RA, Montine TJ, Duggan C, Blount PL, Reid BJ, Anderson GL, Kratz M, White E, et al. Obesity and inflammation markers in relation to leukocyte telomere length in a cross-sectional study of persons with Barrett’s esophagus. BMC Obesity. 2015;2:32.
Alegria-Torres JA, Baccarelli A, Bollati V. Epigenetics and lifestyle. Epigenomics. 2011;3:267–77.
Drummond EM, Gibney ER. Epigenetic regulation in obesity. Curr Opin Clin Nutr Metab Care. 2013;16:392–7.
van Dijk SJ, Molloy PL, Varinli H, Morrison JL, Muhlhausler BS, Members of Epi, S. Epigenetics and human obesity. Int J Obes. 2015;39:85–97.
Ding X, Zheng D, Fan C, Liu Z, Dong H, Lu Y, Qi K. Genome-wide screen of DNA methylation identifies novel markers in childhood obesity. Gene. 2015;566:74–83.
Hair BY, Troester MA, Edmiston SN, Parrish EA, Robinson WR, Wu MC, Olshan AF, Swift-Scanlan T, Conway K. Body mass index is associated with gene methylation in estrogen receptor-positive breast tumors. Cancer Epidemiol Biomarkers Prev. 2015;24:580–6.
Enokida H, Shiina H, Urakami S, Terashima M, Ogishima T, Li LC, Kawahara M, Nakagawa M, Kane CJ, Carroll PR, et al. Smoking influences aberrant CpG hypermethylation of multiple genes in human prostate carcinoma. Cancer. 2006;106:79–86.
Belinsky SA, Palmisano WA, Gilliland FD, Crooks LA, Divine KK, Winters SA, Grimes MJ, Harms HJ, Tellez CS, Smith TM, et al. Aberrant promoter methylation in bronchial epithelium and sputum from current and former smokers. Cancer Res. 2002;62:2370–7.
Wan ES, Qiu W, Baccarelli A, Carey VJ, Bacherman H, Rennard SI, Agusti A, Anderson W, Lomas DA, Demeo DL. Cigarette smoking behaviors and time since quitting are associated with differential DNA methylation across the human genome. Hum Mol Genet. 2012;21:3073–82.
Zeilinger S, Kuhnel B, Klopp N, Baurecht H, Kleinschmidt A, Gieger C, Weidinger S, Lattka E, Adamski J, Peters A, et al. Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS One. 2013;8:e63812.
Chandar AK, Iyer PG. Role of obesity in the pathogenesis and progression of Barrett’s esophagus. Gastroenterol Clin N Am. 2015;44:249–64.
Rubenstein JH, Shaheen NJ. Epidemiology, diagnosis, and management of esophageal adenocarcinoma. Gastroenterology. 2015;149(302–317):e301.
Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, et al. The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009;41:178–86.
Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, et al. High density DNA methylation array with single CpG site resolution. Genomics. 2011;98:288–95.
Yuan R, Wang K, Hu J, Yan C, Li M, Yu X, Liu X, Lei J, Guo W, Wu L, et al. Ubiquitin-like protein FAT10 promotes the invasion and metastasis of hepatocellular carcinoma by modifying beta-catenin degradation. Cancer Res. 2014;74:5287–300.
Howe LR, Subbaramaiah K, Hudis CA, Dannenberg AJ. Molecular pathways: adipose inflammation as a mediator of obesity-associated cancer. Clin Cancer Res. 2013;19:6074–83.
Zeng L, Perks CM, Holly JM. IGFBP-2/PTEN: a critical interaction for tumours and for general physiology? Growth Horm IGF Res. 2015;25:103–7.
Mauro L, Naimo GD, Ricchio E, Panno ML, Ando S. Cross-talk between adiponectin and IGF-IR in breast cancer. Front Oncol. 2015;5:157.
Jiang N, Sun R, Sun Q. Leptin signaling molecular actions and drug target in hepatocellular carcinoma. Drug Des Devel Ther. 2014;8:2295–302.
Fan H, Zhao H, Pang L, Liu L, Zhang G, Yu F, Liu T, Xu C, Xiao Y, Li X. Systematically prioritizing functional differentially methylated regions (fDMRs) by integrating multi-omics data in colorectal cancer. Sci Rep. 2015;5:12789.
Kretzmer H, Bernhart SH, Wang W, Haake A, Weniger MA, Bergmann AK, Betts MJ, Carrillo-de-Santa-Pau E, Doose G, Gutwein J. et al. DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control. Nat Genet. 2015; in press.
El-Maarri O, Becker T, Junen J, Manzoor SS, Diaz-Lacava A, Schwaab R, Wienker T, Oldenburg J. Gender specific differences in levels of DNA methylation at selected loci from human total blood: a tendency toward higher methylation levels in males. Hum Genet. 2007;122:505–14.
Fuke C, Shimabukuro M, Petronis A, Sugimoto J, Oda T, Miura K, Miyazaki T, Ogura C, Okazaki Y, Jinno Y. Age related changes in 5-methylcytosine content in human peripheral leukocytes and placentas: an HPLC-based study. Ann Hum Genet. 2004;68:196–204.
Sarter B, Long TI, Tsong WH, Koh WP, Yu MC, Laird PW. Sex differential in methylation patterns of selected genes in Singapore Chinese. Hum Genet. 2005;117:402–3.
Breitling LP, Yang R, Korn B, Burwinkel B, Brenner H. Tobacco-smoking-related differential DNA methylation: 27 K discovery and replication. Am J Hum Genet. 2011;88:450–7.
Levy P, Ripoche H, Laurendeau I, Lazar V, Ortonne N, Parfait B, Leroy K, Wechsler J, Salmon I, Wolkenstein P, et al. Microarray-based identification of tenascin C and tenascin XB, genes possibly involved in tumorigenesis associated with neurofibromatosis type 1. Clin Cancer Res. 2007;13:398–407.
Hsu MK, Wu IC, Cheng CC, Su JL, Hsieh CH, Lin YS, Chen FC. Triple-layer dissection of the lung adenocarcinoma transcriptome—regulation at the gene, transcript, and exon levels. Oncotarget. 2015; in press.
Klausen C, Leung PC, Auersperg N. Cell motility and spreading are suppressed by HOXA4 in ovarian cancer cells: possible involvement of beta1 integrin. Mol Cancer Res. 2009;7:1425–37.
Musialik E, Bujko M, Kober P, Grygorowicz MA, Libura M, Przestrzelska M, Juszczynski P, Borg K, Florek I, Jakobczyk M, et al. Promoter DNA methylation and expression levels of HOXA4, HOXA5 and MEIS1 in acute myeloid leukemia. Mol Med Rep. 2015;11:3948–54.
Phelan JD, Saba I, Zeng H, Kosan C, Messer MS, Olsson HA, Fraszczak J, Hildeman DA, Aronow BJ, Moroy T, et al. Growth factor independent-1 maintains Notch1-dependent transcriptional programming of lymphoid precursors. PLoS Genet. 2013;9:e1003713.
Khandanpour C, Moroy T. Growth factor independence 1 (Gfi1) as a regulator of p53 activity and a new therapeutical target for ALL. Oncotarget. 2013;4:374–5.
Agarwal R, Mori Y, Cheng Y, Jin Z, Olaru AV, Hamilton JP, David S, Selaru FM, Yang J, Abraham JM, et al. Silencing of claudin-11 is associated with increased invasiveness of gastric cancer cells. PLoS One. 2009;4:e8002.
Agarwal A, Polineni R, Hussein Z, Vigoda I, Bhagat TD, Bhattacharyya S, Maitra A, Verma A. Role of epigenetic alterations in the pathogenesis of Barrett’s esophagus and esophageal adenocarcinoma. Int J Clin Exp Pathol. 2012;5:382–96.
Tang JT, Wang ZH, Fang JY. Assessing the potential value of long interspersed element-1 hypomethylation in colorectal cancer: evidence from retrospective studies. OncoTargets Ther. 2015;8:3265–76.
Ziller MJ, Gu H, Muller F, Donaghey J, Tsai LT, Kohlbacher O, De Jager PL, Rosen ED, Bennett DA, Bernstein BE, et al. Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013;500:477–81.
Baylin SB, Herman JG, Graff JR, Vertino PM, Issa JP. Alterations in DNA methylation: a fundamental aspect of neoplasia. Adv Cancer Res. 1998;72:141–96.
Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, Fan G, Sun YE. DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science. 2003;302:890–3.
Tsai MM, Wang CS, Tsai CY, Chen CY, Chi HC, Tseng YH, Chung PJ, Lin YH, Chung IH, Chen CY, et al. MicroRNA-196a/-196b promote cell metastasis via negative regulation of radixin in human gastric cancer. Cancer Lett. 2014;351:222–31.
Reid BJ. Early events during neoplastic progression in Barrett’s esophagus. Cancer Biomark. 2010;9:307–24.
Reid BJ. p53 and neoplastic progression in Barrett’s esophagus. Am J Gastroenterol. 2001;96:1321–3.
Dolan K, Morris AI, Gosney JR, Field JK, Sutton R. Loss of heterozygosity on chromosome 17p predicts neoplastic progression in Barrett’s esophagus. J Gastroenterol Hepatol. 2003;18:683–9.
Feinberg AP, Irizarry RA, Fradin D, Aryee MJ, Murakami P, Aspelund T, Eiriksdottir G, Harris TB, Launer L, Gudnason V, et al. Personalized epigenomic signatures that are stable over time and covary with body mass index. Sci Transl Med. 2010;2:49ra67.
Demerath EW, Guan W, Grove ML, Aslibekyan S, Mendelson M, Zhou YH, Hedman AK, Sandling JK, Li LA, Irvin MR, et al. Epigenome-wide association study (EWAS) of BMI, BMI change and waist circumference in African American adults identifies multiple replicated loci. Hum Mol Genet. 2015;24:4464–79.
Dick KJ, Nelson CP, Tsaprouni L, Sandling JK, Aissi D, Wahl S, Meduri E, Morange PE, Gagnon F, Grallert H, et al. DNA methylation and body-mass index: a genome-wide analysis. Lancet. 2014;383:1990–8.
Kaz AM, Wong CJ, Dzieciatkowski S, Luo Y, Schoen RE, Grady WM. Patterns of DNA methylation in the normal colon vary by anatomical location, gender, and age. Epigenetics. 2014;9:492–502.
Inoshita M, Numata S, Tajima A, Kinoshita M, Umehara H, Yamamori H, Hashimoto R, Imoto I, Ohmori T. Sex differences of leukocytes DNA methylation adjusted for estimated cellular proportions. Biol Sex Differ. 2015;6:11.
Liu J, Morgan M, Hutchison K, Calhoun VD. A study of the influence of sex on genome wide methylation. PLoS One. 2010;5:e10028.
Runge TM, Abrams JA, Shaheen NJ. Epidemiology of Barrett’s esophagus and esophageal adenocarcinoma. Gastroenterol Clin N Am. 2015;44:203–31.
Toh Y, Oki E, Ohgaki K, Sakamoto Y, Ito S, Egashira A, Saeki H, Kakeji Y, Morita M, Sakaguchi Y, et al. Alcohol drinking, cigarette smoking, and the development of squamous cell carcinoma of the esophagus: molecular mechanisms of carcinogenesis. Int J Clin Oncol. 2010;15:135–44.
Xu XC. Risk factors and gene expression in esophageal cancer. Methods Mol Biol. 2009;471:335–60.
Luo Y, Kaz AM, Kanngurn S, Welsch P, Morris SM, Wang J, Lutterbaugh JD, Markowitz SD, Grady WM. NTRK3 is a potential tumor suppressor gene commonly inactivated by epigenetic mechanisms in colorectal cancer. PLoS Genet. 2013;9:e1003552.
Ahmed D, Danielsen SA, Aagesen TH, Bretthauer M, Thiis-Evensen E, Hoff G, Rognum TO, Nesbakken A, Lothe RA, Lind GE. A tissue-based comparative effectiveness analysis of biomarkers for early detection of colorectal tumors. Clin Transl Gastroenterol. 2012;3:e27.
Yamada Y, Toyota M, Hirokawa Y, Suzuki H, Takagi A, Matsuzaki T, Sugimura Y, Yatani R, Shiraishi T, Watanabe M. Identification of differentially methylated CpG islands in prostate cancer. Int J Cancer. 2004;112:840–5.
Chen J, Guo L, Peiffer DA, Zhou L, Chan OT, Bibikova M, Wickham-Garcia E, Lu SH, Zhan Q, Wang-Rodriguez J, et al. Genomic profiling of 766 cancer-related genes in archived esophageal normal and carcinoma tissues. Int J Cancer. 2008;122:2249–54.
Moinova H, Leidner RS, Ravi L, Lutterbaugh J, Barnholtz-Sloan JS, Chen Y, Chak A, Markowitz SD, Willis JE. Aberrant vimentin methylation is characteristic of upper gastrointestinal pathologies. Cancer Epidemiol Biomarkers Prev. 2012;21:594–600.
Luo Y, Wong CJ, Kaz AM, Dzieciatkowski S, Carter KT, Morris SM, Wang J, Willis JE, Makar KW, Ulrich CM, et al. Differences in DNA methylation signatures reveal multiple pathways of progression from adenoma to colorectal cancer. Gastroenterology. 2014;147(418–429):e418.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8:118–27.
Leek JT, Scharpf RB, Bravo HC, Simcha D, Langmead B, Johnson WE, Geman D, Baggerly K, Irizarry RA. Tackling the widespread and critical impact of batch effects in high-throughput data. Nat Rev Genet. 2010;11:733–9.
Storey JD. The positive false discovery rate: a Bayesian interpretation and the q-value. Ann Stat. 2003;31:2013–35.
Jaffe AE, Murakami P, Lee H, Leek JT, Fallin MD, Feinberg AP, Irizarry RA. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int J Epidemiol. 2012;41:200–9.
Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, Buetow KH. PID: the Pathway Interaction Database. Nucleic Acids Res. 2009;37:D674–9.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125:279–84.
Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23:257–8.
The authors would like to acknowledge the invaluable work of members of the Genomics Shared Resource at the FHCRC.
This work was supported by the National Institutes of Health (U54CA163060 to AK, AC, VV, and WG, U01CA152756, 5P30CA015704, U01CA086402, and U01CA182940 to WG, and P50CA150964 to VV); the Listwin Family Foundation to WG, and the Mercer Island Rotary Club to WG. This material is the result of work supported in part by resources from the VA Puget Sound Health Care System. The views expressed in this article are those of the authors and do not necessarily represent the views of the Department of Veterans Affairs.
Availability of data and materials
Upon acceptance of this manuscript, the dataset supporting the conclusions of this article will be available in the Gene Expression Omnibus (GEO) repository, [GSE89181].
AK designed the study, interpreted the data, and prepared the manuscript. CW analyzed the data and performed statistical analyses. VV analyzed the data and performed pathway analyses. JW identified the primary tissue samples and confirmed the pathological diagnosis. AC reviewed the demographic information and interpreted the data. WG helped design the study and interpreted the data. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Protocols were approved by the institutional review board at University Hospitals Case Western (IRB number 10-87-17). The approval for this study included a waiver of consent for the acquisition of human material as the study was deemed low risk and all personalized health information was disassociated from the samples.
An erratum to this article is available at http://dx.doi.org/10.1186/s13148-017-0324-8.
GO terms represented in DM genes comparing BE high vs. low BMI. (DOC 28 kb)
GO terms represented in DM genes comparing HGD/EAC high vs. low BMI. (DOC 28 kb)
GO terms represented in DM genes comparing BE smokers vs. nonsmokers. (DOC 28 kb)
GO terms represented in DM genes comparing HGD/EAC smokers vs. nonsmokers. (DOC 28 kb)
Samples analyzed on HM450 arrays. (DOC 122 kb)