A transcriptome-wide association study among 97,898 women to identify candidate susceptibility genes for epithelial ovarian cancer risk

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in non-coding regions, and causal genes underlying these associations remain largely unknown. Here we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian-tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P<2.2×10-6, we identified 35 genes including FZD4 at 11q14.2 (Z=5.08, P=3.83×10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly-associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and 3 genes remained (P<1.47 x 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis

Evaluation of copy-number variants as modifiers of breast and ovarian cancer risk for BRCA1 pathogenic variant carriers.

Genome-wide studies of patients carrying pathogenic variants (mutations) in BRCA1 or BRCA2 have reported strong associations between single-nucleotide polymorphisms (SNPs) and cancer risk. To conduct the first genome-wide association analysis of copy-number variants (CNVs) with breast or ovarian cancer risk in a cohort of 2500 BRCA1 pathogenic variant carriers, CNV discovery was performed using multiple calling algorithms and Illumina 610k SNP array data from a previously published genome-wide association study. Our analysis, which focused on functionally disruptive genomic deletions overlapping gene regions, identified a number of loci associated with risk of breast or ovarian cancer for BRCA1 pathogenic variant carriers. Despite only including putative deletions called by at least two or more algorithms, detection of selected CNVs by ancillary molecular technologies only confirmed 40% of predicted common (>1% allele frequency) variants. These include four loci that were associated (unadjusted P<0.05) with breast cancer (GTF2H2, ZNF385B, NAALADL2 and PSG5), and two loci associated with ovarian cancer (CYP2A7 and OR2A1). An interesting finding from this study was an association of a validated CNV deletion at the CYP2A7 locus (19q13.2) with decreased ovarian cancer risk (relative risk=0.50, P=0.007). Genomic analysis found this deletion coincides with a region displaying strong regulatory potential in ovarian tissue, but not in breast epithelial cells. This study highlighted the need to verify CNVs in vitro, but also provides evidence that experimentally validated CNVs (with plausible biological consequences) can modify risk of breast or ovarian cancer in BRCA1 pathogenic variant carriers.European Journal of Human Genetics advance online publication, 1 February 2017; doi:10.1038/ejhg.2016.203.Multiple funders listed in the article

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

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.Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of similar to 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.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 x 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.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.</p

Polygenic risk modeling for prediction of epithelial ovarian cancer risk

Polygenic risk scores (PRS) for epithelial ovarian cancer (EOC) have the potential to improve risk stratification. Joint estimation of Single Nucleotide Polymorphism (SNP) effects in models could improve predictive performance over standard approaches of PRS construction. Here, we implemented computationally efficient, penalized, logistic regression models (lasso, elastic net, stepwise) to individual level genotype data and a Bayesian framework with continuous shrinkage, "select and shrink for summary statistics" (S4), to summary level data for epithelial non-mucinous ovarian cancer risk prediction. We developed the models in a dataset consisting of 23,564 non-mucinous EOC cases and 40,138 controls participating in the Ovarian Cancer Association Consortium (OCAC) and validated the best models in three populations of different ancestries: prospective data from 198,101 women of European ancestries; 7,669 women of East Asian ancestries; 1,072 women of African ancestries, and in 18,915 BRCA1 and 12,337 BRCA2 pathogenic variant carriers of European ancestries. In the external validation data, the model with the strongest association for non-mucinous EOC risk derived from the OCAC model development data was the S4 model (27,240 SNPs) with odds ratios (OR) of 1.38 (95% CI: 1.28-1.48, AUC: 0.588) per unit standard deviation, in women of European ancestries; 1.14 (95% CI: 1.08-1.19, AUC: 0.538) in women of East Asian ancestries; 1.38 (95% CI: 1.21-1.58, AUC: 0.593) in women of African ancestries; hazard ratios of 1.36 (95% CI: 1.29-1.43, AUC: 0.592) in BRCA1 pathogenic variant carriers and 1.49 (95% CI: 1.35-1.64, AUC: 0.624) in BRCA2 pathogenic variant carriers. Incorporation of the S4 PRS in risk prediction models for ovarian cancer may have clinical utility in ovarian cancer prevention programs

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.Peer reviewe

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. Weevaluated whether PSVs inBRCA1/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 30 region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001c.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.Peer reviewe

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer.

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC

Male breast cancer in BRCA1 and BRCA2 mutation carriers: Pathology data from the Consortium of Investigators of Modifiers of BRCA1/2

Background: BRCA1 and, more commonly, BRCA2 mutations are associated with increased risk of male breast cancer (MBC). However, only a paucity of data exists on the pathology of breast cancers (BCs) in men with BRCA1/2 mutations. Using the largest available dataset, we determined whether MBCs arising in BRCA1/2 mutation carriers display specific pathologic features and whether these features differ from those of BRCA1/2 female BCs (FBCs). Methods: We characterised the pathologic features of 419 BRCA1/2 MBCs and, using logistic regression analysis, contrasted those with data from 9675 BRCA1/2 FBCs and with population-based data from 6351 MBCs in the Surveillance, Epidemiology, and End Results (SEER) database. Results: Among BRCA2 MBCs, grade significantly decreased with increasing age at diagnosis (P = 0.005). Compared with BRCA2 FBCs, BRCA2 MBCs were of significantly higher stage (P for trend = 2 × 10-5) and higher grade (P for trend = 0.005) and were more likely to be oestrogen receptor-positive [odds ratio (OR) 10.59; 95 % confidence interval (CI) 5.15-21.80] and progesterone receptor-positive (OR 5.04; 95 % CI 3.17-8.04). With the exception of grade, similar patterns of associations emerged when we compared BRCA1 MBCs and FBCs. BRCA2 MBCs also presented with higher grade than MBCs from the SEER database (P for trend = 4 × 10-12). Conclusions: On the basis of the largest series analysed to date, our results show that BRCA1/2 MBCs display distinct pathologic characteristics compared with BRCA1/2 FBCs, and we identified a specific BRCA2-associated MBC phenotype characterised by a variable suggesting greater biological aggressiveness (i.e., high histologic grade). These findings could lead to the development of gender-specific risk prediction models and guide clinical strategies appropriate for MBC management

Mendelian randomisation study of height and body mass index as modifiers of ovarian cancer risk in 22,588 BRCA1 and BRCA2 mutation carriers.

BackgroundHeight and body mass index (BMI) are associated with higher ovarian cancer risk in the general population, but whether such associations exist among BRCA1/2 mutation carriers is unknown.MethodsWe applied a Mendelian randomisation approach to examine height/BMI with ovarian cancer risk using the Consortium of Investigators for the Modifiers of BRCA1/2 (CIMBA) data set, comprising 14,676 BRCA1 and 7912 BRCA2 mutation carriers, with 2923 ovarian cancer cases. We created a height genetic score (height-GS) using 586 height-associated variants and a BMI genetic score (BMI-GS) using 93 BMI-associated variants. Associations were assessed using weighted Cox models.ResultsObserved height was not associated with ovarian cancer risk (hazard ratio [HR]: 1.07 per 10-cm increase in height, 95% confidence interval [CI]: 0.94-1.23). Height-GS showed similar results (HR = 1.02, 95% CI: 0.85-1.23). Higher BMI was significantly associated with increased risk in premenopausal women with HR = 1.25 (95% CI: 1.06-1.48) and HR = 1.59 (95% CI: 1.08-2.33) per 5-kg/m2 increase in observed and genetically determined BMI, respectively. No association was found for postmenopausal women. Interaction between menopausal status and BMI was significant (Pinteraction ConclusionOur observation of a positive association between BMI and ovarian cancer risk in premenopausal BRCA1/2 mutation carriers is consistent with findings in the general population

Male breast cancer in BRCA1 and BRCA2 mutation carriers: pathology data from the Consortium of Investigators of Modifiers of BRCA1/2.

Background BRCA1 and, more commonly, BRCA2 mutations are associated with increased risk of male breast cancer (MBC). However, only a paucity of data exists on the pathology of breast cancers (BCs) in men with BRCA1/2 mutations. Using the largest available dataset, we determined whether MBCs arising in BRCA1/2 mutation carriers display specific pathologic features and whether these features differ from those of BRCA1/2 female BCs (FBCs).Methods We characterised the pathologic features of 419 BRCA1/2 MBCs and, using logistic regression analysis, contrasted those with data from 9675 BRCA1/2 FBCs and with population-based data from 6351 MBCs in the Surveillance, Epidemiology, and End Results (SEER) database.Results Among BRCA2 MBCs, grade significantly decreased with increasing age at diagnosis (P = 0.005). Compared with BRCA2 FBCs, BRCA2 MBCs were of significantly higher stage (P for trend = 2 × 10(-5)) and higher grade (P for trend = 0.005) and were more likely to be oestrogen receptor-positive [odds ratio (OR) 10.59; 95 % confidence interval (CI) 5.15-21.80] and progesterone receptor-positive (OR 5.04; 95 % CI 3.17-8.04). With the exception of grade, similar patterns of associations emerged when we compared BRCA1 MBCs and FBCs. BRCA2 MBCs also presented with higher grade than MBCs from the SEER database (P for trend = 4 × 10(-12)).Conclusions On the basis of the largest series analysed to date, our results show that BRCA1/2 MBCs display distinct pathologic characteristics compared with BRCA1/2 FBCs, and we identified a specific BRCA2-associated MBC phenotype characterised by a variable suggesting greater biological aggressiveness (i.e., high histologic grade). These findings could lead to the development of gender-specific risk prediction models and guide clinical strategies appropriate for MBC management