Genetic variants within the hTERT gene and the risk of colorectal cancer in Lynch syndrome

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

ABSTRACT

Lynch syndrome is an inherited cancer-predisposing disorder caused by germline mutations in the DNA mismatch repair (MMR) genes but there is a high degree of variability in cancer risk observed among carriers, suggesting the existence of modifying factors. Our aim was to investigate variants within the hTERT gene as a potential colorectal cancer (CRC) risk modifier for MMR gene mutation carriers. We identified 1098 MMR gene mutation carriers (420 MLH1, 481 MSH2, 126 MSH6, 53 PMS2 and 18 EPCAM) from 330 families recruited from either family cancer clinics or population cancer registries of the Australasian Colorectal Cancer Family Registry between 1997 and 2012. Using weighted Cox regression after adjusting for ascertainment bias, we estimated associations between 23 SNPs within the hTERT gene and CRC risk. During 46,836 person-years observation, 392 (36%) carriers were diagnosed with CRC at a mean age of 42.2 (standard deviation 11.4) years. There was no evidence of association between any of the hTERT SNPs and CRC risk, overall and stratified by sex and MMR gene mutated, after adjustment for multiple testing. Our findings suggest no evidence for clinical utility of the SNPs within the hTERT gene in Lynch syndrome.

INTRODUCTION

Lynch syndrome, formerly known as hereditary non-polyposis colorectal cancer, is an autosomal dominant disorder caused by germline mutations in one of the DNA mismatch repair (MMR) genes MLH1, MSH2, MSH6, and PMS2, or EPCAM. MMR gene mutation carriers are at increased risk of developing cancers of the colorectum and endometrium, as well as cancers of the ovary, kidney, pancreas, stomach, and urinary bladder. The risk of colorectal cancer (CRC) to age 70 years for MMR gene mutation carriers is reported to be between 12 and 50% [1, 2]. In addition to environmental factors such as obesity and smoking [3], there is evidence suggesting the existence of genetic factors that may contribute to the variability in cancer risks [4, 5].

The hTERT gene (MIM 187270) encodes telomerase reverse transcriptase, the catalytic subunit of telomerase, an important protein for maintaining telomere length. Genetic variation within this gene has been shown to affect telomerase activity and telomere length and is thought to underlie an increased risk for cancer [6]. Genetic association studies in breast, ovarian, endometrial, lung, glioma, and pancreatic cancers have identified multiple single nucleotide polymorphisms (SNPs) in the hTERT gene that are associated with an increased risk of cancer [7]. Several studies have also identified an association between CRC risk and SNPs within hTERT at 5p15.33 [8, 9]. We hypothesized that genetic variation within hTERT may act as a genetic modifier of CRC risk for MMR gene mutation carriers.

RESULTS

A total of 1098 carriers of a germline mutation in a MMR gene (18 in EPCAM, 420 in MLH1, 481 in MSH2, 126 in MSH6, and 53 in PMS2) from 330 independent families was included in this study. During 46,757 person-years observation, 392 (36%) carriers were diagnosed with CRC at a mean age of 42.2 (standard deviation 11.4) years (Table 1). Of these, 90% of cancer diagnoses were confirmed using pathology reports, medical records, cancer registry reports, and death certificates.

Of the 23 hTERT SNPs that we measured, five SNPs (rs13361701, rs2853691, rs4075202, rs2736100, rs4975612) deviated from Hardy–Weinberg equilibrium in CRC-unaffected carriers (Supplementary Table 1). After adjusting for multiple testing and age and sex of the carriers, there was no evidence of an association between CRC risk and any of these 23 SNPs per allele (Figure 1) or for homozygous or heterozygous carriers of minor allele (Supplementary Table 2). When we stratified by sex of the carriers and the MMR gene mutated, there was also no evidence of any association (detail results not shown). When we censored carriers by age 45 years, we found no evidence of any association (detail results not shown).

We combined the 23 SNPs to examine the effect of haplotypes on CRC risk. Of the 88 haplotypes with frequency greater than 0.002, we estimated the associations between five most common haplotypes and CRC risk. Again, there was no evidence for any association (detail results not shown).

DISCUSSION

In this study, we examined whether SNPs within the hTERT gene act as genetic modifiers of CRC risk for MMR gene mutation carriers. We found no evidence of association between any of the 23 SNPs within the hTERT gene and CRC risk either overall or when stratified by sex or specific MMR gene mutation, suggesting that hTERT SNPs do not modify the risk of CRC for MMR gene mutation carriers.

A recent study of 930 MMR gene mutation carriers reported no overall difference in genotype frequencies between carriers affected and unaffected with CRC for the rs2075786 SNP in hTERT (RR=2.46; 95% CI=0.78–7.82) [10], a finding that is consistent with ours. This study also reported marginally significant evidence of association between the AA genotype of the rs2075786 SNP and cancer risk for carriers diagnosed with Lynch syndrome related-cancer diagnosed before the age of 45 years (RR=2.90; 95% CI=1.02–8.26; p=0.05) but we found no evidence of an association for that SNP and CRC risk before age 45 years.

To date, the search for genetic modifiers of CRC risk for MMR gene mutation carriers has provided inconsistent results. We have shown previously that the 11 SNPs identified from genome-wide association studies (GWAS) of CRC for the general population are not associated with CRC risk for MMR gene mutation carriers [11]. Three studies observed two variants, 8q23.3 (rs16892766) and 11q23.1 (rs3802842), to be associated with increased risk of CRC in MMR gene mutation carriers especially for females only [12, 13] or MLH1 mutation carriers only [5, 13]. In contrast, another study of 748 mutation carriers did not report this association[14].

Additional SNPs have been identified from CRC GWAS, including the rs2853668 SNP within the TERT-CLPTM1L locus on 5p15.33 [9], but it appears that some of these have no relevance to the risk of CRC for MMR gene mutation carriers. An association between the rs2736100 SNP in hTERT and CRC has also been reported in sporadic setting [8] but we found no evidence for an association between this SNP and CRC risk for MMR gene mutation carriers.

Our study has the largest number of MMR gene mutation carriers assessed to date for genetic modifiers of the hTERT gene on CRC risk. A possible limitation of our study is its generalizability only to MMR gene mutation carriers with substantial survival, as to be included in the analysis cases had to survive long enough to provide a blood sample for DNA testing. There were 5 SNPs deviated from Hardy–Weinberg equilibrium in CRC-unaffected carriers in our sample. This might be explained by small sample size and non-random sampling of carriers.

In conclusion, our findings suggest that common genetic variants within the hTERT gene do not contribute to CRC risk for carriers of MMR gene mutation carriers. More information is important clinically, to provide more accurate risk assessment and risk management. Subsequent studies should focus on testing large cohorts of MMR gene mutation carriers using agnostic genome-wide approaches in order to identify genetic modifiers of risk.

MATERIALS AND METHODS

Study recruitment and data collection

The study sample comprised confirmed heterozygote carriers of pathogenic mutations in MMR genes who were participants in the Australasian Colorectal Cancer Family Registry [15]. Between 1997 and 2012, the Australasian Colorectal Cancer Family Registry recruited families via: population-based probands who were recently diagnosed colorectal cancer cases from the Victorian Cancer Registry; or clinic-based probands who were enrolled from multiple-case families referred to family cancer clinics in Melbourne, Adelaide, Perth, Brisbane, Sydney, and Auckland. Written informed consent was obtained from all participants, and the study protocol was approved by the institutional human research ethics committee at each center.

Information on demographics, personal characteristics, personal and family history of cancer, cancer-screening history, and history of polyps, polypectomy, and other surgeries was obtained by questionnaires from all probands and participating relatives. Participants were followed up approximately every five years after baseline to update this information. The present study was based on all available baseline and follow-up data.

MMR gene mutation testing

Testing for MMR germline mutations was performed for all population-based probands with CRC displaying evidence of impaired MMR function as evidenced by either tumor microsatellite instability (MSI) and/or a lack of MMR protein expression by immunohistochemistry. The youngest-onset colorectal case participants from each clinic-based family were also genotyped, regardless of tumor MSI or MMR protein expression status. Mutation testing for the MLH1, MSH2 and MSH6 genes was performed by Sanger sequencing or denaturing high performance liquid chromatography, followed by confirmatory DNA sequencing. Large duplication and deletion mutations including those involving EPCAM, which lead to MSH2 methylation, were detected by Multiplex Ligation Dependent Probe Amplification (MLPA). PMS2 mutation testing involved a modified protocol from Senter et al.[16] where exons 1–5, 9 and 11–15 were amplified in three long range PCRs followed by nested exon specific PCR/sequencing while the remaining exons (6, 7, 8 and 10) were amplified and sequenced direct from genomic DNA. Large-scale deletions in PMS2 were detected using the P008-A1 MLPA kit (MRC Holland). The relatives of probands with a pathogenic MMR germline mutation, who provided a blood sample, underwent testing for the specific mutation identified in the proband.

Genotyping of the SNPs

Using data from the International HapMap project (www.HapMap.org) and Haploview program (version 3.12) and a minimum r² threshold of 0.8 we identified a set of tagging SNPs to capture the genetic variation in the hTERT gene. SNPs previously associated with CRC risk [8, 9, 17] or from the single study of CRC risk in MMR gene mutation carriers [10] were also genotyped. The selected SNPs (Supplementary Table 1) were genotyped using Sequenom’s iPLEX Gold. Briefly, PCR and extension primers for these SNPs were designed using the MassARRAY Assay Design 3.0 software (Sequenom, Inc.). PCR amplification and single base extension reactions were performed according to the manufacturer’s instructions. Extension product sizes were determined by mass spectrometry using Sequenom’s Compact matrix-assisted laser desorption ionization-time of flight mass spectrometer. The resulting mass spectra were converted to genotype data using SpectroTYPER-RT software.

Statistical analysis

Time-at-risk for each MMR gene mutation carrier started at birth and ended at age at diagnosis of colorectal cancer (n=392), any other cancer (n=141), polypectomy (n=282), death (n=4) or last contact (n=279), whichever occurred first.

Since some carriers were ascertained because they were diagnosed with CRC or they had a strong family history, the identification of MMR gene mutation carriers was not random with respect to the disease status. To adjust for this non-random ascertainment, we used the weighted cohort approach.[18] Age-specific incidence rates of CRC for MMR gene mutation carriers [19] were used to calculate sampling fractions to weight the proportion of affected and unaffected carriers in five-year age stratum so the proportion of affected carriers in each age group equalled that expected for mutation carriers in the population.

Cox proportional hazards regression analysis was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for associations between each of 23 hTERT SNPs and CRC risk for MMR gene mutation carriers. We estimated HRs separately for homozygous carriers of the minor allele (2 minor alleles) and heterozygous carriers of the minor allele (1 minor allele) versus non-carriers (0 minor allele); and we estimated HRs per minor allele, i.e. a linear association on the log scale. As a subgroup analysis to investigate the associations between hTERT SNPs and early-onset CRC risk, we censored carriers by age of 45 years and conducted regressions.

The proportional hazards assumption was tested by examining the relationship between the scaled Schoenfeld residuals and survival time. To allow for any correlation of risk between family members, the Huber-White robust variance correction was applied by clustering on family membership [20]. To reduce the false discovery rate expected from the large number of associations investigated, the P value cut-off for classifying a HR as statistically significant was determined using methods by Benjamini and Hochberg [21]. We put the individual P-values in order, from smallest to largest. The smallest P value has a rank of i=1, then next smallest has i=2, etc. Then we compared each individual P value to its Benjamini-Hochberg critical value, (i/m)Q, where i is the rank, m is the total number of tests, and Q is the false discovery rate of 0.05. The largest P value that has P<(i/m)Q is significant, and all of the P values smaller than it are also significant. This method controls the expected high false discovery rate and can result in significant gains in power over traditional multiplicity ‘correction’ methods such as the commonly used Bonferroni procedure [22]. All statistical analyses were performed using STATA 13.0 (College Station, TX: StataCorp LP, 2013).

Haplotypes for the 23 SNPs within the hTERT gene were estimated using the haplologit command in STATA that implements the retrospective profile-likelihood methods of Spinka et al [23]. Logistic regression was used to examine associations between the haplotypes and CRC.

FUNDING

This work was supported by grant UM1 CA167551 from the National Cancer Institute and through cooperative agreements with Australasian Colorectal Cancer Family Registry (U01 CA074778 and U01/U24 CA097735). Aung Ko Win is an Australian National Health and Medical Council (NHMRC) Early Career Fellow. Melissa C. Southey is a NHMRC Senior Research Fellow. Mark A. Jenkins is a NHMRC Senior Research Fellow. John L. Hopper is a NHMRC Senior Principal Research Fellow and Distinguished Visiting Professor at Seoul National University, Korea. Christophe Rosty is the Jass Pathology Fellow. Daniel D. Buchanan is a University of Melbourne Research at Melbourne Accelerator Program (R@MAP) Senior Research Fellow. This project was conducted under Colon-CFR approval C-AU-0208-01.

Disclaimer

The content of this manuscript does not necessarily reflect the views or policies of the National Cancer Institute or any of the collaborating centers in the CFRs, nor does mention of trade names, commercial products, or organizations imply endorsement by the US Government or the CFR. Authors had full responsibility for the design of the study, the collection of the data, the analysis and interpretation of the data, the decision to submit the manuscript for publication, and the writing of the manuscript.

CONFLICT OF INTEREST

The authors have no conflict of interest to declare with respect to this manuscript.

AUTHORS’ CONTRIBUTIONS

Aung Ko Win, Daniel D. Buchanan: study concept and design; acquisition of data; statistical analysis; interpretation of data; drafting and critical review of the manuscript for important intellectual content; approval of the final version of the manuscript

Mark Clendenning, William Crawford, Christophe Rosty, Susan G. Preston: acquisition of data; interpretation of data; critical review of the manuscript for important intellectual content; approval of the final version of the manuscript

Melissa C. Southey, Susan Parry, Graham G. Giles, Finlay A. Macrae, Ingrid M. Winship, John A. Baron, John L. Hopper, Mark A. Jenkins: interpretation of data; critical review of the manuscript for important intellectual content; approval of the final version of the manuscript.

• 1. ten Broeke SW, Brohet RM, Tops CM, van der Klift HM, Velthuizen ME, Bernstein I, Capella Munar G, Gomez Garcia E, Hoogerbrugge N, Letteboer TG, Menko FH, Lindblom A, Mensenkamp AR, Moller P, van Os TA, Rahner N, et al. Lynch syndrome caused by germline PMS2 mutations: delineating the cancer risk. J Clin Oncol. 2015; 33(4):319-325. [PubMed]
• 2. Bonadona V, Bonaiti B, Olschwang S, Grandjouan S, Huiart L, Longy M, Guimbaud R, Buecher B, Bignon YJ, Caron O, Colas C, Nogues C, Lejeune-Dumoulin S, Olivier-Faivre L, Polycarpe-Osaer F, Nguyen TD, et al. Cancer risks associated with germline mutations in MLH1, MSH2, and MSH6 genes in Lynch syndrome. JAMA. 2011; 305(22):2304-2310.
• 3. van Duijnhoven FJ, Botma A, Winkels R, Nagengast FM, Vasen HF and Kampman E. Do lifestyle factors influence colorectal cancer risk in Lynch syndrome? Fam Cancer. 2013.. 2011; 305(22):2304-2310. [PubMed]
• 4. Talseth-Palmer BA, Wijnen JT, Brenne IS, Jagmohan-Changur S, Barker D, Ashton KA, Tops CM, Evans TJ, McPhillips M, Groombridge C, Suchy J, Kurzawski G, Spigelman A, Moller P, Morreau HM, Van Wezel T, et al. Combined analysis of three lynch syndrome cohorts confirms the modifying effects of 8q23.3 and 11q23.1 in MLH1 mutation carriers. Int J Cancer. 2013; 132(7):15561564. [PubMed]
• 5. Wijnen JT, Brohet RM, van Eijk R, Jagmohan-Changur S, Middeldorp A, Tops CM, van Puijenbroek M, Ausems MG, Gomez Garcia E, Hes FJ, Hoogerbrugge N, Menko FH, van Os TA, Sijmons RH, Verhoef S, Wagner A, et al. Chromosome 8q23.3 and 11q23.1 variants modify colorectal cancer risk in Lynch syndrome. Gastroenterology. 2009; 136(1):131-137. [PubMed]
• 6. Baird DM. Variation at the TERT locus and predisposition for cancer. Expert Rev Mol Med. 2010; 12:e16. [PubMed]
• 7. 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, 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, 384e371-372. [PubMed] https://doi.org/10.1038/ng.2566.
• 8. Kinnersley B, Migliorini G, Broderick P, Whiffin N, Dobbins SE, Casey G, Hopper J, Sieber O, Lipton L, Kerr DJ, Dunlop MG, Tomlinson IP, Houlston RS and Colon Cancer Family R. The TERT variant rs2736100 is associated with colorectal cancer risk. Br J Cancer. 2012; 107(6):1001-1008.
• 9. Peters U, Hutter CM, Hsu L, Schumacher FR, Conti DV, Carlson CS, Edlund CK, Haile RW, Gallinger S, Zanke BW, Lemire M, Rangrej J, Vijayaraghavan R, Chan AT, Hazra A, Hunter DJ, et al. Meta-analysis of new genome-wide association studies of colorectal cancer risk. Hum Genet. 2012; 131(2):217-234.
• 10. Bellido F, Guino E, Jagmohan-Changur S, Segui N, Pineda M, Navarro M, Lazaro C, Blanco I, Vasen HF, Moreno V, Capella G, Wijnen JT and Valle L. Genetic variant in the telomerase gene modifies cancer risk in Lynch syndrome. Eur J Hum Genet. 2013; 21(5):511-516. [PubMed] https://doi.org/10.1038/ejhg.2012.204.
• 12. Talseth-Palmer BA, Brenne IS, Ashton KA, Evans TJ, McPhillips M, Groombridge C, Suchy J, Kurzawski G, Spigelman A, Lubinski J and Scott RJ. Colorectal cancer susceptibility loci on chromosome 8q23.3 and 11q23.1 as modifiers for disease expression in Lynch syndrome. Journal of medical genetics. 2011; 48(4):279-284. [PubMed]
• 13. Talseth-Palmer BA, Scott RJ, Vasen HF and Wijnen JT. 8q23.3 and 11q23.1 as modifying loci influencing the risk for CRC in Lynch syndrome. European journal of human genetics : EJHG. 2012; 20(5):487-488; author reply 488. [PubMed] https://doi.org/10.1038/ejhg.2011.232.
• 14. Houlle S, Charbonnier F, Houivet E, Tinat J, Buisine MP, Caron O, Benichou J, Baert-Desurmont S and Frebourg T. Evaluation of Lynch syndrome modifier genes in 748 MMR mutation carriers. Eur J Hum Genet. 2011; 19(8):887-892. [PubMed] https://doi.org/10.1038/ejhg.2011.44.
• 15. Winship I and Win AK. The Australasian Colorectal Cancer Family Registry. Med J Aust. 2012; 197(9):480-481. [PubMed]
• 16. Senter L, Clendenning M, Sotamaa K, Hampel H, Green J, Potter JD, Lindblom A, Lagerstedt K, Thibodeau SN, Lindor NM, Young J, Winship I, Dowty JG, White DM, Hopper JL, Baglietto L, et al. The Clinical Phenotype of Lynch Syndrome Due to Germ-Line PMS2 Mutations. Gastroenterology. 2008; 135(2):419-428.
• 17. Pellatt AJ, Wolff RK, Herrick J, Lundgreen A and Slattery ML. TERT’s role in colorectal carcinogenesis. Mol Carcinog. 2013; 52(7):507-513.
• 18. Antoniou AC, Goldgar DE, Andrieu N, Chang-Claude J, Brohet R, Rookus MA and Easton DF. A weighted cohort approach for analysing factors modifying disease risks in carriers of high-risk susceptibility genes. Genet Epidemiol. 2005; 29(1):1-11. [PubMed]
• 19. Pande M, Lynch PM, Hopper JL, Jenkins MA, Gallinger S, Haile RW, LeMarchand L, Lindor NM, Campbell PT, Newcomb PA, Potter JD, Baron JA, Frazier ML and Amos CI. Smoking and Colorectal Cancer in Lynch Syndrome: Results from the Colon Cancer Family Registry and The University of Texas M.D. Anderson Cancer Center. Clin Cancer Res. 2010; 16(4):1331-1339. [PubMed] https://doi.org/10.1158/1078-0432.CCR-09-1877.
• 20. Williams RL. A Note on Robust Variance Estimation for Cluster-Correlated Data. Biometrics. 2000; 56(2):645-646. [PubMed]
• 21. Benjamini Y and Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc Ser B. 1995; 57(1):289-300.
• 22. Miller RG. (1981). Simultaneous statistical inference. (New York: Springer-Verlag).
• 23. Spinka C, Carroll RJ and Chatterjee N. Analysis of case-control studies of genetic and environmental factors with missing genetic information and haplotype-phase ambiguity. Genet Epidemiol. 2005; 29(2):108-127.