Gene set methods aim to assess the overall evidence of association of a set of genes with a phenotype, such as disease or a quantitative trait. Multiple approaches for gene set analysis of expression data have been proposed. They can be divided into two types: competitive and self-contained. Benefits of self-contained methods include that they can be used for genome-wide, candidate gene, or pathway studies, and have been reported to be more powerful than competitive methods. We therefore investigated ten self-contained methods that can be used for continuous, discrete and time-to-event phenotypes. To assess the power and type I error rate for the various previously proposed and novel approaches, an extensive simulation study was completed in which the scenarios varied according to: number of genes in a gene set, number of genes associated with the phenotype, effect sizes, correlation between expression of genes within a gene set, and the sample size. In addition to the simulated data, the various methods were applied to a pharmacogenomic study of the drug gemcitabine. Simulation results demonstrated that overall Fisher's method and the global model with random effects have the highest power for a wide range of scenarios, while the analysis based on the first principal component and Kolmogorov-Smirnov test tended to have lowest power. The methods investigated here are likely to play an important role in identifying pathways that contribute to complex traits.

Genome-wide association studies (GWASs) assess correlation between traits and DNA sequence variation using large numbers of genetic variants such as single nucleotide polymorphisms (SNPs) distributed across the genome. A GWAS produces many trait-SNP associations with low p-values, but few are replicated in subsequent studies. We sought to determine if characteristics of the genomic loci associated with a trait could be used to identify initial associations with a higher chance of replication in a second cohort. Data from the age-related eye disease study (AREDS) of 100,000 SNPs on 395 subjects with and 198 without age-related macular degeneration (AMD) were employed. Loci highly associated with AMD were characterized based on the distribution of genotypes, level of significance, and clustering of adjacent SNPs also associated with AMD suggesting linkage disequilibrium or multiple effects. Forty nine loci were highly associated with AMD, including 3 loci (CFH, C2/BF, LOC387715/HTRA1) already known to contain important genetic risks for AMD. One additional locus (C3) reported during the course of this study was identified and replicated in an additional study group. Tag-SNPs and haplotypes for each locus were evaluated for association with AMD in additional cohorts to account for population differences between discovery and replication subjects, but no additional clearly significant associations were identified. Relying on a significant genotype tests using a log-additive model would have excluded 57% of the non-replicated and none of the replicated loci, while use of other SNP features and clustering might have missed true associations.

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

One difficult question facing researchers is how to prioritize SNPs detected from genetic association studies for functional studies. Often a list of the top M SNPs is determined based on solely the p-value from an association analysis, where M is determined by financial/time constraints. For many studies of complex diseases, multiple analyses have been completed and integrating these multiple sets of results may be difficult. One may also wish to incorporate biological knowledge, such as whether the SNP is in the exon of a gene or a regulatory region, into the selection of markers to follow-up. In this manuscript, we propose a Bayesian latent variable model (BLVM) for incorporating "features" about a SNP to estimate a latent "quality score", with SNPs prioritized based on the posterior probability distribution of the rankings of these quality scores. We illustrate the method using data from an ovarian cancer genome-wide association study (GWAS). In addition to the application of the BLVM to the ovarian GWAS, we applied the BLVM to simulated data which mimics the setting involving the prioritization of markers across multiple GWAS for related diseases/traits. The top ranked SNP by BLVM for the ovarian GWAS, ranked 2(nd) and 7(th) based on p-values from analyses of all invasive and invasive serous cases. The top SNP based on serous case analysis p-value (which ranked 197(th) for invasive case analysis), was ranked 8(th) based on the posterior probability of being in the top 5 markers (0.13). In summary, the application of the BLVM allows for the systematic integration of multiple SNP "features" for the prioritization of loci for fine-mapping or functional studies, taking into account the uncertainty in ranking.

We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk.

Defective cellular transport processes can lead to aberrant accumulation of trace elements, iron, small molecules and hormones in the cell, which in turn may promote the formation of reactive oxygen species, promoting DNA damage and aberrant expression of key regulatory cancer genes. As DNA damage and uncontrolled proliferation are hallmarks of cancer, including epithelial ovarian cancer (EOC), we hypothesized that inherited variation in the cellular transport genes contributes to EOC risk.

A substantial proportion of familial colorectal cancer (CRC) is not a consequence of known susceptibility loci, such as mismatch repair (MMR) genes, supporting the existence of additional loci. To identify novel CRC loci, we conducted a genome-wide linkage scan in 356 white families with no evidence of defective MMR (i.e., no loss of tumor expression of MMR proteins, no microsatellite instability (MSI)-high tumors, or no evidence of linkage to MMR genes). Families were ascertained via the Colon Cancer Family Registry multi-site NCI-supported consortium (Colon CFR), the City of Hope Comprehensive Cancer Center, and Memorial University of Newfoundland. A total of 1,612 individuals (average 5.0 per family including 2.2 affected) were genotyped using genome-wide single nucleotide polymorphism linkage arrays; parametric and non-parametric linkage analysis used MERLIN in a priori-defined family groups. Five lod scores greater than 3.0 were observed assuming heterogeneity. The greatest were among families with mean age of diagnosis less than 50 years at 4q21.1 (dominant HLOD = 4.51, α = 0.84, 145.40 cM, rs10518142) and among all families at 12q24.32 (dominant HLOD = 3.60, α = 0.48, 285.15 cM, rs952093). Among families with four or more affected individuals and among clinic-based families, a common peak was observed at 15q22.31 (101.40 cM, rs1477798; dominant HLOD = 3.07, α = 0.29; dominant HLOD = 3.03, α = 0.32, respectively). Analysis of families with only two affected individuals yielded a peak at 8q13.2 (recessive HLOD = 3.02, α = 0.51, 132.52 cM, rs1319036). These previously unreported linkage peaks demonstrate the continued utility of family-based data in complex traits and suggest that new CRC risk alleles remain to be elucidated.

FKBP51, (FKBP5), is a negative regulator of Akt. Variability in FKBP5 expression level is a major factor contributing to variation in response to chemotherapeutic agents including gemcitabine, a first line treatment for pancreatic cancer. Genetic variation in FKBP5 could influence its function and, ultimately, treatment response of pancreatic cancer.

In recent years, capabilities for genotyping large sets of single nucleotide polymorphisms (SNPs) has increased considerably with the ability to genotype over 1 million SNP markers across the genome. This advancement in technology has led to an increase in the number of genome-wide association studies (GWAS) for various complex traits. These GWAS have resulted in the implication of over 1500 SNPs associated with disease traits. However, the SNPs identified from these GWAS are not necessarily the functional variants. Therefore, the next phase in GWAS will involve the refining of these putative loci.

During the past decade, the age-adjusted mortality rate for endometrial cancer (EC) increased 1.9% annually with TP53 mutant (TP53-mu) EC disproportionally represented in advanced disease and deaths. Therefore, we aimed to assess pivotal molecular parameters that differentiate clinical outcomes in high- and low-risk EC. Using the Cancer Genome Atlas, we analyzed EC specimens with available DNA sequences and quantitative gene-specific RNA expression data. After polymerase ɛ (POLE)-mutant specimens were excluded, differential gene-specific mutations and mRNA expressions were annotated and integrated. Consequent to TP53-mu failure to induce p21, derepression of multiple oncogenes harboring promoter p21 repressive sites was observed, including CCNA2 and FOXM1 (P < .001 compared with TP53 wild type [TP53-wt]). TP53-wt EC with high CCNA2 expression (CCNA2-H) had a targeted transcriptomic profile similar to that of TP53-mu EC, suggesting CCNA2 is a seminal determinant for both TP53-wt and TP53-mu EC. CCNA2 enhances E2F1 function, upregulating FOXM1 and CIP2A, as observed in TP53-mu and CCNA2-H TP53-wt EC (P < .001). CIP2A inhibits protein phosphatase 2A, leading to AKT inactivation of GSK3β and restricted oncoprotein degradation; PPP2R1A and FBXW7 mutations yield similar results. Upregulation of FOXM1 and failed degradation of FOXM1 is evidenced by marked upregulation of multiple homologous recombination genes (P < .001). Integrating these molecular aberrations generated a molecular biomarker panel with significant prognostic discrimination (P = 5.8×10-7); adjusting for age, histology, grade, myometrial invasion, TP53 status, and stage, only CCNA2-H/E2F1-H (P = .0003), FBXW7-mu/PPP2R1A-mu (P = .0002), and stage (P = .017) were significant. The generated prognostic molecular classification system identifies dissimilar signaling aberrations potentially amenable to targetable therapeutic options.

Ovarian carcinosarcoma is a rare subtype of ovarian cancer with poor clinical outcomes. The low incidence of this disease makes accrual to large clinical trials challenging. However, studies have shown that treatment responses in patient-derived xenograft (PDX) models correlate with matched-patient responses in the clinic, supporting their use for preclinical testing of standard and novel therapies. An ovarian carcinosarcoma PDX is presented herein and showed resistance to carboplatin and paclitaxel (similar to the patient) but exhibited significant sensitivity to ifosfamide and paclitaxel. The PDX demonstrated overexpression of EGFR mRNA and gene amplification by array comparative genomic hybridization (log2 ratio 0.399). EGFR phosphorylation was also detected. Angiogensis and insulin-like growth factor pathways were also implicated by overexpression of VEGFC and IRS1. In order to improve response to chemotherapy, the PDX was treated with carboplatin/paclitaxel with or without a pan-HER and VEGF inhibitor (BMS-690514) but there was no tumor growth inhibition or improved animal survival, which may be explained by a KRAS mutation. Resistance was also observed when the IGF-1R inhibitor BMS-754807 was combined with carboplatin/paclitaxel. Because poly (ADP-ribose) polymerase inhibitors have activity in ovarian cancer patients, with and without BRCA mutations, ABT-888 was also tested but found to have no activity. Pathogenic mutations were also detected in TP53 and PIK3CA. In conclusion, ifosfamide/paclitaxel was superior to carboplatin/paclitaxel in this ovarian carcinosarcoma PDX and gene overexpression or amplification alone was not sufficient to predict response to targeted therapy. Better predictive markers of response are needed.

Recently, structural variation in the genome has been implicated in many complex diseases. Using genomewide single nucleotide polymorphism (SNP) arrays, researchers are able to investigate the impact not only of SNP variation, but also of copy-number variants (CNVs) on the phenotype. The most common analytic approach involves estimating, at the level of the individual genome, the underlying number of copies present at each location. Once this is completed, tests are performed to determine the association between copy number state and phenotype. An alternative approach is to carry out association testing first, between phenotype and raw intensities from the SNP array at the level of the individual marker, and then aggregate neighboring test results to identify CNVs associated with the phenotype. Here, we explore the strengths and weaknesses of these two approaches using both simulations and real data from a pharmacogenomic study of the chemotherapeutic agent gemcitabine. Our results indicate that pooled marker-level testing is capable of offering a dramatic increase in power (> 12-fold) over CNV-level testing, particularly for small CNVs. However, CNV-level testing is superior when CNVs are large and rare; understanding these tradeoffs is an important consideration in conducting association studies of structural variation.

We introduce a novel per-gene measure of intra-gene DNA methylation variability (IGV) based on the Illumina Infinium HumanMethylation450 platform, which is prognostic independently of well-known predictors of clinical outcome. Using IGV, we derive a robust gene-panel prognostic signature for ovarian cancer (OC, n = 221), which validates in two independent data sets from Mayo Clinic (n = 198) and TCGA (n = 358), with significance of p = 0.004 in both sets. The OC prognostic signature gene-panel is comprised of four gene groups, which represent distinct biological processes. We show the IGV measurements of these gene groups are most likely a reflection of a mixture of intra-tumour heterogeneity and transcription factor (TF) binding/activity. IGV can be used to predict clinical outcome in patients individually, providing a surrogate read-out of hard-to-measure disease processes.

Despite low incidence, ovarian cancer is the fifth leading cause of cancer deaths and it has the highest mortality rate of all gynecologic malignancies among US women. The mortality rate would be reduced with an early detection marker. The folate receptor alpha (FRalpha) is one logical choice for a biomarker because of its prevalent overexpression in ovarian cancer and its exclusive expression in only a few normal tissues. In prior work, it was observed that patients with ovarian cancer had elevated serum levels of a protein that bound to a FRalpha-specific monoclonal antibody relative to healthy individuals. However, it was not shown that the protein detected was intact functional FRalpha. In the current study, the goal was to determine whether ovarian cancer patients (n = 30) had elevated serum levels of a fully functional intact FRalpha compared to matched healthy controls (n = 30).

Although ovarian cancer is the most lethal of gynecologic malignancies, wide variation in outcome following conventional therapy continues to exist. The presence of tumor-infiltrating regulatory T cells (Tregs) has a role in outcome of this disease, and a growing body of data supports the existence of inherited prognostic factors. However, the role of inherited variants in genes encoding Treg-related immune molecules has not been fully explored. We analyzed expression quantitative trait loci (eQTL) and sequence-based tagging single nucleotide polymorphisms (tagSNPs) for 54 genes associated with Tregs in 3,662 invasive ovarian cancer cases. With adjustment for known prognostic factors, suggestive results were observed among rarer histological subtypes; poorer survival was associated with minor alleles at SNPs in RGS1 (clear cell, rs10921202, p=2.7×10(-5)), LRRC32 and TNFRSF18/TNFRSF4 (mucinous, rs3781699, p=4.5×10(-4), and rs3753348, p=9.0×10(-4), respectively), and CD80 (endometrioid, rs13071247, p=8.0×10(-4)). Fo0r the latter, correlative data support a CD80 rs13071247 genotype association with CD80 tumor RNA expression (p=0.006). An additional eQTL SNP in CD80 was associated with shorter survival (rs7804190, p=8.1×10(-4)) among all cases combined. As the products of these genes are known to affect induction, trafficking, or immunosuppressive function of Tregs, these results suggest the need for follow-up phenotypic studies.

Cancer patients show large individual variation in their response to chemotherapeutic agents. Gemcitabine (dFdC) and AraC, two cytidine analogues, have shown significant activity against a variety of tumors. We previously used expression data from a lymphoblastoid cell line-based model system to identify genes that might be important for the two drug cytotoxicity. In the present study, we used that same model system to perform a genome-wide association (GWA) study to test the hypothesis that common genetic variation might influence both gene expression and response to the two drugs. Specifically, genome-wide single nucleotide polymorphisms (SNPs) and mRNA expression data were obtained using the Illumina 550K(R) HumanHap550 SNP Chip and Affymetrix U133 Plus 2.0 GeneChip, respectively, for 174 ethnically-defined "Human Variation Panel" lymphoblastoid cell lines. Gemcitabine and AraC cytotoxicity assays were performed to obtain IC(50) values for the cell lines. We then performed GWA studies with SNPs, gene expression and IC(50) of these two drugs. This approach identified SNPs that were associated with gemcitabine or AraC IC(50) values and with the expression regulation for 29 genes or 30 genes, respectively. One SNP in IQGAP2 (rs3797418) was significantly associated with variation in both the expression of multiple genes and gemcitabine and AraC IC(50). A second SNP in TGM3 (rs6082527) was also significantly associated with multiple gene expression and gemcitabine IC50. To confirm the association results, we performed siRNA knock down of selected genes with expression that was associated with rs3797418 and rs6082527 in tumor cell and the knock down altered gemcitabine or AraC sensitivity, confirming our association study results. These results suggest that the application of GWA approaches using cell-based model systems, when combined with complementary functional validation, can provide insights into mechanisms responsible for variation in cytidine analogue response.

Genetic variation underlying the regulation of mRNA gene expression in humans may provide key insights into the molecular mechanisms of human traits and complex diseases. Current statistical methods to map genetic variation associated with mRNA gene expression have typically applied standard linkage and/or association methods; however, when genome-wide SNP and mRNA expression data are available performing all pair wise comparisons is computationally burdensome and may not provide optimal power to detect associations. Consideration of different approaches to account for the high dimensionality and multiple testing issues may provide increased efficiency and statistical power. Here we present a novel approach to model and test the association between genetic variation and mRNA gene expression levels in the context of gene sets (GSs) and pathways, referred to as gene set - expression quantitative trait loci analysis (GS-eQTL). The method uses GSs to initially group SNPs and mRNA expression, followed by the application of principal components analysis (PCA) to collapse the variation and reduce the dimensionality within the GSs. We applied GS-eQTL to assess the association between SNP and mRNA expression level data collected from a cell-based model system using PharmGKB and KEGG defined GSs. We observed a large number of significant GS-eQTL associations, in which the most significant associations arose between genetic variation and mRNA expression from the same GS. However, a number of associations involving genetic variation and mRNA expression from different GSs were also identified. Our proposed GS-eQTL method effectively addresses the multiple testing limitations in eQTL studies and provides biological context for SNP-expression associations.

Alterations in stromal tissue components can inhibit or promote epithelial tumorigenesis. Decorin (DCN) and lumican (LUM) show reduced stromal expression in serous epithelial ovarian cancer (sEOC). We hypothesized that common variants in these genes associate with risk. Associations with sEOC among Caucasians were estimated with odds ratios (OR) among 397 cases and 920 controls in two U.S.-based studies (discovery set), 436 cases and 1,098 controls in Australia (replication set 1) and a consortium of 15 studies comprising 1,668 cases and 4,249 controls (replication set 2). The discovery set and replication set 1 (833 cases and 2,013 controls) showed statistically homogeneous (P(heterogeneity)≥0.48) decreased risks of sEOC at four variants: DCN rs3138165, rs13312816 and rs516115, and LUM rs17018765 (OR = 0.6 to 0.9; P(trend) = 0.001 to 0.03). Results from replication set 2 were statistically homogeneous (P(heterogeneity)≥0.13) and associated with increased risks at DCN rs3138165 and rs13312816, and LUM rs17018765: all ORs = 1.2; P(trend)≤0.02. The ORs at the four variants were statistically heterogeneous across all 18 studies (P(heterogeneity)≤0.03), which precluded combining. In post-hoc analyses, interactions were observed between each variant and recruitment period (P(interaction)≤0.003), age at diagnosis (P(interaction) = 0.04), and year of diagnosis (P(interaction) = 0.05) in the five studies with available information (1,044 cases, 2,469 controls). We conclude that variants in DCN and LUM are not directly associated with sEOC, and that confirmation of possible effect modification of the variants by non-genetic factors is required.

Genetic variation at the TERT-CLPTM1L locus at 5p15.33 is associated with susceptibility to several cancers, including epithelial ovarian cancer (EOC). We have carried out fine-mapping of this region in EOC which implicates an association with a single nucleotide polymorphism (SNP) within the TERT promoter. We demonstrate that the minor alleles at rs2736109, and at an additional TERT promoter SNP, rs2736108, are associated with decreased breast cancer risk, and that the combination of both SNPs substantially reduces TERT promoter activity.

Recent studies have identified several single nucleotide polymorphisms (SNPs) in the population that are associated with variations in the risks of many different diseases including cancers such as breast, prostate and colorectal. For ovarian cancer, the known highly penetrant susceptibility genes (BRCA1 and BRCA2) are probably responsible for only 40% of the excess familial ovarian cancer risks, suggesting that other susceptibility genes of lower penetrance exist.