Nifang Niu

Radiation pharmacogenomics: A genome-wide association approach to identify radiation response biomarkers using human lymphoblastoid cell lines

Radiation therapy is used to treat half of all cancer patients. Response to radiation therapy varies widely among patients. Therefore, we performed a genome-wide association study (GWAS) to identify biomarkers to help predict radiation response using 277 ethnically defined human lymphoblastoid cell lines (LCLs). Basal gene expression levels and 1.3 million genome-wide single nucleotide polymorphism (SNP) markers from both Affymetrix and Illumina platforms were assayed for all 277 human LCLs. MTS [3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium] assays for radiation cytotoxicity were also performed to obtain area under the curve (AUC) as a radiation response phenotype for use in the association studies. Functional validation of candidate genes, selected from an integrated analysis that used SNP, expression, and AUC data, was performed with multiple cancer cell lines using specific siRNA knockdown, followed by MTS and colony-forming assays. A total of 27 loci, each containing at least two SNPs within 50 kb with P-values less than 10(-4) were associated with radiation AUC. A total of 270 expression probe sets were associated with radiation AUC with P < 10(-3). The integrated analysis identified 50 SNPs in 14 of the 27 loci that were associated with both AUC and the expression of 39 genes, which were also associated with radiation AUC (P < 10(-3)). Functional validation using siRNA knockdown in multiple tumor cell lines showed that C13orf34, MAD2L1, PLK4, TPD52, and DEPDC1B each significantly altered radiation sensitivity in at least two cancer cell lines. Studies performed with LCLs can help to identify novel biomarkers that might contribute to variation in response to radiation therapy and enhance our understanding of mechanisms underlying that variation.

Genetic variation predicting cisplatin cytotoxicity associated with overall survival in lung cancer patients receiving platinum-based chemotherapy ,

Purpose: Inherited variability in the prognosis of lung cancer patients treated with platinum-based chemotherapy has been widely investigated. However, the overall contribution of genetic variation to platinum response is not well established. To identify novel candidate single nucleotide polymorphisms (SNP)/genes, we carried out a genome-wide association study (GWAS) for cisplatin cytotoxicity by using lymphoblastoid cell lines (LCL), followed by an association study of selected SNPs from the GWAS with overall survival (OS) in lung cancer patients. Experimental Design: A GWAS for cisplatin was conducted with 283 ethnically diverse LCLs. A total of 168 top SNPs were genotyped in 222 small cell lung cancer (SCLC) and 961 non-SCLC (NSCLC) patients treated with platinum-based therapy. Association of the SNPs with OS was determined by using the Cox regression model. Selected candidate genes were functionally validated by siRNA knockdown in human lung cancer cells. Results: Among 157 successfully genotyped SNPs, 9 and 10 SNPs were top SNPs associated with OS for patients with NSCLC and SCLC, respectively, although they were not significant after adjusting for multiple testing. Fifteen genes, including 7 located within 200 kb up or downstream of the 4 top SNPs and 8 genes for which expression was correlated with 3 SNPs in LCLs were selected for siRNA screening. Knockdown of DAPK3 and METTL6, for which expression levels were correlated with the rs11169748 and rs2440915 SNPs, significantly decreased cisplatin sensitivity in lung cancer cells. Conclusions: This series of clinical and complementary laboratory-based functional studies identified several candidate genes/SNPs that might help predict treatment outcomes for platinum-based therapy of lung cancer. Clin Cancer Res; 17(17); 5801–11. ©2011 AACR.

Discovery of genetic biomarkers contributing to variation in drug response of cytidine analogues using human lymphoblastoid cell lines

BackgroundTwo cytidine analogues, gemcitabine and cytosine arabinoside (AraC), are widely used in the treatment of a variety of cancers with a large individual variation in response. To identify potential genetic biomarkers associated with response to these two drugs, we used a human lymphoblastoid cell line (LCL) model system with extensive genomic data, including 1.3 million SNPs and 54,000 basal expression probesets to perform genome-wide association studies (GWAS) with gemcitabine and AraC IC50 values.ResultsWe identified 11 and 27 SNP loci significantly associated with gemcitabine and AraC IC50 values, respectively. Eleven candidate genes were functionally validated using siRNA knockdown approach in multiple cancer cell lines. We also characterized the potential mechanisms of genes by determining their influence on the activity of 10 cancer-related signaling pathways using reporter gene assays. Most SNPs regulated gene expression in a trans manner, except 7 SNPs in the PIGB gene that were significantly associated with both the expression of PIGB and gemcitabine cytotoxicity.ConclusionThese results suggest that genetic variation might contribute to drug response via either cis- or trans- regulation of gene expression. GWAS analysis followed by functional pharmacogenomics studies might help identify novel biomarkers contributing to variation in response to these two drugs and enhance our understanding of underlying mechanisms of drug action.

The eSNV-detect: a computational system to identify expressed single nucleotide variants from transcriptome sequencing data

Rapid development of next generation sequencing technology has enabled the identification of genomic alterations from short sequencing reads. There are a number of software pipelines available for calling single nucleotide variants from genomic DNA but, no comprehensive pipelines to identify, annotate and prioritize expressed SNVs (eSNVs) from non-directional paired-end RNA-Seq data. We have developed the eSNV-Detect, a novel computational system, which utilizes data from multiple aligners to call, even at low read depths, and rank variants from RNA-Seq. Multi-platform comparisons with the eSNV-Detect variant candidates were performed. The method was first applied to RNA-Seq from a lymphoblastoid cell-line, achieving 99.7% precision and 91.0% sensitivity in the expressed SNPs for the matching HumanOmni2.5 BeadChip data. Comparison of RNA-Seq eSNV candidates from 25 ER+ breast tumors from The Cancer Genome Atlas (TCGA) project with whole exome coding data showed 90.6–96.8% precision and 91.6–95.7% sensitivity. Contrasting single-cell mRNA-Seq variants with matching traditional multicellular RNA-Seq data for the MD-MB231 breast cancer cell-line delineated variant heterogeneity among the single-cells. Further, Sanger sequencing validation was performed for an ER+ breast tumor with paired normal adjacent tissue validating 29 out of 31 candidate eSNVs. The source code and user manuals of the eSNV-Detect pipeline for Sun Grid Engine and virtual machine are available at http://bioinformaticstools.mayo.edu/research/esnv-detect/.

HEATR1 negatively regulates akt to help sensitize pancreatic cancer cells to chemotherapy

Elucidating mechanisms of chemoresistance is critical to improve cancer therapy, especially for the treatment of pancreatic ductal adenocarcinoma (PDAC). Genome-wide association studies have suggested the less studied gene HEAT repeat-containing protein 1 (HEATR1) as a possible determinant of cellular sensitivity to different chemotherapeutic drugs. In this study, we assessed this hypothesized link in PDAC, where HEATR1 expression is downregulated significantly. HEATR1 silencing in PDAC cells increased resistance to gemcitabine and other chemotherapeutics, where this effect was associated with increased AKT kinase phosphorylation at the Thr308 regulatory site. Mechanistically, HEATR1 enhanced cell responsiveness to gemcitabine by acting as a scaffold to facilitate interactions between AKT and the protein phosphatase PP2A, thereby promoting Thr308 dephosphorylation. Consistent with these findings, treatment with the AKT inhibitor triciribine sensitized HEATR1-depleted PDAC cells to gemcitabine, suggesting that this therapeutic combination may overcome gemcitabine resistance in patients with low HEATR1 expression. Clinically, we found that HEATR1 downregulation in PDAC patients was associated with increased AKT phosphorylation, poor response to tumor resection plus gemcitabine standard-of-care treatment, and shorter overall survival. Collectively, our findings establish HEATR1 as a novel regulator of AKT and a candidate predictive and prognostic indicator of drug responsiveness and outcome in PDAC patients.

Inflammation‐Related Genetic Variants Predict Toxicity Following Definitive Radiotherapy for Lung Cancer

Definitive radiotherapy improves locoregional control and survival in inoperable non–small cell lung cancer patients. However, radiation‐induced toxicities (pneumonitis/esophagitis) are common dose‐limiting inflammatory conditions. We therefore conducted a pathway‐based analysis to identify inflammation‐related single‐nucleotide polymorphisms associated with radiation‐induced pneumonitis or esophagitis. A total of 11,930 single‐nucleotide polymorphisms were genotyped in 201 stage I–III non–small cell lung cancer patients treated with definitive radiotherapy. Validation was performed in an additional 220 non–small cell lung cancer cases. After validation, 19 single‐nucleotide polymorphisms remained significant. A polygenic risk score was generated to summarize the effect from validated single‐nucleotide polymorphisms. Significant improvements in discriminative ability were observed when the polygenic risk score was added into the clinical/epidemiological variable‐based model. We then used 277 lymphoblastoid cell lines to assess radiation sensitivity and expression quantitative trait loci (eQTL) relationships of the identified single‐nucleotide polymorphisms. Three genes (PRKCE, DDX58, and TNFSF7) were associated with radiation sensitivity. We concluded that inflammation‐related genetic variants could contribute to the development of radiation‐induced toxicities.

Metformin pharmacogenomics: a genome-wide association study to identify genetic and epigenetic biomarkers involved in metformin anticancer response using human lymphoblastoid cell lines

Metformin is currently considered as a promising anticancer agent in addition to its anti-diabetic effect. To better individualize metformin therapy and explore novel molecular mechanisms in cancer treatment, we conducted a pharmacogenomic study using 266 lymphoblastoid cell lines (LCLs). Metformin cytotoxicity assay was performed using the MTS assay. Genome-wide association (GWA) analyses were performed in LCLs using 1.3 million SNPs, 485k DNA methylation probes, 54k mRNA expression probe sets, and metformin cytotoxicity (IC50s). Top candidate genes were functionally validated using siRNA screening, followed by MTS assay in breast cancer cell lines. Further study of one top candidate, STUB1, was performed to elucidate the mechanisms by which STUB1 might contribute to metformin action. GWA analyses in LCLs identified 198 mRNA expression probe sets, 12 SNP loci, and 5 DNA methylation loci associated with metformin IC50 with P-values <10−4 or <10−5. Integrated SNP/methylation loci-expression-IC50 analyses found 3 SNP loci or 5 DNA methylation loci associated with metformin IC50 through trans-regulation of expression of 11 or 26 genes with P-value <10−4. Functional validation of top 61 candidate genes in 4 IPA networks indicated down regulation of 14 genes significantly altered metformin sensitivity in two breast cancer cell lines. Mechanistic studies revealed that the E3 ubiquitin ligase, STUB1, could influence metformin response by facilitating proteasome-mediated degradation of cyclin A. GWAS using a genomic data-enriched LCL model system, together with functional and mechanistic studies using cancer cell lines, help us to identify novel genetic and epigenetic biomarkers involved in metformin anticancer response.

Abstract 2271: Metformin pharmacogenomics: A genome-wide associate study to identify genetic and epigenetic biomarkers involved in metformin response.

Proceedings: AACR 104th Annual Meeting 2013; Apr 6-10, 2013; Washington, DC BACKGROUND Metformin, a widely used anti-diabetic drug, is being considered as a highly promising agent for treatment and prevention of many types of cancer, including breast cancer. Although the exact mechanism of action of metformin is unknown, it is thought that metformin is an AMPK activator. To better individualize metformin therapy and explore additional underlying molecular mechanism associated with metformin action, we conducted a metformin pharmacogenomic study using 266 genomic data-enriched lymphoblastoid cell lines (LCLs). METHODS Genome-wide microarray data were generated for all LCLs, including 1.3 million SNPs, 485K DNA methylation probes and 54K mRNA expression probe sets. Metformin cytotoxicity assay was performed using MTS assay. Genome-wide association (GWA) analyses and integrated analyses using all the genomic data and metformin cytotoxicity (IC50 value) were performed to identify top candidate genes for functional validation in breast cancer cell lines, MDA-MB 231 and Hs578T, using siRNA screening approach. Further mechanistic study was performed to investigate the role of candidates in regulation of AMPK activity. RESULTS GWA analyses in LCLs identified 198 mRNA expression probe sets and 210 DNA methylation probes associated with metformin IC50 with p-value <10-4, respectively. Analysis of 1.3 million SNPs found 12 loci (a region containing at least 1 SNP of p-value <10-5 and 1 SNP of p-value <10-3 within 50kb) that were associated with metformin IC50. Integrated SNP loci-mRNA expression-IC50 analysis indicated that the SNP loci on chromosome 16 was associated with metformin IC50 through a trans-regulation of expression of 7 genes with p-value <10-4. Integrated methylation-mRNA expression-IC50 analysis showed that 15 DNA methylation probes in 6 genes were associated with metformin IC50 through both cis- and trans-regulation of expression of 48 genes with p-value <10-4, including 4 methylation probes in 3 of those 6 genes that were cis-correlated with its own gene expression. Ingenuity pathway analysis of top 62 candidate genes identified four major pathways and 55 genes within those pathways were selected for functional validation in breast cancer cell lines. The knockdown experiment showed that down regulation of 25 genes significantly altered metformin sensitivity in MDA-MB231, and 14 of them also showed the same effect in Hs578T. Preliminary mechanistic experiment indicated that knockdown of 7 genes significantly reduced AMPK activity in MDA-MB231 and 3 of them had same effect in Hs578T. CONCLUSIONS GWAS using a genomic data-enriched LCL model system, together with functional validation in breast cancer cell lines, could help us to identify novel genetic and epigenetic biomarkers involved in metformin response and help us to better understand the mechanisms of metformin in cancer treatment. Citation Format: Nifang Niu, Xianglin Tan, Brooke L. Fridley, Daniel J. Schaid, Ryan P. Abo, Anthony Batzler, Erin E. Carlson, Gregory Jenkins, Sebastian Moran, Holger Andrea Heyn, Manel Esteller Badosa, Liewei Wang. Metformin pharmacogenomics: A genome-wide associate study to identify genetic and epigenetic biomarkers involved in metformin response. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2271. doi:10.1158/1538-7445.AM2013-2271

TSPYL Family Regulates CYP17A1 and CYP3A4 Expression: Potential Mechanism Contributing to Abiraterone Response in Metastatic Castration-Resistant Prostate Cancer

The testis‐specific Y‐encoded‐like protein (TSPYL) gene family includes TSPYL1 to TSPYL6. We previously reported that TSPYL5 regulates cytochrome P450 (CYP) 19A1 expression. Here we show that TSPYLs, especially TSPYL 1, 2, and 4, can regulate the expression of many CYP genes, including CYP17A1, a key enzyme in androgen biosynthesis, and CYP3A4, an enzyme that catalyzes the metabolism of abiraterone, a CYP17 inhibitor. Furthermore, a common TSPYL1 single nucleotide polymorphism (SNP), rs3828743 (G/A) (Pro62Ser), abolishes TSPYL1s ability to suppress CYP3A4 expression, resulting in reduced abiraterone concentrations and increased cell proliferation. Data from a prospective clinical trial of 87 metastatic castration‐resistant prostate cancer patients treated with abiraterone acetate/prednisone showed that the variant SNP genotype (A) was significantly associated with worse response and progression‐free survival. In summary, TSPYL genes are novel CYP gene transcription regulators, and genetic alteration within these genes significantly influences response to drug therapy through transcriptional regulation of CYP450 genes.

Identification of genetic variants or genes that are associated with Homoharringtonine (HHT) response through a genome-wide association study in human lymphoblastoid cell lines (LCLs)

Homoharringtonine (HHT) has been widely used in China to treat patients with acute and chronic myeloid leukemia for decades. Since response to HHT varies among patients, our study aimed to identify biomarkers that might influence the response to HHT using a panel of various human lymphoblastoid cell lines (LCLs). Genome-wide association (GWA) analysis using single nucleotide polymorphism (SNP) and mRNA expression data was assessed for association with cytotoxicity to HHT in LCLs. Integrated analysis among SNPs, expression, AUC value was also performed to help select candidate genes for further functional characterization. Functional validation of candidate genes was performed using leukemia cell lines (U937, K562). Candidate genes were knocked down using specific siRNA and its response to HHT was assessed using MTS assay. We found that 15 expression probes were associated with HHT AUC with P < 10−4, and 96 individual probe sets with P < 10−3. Eighteen SNPs were associated with HHT AUC with P < 10−5 and 281 SNPs with P < 10−4. The integrated analysis identified 4 unique SNPs that were associated with both expression and AUC. Functional validation using siRNA knockdown in leukemia cell lines showed that knocking down CCDC88A, CTBP2, SOCS4 genes in U937 and K562 cells significantly altered HHT cytotoxicity. In summary, this study performed with LCLs can help to identify novel biomarker that might contribute to variation in response to HHT therapy.