Identification of germline alterations of the mad homology 2 domain of SMAD3 and SMAD4 from the Ontario site of the breast cancer family registry (CFR)

Abstract Introduction A common feature of neoplastic cells is that mutations in SMADs can contribute to the loss of sensitivity to the anti-tumor effects of transforming growth factor-β (TGF-β). However, germline mutation analysis of SMAD3 and SMAD4, the principle substrates of the TGF-β signaling pathway, has not yet been conducted in breast cancer. Thus, it is currently unknown whether germline SMAD3 and SMAD4 mutations are involved in breast cancer predisposition. Methods We performed mutation analysis of the highly conserved mad-homology 2 (MH2) domains for both genes in genomic DNA from 408 non-BRCA1/BRCA2 breast cancer cases and 710 population controls recruited by the Ontario site of the breast cancer family registry (CFR) using denaturing high-performance liquid chromatography (DHPLC) and direct DNA sequencing. The results were interpreted in several ways. First, we adapted nucleotide diversity analysis to quantitatively assess whether the frequency of alterations differ between the two genes. Next, in silico tools were used to predict variants' effect on domain function and mRNA splicing. Finally, 37 cases or controls harboring alterations were tested for aberrant splicing using reverse-transcription polymerase chain reaction (PCR) and real-time PCR statistical comparison of germline expressions by non-parametric Mann-Whitney test of independent samples. Results We identified 27 variants including 2 novel SMAD4 coding variants c.1350G > A (p.Gln450Gln), and c.1701A > G (p.Ile525Val). There were no inactivating mutations even though c.1350G > A was predicted to affect exonic splicing enhancers. However, several additional findings were of note: 1) nucleotide diversity estimate for SMAD3 but not SMAD4 indicated that coding variants of the MH2 domain were more infrequent than expected; 2) in breast cancer cases SMAD3 was significantly over-expressed relative to controls (P A was associated with elevated germline expression (> 5-fold); 3) separate analysis using tissue expression data showed statistically significant over-expression of SMAD3 and SMAD4 in breast carcinomas. Conclusions This study shows that inactivating germline alterations in SMAD3 and SMAD4 are rare, suggesting a limited role in driving tumorigenesis. Nevertheless, aberrant germline expressions of SMAD3 and SMAD4 may be more common in breast cancer than previously suspected and offer novel insight into their roles in predisposition and/or progression of breast cancer

NCI60 Cancer Cell Line Panel Data and RNAi Analysis Help Identify EAF2 as a Modulator of Simvastatin and Lovastatin Response in HCT-116 Cells

Simvastatin and lovastatin are statins traditionally used for lowering serum cholesterol levels. However, there exists evidence indicating their potential chemotherapeutic characteristics in cancer. In this study, we used bioinformatic analysis of publicly available data in order to systematically identify the genes involved in resistance to cytotoxic effects of these two drugs in the NCI60 cell line panel. We used the pharmacological data available for all the NCI60 cell lines to classify simvastatin or lovastatin resistant and sensitive cell lines, respectively. Next, we performed whole-genome single marker case-control association tests for the lovastatin and simvastatin resistant and sensitive cells using their publicly available Affymetrix 125K SNP genomic data. The results were then evaluated using RNAi methodology. After correction of the p-values for multiple testing using False Discovery Rate, our results identified three genes (NRP1, COL13A1, MRPS31) and six genes (EAF2, ANK2, AKAP7, STEAP2, LPIN2, PARVB) associated with resistance to simvastatin and lovastatin, respectively. Functional validation using RNAi confirmed that silencing of EAF2 expression modulated the response of HCT-116 colon cancer cells to both statins. In summary, we have successfully utilized the publicly available data on the NCI60 cell lines to perform whole-genome association studies for simvastatin and lovastatin. Our results indicated genes involved in the cellular response to these statins and siRNA studies confirmed the role of the EAF2 in response to these drugs in HCT-116 colon cancer cells

A Whole-Genome SNP Association Study of NCI60 Cell Line Panel Indicates a Role of Ca2+ Signaling in Selenium Resistance

Epidemiological studies have suggested an association between selenium intake and protection from a variety of cancer. Considering this clinical importance of selenium, we aimed to identify the genes associated with resistance to selenium treatment. We have applied a previous methodology developed by our group, which is based on the genetic and pharmacological data publicly available for the NCI60 cancer cell line panel. In short, we have categorized the NCI60 cell lines as selenium resistant and sensitive based on their growth inhibition (GI50) data. Then, we have utilized the Affymetrix 125K SNP chip data available and carried out a genome-wide case-control association study for the selenium sensitive and resistant NCI60 cell lines. Our results showed statistically significant association of four SNPs in 5q33–34, 10q11.2, 10q22.3 and 14q13.1 with selenium resistance. These SNPs were located in introns of the genes encoding for a kinase-scaffolding protein (AKAP6), a membrane protein (SGCD), a channel protein (KCNMA1), and a protein kinase (PRKG1). The knock-down of KCNMA1 by siRNA showed increased sensitivity to selenium in both LNCaP and PC3 cell lines. Furthermore, SNP-SNP interaction (epistasis) analysis indicated the interactions of the SNPs in AKAP6 with SGCD as well as SNPs in AKAP6 with KCNMA1 with each other, assuming additive genetic model. These genes were also all involved in the Ca2+ signaling, which has a direct role in induction of apoptosis and induction of apoptosis in tumor cells is consistent with the chemopreventive action of selenium. Once our findings are further validated, this knowledge can be translated into clinics where individuals who can benefit from the chemopreventive characteristics of the selenium supplementation will be easily identified using a simple DNA analysis

Bioinformatic analyses identifies novel protein-coding pharmacogenomic markers associated with paclitaxel sensitivity in NCI60 cancer cell lines

<p>Abstract</p> <p>Background</p> <p>Paclitaxel is a microtubule-stabilizing drug that has been commonly used in treating cancer. Due to genetic heterogeneity within patient populations, therapeutic response rates often vary. Here we used the NCI60 panel to identify SNPs associated with paclitaxel sensitivity. Using the panel's GI50 response data available from Developmental Therapeutics Program, cell lines were categorized as either sensitive or resistant. PLINK software was used to perform a genome-wide association analysis of the cellular response to paclitaxel with the panel's SNP-genotype data on the Affymetrix 125 k SNP array. FastSNP software helped predict each SNP's potential impact on their gene product. mRNA expression differences between sensitive and resistant cell lines was examined using data from BioGPS. Using Haploview software, we investigated for haplotypes that were more strongly associated with the cellular response to paclitaxel. Ingenuity Pathway Analysis software helped us understand how our identified genes may alter the cellular response to paclitaxel.</p> <p>Results</p> <p>43 SNPs were found significantly associated (FDR < 0.005) with paclitaxel response, with 10 belonging to protein-coding genes (<it>CFTR</it>, <it>ROBO1</it>, <it>PTPRD</it>, <it>BTBD12</it>, <it>DCT</it>, <it>SNTG1</it>, <it>SGCD</it>, <it>LPHN2</it>, <it>GRIK1</it>, <it>ZNF607</it>). SNPs in <it>GRIK1</it>, <it>DCT</it>, <it>SGCD </it>and <it>CFTR </it>were predicted to be intronic enhancers, altering gene expression, while SNPs in <it>ZNF607 </it>and <it>BTBD12 </it>cause conservative missense mutations. mRNA expression analysis supported these findings as <it>GRIK1</it>, <it>DCT</it>, <it>SNTG1</it>, <it>SGCD </it>and <it>CFTR </it>showed significantly (p < 0.05) increased expression among sensitive cell lines. Haplotypes found in <it>GRIK1, SGCD, ROBO1, LPHN2</it>, and <it>PTPRD </it>were more strongly associated with response than their individual SNPs.</p> <p>Conclusions</p> <p>Our study has taken advantage of available genotypic data and its integration with drug response data obtained from the NCI60 panel. We identified 10 SNPs located within protein-coding genes that were not previously shown to be associated with paclitaxel response. As only five genes showed differential mRNA expression, the remainder would not have been detected solely based on expression data. The identified haplotypes highlight the role of utilizing SNP combinations within genomic loci of interest to improve the risk determination associated with drug response. These genetic variants represent promising biomarkers for predicting paclitaxel response and may play a significant role in the cellular response to paclitaxel.</p

Genetic Variations Associated with Resistance to Doxorubicin and Paclitaxel in Breast Cancer

Anthracycline- and taxane-based regimens have been the mainstay in treating breast cancer patients using chemotherapy. Yet, the genetic make-up of patients and their tumors may have a strong impact on tumor sensitivity to these agents and to treatment outcome. This study represents a new paradigm assimilating bioinformatic tools with in vitro model systems to discover novel genetic variations that may be associated with chemotherapy response in breast cancer. This innovative paradigm integrates drug response data for the NCI60 cell line panel with genome-wide Affymetrix SNP data in order to identify genetic variations associated with drug resistance. This genome wide association study has led to the discovery of 59 candidate loci that may play critical roles in breast tumor sensitivity to doxorubicin and paclitaxel. 16 of them were mapped within well-characterized genes (three related to doxorubicin and 13 to paclitaxel). Further in silico characterization and in vitro functional analysis validated their differential expression in resistant cancer cell lines treated with the drug of interest (over-expression of RORA and DSG1, and under-expression of FRMD6, SGCD, SNTG1, LPHN2 and DCT). Interestingly, three and six genes associated with doxorubicin and paclitaxel resistance, respectively, are involved in the apoptotic process in cells. A constructed interactome suggested that there is cross-talk at the Nrf-2 oxidative stress pathway between genes associated with resistance to doxorubicin and paclitaxel. This unique GWA approach serves as a proof-of-principle study and systematically investigates targets responsible for variable response to chemotherapy in breast tumor cells and possibly the tumors of breast cancer patients. Overall, the model discovered novel candidate genes that have not been previously associated with doxorubicin and paclitaxel cytotoxicity. Future studies will be directed at illustrating a causative relationship between the observed genomic changes and drug resistance in breast cancer patients undergoing doxorubicin and paclitaxel chemotherapy.Ph

KCNMA1 knock down in LNCaP and PC3 cell lines using siRNA methodology.

<p>The percent reduction in response is plotted as a bar graft, the first bar in each set indicating the absence and the second bar representing the presence of selenium treatment of both cell lines. The siRNA experiments are done in triplicates and the average is taken. The (*) indicates the statistically significant (p<0.05) differences in response to selenium treatment.</p

Summary of the whole-genome case-control association study for the Selenium resistant and sensitive NCI60 cell line panel.

<p>C<b>hr</b>: chromosome; <b>All</b>: alleles; <b>CHISQ</b>: X²; <b>OR</b>: odds ratio; <b>L95</b>: lower confidence interval; <b>U95</b>: upper confidence interval; <b>FDR_BH</b>: False discovery rate by Benjamini and Hochberg; <b>BONF:</b> Bonferroni method.</p