Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk.

Background: Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS). Methods: We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan. Results: Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62). Conclusions: This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

A Breast Cancer Candidate Locus at 6q Narrowed to 6q15-q21.

Although a number of high-risk breast cancer genes have been identified, including BRCA1 and BRCA2, the risk profile of many high-risk families cannot be explained using known breast cancer genes. Previously, we have shown strong indications of new breast cancer risk loci at chromosomes 2p, 6q, and 14q in a family of six generations including 10 breast cancer cases. In this study, we identified and traced four new family branches descending from siblings of the parents in the top generation of the studied family. One distantly related branch included four breast cancer cases, two of whom were diagnosed at age < 45 years. DNA samples from the cases were typed at selected polymorphic markers from all three chromosome loci, to test identical origin of the haplotypes. All four cases were shown to segregate a common 6q haplotype with a region identical to the previously identified 6q haplotype. The data strongly support a new breast cancer locus at 6q, and narrow it down to a 17 MB interval at 6q15-q21.

High Atlastin 2-2 (ATL2-2) Expression Associates with Worse Prognosis in Estrogen-Receptor-Positive Breast Cancer.

The disruption of endoplasmic reticulum (ER) homeostasis occurs in many human diseases. Atlastins (ATLs) maintain the branched network of the ER. The dysregulation of ATL2, located at ER network junctions, has been associated with cancer. ATL2 is necessary for lipid droplet formation in murine breast tissue. Thus, we analyzed whether ATL2 has a role in human breast cancer (BC) pathology. The expression of ATL2 variant ATL2-2 was analyzed in breast tumors from the BC cohorts of the TCGA, METABRIC, and two independent Icelandic cohorts, Cohort 1 and 2; its association with clinical, pathological, survival, and cellular pathways was explored. ATL2-2 mRNA and protein expression were higher in breast tumors than in normal tissue. ATL2-2 mRNA associated with tumor characteristics that indicate a worse prognosis. In METABRIC, high ATL2-2 mRNA levels were associated with shorter BC-specific survival (BCSS) in patients with estrogen-receptor-positive luminal breast tumors, which remained significant after correction for grade and tumor size (HR 1.334, CI 1.063-1.673). Tumors with high ATL2 mRNA showed an upregulation of hallmark pathways MYC targets v1, E2F targets, and G2M checkpoint genes. Taken together, the results suggest that high levels of ATL2-2 may support BC progression through key cancer driver pathways.

Ovarian cancer pathology characteristics as predictors of variant pathogenicity in BRCA1 and BRCA2.

BACKGROUND: The distribution of ovarian tumour characteristics differs between germline BRCA1 and BRCA2 pathogenic variant carriers and non-carriers. In this study, we assessed the utility of ovarian tumour characteristics as predictors of BRCA1 and BRCA2 variant pathogenicity, for application using the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) variant classification system. METHODS: Data for 10,373 ovarian cancer cases, including carriers and non-carriers of BRCA1 or BRCA2 pathogenic variants, were collected from unpublished international cohorts and consortia and published studies. Likelihood ratios (LR) were calculated for the association of ovarian cancer histology and other characteristics, with BRCA1 and BRCA2 variant pathogenicity. Estimates were aligned to ACMG/AMP code strengths (supporting, moderate, strong). RESULTS: No histological subtype provided informative ACMG/AMP evidence in favour of BRCA1 and BRCA2 variant pathogenicity. Evidence against variant pathogenicity was estimated for the mucinous and clear cell histologies (supporting) and borderline cases (moderate). Refined associations are provided according to tumour grade, invasion and age at diagnosis. CONCLUSIONS: We provide detailed estimates for predicting BRCA1 and BRCA2 variant pathogenicity based on ovarian tumour characteristics. This evidence can be combined with other variant information under the ACMG/AMP classification system, to improve classification and carrier clinical management.

Breast and Prostate Cancer Risks for Male BRCA1 and BRCA2 Pathogenic Variant Carriers Using Polygenic Risk Scores.

BACKGROUND: Recent population-based female breast cancer and prostate cancer polygenic risk scores (PRS) have been developed. We assessed the associations of these PRS with breast and prostate cancer risks for male BRCA1 and BRCA2 pathogenic variant carriers. METHODS: 483 BRCA1 and 1318 BRCA2 European ancestry male carriers were available from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). A 147-single nucleotide polymorphism (SNP) prostate cancer PRS (PRSPC) and a 313-SNP breast cancer PRS were evaluated. There were 3 versions of the breast cancer PRS, optimized to predict overall (PRSBC), estrogen receptor (ER)-negative (PRSER-), or ER-positive (PRSER+) breast cancer risk. RESULTS: PRSER+ yielded the strongest association with breast cancer risk. The odds ratios (ORs) per PRSER+ standard deviation estimates were 1.40 (95% confidence interval [CI] =1.07 to 1.83) for BRCA1 and 1.33 (95% CI = 1.16 to 1.52) for BRCA2 carriers. PRSPC was associated with prostate cancer risk for BRCA1 (OR = 1.73, 95% CI = 1.28 to 2.33) and BRCA2 (OR = 1.60, 95% CI = 1.34 to 1.91) carriers. The estimated breast cancer odds ratios were larger after adjusting for female relative breast cancer family history. By age 85 years, for BRCA2 carriers, the breast cancer risk varied from 7.7% to 18.4% and prostate cancer risk from 34.1% to 87.6% between the 5th and 95th percentiles of the PRS distributions. CONCLUSIONS: Population-based prostate and female breast cancer PRS are associated with a wide range of absolute breast and prostate cancer risks for male BRCA1 and BRCA2 carriers. These findings warrant further investigation aimed at providing personalized cancer risks for male carriers and informing clinical management.

Hsa-miR-21-3p associates with breast cancer patient survival and targets genes in tumor suppressive pathways.

Breast cancer is the cancer most often diagnosed in women. MicroRNAs (MIRs) are short RNA molecules that bind mRNA resulting in their downregulation. MIR21 has been shown to be an oncomiR in most cancer types, including breast cancer. Most of the effects of miR-21 have been attributed to hsa-miR-21-5p that is transcribed from the leading strand of MIR21, but hsa-miR-21-3p (miR-21-3p), transcribed from the lagging strand, is much less studied. The aim of the study is to analyze whether expression of miR-21-3p is prognostic for breast cancer. MiR-21-3p association with survival, clinical and pathological characteristics was analyzed in a large breast cancer cohort and validated in three separate cohorts, including TCGA and METABRIC. Analytical tools were also used to infer miR-21-3p function and to identify potential target genes and functional pathways. The results showed that in the exploration cohort, high miR-21-3p levels associated with shorter survival and lymph node positivity. In the three validation cohorts, high miR-21-3p levels associated with pathological characteristics that predict worse prognosis. Specifically, in the largest validation cohort, METABRIC (n = 1174), high miR-21-3p levels associated with large tumors, a high grade, lymph node and HER2 positivity, and shorter breast-cancer-specific survival (HR = 1.38, CI 1.13-1.68). This association remained significant after adjusting for confounding factors. The genes with expression levels that correlated with miR-21-3p were enriched in particular pathways, including the epithelial-to-mesenchymal transition and proliferation. Among the most significantly downregulated targets were MAT2A and the tumor suppressive genes STARD13 and ZNF132. The results from this study emphasize that both 3p- and 5p-arms from a MIR warrant independent study. The data show that miR-21-3p overexpression in breast tumors is a marker of worse breast cancer progression and it affects genes in pathways that drive breast cancer by down-regulating tumor suppressor genes. The results suggest miR-21-3p as a potential biomarker.

The predictive ability of the 313 variant-based polygenic risk score for contralateral breast cancer risk prediction in women of European ancestry with a heterozygous BRCA1 or BRCA2 pathogenic variant.

PURPOSE: To evaluate the association between a previously published 313 variant-based breast cancer (BC) polygenic risk score (PRS313) and contralateral breast cancer (CBC) risk, in BRCA1 and BRCA2 pathogenic variant heterozygotes. METHODS: We included women of European ancestry with a prevalent first primary invasive BC (BRCA1 = 6,591 with 1,402 prevalent CBC cases; BRCA2 = 4,208 with 647 prevalent CBC cases) from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA), a large international retrospective series. Cox regression analysis was performed to assess the association between overall and ER-specific PRS313 and CBC risk. RESULTS: For BRCA1 heterozygotes the estrogen receptor (ER)-negative PRS313 showed the largest association with CBC risk, hazard ratio (HR) per SD = 1.12, 95% confidence interval (CI) (1.06-1.18), C-index = 0.53; for BRCA2 heterozygotes, this was the ER-positive PRS313, HR = 1.15, 95% CI (1.07-1.25), C-index = 0.57. Adjusting for family history, age at diagnosis, treatment, or pathological characteristics for the first BC did not change association effect sizes. For women developing first BC < age 40 years, the cumulative PRS313 5th and 95th percentile 10-year CBC risks were 22% and 32% for BRCA1 and 13% and 23% for BRCA2 heterozygotes, respectively. CONCLUSION: The PRS313 can be used to refine individual CBC risks for BRCA1/2 heterozygotes of European ancestry, however the PRS313 needs to be considered in the context of a multifactorial risk model to evaluate whether it might influence clinical decision-making.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status.

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes.

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Association of Genomic Domains in BRCA1 and BRCA2 with Prostate Cancer Risk and Aggressiveness.

Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3' region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25-2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63-5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71-4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12-11.54; P = 0.0002). No genotype-phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. SIGNIFICANCE: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.

The BRCA1 c.4096+3A>G Variant Displays Classical Characteristics of Pathogenic BRCA1 Mutations in Hereditary Breast and Ovarian Cancers, But Still Allows Homozygous Viability.

Mutations in BRCA1 result in predisposal to breast and ovarian cancers, but many variants exist with unknown clinical significance (VUS). One is BRCA1 c.4096+3A>G, which affects production of the full-length BRCA1 transcript, while augmenting transcripts lacking most or all of exon 11. Nonetheless, homozygosity of this variant has been reported in a healthy woman. We saw this variant cosegregate with breast and ovarian cancer in several family branches of four Icelandic pedigrees, with instances of phenocopies and a homozygous woman with lung cancer. We found eight heterozygous carriers (0.44%) in 1820 unselected breast cancer cases, and three (0.15%) in 1968 controls (p = 0.13). Seeking conclusive evidence, we studied tumors from carriers in the pedigrees for wild-type-loss of heterozygosity (wtLOH) and BRCA1-characteristic prevalence of estrogen receptor (ER) negativity. Of 15 breast and six ovarian tumors, wtLOH occurred in nine breast and all six ovarian tumours, and six of the nine breast tumors with wtLOH were ER-negative. These data accord with a pathogenic BRCA1-mutation. Our findings add to the current knowledge of BRCA1, and the role of its exon 11 in cancer pathogenicity, and will be of use in clinical genetic counselling.

High expression of the vacuole membrane protein 1 (VMP1) is a potential marker of poor prognosis in HER2 positive breast cancer.

BACKGROUND: Fusion genes result from genomic structural changes, which can lead to alterations in gene expression that supports tumor development. The aim of the study was to use fusion genes as a tool to identify new breast cancer (BC) genes with a role in BC progression. METHODS: Fusion genes from breast tumors and BC cell lines were collected from publications. RNA-Seq data from tumors and cell lines were retrieved from databanks and analyzed for fusions with SOAPfuse or the analysis was purchased. Fusion genes identified in both tumors (n = 1724) and cell lines (n = 45) were confirmed by qRT-PCR and sequencing. Their individual genes were ranked by selection criteria that included correlation of their mRNA level with copy number. The expression of the top ranked gene was measured by qRT-PCR in normal tissue and in breast tumors from an exploratory cohort (n = 141) and a validation cohort (n = 277). Expression levels were correlated with clinical and pathological factors as well as the patients' survival. The results were followed up in BC cohorts from TCGA (n = 818) and METABRIC (n = 2509). RESULTS: Vacuole membrane protein 1 (VMP1) was the most promising candidate based on specific selection criteria. Its expression was higher in breast tumor tissue than normal tissue (p = 1x10-4), and its expression was significantly higher in HER2 positive than HER2 negative breast tumors in all four cohorts analyzed. High expression of VMP1 associated with breast cancer specific survival (BCSS) in cohort 1 (hazard ratio (HR) = 2.31, CI 1.27-4.18) and METABRIC (HR = 1.26, CI 1.02-1.57), and also after adjusting for HER2 expression in cohort 1 (HR = 2.03, CI 1.10-3.72). BCSS was not significant in cohort 2 or TCGA cohort, which may be due to differences in treatment regimens. CONCLUSIONS: The results suggest that high VMP1 expression is a potential marker of poor prognosis in HER2 positive BC. Further studies are needed to elucidate how VMP1 could affect pathways supportive of tumorigenesis.

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer.

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Mutational spectrum in a worldwide study of 29,700 families with BRCA1 or BRCA2 mutations.

The prevalence and spectrum of germline mutations in BRCA1 and BRCA2 have been reported in single populations, with the majority of reports focused on White in Europe and North America. The Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) has assembled data on 18,435 families with BRCA1 mutations and 11,351 families with BRCA2 mutations ascertained from 69 centers in 49 countries on six continents. This study comprehensively describes the characteristics of the 1,650 unique BRCA1 and 1,731 unique BRCA2 deleterious (disease-associated) mutations identified in the CIMBA database. We observed substantial variation in mutation type and frequency by geographical region and race/ethnicity. In addition to known founder mutations, mutations of relatively high frequency were identified in specific racial/ethnic or geographic groups that may reflect founder mutations and which could be used in targeted (panel) first pass genotyping for specific populations. Knowledge of the population-specific mutational spectrum in BRCA1 and BRCA2 could inform efficient strategies for genetic testing and may justify a more broad-based oncogenetic testing in some populations.

Association of breast cancer risk in BRCA1 and BRCA2 mutation carriers with genetic variants showing differential allelic expression: identification of a modifier of breast cancer risk at locus 11q22.3.

PURPOSE: Cis-acting regulatory SNPs resulting in differential allelic expression (DAE) may, in part, explain the underlying phenotypic variation associated with many complex diseases. To investigate whether common variants associated with DAE were involved in breast cancer susceptibility among BRCA1 and BRCA2 mutation carriers, a list of 175 genes was developed based of their involvement in cancer-related pathways. METHODS: Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of ~320 SNPs located in the vicinity of these genes with breast and ovarian cancer risks in 15,252 BRCA1 and 8211 BRCA2 mutation carriers ascertained from 54 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2. RESULTS: We identified a region on 11q22.3 that is significantly associated with breast cancer risk in BRCA1 mutation carriers (most significant SNP rs228595 p = 7 × 10-6). This association was absent in BRCA2 carriers (p = 0.57). The 11q22.3 region notably encompasses genes such as ACAT1, NPAT, and ATM. Expression quantitative trait loci associations were observed in both normal breast and tumors across this region, namely for ACAT1, ATM, and other genes. In silico analysis revealed some overlap between top risk-associated SNPs and relevant biological features in mammary cell data, which suggests potential functional significance. CONCLUSION: We identified 11q22.3 as a new modifier locus in BRCA1 carriers. Replication in larger studies using estrogen receptor (ER)-negative or triple-negative (i.e., ER-, progesterone receptor-, and HER2-negative) cases could therefore be helpful to confirm the association of this locus with breast cancer risk.

Oestrogen receptor status, treatment and breast cancer prognosis in Icelandic BRCA2 mutation carriers.

BACKGROUND: The impact of an inherited BRCA2 mutation on the prognosis of women with breast cancer has not been well documented. We studied the effects of oestrogen receptor (ER) status, other prognostic factors and treatments on survival in a large cohort of BRCA2 mutation carriers. METHODS: We identified 285 breast cancer patients with a 999del5 BRCA2 mutation and matched them with 570 non-carrier patients. Clinical information was abstracted from patient charts and pathology records and supplemented by evaluation of tumour grade and ER status using archived tissue specimens. Univariate and multivariate hazard ratios (HR) were estimated for breast cancer-specific survival using Cox regression. The effects of various therapies were studied in patients treated from 1980 to 2012. RESULTS: Among mutation carriers, positive ER status was associated with higher risk of death than negative ER status (HR=1.94; 95% CI=1.22-3.07, P=0.005). The reverse association was seen for non-carriers (HR=0.71; 95% CI: 0.51-0.97; P=0.03). CONCLUSIONS: Among BRCA2 carriers, ER-positive status is an adverse prognostic factor. BRCA2 carrier status should be known at the time when treatment decisions are made.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus.

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.

High expression of ZNF703 independent of amplification indicates worse prognosis in patients with luminal B breast cancer.

Amplification of 8p12-p11 is relatively common in breast cancer and several genes within the region have been suggested to affect breast tumor progression. The aim of the study was to map the amplified 8p12-p11 region in a large set of breast tumors in an effort to identify the genetic driver and to explore its impact on tumor progression and prognosis. Copy number alterations (CNAs) were mapped in 359 tumors, and gene expression data from 577 tumors (359 tumors included) were correlated with CNA, clinical-pathological factors, and protein expression (39 tumors). 8p12-p11 was amplified in 11.4% of tumors. The smallest region of amplification harbored one full-length gene, ZNF703. ZNF703 mRNA expression was significantly higher in estrogen receptor (ER)-positive than ER-negative tumors (P = 2 × 10(-16)), a reflection of high expression in luminal tumors. Forty-eight percent of tumors with ZNF703 amplification were luminal B tumors in which the best correlation between DNA copy number and mRNA was seen (P = 1.2 × 10(-7)) as well as correlation between mRNA and protein expression (P = 0.02). High ZNF703 mRNA correlated with poor survival in patients with ER-positive luminal B tumors (log rank P = 0.04). Furthermore, high ZNF703 mRNA expression correlated with poor outcome in patients with ZNF703 copy number neutral, ER-positive, luminal B tumors (log rank P = 0.004). The results support ZNF703 as the driver gene of the 8p12 amplification and suggest that independent of amplification, high expression of the gene affects prognosis in luminal B tumors.

Identification of a BRCA2-specific modifier locus at 6p24 related to breast cancer risk.

Common genetic variants contribute to the observed variation in breast cancer risk for BRCA2 mutation carriers; those known to date have all been found through population-based genome-wide association studies (GWAS). To comprehensively identify breast cancer risk modifying loci for BRCA2 mutation carriers, we conducted a deep replication of an ongoing GWAS discovery study. Using the ranked P-values of the breast cancer associations with the imputed genotype of 1.4 M SNPs, 19,029 SNPs were selected and designed for inclusion on a custom Illumina array that included a total of 211,155 SNPs as part of a multi-consortial project. DNA samples from 3,881 breast cancer affected and 4,330 unaffected BRCA2 mutation carriers from 47 studies belonging to the Consortium of Investigators of Modifiers of BRCA1/2 were genotyped and available for analysis. We replicated previously reported breast cancer susceptibility alleles in these BRCA2 mutation carriers and for several regions (including FGFR2, MAP3K1, CDKN2A/B, and PTHLH) identified SNPs that have stronger evidence of association than those previously published. We also identified a novel susceptibility allele at 6p24 that was inversely associated with risk in BRCA2 mutation carriers (rs9348512; per allele HR = 0.85, 95% CI 0.80-0.90, P = 3.9 × 10(-8)). This SNP was not associated with breast cancer risk either in the general population or in BRCA1 mutation carriers. The locus lies within a region containing TFAP2A, which encodes a transcriptional activation protein that interacts with several tumor suppressor genes. This report identifies the first breast cancer risk locus specific to a BRCA2 mutation background. This comprehensive update of novel and previously reported breast cancer susceptibility loci contributes to the establishment of a panel of SNPs that modify breast cancer risk in BRCA2 mutation carriers. This panel may have clinical utility for women with BRCA2 mutations weighing options for medical prevention of breast cancer.

The risk allele of SNP rs3803662 and the mRNA level of its closest genes TOX3 and LOC643714 predict adverse outcome for breast cancer patients.

BACKGROUND: The minor allele of SNP rs3803662 has been shown to correlate with increased breast cancer risk and with lower expression of TOX3. The SNP is closely located to TOX3 residing within an uncharacterised gene LOC643714. The aim of the study was to examine the association of the risk allele with expression of TOX3 and LOC643714, and of mRNA levels and genotype with clinical and pathological characteristics. METHODS: The SNP was genotyped in DNA isolated from blood and normal tissue from 160 breast cancer patients and mRNA levels were measured by microarrays and quantitative real-time (qRT)-PCR in breast tumours. Association with clinical and pathological characteristics was analysed by parametric tests. RESULTS: An association of the risk allele of rs3803662 with lower TOX3 expression was confirmed in oestrogen receptor (ER) positive tumours. It was more often observed in lobular tumours (p = 0.04), and carriers of the risk allele who had been diagnosed with luminal A tumours had shorter overall survival (OS) than carriers of the non-risk allele (p = 0.01). Positive correlation between the mRNA levels of TOX3 and LOC643714 was observed (r = 0.44 and p < 0.001). Association analysis with tumour pathology showed that low TOX3 and LOC643714 expression correlated with high Ki67 levels (p = 0.026 and p = 0.002) and the basal subtype (p < 0.001 and p < 0.001), whereas high expression correlated with ER (p = 0.004 and p < 0.001) and progesterone receptor (PgR) (p = 0.005 and p < 0.001) expression. Furthermore, high TOX3 and LOC643714 correlated with positive lymph nodes (p < 0.001 and p = 0.01). Patients with ER positive tumours and high levels of TOX3 mRNA had shorter overall- and distant metastasis free-survival (p = 0.017 and p = 0.021), an effect mostly attributable to patients with luminal B tumours. CONCLUSIONS: The results suggest that the effect of the risk allele of rs3803662 is strongest in luminal A tumours and that the expression levels of TOX3 and/or LOC643714 affect the progression of breast cancer. The effect may vary depending on the subtype and developmental stage of the tumour.