Telomere length and common disease: study design and analytical challenges.

Telomeres, the repetitive sequences that protect the ends of chromosomes, help to maintain genomic integrity and are of key importance to human health. The aim here is to give an overview of the evidence for the importance of telomere length (TL) to the risk of common disease, considering the strengths and weaknesses of different epidemiological study designs. Methods for measuring TL are described, all of which are subject to considerable measurement error. TL declines with age and varies in relation to factors such as smoking and obesity. It is also highly heritable (estimated heritability of ~40 to 50%), and genome-wide studies have identified a number of associated genetic variants. Epidemiological studies have shown shorter TL to be associated with risk of a number of common diseases, including cardiovascular disease and some cancers. The relationship with cancer appears complex, in that longer telomeres are associated with higher risk of some cancers. Prospective studies of the relationship between TL and disease, where TL is measured before diagnosis, have numerous advantages over retrospective studies, since they avoid the problems of reverse causality and differences in sample handling, but they are still subject to potential confounding. Studies of the genetic predictors of TL in relation to disease risk avoid these drawbacks, although they are not without limitations. Telomere biology is of major importance to the risk of common disease, but the complexities of the relationship are only now beginning to be understood.This research was supported by Cancer Research UK Programme Awards C588/A10589 and C588/A19167 (MMI and JHB) and C8197/A16565 (AMD and KAP) and the Isaac Newton Trust.This is the final version of the article. It first appeared from Springer via http://dx.doi.org/10.1007/s00439-015-1563-

A genome-wide association study of radiotherapy induced toxicity in head and neck cancer patients identifies a susceptibility locus associated with mucositis

PURPOSE: A two-stage genome-wide association study was carried out in head and neck cancer (HNC) patients aiming to identify genetic variants associated with either specific radiotherapy-induced (RT) toxicity endpoints or a general proneness to develop toxicity after RT.MATERIALS AND METHODS: The analysis included 1780 HNC patients treated with primary RT for laryngeal or oro/hypopharyngeal cancers. In a non-hypothesis-driven explorative discovery study, associations were tested in 1183 patients treated within The Danish Head and Neck Cancer Group. Significant associations were later tested in an independent Dutch cohort of 597 HNC patients and if replicated, summary data obtained from discovery and replication studies were meta-analysed. Further validation of significantly replicated findings was pursued in an Asian cohort of 235 HNC patients with nasopharynx as the primary tumour site.RESULTS: We found and replicated a significant association between a locus on chromosome 5 and mucositis with a pooled OR for rs1131769*C in meta-analysis = 1.95 (95% CI 1.48-2.41; ppooled = 4.34 × 10-16).CONCLUSION: This first exploratory GWAS in European cohorts of HNC patients identified and replicated a risk locus for mucositis. A larger Meta-GWAS to identify further risk variants for RT-induced toxicity in HNC patients is warranted.</p

Seq4SNPs: new software for retrieval of multiple, accurately annotated DNA sequences, ready formatted for SNP assay design.

BACKGROUND: In moderate-throughput SNP genotyping there was a gap in the workflow, between choosing a set of SNPs and submitting their sequences to proprietary assay design software, which was not met by existing software. Retrieval and formatting of sequences flanking each SNP, prior to assay design, becomes rate-limiting for more than about ten SNPs, especially if annotated for repetitive regions and adjacent variations. We routinely process up to 50 SNPs at once. IMPLEMENTATION: We created Seq4SNPs, a web-based, walk-away software that can process one to several hundred SNPs given rs numbers as input. It outputs a file of fully annotated sequences formatted for one of three proprietary design softwares: TaqMan's Primer-By-Design FileBuilder, Sequenom's iPLEX or SNPstream's Autoprimer, as well as unannotated fasta sequences. We found genotyping assays to be inhibited by repetitive sequences or the presence of additional variations flanking the SNP under test, and in multiplexes, repetitive sequence flanking one SNP adversely affects multiple assays. Assay design software programs avoid such regions if the input sequences are appropriately annotated, so we used Seq4SNPs to provide suitably annotated input sequences, and improved our genotyping success rate. Adjacent SNPs can also be avoided, by annotating sequences used as input for primer design. CONCLUSION: The accuracy of annotation by Seq4SNPs is significantly better than manual annotation (P < 1e-5).Using Seq4SNPs to incorporate all annotation for additional SNPs and repetitive elements into sequences, for genotyping assay designer software, minimizes assay failure at the design stage, reducing the cost of genotyping. Seq4SNPs provides a rapid route for replacement of poor test SNP sequences. We routinely use this software for assay sequence preparation. Seq4SNPs is available as a service at (http://moya.srl.cam.ac.uk/oncology/bio/s4shome.html) and (http://moya.srl.cam.ac.uk/cgi-bin/oncology/srl/ncbi/seq4snp1.pl), currently for human SNPs, but easily extended to include any species in dbSNP.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Validation of loci at 2q14.2 and 15q21.3 as risk factors for testicular cancer.

Testicular germ cell tumor (TGCT), the most common cancer in men aged 18 to 45 years, has a strong heritable basis. Genome-wide association studies (GWAS) have proposed single nucleotide polymorphisms (SNPs) at a number of loci influencing TGCT risk. To further evaluate the association of recently proposed risk SNPs with TGCT at 2q14.2, 3q26.2, 7q36.3, 10q26.13 and 15q21.3, we analyzed genotype data on 3,206 cases and 7,422 controls. Our analysis provides independent replication of the associations for risk SNPs at 2q14.2 (rs2713206 at P = 3.03 × 10-2; P-meta = 3.92 × 10-8; nearest gene, TFCP2L1) and rs12912292 at 15q21.3 (P = 7.96 × 10-11; P-meta = 1.55 × 10-19; nearest gene PRTG). Case-only analyses did not reveal specific associations with TGCT histology. TFCP2L1 joins the growing list of genes located within TGCT risk loci with biologically plausible roles in developmental transcriptional regulation, further highlighting the importance of this phenomenon in TGCT oncogenesis

Association between Common Variation in 120 Candidate Genes and Breast Cancer Risk

Association studies in candidate genes have been widely used to search for common low penetrance susceptibility alleles, but few definite associations have been established. We have conducted association studies in breast cancer using an empirical single nucleotide polymorphism (SNP) tagging approach to capture common genetic variation in genes that are candidates for breast cancer based on their known function. We genotyped 710 SNPs in 120 candidate genes in up to 4,400 breast cancer cases and 4,400 controls using a staged design. Correction for population stratification was done using the genomic control method, on the basis of data from 280 genomic control SNPs. Evidence for association with each SNP was assessed using a Cochran–Armitage trend test (p-trend) and a two-degrees of freedom χ(2) test for heterogeneity (p-het). The most significant single SNP (p-trend = 8 × 10(−5)) was not significant at a nominal 5% level after adjusting for population stratification and multiple testing. To evaluate the overall evidence for an excess of positive associations over the proportion expected by chance, we applied two global tests: the admixture maximum likelihood (AML) test and the rank truncated product (RTP) test corrected for population stratification. The admixture maximum likelihood experiment-wise test for association was significant for both the heterogeneity test (p = 0.0031) and the trend test (p = 0.017), but no association was observed using the rank truncated product method for either the heterogeneity test or the trend test (p = 0.12 and p = 0.24, respectively). Genes in the cell-cycle control pathway and genes involved in steroid hormone metabolism and signalling were the main contributors to the association. These results suggest that a proportion of SNPs in these candidate genes are associated with breast cancer risk, but that the effects of individual SNPs is likely to be small. Large sample sizes from multicentre collaboration will be needed to identify associated SNPs with certainty

Common variants in the ATM, BRCA1, BRCA2, CHEK2 and TP53 cancer susceptibility genes are unlikely to increase breast cancer risk.

INTRODUCTION: Certain rare, familial mutations in the ATM, BRCA1, BRCA2, CHEK2 or TP53 genes increase susceptibility to breast cancer but it has not, until now, been clear whether common polymorphic variants in the same genes also increase risk. METHODS: We have attempted a comprehensive, single nucleotide polymorphism (SNP)- and haplotype-tagging association study on each of these five genes in up to 4,474 breast cancer cases from the British, East Anglian SEARCH study and 4,560 controls from the EPIC-Norfolk study, using a two-stage study design. Nine tag SNPs were genotyped in ATM, together with five in BRCA1, sixteen in BRCA2, ten in CHEK2 and five in TP53, with the aim of tagging all other known, common variants. SNPs generating the common amino acid substitutions were specifically forced into the tagging set for each gene. RESULTS: No significant breast cancer associations were detected with any individual or combination of tag SNPs. CONCLUSION: It is unlikely that there are any other common variants in these genes conferring measurably increased risks of breast cancer in our study population.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Targeted Resequencing of the Coding Sequence of 38 Genes Near Breast Cancer GWAS Loci in a Large Case-Control Study.

BACKGROUND: Genes regulated by breast cancer risk alleles identified through genome-wide association studies (GWAS) may harbor rare coding risk alleles. METHODS: We sequenced the coding regions for 38 genes within 500 kb of 38 lead GWAS SNPs in 13,538 breast cancer cases and 5,518 controls. RESULTS: Truncating variants in these genes were rare, and were not associated with breast cancer risk. Burden testing of rare missense variants highlighted 5 genes with some suggestion of an association with breast cancer, although none met the multiple testing thresholds: MKL1, FTO, NEK10, MDM4, and COX11. Six common alleles in COX11, MAP3K1 (two), and NEK10 (three) were associated at the P < 0.0001 significance level, but these likely reflect linkage disequilibrium with causal regulatory variants. CONCLUSIONS: There was no evidence that rare coding variants in these genes confer substantial breast cancer risks. However, more modest effect sizes could not be ruled out. IMPACT: We tested the hypothesis that rare variants in 38 genes near breast cancer GWAS loci may mediate risk. These variants do not appear to play a major role in breast cancer heritability

Germline TERT promoter mutations are rare in familial melanoma.

Germline CDKN2A mutations occur in 40 % of 3-or-more case melanoma families while mutations of CDK4, BAP1, and genes involved in telomere function (ACD, TERF2IP, POT1), have also been implicated in melanomagenesis. Mutation of the promoter of the telomerase reverse transcriptase (TERT) gene (c.-57 T>G variant) has been reported in one family. We tested for the TERT promoter variant in 675 multicase families wild-type for the known high penetrance familial melanoma genes, 1863 UK population-based melanoma cases and 529 controls. Germline lymphocyte telomere length was estimated in carriers. The c.-57 T>G TERT promoter variant was identified in one 7-case family with multiple primaries and early age of onset (earliest, 15 years) but not among population cases or controls. One family member had multiple primary melanomas, basal cell carcinomas and a bladder tumour. The blood leukocyte telomere length of a carrier was similar to wild-type cases. We provide evidence confirming that a rare promoter variant of TERT (c.-57 T>G) is associated with high penetrance, early onset melanoma and potentially other cancers, and explains <1 % of UK melanoma multicase families. The identification of POT1 and TERT germline mutations highlights the importance of telomere integrity in melanoma biology.The authors would like to thank the families for their willingness to participate; and Rajiv Kumar for the provision of mutation positive samples. The collection of samples in the UK population-ascertained sample set was funded by Cancer Research UK (awards C588/A19167 and C8216/A6129) and by the NIH (CA83115). The work of N.A.G. and R.v.D was supported by the Dutch Cancer Society (UL 2012-5489). D.J.A and C.D.R.E are supported by Cancer Research UK, ERC Combat Cancer and the Wellcome Trust. N.K.H is supported by a fellowship from the National Health and Medical Research Council of Australia. A.M.D. and K.A.P. were supported by CRUK grant (C8197/A16565) and The Isaac Newton Trust. K.M.B. is supported by the Intramural Research Program of the Division of Cancer Epidemiology and Genetics; National Cancer Institute; National Institutes of Health.This is the final version of the article. It first appeared from Springer via http://dx.doi.org/10.1007/s10689-015-9841-

Genetic variation in stromal proteins decorin and lumican with breast cancer: investigations in two case-control studies.

INTRODUCTION: The stroma is the supportive framework of biologic tissue in the breast, consisting of various proteins such as the proteoglycans, decorin and lumican. Altered expression of decorin and lumican is associated with breast tumors. We hypothesized that genetic variation in the decorin (DCN) and lumican (LUM) genes may contribute to breast cancer. METHODS: We investigated associations of 14 common polymorphisms in the DCN and LUM genes with 798 breast cancer cases and 843 controls from Mayo Clinic, MN, USA. One polymorphism per gene with the strongest risk association in the Mayo Clinic sample was genotyped in 4,470 breast cancer cases and 4,560 controls from East Anglia, England (Studies of Epidemiology and Risk Factors in Cancer Heredity (SEARCH)). RESULTS: In the Mayo Clinic sample, six polymorphisms were associated with breast cancer risk (P trend <or= 0.05). The association with LUM rs2268578, evaluated further in SEARCH, was positive, although the odds ratios (OR) were weaker and not statistically significant. ORs were 1.4 (95% confidence interval [CI], 1.1 to 1.8) for heterozygotes and 2.2 (95% CI, 1.1 to 4.3; P2 df = 0.002) for homozygotes in the Mayo Clinic sample, and were 1.1 (95% CI, 0.9 to 1.2) for heterozygotes and 1.4 (95% CI, 1.0 to 2.1; P2 df = 0.13) for homozygotes in the SEARCH sample. In combined analyses, the ORs were 1.1 (95% CI, 1.0 to 1.2) for heterozygotes and 1.6 (95% CI, 1.2 to 2.3; P2 df = 0.005) for homozygotes. Positive associations for this polymorphism were observed for estrogen receptor-positive tumors in both the Mayo Clinic sample (OR for heterozygotes = 1.5, 1.1 to 1.9 and OR for homozygotes = 2.5, 1.2 to 5.3;P2 df = 0.001) and the SEARCH sample (OR for heterozygotes = 1.0, 0.9 to 1.1 and OR for homozygotes = 1.6, 1.0 to 2.5; P2 df = 0.10). In combined analyses, the ORs were 1.1 (95% CI, 0.9 to 1.2) for heterozygotes and 1.9 (95% CI, 1.3 to 2.8; P2 df = 0.001) for homozygotes. CONCLUSIONS: Although LUM rs2268578 was associated with breast cancer in the Mayo Clinic study, particularly estrogen receptor-positive breast cancer, weaker and modest associations were observed in the SEARCH sample. These modest associations will require larger samples to adequately assess the importance of this polymorphism in breast cancer

The scale of population structure in Arabidopsis thaliana

The population structure of an organism reflects its evolutionary history and influences its evolutionary trajectory. It constrains the combination of genetic diversity and reveals patterns of past gene flow. Understanding it is a prerequisite for detecting genomic regions under selection, predicting the effect of population disturbances, or modeling gene flow. This paper examines the detailed global population structure of Arabidopsis thaliana. Using a set of 5,707 plants collected from around the globe and genotyped at 149 SNPs, we show that while A. thaliana as a species self-fertilizes 97% of the time, there is considerable variation among local groups. This level of outcrossing greatly limits observed heterozygosity but is sufficient to generate considerable local haplotypic diversity. We also find that in its native Eurasian range A. thaliana exhibits continuous isolation by distance at every geographic scale without natural breaks corresponding to classical notions of populations. By contrast, in North America, where it exists as an exotic species, A. thaliana exhibits little or no population structure at a continental scale but local isolation by distance that extends hundreds of km. This suggests a pattern for the development of isolation by distance that can establish itself shortly after an organism fills a new habitat range. It also raises questions about the general applicability of many standard population genetics models. Any model based on discrete clusters of interchangeable individuals will be an uneasy fit to organisms like A. thaliana which exhibit continuous isolation by distance on many scales