A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer

We conducted a genome-wide association study (GWAS) of breast cancer by genotyping 528,173 SNPs in 1,145 postmenopausal women of European ancestry with invasive breast cancer and 1,142 controls. We identified four SNPs in intron 2 of FGFR2 (which encodes a receptor tyrosine kinase and is amplified or overexpressed in some breast cancers) that were highly associated with breast cancer and confirmed this association in 1,776 affected individuals and 2,072 controls from three additional studies. Across the four studies, the association with all four SNPs was highly statistically significant (P(trend) for the most strongly associated SNP (rs1219648) = 1.1 x 10(-10); population attributable risk = 16%). Four SNPs at other loci most strongly associated with breast cancer in the initial GWAS were not associated in the replication studies. Our summary results from the GWAS are available online in a form that should speed the identification of additional risk loci.     

Association analyses identify six new psoriasis susceptibility loci in the Chinese population

We extended our previous genome-wide association study for psoriasis with a multistage replication study including 8,312 individuals with psoriasis (cases) and 12,919 controls from China as well as 3,293 cases and 4,188 controls from Germany and the United States and 254 nuclear families from the United States. We identified six new susceptibility loci associated with psoriasis in the Chinese study containing the candidate genes ERAP1, PTTG1, CSMD1, GJB2, SERPINB8 and ZNF816A (combined P < 5 × 10??) and replicated one locus, 5q33.1 (TNIP1-ANXA6), previously reported (combined P = 3.8 × 10?21) in the European studies. Two of these loci showed evidence for association in the German study at ZNF816A and GJB2 with P = 3.6 × 10?3 and P = 7.9 × 10?3, respectively. ERAP1 and ZNF816A were associated with type 1 (early onset) psoriasis in the Chinese Han population (test for heterogeneity P = 6.5 × 10?3 and P = 1.5 × 10?3, respectively). Comparisons with the results of previous GWAS of psoriasis highlight the heterogeneity of disease susceptibility between the Chinese and European populations. Our study identifies new genetic susceptibility factors and suggests new biological pathways in psoriasis.     

In situ analyses of genome instability in breast cancer

Transition through telomere crisis is thought to be a crucial event in the development of most breast carcinomas. Our goal in this study was to determine where this occurs in the context of histologically defined breast cancer progression. To this end, we assessed genome instability (using fluorescence in situ hybridization) and other features associated with telomere crisis in normal ductal epithelium, usual ductal hyperplasia, ductal carcinoma in situ and invasive cancer. We modeled this process in vitro by measuring these same features in human mammary epithelial cell cultures during ZNF217-mediated transition through telomere crisis and immortalization. Taken together, the data suggest that transition through telomere crisis and immortalization in breast cancer occurs during progression from usual ductal hyperplasia to ductal carcinoma in situ.     

Genome-wide association analysis identifies three new breast cancer susceptibility loci

Breast cancer is the most common cancer among women. To date, 22 common breast cancer susceptibility loci have been identified accounting for ～8% of the heritability of the disease. We attempted to replicate 72 promising associations from two independent genome-wide association studies (GWAS) in ～70,000 cases and ～68,000 controls from 41 case-control studies and 9 breast cancer GWAS. We identified three new breast cancer risk loci at 12p11 (rs10771399; P = 2.7 × 10(-35)), 12q24 (rs1292011; P = 4.3 × 10(-19)) and 21q21 (rs2823093; P = 1.1 × 10(-12)). rs10771399 was associated with similar relative risks for both estrogen receptor (ER)-negative and ER-positive breast cancer, whereas the other two loci were associated only with ER-positive disease. Two of the loci lie in regions that contain strong plausible candidate genes: PTHLH (12p11) has a crucial role in mammary gland development and the establishment of bone metastasis in breast cancer, and NRIP1 (21q21) encodes an ER cofactor and has a role in the regulation of breast cancer cell growth.     

A common variant at the TERT-CLPTM1L locus is associated with estrogen receptor-negative breast cancer

Estrogen receptor (ER)-negative breast cancer shows a higher incidence in women of African ancestry compared to women of European ancestry. In search of common risk alleles for ER-negative breast cancer, we combined genome-wide association study (GWAS) data from women of African ancestry (1,004 ER-negative cases and 2,745 controls) and European ancestry (1,718 ER-negative cases and 3,670 controls), with replication testing conducted in an additional 2,292 ER-negative cases and 16,901 controls of European ancestry. We identified a common risk variant for ER-negative breast cancer at the TERT-CLPTM1L locus on chromosome 5p15 (rs10069690: per-allele odds ratio (OR) = 1.18 per allele, P = 1.0 × 10(-10)). The variant was also significantly associated with triple-negative (ER-negative, progesterone receptor (PR)-negative and human epidermal growth factor-2 (HER2)-negative) breast cancer (OR = 1.25, P = 1.1 × 10(-9)), particularly in younger women (<50 years of age) (OR = 1.48, P = 1.9 × 10(-9)). Our results identify a genetic locus associated with estrogen receptor negative breast cancer subtypes in multiple populations.     

Identification of loci associated with schizophrenia by genome-wide association and follow-up

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10(-5) in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 x 10(-4)), and the overall pattern of replication was unlikely to occur by chance (P = 9 x 10(-8)). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 x 10(-7)) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 x 10(-9)).     

Common variants on chromosome 5p12 confer susceptibility to estrogen receptor-positive breast cancer

We carried out a genome-wide association study of breast cancer predisposition with replication and refinement studies involving 6,145 cases and 33,016 controls and identified two SNPs (rs4415084 and rs10941679) on 5p12 that confer risk, preferentially for estrogen receptor (ER)-positive tumors (OR = 1.27, P = 2.5 x 10(-12) for rs10941679). The nearest gene, MRPS30, was previously implicated in apoptosis, ER-positive tumors and favorable prognosis. A recently reported signal in FGFR2 was also found to associate specifically with ER-positive breast cancer.     

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; P(trend) = 1.1 x 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; P(trend) = 2.8 x 10(-5)). These results demonstrate that common breast cancer susceptibility alleles with small effects on risk can be identified, given sufficiently powerful studies.     

Quantitative mapping of amplicon structure by array CGH identifies CYP24 as a candidate oncogene

We show here that quantitative measurement of DNA copy number across amplified regions using array comparative genomic hybridization (CGH) may facilitate oncogene identification by providing precise information on the locations of both amplicon boundaries and amplification maxima. Using this analytical capability, we resolved two regions of amplification within an approximately 2-Mb region of recurrent aberration at 20q13.2 in breast cancer. The putative oncogene ZNF217 (ref. 5) mapped to one peak, and CYP24 (encoding vitamin D 24 hydroxylase), whose overexpression is likely to lead to abrogation of growth control mediated by vitamin D, mapped to the other.     

Germline mutations in RAD51D confer susceptibility to ovarian cancer

Recently, RAD51C mutations were identified in families with breast and ovarian cancer. This observation prompted us to investigate the role of RAD51D in cancer susceptibility. We identified eight inactivating RAD51D mutations in unrelated individuals from 911 breast-ovarian cancer families compared with one inactivating mutation identified in 1,060 controls (P = 0.01). The association found here was principally with ovarian cancer, with three mutations identified in the 59 pedigrees with three or more individuals with ovarian cancer (P = 0.0005). The relative risk of ovarian cancer for RAD51D mutation carriers was estimated to be 6.30 (95% CI 2.86-13.85, P = 4.8 × 10(-6)). By contrast, we estimated the relative risk of breast cancer to be 1.32 (95% CI 0.59-2.96, P = 0.50). These data indicate that RAD51D mutation testing may have clinical utility in individuals with ovarian cancer and their families. Moreover, we show that cells deficient in RAD51D are sensitive to treatment with a PARP inhibitor, suggesting a possible therapeutic approach for cancers arising in RAD51D mutation carriers.     

A germline variant in the TP53 polyadenylation signal confers cancer susceptibility

To identify new risk variants for cutaneous basal cell carcinoma, we performed a genome-wide association study of 16 million SNPs identified through whole-genome sequencing of 457 Icelanders. We imputed genotypes for 41,675 Illumina SNP chip-typed Icelanders and their relatives. In the discovery phase, the strongest signal came from rs78378222[C] (odds ratio (OR) = 2.36, P = 5.2 × 10(-17)), which has a frequency of 0.0192 in the Icelandic population. We then confirmed this association in non-Icelandic samples (OR = 1.75, P = 0.0060; overall OR = 2.16, P = 2.2 × 10(-20)). rs78378222 is in the 3' untranslated region of TP53 and changes the AATAAA polyadenylation signal to AATACA, resulting in impaired 3'-end processing of TP53 mRNA. Investigation of other tumor types identified associations of this SNP with prostate cancer (OR = 1.44, P = 2.4 × 10(-6)), glioma (OR = 2.35, P = 1.0 × 10(-5)) and colorectal adenoma (OR = 1.39, P = 1.6 × 10(-4)). However, we observed no effect for breast cancer, a common Li-Fraumeni syndrome tumor (OR = 1.06, P = 0.57, 95% confidence interval 0.88-1.27).     

Sequence variant on 8q24 confers susceptibility to urinary bladder cancer

We conducted a genome-wide SNP association study on 1,803 urinary bladder cancer (UBC) cases and 34,336 controls from Iceland and The Netherlands and follow up studies in seven additional case-control groups (2,165 cases and 3,800 controls). The strongest association was observed with allele T of rs9642880 on chromosome 8q24, 30 kb upstream of MYC (allele-specific odds ratio (OR) = 1.22; P = 9.34 x 10(-12)). Approximately 20% of individuals of European ancestry are homozygous for rs9642880[T], and their estimated risk of developing UBC is 1.49 times that of noncarriers. No association was observed between UBC and the four 8q24 variants previously associated with prostate, colorectal and breast cancers, nor did rs9642880 associate with any of these three cancers. A weaker signal, but nonetheless of genome-wide significance, was captured by rs710521[A] located near TP63 on chromosome 3q28 (allele-specific OR = 1.19; P = 1. 15 x 10(-7)).     

Estrogen receptor alpha (ESR1) gene amplification is frequent in breast cancer

Using an Affymetrix 10K SNP array to screen for gene copy number changes in breast cancer, we detected a single-gene amplification of the ESR1 gene, which encodes estrogen receptor alpha, at 6q25. A subsequent tissue microarray analysis of more than 2,000 clinical breast cancer samples showed ESR1 amplification in 20.6% of breast cancers. Ninety-nine percent of tumors with ESR1 amplification showed estrogen receptor protein overexpression, compared with 66.6% cancers without ESR1 amplification (P < 0.0001). In 175 women who had received adjuvant tamoxifen monotherapy, survival was significantly longer for women with cancer with ESR1 amplification than for women with estrogen receptor-expressing cancers without ESR1 amplification (P = 0.023). Notably, we also found ESR1 amplification in benign and precancerous breast diseases, suggesting that ESR1 amplification may be a common mechanism in proliferative breast disease and a very early genetic alteration in a large subset of breast cancers.     

PALB2, which encodes a BRCA2-interacting protein, is a breast cancer susceptibility gene

PALB2 interacts with BRCA2, and biallelic mutations in PALB2 (also known as FANCN), similar to biallelic BRCA2 mutations, cause Fanconi anemia. We identified monoallelic truncating PALB2 mutations in 10/923 individuals with familial breast cancer compared with 0/1,084 controls (P = 0.0004) and show that such mutations confer a 2.3-fold higher risk of breast cancer (95% confidence interval (c.i.) = 1.4-3.9, P = 0.0025). The results show that PALB2 is a breast cancer susceptibility gene and further demonstrate the close relationship of the Fanconi anemia-DNA repair pathway and breast cancer predisposition.     

Mutations in BRIP1 confer high risk of ovarian cancer

Ovarian cancer causes more deaths than any other gynecologic malignancy in developed countries. Sixteen million sequence variants, identified through whole-genome sequencing of 457 Icelanders, were imputed to 41,675 Icelanders genotyped using SNP chips, as well as to their relatives. Sequence variants were tested for association with ovarian cancer (N of affected individuals = 656). We discovered a rare (0.41% allelic frequency) frameshift mutation, c.2040_2041insTT, in the BRIP1 (FANCJ) gene that confers an increase in ovarian cancer risk (odds ratio (OR) = 8.13, P = 2.8 × 10(-14)). The mutation was also associated with increased risk of cancer in general and reduced lifespan by 3.6 years. In a Spanish population, another frameshift mutation in BRIP1, c.1702_1703del, was seen in 2 out of 144 subjects with ovarian cancer and 1 out of 1,780 control subjects (P = 0.016). This allele was also associated with breast cancer (seen in 6/927 cases; P = 0.0079). Ovarian tumors from heterozygous carriers of the Icelandic mutation show loss of the wild-type allele, indicating that BRIP1 behaves like a classical tumor suppressor gene in ovarian cancer.     

NAD(P)H:quinone oxidoreductase 1 NQO1*2 genotype (P187S) is a strong prognostic and predictive factor in breast cancer

NQO1 guards against oxidative stress and carcinogenesis and stabilizes p53. We find that a homozygous common missense variant (NQO1(*)2, rs1800566(T), NM_000903.2:c.558C>T) that disables NQO1 strongly predicts poor survival among two independent series of women with breast cancer (P = 0.002, N = 1,005; P = 0.005, N = 1,162), an effect particularly evident after anthracycline-based adjuvant chemotherapy with epirubicin (P = 7.52 x 10(-6)) and in p53-aberrant tumors (P = 6.15 x 10(-5)). Survival after metastasis was reduced among NQO1(*)2 homozygotes, further implicating NQO1 deficiency in cancer progression and treatment resistance. Consistently, response to epirubicin was impaired in NQO1(*)2-homozygous breast carcinoma cells in vitro, reflecting both p53-linked and p53-independent roles of NQO1. We propose a model of defective anthracycline response in NQO1-deficient breast tumors, along with increased genomic instability promoted by elevated reactive oxygen species (ROS), and suggest that the NQO1 genotype is a prognostic and predictive marker for breast cancer.     

Genome-wide association study identifies a common variant associated with risk of endometrial cancer

Endometrial cancer is the most common malignancy of the female genital tract in developed countries. To identify genetic variants associated with endometrial cancer risk, we performed a genome-wide association study involving 1,265 individuals with endometrial cancer (cases) from Australia and the UK and 5,190 controls from the Wellcome Trust Case Control Consortium. We compared genotype frequencies in cases and controls for 519,655 SNPs. Forty seven SNPs that showed evidence of association with endometrial cancer in stage 1 were genotyped in 3,957 additional cases and 6,886 controls. We identified an endometrial cancer susceptibility locus close to HNF1B at 17q12 (rs4430796, P = 7.1 × 10(-10)) that is also associated with risk of prostate cancer and is inversely associated with risk of type 2 diabetes.     

Common variants at 11q12, 10q26 and 3p11.2 are associated with prostate cancer susceptibility in Japanese

We have previously reported multiple loci associated with prostate cancer susceptibility in a Japanese population using a genome-wide association study (GWAS). To identify additional prostate cancer susceptibility loci, we genotyped nine SNPs that were nominally associated with prostate cancer (P < 1 × 10(-4)) in our previous GWAS in three independent studies of prostate cancer in Japanese men (2,557 individuals with prostate cancer (cases) and 3,003 controls). In a meta-analysis of our previous GWAS and the replication studies, which included a total of 7,141 prostate cancer cases and 11,804 controls from a single ancestry group, three new loci reached genome-wide significance on chromosomes 11q12 (rs1938781; P = 1.10 × 10(-10); FAM111A-FAM111B), 10q26 (rs2252004; P = 1.98 × 10(-8)) and 3p11.2 (rs2055109; P = 3.94 × 10(-8)). We also found suggestive evidence of association at a previously reported prostate cancer susceptibility locus at 2p11 (rs2028898; P = 1.08 × 10(-7)). The identification of three new susceptibility loci should provide additional insight into the pathogenesis of prostate cancer and emphasizes the importance of conducting GWAS in diverse populations.     

Discovery of common variants associated with low TSH levels and thyroid cancer risk

To search for sequence variants conferring risk of nonmedullary thyroid cancer, we focused our analysis on 22 SNPs with a P < 5 × 10(-8) in a genome-wide association study on levels of thyroid stimulating hormone (TSH) in 27,758 Icelanders. Of those, rs965513 has previously been shown to associate with thyroid cancer. The remaining 21 SNPs were genotyped in 561 Icelandic individuals with thyroid cancer (cases) and up to 40,013 controls. Variants suggestively associated with thyroid cancer (P < 0.05) were genotyped in an additional 595 non-Icelandic cases and 2,604 controls. After combining the results, three variants were shown to associate with thyroid cancer: rs966423 on 2q35 (OR = 1.34; P(combined) = 1.3 × 10(-9)), rs2439302 on 8p12 (OR = 1.36; P(combined) = 2.0 × 10(-9)) and rs116909374 on 14q13.3 (OR = 2.09; P(combined) = 4.6 × 10(-11)), a region previously reported to contain an uncorrelated variant conferring risk of thyroid cancer. A strong association (P = 9.1 × 10(-91)) was observed between rs2439302 on 8p12 and expression of NRG1, which encodes the signaling protein neuregulin 1, in blood.     

Genome-wide association study of prostate cancer identifies a second risk locus at 8q24

Recently, common variants on human chromosome 8q24 were found to be associated with prostate cancer risk. While conducting a genome-wide association study in the Cancer Genetic Markers of Susceptibility project with 550,000 SNPs in a nested case-control study (1,172 cases and 1,157 controls of European origin), we identified a new association at 8q24 with an independent effect on prostate cancer susceptibility. The most significant signal is 70 kb centromeric to the previously reported SNP, rs1447295, but shows little evidence of linkage disequilibrium with it. A combined analysis with four additional studies (total: 4,296 cases and 4,299 controls) confirms association with prostate cancer for rs6983267 in the centromeric locus (P = 9.42 x 10(-13); heterozygote odds ratio (OR): 1.26, 95% confidence interval (c.i.): 1.13-1.41; homozygote OR: 1.58, 95% c.i.: 1.40-1.78). Each SNP remained significant in a joint analysis after adjusting for the other (rs1447295 P = 1.41 x 10(-11); rs6983267 P = 6.62 x 10(-10)). These observations, combined with compelling evidence for a recombination hotspot between the two markers, indicate the presence of at least two independent loci within 8q24 that contribute to prostate cancer in men of European ancestry. We estimate that the population attributable risk of the new locus, marked by rs6983267, is higher than the locus marked by rs1447295 (21% versus 9%).