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.

Structural variation (SV) represents a significant, yet poorly understood contribution to an individual's genetic makeup. Advanced next-generation sequencing technologies are widely used to discover such variations, but there is no single detection tool that is considered a community standard. In an attempt to fulfil this need, we developed an algorithm, SoftSearch, for discovering structural variant breakpoints in Illumina paired-end next-generation sequencing data. SoftSearch combines multiple strategies for detecting SV including split-read, discordant read-pair, and unmated pairs. Co-localized split-reads and discordant read pairs are used to refine the breakpoints.

ATP Binding Cassette B1 (ABCB1) is a transporter with a broad substrate specificity involved in the elimination of several carcinogens from the gut. Several polymorphic variants within the ABCB1 gene have been reported as modulators of ABCB1-mediated transport. We investigated the impact of ABCB1 genetic variants on colorectal cancer (CRC) risk. A hybrid tagging/functional approach was performed to select 28 single nucleotide polymorphisms (SNPs) that were genotyped in 1,321 Czech subjects, 699 CRC cases and 622 controls. In addition, six potentially functional SNPs were genotyped in 3,662 German subjects, 1,809 cases and 1,853 controls from the DACHS study. We found that three functional SNPs (rs1202168, rs1045642 and rs868755) were associated with CRC risk in the German population. Carriers of the rs1202168_T and rs868755_T alleles had an increased risk for CRC (P(trend) = 0.016 and 0.029, respectively), while individuals bearing the rs1045642_C allele showed a decreased risk of CRC (P(trend) = 0.022). We sought to replicate the most significant results in an independent case-control study of 3,803 subjects, 2,169 cases and 1,634 controls carried out in the North of Germany. None of the SNPs tested were significantly associated with CRC risk in the replication study. In conclusion, in this study of about 8,800 individuals we show that ABCB1 gene polymorphisms play at best a minor role in the susceptibility to CRC.

Nod-like receptors (NLRs) are important innate pattern recognition receptors and regulators of inflammation or play a role during development. We systematically analysed 41 non-synonymous single nucleotide polymorphisms (SNPs) in 21 NLR genes in a Czech discovery cohort of sporadic colorectal cancer (CRC) (1237 cases, 787 controls) for their association with CRC risk and survival. Five SNPs were found to be associated with CRC risk and eight with survival at 5% significance level. In a replication analysis using data of two large genome-wide association studies (GWASs) from Germany (DACHS: 1798 cases and 1810 controls) and Scotland (2210 cases and 9350 controls) the associations found in the Czech discovery set were not confirmed. However, expression analysis in human gut-related tissues and immune cells revealed that the NLRs associated with CRC risk or survival in the discovery set were expressed in primary human colon or rectum cells, CRC tissue and/or cell lines, providing preliminary evidence for a potential involvement of NLRs in general in CRC development and/or progression. Most interesting was the finding that the enigmatic development-related NLRP5 (also known as MATER) was not expressed in normal colon tissue but in colon cancer tissue and cell lines. Future studies may show whether regulatory variants instead of coding variants might affect the expression of NLRs and contribute to CRC risk and survival.

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

In the last four years, Genome-Wide Association Studies (GWAS) have identified sixteen low-penetrance polymorphisms on fourteen different loci associated with colorectal cancer (CRC). Due to the low risks conferred by known common variants, most of the 35% broad-sense heritability estimated by twin studies remains unexplained. Recently our group performed a case-control study for eight Single Nucleotide Polymorphisms (SNPs) in 4 CRC genes. The present investigation is a follow-up of that study. We have genotyped six SNPs that showed a positive association and carried out a meta-analysis based on eight additional studies comprising in total more than 8000 cases and 6000 controls. The estimated recessive odds ratio for one of the SNPs, rs3219489 (MUTYH Q338H), decreased from 1.52 in the original Swedish study, to 1.18 in the Swedish replication, and to 1.08 in the initial meta-analysis. Since the corresponding summary probability value was 0.06, we decided to retrieve additional information for this polymorphism. The incorporation of six further studies resulted in around 13000 cases and 13000 controls. The newly updated OR was 1.03. The results from the present large, multicenter study illustrate the possibility of decreasing effect sizes with increasing samples sizes. Phenotypic heterogeneity, differential environmental exposures, and population specific linkage disequilibrium patterns may explain the observed difference of genetic effects between Sweden and the other investigated cohorts.

Basal-like and triple negative breast cancer (TNBC) share common molecular features, poor prognosis and a propensity for metastasis to the brain. Amplification of epidermal growth factor receptor (EGFR) occurs in ~50% of basal-like breast cancer, and mutations in the epidermal growth factor receptor (EGFR) have been reported in up to ~ 10% of Asian TNBC patients. In non-small cell lung cancer several different mutations in the EGFR tyrosine kinase domain confer sensitivity to receptor tyrosine kinase inhibitors, but the tumourigenic potential of EGFR mutations in breast cells and their potential for targeted therapy is unknown.

We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk.

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.

Defective cellular transport processes can lead to aberrant accumulation of trace elements, iron, small molecules and hormones in the cell, which in turn may promote the formation of reactive oxygen species, promoting DNA damage and aberrant expression of key regulatory cancer genes. As DNA damage and uncontrolled proliferation are hallmarks of cancer, including epithelial ovarian cancer (EOC), we hypothesized that inherited variation in the cellular transport genes contributes to EOC risk.

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.

Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) associated with colorectal cancer risk. These SNPs may tag correlated variants with biological importance. Fine-mapping around GWAS loci can facilitate detection of functional candidates and additional independent risk variants. We analyzed 11,900 cases and 14,311 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. To fine-map genomic regions containing all known common risk variants, we imputed high-density genetic data from the 1000 Genomes Project. We tested single-variant associations with colorectal tumor risk for all variants spanning genomic regions 250-kb upstream or downstream of 31 GWAS-identified SNPs (index SNPs). We queried the University of California, Santa Cruz Genome Browser to examine evidence for biological function. Index SNPs did not show the strongest association signals with colorectal tumor risk in their respective genomic regions. Bioinformatics analysis of SNPs showing smaller P-values in each region revealed 21 functional candidates in 12 loci (5q31.1, 8q24, 11q13.4, 11q23, 12p13.32, 12q24.21, 14q22.2, 15q13, 18q21, 19q13.1, 20p12.3, and 20q13.33). We did not observe evidence of additional independent association signals in GWAS-identified regions. Our results support the utility of integrating data from comprehensive fine-mapping with expanding publicly available genomic databases to help clarify GWAS associations and identify functional candidates that warrant more onerous laboratory follow-up. Such efforts may aid the eventual discovery of disease-causing variant(s).

Genome-wide association studies (GWAS) have successfully identified a number of single-nucleotide polymorphisms (SNPs) associated with colorectal cancer (CRC) risk. However, these susceptibility loci known today explain only a small fraction of the genetic risk. Gene-gene interaction (GxG) is considered to be one source of the missing heritability. To address this, we performed a genome-wide search for pair-wise GxG associated with CRC risk using 8,380 cases and 10,558 controls in the discovery phase and 2,527 cases and 2,658 controls in the replication phase. We developed a simple, but powerful method for testing interaction, which we term the Average Risk Due to Interaction (ARDI). With this method, we conducted a genome-wide search to identify SNPs showing evidence for GxG with previously identified CRC susceptibility loci from 14 independent regions. We also conducted a genome-wide search for GxG using the marginal association screening and examining interaction among SNPs that pass the screening threshold (p<10(-4)). For the known locus rs10795668 (10p14), we found an interacting SNP rs367615 (5q21) with replication p = 0.01 and combined p = 4.19×10(-8). Among the top marginal SNPs after LD pruning (n = 163), we identified an interaction between rs1571218 (20p12.3) and rs10879357 (12q21.1) (nominal combined p = 2.51×10(-6); Bonferroni adjusted p = 0.03). Our study represents the first comprehensive search for GxG in CRC, and our results may provide new insight into the genetic etiology of CRC.

During the past decade, the age-adjusted mortality rate for endometrial cancer (EC) increased 1.9% annually with TP53 mutant (TP53-mu) EC disproportionally represented in advanced disease and deaths. Therefore, we aimed to assess pivotal molecular parameters that differentiate clinical outcomes in high- and low-risk EC. Using the Cancer Genome Atlas, we analyzed EC specimens with available DNA sequences and quantitative gene-specific RNA expression data. After polymerase ɛ (POLE)-mutant specimens were excluded, differential gene-specific mutations and mRNA expressions were annotated and integrated. Consequent to TP53-mu failure to induce p21, derepression of multiple oncogenes harboring promoter p21 repressive sites was observed, including CCNA2 and FOXM1 (P < .001 compared with TP53 wild type [TP53-wt]). TP53-wt EC with high CCNA2 expression (CCNA2-H) had a targeted transcriptomic profile similar to that of TP53-mu EC, suggesting CCNA2 is a seminal determinant for both TP53-wt and TP53-mu EC. CCNA2 enhances E2F1 function, upregulating FOXM1 and CIP2A, as observed in TP53-mu and CCNA2-H TP53-wt EC (P < .001). CIP2A inhibits protein phosphatase 2A, leading to AKT inactivation of GSK3β and restricted oncoprotein degradation; PPP2R1A and FBXW7 mutations yield similar results. Upregulation of FOXM1 and failed degradation of FOXM1 is evidenced by marked upregulation of multiple homologous recombination genes (P < .001). Integrating these molecular aberrations generated a molecular biomarker panel with significant prognostic discrimination (P = 5.8×10-7); adjusting for age, histology, grade, myometrial invasion, TP53 status, and stage, only CCNA2-H/E2F1-H (P = .0003), FBXW7-mu/PPP2R1A-mu (P = .0002), and stage (P = .017) were significant. The generated prognostic molecular classification system identifies dissimilar signaling aberrations potentially amenable to targetable therapeutic options.

The Cancer Genetic Markers of Susceptibility genome-wide association study (GWAS) originally identified a single nucleotide polymorphism (SNP) rs11249433 at 1p11.2 associated with breast cancer risk. To fine-map this locus, we genotyped 92 SNPs in a 900kb region (120,505,799-121,481,132) flanking rs11249433 in 45,276 breast cancer cases and 48,998 controls of European, Asian and African ancestry from 50 studies in the Breast Cancer Association Consortium. Genotyping was done using iCOGS, a custom-built array. Due to the complicated nature of the region on chr1p11.2: 120,300,000-120,505,798, that lies near the centromere and contains seven duplicated genomic segments, we restricted analyses to 429 SNPs excluding the duplicated regions (42 genotyped and 387 imputed). Per-allelic associations with breast cancer risk were estimated using logistic regression models adjusting for study and ancestry-specific principal components. The strongest association observed was with the original identified index SNP rs11249433 (minor allele frequency (MAF) 0.402; per-allele odds ratio (OR) = 1.10, 95% confidence interval (CI) 1.08-1.13, P = 1.49 x 10-21). The association for rs11249433 was limited to ER-positive breast cancers (test for heterogeneity P≤8.41 x 10-5). Additional analyses by other tumor characteristics showed stronger associations with moderately/well differentiated tumors and tumors of lobular histology. Although no significant eQTL associations were observed, in silico analyses showed that rs11249433 was located in a region that is likely a weak enhancer/promoter. Fine-mapping analysis of the 1p11.2 breast cancer susceptibility locus confirms this region to be limited to risk to cancers that are ER-positive.

Ovarian carcinosarcoma is a rare subtype of ovarian cancer with poor clinical outcomes. The low incidence of this disease makes accrual to large clinical trials challenging. However, studies have shown that treatment responses in patient-derived xenograft (PDX) models correlate with matched-patient responses in the clinic, supporting their use for preclinical testing of standard and novel therapies. An ovarian carcinosarcoma PDX is presented herein and showed resistance to carboplatin and paclitaxel (similar to the patient) but exhibited significant sensitivity to ifosfamide and paclitaxel. The PDX demonstrated overexpression of EGFR mRNA and gene amplification by array comparative genomic hybridization (log2 ratio 0.399). EGFR phosphorylation was also detected. Angiogensis and insulin-like growth factor pathways were also implicated by overexpression of VEGFC and IRS1. In order to improve response to chemotherapy, the PDX was treated with carboplatin/paclitaxel with or without a pan-HER and VEGF inhibitor (BMS-690514) but there was no tumor growth inhibition or improved animal survival, which may be explained by a KRAS mutation. Resistance was also observed when the IGF-1R inhibitor BMS-754807 was combined with carboplatin/paclitaxel. Because poly (ADP-ribose) polymerase inhibitors have activity in ovarian cancer patients, with and without BRCA mutations, ABT-888 was also tested but found to have no activity. Pathogenic mutations were also detected in TP53 and PIK3CA. In conclusion, ifosfamide/paclitaxel was superior to carboplatin/paclitaxel in this ovarian carcinosarcoma PDX and gene overexpression or amplification alone was not sufficient to predict response to targeted therapy. Better predictive markers of response are needed.

Molecular events leading to epithelial ovarian cancer are poorly understood but ovulatory hormones and a high number of life-time ovulations with concomitant proliferation, apoptosis, and inflammation, increases risk. We identified genes that are regulated during the estrous cycle in murine ovarian surface epithelium and analysed these profiles to identify genes dysregulated in human ovarian cancer, using publically available datasets. We identified 338 genes that are regulated in murine ovarian surface epithelium during the estrous cycle and dysregulated in ovarian cancer. Six of seven candidates selected for immunohistochemical validation were expressed in serous ovarian cancer, inclusion cysts, ovarian surface epithelium and in fallopian tube epithelium. Most were overexpressed in ovarian cancer compared with ovarian surface epithelium and/or inclusion cysts (EpCAM, EZH2, BIRC5) although BIRC5 and EZH2 were expressed as highly in fallopian tube epithelium as in ovarian cancer. We prioritised the 338 genes for those likely to be important for ovarian cancer development by in silico analyses of copy number aberration and mutation using publically available datasets and identified genes with established roles in ovarian cancer as well as novel genes for which we have evidence for involvement in ovarian cancer. Chromosome segregation emerged as an important process in which genes from our list of 338 were over-represented including two (BUB1, NCAPD2) for which there is evidence of amplification and mutation. NUAK2, upregulated in ovarian surface epithelium in proestrus and predicted to have a driver mutation in ovarian cancer, was examined in a larger cohort of serous ovarian cancer where patients with lower NUAK2 expression had shorter overall survival. In conclusion, defining genes that are activated in normal epithelium in the course of ovulation that are also dysregulated in cancer has identified a number of pathways and novel candidate genes that may contribute to the development of ovarian cancer.

Approximately 15-30% of all breast cancer tumors are estrogen receptor negative (ER-). Compared with ER-positive (ER+) disease they have an earlier age at onset and worse prognosis. Despite the vast number of risk variants identified for numerous cancer types, only seven loci have been unambiguously identified for ER-negative breast cancer. With the aim of identifying new susceptibility SNPs for this disease we performed a pleiotropic genome-wide association study (GWAS). We selected 3079 SNPs associated with a human complex trait or disease at genome-wide significance level (P<5 × 10(-8)) to perform a secondary analysis of an ER-negative GWAS from the National Cancer Institute's Breast and Prostate Cancer Cohort Consortium (BPC3), including 1998 cases and 2305 controls from prospective studies. We then tested the top ten associations (i.e. with the lowest P-values) using three additional populations with a total sample size of 3509 ER+ cases, 2543 ER- cases and 7031 healthy controls. None of the 3079 selected variants in the BPC3 ER-GWAS were significant at the adjusted threshold. 186 variants were associated with ER- breast cancer risk at a conventional threshold of P<0.05, with P-values ranging from 0.049 to 2.3 × 10(-4). None of the variants reached statistical significance in the replication phase. In conclusion, this study did not identify any novel susceptibility loci for ER-breast cancer using a "pleiotropic approach".

Common variation on 14q24.1, close to RAD51B, has been associated with breast cancer: rs999737 and rs2588809 with the risk of female breast cancer and rs1314913 with the risk of male breast cancer. The aim of this study was to investigate the role of RAD51B variants in breast cancer predisposition, particularly in the context of familial breast cancer in Finland. We sequenced the coding region of RAD51B in 168 Finnish breast cancer patients from the Helsinki region for identification of possible recurrent founder mutations. In addition, we studied the known rs999737, rs2588809, and rs1314913 SNPs and RAD51B haplotypes in 44,791 breast cancer cases and 43,583 controls from 40 studies participating in the Breast Cancer Association Consortium (BCAC) that were genotyped on a custom chip (iCOGS). We identified one putatively pathogenic missense mutation c.541C>T among the Finnish cancer patients and subsequently genotyped the mutation in additional breast cancer cases (n = 5259) and population controls (n = 3586) from Finland and Belarus. No significant association with breast cancer risk was seen in the meta-analysis of the Finnish datasets or in the large BCAC dataset. The association with previously identified risk variants rs999737, rs2588809, and rs1314913 was replicated among all breast cancer cases and also among familial cases in the BCAC dataset. The most significant association was observed for the haplotype carrying the risk-alleles of all the three SNPs both among all cases (odds ratio (OR): 1.15, 95% confidence interval (CI): 1.11-1.19, P = 8.88 x 10-16) and among familial cases (OR: 1.24, 95% CI: 1.16-1.32, P = 6.19 x 10-11), compared to the haplotype with the respective protective alleles. Our results suggest that loss-of-function mutations in RAD51B are rare, but common variation at the RAD51B region is significantly associated with familial breast cancer risk.