Rare and low frequency variants are not well covered in most germline genotyping arrays and are understudied in relation to epithelial ovarian cancer (EOC) risk. To address this gap, we used genotyping arrays targeting rarer protein-coding variation in 8,165 EOC cases and 11,619 controls from the international Ovarian Cancer Association Consortium (OCAC). Pooled association analyses were conducted at the variant and gene level for 98,543 variants directly genotyped through two exome genotyping projects. Only common variants that represent or are in strong linkage disequilibrium (LD) with previously-identified signals at established loci reached traditional thresholds for exome-wide significance (P < 5.0 × 10 ^{-  7}). One of the most significant signals (P_{all histologies}=1.01 × 10 ^{-  13};P_serous=3.54 × 10 ^{-  14}) occurred at 3q25.31 for rs62273959, a missense variant mapping to the LEKR1 gene that is in LD (r²=0.90) with a previously identified 'best hit' (rs7651446) mapping to an intron of TIPARP. Suggestive associations (5.0 × 10 ^{-  5}>P≥5.0 ×10 ^{-  7}) were detected for rare and low-frequency variants at 16 novel loci. Four rare missense variants were identified (ACTBL2 rs73757391 (5q11.2), BTD rs200337373 (3p25.1), KRT13 rs150321809 (17q21.2) and MC2R rs104894658 (18p11.21)), but only MC2R rs104894668 had a large effect size (OR = 9.66). Genes most strongly associated with EOC risk included ACTBL2 (P_AML=3.23 × 10 ^{-  5}; P_SKAT-o=9.23 × 10 ^{-  4}) and KRT13 (P_AML=1.67 × 10 ^{-  4}; P_SKAT-o=1.07 × 10 ^{-  5}), reaffirming variant-level analysis. In summary, this large study identified several rare and low-frequency variants and genes that may contribute to EOC susceptibility, albeit with possible small effects. Future studies that integrate epidemiology, sequencing, and functional assays are needed to further unravel the unexplained heritability and biology of this disease.

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.

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 noncoding 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 three genes remained (P < 1.47 × 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.Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419-30. ©2018 AACR.