Im, Hae Kyung

Abstract RATIONALE: Chronic Obstructive Pulmonary Disease (COPD) is a leading cause of mortality worldwide. We previously published the polygenic transcriptome risk score (PTRS) which uses the cumulative effect of predicted gene expression to construct genetic predictors for complex diseases, and we demonstrated the value of PTRS for COPD built on Genotype-Tissue Expression (GTEx) Lung tissue with significantly improved cross-ancestry portability. However, the postmortem collection and lack of lung disease status on GTEx samples likely affect results. We hypothesized that the performance of PTRS will be improved by using disease-relevant expression quantitative trait loci (eQTLs) to construct the score. Here, we aimed to improve performance of PTRS by leveraging disease-relevant eQTLs. METHODS: We constructed two transcriptome prediction models, Elastic Net (EN) and Prediction Using Models Informed by Chromatin conformation and Epigenomics (PUMICE), using TOPMed lung RNA-seq data from the Lung Tissue Research Consortium (LTRC), which included COPD patients undergoing lung surgery. We compared their performance with two published models built on RNA-seq from GTEx Lung tissue, GTEx-EN and GTEx-PUMICE. We first integrated each of four transcriptome models with multi-ancestry GWAS of FEV1/FVC ratio (Shrine et al. 2023) to generate transcriptome-wide association study (TWAS) results, which produced trait-associated genes. Pathway analysis was then conducted on significant TWAS genes. Finally, the PTRS was computed on TWAS genes that were included in the top 20 significantly enriched pathways. We tested performance of PTRS in race-stratified analysis for moderate-to-severe and severe COPD in COPDGene. RESULTS: Our models produced more significant TWAS genes (FDR<0.05) and had more overlaps with genes identified for FEV1/FVC ratio from Shrine et al. 2023 compared to two GTEx models (Fig. A). Additionally, our models had significantly higher AUC for predicting both moderate-to-severe and severe COPD in race-stratified analysis. The mean AUC from LTRC-EN for predicting severe COPD was 0.57 and 0.55 respectively for EUR and AFA, which was significantly higher than mean AUC of 0.54 and 0.51 from GTEx-EN (Delong p-value=8.15x10-10 for EUR and 1.42x10-3 for AFA) (Fig. B). Similarly, the PTRS derived from our models produced higher Odds Ratios (OR) per standard deviation. The mean OR from LTRC-PUMICE for moderate-to-severe COPD was 1.25 and 1.18 respectively for EUR and AFA, while it was 1.17 and 1.12 from GTEx-PUMICE (Fig. C). CONCLUSIONS: Our study demonstrates the value of leveraging disease-relevant RNA-seq to construct transcriptome prediction models and their improvement on identifying lung function genes and predicting COPD across ancestries.

Abstract Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10−6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10−7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10−3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis. Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419–30. ©2018 AACR.

Abstract Background: Large-scale genome-wide association studies (GWAS) have identified ~35 loci associated with epithelial ovarian cancer risk. However, these genetic loci together explain only a small portion of the heritability of this malignancy. The large majority of the GWAS-identified variants are located in non-coding regions, thus possible causal genes in these loci remain largely unknown. We performed a transcriptome-wide association study (TWAS) to search for novel genetic loci for ovarian cancer risk and plausible causal genes at GWAS-identified loci. Method: Transcriptome data from normal ovarian tissue samples (n=68) and all tissue samples (n = 369), along with their high-density genotyping data, obtained from only European descendants included in the Genotype-Tissue Expression Project (GTEx), were used to build ovarian linear prediction models and cross-tissue models (to increase statistical power) using the elastic net method. Based on model performance, we evaluated 17,121 genes for their cis-predicted gene expressions in relation to ovarian cancer risk using summary statistics data generated in GWAS of ovarian cancer from 25,509 cases and 40,941 controls. MetaXcan was used to integrate gene expression prediction model with summary statistics. Results: We identified 35 genes with predicted expression levels associated with ovarian cancer risk at P value < 2.2 × 10-6, the Bonferroni corrected significance level for multiple comparisons. Of these, 12 genes at 4 genetic loci are located ≥500kb away from risk SNPs previously reported in GWAS, representing potential novel genetic loci for ovarian cancer risk. The remaining 23 genes at 12 loci are located within known ovarian cancer risk loci. Fifteen of these genes at the 12 loci have not been reported in previous studies. Analyses by ovarian cancer histological subtypes showed that the majority of these 35 genes are associated with serous invasive carcinoma. Several new associations were identified in histological subtype analyses. Conclusion: In this TWAS we identified multiple genes with predicted expressions related to ovarian cancer risk and provide substantial new information to enhance the understanding of ovarian cancer biology and genetics. Citation Format: Yingchang Lu, Joellen M. Schildkraut, Thomas A. Sellers, Lang Wu, Xingyi Guo, Bingshan Li, Y. Ann Chen, Jennifer B. Doherty, Simon Gayther, Ellen L. Goode, Hae Kyung Im, Siddhartha Kar, Kate Lawrenson, Ani W. Manichaikul, Jennifer B. Permuth, Brett M. Reid, Jamie K. Teer, Paul Pharoah, Wei Zheng, Jirong Long. Transcriptome-wide association study among 66,450 women to identify candidate susceptible genes for ovarian cancer risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1308. doi:10.1158/1538-7445.AM2017-1308

Abstract Purpose: Cell-based approaches were used to identify genetic markers predictive of patients' risk for poor response prior to chemotherapy. Experimental Design: We conducted genome-wide association studies (GWAS) to identify single-nucleotide polymorphisms (SNP) associated with cellular sensitivity to carboplatin through their effects on mRNA expression using International HapMap lymphoblastoid cell lines (LCL) and replicated them in additional LCLs. SNPs passing both stages of the cell-based study were tested for association with progression-free survival (PFS) in patients. Phase 1 validation was based on 377 ovarian cancer patients receiving at least four cycles of carboplatin and paclitaxel from the Australian Ovarian Cancer Study (AOCS). Positive associations were then assessed in phase 2 validation analysis of 1,326 patients from the Ovarian Cancer Association Consortium and The Cancer Genome Atlas. Results: In the initial GWAS, 342 SNPs were associated with carboplatin-induced cytotoxicity, of which 18 unique SNPs were retained after assessing their association with gene expression. One SNP (rs1649942) was replicated in an independent LCL set (Bonferroni adjusted P < 0.05). It was found to be significantly associated with decreased PFS in phase 1 AOCS patients (Pper-allele = 2 × 10−2), with a stronger effect in the subset of women with optimally debulked tumors (Pper-allele = 4 × 10−3). rs1649942 was also associated with poorer overall survival in women with optimally debulked tumors (Pper-allele = 9 × 10−3). However, this SNP was not significant in phase 2 validation analysis with patients from numerous cohorts. Conclusion: This study shows the potential of cell-based, genome-wide approaches to identify germline predictors of treatment outcome and highlights the need for extensive validation in patients to assess their clinical effect. Clin Cancer Res; 17(16); 5490–500. ©2011 AACR.