Newborn Screening for SCID and Other Severe Primary Immunodeficiency in the Polish-German Transborder Area: Experience From the First 14 Months of Collaboration

In 2017, in the Polish-German transborder area of West Pomerania, Mecklenburg-Western Pomerania, and Brandenburg, in collaboration with two centers in Warsaw, a partnership in the field of newborn screening (NBS) for severe primary immunodeficiency diseases (PID), mainly severe combined immunodeficiency (SCID), was initiated. SCID, but also some other severe PID, is a group of disorders characterized by the absence of T and/or B and NK cells. Affected infants are susceptible to life-threatening infections, but early detection gives a chance for effective treatment. The prevalence of SCID in the Polish and German populations is unknown but can be comparable to other countries (1:50,000–100,000). SCID NBS tests are based on real-time polymerase chain reaction (qPCR) and the measurement of a number of T cell receptor excision circles (TREC), kappa-deleting recombination excision circles (KREC), and beta-actin (ACTB) as a quality marker of DNA. This method can also be effective in NBS for other severe PID with T- and/or B-cell lymphopenia, including combined immunodeficiency (CID) or agammaglobulinemia. During the 14 months of collaboration, 44,287 newborns were screened according to the ImmunoIVD protocol. Within 65 positive samples, seven were classified to immediate recall and 58 requested a second sample. Examination of the 58 second samples resulted in recalling one newborn. Confirmatory tests included immunophenotyping of lymphocyte subsets with extension to TCR repertoire, lymphoproliferation tests, radiosensitivity tests, maternal engraftment assays, and molecular tests. Final diagnosis included: one case of T-B^lowNK+ SCID, one case of atypical T^low B^lowNK+ CID, one case of autosomal recessive agammaglobulinemia, and one case of Nijmegen breakage syndrome. Among four other positive results, three infants presented with T- and/or B-cell lymphopenia due to either the mother's immunosuppression, prematurity, or unknown reasons, which resolved or almost normalized in the first months of life. One newborn was classified as truly false positive. The overall positive predictive value (PPV) for the diagnosis of severe PID was 50.0%. This is the first population screening study that allowed identification of newborns with T and/or B immunodeficiency in Central and Eastern Europe.

Polymorphisms in MMP-1, MMP-2, MMP-7, MMP-13 and MT2A do not contribute to breast, lung and colon cancer risk in polish population

Background

Matrix metalloproteinases (MMPs) and metallothioneins (MTs) are Zinc-related proteins which are involved in processes crucial for carcinogenesis such as angiogenesis, proliferation and apoptosis. Several single nucleotide polymorphisms (SNPs) in MMPs and MTs that affect genes expression have been associated with cancer risk, including breast, lung and colon.

Methods

The study group consisted of 648 unselected patients (299 with breast cancer, 199 with lung cancer, 150 with colon cancer) and 648 unaffected individuals. Five SNPs, rs1799750 in MMP-1, rs243865 in MMP-2, rs11568818 in MMP-7, rs2252070 in MMP-13 and rs28366003 in MT2A were genotyped and serum zinc (Zn) level was measured. The cancer risk was calculated using multivariable logistic regression with respect to Zn.

Results

None of the 5 tested polymorphisms showed a correlation with cancer risk in studied groups, although for MMP-2, MMP-7 and MT2A non-significant differences in genotypes frequencies among cases and controls were observed.

Conclusions

Analyses of polymorphisms, rs1799750 in MMP-1, rs243865 in MMP-2, rs11568818 in MMP-7, rs2252070 in MMP-13 and rs28366003 in MT2A in relation to serum Zn level did not show significant association with breast, lung and colon cancer risk among polish patients. Further studies are needed to verify this observation.

The FANCM:p.Arg658* truncating variant is associated with risk of triple-negative breast cancer

Breast cancer is a common disease partially caused by genetic risk factors. Germline pathogenic variants in DNA repair genes BRCA1, BRCA2, PALB2, ATM, and CHEK2 are associated with breast cancer risk. FANCM, which encodes for a DNA translocase, has been proposed as a breast cancer predisposition gene, with greater effects for the ER-negative and triple-negative breast cancer (TNBC) subtypes. We tested the three recurrent protein-truncating variants FANCM:p.Arg658*, p.Gln1701*, and p.Arg1931* for association with breast cancer risk in 67,112 cases, 53,766 controls, and 26,662 carriers of pathogenic variants of BRCA1 or BRCA2. These three variants were also studied functionally by measuring survival and chromosome fragility in FANCM -/- patient-derived immortalized fibroblasts treated with diepoxybutane or olaparib. We observed that FANCM:p.Arg658* was associated with increased risk of ER-negative disease and TNBC (OR = 2.44, P = 0.034 and OR = 3.79; P = 0.009, respectively). In a country-restricted analysis, we confirmed the associations detected for FANCM:p.Arg658* and found that also FANCM:p.Arg1931* was associated with ER-negative breast cancer risk (OR = 1.96; P = 0.006). The functional results indicated that all three variants were deleterious affecting cell survival and chromosome stability with FANCM:p.Arg658* causing more severe phenotypes. In conclusion, we confirmed that the two rare FANCM deleterious variants p.Arg658* and p.Arg1931* are risk factors for ER-negative and TNBC subtypes. Overall our data suggest that the effect of truncating variants on breast cancer risk may depend on their position in the gene. Cell sensitivity to olaparib exposure, identifies a possible therapeutic option to treat FANCM-associated tumors.

Iron levels, genes involved in iron metabolism and antioxidative processes and lung cancer incidence

Background

Lung cancer is the most common adult malignancy accounting for the largest proportion of cancer related deaths. Iron (Fe) is an essential trace element and is a component of several major metabolic pathways playing an important role in many physiological processes. In this study we evaluated the association between Fe concentration in serum, iron metabolism parameters and genetic variaton in 7 genes involved in iron metabolism and anti-oxidative processes with the incidence of lung cancer in Poland.

Materials and methods

The study included 200 lung cancer patients and 200 matched healthy control subjects. We analyzed serum iron concentration and iron metabolism parameters (TIBC, UIBC, serum ferritin and transferrin saturation), and genotyped seven variants in seven genes: HFE, TFR1, HAMP, TF, SOD2, CAT and GPX1.

Results

Lung cancer patients compared to their matched controls had significantly higher mean serum iron level (p = 0.01), ferritin level (p = 0.007) and TIBC (p = 0.006). Analysis revealed that higher concentration of iron and ferritin (IVth quartile) compared to the lower concentration (Ist quartile) was associated with over 2-fold increased lung cancer incidence. We also found that higher transferrin saturation (p = 0.01) and lower TIBC (p<0.01) are associated with better survival of lung cancer patients. The analysis of polymorphisms in iron related genes did not reveal a significant difference between lung cancer patients and controls. However, rs10421768 in HAMP showed a borderline statistically significant correlation with lung cancer risk (OR = 2.83, p = 0.05).

Conclusions

The results of this case control study indicate that higher body iron represented by higher Fe and ferritin levels may be associated with lung cancer incidence. Rs10421768 in HAMP may be associated with about 3-times higher lung cancer risk. Higher Fe body content may be associated with better survival of lung cancer patients.

Association of zinc level and polymorphism in MMP-7 gene with prostate cancer in Polish population

INTRODUCTION

Prostate cancer is one of the most commonly diagnosed malignancies among men in Western populations. Evidence reported in the literature suggests that zinc may be related to prostate cancer. In this study we evaluated the association of serum zinc levels and polymorphisms in genes encoding zinc-dependent proteins with prostate cancer in Poland.

METHODS

The study group consisted of 197 men affected with prostate cancer and 197 healthy men. Serum zinc levels were measured and 5 single nucleotide polymorphisms in MMP-1, MMP-2, MMP-7, MMP-13, MT2A genes were genotyped.

RESULTS

The mean serum zinc level was higher in prostate cancer patients than in healthy controls (898.9±12.01 μg/l vs. 856.6±13.05 μg/l, p<0.01). When compared in quartiles a significant association of higher zinc concentration with the incidence of prostate cancer was observed. The highest OR (OR = 4.41, 95%CI 2.07-9.37, p<0.01) was observed in 3rd quartile (>853.0-973.9 μg/l). Among five analyzed genetic variants, rs11568818 in MMP-7 appeared to be correlated with 2-fold increased prostate cancer risk (OR = 2.39, 95% CI = 1.19-4.82, p = 0.015).

CONCLUSION

Our results suggest a significant correlation of higher serum zinc levels with the diagnosis of prostate cancer. The polymorphism rs11568818 in MMP-7 gene was also associated with an increased prostate cancer risk in Poland.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

Serum folate concentration and the incidence of lung cancer

Background

Lung cancer is a leading cause of cancer-related mortality globally. Folate helps to maintain DNA integrity and to regulate gene expression. Serum folate levels may affect the risk of several cancers, including lung cancer. In this study we evaluated the association between serum folate concentration and variations in genes involved in folate metabolism with lung cancer incidence in Poland.

Methods

The study included 366 lung cancer patients and 366 control subjects. We measured serum folate concentration and genotyped six variants in MTHFR, MTR and MTRR genes. The odds ratios of being diagnosed with lung cancer were calculated using conditional univariable and multivariable logistic regression with respect to folate level and genotypes.

Results

The mean serum folate level was lower in lung cancer cases than in control group (20.07 nmol/l vs. 22.52 nmol/l, p = 0.002). The odds ratio for lung cancer declined with increasing serum content of the folate. The folate concentration of >25.71 nmol/l (IVth quartile) in comparison to <15.92 nmol/l (Ist quartile) was associated with an odds ratio of 0.61 (95%CI 0.40–0.95, p = 0.03). The analysis of variations in MTHFR, MTR and MTRR genes did not reveal any significant difference between lung cancer cases and controls in univariable and multivariable analyses.

Conclusion

In this case-control study, lower serum folate concentrations were associated with a higher risk of lung cancer diagnosis. Although previous findings have been somewhat mixed, our results add to the evidence that circulating folate levels may be an indicator of lung cancer risk.

Serum 25(OH)D concentration, common variants of the VDR gene and lung cancer occurrence

The first aim of our study was to examine the association between common variants in VDR [rs2228570 (FokI), rs1544410 (BsmI), rs7975232 (ApaI), rs731236 (TaqI) and rs11568820 (Cdx2)] and lung cancer risk in the Polish population. Genotyping and statistical analysis which included Chi-square test with Yates correction and haplotype frequency analysis were performed on a series of 840 consecutively collected lung cancer patients and 920 healthy controls. The second aim was to evaluate the link between serum 25(OH)D concentration and the number of lung cancers in a subgroup of 200 patients. A separate control group that consisted of 400 matched (by age, sex, smoking habits and the season of blood collection) healthy individuals was used to avoid posterior adjustment on the matched variables. Statistical analysis with the use of Chi-square test with Yates was performed. We found no statistically significant difference in the distribution of the allels of studied VDR variants among cases and controls. A statistically significant over-representation of VDR haplotypes: rs731236_A + rs1544410_T [odds ratio (OR) = 2.43, 95% confidence interval (CI) = 1.11-5.32, p < 0.001], rs731236_G + rs1544410_T (OR = 1.54, 95% CI = 1.31-1.81, p < 0.001) and rs731236_G + rs1544410_C (OR = 0.04, 95% CI = 0.03-0.07, p < 0.001) was detected. We found a tendency toward an increased number of lung cancers among individuals with low serum levels of 25(OH)D. To answer the question, whether VDR can be regarded as lung cancer susceptibility gene and low 25(OH)D serum levels is associated with lung cancer occurrences, additional, multicenter study needs to be performed.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Inheritance of deleterious mutations at both BRCA1 and BRCA2 in an international sample of 32,295 women

BACKGROUND

Most BRCA1 or BRCA2 mutation carriers have inherited a single (heterozygous) mutation. Transheterozygotes (TH) who have inherited deleterious mutations in both BRCA1 and BRCA2 are rare, and the consequences of transheterozygosity are poorly understood.

METHODS

From 32,295 female BRCA1/2 mutation carriers, we identified 93 TH (0.3 %). "Cases" were defined as TH, and "controls" were single mutations at BRCA1 (SH1) or BRCA2 (SH2). Matched SH1 "controls" carried a BRCA1 mutation found in the TH "case". Matched SH2 "controls" carried a BRCA2 mutation found in the TH "case". After matching the TH carriers with SH1 or SH2, 91 TH were matched to 9316 SH1, and 89 TH were matched to 3370 SH2.

RESULTS

The majority of TH (45.2 %) involved the three common Jewish mutations. TH were more likely than SH1 and SH2 women to have been ever diagnosed with breast cancer (BC; p = 0.002). TH were more likely to be diagnosed with ovarian cancer (OC) than SH2 (p = 0.017), but not SH1. Age at BC diagnosis was the same in TH vs. SH1 (p = 0.231), but was on average 4.5 years younger in TH than in SH2 (p < 0.001). BC in TH was more likely to be estrogen receptor (ER) positive (p = 0.010) or progesterone receptor (PR) positive (p = 0.013) than in SH1, but less likely to be ER positive (p < 0.001) or PR positive (p = 0.012) than SH2. Among 15 tumors from TH patients, there was no clear pattern of loss of heterozygosity (LOH) for BRCA1 or BRCA2 in either BC or OC.

CONCLUSIONS

Our observations suggest that clinical TH phenotypes resemble SH1. However, TH breast tumor marker characteristics are phenotypically intermediate to SH1 and SH2.

Fine-Scale Mapping at 9p22.2 Identifies Candidate Causal Variants That Modify Ovarian Cancer Risk in BRCA1 and BRCA2 Mutation Carriers

Population-based genome wide association studies have identified a locus at 9p22.2 associated with ovarian cancer risk, which also modifies ovarian cancer risk in BRCA1 and BRCA2 mutation carriers. We conducted fine-scale mapping at 9p22.2 to identify potential causal variants in BRCA1 and BRCA2 mutation carriers. Genotype data were available for 15,252 (2,462 ovarian cancer cases) BRCA1 and 8,211 (631 ovarian cancer cases) BRCA2 mutation carriers. Following genotype imputation, ovarian cancer associations were assessed for 4,873 and 5,020 SNPs in BRCA1 and BRCA 2 mutation carriers respectively, within a retrospective cohort analytical framework. In BRCA1 mutation carriers one set of eight correlated candidate causal variants for ovarian cancer risk modification was identified (top SNP rs10124837, HR: 0.73, 95%CI: 0.68 to 0.79, p-value 2× 10-16). These variants were located up to 20 kb upstream of BNC2. In BRCA2 mutation carriers one region, up to 45 kb upstream of BNC2, and containing 100 correlated SNPs was identified as candidate causal (top SNP rs62543585, HR: 0.69, 95%CI: 0.59 to 0.80, p-value 1.0 × 10-6). The candidate causal in BRCA1 mutation carriers did not include the strongest associated variant at this locus in the general population. In sum, we identified a set of candidate causal variants in a region that encompasses the BNC2 transcription start site. The ovarian cancer association at 9p22.2 may be mediated by different variants in BRCA1 mutation carriers and in the general population. Thus, potentially different mechanisms may underlie ovarian cancer risk for mutation carriers and the general population.

Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization

BACKGROUND

Epidemiological studies have linked adult height with breast cancer risk in women. However, the magnitude of the association, particularly by subtypes of breast cancer, has not been established. Furthermore, the mechanisms of the association remain unclear.

METHODS

We performed a meta-analysis to investigate associations between height and breast cancer risk using data from 159 prospective cohorts totaling 5216302 women, including 113178 events. In a consortium with individual-level data from 46325 case patients and 42482 control patients, we conducted a Mendelian randomization analysis using a genetic score that comprised 168 height-associated variants as an instrument. This association was further evaluated in a second consortium using summary statistics data from 16003 case patients and 41335 control patients.

RESULTS

The pooled relative risk of breast cancer was 1.17 (95% confidence interval [CI] = 1.15 to 1.19) per 10cm increase in height in the meta-analysis of prospective studies. In Mendelian randomization analysis, the odds ratio of breast cancer per 10cm increase in genetically predicted height was 1.22 (95% CI = 1.13 to 1.32) in the first consortium and 1.21 (95% CI = 1.05 to 1.39) in the second consortium. The association was found in both premenopausal and postmenopausal women but restricted to hormone receptor-positive breast cancer. Analyses of height-associated variants identified eight new loci associated with breast cancer risk after adjusting for multiple comparisons, including three loci at 1q21.2, DNAJC27, and CCDC91 at genome-wide significance level P < 5×10(-8).

CONCLUSIONS

Our study provides strong evidence that adult height is a risk factor for breast cancer in women and certain genetic factors and biological pathways affecting adult height have an important role in the etiology of breast cancer.

Prevalence of Germline Mutations in Genes Engaged in DNA Damage Repair by Homologous Recombination in Patients with Triple-Negative and Hereditary Non-Triple-Negative Breast Cancers

Purpose

This study sought to assess the prevalence of common germline mutations in several genes engaged in the repair of DNA double-strand break by homologous recombination in patients with triple-negative breast cancers and hereditary non-triple-negative breast cancers. Tumors deficient in this type of DNA damage repair are known to be especially sensitive to DNA cross-linking agents (e.g., platinum drugs) and to poly(ADP-ribose) polymerase (PARP) inhibitors.

Methods

Genetic testing was performed for 36 common germline mutations in genes engaged in the repair of DNA by homologous recombination, i.e., BRCA1, BRCA2, CHEK2, NBN, ATM, PALB2, BARD1, and RAD51D, in 202 consecutive patients with triple-negative breast cancers and hereditary non-triple-negative breast cancers.

Results

Thirty five (22.2%) of 158 patients in the triple-negative group carried mutations in genes involved in DNA repair by homologous recombination, while 10 (22.7%) of the 44 patients in the hereditary non-triple-negative group carried such mutations. Mutations in BRCA1 were most frequent in patients with triple-negative breast cancer (18.4%), and mutations in CHEK2 were most frequent in patients with hereditary non-triple-negative breast cancers (15.9%). In addition, in the triple-negative group, mutations in CHEK2, NBN, and ATM (3.8% combined) were found, while mutations in BRCA1, NBN, and PALB2 (6.8% combined) were identified in the hereditary non-triple-negative group.

Conclusions

Identifying mutations in genes engaged in DNA damage repair by homologous recombination other than BRCA1/2 can substantially increase the proportion of patients with triple-negative breast cancer and hereditary non-triple-negative breast cancer who may be eligible for therapy using PARP inhibitors and platinum drugs.

Fine-mapping identifies two additional breast cancer susceptibility loci at 9q31.2

We recently identified a novel susceptibility variant, rs865686, for estrogen-receptor positive breast cancer at 9q31.2. Here, we report a fine-mapping analysis of the 9q31.2 susceptibility locus using 43 160 cases and 42 600 controls of European ancestry ascertained from 52 studies and a further 5795 cases and 6624 controls of Asian ancestry from nine studies. Single nucleotide polymorphism (SNP) rs676256 was most strongly associated with risk in Europeans (odds ratios [OR] = 0.90 [0.88-0.92]; P-value = 1.58 × 10(-25)). This SNP is one of a cluster of highly correlated variants, including rs865686, that spans ∼14.5 kb. We identified two additional independent association signals demarcated by SNPs rs10816625 (OR = 1.12 [1.08-1.17]; P-value = 7.89 × 10(-09)) and rs13294895 (OR = 1.09 [1.06-1.12]; P-value = 2.97 × 10(-11)). SNP rs10816625, but not rs13294895, was also associated with risk of breast cancer in Asian individuals (OR = 1.12 [1.06-1.18]; P-value = 2.77 × 10(-05)). Functional genomic annotation using data derived from breast cancer cell-line models indicates that these SNPs localise to putative enhancer elements that bind known drivers of hormone-dependent breast cancer, including ER-α, FOXA1 and GATA-3. In vitro analyses indicate that rs10816625 and rs13294895 have allele-specific effects on enhancer activity and suggest chromatin interactions with the KLF4 gene locus. These results demonstrate the power of dense genotyping in large studies to identify independent susceptibility variants. Analysis of associations using subjects with different ancestry, combined with bioinformatic and genomic characterisation, can provide strong evidence for the likely causative alleles and their functional basis.

Identification of novel genetic markers of breast cancer survival

BACKGROUND

Survival after a diagnosis of breast cancer varies considerably between patients, and some of this variation may be because of germline genetic variation. We aimed to identify genetic markers associated with breast cancer-specific survival.

METHODS

We conducted a large meta-analysis of studies in populations of European ancestry, including 37954 patients with 2900 deaths from breast cancer. Each study had been genotyped for between 200000 and 900000 single nucleotide polymorphisms (SNPs) across the genome; genotypes for nine million common variants were imputed using a common reference panel from the 1000 Genomes Project. We also carried out subtype-specific analyses based on 6881 estrogen receptor (ER)-negative patients (920 events) and 23059 ER-positive patients (1333 events). All statistical tests were two-sided.

RESULTS

We identified one new locus (rs2059614 at 11q24.2) associated with survival in ER-negative breast cancer cases (hazard ratio [HR] = 1.95, 95% confidence interval [CI] = 1.55 to 2.47, P = 1.91 x 10(-8)). Genotyping a subset of 2113 case patients, of which 300 were ER negative, provided supporting evidence for the quality of the imputation. The association in this set of case patients was stronger for the observed genotypes than for the imputed genotypes. A second locus (rs148760487 at 2q24.2) was associated at genome-wide statistical significance in initial analyses; the association was similar in ER-positive and ER-negative case patients. Here the results of genotyping suggested that the finding was less robust.

CONCLUSIONS

This is currently the largest study investigating genetic variation associated with breast cancer survival. Our results have potential clinical implications, as they confirm that germline genotype can provide prognostic information in addition to standard tumor prognostic factors.

Prediction of breast cancer risk based on profiling with common genetic variants

BACKGROUND

Data for multiple common susceptibility alleles for breast cancer may be combined to identify women at different levels of breast cancer risk. Such stratification could guide preventive and screening strategies. However, empirical evidence for genetic risk stratification is lacking.

METHODS

We investigated the value of using 77 breast cancer-associated single nucleotide polymorphisms (SNPs) for risk stratification, in a study of 33 673 breast cancer cases and 33 381 control women of European origin. We tested all possible pair-wise multiplicative interactions and constructed a 77-SNP polygenic risk score (PRS) for breast cancer overall and by estrogen receptor (ER) status. Absolute risks of breast cancer by PRS were derived from relative risk estimates and UK incidence and mortality rates.

RESULTS

There was no strong evidence for departure from a multiplicative model for any SNP pair. Women in the highest 1% of the PRS had a three-fold increased risk of developing breast cancer compared with women in the middle quintile (odds ratio [OR] = 3.36, 95% confidence interval [CI] = 2.95 to 3.83). The ORs for ER-positive and ER-negative disease were 3.73 (95% CI = 3.24 to 4.30) and 2.80 (95% CI = 2.26 to 3.46), respectively. Lifetime risk of breast cancer for women in the lowest and highest quintiles of the PRS were 5.2% and 16.6% for a woman without family history, and 8.6% and 24.4% for a woman with a first-degree family history of breast cancer.

CONCLUSIONS

The PRS stratifies breast cancer risk in women both with and without a family history of breast cancer. The observed level of risk discrimination could inform targeted screening and prevention strategies. Further discrimination may be achievable through combining the PRS with lifestyle/environmental factors, although these were not considered in this report.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

INTRODUCTION

Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

METHODS

We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

RESULTS

We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

CONCLUSIONS

This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Common germline polymorphisms associated with breast cancer-specific survival

INTRODUCTION

Previous studies have identified common germline variants nominally associated with breast cancer survival. These associations have not been widely replicated in further studies. The purpose of this study was to evaluate the association of previously reported SNPs with breast cancer-specific survival using data from a pooled analysis of eight breast cancer survival genome-wide association studies (GWAS) from the Breast Cancer Association Consortium.

METHODS

A literature review was conducted of all previously published associations between common germline variants and three survival outcomes: breast cancer-specific survival, overall survival and disease-free survival. All associations that reached the nominal significance level of P value <0.05 were included. Single nucleotide polymorphisms that had been previously reported as nominally associated with at least one survival outcome were evaluated in the pooled analysis of over 37,000 breast cancer cases for association with breast cancer-specific survival. Previous associations were evaluated using a one-sided test based on the reported direction of effect.

RESULTS

Fifty-six variants from 45 previous publications were evaluated in the meta-analysis. Fifty-four of these were evaluated in the full set of 37,954 breast cancer cases with 2,900 events and the two additional variants were evaluated in a reduced sample size of 30,000 samples in order to ensure independence from the previously published studies. Five variants reached nominal significance (P <0.05) in the pooled GWAS data compared to 2.8 expected under the null hypothesis. Seven additional variants were associated (P <0.05) with ER-positive disease.

CONCLUSIONS

Although no variants reached genome-wide significance (P <5 x 10(-8)), these results suggest that there is some evidence of association between candidate common germline variants and breast cancer prognosis. Larger studies from multinational collaborations are necessary to increase the power to detect associations, between common variants and prognosis, at more stringent significance levels.

Association of type and location of BRCA1 and BRCA2 mutations with risk of breast and ovarian cancer

IMPORTANCE

Limited information about the relationship between specific mutations in BRCA1 or BRCA2 (BRCA1/2) and cancer risk exists.

OBJECTIVE

To identify mutation-specific cancer risks for carriers of BRCA1/2.

DESIGN, SETTING, AND PARTICIPANTS

Observational study of women who were ascertained between 1937 and 2011 (median, 1999) and found to carry disease-associated BRCA1 or BRCA2 mutations. The international sample comprised 19,581 carriers of BRCA1 mutations and 11,900 carriers of BRCA2 mutations from 55 centers in 33 countries on 6 continents. We estimated hazard ratios for breast and ovarian cancer based on mutation type, function, and nucleotide position. We also estimated RHR, the ratio of breast vs ovarian cancer hazard ratios. A value of RHR greater than 1 indicated elevated breast cancer risk; a value of RHR less than 1 indicated elevated ovarian cancer risk.

EXPOSURES

Mutations of BRCA1 or BRCA2.

MAIN OUTCOMES AND MEASURES

Breast and ovarian cancer risks.

RESULTS

Among BRCA1 mutation carriers, 9052 women (46%) were diagnosed with breast cancer, 2317 (12%) with ovarian cancer, 1041 (5%) with breast and ovarian cancer, and 7171 (37%) without cancer. Among BRCA2 mutation carriers, 6180 women (52%) were diagnosed with breast cancer, 682 (6%) with ovarian cancer, 272 (2%) with breast and ovarian cancer, and 4766 (40%) without cancer. In BRCA1, we identified 3 breast cancer cluster regions (BCCRs) located at c.179 to c.505 (BCCR1; RHR = 1.46; 95% CI, 1.22-1.74; P = 2 × 10(-6)), c.4328 to c.4945 (BCCR2; RHR = 1.34; 95% CI, 1.01-1.78; P = .04), and c. 5261 to c.5563 (BCCR2', RHR = 1.38; 95% CI, 1.22-1.55; P = 6 × 10(-9)). We also identified an ovarian cancer cluster region (OCCR) from c.1380 to c.4062 (approximately exon 11) with RHR = 0.62 (95% CI, 0.56-0.70; P = 9 × 10(-17)). In BRCA2, we observed multiple BCCRs spanning c.1 to c.596 (BCCR1; RHR = 1.71; 95% CI, 1.06-2.78; P = .03), c.772 to c.1806 (BCCR1'; RHR = 1.63; 95% CI, 1.10-2.40; P = .01), and c.7394 to c.8904 (BCCR2; RHR = 2.31; 95% CI, 1.69-3.16; P = .00002). We also identified 3 OCCRs: the first (OCCR1) spanned c.3249 to c.5681 that was adjacent to c.5946delT (6174delT; RHR = 0.51; 95% CI, 0.44-0.60; P = 6 × 10(-17)). The second OCCR spanned c.6645 to c.7471 (OCCR2; RHR = 0.57; 95% CI, 0.41-0.80; P = .001). Mutations conferring nonsense-mediated decay were associated with differential breast or ovarian cancer risks and an earlier age of breast cancer diagnosis for both BRCA1 and BRCA2 mutation carriers.

CONCLUSIONS AND RELEVANCE

Breast and ovarian cancer risks varied by type and location of BRCA1/2 mutations. With appropriate validation, these data may have implications for risk assessment and cancer prevention decision making for carriers of BRCA1 and BRCA2 mutations.

Assessing associations between the AURKA-HMMR-TPX2-TUBG1 functional module and breast cancer risk in BRCA1/2 mutation carriers

While interplay between BRCA1 and AURKA-RHAMM-TPX2-TUBG1 regulates mammary epithelial polarization, common genetic variation in HMMR (gene product RHAMM) may be associated with risk of breast cancer in BRCA1 mutation carriers. Following on these observations, we further assessed the link between the AURKA-HMMR-TPX2-TUBG1 functional module and risk of breast cancer in BRCA1 or BRCA2 mutation carriers. Forty-one single nucleotide polymorphisms (SNPs) were genotyped in 15,252 BRCA1 and 8,211 BRCA2 mutation carriers and subsequently analyzed using a retrospective likelihood approach. The association of HMMR rs299290 with breast cancer risk in BRCA1 mutation carriers was confirmed: per-allele hazard ratio (HR) = 1.10, 95% confidence interval (CI) 1.04-1.15, p = 1.9 x 10(-4) (false discovery rate (FDR)-adjusted p = 0.043). Variation in CSTF1, located next to AURKA, was also found to be associated with breast cancer risk in BRCA2 mutation carriers: rs2426618 per-allele HR = 1.10, 95% CI 1.03-1.16, p = 0.005 (FDR-adjusted p = 0.045). Assessment of pairwise interactions provided suggestions (FDR-adjusted pinteraction values > 0.05) for deviations from the multiplicative model for rs299290 and CSTF1 rs6064391, and rs299290 and TUBG1 rs11649877 in both BRCA1 and BRCA2 mutation carriers. Following these suggestions, the expression of HMMR and AURKA or TUBG1 in sporadic breast tumors was found to potentially interact, influencing patients' survival. Together, the results of this study support the hypothesis of a causative link between altered function of AURKA-HMMR-TPX2-TUBG1 and breast carcinogenesis in BRCA1/2 mutation carriers.

Inherited variants in the inner centromere protein (INCENP) gene of the chromosomal passenger complex contribute to the susceptibility of ER-negative breast cancer

The chromosomal passenger complex (CPC) plays a pivotal role in the regulation of cell division. Therefore, inherited CPC variability could influence tumor development. The present candidate gene approach investigates the relationship between single nucleotide polymorphisms (SNPs) in genes encoding key CPC components and breast cancer risk. Fifteen SNPs in four CPC genes (INCENP, AURKB, BIRC5 and CDCA8) were genotyped in 88 911 European women from 39 case-control studies of the Breast Cancer Association Consortium. Possible associations were investigated in fixed-effects meta-analyses. The synonymous SNP rs1675126 in exon 7 of INCENP was associated with overall breast cancer risk [per A allele odds ratio (OR) 0.95, 95% confidence interval (CI) 0.92-0.98, P = 0.007] and particularly with estrogen receptor (ER)-negative breast tumors (per A allele OR 0.89, 95% CI 0.83-0.95, P = 0.0005). SNPs not directly genotyped were imputed based on 1000 Genomes. The SNPs rs1047739 in the 3' untranslated region and rs144045115 downstream of INCENP showed the strongest association signals for overall (per T allele OR 1.03, 95% CI 1.00-1.06, P = 0.0009) and ER-negative breast cancer risk (per A allele OR 1.06, 95% CI 1.02-1.10, P = 0.0002). Two genotyped SNPs in BIRC5 were associated with familial breast cancer risk (top SNP rs2071214: per G allele OR 1.12, 95% CI 1.04-1.21, P = 0.002). The data suggest that INCENP in the CPC pathway contributes to ER-negative breast cancer susceptibility in the European population. In spite of a modest contribution of CPC-inherited variants to the total burden of sporadic and familial breast cancer, their potential as novel targets for breast cancer treatment should be further investigated.

Fine-scale mapping of the 5q11.2 breast cancer locus reveals at least three independent risk variants regulating MAP3K1

Genome-wide association studies (GWASs) have revealed SNP rs889312 on 5q11.2 to be associated with breast cancer risk in women of European ancestry. In an attempt to identify the biologically relevant variants, we analyzed 909 genetic variants across 5q11.2 in 103,991 breast cancer individuals and control individuals from 52 studies in the Breast Cancer Association Consortium. Multiple logistic regression analyses identified three independent risk signals: the strongest associations were with 15 correlated variants (iCHAV1), where the minor allele of the best candidate, rs62355902, associated with significantly increased risks of both estrogen-receptor-positive (ER(+): odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.21-1.27, ptrend = 5.7 × 10(-44)) and estrogen-receptor-negative (ER(-): OR = 1.10, 95% CI = 1.05-1.15, ptrend = 3.0 × 10(-4)) tumors. After adjustment for rs62355902, we found evidence of association of a further 173 variants (iCHAV2) containing three subsets with a range of effects (the strongest was rs113317823 [pcond = 1.61 × 10(-5)]) and five variants composing iCHAV3 (lead rs11949391; ER(+): OR = 0.90, 95% CI = 0.87-0.93, pcond = 1.4 × 10(-4)). Twenty-six percent of the prioritized candidate variants coincided with four putative regulatory elements that interact with the MAP3K1 promoter through chromatin looping and affect MAP3K1 promoter activity. Functional analysis indicated that the cancer risk alleles of four candidates (rs74345699 and rs62355900 [iCHAV1], rs16886397 [iCHAV2a], and rs17432750 [iCHAV3]) increased MAP3K1 transcriptional activity. Chromatin immunoprecipitation analysis revealed diminished GATA3 binding to the minor (cancer-protective) allele of rs17432750, indicating a mechanism for its action. We propose that the cancer risk alleles act to increase MAP3K1 expression in vivo and might promote breast cancer cell survival.

Candidate genetic modifiers for breast and ovarian cancer risk in BRCA1 and BRCA2 mutation carriers

BACKGROUND

BRCA1 and BRCA2 mutation carriers are at substantially increased risk for developing breast and ovarian cancer. The incomplete penetrance coupled with the variable age at diagnosis in carriers of the same mutation suggests the existence of genetic and nongenetic modifying factors. In this study, we evaluated the putative role of variants in many candidate modifier genes.

METHODS

Genotyping data from 15,252 BRCA1 and 8,211 BRCA2 mutation carriers, for known variants (n = 3,248) located within or around 445 candidate genes, were available through the iCOGS custom-designed array. Breast and ovarian cancer association analysis was performed within a retrospective cohort approach.

RESULTS

The observed P values of association ranged between 0.005 and 1.000. None of the variants was significantly associated with breast or ovarian cancer risk in either BRCA1 or BRCA2 mutation carriers, after multiple testing adjustments.

CONCLUSION

There is little evidence that any of the evaluated candidate variants act as modifiers of breast and/or ovarian cancer risk in BRCA1 or BRCA2 mutation carriers.

IMPACT

Genome-wide association studies have been more successful at identifying genetic modifiers of BRCA1/2 penetrance than candidate gene studies.

Identification and characterization of novel associations in the CASP8/ALS2CR12 region on chromosome 2 with breast cancer risk

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.

Associations of common breast cancer susceptibility alleles with risk of breast cancer subtypes in BRCA1 and BRCA2 mutation carriers

INTRODUCTION

More than 70 common alleles are known to be involved in breast cancer (BC) susceptibility, and several exhibit significant heterogeneity in their associations with different BC subtypes. Although there are differences in the association patterns between BRCA1 and BRCA2 mutation carriers and the general population for several loci, no study has comprehensively evaluated the associations of all known BC susceptibility alleles with risk of BC subtypes in BRCA1 and BRCA2 carriers.

METHODS

We used data from 15,252 BRCA1 and 8,211 BRCA2 carriers to analyze the associations between approximately 200,000 genetic variants on the iCOGS array and risk of BC subtypes defined by estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and triple-negative- (TN) status; morphologic subtypes; histological grade; and nodal involvement.

RESULTS

The estimated BC hazard ratios (HRs) for the 74 known BC alleles in BRCA1 carriers exhibited moderate correlations with the corresponding odds ratios from the general population. However, their associations with ER-positive BC in BRCA1 carriers were more consistent with the ER-positive associations in the general population (intraclass correlation (ICC) = 0.61, 95% confidence interval (CI): 0.45 to 0.74), and the same was true when considering ER-negative associations in both groups (ICC = 0.59, 95% CI: 0.42 to 0.72). Similarly, there was strong correlation between the ER-positive associations for BRCA1 and BRCA2 carriers (ICC = 0.67, 95% CI: 0.52 to 0.78), whereas ER-positive associations in any one of the groups were generally inconsistent with ER-negative associations in any of the others. After stratifying by ER status in mutation carriers, additional significant associations were observed. Several previously unreported variants exhibited associations at P <10(-6) in the analyses by PR status, HER2 status, TN phenotype, morphologic subtypes, histological grade and nodal involvement.

CONCLUSIONS

Differences in associations of common BC susceptibility alleles between BRCA1 and BRCA2 carriers and the general population are explained to a large extent by differences in the prevalence of ER-positive and ER-negative tumors. Estimates of the risks associated with these variants based on population-based studies are likely to be applicable to mutation carriers after taking ER status into account, which has implications for risk prediction.

Refined histopathological predictors of BRCA1 and BRCA2 mutation status: a large-scale analysis of breast cancer characteristics from the BCAC, CIMBA, and ENIGMA consortia

INTRODUCTION

The distribution of histopathological features of invasive breast tumors in BRCA1 or BRCA2 germline mutation carriers differs from that of individuals with no known mutation. Histopathological features thus have utility for mutation prediction, including statistical modeling to assess pathogenicity of BRCA1 or BRCA2 variants of uncertain clinical significance. We analyzed large pathology datasets accrued by the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC) to reassess histopathological predictors of BRCA1 and BRCA2 mutation status, and provide robust likelihood ratio (LR) estimates for statistical modeling.

METHODS

Selection criteria for study/center inclusion were estrogen receptor (ER) status or grade data available for invasive breast cancer diagnosed younger than 70 years. The dataset included 4,477 BRCA1 mutation carriers, 2,565 BRCA2 mutation carriers, and 47,565 BCAC breast cancer cases. Country-stratified estimates of the likelihood of mutation status by histopathological markers were derived using a Mantel-Haenszel approach.

RESULTS

ER-positive phenotype negatively predicted BRCA1 mutation status, irrespective of grade (LRs from 0.08 to 0.90). ER-negative grade 3 histopathology was more predictive of positive BRCA1 mutation status in women 50 years or older (LR = 4.13 (3.70 to 4.62)) versus younger than 50 years (LR = 3.16 (2.96 to 3.37)). For BRCA2, ER-positive grade 3 phenotype modestly predicted positive mutation status irrespective of age (LR = 1.7-fold), whereas ER-negative grade 3 features modestly predicted positive mutation status at 50 years or older (LR = 1.54 (1.27 to 1.88)). Triple-negative tumor status was highly predictive of BRCA1 mutation status for women younger than 50 years (LR = 3.73 (3.43 to 4.05)) and 50 years or older (LR = 4.41 (3.86 to 5.04)), and modestly predictive of positive BRCA2 mutation status in women 50 years or older (LR = 1.79 (1.42 to 2.24)).

CONCLUSIONS

These results refine likelihood-ratio estimates for predicting BRCA1 and BRCA2 mutation status by using commonly measured histopathological features. Age at diagnosis is an important variable for most analyses, and grade is more informative than ER status for BRCA2 mutation carrier prediction. The estimates will improve BRCA1 and BRCA2 variant classification and inform patient mutation testing and clinical management.

Genetic variation in mitotic regulatory pathway genes is associated with breast tumor grade

Mitotic index is an important component of histologic grade and has an etiologic role in breast tumorigenesis. Several small candidate gene studies have reported associations between variation in mitotic genes and breast cancer risk. We measured associations between 2156 single nucleotide polymorphisms (SNPs) from 194 mitotic genes and breast cancer risk, overall and by histologic grade, in the Breast Cancer Association Consortium (BCAC) iCOGS study (n = 39 067 cases; n = 42 106 controls). SNPs in TACC2 [rs17550038: odds ratio (OR) = 1.24, 95% confidence interval (CI) 1.16-1.33, P = 4.2 × 10(-10)) and EIF3H (rs799890: OR = 1.07, 95% CI 1.04-1.11, P = 8.7 × 10(-6)) were significantly associated with risk of low-grade breast cancer. The TACC2 signal was retained (rs17550038: OR = 1.15, 95% CI 1.07-1.23, P = 7.9 × 10(-5)) after adjustment for breast cancer risk SNPs in the nearby FGFR2 gene, suggesting that TACC2 is a novel, independent genome-wide significant genetic risk locus for low-grade breast cancer. While no SNPs were individually associated with high-grade disease, a pathway-level gene set analysis showed that variation across the 194 mitotic genes was associated with high-grade breast cancer risk (P = 2.1 × 10(-3)). These observations will provide insight into the contribution of mitotic defects to histological grade and the etiology of breast cancer.

MicroRNA related polymorphisms and breast cancer risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology.

Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study

INTRODUCTION

We have previously shown that a tag single nucleotide polymorphism (rs10235235), which maps to the CYP3A locus (7q22.1), was associated with a reduction in premenopausal urinary estrone glucuronide levels and a modest reduction in risk of breast cancer in women age ≤50 years.

METHODS

We further investigated the association of rs10235235 with breast cancer risk in a large case control study of 47,346 cases and 47,570 controls from 52 studies participating in the Breast Cancer Association Consortium. Genotyping of rs10235235 was conducted using a custom Illumina Infinium array. Stratified analyses were conducted to determine whether this association was modified by age at diagnosis, ethnicity, age at menarche or tumor characteristics.

RESULTS

We confirmed the association of rs10235235 with breast cancer risk for women of European ancestry but found no evidence that this association differed with age at diagnosis. Heterozygote and homozygote odds ratios (ORs) were OR = 0.98 (95% CI 0.94, 1.01; P = 0.2) and OR = 0.80 (95% CI 0.69, 0.93; P = 0.004), respectively (P(trend) = 0.02). There was no evidence of effect modification by tumor characteristics. rs10235235 was, however, associated with age at menarche in controls (P(trend) = 0.005) but not cases (P(trend) = 0.97). Consequently the association between rs10235235 and breast cancer risk differed according to age at menarche (P(het) = 0.02); the rare allele of rs10235235 was associated with a reduction in breast cancer risk for women who had their menarche age ≥15 years (OR(het) = 0.84, 95% CI 0.75, 0.94; OR(hom) = 0.81, 95% CI 0.51, 1.30; P(trend) = 0.002) but not for those who had their menarche age ≤11 years (OR(het) = 1.06, 95% CI 0.95, 1.19, OR(hom) = 1.07, 95% CI 0.67, 1.72; P(trend) = 0.29).

CONCLUSIONS

To our knowledge rs10235235 is the first single nucleotide polymorphism to be associated with both breast cancer risk and age at menarche consistent with the well-documented association between later age at menarche and a reduction in breast cancer risk. These associations are likely mediated via an effect on circulating hormone levels.

DNA glycosylases involved in base excision repair may be associated with cancer risk in BRCA1 and BRCA2 mutation carriers

Single Nucleotide Polymorphisms (SNPs) in genes involved in the DNA Base Excision Repair (BER) pathway could be associated with cancer risk in carriers of mutations in the high-penetrance susceptibility genes BRCA1 and BRCA2, given the relation of synthetic lethality that exists between one of the components of the BER pathway, PARP1 (poly ADP ribose polymerase), and both BRCA1 and BRCA2. In the present study, we have performed a comprehensive analysis of 18 genes involved in BER using a tagging SNP approach in a large series of BRCA1 and BRCA2 mutation carriers. 144 SNPs were analyzed in a two stage study involving 23,463 carriers from the CIMBA consortium (the Consortium of Investigators of Modifiers of BRCA1 and BRCA2). Eleven SNPs showed evidence of association with breast and/or ovarian cancer at p<0.05 in the combined analysis. Four of the five genes for which strongest evidence of association was observed were DNA glycosylases. The strongest evidence was for rs1466785 in the NEIL2 (endonuclease VIII-like 2) gene (HR: 1.09, 95% CI (1.03-1.16), p = 2.7 × 10(-3)) for association with breast cancer risk in BRCA2 mutation carriers, and rs2304277 in the OGG1 (8-guanine DNA glycosylase) gene, with ovarian cancer risk in BRCA1 mutation carriers (HR: 1.12 95%CI: 1.03-1.21, p = 4.8 × 10(-3)). DNA glycosylases involved in the first steps of the BER pathway may be associated with cancer risk in BRCA1/2 mutation carriers and should be more comprehensively studied.

A large-scale assessment of two-way SNP interactions in breast cancer susceptibility using 46,450 cases and 42,461 controls from the breast cancer association consortium

Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70,917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46,450 breast cancer cases and 42,461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10(-4)) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10(-8). Results from the second analytic approach were consistent with those from the first (P > 10(-10)). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

Genetic predisposition to in situ and invasive lobular carcinoma of the breast

Invasive lobular breast cancer (ILC) accounts for 10-15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09-1.18), P = 6.0 × 10(-10); P-het for ILC vs IDC ER+ tumors = 1.8 × 10(-4)). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes.

FGF receptor genes and breast cancer susceptibility: results from the Breast Cancer Association Consortium

BACKGROUND

Breast cancer is one of the most common malignancies in women. Genome-wide association studies have identified FGFR2 as a breast cancer susceptibility gene. Common variation in other fibroblast growth factor (FGF) receptors might also modify risk. We tested this hypothesis by studying genotyped single-nucleotide polymorphisms (SNPs) and imputed SNPs in FGFR1, FGFR3, FGFR4 and FGFRL1 in the Breast Cancer Association Consortium.

METHODS

Data were combined from 49 studies, including 53 835 cases and 50 156 controls, of which 89 050 (46 450 cases and 42 600 controls) were of European ancestry, 12 893 (6269 cases and 6624 controls) of Asian and 2048 (1116 cases and 932 controls) of African ancestry. Associations with risk of breast cancer, overall and by disease sub-type, were assessed using unconditional logistic regression.

RESULTS

Little evidence of association with breast cancer risk was observed for SNPs in the FGF receptor genes. The strongest evidence in European women was for rs743682 in FGFR3; the estimated per-allele odds ratio was 1.05 (95% confidence interval=1.02-1.09, P=0.0020), which is substantially lower than that observed for SNPs in FGFR2.

CONCLUSION

Our results suggest that common variants in the other FGF receptors are not associated with risk of breast cancer to the degree observed for FGFR2.

Fine-scale mapping of the FGFR2 breast cancer risk locus: putative functional variants differentially bind FOXA1 and E2F1

The 10q26 locus in the second intron of FGFR2 is the locus most strongly associated with estrogen-receptor-positive breast cancer in genome-wide association studies. We conducted fine-scale mapping in case-control studies genotyped with a custom chip (iCOGS), comprising 41 studies (n = 89,050) of European ancestry, 9 Asian ancestry studies (n = 13,983), and 2 African ancestry studies (n = 2,028) from the Breast Cancer Association Consortium. We identified three statistically independent risk signals within the locus. Within risk signals 1 and 3, genetic analysis identified five and two variants, respectively, highly correlated with the most strongly associated SNPs. By using a combination of genetic fine mapping, data on DNase hypersensitivity, and electrophoretic mobility shift assays to study protein-DNA binding, we identified rs35054928, rs2981578, and rs45631563 as putative functional SNPs. Chromatin immunoprecipitation showed that FOXA1 preferentially bound to the risk-associated allele (C) of rs2981578 and was able to recruit ERα to this site in an allele-specific manner, whereas E2F1 preferentially bound the risk variant of rs35054928. The risk alleles were preferentially found in open chromatin and bound by Ser5 phosphorylated RNA polymerase II, suggesting that the risk alleles are associated with changes in transcription. Chromatin conformation capture demonstrated that the risk region was able to interact with the promoter of FGFR2, the likely target gene of this risk region. A role for FOXA1 in mediating breast cancer susceptibility at this locus is consistent with the finding that the FGFR2 risk locus primarily predisposes to estrogen-receptor-positive disease.

Can selenium levels act as a marker of colorectal cancer risk?

Background

Selenium has attracted attention because of its antioxidant properties. Antioxidants protects cells from damage. Certain breakdown products of selenium are believed to prevent tumor growth by enhancing the immune cell activity and suppressing the development of tumor blood vessels. In this observational study, selenium level was measured in a series of patients from Poland and Estonia to determine a correlation between levels of this microelement and colorectal cancer risk.

Methods

A total of 169 colorectal cancer patients and 169 healthy controls were enrolled in the study after obtaining their informed consent. Selenium level in the blood serum was measured using Graphite Furnace Atomic Absorption Spectrometry (GFAAS). The statistical analysis was performed by Fisher’s exact test.

Results

The threshold point of selenium level was 55 μg/l and 65 μg/l for Poland and Estonia respectively, for an increase in cancer risk. The lower levels of selenium were associated with greater risk of colorectal cancer.

Conclusions

The result reveals a significant strong association between low selenium level and the colorectal cancer risk in both Estonian and Polish populations.

Functional variants at the 11q13 risk locus for breast cancer regulate cyclin D1 expression through long-range enhancers

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1.

Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer

TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ∼480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10(-7)), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10(-8)) and BRCA1 mutation carrier (P = 1.1 × 10(-5)) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10(-14)), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10(-15)) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10(-12)) and BRCA1 mutation carrier (P = 1.6 × 10(-14)) breast and invasive ovarian (P = 1.3 × 10(-11)) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant.

Large-scale genotyping identifies 41 new loci associated with breast cancer risk

Breast cancer is the most common cancer among women. Common variants at 27 loci have been identified as associated with susceptibility to breast cancer, and these account for ∼9% of the familial risk of the disease. We report here a meta-analysis of 9 genome-wide association studies, including 10,052 breast cancer cases and 12,575 controls of European ancestry, from which we selected 29,807 SNPs for further genotyping. These SNPs were genotyped in 45,290 cases and 41,880 controls of European ancestry from 41 studies in the Breast Cancer Association Consortium (BCAC). The SNPs were genotyped as part of a collaborative genotyping experiment involving four consortia (Collaborative Oncological Gene-environment Study, COGS) and used a custom Illumina iSelect genotyping array, iCOGS, comprising more than 200,000 SNPs. We identified SNPs at 41 new breast cancer susceptibility loci at genome-wide significance (P < 5 × 10(-8)). Further analyses suggest that more than 1,000 additional loci are involved in breast cancer susceptibility.

Genome-wide association studies identify four ER negative-specific breast cancer risk loci

Estrogen receptor (ER)-negative tumors represent 20-30% of all breast cancers, with a higher proportion occurring in younger women and women of African ancestry. The etiology and clinical behavior of ER-negative tumors are different from those of tumors expressing ER (ER positive), including differences in genetic predisposition. To identify susceptibility loci specific to ER-negative disease, we combined in a meta-analysis 3 genome-wide association studies of 4,193 ER-negative breast cancer cases and 35,194 controls with a series of 40 follow-up studies (6,514 cases and 41,455 controls), genotyped using a custom Illumina array, iCOGS, developed by the Collaborative Oncological Gene-environment Study (COGS). SNPs at four loci, 1q32.1 (MDM4, P = 2.1 × 10(-12) and LGR6, P = 1.4 × 10(-8)), 2p24.1 (P = 4.6 × 10(-8)) and 16q12.2 (FTO, P = 4.0 × 10(-8)), were associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancers.

Genome-wide association study in BRCA1 mutation carriers identifies novel loci associated with breast and ovarian cancer risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7 × 10(-8), HR = 1.14, 95% CI: 1.09-1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4 × 10(-8), HR = 1.27, 95% CI: 1.17-1.38) and 4q32.3 (rs4691139, P = 3.4 × 10(-8), HR = 1.20, 95% CI: 1.17-1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10(-4)). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%-50% compared to 81%-100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.

Comparison of 6q25 breast cancer hits from Asian and European Genome Wide Association Studies in the Breast Cancer Association Consortium (BCAC)

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER-) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26-1.48), p = 7.6 × 10(-14) in Asians and 1.09 (95% CI 1.07-1.11), p = 6.8 × 10(-18) in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19-1.41), p = 1.2 × 10(-9) in Asians and 1.12 (95% CI 1.08-1.17), p = 3.8 × 10(-9) in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER- than ER+ tumours in Europeans [OR (ER-) = 1.20 (95% CI 1.15-1.25), p = 1.8 × 10(-17) versus OR (ER+) = 1.07 (95% CI 1.04-1.1), p = 1.3 × 10(-7), p(heterogeneity) = 5.1 × 10(-6)]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER- tumours.

A nonsynonymous polymorphism in IRS1 modifies risk of developing breast and ovarian cancers in BRCA1 and ovarian cancer in BRCA2 mutation carriers

BACKGROUND

We previously reported significant associations between genetic variants in insulin receptor substrate 1 (IRS1) and breast cancer risk in women carrying BRCA1 mutations. The objectives of this study were to investigate whether the IRS1 variants modified ovarian cancer risk and were associated with breast cancer risk in a larger cohort of BRCA1 and BRCA2 mutation carriers.

METHODS

IRS1 rs1801123, rs1330645, and rs1801278 were genotyped in samples from 36 centers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Data were analyzed by a retrospective cohort approach modeling the associations with breast and ovarian cancer risks simultaneously. Analyses were stratified by BRCA1 and BRCA2 status and mutation class in BRCA1 carriers.

RESULTS

Rs1801278 (Gly972Arg) was associated with ovarian cancer risk for both BRCA1 (HR, 1.43; 95% confidence interval (CI), 1.06-1.92; P = 0.019) and BRCA2 mutation carriers (HR, 2.21; 95% CI, 1.39-3.52, P = 0.0008). For BRCA1 mutation carriers, the breast cancer risk was higher in carriers with class II mutations than class I mutations (class II HR, 1.86; 95% CI, 1.28-2.70; class I HR, 0.86; 95%CI, 0.69-1.09; P(difference), 0.0006). Rs13306465 was associated with ovarian cancer risk in BRCA1 class II mutation carriers (HR, 2.42; P = 0.03).

CONCLUSION

The IRS1 Gly972Arg single-nucleotide polymorphism, which affects insulin-like growth factor and insulin signaling, modifies ovarian cancer risk in BRCA1 and BRCA2 mutation carriers and breast cancer risk in BRCA1 class II mutation carriers.

IMPACT

These findings may prove useful for risk prediction for breast and ovarian cancers in BRCA1 and BRCA2 mutation carriers.