This article has Open Peer Review reports available.
Methylation of leukocyte DNA and ovarian cancer: relationships with disease status and outcome
© Fridley et al.; licensee BioMed Central Ltd. 2014
Received: 13 December 2013
Accepted: 16 April 2014
Published: 28 April 2014
Genome-wide interrogation of DNA methylation (DNAm) in blood-derived leukocytes has become feasible with the advent of CpG genotyping arrays. In epithelial ovarian cancer (EOC), one report found substantial DNAm differences between cases and controls; however, many of these disease-associated CpGs were attributed to differences in white blood cell type distributions.
We examined blood-based DNAm in 336 EOC cases and 398 controls; we included only high-quality CpG loci that did not show evidence of association with white blood cell type distributions to evaluate association with case status and overall survival.
Of 13,816 CpGs, no significant associations were observed with survival, although eight CpGs associated with survival at p < 10-3, including methylation within a CpG island located in the promoter region of GABRE (p = 5.38 x 10-5, HR = 0.95). In contrast, 53 CpG methylation sites were significantly associated with EOC risk (p <5 x10-6). The top association was observed for the methylation probe cg04834572 located approximately 315 kb upstream of DUSP13 (p = 1.6 x10-14). Other disease-associated CpGs included those near or within HHIP (cg14580567; p =5.6x10-11), HDAC3 (cg10414058; p = 6.3x10-12), and SCR (cg05498681; p = 4.8x10-7).
We have identified several CpGs in leukocytes that are differentially methylated by case-control status. Since a retrospective study design was used, we cannot differentiate whether DNAm was etiologic or resulting from EOC; thus, prospective studies of EOC-associated loci are the critical next step.
The role of DNA methylation (DNAm) in ovarian cancer is multi-faceted. While tumor tissue shows clear methylation patterns associating with histopathology, the role of blood-based DNAm patterns on disease etiology and outcome has been a subject of growing interest [1–4]. This includes study of variation in inherent global methylation levels, the relationship between exogenous exposures and leukocyte methylation, and the role of inherited variants on leukocyte methylation (mQTL) [5–8]. Five of the eleven confirmed ovarian cancer susceptibility variants and an endometriosis locus are located in homeobox gene clusters (HOXA, HOXB, and HOXD), homeobox related genes (HNF1B), or genes expressed in early progenitor cells (TERT) [9–13]. Thus, we hypothesize that DNAm levels in circulating systemic leukocytes of ovarian cancer cases and controls may differ, and that among cases, leukocyte methylation may vary by disease outcome.
Previous work by Teschendorff et al. (2009)  identified peripheral blood methylation signatures that predicted ovarian cancer case-control status using methylation measurements at more than 27,000 CpGs in 113 cases and 148 controls. However, as pointed out in the discussion by Teschendorff et al. and subsequently by Koestler et al.(2009)  and Houseman et al. (2012) , blood-based methylation measurements are dependent on distribution of white blood cell (leukocyte) types and the distribution of cell types is also related to disease status (i.e., confounding). Therefore, in order to minimize confounding by distribution of cell types, we performed case-control and survival analyses using 336 EOC cases and 398 controls, accounting for cell type associations to better understand the role of blood-based DNAm in ovarian cancer risk and survival.
Characteristics of study participants
Age at case diagnosis
Age at control enrollment
EOC case vital status
Follow-up time (years)
Histology, cases only
Grade, cases only
Grade 1 or 2
Grade 3 or 4
Stage, cases only
Stage I or II
Stage III or IV
Surgical debulking, cases only
Optimal (<1 cm)
Sub-optimal (>1 cm)
Ascites, cases only
DNA methylation assays and arrays
Peripheral blood (leukocytes) was used as the source of DNA. DNA was extracted from four milliliters of fresh peripheral blood using the Autogen Flexstar instrument utilizing Flexigene chemistry (salting out methodology). Blood is aliquoted for DNA extraction using an automated liquid handler with barcoding to ensure proper sample placement. Post-DNA extraction, DNA is aliquoted into a permanent storage tube utilizing an automated liquid handler with barcoding, again, to ensure proper sample placement. DNA samples are assessed for quality and concentration using a Trinean DropSense 96 spectrophotometer and DNA is then stored longterm at -80°C. The leukocyte-derived DNA (1 ug) was bisulfite modified (BSM) using the Zymo EZ96 DNA Methylation Kit (Zymo Research, Orange, CA) according to the manufacturer’s protocol. BSM DNA (250 ng) was then assayed on 96 well plates in three batches at the Mayo Clinic Molecular Genome Facility (Rochester, MN): Batch 1 used the Infinium HumanMethylation27 BeadChip on 84 cases and 91 controls, Batch 2 used this array on 172 cases and 176 controls and Batch 3 used Illumina Infinium HumanMethylation450 BeadChip on 156 cases and 157 controls. Methylation status at the target CpG sites was determined by comparing the ratio of fluorescent signal from the methylated allele to the sum from the fluorescent signal from both methylated and unmethylated alleles (i.e., the beta value).
To assess the quality of the DNAm data produced from the Illumina arrays, Centre d'Etudes du Polymorphisme Humain (CEPH) DNA, positive BSM controls (placental DNA) and negative BSM controls (whole genome amplified [WGA] DNA) were assayed within each batch. For the HumanMethylation27 BeadChips (Batch 1 and Batch 2), 9 CEPH DNA, 12 positive control DNA samples and 8 negative control DNA samples were also assayed, in addition to 12 replicate/duplicate samples. Similarly, for the HumanMethylation450 BeadChip batch (Batch 3), 6 CEPH samples, 11 positive control samples, 6 negative control samples and 6 replicate samples were assayed. Lastly, twenty duplicate samples were assayed using both Illumina Infinium HumanMethylation27 and HumanMethylation450 BeadChip in order to compare the methylation levels between the two arrays.
Quality control and normalization
Using Illumina GenomeStudio software, DNAm values from the HumanMethylation27 BeadChip assays were scored as beta values, ranging from 0 (unmethylated) to 1 (methylated), resulting in methylation beta values for 27,578 probes. Quality control was done for Batch 1 and Batch 2 combined and then separately for Batch 3. Probes were then excluded if they were on the Y chromosome, positioned at a single nucleotide polymorphism (dbSNP build 137), had high beta values in BSM negative controls (beyond four standard deviations of mean), or were detected in less than 70% of samples. Quality control was also conducted at the sample level, based on the bisulfite conversion ratio and call rate rates (based on a detection p-value of 0.05). Histograms and scatterplots of these statistics were used to determine which samples to exclude (i.e., “outliers”). Similar quality control steps were completed for the samples assayed using the HumanMethylation450 BeadChips, which contained 485,577 CpG site-specific probes.
For the HumanMethylation27 BeadChip arrays, 25,922 (94%) methylation probes passed quality control; for the HumanMethylation450, 441,716 (91%) methylation probes passed quality control. The pairwise correlations for beta values among CEPH replicates were excellent (≥0.97 for Batches 1 and 2, and >0.99 for Batch 3), as were the intra-class correlations of beta values among CEPH replicates (>0.98 for Batches 1 and 2, and >0.99 for Batch 3) and among duplicated study participant samples (>0.93 for Batches 1 and 2, and >0.81 for Batch 3). For 20 samples assessed across batches, the intra-class correlation for beta values of the 24,520 overlapping probes in the HumanMethylation27 and HumanMethylation450 BeadChips was > 0.88. Of samples in Batches 1 and 2, 6 were excluded based on call rates, and one failed bisulfite conversion; in Batch 3, 10 samples were removed following quality control (9 samples failed the bisulfite conversion, one sample with low mean methylation beta value across probes). Following exclusions, we included 69 cases and 87 controls in Batch 1, 146 cases and 176 controls in Batch 2, and 121 cases and 135 controls in Batch 3.
We assessed possible differences by plate and chips within plates (8 BeadChips per plate were assessed with 12 DNAs each) through principal component analyses. Based on the assessment of technical artifacts using principal component analyses, a plate effect was observed within each of the three batches and a chip within batch effect for the HumanMethylation27 data (Additional file 1: Figure S1). Therefore, we adjusted for a plate effect for batch 3 and for chip within plate effect for batches 1 and 2 using a linear model of the logit-transformed beta value for each CpG site, with the unstandardized residuals saved. The logit-transformed locus mean was added back onto the residuals followed by the transformation of the residual to the 0 to 1 scale, producing an “adjusted beta” value for all CpG sites.
Finally, we restricted analyses to probes in common between the DNAm arrays following quality control, excluding 9,341 CpG probes on the Illumina Infinium HumanMethylation27 shown to associate with cell type distribution at q-value < 0.05 [15, 16], as well as 1,363 CpG probes found by Chen et al. to be non-specific (i.e., mapped to multiple places along the genome) . Thus, analyses focused on the remaining 13,816 CpG probes (i.e., 24,520 probes in common between the two panels following quality control minus 9,341 probes associated with cell type distribution minus 1,363 non-specific probes).
Statistical association analysis
We analyzed each batch separately using Van der Waerden rank, or rank-based inverse Gaussian, transformed beta values and combined results across batch using meta-analysis techniques. This allowed us to examine similarity of effects across batches and to estimate the combined effect. Meta-analysis was completed using a random effect meta-analysis. A Woolf’s test of homogeneity of regression coefficients across batches was performed, i.e. the distribution of regression estimates across batches for each probe is compatible with that expected given a common regression estimate. All statistical tests were 2-sided, and analyses of individual batches were carried out using SAS (version 9.3; SAS Institute Inc., Cary, NC) and R (version 2.14.0). Meta-analyses were carried out using the R package rmeta (http://CRAN.R-project.org/package=rmeta). To control for multiple testing, associations with p < 5 × 10-6 were considered statistically significant (e.g., Bonferroni adjustment based on number of independent tests). Pathway analysis used Ingenuity Pathway Analysis (IPA) (Ingenuity® Systems, http://www.ingenuity.com) for genes closest to CpG probes associated with disease status or outcome at p < 0.0001.
The following linear model was used to determine if DNAm levels differ between EOC cases and matched controls for each CpG site. Let, where Y ij represents the adjusted methylation beta value for subject i and CpG probe j (j = 1…, 13816), X i represents disease status for subject i (1 if case and 0 if control), Z i represents covariates for subject i and To identify covariates that differ between EOC cases and controls to include in the model (i.e., potential confounders), potential covariates were examined for association with disease status within a stepwise logistic regression model, resulting in the inclusion of parity/age at first live birth combination (nulliparous, 1-2 and age < = 20 years, 1-2 and age > 20 years, 3+ and age < = 20 years, 3+ and age > =20 years, missing), current alcohol use (never, former, current, missing), current smoking status (never or former, current, missing), enrollment year, and recruitment state (MN vs. non-MN). For each CpG probe j, the disease status parameter () was estimated using the rank-transformed adjusted beta methylation values, along with a 95% confidence interval (CI).
We assessed associations of methylation beta values with overall survival (OS) using Cox proportional hazards regression analyses, adjusted for age at diagnosis, tumor stage (III/IV, I/ II), presence of ascites (yes, no, missing) and volume of residual tumor following debulking surgery (<1 cm, >1 cm, missing) based on stepwise Cox regression analysis. The proportionality assumption was assessed by the analysis of scaled Schoenfeld residuals for all covariates included in the statistical analysis and found to be upheld . We accounted for left truncation using start-stop counting process style of input and estimated hazards ratios (HR) and 95% CIs .
Disease status and DNA methylation
CpG sites associated with disease status (p ≤ 5×10 -7 )
Location of nearest gene (bp)#
Location to Island^
To identify any commonality of highlighted genes within biological pathways, pathway analysis using Ingenuity Pathway Analysis (IPA) was completed for the 155 genes closest to the CpG probes (based on Illumina provided annotation) that were associated with disease status based on a liberal threshold of p < 0.0001. The top pathways enriched for these 155 genes were the telomerase signaling (five genes in our top 155 were in the list of 99 genes within the telomerase signaling pathway; p = 1.24×10-3 for enrichment of pathway) and the paxillin signaling (five genes in our top 155 were in the list of 110 genes within the paxillin signaling pathway; p = 1.42×10-3). The five genes in the telomerase signaling pathway with methylation associated with disease status at p < 0.0001 were HDAC3 (p = 6.33×10-12), IL2RG (p = 4.33×10-6), PIK3C2B (p = 1.97×10-5), PIK3R1 (p = 5.19×10-5), and POT1 (p = 1.38×10-6). PIK3C2B has been implicated in development of glioblastoma multiforme, while mutations in PIK3R1 have been seen in ovarian tumors and cancer cell lines and endometrial cancer [21–23]. POT1 has been found to be associated with tumor stage and telomere length in gastric cancer [24–26]. For the paxillin signaling pathway, the five differentially methylated CpGs were near ARFIP2 (p = 4.60×10-5), ITGB6 (p = 3.95×10-5), PIK3C2B, PIK3R1 and SRC, with some overlap between the top two pathways (PIK3R1 and PIK3C2B).
Survival following EOC and DNA methylation
Similar to the analysis of the disease-associated genes, an exploratory pathway analysis using IPA was completed for the 61 genes closest to the CpG probes most associated with OS (meta-analysis p < 0.01). The top canonical pathways enriched for these 61 genes were relaxin signaling (five genes out of 147; GNA12, GNB1, PIK3R4, RAP1A, TDP2; p = 7.09×10-5 for enrichment of pathway) and CXCR4 signaling (five genes out of 160; GNA12, GNB1, ITPR1, PIK3R4, ROCK1; p = 1.25×10-4) and IL-8 signaling (five genes out of 192; ARRB2, GNA12, GNB1, PIK3R4, ROCK1; p = 3.05×10-4). Three genes (GNB1 (p = 0.006), GNA12 (p = 0.009), and PIK3R4 (p = 0.002)) are part of all three of these canonical pathways.
Via a CpG-by-CpG approach excluding CpGs known to correlate with potentially confounding white blood cell types, we identified methylation CpG sites (and corresponding neighboring genes) with differential specific hyper- or hypo-methylation signals by case-control status and by survival time. To increase power to detect differentially methylated CpG sites, we completed a meta-analysis of results from three DNAm experiments using two genome-wide methylation arrays, restricting focus to high quality probes on both arrays.
A number of CpG sites were found to be differentially methylated between EOC cases and age-matched controls (Table 2). The CpG site that was most differentially methylated between EOC cases and controls was cg04834572 located at the splice region of exon 1 and intron 1 of DUSP13 (p = 1.6×10-14). The blood-based methylation of this CpG site was also reported to be associated with EOC risk in a previous study with p = 0.002 (Figure 2A) . In addition to the replication of the association for the methylation site at DUSP13, four regions identified in this study were replicated for association EOC risk with a p ≤ 0.001, as reported in Teschendorff et al. : cg02449608 (C19orf18, p = 0.0002), cg19280776 (PAG1, p = 8×10-6); cg17271365 (MARCH1, p = 2×10-5), cg10414058 (HDAC3, p = 0.001), and cg24959428 (GBP6, p = 0.001).
Many of the genes neighboring the top associated CpG sites have biological relevance to cancer development. Methylation at a CpG site on chromosome 20 at bp 35973318 (cg0549868), located within the splicing region of exon 1 and intron 1 of gene SRC (35973088-36033821 bp), was found to be associated with EOC risk (p = 4.8×10-7) (Figure 2B). SRC is a proto-oncogene which regulates EGFR, Akt, MAPK1 and NF-κB. SRC is a target for many anticancer drugs . A CpG island (cg14580567, bp 145567271) within HHIP (145567148-145659881 bp) was also found to be associated with EOC risk (p = 5.6×10-11). The genomic region surrounding HHIP (hedgehog-interacting protein) (Figure 2C) has been implicated in many cancers, with hypermethylation of the promoter region found to down-regulate the expression of HHIP found in many tumors, such as gastric and pancreatic cancer . The hedgehog proteins are evolutionarily conserved and are important for a wide range of developmental processes; members of this family control cell proliferation and differentiation, thus linking them with many cancers, including basal-cell carcinoma, small cell lung cancer and pancreatic cancer .
The methylation at a CpG site in the shore of a CpG island, approximately 2.5 kb upstream of HDAC3, was observed to be associated with EOC risk (p =6.3×10-13) (Figure 2D). This association was also observed in a previous study (p = 0.001) . Other studies looking at the role of histone deacetylases (HDACs) found that the expression of HDAC1, along with the expression levels of HDAC2 and HDAC3, to be increased in ovarian tumors compared to levels in benign tumors and normal tissue, suggesting the oncogenic potential of HDACs in ovarian tumors [33–35]. Lastly, a CpG site near CUL7 (cg18159180, p = 2.2×10-8) was differentially methylated between EOC cases and controls; CUL7 has been shown to block Myc-induced apoptosis in a p53-dependent manner (Figure 2E) [36, 37].
In addition to individual CpG sites associated with EOC risk, we also found the telomerase signaling and paxillin signaling pathways to be enriched for genes with CpGs that were differentially methylated between cases and controls. The telomerase signaling pathway and inherited variation in TERT have been found to be associated with the development of EOC and other cancers . The maintenance of functional telomeres is critical in that telomeres that become too short are unable to protect the chromosome from DNA damage. TERT plays an extensive role in the maintenance of functional telomeres, and TERT can be activated by AKT and HSP90 and inhibited by c-Abl. One gene identified to be moderately associated with EOC risk was PIK3R1, which is also a member of the telomerase signaling pathway. PIK3R1 is involved in ATPase and estrogen receptor binding and regulates numerous genes, such as AKT, NFKB, TNF, and is involved in apoptosis, proliferation and differentiation. PIK3R1 has also been linked to epithelial neoplasia and cancer, endometrial ovarian cancer, and endometrioid carcinoma [20, 22, 23, 38]. In contrast, the paxillin signaling pathway is involved in the recruitment of signaling and structural proteins to paxillin required for regulation of cell motility, with many of the paxillin-binding proteins having oncogenic equivalents.
CpG sites associated with overall survival following EOC (p < 10 -3 )
Location of nearest gene (bp)#
Relation to Island^
Single CpG probe analysis of the association of blood-based DNAm with survival following EOC, followed by pathway analyses found the top pathways to all contain three genes (GNB1, GNA12, and PIK3R4), although individual CpG evidence for these three genes were modest. The standard chemotherapy regimen for EOC patients following surgery is a combination therapy involving a taxane (e.g., paclitaxel) and platinum (e.g., cisplatin, carboplatin) agent, increasing our interest in the gene GNA12 found to be associated with response to cisplatin/paclitaxel . Guanine nucleotide binding protein (G protein), beta polypeptide 1 (GNB1) has been recently found to be associated with breast cancer outcomes and clinical and pathological measurements . PIK3R4 is a member of the phosphoinositide 3-kinases (PI3Ks) family that is involved in multiple cell functions (e.g., proliferation, cell survival, degranulation), and this gene is a novel candidate for outcome following EOC.
In summary, we have identified several methylation CpGs sites, using blood-based or leukocyte DNAm, which are differentially methylated by case-control status. Of these CpGs, four CpGs and two genes containing significant CpGs were replicated in an independent study of DNAm and EOC risk. Strengths of our study are large sample size, exclusion of CpGs associated with white blood cell types, and inclusion of relevant covariates. Prior work in a smaller set of cases and controls showed that blood-based DNAm associated with case-control status , thus providing additional evidence to “confirmed” CpG regions associated with EOC. To ensure that none of our findings could be attributed to confounding due to cell type distribution, we removed of probes associating with cell types (which in fact showed very strong associations with case-control status; data not shown). In addition to these strengths, there are also limitations to this study. First, this study was limited to CpG sites assayed on the Illumina array; future application of genome-wide DNA methylation sequencing (i.e., methyl-seq) will enable additional EOC related methylation marks to be discovered. Secondly, the retrospective case-control design used in this study precludes interpretation of these results as indicators of EOC risk. As blood was drawn upon diagnosis, we cannot exclude the possibility that the case-control differences resulted from the cancer itself, from its treatment, or from lifestyle changes. Nonetheless, this short list of CpGs should be of high priority for cohort studies with baseline blood draws and follow-up for later EOC. We note that our survival studies were limited primarily by sample size (336 cases), and thus may have been underpowered to detect modest effects; combining this study with other blood-based methylation case studies will be a key next step.
In conclusion, this early examination of blood-based DNAm provides added experience to a relatively nascent field, suggesting that careful pre-processing and consideration of probes associating with distributions of white blood cell types is critical. We also report specific CpGs that associate either with case-control status or outcome, which are worthy of follow-up in prospective cohort and clinical studies.
Funding was provided by the Fred C. and Katherine B. Andersen Foundation and the US National Institute of Health (P30 CA168524, P50 CA136393, R01 CA122443, P20 GM103418, R21 GM86689).
- Cancer Genome Atlas Research Network: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474 (7353): 609-615. 10.1038/nature10166.View ArticleGoogle Scholar
- Cicek MS, Koestler DC, Fridley BL, Kalli KR, Armasu SM, Larson MC, Wang C, Winham SJ, Vierkant RA, Rider DN, Block MS, Klotzle B, Konecny G, Winterhoff BJ, Hamidi H, Shridhar V, Fan JB, Visscher DW, Olson JE, Hartmann LC, Bibikova M, Chien J, Cunningham JM, Goode EL: Epigenome-wide ovarian cancer analysis identifies a methylation profile differentiating clear-cell histology with epigenetic silencing of the HERG K+ channel. Hum Mol Genet. 2013, 22 (15): 3038-3047. 10.1093/hmg/ddt160.View ArticlePubMedPubMed CentralGoogle Scholar
- Rakyan VK, Down TA, Balding DJ, Beck S: Epigenome-wide association studies for common human diseases. Nat Rev Genet. 2011, 12 (8): 529-541. 10.1038/nrg3000.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsai PC, Spector TD, Bell JT: Using epigenome-wide association scans of DNA methylation in age-related complex human traits. Epigenomics. 2012, 4 (5): 511-526. 10.2217/epi.12.45.View ArticlePubMedGoogle Scholar
- Brennan K, Flanagan JM: Epigenetic epidemiology for cancer risk: harnessing germline epigenetic variation. Methods Mol Biol. 2012, 863: 439-465. 10.1007/978-1-61779-612-8_27.View ArticlePubMedGoogle Scholar
- Brennan K, Flanagan JM: Is there a link between genome-wide hypomethylation in blood and cancer risk?. Cancer Prev Res (Phila). 2012, 5 (12): 1345-1357. 10.1158/1940-6207.CAPR-12-0316.View ArticleGoogle Scholar
- Shenker NS, Polidoro S, van Veldhoven K, Sacerdote C, Ricceri F, Birrell MA, Belvisi MG, Brown R, Vineis P, Flanagan JM: Epigenome-wide association study in the European Prospective Investigation into Cancer and Nutrition (EPIC-Turin) identifies novel genetic loci associated with smoking. Hum Mol Genet. 2013, 22 (5): 843-851. 10.1093/hmg/dds488.View ArticlePubMedGoogle Scholar
- Pedersen KS, Bamlet WR, Oberg AL, de Andrade M, Matsumoto ME, Tang H, Thibodeau SN, Petersen GM, Wang L: Leukocyte DNA methylation signature differentiates pancreatic cancer patients from healthy controls. PLoS One. 2011, 6 (3): e18223-10.1371/journal.pone.0018223.View ArticlePubMedPubMed CentralGoogle Scholar
- Goode EL, Chenevix-Trench G, Song H, Ramus SJ, Notaridou M, Lawrenson K, Widschwendter M, Vierkant RA, Larson MC, Kjaer SK, Birrer MJ, Berchuck A, Schildkraut J, Tomlinson I, Kiemeney LA, Cook LS, Gronwald J, Garcia-Closas M, Gore ME, Campbell I, Whittemore AS, Sutphen R, Phelan C, Anton-Culver H, Pearce CL, Lambrechts D, Rossing MA, Chang-Claude J, Moysich KB, Goodman MT, et al: A genome-wide association study identifies susceptibility loci for ovarian cancer at 2q31 and 8q24. Nat Genet. 2010, 42 (10): 874-879. 10.1038/ng.668.View ArticlePubMedPubMed CentralGoogle Scholar
- Song H, Ramus SJ, Tyrer J, Bolton KL, Gentry-Maharaj A, Wozniak E, Anton-Culver H, Chang-Claude J, Cramer DW, DiCioccio R, Dork T, Goode EL, Goodman MT, Schildkraut JM, Sellers T, Baglietto L, Beckmann MW, Beesley J, Blaakaer J, Carney ME, Chanock S, Chen Z, Cunningham JM, Dicks E, Doherty JA, Durst M, Ekici AB, Fenstermacher D, Fridley BL, Giles G, et al: A genome-wide association study identifies a new ovarian cancer susceptibility locus on 9p22.2. Nat Genet. 2009, 41 (9): 996-1000. 10.1038/ng.424.View ArticlePubMedPubMed CentralGoogle Scholar
- Bojesen SE, Pooley KA, Johnatty SE, Beesley J, Michailidou K, Tyrer JP, Edwards SL, Pickett HA, Shen HC, Smart CE, Hillman KM, Mai PL, Lawrenson K, Stutz MD, Lu Y, Karevan R, Woods N, Johnston RL, French JD, Chen X, Weischer M, Nielsen SF, Maranian MJ, Ghoussaini M, Ahmed S, Baynes C, Bolla MK, Wang Q, Dennis J, McGuffog L, et al: Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer. Nat Genet. 2013, 45 (4): 371-384. 10.1038/ng.2566. 384e371-372View ArticlePubMedPubMed CentralGoogle Scholar
- Shen H, Fridley BL, Song H, Lawrenson K, Cunningham JM, Ramus SJ, Cicek MS, Tyrer J, Stram D, Larson MC, Kobel M, Consortium P, Ziogas A, Zheng W, Yang HP, Wu AH, Wozniak EL, Woo YL, Winterhoff B, Wik E, Whittemore AS, Wentzensen N, Weber RP, Vitonis AF, Vincent D, Vierkant RA, Vergote I, Van Den Berg D, Van Altena AM, Tworoger SS, et al: Epigenetic analysis leads to identification of HNF1B as a subtype-specific susceptibility gene for ovarian cancer. Nat Commun. 2013, 4: 1628.View ArticlePubMedGoogle Scholar
- Painter JN, Anderson CA, Nyholt DR, Macgregor S, Lin J, Lee SH, Lambert A, Zhao ZZ, Roseman F, Guo Q, Gordon SD, Wallace L, Henders AK, Visscher PM, Kraft P, Martin NG, Morris AP, Treloar SA, Kennedy SH, Missmer SA, Montgomery GW, Zondervan KT: Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis. Nat Genet. 2011, 43 (1): 51-54. 10.1038/ng.731.View ArticlePubMedGoogle Scholar
- Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Gayther SA, Apostolidou S, Jones A, Lechner M, Beck S, Jacobs IJ, Widschwendter M: An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS One. 2009, 4 (12): e8274-10.1371/journal.pone.0008274.View ArticlePubMedPubMed CentralGoogle Scholar
- Koestler DC, Marsit CJ, Christensen BC, Accomando W, Langevin SM, Houseman EA, Nelson HH, Karagas MR, Wiencke JK, Kelsey KT: Peripheral blood immune cell methylation profiles are associated with nonhematopoietic cancers. Cancer Epidemiol Biomarkers Prev. 2012, 21 (8): 1293-1302. 10.1158/1055-9965.EPI-12-0361.View ArticlePubMedPubMed CentralGoogle Scholar
- Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, Wiencke JK, Kelsey KT: DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012, 13 (1): 86-10.1186/1471-2105-13-86.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen YA, Choufani S, Ferreira JC, Grafodatskaya D, Butcher DT, Weksberg R: Sequence overlap between autosomal and sex-linked probes on the Illumina HumanMethylation27 microarray. Genomics. 2011, 97 (4): 214-222. 10.1016/j.ygeno.2010.12.004.View ArticlePubMedGoogle Scholar
- Grambsch PM, Therneau TM: Proportional hazards tests and diagnostics based on weighted residuals. Biometrika. 1994, 81 (3): 515-526. 10.1093/biomet/81.3.515.View ArticleGoogle Scholar
- Therneau TM, Grambsch PM, Fleming TR: Martingale-based residuals for survival models. Biometrika. 1990, 77: 147-160. 10.1093/biomet/77.1.147.View ArticleGoogle Scholar
- Sjodahl G, Lauss M, Gudjonsson S, Liedberg F, Hallden C, Chebil G, Mansson W, Hoglund M, Lindgren D: A systematic study of gene mutations in urothelial carcinoma; inactivating mutations in TSC2 and PIK3R1. PLoS One. 2011, 6 (4): e18583-10.1371/journal.pone.0018583.View ArticlePubMedPubMed CentralGoogle Scholar
- Rao SK, Edwards J, Joshi AD, Siu IM, Riggins GJ: A survey of glioblastoma genomic amplifications and deletions. J Neurooncol. 2010, 96 (2): 169-179. 10.1007/s11060-009-9959-4.View ArticlePubMedGoogle Scholar
- Philp AJ, Campbell IG, Leet C, Vincan E, Rockman SP, Whitehead RH, Thomas RJ, Phillips WA: The phosphatidylinositol 3’-kinase p85alpha gene is an oncogene in human ovarian and colon tumors. Cancer Res. 2001, 61 (20): 7426-7429.PubMedGoogle Scholar
- Urick ME, Rudd ML, Godwin AK, Sgroi D, Merino M, Bell DW: PIK3R1 (p85alpha) is somatically mutated at high frequency in primary endometrial cancer. Cancer Res. 2011, 71 (12): 4061-4067. 10.1158/0008-5472.CAN-11-0549.View ArticlePubMedPubMed CentralGoogle Scholar
- Kondo T, Oue N, Yoshida K, Mitani Y, Naka K, Nakayama H, Yasui W: Expression of POT1 is associated with tumor stage and telomere length in gastric carcinoma. Cancer Res. 2004, 64 (2): 523-529. 10.1158/0008-5472.CAN-03-1196.View ArticlePubMedGoogle Scholar
- Ning X, Yang S, Wang R, Zhang R, Guo L, Tie J, Cheng Y, Wang G, Wan S, Fang D: POT1 deficiency alters telomere length and telomere-associated gene expression in human gastric cancer cells. Eur J Cancer Prev. 2010, 19 (5): 345-351. 10.1097/CEJ.0b013e32833b4812.View ArticlePubMedGoogle Scholar
- Wan SM, Tie J, Zhang YF, Guo J, Yang LQ, Wang J, Xia SH, Yang SM, Wang RQ, Fang DC: Silencing of the hPOT1 gene by RNA inference promotes apoptosis and inhibits proliferation and aggressive phenotype of gastric cancer cells, likely through up-regulating PinX1 expression. J Clin Pathol. 2011, 64 (12): 1051-1057. 10.1136/jclinpath-2011-200211.View ArticlePubMedGoogle Scholar
- Mohler H, Fritschy JM, Rudolph U: A new benzodiazepine pharmacology. J Pharmacol Exp Ther. 2002, 300 (1): 2-8. 10.1124/jpet.300.1.2.View ArticlePubMedGoogle Scholar
- Riss J, Cloyd J, Gates J, Collins S: Benzodiazepines in epilepsy: pharmacology and pharmacokinetics. Acta Neurol Scand. 2008, 118 (2): 69-86. 10.1111/j.1600-0404.2008.01004.x.View ArticlePubMedGoogle Scholar
- Salzman C, Miyawaki EK, le Bars P, Kerrihard TN: Neurobiologic basis of anxiety and its treatment. Harv Rev Psychiatry. 1993, 1 (4): 197-206. 10.3109/10673229309017080.View ArticlePubMedGoogle Scholar
- Yeatman TJ: A renaissance for SRC. Nat Rev Cancer. 2004, 4 (6): 470-480. 10.1038/nrc1366.View ArticlePubMedGoogle Scholar
- Katoh Y, Katoh M: Hedgehog signaling pathway and gastric cancer. Cancer Biol Ther. 2005, 4 (10): 1050-1054. 10.4161/cbt.4.10.2184.View ArticlePubMedGoogle Scholar
- Pasca di Magliano M, Hebrok M: Hedgehog signalling in cancer formation and maintenance. Nat Rev Cancer. 2003, 3 (12): 903-911. 10.1038/nrc1229.View ArticlePubMedGoogle Scholar
- Hayashi A, Horiuchi A, Kikuchi N, Hayashi T, Fuseya C, Suzuki A, Konishi I, Shiozawa T: Type-specific roles of histone deacetylase (HDAC) overexpression in ovarian carcinoma: HDAC1 enhances cell proliferation and HDAC3 stimulates cell migration with downregulation of E-cadherin. Int J Cancer. 2010, 127 (6): 1332-1346. 10.1002/ijc.25151.View ArticlePubMedGoogle Scholar
- Khabele D, Son DS, Parl AK, Goldberg GL, Augenlicht LH, Mariadason JM, Rice VM: Drug-induced inactivation or gene silencing of class I histone deacetylases suppresses ovarian cancer cell growth: implications for therapy. Cancer Biol Ther. 2007, 6 (5): 795-801. 10.4161/cbt.6.5.4007.View ArticlePubMedGoogle Scholar
- Jin KL, Pak JH, Park JY, Choi WH, Lee JY, Kim JH, Nam JH: Expression profile of histone deacetylases 1, 2 and 3 in ovarian cancer tissues. J Gynecol Oncol. 2008, 19 (3): 185-190. 10.3802/jgo.2008.19.3.185.View ArticlePubMedPubMed CentralGoogle Scholar
- Jung P, Verdoodt B, Bailey A, Yates JR, Menssen A, Hermeking H: Induction of cullin 7 by DNA damage attenuates p53 function. Proc Natl Acad Sci U S A. 2007, 104 (27): 11388-11393. 10.1073/pnas.0609467104.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim SS, Shago M, Kaustov L, Boutros PC, Clendening JW, Sheng Y, Trentin GA, Barsyte-Lovejoy D, Mao DY, Kay R, Jurisica I, Arrowsmith CH, Penn LZ: CUL7 is a novel antiapoptotic oncogene. Cancer Res. 2007, 67 (20): 9616-9622. 10.1158/0008-5472.CAN-07-0644.View ArticlePubMedGoogle Scholar
- Jaiswal BS, Janakiraman V, Kljavin NM, Chaudhuri S, Stern HM, Wang W, Kan Z, Dbouk HA, Peters BA, Waring P, Dela Vega T, Kenski DM, Bowman KK, Lorenzo M, Li H, Wu J, Modrusan Z, Stinson J, Eby M, Yue P, Kaminker JS, de Sauvage FJ, Backer JM, Seshagiri S: Somatic mutations in p85alpha promote tumorigenesis through class IA PI3K activation. Cancer Cell. 2009, 16 (6): 463-474. 10.1016/j.ccr.2009.10.016.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng B, Gu Y, Xiong Y, Zheng G, He Z: Microarray-assisted pathway analysis identifies MT1X & NFkappaB as mediators of TCRP1-associated resistance to cisplatin in oral squamous cell carcinoma. PLoS One. 2012, 7 (12): e51413-10.1371/journal.pone.0051413.View ArticlePubMedPubMed CentralGoogle Scholar
- Chun JH, Kim HK, Kim E, Kim IH, Kim JH, Chang HJ, Choi IJ, Lim HS, Kim IJ, Kang HC, Park JH, Bae JM, Park JG: Increased expression of metallothionein is associated with irinotecan resistance in gastric cancer. Cancer Res. 2004, 64 (14): 4703-4706. 10.1158/0008-5472.CAN-04-1063.View ArticlePubMedGoogle Scholar
- Desmet CJ, Gallenne T, Prieur A, Reyal F, Visser NL, Wittner BS, Smit MA, Geiger TR, Laoukili J, Iskit S, Rodenko B, Zwart W, Evers B, Horlings H, Ajouaou A, Zevenhoven J, van Vliet M, Ramaswamy S, Wessels LF, Peeper DS: Identification of a pharmacologically tractable Fra-1/ADORA2B axis promoting breast cancer metastasis. Proc Natl Acad Sci U S A. 2013, 110 (13): 5139-5144. 10.1073/pnas.1222085110.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma DF, Kondo T, Nakazawa T, Niu DF, Mochizuki K, Kawasaki T, Yamane T, Katoh R: Hypoxia-inducible adenosine A2B receptor modulates proliferation of colon carcinoma cells. Hum Pathol. 2010, 41 (11): 1550-1557. 10.1016/j.humpath.2010.04.008.View ArticlePubMedGoogle Scholar
- Barreto B, Feher O, Carvalho A, Cristo E, Neves E, Kowalski L, Reis L: Molecular classifiers as predictors of responsiveness to concomitant chemoradiotherapy in SCC of the larynx. AACR Meeting Abstracts 2008. 2008, 2008 (1_Annual_Meeting): 4142.Google Scholar
- Wazir U, Jiang WG, Sharma AK, Mokbel K: Guanine nucleotide binding protein beta 1: a novel transduction protein with a possible role in human breast cancer. Cancer Genomics Proteomics. 2013, 10 (2): 69-73.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/7/21/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.