The admixture maximum likelihood test to test for association between rare variants and disease phenotypes.

BACKGROUND: The development of genotyping arrays containing hundreds of thousands of rare variants across the genome and advances in high-throughput sequencing technologies have made feasible empirical genetic association studies to search for rare disease susceptibility alleles. As single variant testing is underpowered to detect associations, the development of statistical methods to combine analysis across variants - so-called "burden tests" - is an area of active research interest. We previously developed a method, the admixture maximum likelihood test, to test multiple, common variants for association with a trait of interest. We have extended this method, called the rare admixture maximum likelihood test (RAML), for the analysis of rare variants. In this paper we compare the performance of RAML with six other burden tests designed to test for association of rare variants. RESULTS: We used simulation testing over a range of scenarios to test the power of RAML compared to the other rare variant association testing methods. These scenarios modelled differences in effect variability, the average direction of effect and the proportion of associated variants. We evaluated the power for all the different scenarios. RAML tended to have the greatest power for most scenarios where the proportion of associated variants was small, whereas SKAT-O performed a little better for the scenarios with a higher proportion of associated variants. CONCLUSIONS: The RAML method makes no assumptions about the proportion of variants that are associated with the phenotype of interest or the magnitude and direction of their effect. The method is flexible and can be applied to both dichotomous and quantitative traits and allows for the inclusion of covariates in the underlying regression model. The RAML method performed well compared to the other methods over a wide range of scenarios. Generally power was moderate in most of the scenarios, underlying the need for large sample sizes in any form of association testing.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

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

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

BRCA1 and BRCA2 mutations in a population-based study of male breast cancer

Background: The contribution of BRCA1 and BRCA2 to the incidence of male breast cancer (MBC) in the United Kingdom is not known, and the importance of these genes in the increased risk of female breast cancer associated with a family history of breast cancer in a male first-degree relative is unclear. Methods: We have carried out a population-based study of 94 MBC cases collected in the UK. We screened genomic DNA for mutations in BRCA1 and BRCA2 and used family history data from these cases to calculate the risk of breast cancer to female relatives of MBC cases. We also estimated the contribution of BRCA1 and BRCA2 to this risk. Results: Nineteen cases (20%) reported a first-degree relative with breast cancer, of whom seven also had an affected second-degree relative. The breast cancer risk in female first-degree relatives was 2.4 times (95% confidence interval [CI] = 1.4–4.0) the risk in the general population. No BRCA1 mutation carriers were identified and five cases were found to carry a mutation in BRCA2. Allowing for a mutation detection sensitivity frequency of 70%, the carrier frequency for BRCA2 mutations was 8% (95% CI = 3–19). All the mutation carriers had a family history of breast, ovarian, prostate or pancreatic cancer. However, BRCA2 accounted for only 15% of the excess familial risk of breast cancer in female first-degree relatives. Conclusion: These data suggest that other genes that confer an increased risk for both female and male breast cancer have yet to be found

Combined effects of single nucleotide polymorphisms TP53 R72P and MDM2 SNP309, and p53 expression on survival of breast cancer patients.

INTRODUCTION: Somatic inactivation of the TP53 gene in breast tumors is a marker for poor outcome, and breast cancer outcome might also be affected by germ-line variation in the TP53 gene or its regulators. We investigated the effects of the germ-line single nucleotide polymorphisms TP53 R72P (215G>C) and MDM2 SNP309 (-410T>G), and p53 protein expression in breast tumors on survival. METHODS: We pooled data from four breast cancer cohorts within the Breast Cancer Association Consortium for which both TP53 R72P and MDM2 SNP309 were genotyped and follow-up was available (n = 3,749). Overall and breast cancer-specific survival analyses were performed using Kaplan-Meier analysis and multivariate Cox's proportional hazards regression models. RESULTS: Survival of patients did not differ by carriership of either germ-line variant, R72P (215G>C) or SNP309 (-410G>T) alone. Immunohistochemical p53 staining of the tumor was available for two cohorts (n = 1,109 patients). Survival was worse in patients with p53-positive tumors (n = 301) compared to patients with p53-negative tumors (n = 808); breast cancer-specific survival: HR 1.6 (95% CI 1.2 to 2.1), P = 0.001. Within the patient group with p53-negative tumors, TP53 rare homozygous (CC) carriers had a worse survival than G-allele (GG/GC) carriers; actuarial breast cancer-specific survival 71% versus 80%, P = 0.07; HR 1.8 (1.1 to 3.1), P = 0.03. We also found a differential effect of combinations of the two germ-line variants on overall survival; homozygous carriers of the G-allele in MDM2 had worse survival only within the group of TP53 C-allele carriers; actuarial overall survival (GG versus TT/TG) 64% versus 75%, P = 0.001; HR (GG versus TT) 1.5 (1.1 to 2.0), P = 0.01. We found no evidence for a differential effect of MDM2 SNP309 by p53 protein expression on survival. CONCLUSIONS: The TP53 R72P variant may be an independent predictor for survival of patients with p53-negative tumors. The combined effect of TP53 R72P and MDM2 SNP309 on survival is in line with our a priori biologically-supported hypothesis, that is, the role of enhanced DNA repair function of the TP53 Pro-variant, combined with increased expression of the Mdm2 protein, and thus overall attenuation of the p53 pathway in the tumor cells.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

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

Evidence of a Causal Association Between Insulinemia and Endometrial Cancer: A Mendelian Randomization Analysis.

BACKGROUND: Insulinemia and type 2 diabetes (T2D) have been associated with endometrial cancer risk in numerous observational studies. However, the causality of these associations is uncertain. Here we use a Mendelian randomization (MR) approach to assess whether insulinemia and T2D are causally associated with endometrial cancer. METHODS: We used single nucleotide polymorphisms (SNPs) associated with T2D (49 variants), fasting glucose (36 variants), fasting insulin (18 variants), early insulin secretion (17 variants), and body mass index (BMI) (32 variants) as instrumental variables in MR analyses. We calculated MR estimates for each risk factor with endometrial cancer using an inverse-variance weighted method with SNP-endometrial cancer associations from 1287 case patients and 8273 control participants. RESULTS: Genetically predicted higher fasting insulin levels were associated with greater risk of endometrial cancer (odds ratio [OR] per standard deviation = 2.34, 95% confidence internal [CI] = 1.06 to 5.14, P = .03). Consistently, genetically predicted higher 30-minute postchallenge insulin levels were also associated with endometrial cancer risk (OR = 1.40, 95% CI = 1.12 to 1.76, P = .003). We observed no associations between genetic risk of type 2 diabetes (OR = 0.91, 95% CI = 0.79 to 1.04, P = .16) or higher fasting glucose (OR = 1.00, 95% CI = 0.67 to 1.50, P = .99) and endometrial cancer. In contrast, endometrial cancer risk was higher in individuals with genetically predicted higher BMI (OR = 3.86, 95% CI = 2.24 to 6.64, P = 1.2x10(-6)). CONCLUSION: This study provides evidence to support a causal association of higher insulin levels, independently of BMI, with endometrial cancer risk.This study was supported by MRC grant MC_UU_12015/1 and by the Innovative Medicines Initiative Joint Undertaking under EMIF grant agreement n° 115372 (contributions from the European Union's Seventh Framework Programme (FP7/2007-2013) and EFPIA companies). ANECS recruitment was supported by project grants from the National Health and Medical Research Council of Australia (ID#339435), The Cancer Council Queensland (ID#4196615) and Cancer Council Tasmania (ID#403031 and ID#457636). SEARCH recruitment was funded by a programme grant from Cancer Research UK [C490/A10124]. Case genotyping was supported by the National Health and Medical Research Council (ID#552402). Control data was generated by the Wellcome Trust Case Control Consortium (WTCCC), and a full list of the investigators who contributed to the generation of the data is available from the WTCCC website. We acknowledge use of DNA from the British 1958 Birth Cohort collection, funded by the Medical Research Council grant G0000934 and the Wellcome Trust grant 068545/Z/02. Funding for this project was provided by the Wellcome Trust under award 085475. Recruitment of the QIMR controls was supported by the National Health and Medical Research Council of Australia (NHMRC). The University of Newcastle, the Gladys M Brawn Senior Research Fellowship scheme, The Vincent Fairfax Family Foundation, the Hunter Medical Research Institute and the Hunter Area Pathology Service all contributed towards the costs of establishing the Hunter Community Study. K.T.N. was supported by the Gates Cambridge Trust. R.K.S. is supported by the Wellcome Trust (grant number WT098498). A.B.S. is supported by the National Health and Medical Research Council (NHMRC) Fellowship Scheme. D.F.E. is a Principal Research Fellow of Cancer Research UK. A.M.D is supported by the Joseph Mitchell Trust.This is the final version of the article. It first appeared from Oxford University Press via http://dx.doi.org/10.1093/jnci/djv17

A genome-wide association scan on estrogen receptor-negative breast cancer.

INTRODUCTION: Breast cancer is a heterogeneous disease and may be characterized on the basis of whether estrogen receptors (ER) are expressed in the tumour cells. ER status of breast cancer is important clinically, and is used both as a prognostic indicator and treatment predictor. In this study, we focused on identifying genetic markers associated with ER-negative breast cancer risk. METHODS: We conducted a genome-wide association analysis of 285,984 single nucleotide polymorphisms (SNPs) genotyped in 617 ER-negative breast cancer cases and 4,583 controls. We also conducted a genome-wide pathway analysis on the discovery dataset using permutation-based tests on pre-defined pathways. The extent of shared polygenic variation between ER-negative and ER-positive breast cancers was assessed by relating risk scores, derived using ER-positive breast cancer samples, to disease state in independent, ER-negative breast cancer cases. RESULTS: Association with ER-negative breast cancer was not validated for any of the five most strongly associated SNPs followed up in independent studies (1,011 ER-negative breast cancer cases, 7,604 controls). However, an excess of small P-values for SNPs with known regulatory functions in cancer-related pathways was found (global P = 0.052). We found no evidence to suggest that ER-negative breast cancer shares a polygenic basis to disease with ER-positive breast cancer. CONCLUSIONS: ER-negative breast cancer is a distinct breast cancer subtype that merits independent analyses. Given the clinical importance of this phenotype and the likelihood that genetic effect sizes are small, greater sample sizes and further studies are required to understand the etiology of ER-negative breast cancers.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