Abstract
Background: Telomeres protect chromosomal ends, shorten with cellular division, and signal cellular senescence, but unchecked telomere attrition can lead to telomere dysfunction, upregulation of telomerase, and carcinogenesis. Shorter telomeres in peripheral blood leukocytes (PBL) have been associated with elevated cancer risk. Furthermore, genetic variants in and around the TERT gene have been implicated in carcinogenesis.
Methods: We measured relative telomere length (RTL) in PBLs of 911 cases and 948 controls from the New England case–control (NECC) study, a population-based study of ovarian cancer. In addition, we assessed germ line genetic variation in five telomere maintenance genes among 2,112 cases and 2,456 controls from the NECC study and the Nurses' Health Study, a prospective cohort study. ORs and 95% CIs were estimated by logistic regression.
Results: Overall, we observed no differences in telomere length between cases and controls. Compared with women with RTL in the longest tertile, women with RTL in the shortest tertile had no increase in risk (OR = 1.01, 95% CI: 0.80–1.28). However, several SNPs in the TERT gene, including rs2736122, rs4246742, rs4975605, rs10069690, rs2736100, rs2853676, and rs7726159, were significantly associated with ovarian cancer risk. We observed a significant gene-level association between TERT and ovarian cancer risk (P = 0.00008).
Conclusion: Our observations suggest that genetic variation in the TERT gene may influence ovarian cancer risk, but the association between average telomere length in PBLs and ovarian cancer remains unclear.
Impact: The role of telomeres in ovarian carcinogenesis remains unsettled and warrants further investigation. Cancer Epidemiol Biomarkers Prev; 21(3); 504–12. ©2012 AACR.
Background
Telomeres are repeated DNA sequences that protect chromosomal ends, shorten with each cellular division, and signal cellular senescence. However, inactivation of tumor suppressors such as p53 allow cells to circumvent senescence and telomeres erode until a crisis point is reached (1). At this stage, cells are characterized by chromosomal aberrations and the majority undergoes apoptosis. A rare cell may reactivate telomerase expression allowing the cell to escape apoptosis and become immortalized (2). Shorter telomeres in peripheral blood leukocytes (PBL) have been associated with a variety of cancers, including ovarian cancer (3–5).
Telomerase activity has been identified in some normal ovarian and fallopian tube tissue and is thought to contribute in a controlled manner to oogenesis and fertility (6, 7). Conversely, telomerase is not detectable in premalignant lesions but is upregulated in 90% to 97% of ovarian cancers, suggesting deregulation of telomerase is a step in ovarian carcinogenesis (6). Twin studies suggest that the majority (40%–80% heritability) of interindividual variation in telomere length is genetic (8, 9). Telomere length and stability is maintained by several proteins, including telomerase (encoded in part by the TERT gene), TRF1, TRF2, and POT1 that act directly on telomeres, as well as TRF2IP, TINF2, and TNKS that interact directly with components of the telomere complex (10). In addition, the importance of genetic variation in TERT to cancer susceptibility has emerged from recent genome-wide association studies (GWAS) and candidate gene studies (11–17). Here, we report the results of relative telomere length (RTL) measurements in PBLs of 911 ovarian cancer cases and 948 controls and germ line genetic analyses in 5 telomere maintenance genes among 2,112 ovarian cancer cases and 2,456 controls.
Methods
Study populations
New England case–control study.
Data and specimens from the New England case–control (NECC) study of ovarian cancer come from 3 enrollment phases (1992–1997, 1998–2003, and 2003–2008). Briefly, 3,957 women residing in eastern Massachusetts or New Hampshire with a diagnosis of incident ovarian cancer were identified through hospital tumor boards and statewide cancer registries over 3 study phases. Women were not eligible to be a case in the study if they had died, moved, lived outside the study area, had no phone, did not speak English, or had a nonovarian primary tumor after review. Of the 3,083 eligible cases, 2,203 (71% of the eligible cases, 2,076 epithelial cases) agreed to participate. All pathology reports were reviewed by a study pathologist for consistency. Controls were identified through a combination of random digit dialing, drivers' license lists, and town resident lists. In the first phase, 420 (72%) of the eligible women identified through random digit dialing agreed to participate, and 102 (51%) of the eligible women identified through town resident lists agreed to participate. In the second and third phases, 4,366 potential controls were identified, 2,940 were eligible, 1,362 declined to participate by phone or by mail via an “opt-out” postcard, and 1,578 (54%) were enrolled. Potential controls were required to meet the same eligibility criteria as cases. In addition, controls were required to have at least 1 ovary. Controls were frequency matched to cases on age and state of residence. All study participants were interviewed at the time of enrollment about known and suspected ovarian cancer risk factors and asked to complete a food frequency questionnaire. Over 97% of the participants provided a blood specimen. In 558 cases diagnosed at Brigham and Women's Hospital or Massachusetts General Hospital, we abstracted chemotherapy data from medical records. Of the 434 cases with RTL and abstracted clinical data, 302 received chemotherapy. Of these, 51 had blood drawn before her chemotherapy. For the remainder, the mean time between chemotherapy and blood draw was 35 days with an interquartile range of 20 to 165 days.
Nurses' Health Study.
The Nurses' Health Study (NHS) cohort was established in 1976 among 121,700 U.S. female registered nurses, aged 30 to 55 years. Women completed an initial questionnaire, were followed biennially by questionnaire to update exposure status and disease diagnoses, and completed a food frequency questionnaire every 4 years. In 1989 to 1990, 32,826 participants submitted a blood sample (18). Follow-up of the NHS blood study cohort was 96% in 2006. Self-reported diagnoses and deaths reported by family members or identified through the U.S. National Death Index or U.S. Postal Service were followed up for confirmation. A gynecologic pathologist confirmed the diagnosis of all cases through a record review (19). For this analysis, we included all epithelial cases who submitted a blood sample prior to diagnosis or within 4 years after diagnosis (n = 162). All cases were diagnosed before June 1, 2006 and had no history of a prior cancer other than nonmelanoma skin cancer. We randomly selected 3 controls per case from the study participants with DNA available, no prior bilateral oophorectomy, and no history of cancer, other than nonmelanoma skin cancer, at the questionnaire cycle coinciding with diagnosis of the case, matching on age, menopausal status, postmenopausal hormone use, and timing of blood draw.
Telomere length measurements
Genomic DNA from 990 cases and 1,021 controls in the NECC study were standardized and concentrations were measured with PicoGreen quantification by 96-well spectrophotometer (Molecular Devices). We used quantitative PCR as described previously (20, 21) to estimate the ratio (T/S) of telomeric DNA to that of a single gene (36B4) copy number derived from PBLs. Briefly, 5 ng of buffy coat–derived genomic DNA was dried down in a 384-well plate and resuspended in 10 μL of either the telomere or 36B4 PCR reaction mixture for 2 hours at 4°C. The telomere reaction mixture consisted of 1× QIAGEN QuantiTect SYBR Green Master Mix, 2.5 mmol/L of DTT, 270 nmol/L of Tel-1b primer (GGTTTTTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGT), and 900 nmol/L of Tel-2b primer (TCCCGACTATCCCTATCCCTATCCCTATCCCTATCCCTA). The reaction proceeded for 1 cycle at 95°C for 5 minutes, followed by 30 cycles of at 95°C for 15 seconds, and 54°C for 2 minutes. The 36B4 reaction consisted of 1× QIAGEN QuantiTect SYBR Green Master Mix, 300 nmol/L of 36B4U primer (CAGCAAGTGGGAAGGTGTAATCC), and 500 nmol/L of 36B4D primer (CCCATTCTATCATCAACGGGTACAA). The 36B4 reaction proceeded for 1 cycle at 95°C for 5 minutes, followed by 30 cycles at 95°C for 15 seconds, and 58°C for 1 minute 10 seconds. Each 384-well plate contained a 6-point standard curve from 0.6 to 20 ng using pooled buffy coat–derived genomic DNA. By obtaining slopes within −3.33 ± 0.15 for each standard curve, we were able to achieve 94% PCR efficiency. The T/S ratio (-dCt) for each sample was calculated by subtracting the average 36B4 Ct value from the average telomere Ct value. The relative T/S ratio (-ddCt) was determined by subtracting the T/S ratio value of the 5 ng standard curve point from the T/S ratio of each unknown sample. Assays were conducted in triplicate. A total of 926 cases and 961 controls were successfully assayed. Telomere and single gene coefficients of variation (CV) for triplicates were 1.7 and 1.6, respectively. Cases and controls were evenly distributed across plates. Each plate contained its own unique set of 8 quality control samples, generated by aliquotting 4 replicates from 2 designated samples. The intraplate CV for T/S ratio of QC samples was 10%. Our assay was validated with cancer cell lines (SAOS, CYS-SV, and ATS-V1), which showed a strong correlation between T/S ratio and telomere length determined by Southern blot with correlations ranging from 0.82 to 0.99. Cell lines were kindly provided by the laboratory of Dr. Nikhil Munshi of the Dana Farber Cancer Institute. The cell lines were obtained more than 6 months prior to the assay from American Type Culture Collection and were authenticated by the repository with standard cytogenetic, morphologic, antigen expression, and DNA profiling criteria.
SNP selection and genotyping
We identified 53 tagging SNPs in 7 genes involved in telomere maintenance (TERT, POT1, TNKS, TRF1, TINF2, TRF2, and TRF2IP) using publicly available data from the HapMap phase II. SNPs in TINF2 and TRF2IP were excluded because they were monomorphic in whites and POT1 SNP rs6946757 was excluded because we were unable to design an assay for this SNP, resulting in 39 SNPs genotyped on 1,173 cases, 1,200 controls from the first 2 phases of the NECC study. In phase 3 NECC and NHS, we used an alternative platform that only allows multiples of 16 SNPs, therefore, we limited our SNP set to the top 31 SNPs and added a promising candidate SNP from the TERT gene (rs7726159; ref. 11).
DNA was extracted from buffy coat samples with QIAmp (QIAGEN); samples were amplified with Genomiphi (GE Healthcare). All genotyping was carried out at the DF/HCC High Throughput Polymorphism Core. First, we genotyped 39 SNPs on the first 2 phases of the NECC study using 5′ nuclease assays (Taqman) on the Applied Biosystems Prism 7900HT Sequence Detection System (Applied Biosystems). Then, we genotyped 32 SNPs on phase 3 NECC samples and the NHS samples using the Applied Biosystems Taqman OpenArray genotyping platform. Replicate samples (approximately 10%) included for quality control had 100% concordance. Genotyping was carried out by laboratory personnel blinded to case–control status and quality control replicates. Over 95% of the samples were successfully genotyped for each polymorphism.
Statistical analysis
RTL analysis.
RTL values were adjusted for 384-well plate using batch adjustment methods described previously (22). Outliers were identified by the extreme studentized deviate many-outlier procedure (23) and excluded from further analysis (n = 29). Tertiles were determined based on the distribution in controls. Among controls, we tested for differences in the distribution of RTL by categories of epidemiologic characteristics using the Wilcoxon rank-sum and Kruskal–Wallis tests as well as linear regression for multivariate adjusted analyses. Unconditional logistic regression adjusted for the matching factors and study phase were used to calculate ORs and 95% CI for the association between tertiles of RTL and ovarian cancer risk. Multivariate models were additionally adjusted for oral contraceptive use (<3 months, 3 months–2 years, 2–5 years, and >5 years), parity (0, 1, 2, 3, and 4+), tubal ligation, body mass index (BMI; <21.9, 22–24.6, 25.7–28.3, and >28.3), smoking (never, past, and current), family history of breast or ovarian cancer, and race (white and nonwhite).
Histology-specific analyses were conducted comparing RTL in a particular histologic group with all controls. We evaluated the association between RTL and ovarian cancer risk within subgroups of ovulatory cycles, genital talc use, smoking, and BMI. We conducted tests for interaction using a likelihood ratio test to compare models with and without the interaction terms.
Among the subset of cases with telomere length measurements and date of last chemotherapy before blood draw (n = 256), we conducted a sensitivity analysis to assess the influence of recent chemotherapy on the association between RTL in PBLs and ovarian cancer by excluding cases with blood drawn within 3 or 6 months of chemotherapy. Further, we used polytomous logistic regression (24) to assess whether the association between RTL and risk of serous invasive ovarian cancer differed by time between diagnosis and blood draw, as this may also reflect chemotherapy status. We calculated estimates for the association between RTL and ovarian cancer risk for cases with blood draw within 283 days of diagnosis (population median) versus more than 283 after diagnosis.
SNP analysis.
There were no SNPs out of Hardy–Weinberg Equilibrium at P < 0.001 in any data set. SNP associations were evaluated by a log-additive model in which the OR represents the incremental increase or decrease in risk with each additional allele. Unconditional logistic regression adjusted for the matching factors (age and study center) was used for the NECC analyses and conditional logistic regression was used in the NHS analysis. We restricted to whites because we observed that several SNP frequencies varied by race. Multivariate models were additionally adjusted for oral contraceptive use, parity, tubal ligation, and family history of breast or ovarian cancer. Estimates from NECC phases 1 and 2, NECC phase 3, and NHS were combined by meta-analysis using fixed effects models to calculate summary ORs and 95% CIs. We used PACT to account for multiple testing in correlated tests (25). Variables evaluated for effect modification included those associated with ovarian cancer and/or thought to influence telomere length, including estimated lifetime number of ovulatory cycles, talc use, BMI, exercise, dietary folate, and total folate (which includes vitamins and supplements).
To evaluate gene level associations, we employed a principal components (PC) approach described previously by Gauderman and colleagues that accounts for linkage disequilibrium between SNPs (26). Briefly, we estimated the combinations of SNPs, grouped as PCs that represent the genetic variation across the gene. Then, we included the fewest number of PCs that together describe at least 80% of the variation in a logistic regression with ovarian cancer as the outcome. Using a likelihood ratio test, we compared models with and without selected PCs to determine the association between the gene of interest and ovarian cancer risk. All analyses were conducted with SAS (version 9.2; SAS Institute) or Intercooled Stata 9.2 (StataCorp).
Results
Population characteristics for the first 2 phases of the NECC study and the NHS-nested case–control study have been reported previously (27). Participants in NECC phase 3 were similar to those in the previous 2 phases. However, phase 3 participants were slightly older, more likely to have taken oral contraceptives, and cases were more likely to be invasive. Among controls with telomere length measurements (Table 1), RTL was inversely associated with age (r = −0.21, P < 0.0001) and BMI (r = −0.12, P = 0.0002). BMI remained independently, although modestly, associated with RTL after adjustment for age (r = −0.08, P = 0.02). Former smokers had significantly shorter telomeres on average than women who never smoked (age-adjusted P = 0.04), but there were no differences between current and never smokers. Oral contraceptive use, parity, tubal ligation, family history of breast or ovarian cancer, estimated lifetime number of ovulatory cycles, and Jewish ancestry were not associated with telomere length.
Characteristic . | Number of controls . | Batch adjusted telomere length . | Pregressiona . | Correlation . | Pcorrelation . |
---|---|---|---|---|---|
All controls | 947 | 0.77 (0.51) | |||
Age, y | |||||
≤40 | 230 | 0.84 (0.53) | ref. | –0.21 | <0.0001 |
41–50 | 230 | 0.83 (0.56) | 0.62 | ||
51–60 | 234 | 0.75 (0.43) | 0.01 | ||
61–70 | 210 | 0.68 (0.48) | <0.0001 | ||
>70 | 43 | 0.79 (0.63) | 0.19 | ||
Smoking | |||||
Never | 439 | 0.81 (0.54) | ref. | 0.08c | |
Former | 366 | 0.73 (0.42) | 0.04 | ||
Current | 142 | 0.79 (0.61) | 0.16 | ||
Body mass indexb | |||||
Quartile 1 (≤21.9) | 225 | 0.82 (0.55) | ref. | –0.12 | 0.0002 |
Quartile 2 (22–24.6) | 251 | 0.83 (0.62) | 0.69 | –0.08d | |
Quartile 3 (24.7–28.3) | 230 | 0.73 (0.34) | 0.31 | ||
Quartile 4 (≥28.4) | 241 | 0.71 (0.47) | 0.06 |
Characteristic . | Number of controls . | Batch adjusted telomere length . | Pregressiona . | Correlation . | Pcorrelation . |
---|---|---|---|---|---|
All controls | 947 | 0.77 (0.51) | |||
Age, y | |||||
≤40 | 230 | 0.84 (0.53) | ref. | –0.21 | <0.0001 |
41–50 | 230 | 0.83 (0.56) | 0.62 | ||
51–60 | 234 | 0.75 (0.43) | 0.01 | ||
61–70 | 210 | 0.68 (0.48) | <0.0001 | ||
>70 | 43 | 0.79 (0.63) | 0.19 | ||
Smoking | |||||
Never | 439 | 0.81 (0.54) | ref. | 0.08c | |
Former | 366 | 0.73 (0.42) | 0.04 | ||
Current | 142 | 0.79 (0.61) | 0.16 | ||
Body mass indexb | |||||
Quartile 1 (≤21.9) | 225 | 0.82 (0.55) | ref. | –0.12 | 0.0002 |
Quartile 2 (22–24.6) | 251 | 0.83 (0.62) | 0.69 | –0.08d | |
Quartile 3 (24.7–28.3) | 230 | 0.73 (0.34) | 0.31 | ||
Quartile 4 (≥28.4) | 241 | 0.71 (0.47) | 0.06 |
aBased on age-adjusted linear regression comparing batch-adjusted RTL of each category to reference group.
bBMI missing for 3 controls.
cBased on the Kruskal–Wallis test.
dCorrelation between BMI and RTL adjusted for age.
Overall, we observed no significant difference in mean RTL between 911 women with epithelial ovarian cancer (mean = 0.75, SD = 0.40) versus 948 controls (mean = 0.77, SD = 0.51). Compared with women with telomeres in the longest tertile of length, those with telomeres in the shortest tertile of length had no increased risk of ovarian cancer overall (OR = 1.01, 95% CI: 0.80–1.28; Table 2) or serous invasive cancer (OR = 1.11, 95% CI: 0.82–1.51). Associations did not significantly differ between high-grade (OR = 1.09, 95% CI: 0.80–1.50) serous invasive ovarian cancer and low-grade (OR = 0.88, 95% CI: 0.30–2.64) serous invasive ovarian cancer. Results were similar for other histologic subtypes and when considering quartile or quintile cutpoints. Evaluating the RTL-ovarian cancer association within subgroups did not reveal any interactions with ovulatory cycles, genital talc use, endometriosis, or smoking. However, the association between RTL and ovarian cancer differed significantly by quartiles of BMI (Pheterogeneity = 0.01). Among overweight women (BMI > 28.3), women with the shortest telomeres had a borderline decreased risk compared with women with the longest telomeres (OR = 0.62, 95% CI: 0.38–1.01), but among women with a BMI less than 21.9 (OR = 1.29, 95% CI: 0.79–2.10) and women with a BMI between 22 and 24.6 (OR = 1.06, 95% CI: 0.67–1.66) shorter telomeres were not associated with ovarian cancer risk.
. | n . | RTL mean (SD) . | Tertile of RTL . | ||
---|---|---|---|---|---|
. | . | . | T1 (shortest) OR (95% CI)a . | T2 OR (95% CI)a . | T3 (longest) . |
Controls, n (%) | 947 | 0.77 (0.51) | 314 (33) | 305 (32) | 328 (35) |
Cases, n (%) | 911 | 0.75 (0.40) | 319 (35) | 275 (30) | 317 (35) |
Overall | 1.01 (0.80–1.28) | 0.95 (0.75–1.20) | 1.00 | ||
Serous | 482 | 0.73 (0.37) | 1.05 (0.79–1.39) | 0.94 (0.71–1.25) | 1.00 |
Mucinous | 120 | 0.79 (0.39) | 0.77 (0.46–1.32) | 1.41 (0.89–2.24) | 1.00 |
Endometrioid | 132 | 0.77 (0.46) | 0.88 (0.55–1.42) | 0.89 (0.55–1.43) | 1.00 |
Clear cell | 121 | 0.72 (0.41) | 1.07 (0.66–1.74) | 0.84 (0.51–1.39) | 1.00 |
. | n . | RTL mean (SD) . | Tertile of RTL . | ||
---|---|---|---|---|---|
. | . | . | T1 (shortest) OR (95% CI)a . | T2 OR (95% CI)a . | T3 (longest) . |
Controls, n (%) | 947 | 0.77 (0.51) | 314 (33) | 305 (32) | 328 (35) |
Cases, n (%) | 911 | 0.75 (0.40) | 319 (35) | 275 (30) | 317 (35) |
Overall | 1.01 (0.80–1.28) | 0.95 (0.75–1.20) | 1.00 | ||
Serous | 482 | 0.73 (0.37) | 1.05 (0.79–1.39) | 0.94 (0.71–1.25) | 1.00 |
Mucinous | 120 | 0.79 (0.39) | 0.77 (0.46–1.32) | 1.41 (0.89–2.24) | 1.00 |
Endometrioid | 132 | 0.77 (0.46) | 0.88 (0.55–1.42) | 0.89 (0.55–1.43) | 1.00 |
Clear cell | 121 | 0.72 (0.41) | 1.07 (0.66–1.74) | 0.84 (0.51–1.39) | 1.00 |
aORs and 95% CIs are age adjusted, study center, study, oral contraceptive use (<3 mo, 3 mo–2 y, 2–5 y, and >5 y), parity (0, 1, 2, 3, and 4+), tubal ligation, BMI (<21.9, 22–24.6, 25.7–28.3, and >28.3), smoking (never, past, and current), family history of breast or ovarian cancer, and race.
The association between RTL and risk of serous invasive ovarian cancer did not differ significantly between cases with blood drawn less than or equal to median time between diagnosis and blood draw (OR = 0.89, 95% CI: 0.67–1.18) compared with those with blood drawn longer after diagnosis (OR = 1.21, 95% CI: 0.90–1.62, Pheterogeneity = 0.13). Although nonsignificant, the association between RTL and ovarian cancer risk tended to increase as cases that had recently had chemotherapy were excluded from the analyses. For example, ORs for ovarian cancer, comparing shortest with longest tertiles of RTL, were 1.20 (95% CI: 0.66–2.19) when restricting to cases with blood drawn more than 3 months after chemotherapy (n = 82) and 1.51 (95% CI: 0.67–3.40) when restricting to cases with blood drawn more than 6 months after chemotherapy (n = 53).
We genotyped 40 tagging SNPs in 5 telomere maintenance genes (TERT, TRF2, TRF1, TNKS, and POT1) in the NECC and NHS studies. Minor allele frequencies were similar across all 3 data sets. Overall, we observed 7 significant associations for TERT SNPs, but no significant associations for other SNPs. TERT SNPs rs4246742, rs4975605, rs10069690, rs2736100, rs2853676, and rs7726159 were associated with 11% to 19% increase in ovarian cancer risk with each additional variant allele, whereas rs2736122 was associated with an 11% decrease in risk with each allele (Table 3), but none of these associations remained significant after permutation testing. On a gene level, estimated PCs accounting for 88% of the genetic variation in TERT were significantly associated with ovarian cancer (P = 0.00008), but PCs of other genes were not associated with risk.
. | . | Per allele . | . |
---|---|---|---|
Gene . | rs number . | ORb (95% CI) . | Pgened . |
TERT | rs2736122 | 0.89 (0.81–0.99) | 0.00008 |
TERT | rs2075786 | 1.06 (0.97–1.17) | |
TERT | rs4246742 | 1.19 (1.05–1.35) | |
TERT | rs6882077c | 1.31 (0.27–6.29) | |
TERT | rs4975605 | 1.11 (1.02–1.22) | |
TERT | rs10069690 | 1.11 (1.00–1.23) | |
TERT | rs2242652c | 1.06 (0.91–1.24) | |
TERT | rs2736100 | 1.12 (1.02–1.23) | |
TERT | rs2853676 | 1.13 (1.02–1.25) | |
TERT | rs2736098 | 0.99 (0.89–1.09) | |
TERT | rs7726159c | 1.17 (1.02–1.34) | |
TRF2 | rs251796 | 0.94 (0.85–1.04) | 0.59 |
TRF2 | rs153045 | 0.94 (0.86–1.04) | |
TRF2 | rs3785074 | 1.06 (0.96–1.18) | |
TRF2 | rs166134c | 0.79 (0.48–1.28) | |
TRF2 | rs8061382c | 1.07 (0.74–1.55) | |
TRF1 | rs2975842 | 1.01 (0.93–1.11) | 0.99 |
TRF1 | rs2975852 | 1.03 (0.94–1.14) | |
TRF1 | rs6989159c | 0.95 (0.40–2.25) | |
TRF1 | rs6989493c | 1.19 (0.44–3.26) | |
TRF1 | rs12334686 | 1.02 (0.93–1.12) | |
TRF1 | rs6982126 | 0.97 (0.87–1.08) | |
TRF1 | rs2981096c | 0.94 (0.69–1.26) | |
TRF1 | rs10107605 | 1.00 (0.87–1.15) | |
TRF1 | rs1545827 | 1.01 (0.93–1.10) | |
TNKS | rs1539041 | 0.98 (0.89–1.08) | 0.34 |
TNKS | rs3802650 | 1.02 (0.94–1.12) | |
TNKS | rs10509637 | 1.01 (0.89–1.15) | |
TNKS | rs1772180 | 0.93 (0.85–1.02) | |
TNKS | rs10509639 | 1.11 (0.94–1.30) | |
TNKS | rs1772186 | 0.98 (0.86–1.11) | |
TNKS | rs10881982c | 1.03 (0.72–1.49) | |
TNKS | rs12412538 | 1.07 (0.97–1.18) | |
TNKS | rs7087365 | 1.06 (0.96–1.17) | |
POT1 | rs929365 | 0.96 (0.79–1.18) | 0.79 |
POT1 | rs7801661 | 1.02 (0.92–1.13) | |
POT1 | rs11972248 | 0.96 (0.86–1.07) | |
POT1 | rs12532038 | 1.01 (0.91–1.11) | |
POT1 | rs4360236 | 1.14 (0.99–1.32) | |
POT1 | rs2896361 | 0.98 (0.90–1.08) |
. | . | Per allele . | . |
---|---|---|---|
Gene . | rs number . | ORb (95% CI) . | Pgened . |
TERT | rs2736122 | 0.89 (0.81–0.99) | 0.00008 |
TERT | rs2075786 | 1.06 (0.97–1.17) | |
TERT | rs4246742 | 1.19 (1.05–1.35) | |
TERT | rs6882077c | 1.31 (0.27–6.29) | |
TERT | rs4975605 | 1.11 (1.02–1.22) | |
TERT | rs10069690 | 1.11 (1.00–1.23) | |
TERT | rs2242652c | 1.06 (0.91–1.24) | |
TERT | rs2736100 | 1.12 (1.02–1.23) | |
TERT | rs2853676 | 1.13 (1.02–1.25) | |
TERT | rs2736098 | 0.99 (0.89–1.09) | |
TERT | rs7726159c | 1.17 (1.02–1.34) | |
TRF2 | rs251796 | 0.94 (0.85–1.04) | 0.59 |
TRF2 | rs153045 | 0.94 (0.86–1.04) | |
TRF2 | rs3785074 | 1.06 (0.96–1.18) | |
TRF2 | rs166134c | 0.79 (0.48–1.28) | |
TRF2 | rs8061382c | 1.07 (0.74–1.55) | |
TRF1 | rs2975842 | 1.01 (0.93–1.11) | 0.99 |
TRF1 | rs2975852 | 1.03 (0.94–1.14) | |
TRF1 | rs6989159c | 0.95 (0.40–2.25) | |
TRF1 | rs6989493c | 1.19 (0.44–3.26) | |
TRF1 | rs12334686 | 1.02 (0.93–1.12) | |
TRF1 | rs6982126 | 0.97 (0.87–1.08) | |
TRF1 | rs2981096c | 0.94 (0.69–1.26) | |
TRF1 | rs10107605 | 1.00 (0.87–1.15) | |
TRF1 | rs1545827 | 1.01 (0.93–1.10) | |
TNKS | rs1539041 | 0.98 (0.89–1.08) | 0.34 |
TNKS | rs3802650 | 1.02 (0.94–1.12) | |
TNKS | rs10509637 | 1.01 (0.89–1.15) | |
TNKS | rs1772180 | 0.93 (0.85–1.02) | |
TNKS | rs10509639 | 1.11 (0.94–1.30) | |
TNKS | rs1772186 | 0.98 (0.86–1.11) | |
TNKS | rs10881982c | 1.03 (0.72–1.49) | |
TNKS | rs12412538 | 1.07 (0.97–1.18) | |
TNKS | rs7087365 | 1.06 (0.96–1.17) | |
POT1 | rs929365 | 0.96 (0.79–1.18) | 0.79 |
POT1 | rs7801661 | 1.02 (0.92–1.13) | |
POT1 | rs11972248 | 0.96 (0.86–1.07) | |
POT1 | rs12532038 | 1.01 (0.91–1.11) | |
POT1 | rs4360236 | 1.14 (0.99–1.32) | |
POT1 | rs2896361 | 0.98 (0.90–1.08) |
aRestricted to white women.
bAdjusted for age (continuous), center, oral contraceptive use (never, 1–24, 25–72, and >72 months), parity (0, 1, and 2+), family history of breast or ovarian cancer, tubal ligation, and genotyping batch. Meta-analysis of NECC phases 1 and 2, NECC phase 3, and NHS using fixed effects model. We observed significant heterogeneity between data sets for TRF1 SNP rs12334686 (Pheterogeneity = 0.04) and TNKS SNP rs1772180 (Pheterogeneity = 0.02).
crs6882077, rs2242652, rs166134, rs8061382, rs6989159, rs6989493, rs2981096, and rs10881982 were not genotyped in NECC phase 3 and NHS data sets and rs7726159 was not genotyped in NECC phases 1 and 2.
dPC analyses were used to determine gene-level associations with ovarian cancer risk accounting for linkage disequilibrium between SNPs.
Interestingly, histology-specific associations revealed that TERT SNPs rs4246742 (pooled OR = 1.43, 95% CI: 1.12–1.83), rs4975605 (pooled OR = 1.26, 95% CI: 1.05–1.51), and rs10069690 (pooled OR = 1.28, 95% CI: 1.04–1.57) were most strongly associated with endometrioid tumors, whereas rs2853676 was more strongly associated with serous invasive ovarian cancer (OR = 1.14, 95% CI: 1.00–1.30). rs2736100 was associated with a significant increase in borderline tumors (pooled OR = 1.24, 95% CI: 1.05–1.46). In stratified analyses, we evaluated the association between telomere maintenance SNPs and ovarian cancer risk by estimated lifetime number of ovulatory cycles, BMI, exercise, genital talc use, folate intake, and smoking; however, we did not observe any significant interactions after adjustment for multiple testing.
Among 918 white controls from the NECC study, we evaluated the association between SNPs in telomere maintenance genes and RTL in PBLs. We observed significant age-adjusted differences in telomere length by genotype for POT1 SNPs rs929365 (P = 0.02) and rs4360236 (P = 0.04), but these SNPs were not associated with ovarian cancer risk.
Discussion
Our data show an association between genetic variation in the TERT gene and ovarian cancer risk in 3 data sets with significant associations between 6 individual TERT SNPs and ovarian cancer risk and a highly significant gene level association with risk. In our telomere length analysis, we observed shorter telomeres in PBLs with increasing age, former smoking, and higher BMI as previously reported in the literature (28, 29) but no significant association between RTL in PBLs and ovarian cancer risk in the largest study of telomere length and ovarian cancer risk to date.
When considered independently, the magnitude of the association between any individual TERT SNP we genotyped and ovarian cancer risk was small and of borderline significance, and a combined analysis of 4 ovarian cancer GWAS studies showed no significant associations between these individual TERT SNPs and overall ovarian cancer risk (unpublished data). However, when we considered the TERT variants together in a PC analysis that accounts for linkage disequilibrium between the SNPs, we observed a highly significant association (P = 0.00008), suggesting that these individual SNPs do not capture the importance of genetic variation in this region. In the context of other recent reports of significant associations between TERT SNPs and risk of various cancers, including ovarian cancer (11–17), our data suggest that genetic variation in the TERT gene may be important to ovarian cancer susceptibility. In a study of 1,309 SNPs in 173 genes involved in stromal epithelial crosstalk, Johnatty and colleagues observed that TERT rs7726159 was significantly associated with serous invasive ovarian cancer (ORper allele = 1.14, 95% CI: 1.04–1.24) and an inverse association of borderline significance between rs2736100 and serous ovarian cancer (OR = 0.88, 95% CI: 0.77–1.01; ref. 11). Furthermore, rs2736100 has emerged from recent GWAS studies as a potential susceptibility allele for glioma (OR = 1.27, P = 2 × 10−17; ref. 15), lung cancer (OR = 1.14, P = 4 × 10−6; ref. 12), and testicular cancer (OR = 1.33, P = 8 × 10−15; ref. 17). Based on this accumulating evidence, the area around the TERT gene has been hypothesized to be cancer polymorphism “hot spot” that may influence many types of cancer (11, 13, 14).
Given the emerging evidence that p53 mutations are likely precursors to most serous ovarian cancers (36) and p53 mutations are also observed in endometrioid ovarian tumors (37), the telomere maintenance pathway may have particular relevance to these tumors. Interestingly, Kuhn and colleagues found that telomeres measured in tissue from serous tubal intraepithelial carcinoma (STIC) were significantly shorter than those found in concurrent high-grade serous ovarian carcinoma and suggested that the shortened telomeres STICs represent an early carcinogenic event whereas longer telomeres in coexisting high-grade serous tumors are a result of subsequent activation of telomerase and immortalization (38).
Although none of the TERT SNPs we found to be associated with ovarian cancer was associated with telomere length in PBLs of controls, the influence of these SNPs on other telomere attributes cannot be ruled out. Telomere attrition can occur due to a range of insults that increase cellular turnover, accelerate telomeric DNA damage, or result in defective telomere repair. Therefore, the influence of a single SNP on telomere length may be difficult to detect. Furthermore, these SNPs may influence carcinogenesis through a pathway that does not involve telomere length. For instance, Artandi and colleagues have shown that TERT interacts directly with WNT transcription factors and cooperate in activating WNT responsive promoters, whereas Mancke and colleagues have recently shown telomerase interacts with β-catenin (reviewed in ref. 39). Interestingly, WNT and β-catenin have been implicated in ovarian carcinogenesis (40).
Our telomere length and ovarian cancer risk results contrast those of 2 small case–control studies that previously reported significantly increased ovarian cancer risk with shorter telomeres (4, 5). Idei and colleagues observed shorter mean telomere lengths in DNA from plasma from 32 cases (≥stage 3) versus controls in a Japanese population. In addition, long-term (≥10 years) ovarian cancer survivors had significantly longer telomeres than other cases, suggesting a survival benefit with increased telomere length (4). Similarly, Mirabello and colleagues reported shorter leukocyte telomere length in 99 serous ovarian cancer cases compared with 100 controls in a Polish population, with a 3-fold increase in risk for women with telomere length in the shortest tertile (OR = 3.39, 95% CI: 1.54–7.46; ref. 5). Differences in study populations, study design, sample size, timing of chemotherapy, and type of specimen or assay may explain inconsistencies between our own results and previous studies. Telomere length has been shown to vary by ethnicity (30), which varies between these studies. Furthermore, telomere length in DNA obtained from plasma (4) may differ from that derived from buffy coat, and cases in the Polish study were recruited sooner after diagnosis (median, 2.4 months) than in our study (median, 9.3 months), such that time since chemotherapy likely differed between the studies.
Several studies have suggested that chemotherapy may influence measurable telomere length in PBLs, but conclusions are mixed showing both increases and decreases in length as well as changes in telomerase activity after chemotherapy (4, 31–34). In prospectively collected serial samples, one study reported that telomere length postchemotherapy was approximately 550 base pairs shorter than prechemotherapy (34), whereas another reported significant telomere shortening in peripheral blood mononuclear cells after multiple cycles of chemotherapy (33). Cisplatin, a standard therapy for ovarian cancer, forms DNA adducts with telomeric sequences and results in telomere loss in cell lines (31). In our study, telomere lengths were shorter on average as time since chemotherapy increased, suggesting an influence of chemotherapy. Thus, the association between RTL and ovarian cancer risk may be obscured in samples collected after treatment. Yet, Mirabello and colleagues observed no difference in telomere length between samples collected before (telomere length = 0.89, 95% CI: 0.89–0.99) and after chemotherapy (telomere length = 0.91, 95% CI: 0.82–0.99). Regardless, telomere length measurements in postdiagnostic samples can be difficult to interpret because samples collected soon after diagnosis are likely to be influenced by chemotherapy or active disease and those collected later are likely influenced by survival (35).
Given the well documented inverse correlation between BMI and RTL noted in our study and by others (28, 29), our subgroup analyses showing a reduced risk of ovarian cancer for women in the highest quartile of weight is counterintuitive as one would expect the highest cancer risk for obese women with short telomeres. These results are likely a chance finding due to a limited distribution of telomere length in this subgroup.
A limitation of our study is the use of telomere length in PBLs as a surrogate for telomere length in the ovaries or fallopian tubes. A range of exposures including chronic infection or other immune insults could accelerate cellular turnover in white blood cells leading to shorter telomeres compared with tissue. However, telomere length in a variety of tissues are highly synchronous (41), and matched blood and tissue samples correlated (4, 42). Although this is a large study with over 2,112 cases and 2,456 controls in the SNP analysis and 911 cases in 948 controls in the telomere length analysis, power is limited to detect significant associations within subgroups of histologic types, particularly rarer subtypes.
Our findings regarding genetic variation in the TERT gene and ovarian cancer risk add to the growing body of literature suggesting the importance of this genetic region to cancer development. Further research is needed in this area to understand how changes in telomere length over time may influence ovarian carcinogenesis in a prospective setting and interaction with a p53 pathway of development.
Disclosure of Potential Conflicts of Interest
No conflicts of interest were disclosed.
Acknowledgments
The authors thank the ovarian cancer GWAS investigators for access to the combined GWAS analysis results. In addition, the authors thank the participants and staff of the NHS, for their valuable contributions as well as the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, and WY.
Grant Support
This work was supported by NIH grants R01CA54419, P50CA105009, P01CA876969, and R01CA49449; Department of Defense grant W81XWH-10-1-0280; and the Ovarian Cancer Research Fund through the Liz Tilberis Scholarship.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.