Women with epithelial ovarian cancer (EOC) are usually treated with platinum/taxane therapy after cytoreductive surgery but there is considerable inter-individual variation in response. To identify germline single-nucleotide polymorphisms (SNPs) that contribute to variations in individual responses to chemotherapy, we carried out a multi-phase genome-wide association study (GWAS) in 1,244 women diagnosed with serous EOC who were treated with the same first-line chemotherapy, carboplatin and paclitaxel. We identified two SNPs (rs7874043 and rs72700653) in TTC39B (best P=7x10-5, HR=1.90, for rs7874043) associated with progression-free survival (PFS). Functional analyses show that both SNPs lie in a putative regulatory element (PRE) that physically interacts with the promoters of PSIP1, CCDC171 and an alternative promoter of TTC39B. The C allele of rs7874043 is associated with poor PFS and showed increased binding of the Sp1 transcription factor, which is critical for chromatin interactions with PSIP1. Silencing of PSIP1 significantly impaired DNA damage-induced Rad51 nuclear foci and reduced cell viability in ovarian cancer lines. PSIP1 (PC4 and SFRS1 Interacting Protein 1) is known to protect cells from stress-induced apoptosis, and high expression is associated with poor PFS in EOC patients. We therefore suggest that the minor allele of rs7874043 confers poor PFS by increasing PSIP1 expression.

Circulating microRNAs have drawn a great deal of attention as promising novel biomarkers for breast cancer. However, to date, the results are mixed. Here, we performed a three-stage microRNA analysis using plasma samples from breast cancer patients and healthy controls, with efforts taken to address several pitfalls in detection techniques and study design observed in previous studies. In the discovery phase with 122 Caucasian study subjects, we identified 43 microRNAs differentially expressed between breast cancer cases and healthy controls. When those microRNAs were compared with published data from other studies, we identified three microRNAs, including miR-148b, miR-133a and miR-409-3p, whose plasma levels were significantly higher in breast cancer cases than healthy controls and were also significant in previous independent studies. In the validation phase with 50 breast cancer cases and 50 healthy controls, we validated the associations with breast cancer detection for miR-148b and miR-133a (P = 1.5×10-6 and 1.3×10-10, respectively). In the in-vitro study phase, we found that both miR-148b and miR-133a were secreted from breast cancer cell lines, showing their secretory potential and possible tumor origin. Thus, our data suggest that both miR-148b and miR-133a have potential use as biomarkers for breast cancer detection.

Most non-BRCA1/2 breast cancer families have no identified genetic cause. We used linkage and haplotype analyses in familial and sporadic breast cancer cases to identify a susceptibility locus on chromosome 6q. Two independent genome-wide linkage analysis studies suggested a 3 Mb locus on chromosome 6q and two unrelated Swedish families with a LOD >2 together seemed to share a haplotype in 6q14.1. We hypothesized that this region harbored a rare high-risk founder allele contributing to breast cancer in these two families. Sequencing of DNA and RNA from the two families did not detect any pathogenic mutations. Finally, 29 SNPs in the region were analyzed in 44,214 cases and 43,532 controls from BCAC, and the original haplotypes in the two families were suggested as low-risk alleles for European and Swedish women specifically. There was also some support for one additional independent moderate-risk allele in Swedish familial samples. The results were consistent with our previous findings in familial breast cancer and supported a breast cancer susceptibility locus at 6q14.1 around the PHIP gene.

There are significant inter-individual differences in the levels of gene expression. Through modulation of gene expression, cis-acting variants represent an important source of phenotypic variation. Consequently, cis-regulatory SNPs associated with differential allelic expression are functional candidates for further investigation as disease-causing variants. To investigate whether common variants associated with differential allelic expression were involved in breast cancer susceptibility, a list of genes was established on the basis of their involvement in cancer related pathways and/or mechanisms. Thereafter, using data from a genome-wide map of allelic expression associated SNPs, 313 genetic variants were selected and their association with breast cancer risk was then evaluated in 46,451 breast cancer cases and 42,599 controls of European ancestry ascertained from 41 studies participating in the Breast Cancer Association Consortium. The associations were evaluated with overall breast cancer risk and with estrogen receptor negative and positive disease. One novel breast cancer susceptibility locus on 4q21 (rs11099601) was identified (OR = 1.05, P = 5.6x10-6). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, HELQ, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, MRPS18C encoding the Mitochondrial Ribosomal Protein S18C and FAM175A (ABRAXAS), encoding a BRCA1 BRCT domain-interacting protein involved in DNA damage response and double-strand break (DSB) repair. Expression QTL analysis in breast cancer tissue showed rs11099601 to be associated with HELQ (P = 8.28x10-14), MRPS18C (P = 1.94x10-27) and FAM175A (P = 3.83x10-3), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.

DNA methylation is altered in many types of disease, including metastatic colorectal cancer. However, the methylome has not yet been fully described in archival formalin-fixed paraffin embedded (FFPE) samples in the context of matched fresh-frozen (FF) tumor material at base-pair resolution using a targeted approach. Using next-generation sequencing, we investigated three pairs of matched FFPE and FF samples to determine the extent of their similarity. We identified a 'bowing' pattern specific to FFPE samples categorized by a lower CG proportion at the start of sequence reads. We have found no evidence that this affected methylation calling, nor concordance of results. We also found no significant increase in deamination, measured by C>T transitions, previously considered a result of crosslinking DNA by formalin fixation and a barrier to the use of FFPE in methylation studies. The methods used in this study have shown sensitivity of between 60-70% based on positions also methylated in colorectal cancer cell lines. We demonstrate that FFPE material is a useful source of tumor material for methylation studies using targeted sequencing.

Hematotoxicity is one of the major side effects of chemotherapy. The aim of this study was to examine the association between single nucleotide polymorphisms (SNPs) and hematotoxicity in breast cancer patients in a subset of patients of the SUCCESS prospective phase III chemotherapy study. All patients (n = 1678) received three cycles of 5-fluorouracil, epirubicin, and cyclophosphamide (FEC) followed by three cycles of docetaxel or docetaxel/gemcitabine, depending on randomization. Germline DNA was genotyped for 246 SNPs selected from a previous genome-wide association study (GWAS) in a panel of lymphoblastoid cell lines, with gemcitabine toxicity as the phenotype. All SNPs were tested for their value in predicting grade 3 or 4 neutropenic or leukopenic events (NLEs). Their prognostic value in relation to overall survival and disease-free survival was also tested. None of the SNPs was found to be predictive for NLEs during treatment with docetaxel/gemcitabine. Two SNPs in and close to the PIGB gene significantly improved the prediction of NLEs after FEC, in addition to the factors of age and body surface area. The top SNP (rs12050587) had an odds ratio of 1.38 per minor allele (95% confidence interval, 1.17 to 1.62). No associations were identified for predicting disease-free or overall survival. Genetic variance in the PIGB gene may play a role in determining interindividual differences in relation to hematotoxicity after FEC chemotherapy.

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.

Hepatic angiosarcoma is a rare and aggressive vascular neoplasm. Pathogenic driver mutations are largely unknown. We present the case of a patient with recurrent hepatic angiosarcoma, who initially showed good response to sunitinib, followed by progression. Using comprehensive molecular techniques, we explored the potential mechanisms of resistance. By low-read-depth whole-genome sequencing, the comparison of copy number aberrations (CNAs) of the primary tumor to the skin metastatic lesion that developed after progression on sunitinib, revealed high-level amplification of the 4q11-q13.1 region (containing KIT, PDGFRA and VEGFR2 genes) that was sustained in both lesions. Whole exome sequencing on the germline, primary and metastatic tumor DNAs, resulted in 27 confirmed mutations, 19 of which (including TP53 mutation) presented in both primary and metastatic lesions. One mutation, ZNF331 frameshift deletion, was detected only in the primary tumor. In contrast, seven other mutations, including phospholipase C-gamma1 (PLCG1) R707Q mutation, were found only in the metastatic tumor, indicating selection of cells with the resistant genotype under sunitinib pressure. Our study supports the notion that PLCG1-R707Q mutation may confer VEGFR2-independent signaling and may thus cause resistance against VEGF(R)-directed therapies. This case illustrates also the advantages of using next-generation technologies in identifying individualized targeted therapy.

We analyzed whole exome sequencing data in germline DNA from 412 high grade serous ovarian cancer (HGSOC) cases from The Cancer Genome Atlas Project and identified 5,517 genes harboring a predicted deleterious germline coding mutation in at least one HGSOC case. Gene-set enrichment analysis showed enrichment for genes involved in DNA repair (p = 1.8×10-3). Twelve DNA repair genes - APEX1, APLF, ATX, EME1, FANCL, FANCM, MAD2L2, PARP2, PARP3, POLN, RAD54L and SMUG1 - were prioritized for targeted sequencing in up to 3,107 HGSOC cases, 1,491 cases of other epithelial ovarian cancer (EOC) subtypes and 3,368 unaffected controls of European origin. We estimated mutation prevalence for each gene and tested for associations with disease risk. Mutations were identified in both cases and controls in all genes except MAD2L2, where we found no evidence of mutations in controls. In FANCM we observed a higher mutation frequency in HGSOC cases compared to controls (29/3,107 cases, 0.96 percent; 13/3,368 controls, 0.38 percent; P=0.008) with little evidence for association with other subtypes (6/1,491, 0.40 percent; P=0.82). The relative risk of HGSOC associated with deleterious FANCM mutations was estimated to be 2.5 (95% CI 1.3 - 5.0; P=0.006). In summary, whole exome sequencing of EOC cases with large-scale replication in case-control studies has identified FANCM as a likely novel susceptibility gene for HGSOC, with mutations associated with a moderate increase in risk. These data may have clinical implications for risk prediction and prevention approaches for high-grade serous ovarian cancer in the future and a significant impact on reducing disease mortality.

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance (p<1×10-5). Larger patient numbers will be needed to convincingly identify any true associations at these loci.

RNA editing in mammals is a form of post-transcriptional modification in which adenosine is converted to inosine by the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on evidence of altered ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in ADAR genes modify EOC susceptibility, potentially by altering ovarian tissue gene expression. Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 controls, we evaluated the associations of 5,303 SNPs in ADAD1, ADAR, ADAR2, ADAR3, and SND1. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI), with adjustment for European ancestry. We conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for common and rare variants (SKAT-CR). Association analysis revealed top risk-associated SNP rs77027562 (OR (95% CI)= 1.39 (1.17-1.64), P=1.0x10-4) in ADAR3 and rs185455523 in SND1 (OR (95% CI)= 0.68 (0.56-0.83), P=2.0x10-4). When restricting to serous histology (n=6,500), the magnitude of association strengthened for rs185455523 (OR=0.60, P=1.0x10-4). Gene-level analyses revealed that variation in ADAR was associated (P<0.05) with EOC susceptibility, with PAML=0.022 and PSKAT-CR=0.020. Expression quantitative trait locus analysis in EOC tissue revealed significant associations (P<0.05) with ADAR expression for several SNPs in ADAR, including rs1127313 (G/A), a SNP in the 3' untranslated region. In summary, germline variation involving RNA editing genes may influence EOC susceptibility, warranting further investigation of inherited and acquired alterations affecting RNA editing.

Endometrial carcinoma (EnCa) is the most common invasive gynaecologic carcinoma. Over 85% of EnCa are classified as endometrioid, expressing steroid hormone receptors and mostly involving pathological prestages. Human endogenous retroviruses (ERV) are chromosomally integrated genes, account for about 8% of the human genome and are implicated in the etiology of carcinomas. The majority of ERV envelope (env) coding genes are either not present or not consistently represented between common gene expression microarrays. The aim of this study was to analyse the absolute gene expression of all known 21 ERV env genes including 19 codogenic and two env genes with premature stop codons in EnCa, endometrium as well as in hyperplasia and polyps. For EnCa seven env genes had high expression with >200 mol/ng cDNA (e.g. envH1-3, Syncytin-1, envT), two middle >50 mol/ng cDNA (envFc2, erv-3) and 12 low less than 50 mol/ng cDNA (e.g. Syncytin-2, envV2). Regarding tumor parameters, Syncytin-1 and Syncytin-2 were significantly over-expressed in advanced stage pT2 compared to pT1b. In less differentiated EnCa Syncytin-1, erv-3, envT and envFc2 were significantly over-expressed. Syncytin-1, Syncytin-2 and erv-3 were specific to glandular epithelial cells of polyps, hyperplasia and EnCa using immunohistochemistry. An analysis of 10 patient-matched EnCa with endometrium revealed that the ERV-W 5' long terminal repeat regulating Syncytin-1 was hypomethylated, including the ERE and CRE overlapping MeCP2 sites. Functional analyses showed that 10 env genes were regulated by methylation in EnCa using the RL95-2 cell line. In conclusion, over-expressed env genes could serve as indicators for pathological pre-stages and EnCa.

Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci.

An important precursor to lung cancer development is chronic obstructive pulmonary disease (COPD), independent of exposure to tobacco smoke. Both diseases are associated with increased host susceptibility, inflammation, and genomic instability. However, validation of the candidate genes and functional confirmation to test shared genetic contribution and cellular mechanisms to the development of lung cancer in patients with COPD remains underexplored. Here, we show that loss of PARK2 (encoding Parkin) increases the expression of proinflammation factors as well as nuclear NF-κB localization, suggesting a role of PARK2 loss in inflammation. Additional exploration showed that PARK2 deficiency promotes genomic instability and cell transformation. This role of PARK2 in inflammation and chromosome instability provides a potential link among Parkin, COPD and lung cancer. A further comprehensive validation of 114 informative single nucleotide polymorphism (SNP) variants of PARK2, in 2,484 cases and controls with well-defined lung cancer and COPD phenotypes, found rs577876, rs6455728 and rs9346917 (p<0.01) to be significantly associated with lung cancer development in people with COPD. Our findings support the evidence that PARK2 might have a tumor suppressor role in the development of COPD and lung cancer.

The clinical heterogeneity of prostate cancer (PCa) makes it difficult to identify those patients that could benefit from more aggressive treatments. As a contribution to a better understanding of the genomic changes in the primary tumor that are associated with the development of high-risk disease, we performed exome sequencing and copy number determination of a clinically homogeneous cohort of 47 high-risk PCas. We confirmed recurrent mutations in SPOP, PTEN and TP53 among the 850 point mutations we detected. In seven cases, we discovered genomic aberrations in the TET1 (Ten-Eleven Translocation 1) gene which encodes a DNA hydroxymethylase than can modify methylated cytosines in genomic DNA and thus is linked with gene expression changes. TET1 protein levels were reduced in tumor versus non-tumor prostate tissue in 39 of 40 cases. The clinical relevance of changes in TET1 levels was demonstrated in an independent PCa cohort, in which low TET1 mRNA levels were significantly associated with worse metastases-free survival. We also demonstrate a strong reduction in hydroxymethylated DNA in tumor tissue in 27 of 41 cases. Furthermore, we report the first exploratory (h)MeDIP-Seq analyses of eight high-risk PCa samples. This reveals a large heterogeneity in hydroxymethylation changes in tumor versus non-tumor genomes which can be linked with cell polarity.

Gene fusions play a critical role in some cancers and can serve as important clinical targets. In epithelial ovarian cancer (EOC), the contribution of fusions, especially by histological type, is unclear. We therefore screened for recurrent fusions in a histologically diverse panel of 220 EOCs using RNA sequencing. The Pipeline for RNA-Sequencing Data Analysis (PRADA) was used to identify fusions and allow for comparison with The Cancer Genome Atlas (TCGA) tumors. Associations between fusions and clinical prognosis were evaluated using Cox proportional hazards regression models. Nine recurrent fusions, defined as occurring in two or more tumors, were observed. CRHR1-KANSL1 was the most frequently identified fusion, identified in 6 tumors (2.7% of all tumors). This fusion was not associated with survival; other recurrent fusions were too rare to warrant survival analyses. One recurrent in-frame fusion, UBAP1-TGM7, was unique to clear cell (CC) EOC tumors (in 10%, or 2 of 20 CC tumors). We found some evidence that CC tumors harbor more fusions on average than any other EOC histological type, including high-grade serous (HGS) tumors. CC tumors harbored a mean of 7.4 fusions (standard deviation [sd] = 7.4, N = 20), compared to HGS EOC tumors mean of 2.0 fusions (sd = 3.3, N = 141). Few fusion genes were detected in endometrioid tumors (mean = 0.24, sd = 0.74, N = 55) or mucinous tumors (mean = 0.25, sd = 0.5, N = 4) tumors. To conclude, we identify one fusion at 10% frequency in the CC EOC subtype, but find little evidence for common (> 5% frequency) recurrent fusion genes in EOC overall, or in HGS subtype-specific EOC tumors.

Regulatory T (Treg) cells, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and, thus, variants in genes encoding Treg cell immune molecules could be associated with ovarian cancer.

TP53 overexpression is indicative of somatic TP53 mutations and associates with aggressive tumors and poor prognosis in breast cancer. We utilized a two-stage SNP association study to detect variants associated with breast cancer survival in a TP53-dependent manner. Initially, a genome-wide study (n = 575 cases) was conducted to discover candidate SNPs for genotyping and validation in the Breast Cancer Association Consortium (BCAC). The SNPs were then tested for interaction with tumor TP53 status (n = 4,610) and anthracycline treatment (n = 17,828). For SNPs interacting with anthracycline treatment, siRNA knockdown experiments were carried out to validate candidate genes.In the test for interaction between SNP genotype and TP53 status, we identified one locus, represented by rs10916264 (p(interaction) = 3.44 × 10-5; FDR-adjusted p = 0.0011) in estrogen receptor (ER) positive cases. The rs10916264 AA genotype associated with worse survival among cases with ER-positive, TP53-positive tumors (hazard ratio [HR] 2.36, 95% confidence interval [C.I] 1.45 - 3.82). This is a cis-eQTL locus for FBXO28 and TP53BP2; expression levels of these genes were associated with patient survival specifically in ER-positive, TP53-mutated tumors. Additionally, the SNP rs798755 was associated with survival in interaction with anthracycline treatment (p(interaction) = 9.57 × 10-5, FDR-adjusted p = 0.0130). RNAi-based depletion of a predicted regulatory target gene, FAM53A, indicated that this gene can modulate doxorubicin sensitivity in breast cancer cell lines.If confirmed in independent data sets, these results may be of clinical relevance in the development of prognostic and predictive marker panels for breast cancer.