Fridley, Brooke L.

Abstract Variation in drug response results from a combination of factors that include differences in gender, ethnicity, and environment, as well as genetic variation that may result in differences in mRNA and protein expression. This article presents two integrative analytic approaches that make use of both genome-wide SNP and mRNA expression data available on the same set of subjects: a step-wise integrative approach and a comprehensive analysis using sparse canonical correlation analysis (SCCA). In addition to applying standard SCCA, we present a novel modification of SCCA which allows different weighting for the various pair-wise relationships in the SCCA. These integrative approaches are illustrated with both simulated data and data from a pharmacogenomic study of the drug gemcitabine. Results from these analyses found little overlap in terms of genes detected, possibly detecting different biological mechanisms. In addition, we found the proposed weighted SCCA to outperform its unweighted counterpart in detecting associations between the genomic features and phenotype. Further research is needed to develop and assess new integrative methods for pharmacogenomic studies, as these types of analyses may uncover novel insights into the relationship between genomic variation and drug response.

Abstract Introduction: Androstenedione is a precursor for estrogen biosynthesis catalyzed by aromatase. Plasma androstenedione concentrations in postmenopausal women with early stage breast cancer display large interindividual variation (Ingle et al., Cancer Res. 70:3278–86, 2010). We performed a GWAS to identify single nucleotide polymorphisms (SNPs) associated with plasma levels of androstenedione in 776 postmenopausal women with resected early stage breast cancer. Methods: Plasma androstenedione concentrations were measured by GC-MS/MS. GWAS was performed using Illumina Human610-Quad BeadChips. The GWAS results were controlled for population stratification (PS) by Eigenstrat analyses. Plasma androstenedione concentrations were associated with genotype by using a quasi-likelihood model with adjustment for PS, BMI, age and other clinical variables associated with androstenedione concentrations (p<0.01). Imputation using HapMap2 data was performed within 200 kb on either side of regions containing SNPs with p<E-05. Results: 563,945 SNPs were analyzed, followed by imputation. The genotyped SNP with the lowest p value, rs4736317 (p=4.61E-06), was on chromosome 8 near the CYP11B1 and CYP11B2 genes. This SNP had a multiplicative effect of 1.128 (95% Cl: 1.07, 1.19) and a MAF of 0.40. After imputation, 20 SNPs with apparent p<E-06 were observed in this same region, with the lowest p value for rs9773022 (p=8.1E-08). We then performed functional genomic studies, and observed that two of the SNPs with low p values (rs4736312, p=1.49E-05 and rs4736350, p=6.34E-06) created apparent Oct-1 transcription factor binding sites (TFBS). Differential nuclear protein binding of Oct-1 to probes for these SNPs was verified with EMS and supershift assays using nuclear extract from A549 and NCI-H295R cells. Treatment of these two cell lines with 10 nM androstenedione led to the disruption of Oct-1 binding with the A (wild type) allele on rs4736312. rs4736312 is located in the 3'-UTR of CYP11B1 15 bp distant from response elements for the progesterone receptor and the glucocorticoid receptor. We also demonstrated that CYP11B1 and CYP11B2 were induced by estradiol in U2OS cell lines that express ER receptor and by androstenedione in H295R cells that express the androgen receptor. Transfection of a reporter gene construct with the rs4736312 SNP sequences into COS-1 and H295R cells suppressed reporter gene activity for the variant SNP sequence as compared to wild type, an observation that may correlate with the higher plasma levels of androstenedione in women (549× 328 vs 410×208 mean×SD, M/ml) with the variant allele. Conclusion: We have identified a series of SNPs in a region of chromosome 8 in and around CYP11B1 and CYP11B2 that are associated with elevated plasma androstenedione concentrations in postmenopausal women with resected early stage breast cancer. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2011 Nov 12-16; San Francisco, CA. Philadelphia (PA): AACR; Mol Cancer Ther 2011;10(11 Suppl):Abstract nr C114.

Abstract It is believed that multiple genetic variants with small individual effects contribute to the risk of alcohol dependence. Such polygenic effects are difficult to detect in genome-wide association studies that test for association of the phenotype with each single nucleotide polymorphism (SNP) individually. To overcome this challenge, gene-set analysis (GSA) methods that jointly test for the effects of pre-defined groups of genes have been proposed. Rather than testing for association between the phenotype and individual SNPs, these analyses evaluate the global evidence of association with a set of related genes enabling the identification of cellular or molecular pathways or biological processes that play a role in development of the disease. It is hoped that by aggregating the evidence of association for all available SNPs in a group of related genes, these approaches will have enhanced power to detect genetic associations with complex traits. We performed GSA using data from a genome-wide study of 1165 alcohol-dependent cases and 1379 controls from the Study of Addiction: Genetics and Environment (SAGE), for all 200 pathways listed in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Results demonstrated a potential role of the ‘synthesis and degradation of ketone bodies’ pathway. Our results also support the potential involvement of the ‘neuroactive ligand–receptor interaction’ pathway, which has previously been implicated in addictive disorders. These findings demonstrate the utility of GSA in the study of complex disease, and suggest specific directions for further research into the genetic architecture of alcohol dependence.

Abstract Purpose: Glioblastoma is a devastating, incurable disease with few known prognostic factors. Here, we present the first genome-wide survival and validation study for glioblastoma. Experimental Design: Cox regressions for survival with 314,635 inherited autosomal single-nucleotide polymorphisms (SNP) among 315 San Francisco Adult Glioma Study patients for discovery and three independent validation data sets [87 Mayo Clinic, 232 glioma patients recruited from several medical centers in Southeastern United States (GliomaSE), and 115 The Cancer Genome Atlas patients] were used to identify SNPs associated with overall survival for Caucasian glioblastoma patients treated with the current standard of care, resection, radiation, and temozolomide (total n = 749). Tumor expression of the gene that contained the identified prognostic SNP was examined in three separate data sets (total n = 619). Genotype imputation was used to estimate hazard ratios (HR) for SNPs that had not been directly genotyped. Results: From the discovery and validation analyses, we identified a variant in single-stranded DNA-binding protein 2 (SSBP2) on 5q14.1 associated with overall survival in combined analyses (HR, 1.64; P = 1.3 × 10−6). Expression of SSBP2 in tumors from three independent data sets also was significantly related to patient survival (P = 5.3 × 10−4). Using genotype imputation, the SSBP2 SNP rs17296479 had the strongest statistically significant genome-wide association with poorer overall patient survival (HR, 1.79; 95% CI, 1.45-2.22; P = 1.0 × 10−7). Conclusion: The minor allele of SSBP2 SNP rs17296479 and the increased tumor expression of SSBP2 were statistically significantly associated with poorer overall survival among glioblastoma patients. With further confirmation, previously unrecognized inherited variations influencing survival may warrant inclusion in clinical trials to improve randomization. Unaccounted for genetic influence on survival could produce unwanted bias in such studies. Clin Cancer Res; 18(11); 3154–62. ©2012 AACR.

With new technologies, multiple types of genomic data are commonly collected on a single set of samples. However, standard analysis methods concentrate on a single data type at a time and ignore the relationships between genes, proteins, and biochemical reactions that give rise to complex phenotypes. In this paper, we propose a novel integrative model to incorporate multiple types of genomic data into an association analysis with a complex phenotype. The method combines path analysis and stochastic search variable selection into a Bayesian hierarchical model that simultaneously identifies both direct and indirect genomic effects on the phenotype. Results from a simulation study and application of the Bayesian model to a pharmacogenomic study of the drug gemcitabine demonstrate greater sensitivity to detect genomic effects in some simulation scenarios, when compared to the standard single data type analysis. Further research is required to extend and modify this integrative modeling framework to increase computational efficiency to best capitalize on the wealth of information available across multiple genomic data types. Genet. Epidemiol. 36:352–359, 2012. © 2012 Wiley Periodicals, Inc.

Abstract Background: Overexpression of mitotic kinases has been associated with prognosis, histologic grade, and clinical stage in ovarian cancer, but the relationship between inherited variation in these genes and ovarian cancer risk has not been well defined. Methods: We measured associations between 397 single nucleotide polymorphisms (SNPs) from 67 mitotic kinases and invasive epithelial ovarian cancer risk in two case–control studies (n = 671 cases; n = 939 controls). Thirty-six candidate SNPs (P < 0.05) were assessed in a replication analysis consisting of three additional studies (n = 1,094 cases; n = 829 controls). Results: In initial analysis, thirty-six SNPs were suggestive of association with risk of serous ovarian cancer, all subtypes of ovarian cancer, or both (P < 0.05). Replication analyses suggested an association between rs2125846 in the Nemo-like kinase (NLK) gene and ovarian cancer (serous OR = 1.36, 95% CI: 1.11–1.67, P = 1.77 × 10−3; all subtypes OR = 1.30, 95% CI: 1.08–1.56, P = 2.97 × 10−3). Furthermore, rs2125846 was associated with risk in the combined discovery and replication sets (serous OR = 1.33, 95% CI: 1.15–1.54; all subtypes OR = 1.27, 95% CI: 1.12–1.45). Conclusions: Variation in NLK may be associated with risk of invasive epithelial ovarian cancer. Further studies are needed to confirm and understand the biologic relationship between this mitotic kinase and ovarian cancer risk. Impact: An association between SNPs in NLK and ovarian cancer may provide biologic insight into the development of this disease. Cancer Epidemiol Biomarkers Prev; 21(3); 523–8. ©2012 AACR.

Abstract Background: Genome-wide association studies (GWAS) for epithelial ovarian cancer (EOC), the most lethal gynecologic malignancy, have identified novel susceptibility loci. GWAS for survival after EOC have had more limited success. The association of each single-nucleotide polymorphism (SNP) individually may not be well suited to detect small effects of multiple SNPs, such as those operating within the same biologic pathway. Gene set analysis (GSA) overcomes this limitation by assessing overall evidence for association of a phenotype with all measured variation in a set of genes. Methods: To determine gene sets associated with EOC overall survival, we conducted GSA using data from two large GWAS (N cases = 2,813, N deaths = 1,116), with a novel Principal Component-Gamma GSA method. Analysis was completed for all cases and then separately for high-grade serous histologic subtype. Results: Analysis of the high-grade serous subjects resulted in 43 gene sets with P < 0.005 (1.7%); of these, 21 gene sets had P < 0.10 in both GWAS, including intracellular signaling pathway (P = 7.3 × 10−5) and macrolide binding (P = 6.2 × 10−4) gene sets. The top gene sets in analysis of all cases were meiotic mismatch repair (P = 6.3 × 10−4) and macrolide binding (P = 1.0 × 10−3). Of 18 gene sets with P < 0.005 (0.7%), eight had P < 0.10 in both GWAS. Conclusion: This research detected novel gene sets associated with EOC survival. Impact: Novel gene sets associated with EOC survival might lead to new insights and avenues for development of novel therapies for EOC and pharmacogenomic studies. Cancer Epidemiol Biomarkers Prev; 21(3); 529–36. ©2012 AACR.

Abstract The importance of inflammation pathways to the development of many human cancers prompted us to examine the associations between single-nucleotide polymorphisms (SNP) in inflammation-related genes and risk of ovarian cancer. In a multisite case–control study, we genotyped SNPs in a large panel of inflammatory genes in 930 epithelial ovarian cancer cases and 1,037 controls using a custom array and analyzed by logistic regression. SNPs with P < 0.10 were evaluated among 3,143 cases and 2,102 controls from the Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI) post-GWAS collaboration. Combined analysis revealed association with SNPs rs17561 and rs4848300 in the interleukin gene IL1A which varied by histologic subtype (Pheterogeneity = 0.03). For example, IL1A rs17561, which correlates with numerous inflammatory phenotypes, was associated with decreased risk of clear cell, mucinous, and endometrioid subtype, but not with the most common serous subtype. Genotype at rs1864414 in the arachidonate 5-lipoxygenase ALOX5 was also associated with decreased risk. Thus, inherited variation in IL1A and ALOX5 seems to affect ovarian cancer risk which, for IL1A, is limited to rarer subtypes. Given the importance of inflammation in tumorigenesis and growing evidence of subtype-specific features in ovarian cancer, functional investigations will be important to help clarify the importance of inherited variation related to inflammation in ovarian carcinogenesis. Cancer Res; 72(5); 1064–9. ©2012 AACR.

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 Background: Invasive ovarian cancer is a significant cause of gynecologic cancer mortality. Methods: We examined whether this mortality was associated with inherited variation in approximately 170 candidate genes/regions [993 single-nucleotide polymorphisms (SNPs)] in a multistage analysis based initially on 312 Mayo Clinic cases (172 deaths). Additional analyses used The Cancer Genome Atlas (TCGA; 127 cases, 62 deaths). For the most compelling gene, we immunostained Mayo Clinic tissue microarrays (TMA, 326 cases) and conducted consortium-based SNP replication analysis (2,560 cases, 1,046 deaths). Results: The strongest initial mortality association was in HGF (hepatocyte growth factor) at rs1800793 (HR = 1.7, 95% CI = 1.3–2.2, P = 2.0 × 10−5) and with overall variation in HGF (gene-level test, P = 3.7 × 10−4). Analysis of TCGA data revealed consistent associations [e.g., rs5745709 (r2 = 0.96 with rs1800793): TCGA HR = 2.4, CI = 1.4–4.1, P = 2.2 × 10−3; Mayo Clinic + TCGA HR = 1.6, CI = 1.3–1.9, P = 7.0 × 10−5] and suggested genotype correlation with reduced HGF mRNA levels (P = 0.01). In Mayo Clinic TMAs, protein levels of HGF, its receptor MET (C-MET), and phospho-MET were not associated with genotype and did not serve as an intermediate phenotype; however, phospho-MET was associated with reduced mortality (P = 0.01) likely due to higher expression in early-stage disease. In eight additional ovarian cancer case series, HGF rs5745709 was not associated with mortality (HR = 1.0, CI = 0.9–1.1, P = 0.87). Conclusions: We conclude that although HGF signaling is critical to migration, invasion, and apoptosis, it is unlikely that HGF genetic variation plays a major role in ovarian cancer mortality. Furthermore, any minor role is not related to genetically-determined expression. Impact: Our study shows the utility of multiple data types and multiple data sets in observational studies. Cancer Epidemiol Biomarkers Prev; 20(8); 1638–48. ©2011 AACR.