Statistical Genetics of Dose Response Traits

National Institute Of Environmental Health Sciences

Investigators

Linked publications, trials & patents

Abstract

Dissecting the genetic etiology of cancer drug response using a lymphoblastoid cell line model has been a longstanding research goal. The use of ex-vivo model systems to provide a level of forecasting for in-vivo characteristics remains an important need for cancer therapeutics. The use of lymphoblastoid cell lines (LCLs) is an attractive approach for pharmacogenomics and toxicogenomics, due to their scalability, efficiency, and cost-effectiveness. In previous years of the project, we used LCLs to conduct genome-wide association studies for numerous anti-cancer drugs. The GWAS analyses identified a number of significant findings, and this year the project was focused on functional follow up studies to confirm associations and begin to understand the functional consequences of genetic variation in this context. Summaries of these studies are below. Oxaliplatin (OXAL) is a commonly used chemotherapy for treating colorectal cancer (CRC). A recent genome wide association study (GWAS) showed that a genetic variant (rs11006706) in the lncRNA gene MKX-AS1 and partnered sense gene MKX could impact the response of genetically varied cell lines to OXAL treatment. This study found that the expression levels of MKX-AS1 and MKX in lymphocytes (LCLs) and CRC cell lines differed between the rs11006706 genotypes, indicating that this gene pair could play a role in OXAL response. Further analysis of patient survival data from the Cancer Genome Atlas (TCGA) and other sources showed that patients with high MKX-AS1 expression status had significantly worse overall survival (HR = 3.2; 95%CI = (1.179); p = 0.024) compared to cases with low MKX-AS1 expression status. Alternatively, high MKX expression status had significantly better overall survival (HR = 0.22; 95%CI = (0.070.7); p = 0.01) compared to cases with low MKX expression status. These results suggest an association between MKX-AS1 and MKX expression status that could be useful as a prognostic marker of response to OXAL and potential patient outcomes in CRC. Temozolomide (TMZ) chemotherapy is an important tool in the treatment of glioma brain tumors. However, variable patient response and chemo-resistance remain exceptionally challenging. Our previous genome-wide association study (GWAS) identified a suggestively significant association of SNP rs4470517 in the RYK (receptor-like kinase) gene with TMZ drug response. Functional validation of RYK using lymphocytes and glioma cell lines resulted in gene expression analysis indicating differences in expression status between genotypes of the cell lines and TMZ dose response. We conducted univariate and multivariate Cox regression analyses using publicly available TCGA and GEO datasets to investigate the impact of RYK gene expression status on glioma patient overall (OS) and progression-free survival (PFS). Our results indicated that in IDH mutant gliomas, RYK expression and tumor grade were significant predictors of survival. In IDH wildtype glioblastomas (GBM), MGMT status was the only significant predictor. Despite this result, we revealed a potential benefit of RYK expression in IDH wildtype GBM patients. We found that a combination of RYK expression and MGMT status could serve as an additional biomarker for improved survival. Overall, our findings suggest that RYK expression may serve as an important prognostic or predictor of TMZ response and survival for glioma patients. Monoclonal antibody (mAb) therapy directed against CD20 is an important tool in the treatment of B cell disorders. However, variable patient response and acquired resistance remain important clinical challenges. To identify genetic factors that may influence sensitivity to treatment, the cytotoxic activity of three CD20 mAbs: rituximab; ofatumumab; and obinutuzumab, were screened in high-throughput assays using 680 ethnically diverse lymphoblastoid cell lines (LCLs) followed by a pharmacogenomic assessment. GWAS analysis identified several novel gene candidates. The most significant SNP, rs58600101, in the gene MKL1 displayed ethnic stratification, with the variant being significantly more prevalent in the African cohort and resulting in reduced transcript levels as measured by qPCR. Functional validation of MKL1 by shRNA-mediated knockdown of MKL1 resulted in a more resistant phenotype. Gene expression analysis identified the developmentally associated TGFB1I1 as the most significant gene associated with sensitivity. qPCR among a panel of sensitive and resistant LCLs revealed immunoglobulin class-switching as well as differences in the expression of B cell activation markers. Flow cytometry showed heterogeneity within some cell lines relative to surface Ig isotype with a shift to more IgG+ cells among the resistant lines. Pretreatment with prednisolone could partly reverse the resistant phenotype. Results suggest that the efficacy of anti-CD20 mAb therapy may be influenced by B cell developmental status as well as polymorphism in the MKL1 gene