Cheng, Jin Q.

Abstract Background: Genome-wide association studies (GWAS) have identified multiple loci associated with epithelial ovarian cancer (EOC) susceptibility, but further progress requires integration of epidemiology and biology to illuminate true risk loci below genome-wide significance levels (P < 5 × 10−8). Most risk SNPs lie within non–protein-encoding regions, and we hypothesize that long noncoding RNA (lncRNA) genes are enriched at EOC risk regions and represent biologically relevant functional targets. Methods: Using imputed GWAS data from about 18,000 invasive EOC cases and 34,000 controls of European ancestry, the GENCODE (v19) lncRNA database was used to annotate SNPs from 13,442 lncRNAs for permutation-based enrichment analysis. Tumor expression quantitative trait locus (eQTL) analysis was performed for sub-genome-wide regions (1 × 10−5 > P > 5 × 10−8) overlapping lncRNAs. Results: Of 5,294 EOC-associated SNPs (P < 1.0 × 10−5), 1,464 (28%) mapped within 53 unique lncRNAs and an additional 3,484 (66%) SNPs were correlated (r2 > 0.2) with SNPs within 115 lncRNAs. EOC-associated SNPs comprised 130 independent regions, of which 72 (55%) overlapped with lncRNAs, representing a significant enrichment (P = 5.0 × 10−4) that was more pronounced among a subset of 5,401 lncRNAs with active epigenetic regulation in normal ovarian tissue. EOC-associated lncRNAs and their putative promoters and transcription factors were enriched for biologically relevant pathways and eQTL analysis identified five novel putative risk regions with allele-specific effects on lncRNA gene expression. Conclusions: lncRNAs are significantly enriched at EOC risk regions, suggesting a mechanistic role for lncRNAs in driving predisposition to EOC. Impact: lncRNAs represent key candidates for integrative epidemiologic and functional studies. Further research on their biologic role in ovarian cancer is indicated. Cancer Epidemiol Biomarkers Prev; 26(1); 116–25. ©2016 AACR.

Abstract Background: Genome wide association studies (GWAS) have identified 20 loci associated with epithelial ovarian cancer (EOC) risk. Additional EOC risk loci await identification, and biological approaches may represent a useful strategy for overcoming sample size limitations. The ENCODE project recently concluded that only ∼1.2% of the genome encodes proteins, but at least 20% exhibits biological function and over 80% exhibits biochemical indices of function. Post-GWAS follow-up studies have firmly established that some associations are due to inter-individual differences in non-coding RNAs (ncRNAs). To estimate the magnitude of impact of variants in long non-coding RNAs (lncRNAs), we assessed enrichment of EOC-associated SNPs in lncRNA regions by comparing the density of EOC-associated loci between lncRNA regions, non-lncRNA regions, and the whole genome. Methods: The study population included ∼18,000 invasive EOC cases and 34,000 healthy controls of European ancestry from the Ovarian Cancer Association Consortium with GWAS data imputed to 1000 Genomes Project (1KGP) density. Density was compared with two measures: average chromosome length required for one EOC-associated SNP (kb/locus) and average number of tested SNPs containing one EOC-associated SNP (SNPs/locus). Coordinates and annotation for 13,870 lncRNA genes were downloaded from the publicly available GENCODE (v19) database and were used to annotate SNPs from the 1KGP imputed GWAS data, yielding 13,192 lncRNA genes with biallelic variants imputed at rsquare ≥0.25. A SNP was considered EOC-associated if it reached a significance level of P<10−5. Results: Genome-wide, 11.6% of the 15,123,646 tested SNPs fell within lncRNA regions. Of 5,294 EOC-associated SNPs that were identified, 27.7% were located within lncRNAs, a higher frequency than expected (<v:shape style = “WIDTH: 16.5pt; HEIGHT: 18pt” id = _x0000_i1026 equationxml = ‘14pχ2′ type = “#_x0000_t75”> = 9.35×10−22). Additionally, EOC-associated SNPs were located on average every 237.6 kb (or 1204 tested SNPs) in lncRNA regions versus every 657.7 kb (or 3489 SNPS) in non-lncRNA regions and every 541.5 kb (or 2857 SNPs) in the whole genome, clearly demonstrating an enrichment of EOC-associated loci in lncRNA regions. A total of 1,464 lncRNA SNPs from 53 unique lncRNAs associated with EOC risk (P<10−5); 70 of the identified lncRNA SNPs had P-values between 10−20 and 10−32, and 221 had P-values between 10−10 and 10−15, representing 21 unique regions > 500kb apart. Conclusions: These results justify efforts to integrate functional and population level genetic data on lncRNA regions to explore a new paradigm that addresses the problem of the missing heritability for cancer where restrictions in sample size and statistical power limit the ability to discover additional risk loci. Citation Format: Thomas A. Sellers, Brett M. Reid, Y. Ann Chen, Hui-Yi Lin, Edward Richards, Jamie Teer, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Ellen Goode, Edwin Iverson, Leigh Pearce, Paul Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Jin Cheng, Simon Gayther, Jennifer Permuth-Wey, on behalf of the Ovarian Cancer Association Consortium. Evidence that long non-coding RNA variants associate with epithelial ovarian cancer risk. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4633. doi:10.1158/1538-7445.AM2015-4633

Abstract Mi(cro)RNAs are short non-coding RNA molecules that play a key role in carcinogenesis by regulating tumor suppressors and oncogenes. Published data has implicated miRNAs in ovarian cancer (OC) development and progression, and we hypothesize that single nucleotide polymorphisms (SNPs) in sites of miRNA: messenger RNA (mRNA) binding may influence OC risk by altering expression levels of targeted mRNAs. To evaluate associations between SNPs in miRNA binding sites and OC risk, we used bioinformatics tools and public databases to identify miRNA binding SNPs relevant to ovarian cancer. We then evaluated the frequency of approximately 1,100 SNPs in ∼20,000 epithelial OC cases and 20,000 controls represented in the international Collaborative Oncological Gene-environment Study. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI) between genotypes and case status, with adjustment for the first five principal components representing European ancestry. Log-additive genetic models were applied to each SNP, modeling the number of copies of the minor allele. Subgroup analysis was conducted for serous adenocarcinomas, the most predominant histologic subtype of epithelial OC. Preliminary analysis of 11,215 invasive epithelial OCs and 15,910 controls included four SNPs with P<10−6 and 1 SNP with P<10−5. Analysis of 6,215 invasive serous adenocarcinomas revealed two SNPs with P<10−7 and three SNPs with P<10−6. Noteworthy is rs6104808, a C>T SNP predicted to reside in a highly conserved miRNA binding site of BTBD3 (BTB POZ domain containing 3), a transcription factor suggested to be involved in cell proliferation, cancer progression, and response to platinum-based chemotherapy. rs6104808 was associated with a decreased risk of OC overall (OR (95% CI): 0.42 (0.28-0.62), P=1.26 x 10−5)), and the association was stronger among cases with serous adenocarcinoma (OR (95% CI): 0.30 (0.17-0.50), P=9.89 x 10−6)). This represents the largest, most comprehensive epidemiologic study to date to evaluate associations between miRNA binding site SNPs and OC susceptibility. Although based on a preliminary analysis, these data suggest a significant role and putative biological mechanism for common germline variants in miRNA binding sites leading to cancer development. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2927. doi:1538-7445.AM2012-2927

Abstract Background: Single nucleotide polymorphisms (SNP) in microRNA-related genes have been associated with epithelial ovarian cancer (EOC) risk in two reports, yet associated alleles may be inconsistent across studies. Methods: We conducted a pooled analysis of previously identified SNPs by combining genotype data from 3,973 invasive EOC cases and 3,276 controls from the Ovarian Cancer Association Consortium. We also conducted imputation to obtain dense coverage of genes and comparable genotype data for all studies. In total, 226 SNPs within 15 kb of 4 miRNA biogenesis genes (DDX20, DROSHA, GEMIN4, and XPO5) and 23 SNPs located within putative miRNA binding sites of 6 genes (CAV1, COL18A1, E2F2, IL1R1, KRAS, and UGT2A3) were genotyped or imputed and analyzed in the entire dataset. Results: After adjustment for European ancestry, no overall association was observed between any of the analyzed SNPs and EOC risk. Conclusions: Common variants in these evaluated genes do not seem to be strongly associated with EOC risk. Impact: This analysis suggests earlier associations between EOC risk and SNPs in these genes may have been chance findings, possibly confounded by population admixture. To more adequately evaluate the relationship between genetic variants and cancer risk, large sample sizes are needed, adjustment for population stratification should be carried out, and use of imputed SNP data should be considered. Cancer Epidemiol Biomarkers Prev; 20(8); 1793–7. ©2011 AACR.

Abstract Defective microRNA (miRNA) biogenesis contributes to the development and progression of epithelial ovarian cancer (EOC). In this study, we examined the hypothesis that single nucleotide polymorphisms (SNP) in miRNA biogenesis genes may influence EOC risk. In an initial investigation, 318 SNPs in 18 genes were evaluated among 1,815 EOC cases and 1,900 controls, followed up by a replicative joint meta-analysis of data from an additional 2,172 cases and 3,052 controls. Of 23 SNPs from 9 genes associated with risk (empirical P < 0.05) in the initial investigation, the meta-analysis replicated 6 SNPs from the DROSHA, FMR1, LIN28, and LIN28B genes, including rs12194974 (G>A), an SNP in a putative transcription factor binding site in the LIN28B promoter region (summary OR = 0.90, 95% CI: 0.82–0.98; P = 0.015) which has been recently implicated in age of menarche and other phenotypes. Consistent with reports that LIN28B overexpression in EOC contributes to tumorigenesis by repressing tumor suppressor let-7 expression, we provide data from luciferase reporter assays and quantitative RT-PCR to suggest that the inverse association among rs12194974 A allele carriers may be because of reduced LIN28B expression. Our findings suggest that variants in LIN28B and possibly other miRNA biogenesis genes may influence EOC susceptibility. Cancer Res; 71(11); 3896–903. ©2011 AACR.