|Home | About | Journals | Submit | Contact Us | Français|
Telomeres cap chromosome ends and are critical for genomic stability. Many telomere-associated proteins are important for telomere length maintenance. Recent genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) in genes encoding telomere-associated proteins (RTEL1 and TERT-CLPTM1) as markers of cancer risk. We conducted an association study of telomere length and 743 SNPs in 43 telomere biology genes. Telomere length in peripheral blood DNA was determined by Q-PCR in 3,646 participants from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial and Nurses' Health Study. We investigated associations by SNP, gene, and pathway (functional group). We found no associations between telomere length and SNPs in TERT-CLPTM1L or RTEL1. Telomere length was not significantly associated with specific functional groups. Thirteen SNPs from four genes (MEN1, MRE11A, RECQL5, and TNKS) were significantly associated with telomere length. The strongest findings were in MEN1 (Gene-based P=0.006), menin, which associates with the telomerase promoter and may negatively regulate telomerase. This large association study did not find strong associations with telomere length. The combination of limited diversity and evolutionary conservation suggest that these genes may be under selective pressure. More work is needed to explore the role of genetic variants in telomere length regulation.
Telomeres form the ends of chromosomes and help maintain genomic structural integrity (Moon and Jarstfer, 2007). They consist of tandem hexameric (TTAGGG)n nucleotide repeats with a single-stranded overhang and protein complex. The overhang folds back to form a t-loop (Griffith et al., 1999), which prevents the telomere ends from being recognized as break points by the DNA damage repair machinery (Palm and de Lange, 2008). Many proteins bind to or interact with the telomere to maintain telomeric integrity. Shelterin is an ordered associated protein complex that consists of TERF1, TERF2, TINF2, TERF2IP, ACD, and POT1. This complex helps form the t-loop, and protects the telomeres from degradation and inappropriate DNA repair, thereby avoiding end-to-end fusion, atypical recombination, and premature senescence (Palm and de Lange, 2008). The telomerase reverse transcriptase (TERT) and its telomere template-containing RNA component (TERC) are telomere-associated proteins that add telomeric repeats to elongate telomeres (Collins and Mitchell, 2002). Telomerase activity is usually absent in differentiated cells. There are numerous other important telomere-associated proteins and complexes that transiently associates with telomeric DNA, including proteins involved in DNA repair (e.g., ATM and MRE11A) and helicases (e.g., BLM and RECQL) (Aubert and Lansdorp, 2008). A large number of additional proteins interact either directly or indirectly at telomeric ends, and regulate protein-protein and protein-DNA interactions, cellular protein trafficking, and additional telomere-specific functions.
Telomeric repeats range in size from 0.15 to 50 kilobases (kb), and progressively shorten with each cell division (Aubert and Lansdorp, 2008). Telomere length is dependent on many factors including age, replicative history of the cell, chromosome arm, and tissue type (Aubert and Lansdorp, 2008; Wise et al., 2009). There is considerable inter-individual variation in telomere length (Aviv et al., 2009; Nordfjäll et al., 2009), and a genetic influence on length variation. Twin and family studies have estimated the heritability of mean leukocyte telomere length to range from 44% to 84% (Jeanclos et al., 2000; Njajou et al., 2007; Slagboom et al., 1994; Vasa-Nicotera et al., 2004). Quantitative-trait linkage analyses have mapped loci influencing telomere length to chromosomes 12q12.22 and 14q23.2 (Andrew et al., 2006; Vasa-Nicotera et al., 2004). Fine mapping of the 12q12.22 locus identified a polymorphism in the BICD1 gene that was significantly associated with shorter telomere length (Mangino et al., 2008). A genome-wide association study (GWAS) identified two SNPs on chromosome 18q12.2 associated with leukocyte telomere length in the region of the gene VPS34/PIKC3C (Mangino et al., 2009), which has been suggested to be involved in controlling telomere length variation in yeast (Rog et al., 2005). Other recent GWASs have identified associations between leukocyte telomere length and the telomere-associated protein TERC (Codd et al., 2010; Levy et al., 2010) and a gene suggested to be involved with telomere length regulation, OBFC1 (oligonucleotide/oligosaccharide-binding folds containing one) (Levy et al., 2010).
The inverse relationships between telomere length and aging, and its role in age-related and premature aging diseases have been well documented (Aubert and Lansdorp, 2008; Garcia et al., 2007). Telomere attrition has also been associated with inflammatory processes, oxidative stress, and an unhealthy lifestyle (Mirabello et al., 2009; Morlá et al., 2006; von Zglinicki, 2002). There is growing evidence that short telomeres are associated with the initiation and progression of cancer (Blasco MA, 1997; Hackett and Greider, 2002).
Cancer GWAS have shown that SNPs in genes encoding telomere-associated proteins at 5p15.33 (TERT-CLPTM1L locus) and RTEL1 were associated with risk of glioma (Shete et al., 2009; Wrensch et al., 2009), pancreatic (Petersen et al., 2010), and/or lung cancer (Jin et al., 2009; Landi et al., 2009; McKay et al., 2008). In addition, an association study of multiple tumor types suggests that this TERT-CLPTM1L region may contain important markers of overall cancer risk (Rafnar et al., 2009). It is possible that these cancer-associated sequence variants in telomere-associated genes may be associated with shorter telomeres.
Studies of dyskeratosis congenita, an inherited bone marrow failure and cancer predisposition syndrome have also been important in understanding the consequences of telomere dysfunction (Savage and Alter, 2009). Patients with dyskeratosis congenita have extremely short telomeres for their age, very high risk of several cancers, and germline mutations in genes important in the maintenance of telomeres (DKC1, TERC, TERT, NOLA3, TINF2, or NOLA2) (Armanios, 2009; Savage and Alter, 2009). The phenotypic spectrum of telomere biology disorders also includes patients with isolated aplastic anemia (Yamaguchi et al., 2005), acute myelogenous leukemia (Calado et al., 2009), and idiopathic pulmonary fibrosis (Armanios et al., 2007) who may have mutations in TERC or TERT.
The majority of the associated and telomeric complex proteins are highly evolutionarily conserved (de Lange, 2004; Kanoh and Ishikawa, 2003; Li et al., 2000; Mirabello et al., 2008; Nakamura and Cech, 1998; Savage et al., 2005). A recent population genetics study targeting 37 telomere maintenance genes in 53 worldwide populations found that these genes have limited genetic variation (Mirabello et al., 2008). The majority of telomere genes had low diversity, high ancestral allele frequencies, and low population differentiation (Mirabello et al., 2008).
There is little known about how common genetic variation relates to telomere length. We hypothesize that common genetic variation in genes encoding telomere-associated proteins could affect telomere length. We evaluated the association between genetic variation in 43 candidate telomere biology genes and leukocyte telomere length using SNP markers from GWAS of breast and prostate cancers (Hunter et al., 2007; Yeager et al., 2007). These genes encode proteins that are thought to be involved either in telomere length maintenance or with telomere binding proteins necessary for telomere stability and structure. Telomere length data was obtained from prospectively collected blood samples and measured using quantitative-polymerase chain reaction (Q-PCR) by the same laboratory (De Vivo et al., 2009; Mirabello et al., 2009).
The PLCO Cancer Screening Trial is an ongoing randomized trial with 154,942 persons aged 55 to 74 enrolled between September 1993 and July 2001 from 10 screening centers nationwide (Prorok et al., 2000). Detailed questionnaire data was collected from all subjects at baseline. Participants provided blood and tissue samples for etiologic studies of cancer. Detailed information about the PLCO Trial and eligibility criteria are described elsewhere (Gohagan et al., 2000). Institutional review boards at the U.S. National Cancer Institute and the 10 screening centers approved the PLCO protocol, and all participants provided written informed consent. Participants in this study were male subjects selected for the prospective case-control study of telomere length, prostate cancer risk and life-style variables from Mirabello et al. (Mirabello et al., 2009). In brief, men for this study were selected from the screening arm of the PLCO trial who (a) were of non-Hispanic white race/ethnicity; (b) had a prostate cancer screen (PLCO PSA test) prior to October 1, 2003; (c) completed a baseline questionnaire about cancer risk factors; (d) provided a blood sample between one month and three years prior to prostate cancer diagnosis for cases; (e) age 55-74; and, (f) no history of cancer (other than non-melanoma skin cancer) to study entry. Medical and pathology records related to prostate cancer diagnosis were acquired for men with suspected prostate cancer by screening examination or annual questionnaire, and these data were abstracted by trained medical record specialists. All cases had confirmed aggressive prostate cancer and a Gleason score of ≥ 7. There were a total of 616 prostate cancer cases and 1061 matched male controls. Since there was no difference between the telomere length of cases and controls and no significant association with prostate cancer risk in this study (Mirabello et al., 2009), we combined telomere length data from the cases and controls for the current study.
The NHS is a cohort of 121,700 female registered nurses aged 30 to 55 from 11 states in the United States. Detailed questionnaires data was collected at enrollment and biennially thereafter, and blood samples were collected from a subset of 32,826 women. Details of this study are described elsewhere (De Vivo et al., 2009). Participants in the current study were female subjects selected for a case-control study of telomere length and postmenopausal breast cancer risk by De Vivo et al. (De Vivo et al., 2009). In short, eligible cases were (a) postmenopausal woman; (b) had confirmed incident invasive breast cancer diagnosed after blood collection up to June 1, 2004; (d) were of Caucasian ethnicity; and, (c) had no prior diagnosis of cancer. Controls were randomly selected postmenopausal women free of cancer and matched to cases on age, blood collection, and ethnicity. Completion of the questionnaire and submission of the blood sample were considered to imply informed consent. The study protocol was approved by the Committee on Use of Human Subjects of the Brigham and Women's Hospital, Boston, MA. There were 1,122 postmenopausal breast cancer cases and 1,147 controls included in the study. They also found no significant difference between the telomere length of cases and controls and no significant association with breast cancer risk (De Vivo et al., 2009), therefore, we included telomere length data from both cases and controls for the current study.
Details of the assay have been previously described (De Vivo et al., 2009; Mirabello et al., 2009). In brief, genomic DNA was extracted from peripheral blood samples (buffy coats, including all leukocytes) using the QIAmp (Qiagen, Chatsworth, CA) 96-spin blood protocol. The DNA concentration was quantified using a Nanodrop SD-1000 spectrophotometer, and subsequently dried down and re-suspended to ensure accurate and uniform DNA concentrations. Telomere length was measured by a modified version of the quantitative real time polymerase chain reaction (PCR)-based assay (Cawthon, 2002). The ratio of telomere repeat copy number to single-gene (β-globin, 36B4) copy number (T/S) was determined using an Applied Biosystems 7900HT PCR system in a high-throughput 384-well format. Five nanograms of genomic DNA was dried down and then re-suspended in 10 μL of either the telomere or 36B4 PCR reaction mixture. The telomere reaction mixture contained 1X Qiagen Quantitect Sybr Green Master Mix, 2.5 mM of DTT, 270 nM of Tel-1b primer (GGTTTTTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGT), and 900 nM of Tel-2b primer (TCCCGACTATCCCTATCCCTATCCCTATCCCTATCCCTA). The PCR reaction ran for 1 cycle at 95°C for 5 minutes, followed by 40 cycles at 95°C for 15 seconds, and 54°C for 2 minutes. The 36B4 reaction mixture consisted of 1X Qiagen Quantitect Sybr Green Master Mix, 300 nM of 36B4U primer (CAGCAAGTGGGAAGGTGTAATCC), and 500 nM of 36B4D primer (CCCATTCTATCATCAACGGGTACAA). The 36B4 PCR reaction ran for 1 cycle at 95°C for 5 minutes, followed by 40 cycles at 95°C for 15 seconds, and 58°C for 1 minute 10 seconds.
All samples of each assay were performed in triplicate. Each 384-well plate included an 8-point standard curve to assess and compensate for interpolated variations in PCR efficiency. Blinded replicate samples were interspersed with the samples to assess inter-plate and intra-plate variability of threshold cycle (Ct) values. The relative average telomere length was calculated as the ratio of telomere repeat copy number to single-gene copy number (T/S) in the study subjects compared with the reference DNA sample; derived by exponentiating the T/S ratio (−dCt) between the average telomere Ct and average 36B4 Ct values. For the PLCO study samples, the coefficients of variation (CV) within triplicates of the telomere and single-gene assay were 1.11% and 0.77%, respectively. For the NHS samples, the CVs of the telomere and single-gene assay were 1.03% and 0.56%, respectively.
We created a subset of genotype data from the GWAS of NHS (Hunter et al., 2007; Cancer Genetic Markers of Susceptibility project: http://cgems.cancer.gov) and PLCO (Yeager et al., 2007; Cancer Genetic Markers of Susceptibility project: http://cgems.cancer.gov) participants with telomere length data (De Vivo et al., 2009; Mirabello et al., 2009) for 43 genes (743 SNPs), including 20kb upstream and 10kb downstream, that have high a priori probabilities of effecting telomere length. We also included an expanded region around TERT-CLPTM1 due to recent findings in this region (Jin et al., 2009; Landi et al., 2009; McKay et al., 2008; Rafnar et al., 2009; Shete et al., 2009), from chromosome 5 position 1,201,710 to 1,445,536, which includes SLC6A19, SLC6A18, TERT, CLPTM1L, and SLC6A3. A detailed description of the genotyping methods and GWAS platforms are given elsewhere (Hunter et al., 2007; Yeager et al., 2007).
SNPs were excluded if they had less than a 90% genotyping rate in either study population or if they failed the Hardy-Weinberg equilibrium test. Individuals were excluded if they had more than 10% missing genotypes.
Linear regression models were used to estimate the association with telomere length for each SNP independently. Models were adjusted for smoking status (categorized as pack-years of smoking: 0, 1-9 pack-years, 10-29 pack-years, ≥30 pack-years), age (categorized as quartiles), case-control status (for both the PLCO and NHS participants), gender, NHS principal components and PLCO principal components (Patterson et al., 2006). For the principal components analysis, pooled case and control samples were analyzed using a set of 14,111 SNPs from the GWAS (Hunter et al., 2007; Yeager et al., 2007) with low pairwise linkage disequilibrium (LD; r2 < 0.01) according to the procedure described in Patterson et al. (2006). Testing for significance using the Tracy-Widom statistics (Patterson et al., 2006) showed four significant principal components for NHS and three for PLCO at the level of P < 0.05. The common allele was used as the referent category for the additive model to evaluate the additive effect of each minor allele. The gene–dose effects for each SNP were estimated by a linear trend test by coding the genotypes based on the number of variant alleles (0, 1 and 2). FDR corrections (Benjamini and Yekutieli, 2001) were conducted by gene (for all SNPs in a gene) for correction of multiple statistical tests.
We conducted gene-level and pathway-level analyses based on Yu et al. (Yu et al., 2009). The gene-level analysis is a global test for the association between the outcome and any subset of SNPs within a given gene or region. The pathway-level analysis is a global test for the association between the outcome and any subset of genes within a given pathway. P-values for these analyses were estimated with 20,000 permutation steps according to the algorithm.(Yu et al., 2009).
Statistical analyses were performed using SAS software, version 9.1 (SAS Institute, Cary, NC), R language, and PLINK software, version 1.06 (http://pngu.mgh.harvard.edu/purcell/plink/). Haploview version 4.1 (Barrett et al., 2005) was used to determine the degree of LD using SNP data from the HapMap Phase 2 (The International HapMap Consortium, 2003) public database.
Table 1 shows the characteristics of the study subjects, Table 2 shows the genes analyzed by functional pathway, and Supp. Table S1 lists the SNPs in each gene. We analyzed SNP data for 43 gene regions involved in telomere biology, in a dataset consisting of 743 SNPs in 3,646 individuals from the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial (Mirabello et al., 2009) and Nurses' Health Study (NHS) (De Vivo et al., 2009). There were 1635 males and 2011 females; their mean age was 60.8 [standard deviation (SD) of 6.5; range 43-74]. The mean relative telomere length (T/S) among the 3,646 subjects was 2.70 (SD 0.47) (Table 1).
We took three statistical approaches to evaluate genetic variation associated with telomere length, (1) at the individual SNP-level, (2) gene-level, and (3) pathway-level. We used the False Discovery Rate (FDR; Padj) correction to correct for multiple statistical tests at the SNP-level by gene (for all SNPs in a gene). Table 2 shows the most significant SNP in each gene. The SNP approach identified 13 SNPs that were significantly associated with telomere length after correction for multiple tests (Table 3). The direction of the regression coefficient (Beta) represents the effect of each minor allele increasing (+) or decreasing (−) telomere length. All four MEN1 SNPs were inversely associated with telomere length (top SNP, rs670358: Beta −0.046, standard error [SE] 0.017, Padj = 0.014). Two SNPs in MRE11A were positively associated with telomere length (rs12270338 and rs13447720: Beta 0.04, SE 0.012, Padj = 0.025). One SNP in RECQL5 and six SNPs in TNKS were also inversely associated with telomere length (Padj = 0.025-0.049; Table 3). The gene- and pathway-level analyses were global tests used to investigate the association between telomere length and all the SNPs within a gene and all genes within a pathway, respectively. The gene-level analysis identified 4 genes that were significantly associated with telomere length: MRE11A, RECQL5, MEN1, and TNKS (Table 2). MEN1 was the most significant gene (Gene P = 0.006). Genes were grouped by their function in telomere biology for the pathway-level analysis. No single functional pathway showed a significant association with telomere length (Table 2). Within functional pathways, the DNA repair and Other groups appeared to have the most SNPs with significant associations with telomere length, with no strong influence by any one group (Figure 1).
Since recent GWAS have suggested associations between SNPs at the RTEL1 and TERT-CLPTM1L locus and cancer (Jin et al., 2009; Landi et al., 2009; McKay et al., 2008; Petersen et al., 2010; Rafnar et al., 2009; Shete et al., 2009; Wrensch et al., 2009), we evaluated the five SNPs in these regions that were included in our dataset (Table 4). The genotype frequencies of these five SNPS did not differ between the cases and controls in our study population (P > 0.4). None of these SNPs were significantly associated with telomere length in our subjects. We also examined an extended region around the TERT-CLPTM1L locus on 5p13.33 to further investigate this region for any associations with telomere length. No SNPs in this region were significantly associated with telomere length after correction for multiple tests. Eight SNPs in TERT, CLPTM1L or the surrounding region were significantly associated without correction for multiple tests using co-dominant, dominant and/or recessive inheritance models (data not shown). The strongest association was observed for rs31489 in CLPTM1L with a recessive inheritance model (Beta 0.05, SE 0.02, P = 0.005).
There are many genes involved in telomere length maintenance and stability. We have examined the association between telomere length and genetic variation in 43 telomere biology genes among 3,646 individuals. This large dataset allowed us to investigate genetic variation in candidate genes associated with even small changes in telomere length. No specific functional pathway had a strong effect on telomere length. We noted several significant individual SNP associations after correction for multiple tests; however, we did not find compelling evidence that common variants in these pathways were associated with large differences in telomere length. Four genes (MRE11A, RECQL5, MEN1, and TNKS), and 13 SNPs within these genes had statistically significant effects on telomere length after correction at the gene level, although the effects were all relatively small.
We observed that variation in MRE11A (meiotic recombination 11) was positively associated with telomere length, and the two SNPs identified were in strong LD (r2 = 0.92). MRE11A is part of the telomere MRN (MRE11-RAD50-NBN) complex that senses DNA damage and is involved in modulating t-loop formation (Zhu et al., 2000). This complex positively regulates telomerase-dependent telomere elongation through interactions with ATM and TERF1 (Wu et al., 2007).
One SNP in RECQL5 (RECQ protein-like 5) was found to be inversely associated with telomere length. The RECQ protein family of helicases have critical roles in protecting and stabilizing the genome (Bohr, 2008). Deficiencies in these proteins can lead to genomic instability, premature aging and an increased susceptibility to cancer (Bohr, 2008). Detailed molecular functions of RECQL5 are not clear. RECQL5 has a role in preventing inappropriate homologous recombination (Hu et al., 2007), and RECQL5 deficient mice have an increased susceptibility to cancer.
MEN1 (multiple endocrine neoplasia I) was the most significantly implicated gene, with all SNPs included in this region showing similar inverse associations with telomere length. Surprisingly, there is very little LD across this gene region and these SNPs do not appear to be linked. MEN1 is a strong tumor suppressor gene that encodes menin. Menin associates with the promoter region of TERT, and is thought to negatively regulate TERT (Lin and Elledge, 2003). Through TERT, MEN1 may affect telomere length maintenance, since TERT expression has been shown to stimulate telomerase activity (Lin and Elledge, 2003; Rufer et al., 2001).
Six SNPs in TNKS (tankyrase 1) were similarly negatively associated with telomere length. TNKS is a poly (ADP-ribose) polymerase thought to be a positive regulator of telomere length through its modification of TERF1 (PARsylation), which releases TERF1 repression and allows telomerase to bind and elongate telomeric DNA (Cook et al., 2002; Hsiao and Smith, 2008; Smith and de Lange, 2000). In a small study of genetic variation in nine telomere-associated genes (ACD, POT1, TEP1, TERF1, TERF2, TERF2IP, TERT, TNKS, TNKS2) among 100 breast cancer cases, a SNP in TNKS that significantly associated with telomere length was identified (Varadi et al., 2009). Although this TNKS SNP was not included in our analysis, HapMap (The International HapMap Consortium, 2003) data suggests that this SNP (rs6990097) (Varadi et al., 2009) and one of our significant TNKS SNPs (rs11991621) may be in LD (D′ = 0.92 and r2 = 0.40). In addition, a recent GWAS identified two SNPs in TNKS associated with leukocyte telomere length with P-values of 10−2 to 10−4 (Mangino et al., 2009). One of these SNPs, rs11249943, was also statistically significant in our dataset. Our study is therefore the third to have identified significant associations between genetic variation in TNKS and telomere length.
We included an extended region around TERT-CLPTM1L since recent GWAS (Jin et al., 2009; Landi et al., 2009; McKay et al., 2008; Petersen et al., 2010; Rafnar et al., 2009; Shete et al., 2009) suggested that this region may contain important markers of overall cancer risk. We could not confirm a significant effect of these cancer-associated SNPs (Jin et al., 2009; Landi et al., 2009; McKay et al., 2008; Rafnar et al., 2009; Shete et al., 2009; Wrensch et al., 2009) on telomere length in our study population of 3646 subjects. Rafnar et al. (Rafnar et al., 2009) examined the association between rs401681 and rs2736098 and telomere length in DNA from whole blood, and found that only older subjects (n = 276; born 1925-1935) with the risk allele (rs401681 C allele, and rs2736098 A allele) had significantly shorter telomeres (P = 0.017 and 0.027, respectively) compared with younger women (n = 260; born 1940-1950). They suggested that these variants may gradually shorten telomeres over time with an effect becoming apparent at an older age (Rafnar et al., 2009). The mean age of individuals in our study was 60.8 (age range 43-73), so we did not stratify on age further. The functional significance of these SNP associations with cancer risk requires further investigation.
Variants in intron 1 of BICD1 (most notably rs2630578) were associated with telomere length (Mangino et al., 2008). Our dataset did not include rs2630578 but did include rs1798255, rs2668301, rs10771917 which are in LD with rs2630578. These SNPs in BICD1 were not significantly associated with telomere length after correction for multiple tests in our dataset. Another linkage study found significant linkage of telomere length to a region on chromosome 12 including the DNA helicase DDX11; however, after further analysis of the SNPs in this region, they found no significant effects of DDX11 genotypes on telomere length (Vasa-Nicotera et al., 2004); we observed no significant associations between SNPs in DDX11 and telomere length. Our dataset only included one SNP (rs10936599) in TERC, therefore, we could not determine if the SNP (rs12696304) identified by Codd et al. (2010) in TERC was specifically associated with telomere length. However, these two SNPs appear to be in high LD (D′ = 1.0 and r2 = 0.91), and we could not confirm a significant association in our dataset (P = 0.17). Another recent GWAS identified SNPs in the region of the OBFC1 gene, suggested to be involved with the replication and capping of telomeres, associated with leukocyte telomere length, and they also confirmed an association between TERC and telomere length (Levy et al., 2010). The inconsistencies between many of these studies could be due to differences in genetic background of the study populations, limitations in statistical power, differing telomere lengths in the different cell types analyzed, and/or differences in data interpretation methods. In addition, common SNPs may not have large effects on telomere length. It is also important to point out that our findings are for blood leukocyte DNA and may not necessarily apply to other tissue DNA.
A limitation of our study was that the SNPs available for analysis were limited to those in the specific genotyping platform used for the previous GWAS studies from which our individual data were drawn. This resulted in limited representation of SNPs in some regions; e.g., we only had one SNP in TERC. Strengths of our study include our comprehensive analysis of SNPs in the vast majority of known telomere-associated genes using 3 statistical approaches (at the individual SNP-level, gene-level and pathway-level), the analysis of telomere length in peripheral blood samples by Q-PCR from two studies done by the same laboratory, and the large sample size (3646 individuals).
The relatively limited effect of common germline genetic variants on telomere length could be due to the fact that the majority of these genes are highly evolutionarily conserved (de Lange, 2004; Kanoh and Ishikawa, 2003; Li et al., 2000; Mirabello et al., 2008; Nakamura and Cech, 1998; Savage et al., 2005). Our recent population genetic study of the majority of telomere maintenance genes using a similar approach of extracting genotype data from a GWA platform identified little variation in these genes as a group (Mirabello et al., 2008). MEN1 was not included in that gene-set, but the SNPs in TNKS, RECQL5, and MRE11A all had low population differentiation [FST ≤0.06; genome-wide average for autosomal SNPs is 0.10 ~ 0.15 (Akey et al., 2002; Shriver et al., 2004; Shriver et al., 2005; Weir et al., 2005)] and heterozygosity (0.22-0.37) (Mirabello et al., 2008). There was also evidence of evolutionary selection in TNKS (Mirabello et al., 2008; Savage et al., 2005; The International HapMap Consortium, 2007). A large amount of variation in telomere genes may not be tolerated due to their critical roles in telomere maintenance and chromosomal stability. Rare variants (i.e., mutations) may only contribute to rare telomere biology disorders, such as dyskeratosis congenita. Telomere length and the rate of telomere length change are highly variable among individuals, and have been shown to be effected by many factors including race, BMI, smoking, inflammation, and oxidative stress (Aviv et al., 2009; Mirabello et al., 2009; Nordfjäll et al., 2009). It may be that telomere length is controlled by small contributions from many environmental factors and/or polymorphisms.
In summary, this study of common genetic variation in candidate genes and telomere length identified SNPs in four genes that are associated with telomere length after correction for multiple tests by gene. However, the effects on telomere length were found to be small, and we found no association using global tests by functional pathway. Common SNPs, like those used in GWAS, may not have strong effects on telomere length. The combination of limited nucleotide diversity, high ancestral allele frequencies, and evolutionary conservation suggest that these genes may be under selective pressure. Additional studies of rare variants will be helpful in understanding the role of genetic variation in telomere length regulation.
We are grateful to the NHS and PLCO participants for their valuable contributions. The work of Drs. Mirabello, Yu, Chatterjee, Hayes and Savage was supported [in part] by the Intramural Research Program of the National Institutes of Health and the National Cancer Institute.
Supporting Information for this preprint is available from the Human Mutation editorial office upon request (firstname.lastname@example.org)