PhenomeXcan: Mapping the genome to the phenome through the transcriptome.

Large-scale genomic and transcriptomic initiatives offer unprecedented insight into complex traits, but clinical translation remains limited by variant-level associations without biological context and lack of analytic resources. Our resource, PhenomeXcan, synthesizes 8.87 million variants from genome-wide association study summary statistics on 4091 traits with transcriptomic data from 49 tissues in Genotype-Tissue Expression v8 into a gene-based, queryable platform including 22,515 genes. We developed a novel Bayesian colocalization method, fast enrichment estimation aided colocalization analysis (fastENLOC), to prioritize likely causal gene-trait associations. We successfully replicate associations from the phenome-wide association studies (PheWAS) catalog Online Mendelian Inheritance in Man, and an evidence-based curated gene list. Using PhenomeXcan results, we provide examples of novel and underreported genome-to-phenome associations, complex gene-trait clusters, shared causal genes between common and rare diseases via further integration of PhenomeXcan with ClinVar, and potential therapeutic targets. PhenomeXcan (phenomexcan.org) provides broad, user-friendly access to complex data for translational researchers.

Fine-mapping and QTL tissue-sharing information improves the reliability of causal gene identification.

The integration of transcriptomic studies and genome-wide association studies (GWAS) via imputed expression has seen extensive application in recent years, enabling the functional characterization and causal gene prioritization of GWAS loci. However, the techniques for imputing transcriptomic traits from DNA variation remain underdeveloped. Furthermore, associations found when linking eQTL studies to complex traits through methods like PrediXcan can lead to false positives due to linkage disequilibrium between distinct causal variants. Therefore, the best prediction performance models may not necessarily lead to more reliable causal gene discovery. With the goal of improving discoveries without increasing false positives, we develop and compare multiple transcriptomic imputation approaches using the most recent GTEx release of expression and splicing data on 17,382 RNA-sequencing samples from 948 post-mortem donors in 54 tissues. We find that informing prediction models with posterior causal probability from fine-mapping (dap-g) and borrowing information across tissues (mashr) can lead to better performance in terms of number and proportion of significant associations that are colocalized and the proportion of silver standard genes identified as indicated by precision-recall and receiver operating characteristic curves. All prediction models are made publicly available at predictdb.org.

Protein quantitative trait loci identify novel candidates modulating cellular response to chemotherapy.

Annotating and interpreting the results of genome-wide association studies (GWAS) remains challenging. Assigning function to genetic variants as expression quantitative trait loci is an expanding and useful approach, but focuses exclusively on mRNA rather than protein levels. Many variants remain without annotation. To address this problem, we measured the steady state abundance of 441 human signaling and transcription factor proteins from 68 Yoruba HapMap lymphoblastoid cell lines to identify novel relationships between inter-individual protein levels, genetic variants, and sensitivity to chemotherapeutic agents. Proteins were measured using micro-western and reverse phase protein arrays from three independent cell line thaws to permit mixed effect modeling of protein biological replicates. We observed enrichment of protein quantitative trait loci (pQTLs) for cellular sensitivity to two commonly used chemotherapeutics: cisplatin and paclitaxel. We functionally validated the target protein of a genome-wide significant trans-pQTL for its relevance in paclitaxel-induced apoptosis. GWAS overlap results of drug-induced apoptosis and cytotoxicity for paclitaxel and cisplatin revealed unique SNPs associated with the pharmacologic traits (at p<0.001). Interestingly, GWAS SNPs from various regions of the genome implicated the same target protein (p<0.0001) that correlated with drug induced cytotoxicity or apoptosis (p ≤ 0.05). Two genes were functionally validated for association with drug response using siRNA: SMC1A with cisplatin response and ZNF569 with paclitaxel response. This work allows pharmacogenomic discovery to progress from the transcriptome to the proteome and offers potential for identification of new therapeutic targets. This approach, linking targeted proteomic data to variation in pharmacologic response, can be generalized to other studies evaluating genotype-phenotype relationships and provide insight into chemotherapeutic mechanisms.

EPS8 inhibition increases cisplatin sensitivity in lung cancer cells.

Cisplatin, a commonly used chemotherapeutic, is associated with ototoxicity, renal toxicity and neurotoxicity, thus identifying means to increase the therapeutic index of cisplatin may allow for improved outcomes. A SNP (rs4343077) within EPS8, discovered through a genome wide association study of cisplatin-induced cytotoxicity and apoptosis in lymphoblastoid cell lines (LCLs), provided impetus to further study this gene. The purpose of this work was to evaluate the role of EPS8 in cellular susceptibility to cisplatin in cancerous and non-cancerous cells. We used EPS8 RNA interference to determine the effect of decreased EPS8 expression on LCL and A549 lung cancer cell sensitivity to cisplatin. EPS8 knockdown in LCLs resulted in a 7.9% increase in cisplatin-induced survival (P = 1.98 × 10(-7)) and an 8.7% decrease in apoptosis (P = 0.004) compared to control. In contrast, reduced EPS8 expression in lung cancer cells resulted in a 20.6% decrease in cisplatin-induced survival (P = 5.08 × 10(-5)). We then investigated an EPS8 inhibitor, mithramycin A, as a potential agent to increase the therapeutic index of cisplatin. Mithramycin A decreased EPS8 expression in LCLs resulting in decreased cellular sensitivity to cisplatin as evidenced by lower caspase 3/7 activation following cisplatin treatment (42.7% ± 6.8% relative to control P = 0.0002). In 5 non-small-cell lung carcinoma (NSCLC) cell lines, mithramycin A also resulted in decreased EPS8 expression. Adding mithramycin to 4 NSCLC cell lines and a bladder cancer cell line, resulted in increased sensitivity to cisplatin that was significantly more pronounced in tumor cell lines than in LCL lines (p<0.0001). An EGFR mutant NSCLC cell line (H1975) showed no significant change in sensitivity to cisplatin with the addition of mithramycin treatment. Therefore, an inhibitor of EPS8, such as mithramycin A, could improve cisplatin treatment by increasing sensitivity of tumor relative to normal cells.

Neo-adjuvant capecitabine chemotherapy in women with newly diagnosed locally advanced breast cancer in a resource-poor setting (Nigeria): efficacy and safety in a phase II feasibility study.

The majority of clinical trials of neo-adjuvant therapy for breast cancer have been conducted in resource-rich countries. We chose Nigeria, a resource-poor country, as the major site for a phase II feasibility open-label multicenter clinical trial designed to evaluate the efficacy, safety, and tolerability of neo-adjuvant capecitabine in locally advanced breast cancer (LABC). Planned treatment consisted of 24 weeks of capecitabine at a dose of 1,000 mg/m(2) twice daily (2,000 mg/m(2) total per day). The primary endpoints were overall, partial, complete clinical response rate (OCR, PCR, CCR) and complete pathologic response (cPR). A total of 16 patients were recruited from August 2007 to April 2010. The study was terminated early as a result of slow accrual. After the first three cycles of therapy, PCR were seen in five of 16 patients (31%; 95% CI 11-59%). Of the remaining 11 patients, eight had no response (NR) or stable disease (SD), and three had progressive disease (PD). Seven patients proceeded with further therapy of which had SD. OCR at the end of eight cycles was 44% (95% CI 20-70%). Clinical response and radiologic response by ultrasonomammography were highly concordant (spearman correlation 0.70). The most common adverse effect was Grade 1 hand-foot syndrome, which was seen in 75% of patients. Despite several limitations, we successfully carried out this phase II feasibility study of neo-adjuvant capecitabine for LABC in Nigeria. Capecitabine monotherapy showed good overall response rates with minimal toxicity and further studies are warranted.

Comprehensive genetic analysis of cytarabine sensitivity in a cell-based model identifies polymorphisms associated with outcome in AML patients.

A whole-genome approach was used to investigate the genetic determinants of cytarabine-induced cytotoxicity. We performed a meta-analysis of genome-wide association studies involving 523 lymphoblastoid cell lines (LCLs) from individuals of European, African, Asian, and African American ancestry. Several of the highest-ranked single-nucleotide polymorphisms (SNPs) were within the mutated in colorectal cancers (MCC) gene. MCC expression was induced by cytarabine treatment from 1.7- to 26.6-fold in LCLs. A total of 33 SNPs ranked at the top of the meta-analysis (P < 10(-5)) were successfully tested in a clinical trial of patients randomized to receive low-dose or high-dose cytarabine plus daunorubicin and etoposide; of these, 18 showed association (P < .05) with either cytarabine 50% inhibitory concentration in leukemia cells or clinical response parameters (minimal residual disease, overall survival (OS), and treatment-related mortality). This count (n = 18) was significantly greater than expected by chance (P = .016). For rs1203633, LCLs with AA genotype were more sensitive to cytarabine-induced cytotoxicity (P = 1.31 × 10(-6)) and AA (vs GA or GG) genotype was associated with poorer OS (P = .015), likely as a result of greater treatment-related mortality (P = .0037) in patients with acute myeloid leukemia (AML). This multicenter AML02 study trial was registered at www.clinicaltrials.gov as #NCT00136084.