Serena Scollen

A common coding variant in CASP8 is associated with breast cancer risk

The Breast Cancer Association Consortium (BCAC) has been established to conduct combined case-control analyses with augmented statistical power to try to confirm putative genetic associations with breast cancer. We genotyped nine SNPs for which there was some prior evidence of an association with breast cancer: CASP8 D302H (rs1045485), IGFBP3 −202 C → A (rs2854744), SOD2 V16A (rs1799725), TGFB1 L10P (rs1982073), ATM S49C (rs1800054), ADH1B 3′ UTR A → G (rs1042026), CDKN1A S31R (rs1801270), ICAM5 V301I (rs1056538) and NUMA1 A794G (rs3750913). We included data from 9–15 studies, comprising 11,391–18,290 cases and 14,753–22,670 controls. We found evidence of an association with breast cancer for CASP8 D302H (with odds ratios (OR) of 0.89 (95% confidence interval (c.i.): 0.85–0.94) and 0.74 (95% c.i.: 0.62–0.87) for heterozygotes and rare homozygotes, respectively, compared with common homozygotes; Ptrend = 1.1 × 10−7) and weaker evidence for TGFB1 L10P (OR = 1.07 (95% c.i.: 1.02–1.13) and 1.16 (95% c.i.: 1.08–1.25), respectively; Ptrend = 2.8 × 10−5). These results demonstrate that common breast cancer susceptibility alleles with small effects on risk can be identified, given sufficiently powerful studies.NOTE: In the version of this article initially published, there was an error that affected the calculations of the odds ratios, confidence intervals, between-study heterogeneity, trend test and test for association for SNP ICAM5 V301I in Table 1 (ICAM5 V301I); genotype counts in Supplementary Table 2 (ICAM5; ICR_FBCS and Kuopio studies) and minor allele frequencies, trend test and odds ratios for heterozygotes and rare homozygotes in Supplementary Table 3 (ICAM5; ICR_FBCS and Kuopio studies). The errors in Table 1 have been corrected in the PDF version of the article. The errors in supplementary information have been corrected online.

Large-Scale Analysis of Association Between LRP5 and LRP6 Variants and Osteoporosis

CONTEXT Mutations in the low-density lipoprotein receptor-related protein 5 (LRP5) gene cause rare syndromes characterized by altered bone mineral density (BMD). More common LRP5 variants may affect osteoporosis risk in the general population. OBJECTIVE To generate large-scale evidence on whether 2 common variants of LRP5 (Val667Met, Ala1330Val) and 1 variant of LRP6 (Ile1062Val) are associated with BMD and fracture risk. DESIGN AND SETTING Prospective, multicenter, collaborative study of individual-level data on 37,534 individuals from 18 participating teams in Europe and North America. Data were collected between September 2004 and January 2007; analysis of the collected data was performed between February and May 2007. Bone mineral density was assessed by dual-energy x-ray absorptiometry. Fractures were identified via questionnaire, medical records, or radiographic documentation; incident fracture data were available for some cohorts, ascertained via routine surveillance methods, including radiographic examination for vertebral fractures. MAIN OUTCOME MEASURES Bone mineral density of the lumbar spine and femoral neck; prevalence of all fractures and vertebral fractures. RESULTS The Met667 allele of LRP5 was associated with reduced lumbar spine BMD (n = 25,052 [number of participants with available data]; 20-mg/cm2 lower BMD per Met667 allele copy; P = 3.3 x 10(-8)), as was the Val1330 allele (n = 24,812; 14-mg/cm2 lower BMD per Val1330 copy; P = 2.6 x 10(-9)). Similar effects were observed for femoral neck BMD, with a decrease of 11 mg/cm2 (P = 3.8 x 10(-5)) and 8 mg/cm2 (P = 5.0 x 10(-6)) for the Met667 and Val1330 alleles, respectively (n = 25 193). Findings were consistent across studies for both LRP5 alleles. Both alleles were associated with vertebral fractures (odds ratio [OR], 1.26; 95% confidence interval [CI], 1.08-1.47 for Met667 [2001 fractures among 20 488 individuals] and OR, 1.12; 95% CI, 1.01-1.24 for Val1330 [1988 fractures among 20,096 individuals]). Risk of all fractures was also increased with Met667 (OR, 1.14; 95% CI, 1.05-1.24 per allele [7876 fractures among 31,435 individuals)]) and Val1330 (OR, 1.06; 95% CI, 1.01-1.12 per allele [7802 fractures among 31 199 individuals]). Effects were similar when adjustments were made for age, weight, height, menopausal status, and use of hormone therapy. Fracture risks were partly attenuated by adjustment for BMD. Haplotype analysis indicated that Met667 and Val1330 variants both independently affected BMD. The LRP6 Ile1062Val polymorphism was not associated with any osteoporosis phenotype. All aforementioned associations except that between Val1330 and all fractures and vertebral fractures remained significant after multiple-comparison adjustments. CONCLUSIONS Common LRP5 variants are consistently associated with BMD and fracture risk across different white populations. The magnitude of the effect is modest. LRP5 may be the first gene to reach a genome-wide significance level (a conservative level of significance [herein, unadjusted P < 10(-7)] that accounts for the many possible comparisons in the human genome) for a phenotype related to osteoporosis.

Large-Scale Evidence for the Effect of the COLIA1 Sp1 Polymorphism on Osteoporosis Outcomes: The GENOMOS Study

Background Osteoporosis and fracture risk are considered to be under genetic control. Extensive work is being performed to identify the exact genetic variants that determine this risk. Previous work has suggested that a G/T polymorphism affecting an Sp1 binding site in the COLIA1 gene is a genetic marker for low bone mineral density (BMD) and osteoporotic fracture, but there have been no very-large-scale studies of COLIA1 alleles in relation to these phenotypes. Methods and Findings Here we evaluated the role of COLIA1 Sp1 alleles as a predictor of BMD and fracture in a multicenter study involving 20,786 individuals from several European countries. At the femoral neck, the average (95% confidence interval [CI]) BMD values were 25 mg/cm 2 (CI, 16 to 34 mg/cm 2) lower in TT homozygotes than the other genotype groups ( p < 0.001), and a similar difference was observed at the lumbar spine; 21 mg/cm 2 (CI, 1 to 42 mg/cm 2), ( p = 0.039). These associations were unaltered after adjustment for potential confounding factors. There was no association with fracture overall (odds ratio [OR] = 1.01 [CI, 0.95 to 1.08]) in either unadjusted or adjusted analyses, but there was a non-significant trend for association with vertebral fracture and a nominally significant association with incident vertebral fractures in females (OR = 1.33 [CI, 1.00 to 1.77]) that was independent of BMD, and unaltered in adjusted analyses. Conclusions Allowing for the inevitable heterogeneity between participating teams, this study—which to our knowledge is the largest ever performed in the field of osteoporosis genetics for a single gene—demonstrates that the COLIA1 Sp1 polymorphism is associated with reduced BMD and could predispose to incident vertebral fractures in women, independent of BMD. The associations we observed were modest however, demonstrating the importance of conducting studies that are adequately powered to detect and quantify the effects of common genetic variants on complex diseases.

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

IntroductionCertain 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.MethodsWe 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.ResultsNo significant breast cancer associations were detected with any individual or combination of tag SNPs.ConclusionIt is unlikely that there are any other common variants in these genes conferring measurably increased risks of breast cancer in our study population.

Large-scale analysis of association between polymorphisms in the transforming growth factor beta 1 gene (TGFB1) and osteoporosis : The GENOMOS study

INTRODUCTION The TGFB1 gene which encodes transforming growth factor beta 1, is a strong candidate for susceptibility to osteoporosis and several studies have reported associations between bone mineral density (BMD), osteoporotic fractures and polymorphisms of TGFB1, although these studies have yielded conflicting results. METHODS We investigated associations between TGFB1 polymorphisms and BMD and fracture in the GENOMOS study: a prospective multicenter study involving 10 European research studies including a total of 28,924 participants. Genotyping was conducted for known TGFB1 polymorphisms at the following sites: G-1639-A (G-800-A, rs1800468), C-1348-T (C-509-T, rs1800469), T29-C (Leu10Pro, rs1982073), G74-C (Arg25Pro, rs1800471) and C788-T (Thr263Ile, rs1800472). These polymorphisms were genotyped prospectively and methodology was standardized across research centers. Genotypes and haplotypes were related to BMD at the lumbar sine and femoral neck and fractures. RESULTS There were no significant differences in either women or men at either skeletal site for any of the examined polymorphisms with the possible exception of a weak association with reduced BMD (-12 mg/cm2) in men with the T-1348 allele (p<0.05). None of the haplotypes was associated with BMD and none of the polymorphisms or haplotypes significantly affected overall risk of fractures, however, the odds ratio for incident vertebral fracture in carriers of the rare T788 allele was 1.64 (95% CI: 1.09-2.64), p<0.05. CONCLUSIONS This study indicates that polymorphic variation in the TGFB1 gene does not play a major role in regulating BMD or susceptibility to fractures. The weak associations we observed between the C-1348-T and lumbar spine BMD in men and between C788-T and risk of incident vertebral fractures are of interest but could be chance findings and will need replication in future studies.

Association of ESR1 gene tagging SNPs with breast cancer risk

We have conducted a three-stage, comprehensive single nucleotide polymorphism (SNP)-tagging association study of ESR1 gene variants (SNPs) in more than 55,000 breast cancer cases and controls from studies within the Breast Cancer Association Consortium (BCAC). No large risks or highly significant associations were revealed. SNP rs3020314, tagging a region of ESR1 intron 4, is associated with an increase in breast cancer susceptibility with a dominant mode of action in European populations. Carriers of the c-allele have an odds ratio (OR) of 1.05 [95% Confidence Intervals (CI) 1.02-1.09] relative to t-allele homozygotes, P = 0.004. There is significant heterogeneity between studies, P = 0.002. The increased risk appears largely confined to oestrogen receptor-positive tumour risk. The region tagged by SNP rs3020314 contains sequence that is more highly conserved across mammalian species than the rest of intron 4, and it may subtly alter the ratio of two mRNA splice forms.

Phytoestrogen Exposure, Polymorphisms in COMT, CYP19, ESR1, and SHBG Genes, and Their Associations With Prostate Cancer Risk

Abstract: Prospective phytoestrogen exposure was assessed using both biomarkers and estimates of intake in 89 British men recruited into the Norfolk arm of the European Prospective Investigation into Cancer and Nutrition study, men who subsequently developed prostate cancer. Results were compared with those from 178 healthy men matched by age and date of recruitment. Levels of seven phytoestrogens (daidzein, genistein, glycitein, O-desmethylangolensin, equol, enterodiol, and enterolactone) were measured in spot urine and serum samples. Five single-nucleotide polymorphisms in COMT, CYP19, ESR1, and SHBG genes were genotyped. Urinary levels of all phytoestrogens correlated strongly with serum levels. Correlation coefficients ranged from 0.63 (glycitein) to 0.88 (daidzein) (P < 0.001). Urinary and serum levels correlated significantly with isoflavone intake assessed from food diaries (R = 0.15–0.20; P < 0.05) but not with that from a food-frequency questionnaire. Odds ratios for phytoestrogen exposure, as assessed using the four methods, were not significantly associated with prostate cancer risk (P = 0.15–0.94). Men with the CC genotype for the ESRI PvuII polymorphism had significantly higher risk for prostate cancer compared with men with the TT genotype [adjusted odds ratio = 4.65 (1.60–13.49); P = 0.005]. Our results utilizing a combined prospective exposure provide no evidence that phytoestrogens alter prostate cancer risk in British men, whereas the C allele for the PvuII polymorphism may be associated with increased risk.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus

Background Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. Aim To identify CNVs associated with osteoporotic bone fracture risk. Method We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. Results A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p=8.69×10−5). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p=0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. Conclusions These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

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

BackgroundIn 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.ImplementationWe 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: TaqMans Primer-By-Design FileBuilder, Sequenoms iPLEX or SNPstreams 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.ConclusionThe 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.

Association of gene variants in the TGF-beta signalling pathways with invasive breast cancer risk

Many studies of normal cells in vitro, of transgenic mouse models and of somatic mutations in human cancers have provided evidence that the cytokine transforming growth factor beta (TGFβ) acts as a suppressor of primary tumour initiation. However, studies of transformed cells in vitro and of mouse models have also implicated TGFβ as a promoter of the later stages of tumour development. A hypothesis that has been proposed to account for this dual action is that TGFβ acts as a tumour suppressor through the ubiquitous ALK5 receptor signalling via SMAD2 and SMAD3 to inhibit proliferation of primary tumour cells, but acts subsequently through the endothelial-specific ALK1 receptor via SMAD1 and SMAD5 to promote angiogenesis, which is required for tumour progression [1]. In a recent meta-analysis we showed that a single nucleotide polymorphism (SNP) generating a leucine to proline substitution in the signal peptide of the TGFB1 protein is associated with an increase in the risk of invasive breast cancer (OR per additional proline allele = 1.08 (95% CI = 1.04 to 1.11), Ptrend = 2.8 × 10-5) [2]. We have also reported that this SNP increases the amount of TGFB1 protein secreted in vitro by threefold [3]. These data suggest that higher levels of TGFB1 may promote the invasive breast tumour phenotype. To determine the effect of host TGFB1 levels in an in vivo model, Tgfb1+/- and Tgfb1+/+ mice have been compared in which the mice carry one or two TGFB1 alleles with the ancestral SNP form encoding proline in the signal peptide. These studies have revealed major effects of TGFB1 in controlling the site of metastatic seeding and the number of metastases that develop. The TGFB1 SNP association with breast cancer suggested that other genes in the TGFβ signalling pathways might be associated with altered risk. We have conducted association studies with SNPs in 16 further genes encoding proteins directly implicated in TGFβ signalling. LTBP1, LTBP2, LTBP4, TGFB1, TGFB2, TGFB3, ALK1, ALK5, TGFBR2, Endoglin, SMAD1, SMAD2, SMAD3, SMAD4, SMAD5, SMAD6 and SMAD7 were analysed. A comprehensive SNP tagging approach was used to select variants for genotyping in a staged study design using up to 6,900 cases and 6,900 controls, all collected from the East Anglia region of the UK (>98% of northwestern European ancestry). From 1,254 common SNPs (minor allele frequency >0.05) in these genes identified from the International HapMap project data, we defined and genotyped 354 tagging SNPs in the East Anglia cases and controls. Statistically significant associations were followed up by genotyping in a Polish set of 2,215 cases and 2,374 controls. Meta-analysis of these results identified associations with cancer susceptibility of a variant in ALK5 (OR per additional rare G allele = 0.88 (95% CI = 0.81 to 0.95), Ptrend = 0.001) and in TGFBR2 (OR per additional rare G allele = 0.96 (95% CI = 0.92 to 1.00), Ptrend = 0.039). Data from two genome-wide studies have been examined to search further for associations in these genes. The haplotype risks and interactions between two or more loci have been investigated and survival analyses have been conducted. From this comprehensive study we have identified tagging SNPs in two TGFβ receptor genes that are significantly associated with risk of invasive breast cancer. Identification of the causative SNPs within the LD blocks tagged by the two SNPs and their effects on signalling via the ALK1 and ALK5 pathways remain to be determined.