Breast cancer (BC) is a heterogeneous disease, and different tumor characteristics and genetic variation may affect the clinical outcome. The FANCM c.5101C > T nonsense mutation in the Finnish population associates with increased risk of breast cancer, especially for triple-negative breast cancer patients. To investigate the association of the mutation with disease prognosis, we studied tumor phenotype, treatment outcome, and patient survival in 3,933 invasive breast cancer patients, including 101 FANCM c.5101C > T mutation carriers and 3,832 non-carriers. We also examined association of the mutation with nuclear immunohistochemical staining of DNA repair markers in 1,240 breast tumors. The FANCM c.5101C > T mutation associated with poor 10-year breast cancer-specific survival (hazard ratio (HR)=1.66, 95% confidence interval (CI) 1.09-2.52, p = 0.018), with a more pronounced survival effect among familial cases (HR = 2.93, 95% CI 1.5-5.76, p = 1.80 × 10-3 ). Poor disease outcome of the carriers was also found among the estrogen receptor (ER) positive subgroup of patients (HR = 1.8, 95% CI 1.09-2.98, p = 0.021). Reduced survival was seen especially among patients who had not received radiotherapy (HR = 3.43, 95% CI 1.6-7.34, p = 1.50 × 10-3 ) but not among radiotherapy treated patients (HR = 1.35, 95% CI 0.82-2.23, p = 0.237). Significant interaction was found between the mutation and radiotherapy (p = 0.040). Immunohistochemical analyses show that c.5101C > T carriers have reduced PAR-activity. Our results suggest that FANCM c.5101C > T nonsense mutation carriers have a reduced breast cancer survival but postoperative radiotherapy may diminish this survival disadvantage.

Genome-wide association studies have found SNPs at 17q22 to be associated with breast cancer risk. To identify potential causal variants related to breast cancer risk, we performed a high resolution fine-mapping analysis that involved genotyping 517 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of genotypes for 3,134 SNPs in more than 89,000 participants of European ancestry from the Breast Cancer Association Consortium (BCAC). We identified 28 highly correlated common variants, in a 53 Kb region spanning two introns of the STXBP4 gene, that are strong candidates for driving breast cancer risk (lead SNP rs2787486 (OR = 0.92; CI 0.90-0.94; P = 8.96 × 10(-15))) and are correlated with two previously reported risk-associated variants at this locus, SNPs rs6504950 (OR = 0.94, P = 2.04 × 10(-09), r(2) = 0.73 with lead SNP) and rs1156287 (OR = 0.93, P = 3.41 × 10(-11), r(2) = 0.83 with lead SNP). Analyses indicate only one causal SNP in the region and several enhancer elements targeting STXBP4 are located within the 53 kb association signal. Expression studies in breast tumor tissues found SNP rs2787486 to be associated with increased STXBP4 expression, suggesting this may be a target gene of this locus.

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Genome-wide association studies (GWASs) have revealed increased breast cancer risk associated with multiple genetic variants at 5p12. Here, we report the fine mapping of this locus using data from 104,660 subjects from 50 case-control studies in the Breast Cancer Association Consortium (BCAC). With data for 3,365 genotyped and imputed SNPs across a 1 Mb region (positions 44,394,495-45,364,167; NCBI build 37), we found evidence for at least three independent signals: the strongest signal, consisting of a single SNP rs10941679, was associated with risk of estrogen-receptor-positive (ER+) breast cancer (per-g allele OR ER+ = 1.15; 95% CI 1.13-1.18; p = 8.35 × 10-30). After adjustment for rs10941679, we detected signal 2, consisting of 38 SNPs more strongly associated with ER-negative (ER-) breast cancer (lead SNP rs6864776: per-a allele OR ER- = 1.10; 95% CI 1.05-1.14; p conditional = 1.44 × 10-12), and a single signal 3 SNP (rs200229088: per-t allele OR ER+ = 1.12; 95% CI 1.09-1.15; p conditional = 1.12 × 10-05). Expression quantitative trait locus analysis in normal breast tissues and breast tumors showed that the g (risk) allele of rs10941679 was associated with increased expression of FGF10 and MRPS30. Functional assays demonstrated that SNP rs10941679 maps to an enhancer element that physically interacts with the FGF10 and MRPS30 promoter regions in breast cancer cell lines. FGF10 is an oncogene that binds to FGFR2 and is overexpressed in ∼10% of human breast cancers, whereas MRPS30 plays a key role in apoptosis. These data suggest that the strongest signal of association at 5p12 is mediated through coordinated activation of FGF10 and MRPS30, two candidate genes for breast cancer pathogenesis.

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

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.

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.

Genome wide association studies (GWAs) of breast cancer mortality have identified few potential associations. The concordance between these studies is unclear. In this study, we used a meta-analysis of two prognostic GWAs and a replication cohort to identify the strongest associations and to evaluate the loci suggested in previous studies. We attempt to identify those SNPs which could impact overall survival irrespective of the age of onset.

1α,25-Dihydroxyvitamin D3 (1α,25(OH)2D3) had earlier been regarded as the only active hormone. The newly identified actions of 25-hydroxyvitamin D3 (25(OH)D3) and 24R,25-dihydroxyvitamin D3 (24R,25(OH)2D3) broadened the vitamin D3 endocrine system, however, the current data are fragmented and a systematic understanding is lacking. Here we performed the first systematic study of global gene expression to clarify their similarities and differences. Three metabolites at physiologically comparable levels were utilized to treat human and mouse fibroblasts prior to DNA microarray analyses. Human primary prostate stromal P29SN cells (hP29SN), which convert 25(OH)D3 into 1α,25(OH)2D3 by 1α-hydroxylase (encoded by the gene CYP27B1), displayed regulation of 164, 171, and 175 genes by treatment with 1α,25(OH)2D3, 25(OH)D3, and 24R,25(OH)2D3, respectively. Mouse primary Cyp27b1 knockout fibroblasts (mCyp27b1 (-/-)), which lack 1α-hydroxylation, displayed regulation of 619, 469, and 66 genes using the same respective treatments. The number of shared genes regulated by two metabolites is much lower in hP29SN than in mCyp27b1 (-/-). By using DAVID Functional Annotation Bioinformatics Microarray Analysis tools and Ingenuity Pathways Analysis, we identified the agonistic regulation of calcium homeostasis and bone remodeling between 1α,25(OH)2D3 and 25(OH)D3 and unique non-classical actions of each metabolite in physiological and pathological processes, including cell cycle, keratinocyte differentiation, amyotrophic lateral sclerosis signaling, gene transcription, immunomodulation, epigenetics, cell differentiation, and membrane protein expression. In conclusion, there are three distinct vitamin D3 hormones with clearly different biological activities. This study presents a new conceptual insight into the vitamin D3 endocrine system, which may guide the strategic use of vitamin D3 in disease prevention and treatment.

Purpose BRCA1/2 mutations increase the risk of breast and prostate cancer in men. Common genetic variants modify cancer risks for female carriers of BRCA1/2 mutations. We investigated-for the first time to our knowledge-associations of common genetic variants with breast and prostate cancer risks for male carriers of BRCA1/ 2 mutations and implications for cancer risk prediction. Materials and Methods We genotyped 1,802 male carriers of BRCA1/2 mutations from the Consortium of Investigators of Modifiers of BRCA1/2 by using the custom Illumina OncoArray. We investigated the combined effects of established breast and prostate cancer susceptibility variants on cancer risks for male carriers of BRCA1/2 mutations by constructing weighted polygenic risk scores (PRSs) using published effect estimates as weights. Results In male carriers of BRCA1/2 mutations, PRS that was based on 88 female breast cancer susceptibility variants was associated with breast cancer risk (odds ratio per standard deviation of PRS, 1.36; 95% CI, 1.19 to 1.56; P = 8.6 × 10-6). Similarly, PRS that was based on 103 prostate cancer susceptibility variants was associated with prostate cancer risk (odds ratio per SD of PRS, 1.56; 95% CI, 1.35 to 1.81; P = 3.2 × 10-9). Large differences in absolute cancer risks were observed at the extremes of the PRS distribution. For example, prostate cancer risk by age 80 years at the 5th and 95th percentiles of the PRS varies from 7% to 26% for carriers of BRCA1 mutations and from 19% to 61% for carriers of BRCA2 mutations, respectively. Conclusion PRSs may provide informative cancer risk stratification for male carriers of BRCA1/2 mutations that might enable these men and their physicians to make informed decisions on the type and timing of breast and prostate cancer risk management.

Germline genetic variation has been suggested to influence the survival of breast cancer patients independently of tumor pathology. We have studied survival associations of genetic variants in two etiologically unique groups of breast cancer patients, the carriers of germline pathogenic variants in BRCA1 or BRCA2 genes. We found that rs57025206 was significantly associated with the overall survival, predicting higher mortality of BRCA1 carrier patients with estrogen receptor-negative breast cancer, with a hazard ratio 4.37 (95% confidence interval 3.03-6.30, P = 3.1 × 10-9). Multivariable analysis adjusted for tumor characteristics suggested that rs57025206 was an independent survival marker. In addition, our exploratory analyses suggest that the associations between genetic variants and breast cancer patient survival may depend on tumor biological subgroup and clinical patient characteristics.

Insertions and deletions (indels) in human genomes are associated with a wide range of phenotypes, including various clinical disorders. High-throughput, next generation sequencing (NGS) technologies enable the detection of short genetic variants, such as single nucleotide variants (SNVs) and indels. However, the variant calling accuracy for indels remains considerably lower than for SNVs. Here we present a comparative study of the performance of variant calling tools for indel calling, evaluated with a wide repertoire of NGS datasets. While there is no single optimal tool to suit all circumstances, our results demonstrate that the choice of variant calling tool greatly impacts the precision and recall of indel calling. Furthermore, to reliably detect indels, it is essential to choose NGS technologies that offer a long read length and high coverage coupled with specific variant calling tools.

Basal-like breast cancers (BLBC) are characterized by defects in homologous recombination (HR), deficient mitotic checkpoint, and high-proliferation activity. Here, we discover CIP2A as a candidate driver of BLBC. CIP2A was essential for DNA damage-induced initiation of mouse BLBC-like mammary tumors and for survival of HR-defective BLBC cells. CIP2A was dispensable for normal mammary gland development and for unperturbed mitosis, but selectively essential for mitotic progression of DNA damaged cells. A direct interaction between CIP2A and a DNA repair scaffold protein TopBP1 was identified, and CIP2A inhibition resulted in enhanced DNA damage-induced TopBP1 and RAD51 recruitment to chromatin in mammary epithelial cells. In addition to its role in tumor initiation, and survival of BRCA-deficient cells, CIP2A also drove proliferative MYC and E2F1 signaling in basal-like triple-negative breast cancer (BL-TNBC) cells. Clinically, high CIP2A expression was associated with poor patient prognosis in BL-TNBCs but not in other breast cancer subtypes. Small-molecule reactivators of PP2A (SMAP) inhibited CIP2A transcription, phenocopied the CIP2A-deficient DNA damage response (DDR), and inhibited growth of patient-derived BLBC xenograft. In summary, these results demonstrate that CIP2A directly interacts with TopBP1 and coordinates DNA damage-induced mitotic checkpoint and proliferation, thereby driving BLBC initiation and progression. SMAPs could serve as a surrogate therapeutic strategy to inhibit the oncogenic activity of CIP2A in BLBCs. SIGNIFICANCE: These results identify CIP2A as a nongenetic driver and therapeutic target in basal-like breast cancer that regulates DNA damage-induced G2-M checkpoint and proliferative signaling.

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82-0.88]) and ER-negative (OR = 0.87 [0.82-0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91-0.95] and OR = 1.06 [1.03-1.09]) and ER-negative (OR = 0.95 [0.91-0.98] and OR = 1.08 [1.04-1.13]) disease. There was weaker evidence for iCHAV4, located 5' of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90-0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1-4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.

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.

Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

Patient-derived xenografts (PDXs) are considered to better recapitulate the histopathological and molecular heterogeneity of human cancer than other preclinical models. Despite technological advances, PDX models from hormone naïve primary prostate cancer are scarce. We performed a detailed analysis of PDX methodology using a robust subcutaneous model and fresh tissues from patients with primary hormone naïve prostate cancer.

The sequence contexts of genomic variants play important roles in understanding biological significances of variants and potential sequencing related variant calling issues. However, methods for assessing the diverse sequence contexts of genomic variants such as tandem repeats and unambiguous annotations have been limited. Herein, we describe the Variant Sequence Context Annotation Tool (VarSCAT) for annotating the sequence contexts of genomic variants, including breakpoint ambiguities, flanking bases of variants, wildtype/mutated DNA sequences, variant nomenclatures, distances between adjacent variants, tandem repeat regions, and custom annotation with user customizable options. Our analyses demonstrate that VarSCAT is more versatile and customizable than the currently available methods or strategies for annotating variants in short tandem repeat (STR) regions or insertions and deletions (indels) with breakpoint ambiguity. Variant sequence context annotations of high-confidence human variant sets with VarSCAT revealed that more than 75% of all human individual germline and clinically relevant indels have breakpoint ambiguities. Moreover, we illustrate that more than 80% of human individual germline small variants in STR regions are indels and that the sizes of these indels correlated with STR motif sizes. VarSCAT is available from https://github.com/elolab/VarSCAT.

We have utilized a two-stage study design to search for SNPs associated with the survival of breast cancer patients treated with adjuvant chemotherapy. Our initial GWS data set consisted of 805 Finnish breast cancer cases (360 treated with adjuvant chemotherapy). The top 39 SNPs from this stage were analyzed in three independent data sets: iCOGS (n=6720 chemotherapy-treated cases), SUCCESS-A (n=3596), and POSH (n=518). Two SNPs were successfully validated: rs6500843 (any chemotherapy; per-allele HR 1.16, 95% C.I. 1.08-1.26, p=0.0001, p(adjusted)=0.0091), and rs11155012 (anthracycline therapy; per-allele HR 1.21, 95% C.I. 1.08-1.35, p=0.0010, p(adjusted)=0.0270). The SNP rs6500843 was found to specifically interact with adjuvant chemotherapy, independently of standard prognostic markers (p(interaction)=0.0009), with the rs6500843-GG genotype corresponding to the highest hazard among chemotherapy-treated cases (HR 1.47, 95% C.I. 1.20-1.80). Upon trans-eQTL analysis of public microarray data, the rs6500843 locus was found to associate with the expression of a group of genes involved in cell cycle control, notably AURKA, the expression of which also exhibited differential prognostic value between chemotherapy-treated and untreated cases in our analysis of microarray data. Based on previously published information, we propose that the eQTL genes may be connected to the rs6500843 locus via a RBFOX1-FOXM1 -mediated regulatory pathway.