Li, Qiyuan

Abstract Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10−8 seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type–specific expression quantitative trait locus and enhancer–gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10−5 in the three-cancer meta-analysis. Significance: We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052–67. ©2016 AACR. This article is highlighted in the In This Issue feature, p. 932

Abstract Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103–14. ©2016 AACR.

Abstract Genome-wide association studies (GWAS) have identified several susceptibility regions that confer risk of two or more cancers. A region at 19p13 is associated with both breast and ovarian cancer risk in the general population as well as in BRCA1/BRCA2 mutation carriers. Using the iCOGs custom genotyping array, we analyzed 438 SNPs in this region in 7,435 ER-negative breast cancer cases and 42,599 controls from BCAC, and 15,347 ovarian cancer cases and 30,845 controls from OCAC and 15,252 BRCA1 mutation carriers. We identified two significant peaks of association. Peak 1, consisting of 13 highly correlated SNPs (r2>0.95), spans a 19.4kb region harboring the BABAM1, ANKLE1 and ABHD8 genes. Peak 1 SNPs are associated with serous ovarian cancer (top hit rs4808075, P = 9.2×10-20), ER-negative breast cancer (rs61494113, P = 1.1×10-13) and breast cancer in BRCA1 mutation carriers (rs61494113, P = 7.7×10-16). Peak 2 spans over 100kb and contains 7 SNPs associated with breast cancer in BRCA1 mutation carriers only (rs4464206, P = 8.9×10-8). Multiple assays were used to establish the putative functional roles of risk SNPs and candidate susceptibility genes in this region. Expression quantitative trait locus (eQTL) analyses identified significant genotype-gene expression correlations between genotype at the Peak 1 SNP rs480816 and ABHD8 expression in ovarian tumors (P = 3.0×10-5) and normal breast tissues (P = 2.8×10-3). The Peak 2 SNP, rs7246243, was significantly associated with BABAM1 expression (P = 5×10-3) in normal breast tissue only. Overexpression of BABAM1 in normal ovarian and breast cells did not modulate neoplastic phenotypes, whereas overexpression of ABHD8 induced a significant reduction in migration in both cell types in vitro. Differentially expressed genes in ABHD8 models, identified by RNAseq, were associated with cellular movement (P<1.82×10-44). Chromosome conformation capture (3C) identified an interaction between the ABHD8 promoter and a region containing 4 risk SNPs (rs4808075, rs10419397, rs56069439 and rs4808076) in breast and ovarian normal and cancer cells but only rs56069439 coincides with epigenetic enhancer marks in breast and ovarian cells. CRISPR/Cas9 genome editing was used to delete a 57bp region containing rs56069439 and expression of BABAM1, ANKLE1 and ABHD8 were measured. Taken together, these data suggest a common functional mechanism underlying risk of breast and ovarian cancer at this locus, mediated by rs56069439 variants differentially regulating expression of ABHD8, a gene involved in cell motility. Citation Format: Kate Lawrenson, Siddhartha Kar, Karoline Kuckenbaeker, Stacey Edwards, Qiyuan Li, Jonathan Tyrer, Jonathan Beesley, Kyriaki Michailidou, Susan Ramus, Alison Dunning, Fergus Couch, Alvaro Monteiro, Jacques Simard, Matthew Freedman, Douglas Easton, Georgia Chenevix-Trench, Paul Pharoah, Antonis Antoniou, Simon Gayther, OCAC, BCAC and CIMBA. Common functional mechanisms underlying pleiotropy at the 19p13.1 breast and ovarian cancer cusceptibility locus. [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 2783. doi:10.1158/1538-7445.AM2015-2783

Abstract Background: Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations. Methods: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls). Results: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network. Conclusion: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development. Impact: Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization. Cancer Epidemiol Biomarkers Prev; 24(10); 1574–84. ©2015 AACR.

Abstract A genome-wide association study (GWAS) of prostate cancer in Kaiser Permanente health plan members (7,783 cases, 38,595 controls; 80.3% non-Hispanic white, 4.9% African-American, 7.0% East Asian, and 7.8% Latino) revealed a new independent risk indel rs4646284 at the previously identified locus 6q25.3 that replicated in PEGASUS (N = 7,539) and the Multiethnic Cohort (N = 4,679) with an overall P = 1.0 × 10−19 (OR, 1.18). Across the 6q25.3 locus, rs4646284 exhibited the strongest association with expression of SLC22A1 (P = 1.3 × 10−23) and SLC22A3 (P = 3.2 × 10−52). At the known 19q13.33 locus, rs2659124 (P = 1.3 × 10−13; OR, 1.18) nominally replicated in PEGASUS. A risk score of 105 known risk SNPs was strongly associated with prostate cancer (P < 1.0 × 10−8). Comparing the highest to lowest risk score deciles, the OR was 6.22 for non-Hispanic whites, 5.82 for Latinos, 3.77 for African-Americans, and 3.38 for East Asians. In non-Hispanic whites, the 105 risk SNPs explained approximately 7.6% of disease heritability. The entire GWAS array explained approximately 33.4% of heritability, with a 4.3-fold enrichment within DNaseI hypersensitivity sites (P = 0.004). Significance: Taken together, our findings of independent risk variants, ethnic variation in existing SNP replication, and remaining unexplained heritability have important implications for further clarifying the genetic risk of prostate cancer. Our findings also suggest that there may be much promise in evaluating understudied variation, such as indels and ethnically diverse populations. Cancer Discov; 5(8); 878–91. ©2015 AACR. This article is highlighted in the In This Issue feature, p. 783