Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions.

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

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.

Genome-wide association analyses of ovarian cancer patients undergoing primary debulking surgery identify candidate genes for residual disease.

Survival from ovarian cancer depends on the resection status after primary surgery. We performed genome-wide association analyses for resection status of 7705 ovarian cancer patients, including 4954 with high-grade serous carcinoma (HGSOC), to identify variants associated with residual disease. The most significant association with resection status was observed for rs72845444, upstream of MGMT, in HGSOC (p = 3.9 × 10-8). In gene-based analyses, PPP2R5C was the most strongly associated gene in HGSOC after stage adjustment. In an independent set of 378 ovarian tumours from the AGO-OVAR 11 study, variants near MGMT and PPP2R5C correlated with methylation and transcript levels, and PPP2R5C mRNA levels predicted progression-free survival in patients with residual disease. MGMT encodes a DNA repair enzyme, and PPP2R5C encodes the B56γ subunit of the PP2A tumour suppressor. Our results link heritable variation at these two loci with resection status in HGSOC.

Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

Folate Intake and Ovarian Cancer Risk among Women with Endometriosis: A Case-Control Study from the Ovarian Cancer Association Consortium.

BACKGROUND: Although folate intake has not been associated with an increased risk of ovarian cancer overall, studies of other cancer types have suggested that high folate intake may promote carcinogenesis in precancerous lesions. Women with endometriosis (a potential precancerous lesion) have an increased risk of developing ovarian cancer; however, whether high folate intake increases risk in this group is unknown. METHODS: We conducted a pooled analysis of six case-control studies from the Ovarian Cancer Association Consortium to investigate the association between folate intake and risk of ovarian cancer among women with and without self-reported endometriosis. We included 570 cases/558 controls with and 5,171/7,559 without endometriosis. We used logistic regression to estimate odds ratios (OR) and 95% confidence intervals for the association between folate intake (dietary, supplemental, and total) and ovarian cancer risk. Finally, we used Mendelian randomization (MR) to evaluate our results using genetic markers as a proxy for folate status. RESULTS: Higher dietary folate intake was associated with an increased risk of ovarian cancer for women with endometriosis [OR, 1.37 (1.01-1.86)] but not for women without endometriosis. There was no association between supplemental folate intake and ovarian cancer risk for women with or without endometriosis. A similar pattern was seen using MR. CONCLUSIONS: High dietary folate intake may be associated with an increased risk of ovarian cancer among women with endometriosis. IMPACT: Women with endometriosis with high folate diets may be at increased risk of ovarian cancer. Further research is needed on the potential cancer-promoting effects of folate in this group.

Copy Number Variants Are Ovarian Cancer Risk Alleles at Known and Novel Risk Loci.

BACKGROUND: Known risk alleles for epithelial ovarian cancer (EOC) account for approximately 40% of the heritability for EOC. Copy number variants (CNVs) have not been investigated as EOC risk alleles in a large population cohort. METHODS: Single nucleotide polymorphism array data from 13 071 EOC cases and 17 306 controls of White European ancestry were used to identify CNVs associated with EOC risk using a rare admixture maximum likelihood test for gene burden and a by-probe ratio test. We performed enrichment analysis of CNVs at known EOC risk loci and functional biofeatures in ovarian cancer-related cell types. RESULTS: We identified statistically significant risk associations with CNVs at known EOC risk genes; BRCA1 (PEOC = 1.60E-21; OREOC = 8.24), RAD51C (Phigh-grade serous ovarian cancer [HGSOC] = 5.5E-4; odds ratio [OR]HGSOC = 5.74 del), and BRCA2 (PHGSOC = 7.0E-4; ORHGSOC = 3.31 deletion). Four suggestive associations (P < .001) were identified for rare CNVs. Risk-associated CNVs were enriched (P < .05) at known EOC risk loci identified by genome-wide association study. Noncoding CNVs were enriched in active promoters and insulators in EOC-related cell types. CONCLUSIONS: CNVs in BRCA1 have been previously reported in smaller studies, but their observed frequency in this large population-based cohort, along with the CNVs observed at BRCA2 and RAD51C gene loci in EOC cases, suggests that these CNVs are potentially pathogenic and may contribute to the spectrum of disease-causing mutations in these genes. CNVs are likely to occur in a wider set of susceptibility regions, with potential implications for clinical genetic testing and disease prevention.

No Association Between Polygenic Risk Scores for Cancer and Development of Radiation Therapy Toxicity.

PURPOSE: Our aim was to test whether updated polygenic risk scores (PRS) for susceptibility to cancer affect risk of radiation therapy toxicity. METHODS AND MATERIALS: Analyses included 9,717 patients with breast (n=3,078), prostate (n=5,748) or lung (n=891) cancer from Radiogenomics and REQUITE Consortia cohorts. Patients underwent potentially curative radiation therapy and were assessed prospectively for toxicity. Germline genotyping involved genome-wide single nucleotide polymorphism (SNP) arrays with nontyped SNPs imputed. PRS for each cancer were generated by summing literature-identified cancer susceptibility risk alleles: 352 breast, 136 prostate, and 24 lung. Weighted PRS were generated using log odds ratio (ORs) for cancer susceptibility. Standardized total average toxicity (STAT) scores at 2 and 5 years (breast, prostate) or 6 to 12 months (lung) quantified toxicity. Primary analysis tested late STAT, secondary analyses investigated acute STAT, and individual endpoints and SNPs using multivariable regression. RESULTS: Increasing PRS did not increase risk of late toxicity in patients with breast (OR, 1.000; 95% confidence interval [CI], 0.997-1.002), prostate (OR, 0.99; 95% CI, 0.98-1.00; weighted PRS OR, 0.93; 95% CI, 0.83-1.03), or lung (OR, 0.93; 95% CI, 0.87-1.00; weighted PRS OR, 0.68; 95% CI, 0.45-1.03) cancer. Similar results were seen for acute toxicity. Secondary analyses identified rs138944387 associated with breast pain (OR, 3.05; 95% CI, 1.86-5.01; P = 1.09 × 10-5) and rs17513613 with breast edema (OR, 0.94; 95% CI, 0.92-0.97; P = 1.08 × 10-5). CONCLUSIONS: Patients with increased polygenic predisposition to breast, prostate, or lung cancer can safely undergo radiation therapy with no anticipated excess toxicity risk. Some individual SNPs increase the likelihood of a specific toxicity endpoint, warranting validation in independent cohorts and functional studies to elucidate biologic mechanisms.

chromMAGMA: regulatory element-centric interrogation of risk variants.

Candidate causal risk variants from genome-wide association studies reside almost exclusively in noncoding regions of the genome and innovative approaches are necessary to understand their biological function. Multi-marker analysis of genomic annotation (MAGMA) is a widely used program that nominates candidate risk genes by mapping single-nucleotide polymorphism summary statistics from genome-wide association studies to gene bodies. We augmented MAGMA to create chromatin-MAGMA (chromMAGMA), a method to nominate candidate risk genes based on the presence of risk variants within noncoding regulatory elements (REs). We applied chromMAGMA to a genetic susceptibility dataset for epithelial ovarian cancer (EOC), a rare gynecologic malignancy characterized by high mortality. This identified 155 unique candidate EOC risk genes across five EOC histotypes; 83% (105/127) of high-grade serous ovarian cancer risk genes had not previously been implicated in this EOC histotype. Risk genes nominated by chromMAGMA converged on mRNA splicing and transcriptional dysregulation pathways. chromMAGMA is a pipeline that nominates candidate risk genes through a gene regulation-focused approach and helps interpret the biological mechanism of noncoding risk variants for complex diseases.

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.

Rare germline copy number variants (CNVs) and breast cancer risk.

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Comprehensive epithelial tubo-ovarian cancer risk prediction model incorporating genetic and epidemiological risk factors.

BACKGROUND: Epithelial tubo-ovarian cancer (EOC) has high mortality partly due to late diagnosis. Prevention is available but may be associated with adverse effects. A multifactorial risk model based on known genetic and epidemiological risk factors (RFs) for EOC can help identify women at higher risk who could benefit from targeted screening and prevention. METHODS: We developed a multifactorial EOC risk model for women of European ancestry incorporating the effects of pathogenic variants (PVs) in BRCA1, BRCA2, RAD51C, RAD51D and BRIP1, a Polygenic Risk Score (PRS) of arbitrary size, the effects of RFs and explicit family history (FH) using a synthetic model approach. The PRS, PV and RFs were assumed to act multiplicatively. RESULTS: Based on a currently available PRS for EOC that explains 5% of the EOC polygenic variance, the estimated lifetime risks under the multifactorial model in the general population vary from 0.5% to 4.6% for the first to 99th percentiles of the EOC risk distribution. The corresponding range for women with an affected first-degree relative is 1.9%-10.3%. Based on the combined risk distribution, 33% of RAD51D PV carriers are expected to have a lifetime EOC risk of less than 10%. RFs provided the widest distribution, followed by the PRS. In an independent partial model validation, absolute and relative 5-year risks were well calibrated in quintiles of predicted risk. CONCLUSION: This multifactorial risk model can facilitate stratification, in particular among women with FH of cancer and/or moderate-risk and high-risk PVs. The model is available via the CanRisk Tool (www.canrisk.org).

Functional analysis of the 1p34.3 risk locus implicates GNL2 in high-grade serous ovarian cancer.

The high-grade serous ovarian cancer (HGSOC) risk locus at chromosome 1p34.3 resides within a frequently amplified genomic region signifying the presence of an oncogene. Here, we integrate in silico variant-to-function analysis with functional studies to characterize the oncogenic potential of candidate genes in the 1p34.3 locus. Fine mapping of genome-wide association statistics identified candidate causal SNPs local to H3K27ac-demarcated enhancer regions that exhibit allele-specific binding for CTCF in HGSOC and normal fallopian tube secretory epithelium cells (FTSECs). SNP risk associations colocalized with eQTL for six genes (DNALI1, GNL2, RSPO1, SNIP1, MEAF6, and LINC01137) that are more highly expressed in carriers of the risk allele, and three (DNALI1, GNL2, and RSPO1) were upregulated in HGSOC compared to normal ovarian surface epithelium cells and/or FTSECs. Increased expression of GNL2 and MEAF6 was associated with shorter survival in HGSOC with 1p34.3 amplifications. Despite its activation of ß-catenin signaling, RSPO1 overexpression exerted no effects on proliferation or colony formation in our study of ovarian cancer and FTSECs. Instead, GNL2, MEAF6, and SNIP1 silencing impaired in vitro ovarian cancer cell growth. Additionally, GNL2 silencing diminished xenograft tumor formation, whereas overexpression stimulated proliferation and colony formation in FTSECs. GNL2 influences 60S ribosomal subunit maturation and global protein synthesis in ovarian cancer and FTSECs, providing a potential mechanism of how GNL2 upregulation might promote ovarian cancer development and mediate genetic susceptibility of HGSOC.

Genomic risk prediction of coronary artery disease in women with breast cancer: a prospective cohort study.

BACKGROUND: Advancements in cancer therapeutics have resulted in increases in cancer-related survival; however, there is a growing clinical dilemma. The current balancing of survival benefits and future cardiotoxic harms of oncotherapies has resulted in an increased burden of cardiovascular disease in breast cancer survivors. Risk stratification may help address this clinical dilemma. This study is the first to assess the association between a coronary artery disease-specific polygenic risk score and incident coronary artery events in female breast cancer survivors. METHODS: We utilized the Studies in Epidemiology and Research in Cancer Heredity prospective cohort involving 12,413 women with breast cancer with genotype information and without a baseline history of cardiovascular disease. Cause-specific hazard ratios for association of the polygenic risk score and incident coronary artery disease (CAD) were obtained using left-truncated Cox regression adjusting for age, genotype array, conventional risk factors such as smoking and body mass index, as well as other sociodemographic, lifestyle, and medical variables. RESULTS: Over a median follow-up of 10.3 years (IQR: 16.8) years, 750 incident fatal or non-fatal coronary artery events were recorded. A 1 standard deviation higher polygenic risk score was associated with an adjusted hazard ratio of 1.33 (95% CI 1.20, 1.47) for incident CAD. CONCLUSIONS: This study provides evidence that a coronary artery disease-specific polygenic risk score can risk-stratify breast cancer survivors independently of other established cardiovascular risk factors.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals.

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Pleiotropy-guided transcriptome imputation from normal and tumor tissues identifies candidate susceptibility genes for breast and ovarian cancer.

Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Identification of a Locus Near ULK1 Associated With Progression-Free Survival in Ovarian Cancer.

BACKGROUND: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. METHODS: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. RESULTS: We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10-8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10-8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. CONCLUSIONS: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. IMPACT: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options.See related commentary by Peres and Monteiro, p. 1604.

Author Correction: A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23162-4.

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers.

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.