Simultaneous Confidence Intervals to Compare Gene Expression Profiles Using ABC Transporter Taqman Microfluidic Cards
Total Page:16
File Type:pdf, Size:1020Kb
853-860.qxd 20/1/2010 12:25 ÌÌ Page 853 ONCOLOGY REPORTS 23: 853-860, 2010 853 Simultaneous confidence intervals to compare gene expression profiles using ABC transporter TaqMan microfluidic cards SARA PIZZAMIGLIO1*, GIACOMO COSSA2*, LAURA GATTI2, GIOVANNI L. BERETTA2, ELISABETTA CORNA2, STELLA TINELLI2, PAOLO VERDERIO1** and PAOLA PEREGO2** 1Medical Statistics and Biometry Unit, 2Preclinical Chemotherapy and Pharmacology Unit, Department of Experimental Oncology and Molecular Medicine, Fondazione IRCCS Istituto Nazionale Tumori, Via Venezian 1, 20133 Milan, Italy Received September 14, 2009; Accepted November 20, 2009 DOI: 10.3892/or_00000707 Abstract. The rapid evolution of techniques for measuring Introduction gene expression makes available substantial data which require careful analysis. In particular, relative quantification Ovarian carcinoma is a major cause of death from gynaeco- based on microfluidic cards allows performing of rapid large logical cancers in the Western world. This disease is usually scale analyses. In the present study, we employed ovarian diagnosed at advanced stage when surgery cannot be curative carcinoma cell lines resistant to cisplatin (IGROV-1/Pt1) or and, after debulking surgery, patients are treated with platinum- to a camptothecin (IGROV-1CPT/L), both characterized by a taxane combination chemotherapy, which has become a complex pattern of resistance to multiple agents, to examine standard treatment for advanced-stage disease (1). A large the expression of genes of the superfamily of ATP binding- number of cytotoxic agents including camptothecins, the cassette (ABC) transporters by TaqMan microfluidic cards clinically available DNA topoisomerase I inhibitors, are used with the aim of developing an analytical tool to process data in the treatment of recurrence, but most therapies have failed in this particular framework. The transcript quantification to improve survival (2). Drug resistance represents a major was based on the comparative threshold cycle method, which cause of treatment failure, but the cellular bases of resistance compares the expression of a target gene normalized to the to cytotoxic agents remain to be defined. Although extensive expression of one or more reference genes (relative quantifi- cellular studies have underlined that resistance is multi- cation). To process expression of ABC transporters, we applied factorial in nature, a major impact of defence mechanisms a statistical procedure based on multivariate approaches and i.e., factors regulating antitumor drug influx and efflux has re-sampling techniques. The transporters that were signi- been recognized (3,4). Indeed, for all cytotoxic agents, ficantly modulated included members of the ABCA, ABCB, inadequate intra-tumor concentration could explain at least ABCC and ABCG subgroups. A consistent up-regulation in part the ‘pharmacological’ resistance to treatment. of ABCC2 as compared with the parental IGROV-1 cell line Selected components of the super-family of ATP binding- was observed in the IGROV-1/Pt1 cells, whereas down- cassette (ABC) transporters have been implicated in conferring regulation of ABCC6 and ABCG1 was found in IGROV-1/ resistance to structurally unrelated antitumor drugs including CPT-L cells. The use of rigorous analytic tools for gene cisplatin and camptothecins (5,6). The recent achievement expression data in preclinical models may lead to the identi- that ABC transporters are a large family of genes supports fication of signatures to test in ovarian carcinoma clinical the need for novel studies directed at clarifying the relation- samples. Moreover, the developed procedure may be useful ship between less characterized transporters and resistance. in the analysis of relative quantification data obtained with Thus, the availability of techniques allowing the simultaneous microfluidic cards in different experimental settings. analysis of a large number of transcripts with minimal quantity of total RNA may be useful in an attempt to establish pro- cedures to translate to the clinical setting. In this context, it is _________________________________________ important to define rigorous analytical procedures to analyze consistently the gene expression data in an attempt to define Correspondence to: Dr Laura Gatti, Preclinical Chemotherapy tools to translate insights from preclinical models into the and Pharmacology Unit, Department of Experimental Oncology and clinical practise (7). Based on this rationale, we used ovarian Molecular Medicine, Fondazione IRCCS Istituto Nazionale Tumori, carcinoma cell systems including a parental cell line and two Via Venezian 1, 20133 Milan, Italy E-mail: [email protected] drug-resistant sublines exhibiting a variable pattern of response to structurally unrelated antitumor agents to *Contributed equally examine the expression of the whole ABC transporters super- **Senior co-authorship family. The amount of ABC transporter transcripts were quantified by means of TaqMan microfluidic cards, a tool Key words: simultaneous confidence intervals, microfluidic cards, that allows simultaneous real-time PCR analysis of a large ABC transporters, ovarian carcinoma number of mRNAs because TaqMan gene expression assays are carried out on a microfluidic card with minimal quantity 853-860.qxd 20/1/2010 12:25 ÌÌ Page 854 854 PIZZAMIGLIO et al: STATISTICAL ANALYSIS AND GENE EXPRESSION PROFILES of total RNA. The transcript quantification is based on and DNase digestion were carried out with the RNAqueous®- the comparative threshold cycle (Ct) method (8), which 4PCR Kit (Ambion Europe LTD, Huntingdon, UK), according compares the expression of a target gene normalized to the to the manufacturer's instructions. RNA purity and integrity expression of one or more reference genes (relative quanti- were assessed with denaturing gel electrophoresis and the fication). Although relative quantification based on micro- RNA was quantified spectrophotometrically and then stored fluidic cards enables researchers to perform rapidly large at -80˚C. cDNA synthesis was performed with the High- scale analyses, its statistical background needs to be better Capacity cDNA Reverse Transcription Kit (Applied Bio- investigated. Currently, the software programs supplied systems, Foster City, CA, USA) with a Master Mix containing along with the various qPCR instruments usually do not 2.5 U/μl of MultiScribe Reverse Transcriptase and 1 μg provide an adequate solution for the processing of the raw of total RNA. The reaction mixture was incubated at 25˚C data into meaningful results due to their intrinsic restrictions for 10 min, followed by 120 min at 37˚C and then by heat mainly related to a closed software architecture. Consequently, inactivation of the enzyme at 85˚C for 5 sec. We then mixed we applied an ad hoc statistical procedure based on multi- 2 μl of single-stranded cDNA (equivalent to around 100 ng variate approaches and re-sampling techniques specifically of total RNA) with 48 μl of nuclease-free water and 50 μl developed to process gene expression data obtained using of TaqMan Universal PCR Master Mix. After loading TaqMan microfluidic cards. This procedure, exploiting the 100 μl of the sample-specific PCR mixture into one sample comparative Ct method (8), is based on the following steps: port of the microfluidic cards (Human ABC Transporter data cleaning, outlier detection, data normalization and Panel; Applied Biosystems), the cards were centrifuged twice selection of a small subset of genes confidently inferred to for 1 min at 280 g and sealed to prevent well-to-well conta- be involved in the multidrug-resistant phenotype. mination. The cards were placed in the microfluidic card Sample Block of an ABI Prism 7900 HT Sequence Detection Materials and methods System (Applied Biosystems). The thermal cycling conditions were 2 min at 50˚C and 10 min at 95˚C, followed by 40 cycles Cell culture and drugs. The human ovarian carcinoma of 30 sec at 97˚C and 1 min at 59.7˚C. The assay for each IGROV-1 cell line and the cisplatin- and camptothecin- gene on the microfluidic card was carried out in triplicate, resistant sublines, IGROV-1/Pt1 and IGROV-1/CPT-L, were due to the design of this specific panel. The calculation of the maintained in RPMI-1640 medium (BioWhittaker Lonza, threshold cycle (Ct) values was performed using the SDS 2.2 Lonza Milano S.r.l., Italy) supplemented with 10% fetal software (Applied Biosystems), after automatically setting bovine serum (Gibco, Invitrogen, San Giuliano M, Italy). the baseline and the threshold. The IGROV-1/Pt1 variant was generated as described (9). IGROV-1/CPT-L cells were developed by exposure of the Methodological background. Microfluidic card experiments parental cell line to increasing concentrations of the lipo- allow the detection and quantification of a large number philic camptothecin gimatecan. Resistance was stable for of genes by taking advantage of the real-time PCR method. at least 6 months when cells were grown in the absence Usually, each j-th (j=1,2,...J) Ct replicate is processed of selecting agent. Cell cultures were routinely checked for according to the comparative Ct method to compute the being mycoplasma-free and experiments were carried out relative quantity (RQ) of the i-th (i=1,2,...I) transcript in the using cell lines at similar passages following thawing from a k-th (k=1,2,...K) sample, normalized to an endogenous gene frozen stock. Cisplatin (Teva Italia S.r.l., Milan, Italy) was (R) and relative to a calibrator sample