Name: using rna-sequencing to detect novel splice variants related to drug resistance in in vitro cancer models
Uploaded: 2016-12-09T15:00:00
Description: Watch this Scientific Journal Video about Using RNA-sequencing to Detect Novel Splice Variants Related to Drug Resistance in In Vitro Cancer Models at JoVE.com

The authors would like to acknowledge prof. J.J. McGuire, prof. R. Kofler and dr. K. Quint for providing the resistant cell lines used in this work. The study has been founded by grants from Cancer Center Amsterdam (CCA) Foundation (to JC, EG and RS), the KiKa (Children Cancer-free grant for AW) foundation, the Law Offices of Peter G. Angelos Grant from the Mesothelioma Applied Research Foundation (to VEG and EG), Associazione Italiana per la Ricerca sul Cancro (AIRC), Istituto Toscano Tumori (ITT), and Regione Toscana Bando FAS Salute (to EG).

1. Characterization of Drug Resistance Profiles through Cytotoxicity Assays
<ol>
	<li>Leukemic Cell Line Culture
	<ol>
		<li>Maintain the parental T-cell ALL cell line CCRF-CEM (CEM-WT) as well as its drug resistant sublines, including CEM/R30dm, CEM-R5 and CEM-C3, in 25 cm2 in 10 ml RPMI-1640 medium containing 2.3 &#181;M folic acid supplemented with 10% fetal calf serum and 100 units/ml penicillin G and 100 &#181;g/ml streptomycin.</li>
		<li>Culture the cells in a humidified atmosphere at...

<ol>
	<li>Gottesman, M. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mechanisms+of+cancer+drug+resistance.">Mechanisms of cancer drug resistance.</a> Annu Rev Med. 53, 615-627 (2002).</li><li>McDermott, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In+vitro+Development+of+Chemotherapy+and+Targeted+Therapy+Drug-Resistant+Cancer+Cell+Lines:+A+Practical+Guide+with+Case+Studies.">In vitro Development of Chemotherapy and Targeted Therapy Drug-Resistant Cancer Cell Lines: A Practical Guide with Case Studies.</a> Front Oncol. 4, 40 (2014).</li><li>Kaspers, G. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Prednisolone+resistance+in+childhood+acute+lymphoblastic+leukemia:+vitro-vivo+correlations+and+cross-resistance+to+other+drugs.">Prednisolone resistance in childhood acute lymphoblastic leukemia: vitro-vivo correlations and cross-resistance to other drugs.</a> Blood. (1), 259-266 (1998).</li><li>van Meerloo, J., Kaspers, G. J., Cloos, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cell+sensitivity+assays:+the+MTT+assay.">Cell sensitivity assays: the MTT assay.</a> Methods Mol Biol. 731, 237-245 (2011).</li><li>Vichai, V., Kirtikara, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Sulforhodamine+B+colorimetric+assay+for+cytotoxicity+screening.">Sulforhodamine B colorimetric assay for cytotoxicity screening.</a> Nat Protoc. 1 (3), 1112-1116 (2006).</li><li>Wojtuszkiewicz, A., Assaraf, Y. G., Maas, M. J., Kaspers, G. J., Jansen, G., Cloos, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pre-mRNA+splicing+in+cancer:+the+relevance+in+oncogenesis,+treatment+and+drug+resistance.">Pre-mRNA splicing in cancer: the relevance in oncogenesis, treatment and drug resistance.</a> Expert Opin Drug Metab Toxicol. 5, 673-689 (2015).</li><li>Liu, S., Cheng, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Alternative+RNA+splicing+and+cancer.">Alternative RNA splicing and cancer.</a> Wiley Interdiscip Rev RNA. 4 (5), 547-566 (2013).</li><li>Wojtuszkiewicz, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Folylpolyglutamate+synthetase+splicing+alterations+in+acute+lymphoblastic+leukemia+are+provoked+by+methotrexate+and+other+chemotherapeutics+and+mediate+chemoresistance.">Folylpolyglutamate synthetase splicing alterations in acute lymphoblastic leukemia are provoked by methotrexate and other chemotherapeutics and mediate chemoresistance.</a> Int J Cancer. 138 (7), 1645-1656 (2016).</li><li>Shendure, J., Ji, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Next-generation+DNA+sequencing.">Next-generation DNA sequencing.</a> Nat Biotechnol. 10, 1135-1145 (2008).</li><li>Wang, Z., Gerstein, M., Snyder, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RNA-Seq:+a+revolutionary+tool+for+transcriptomics.">RNA-Seq: a revolutionary tool for transcriptomics.</a> Nat Rev Genet. 1, 57-63 (2009).</li><li>Ozsolak, F., Milos, P. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RNA+sequencing:+advances,+challenges+and+opportunities.">RNA sequencing: advances, challenges and opportunities.</a> Nat Rev Genet. 2, 87-98 (2011).</li><li>Han, Y., Gao, S., Muegge, K., Zhang, W., Zhou, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Advanced+Applications+of+RNA+Sequencing+and+Challenges.">Advanced Applications of RNA Sequencing and Challenges.</a> Bioinform Biol Insights. 9 (Suppl 1), 29-46 (2015).</li><li>Khatoon, Z., Figler, B., Zhang, H., Cheng, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Introduction+to+RNA-Seq+and+its+applications+to+drug+discovery+and+development.">Introduction to RNA-Seq and its applications to drug discovery and development.</a> Drug Dev Res. 5, 324-330 (2014).</li><li>Zhao, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+stranded+and+non-stranded+RNAseq+transcriptome+profiling+and+investigation+of+gene+overlap.">Comparison of stranded and non-stranded RNAseq transcriptome profiling and investigation of gene overlap.</a> BMC Genomics. 16, 675 (2015).</li><li>Carrara, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Alternative+splicing+detection+workflow+needs+a+careful+combination+of+sample+prep+and+bioinformatics+analysis.">Alternative splicing detection workflow needs a careful combination of sample prep and bioinformatics analysis.</a> BMC Bioinformatics. 16 (Suppl 9), (2015).</li><li>Dobin, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=STAR:+ultrafast+universal+RNA-seq+aligner.">STAR: ultrafast universal RNA-seq aligner.</a> Bioinformatics. 29, 15-21 (2013).</li><li>Wang, L., Wang, S., Li, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RSeQC:+quality+control+of+RNA-seq+experiments.">RSeQC: quality control of RNA-seq experiments.</a> Bioinformatics. 28, 2184-2185 (2012).</li><li>Hooper, J. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+survey+of+software+for+genome-wide+discovery+of+differential+splicing+in+RNA-Seq+data.">A survey of software for genome-wide discovery of differential splicing in RNA-Seq data.</a> Hum Genomics. 8, (2014).</li><li>Shen, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MATS:+a+Bayesian+framework+for+flexible+detection+of+differential+alternative+splicing+from+RNA-Seq+data.">MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data.</a> Nucleic Acids Res. 8, e61 (2012).</li><li>Vargas, I. M., Vivas-Mej&#236;a, P. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Assessment+of+mRNA+splice+variants+by+qRT-PCR.">Assessment of mRNA splice variants by qRT-PCR.</a> Methods Mol Biol. 1049, 171-186 (2013).</li><li>Hala, M., Hartmann, B. L., B&#246;ck, G., Geley, S., Kofler, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Glucocorticoid-receptor-gene+defects+and+resistance+to+glucocorticoid-induced+apoptosis+in+human+leukemic+cell+lines.">Glucocorticoid-receptor-gene defects and resistance to glucocorticoid-induced apoptosis in human leukemic cell lines.</a> Int J Cancer. 68 (5), 663-668 (1996).</li><li>Schmidt, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Glucocorticoid+resistance+in+two+key+models+of+acute+lymphoblastic+leukemia+occurs+at+the+level+of+theglucocorticoid+receptor.">Glucocorticoid resistance in two key models of acute lymphoblastic leukemia occurs at the level of theglucocorticoid receptor.</a> FASEB J. 20 (14), 2600-2602 (2006).</li><li>McCloskey, D. E., McGuire, J. J., Russell, C. A., Rowan, B. G., Bertino, J. R., Giuseppe Pizzorno, G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Decreased+folylpolyglutamate+synthetase+activity+as+a+mechanism+of+methotrexate+resistance+in+CCRF-CEM+human+leukemia+sublines.">Decreased folylpolyglutamate synthetase activity as a mechanism of methotrexate resistance in CCRF-CEM human leukemia sublines.</a> The Journal of Biological Chemistry. 266 (10), 6181-6187 (1991).</li><li>Quint, K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pancreatic+cancer+cells+surviving+gemcitabine+treatment+express+markers+of+stem+cell+differentiation+and+epithelial-mesenchymal+transition.">Pancreatic cancer cells surviving gemcitabine treatment express markers of stem cell differentiation and epithelial-mesenchymal transition.</a> Int J Oncol. 41 (6), 2093-2102 (2012).</li><li>Lin, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DDX5+is+a+positive+regulator+of+oncogenic+NOTCH1+signaling+in+T+cell+acute+lymphoblastic+leukemia.">DDX5 is a positive regulator of oncogenic NOTCH1 signaling in T cell acute lymphoblastic leukemia.</a> Oncogene. 32 (40), 4845-4853 (2013).</li><li>Mazurek, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Acquired+dependence+of+acute+myeloid+leukemia+on+the+DEAD-box+RNA+helicase+DDX5.">Acquired dependence of acute myeloid leukemia on the DEAD-box RNA helicase DDX5.</a> Cell Rep. 7 (6), 1887-1899 (2014).</li><li>Calabretta, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Modulation+of+PKM+alternative+splicing+by+PTBP1+promotes+gemcitabine+resistance+in+pancreatic+cancer+cells.">Modulation of PKM alternative splicing by PTBP1 promotes gemcitabine resistance in pancreatic cancer cells.</a> Oncogene. , (2016).</li><li>Azoitei, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=PKM2+promotes+tumor+angiogenesis+by+regulating+HIF-1&#945;+through+NF-&#954;B+activation.">PKM2 promotes tumor angiogenesis by regulating HIF-1&#945; through NF-&#954;B activation.</a> Mol Cancer. 15 (1), (2016).</li><li>Samuels, A. L., Heng, J. Y., Beesley, A. H., Kees, U. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Bioenergetic+modulation+overcomes+glucocorticoid+resistance+in+T-lineage+acute+lymphoblastic+leukaemia.">Bioenergetic modulation overcomes glucocorticoid resistance in T-lineage acute lymphoblastic leukaemia.</a> Br J Haematol. 165 (1), 57-66 (2014).</li><li>Rots, M. G., Pieters, R., Kaspers, G. J., Veerman, A. J., Peters, G. J., Jansen, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Classification+of+ex+vivo+methotrexate+resistance+in+acute+lymphoblastic+and+myeloid+leukaemia.">Classification of ex vivo methotrexate resistance in acute lymphoblastic and myeloid leukaemia.</a> Br J Haematol. 110 (4), 791-800 (2000).</li><li>Funel, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Laser+microdissection+and+primary+cell+cultures+improve+pharmacogenetic+analysis+in+pancreatic+adenocarcinoma.">Laser microdissection and primary cell cultures improve pharmacogenetic analysis in pancreatic adenocarcinoma.</a> Lab Invest. 88 (7), 773-784 (2008).</li><li>Gillet, J. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Redefining+the+relevance+of+established+cancer+cell+lines+to+the+study+of+mechanisms+of+clinical+anti-cancer+drug+resistance.">Redefining the relevance of established cancer cell lines to the study of mechanisms of clinical anti-cancer drug resistance.</a> Proc Natl Acad Sci U S A. 108 (46), 18708-18713 (2011).</li><li>Kole, R., Krainer, A. R., Altman, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RNA+therapeutics:+beyond+RNA+interference+and+antisense+oligonucleotides.">RNA therapeutics: beyond RNA interference and antisense oligonucleotides.</a> Nat Rev Drug Discov. 11 (2), 125-140 (2012).</li><li>Shen, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=rMATS:+robust+and+flexible+detection+of+differential+alternative+splicing+from+replicate+RNA-Seq+data.">rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data.</a> Proc Natl Acad Sci U S A. 111 (51), E5593-E5601 (2014).</li></ol>

Here we describe a novel approach that combines well-established cytotoxicity screening techniques and powerful NGS-based transcriptomic analyses to identify differential splicing events in relation to drug resistance. Spectrophotometric assays are convenient and robust high-throughput methods to assess drug sensitivity in in vitro cancer models and represent the first choice for many laboratories performing cytotoxicity screenings. Troubleshooting as well as possible variations for this method were extensively ...

Despite considerable advances in cancer treatment, resistance of malignant cells to chemotherapy, either intrinsic or acquired upon prolonged drug exposure, is the major reason for treatment failure in a wide range of leukemia and solid tumors1.
In order to delineate the mechanisms underlying drug resistance, in vitro cell line models are developed by stepwise selection of cancer cells resistant to chemotherapeutic agents. This procedure mimics the regimes used in the clinical settings and therefore allows in depth investigation of relevant resistance mechanisms. Resistant cells which survive the treatment are then distinguished from parental sensitive cells by using cell viability/cytotoxicity assays2. In vitro drug resistance profiles of primary cells have been shown to be significantly related to clinical response to chemotherapy3.
High-throughput cytotoxicity assays constitute a convenient method to determine drug sensitivity in vitro. Herein, the viability of cells is assessed by for instance the 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyl tetrazolium bromide - MTT assay4, which is based on metabolic conversion of certain substrates (i.e., tetrazolium salts) into colored products, thereby reflecting the mitochondrial activity of cells. Alternatively, the cellular protein content can be quantified using the sulforhodamine B (SRB) assay5. Here, the number of viable cells is proportional to the optical density (OD) measured at an appropriate wavelength in a spectrophotometer, with no need of extensive and time-consuming cell counting procedures. The growth inhibition induced by a certain chemotherapeutic drug can be calculated based on the OD of the wells in which cells were treated with a test agent and compared with the OD of untreated control cells. A dose-response curve is obtained by plotting drug concentrations versus percentages of viable cells relative to control cells. Finally, drug sensitivity can be reported as the concentration that results in 50% of cell growth inhibition as compared to untreated cells (IC50).
The mechanisms underlying drug resistance include many different abnormalities, such as alterations affecting gene expression of determinants of drug activity and cellular metabolism. These molecular lesions, including mutations, aberrations at a transcriptional and post-transcriptional level as well as disturbed epigenetic regulation often affect genes involved either in drug metabolism or apoptosis6.
Alternative pre-mRNA splicing and its intricate regulation have recently received considerable attention as a novel entity that may dictate drug resistance of cancer cells7. Up to 95% of human genes are alternatively spliced in normal cells by means of this tightly regulated process which produces many different protein isoforms from the same gene. Alternative splicing is often deregulated in cancer and several tumors are characterized by altered splicing of a growing number of genes involved in drug metabolism (i.e., deoxycytidine kinase, folylpolyglutamate synthetase, or Multidrug resistance proteins)6,8. However, comprehensive analysis of splicing profiles of drug resistant cells is painfully lacking. Therefore, it is imperative to develop high throughput methods for alternative splicing analysis. This could help to develop more effective therapeutic approaches.
During the last decade, the rapid development of next generation sequencing (NGS) technologies has enriched biomedical research with new insights into the molecular mechanisms governing regulation of genome expression and their role in various biological processes9. RNA-sequencing (RNA-seq) is a powerful sub-application of NGS in the field of transcriptomics. It allows a genome-wide profiling (both qualitatively and quantitatively) of the expression patterns of thousands of genes simultaneously and is well suited for the characterization of novel coding mRNAs as well as long non-coding RNA, miRNA, siRNA, and other small RNA classes (e.g., snRNA and piRNA)10,11.
RNA-Seq has many advantages over previous technologies for transcriptome characterization (e.g., Sanger sequencing and expression microarrays). It is not based on existing genome annotation, it has a single-nucleotide level of resolution and it has a broader dynamic range for expression level estimation. Briefly, the basic experimental workflow of RNA-seq experiments consists of polyadenylated transcript (mRNA) selection and fragmentation, followed by conversion into cDNA, library construction and finally, massively parallel deep sequencing12,13. Due to rapid drop of sequencing costs over the last few years, RNA-seq is gradually replacing other technologies and significant efforts are being made to improve the library preparation protocols. For instance, it is now possible to retain the strand information of mRNA transcripts by marking the second strand cDNA with deoxyuridine triphosphate (dUTP) and, prior to PCR amplification, digesting the marked strand with uracil-DNA-glycosilase (UDG). This process enhances the accuracy of gene annotation and estimation of the expression levels14,15.
The analysis and interpretation of RNA-seq data require complex and powerful computational software packages and processing within bioinformatic pipelines16,17. First, the raw reads undergo quality control by removing technical and biological artefacts and discarding (trimming) the sequences which do not reach stringent quality requirements. Subsequently the reads for each sample are mapped to a reference genome and indexed into gene-level, exon-level, or transcript-level, in order to determine the abundance of each category. Depending on the application, refined data are then computed through statistical models for the identification of allele-specific expression, alternative splicing, gene fusions and single nucleotide polymorphisms (SNPs)12. Finally, differential analysis on selected level (i.e., gene expression or alternative splicing) can be used to compare samples obtained under different conditions.
Differential splicing analysis describes the differences in splice site usage between two samples. An increasing number of software packages devoted to this purpose are available based on different statistical models, performances and user interface18. Among these, MATS (Multivariate Analysis of Transcript Splicing) emerges as a freely available and precise computational tool based on a Bayesian statistical framework and designed to detect differential splicing events from either single or paired end RNA-seq data. Starting from the aligned (.bam) files, MATS can detect all major types of alternative splicing events (exon skipping, alternative 3&#39; splice site, alternative 5&#39; splice site, mutually exclusive exons and intron retention - also see Figure 1).
First, the software identifies reads which support a certain splice event, for instance exon skipping, and classifies them into two types. &#34;Inclusion reads&#34; (for the canonical splice event) map within the investigated exon and span the junctions between that specific exon and the two upstream and downstream flanking exons. &#34;Skipping reads&#34; (for the alternative splice event) span the junction between the two flanking exons. Subsequently, MATS returns the normalized inclusion level for both the canonical and alternative events and compares values between samples or conditions. Ultimately, it calculates P-value and false discovery rate (FDR) assuming that the difference in the variant ratio of a gene between two conditions exceeds a given user-defined threshold for each splicing event19,34.
Following differential splicing analysis in conjunction with RNA-seq, an extensive experimental validation is warranted in order to identify true-positive gene candidates18. Quantitative reverse transcribed-polymerase chain reaction (qRT-PCR) is the most commonly used and optimal method in validation of candidates obtained from RNA-Seq analysis20. The aim of this paper is to provide a robust methodology to investigate drug resistance-related splicing profiles in solid tumors and hematological malignancies. Our approach utilizes RNA-seq-based transcriptome profiling of selected cell line models of drug resistant cancers in combination with an established qRT-PCR method for the validation of candidate genes implicated in drug resistance.
The human leukemia cell line models used in this study included pediatric T-cell acute lymphoblastic leukemia (T-ALL) cell line CCRF-CEM (CEM-WT), its two glucocorticoid (GC)-resistant subclones CEM-C7H2-R5C3 (CEM-C3) and CEM-C7R5 (CEM-R5)21,22 and the methotrexate (MTX)-resistant subline CEM/R30dm23. Although current therapies based on GCs and MTX establish clinical benefit in about 90% of cases, the emergence of GC-resistance still represents an unsolved problem with an unclear molecular mechanism. To isolate GC-resistant sub-clones, CEM-WT cells were cultured in 1 &#181;M dexamethasone (Dex) for 2 to 3 weeks. MTX-resistant subline CEM/R30dm was developed through repeated short- term (24 hr) exposure of CEM-WT cells to 30 &#181;M MTX as a mimic of clinical protocols. Interestingly, this cell line also displayed cross-resistance to Dex (unpublished results) for which the mechanism is not fully understood.
The solid tumor model investigated in the present study is pancreatic ductal adenocarcinoma, notorious for its extraordinary refractoriness to chemotherapy. To this end, we selected Panc-1 cell line and its gemcitabine-resistant sub-clone Panc-1R obtained by continuous incubation with 1 &#181;M of the drug24. Here we describe an approach to discover novel mechanisms underlying in-vitro drug resistance by combining three protocols: colorimetric cytotoxicity assays to assess drug sensitivity in leukemic cells and cancer cells from solid tumors, RNA-seq-based pipeline to identify novel splice variants related to drug sensitivity/resistance and RT-PCR and qRT-PCR analysis to validate potential candidates.

<table border="0" cellpadding="0" cellspacing="0" width="732">
	<tbody>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">Sulforhodamine B</td>
			<td style="width:195px;">Sigma-Aldrich</td>
			<td style="width:119px;">230162</td>
			<td style="width:167px;"></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">Trichloroacetic acid</td>
			<td style="width:195px;">Sigma-Aldrich</td>
			<td style="width:119px;">251399</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">CCRF-CEM</td>
			<td style="width:195px;">ATCC, Manassas, VA, USA</td>
			<td style="width:119px;">ATCC CCL-119</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">Panc-1</td>
			<td style="width:195px;">ATCC, Manassas, VA, USA</td>
			<td style="width:119px;">ATCC CRL-1469</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">DMEM high glucose</td>
			<td style="width:195px;">Lonza, Basel, Switzerland</td>
			<td style="width:119px;">12-604F</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">RPMI-1640&#160;</td>
			<td style="width:195px;">Gibco, Carlsbad, CA, USA</td>
			<td style="width:119px;">11875093</td>
			<td></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:252px;">Fetal bovine and calf serum</td>
			<td style="width:195px;">Greiner Bio-One, Frickenhausen, Germany</td>
			<td style="width:119px;">758093</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">penicillin G streptomycin sulphate</td>
			<td style="width:195px;">Gibco, Carlsbad, CA, USA</td>
			<td style="width:119px;">15140122</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">Tris(hydroxymethyl)-aminomethane</td>
			<td style="width:195px;">Sigma Aldrich</td>
			<td style="width:119px;">252859</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">MTT formazan</td>
			<td style="width:195px;">Sigma Aldrich</td>
			<td style="width:119px;">M2003</td>
			<td></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:252px;">Anthos-Elisa-reader 2001</td>
			<td style="width:195px;">Labtec, Heerhugowaard, Netherlands</td>
			<td></td>
			<td>UV-Vis 96-well plate spectrophotometer</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">Greiner CELLSTAR 96 well plates</td>
			<td style="width:195px;">Greiner/Sigma</td>
			<td>M0812-100EA</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:252px;">Trypsin/EDTA Solution 100 ml</td>
			<td style="width:195px;">Lonza, Basel, Switzerland</td>
			<td>CC-5012</td>
			<td></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;">CELLSTAR Cell Culture Flasks 25 cm2</td>
			<td style="width:195px;">Greiner Bio-One, Frickenhausen, Germany</td>
			<td>82051-074</td>
			<td></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;">CELLSTAR Cell Culture Flasks 75 cm3</td>
			<td style="width:195px;">Greiner Bio-One, Frickenhausen, Germany</td>
			<td>82050-856</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Phosphate Buffered Saline (NaCl 0.9%)</td>
			<td>B.Braun Melsungen AG, Germany</td>
			<td>362 3140</td>
			<td></td>
		</tr>
	</tbody>
</table>

using rna-sequencing to detect novel splice variants related to drug resistance in in vitro cancer models

Drug resistance remains a major problem in the treatment of cancer for both hematological malignancies and solid tumors. Intrinsic or acquired resistance can be caused by a range of mechanisms, including increased drug elimination, decreased drug uptake, drug inactivation and alterations of drug targets. Recent data showed that other than by well-known genetic (mutation, amplification) and epigenetic (DNA hypermethylation, histone post-translational modification) modifications, drug resistance mechanisms might also be regulated by splicing aberrations. This is a rapidly growing field of investigation that deserves future attention in order to plan more effective therapeutic approaches. The protocol described in this paper is aimed at investigating the impact of aberrant splicing on drug resistance in solid tumors and hematological malignancies. To this goal, we analyzed the transcriptomic profiles of several in vitro models through RNA-seq and established a qRT-PCR based method to validate candidate genes. In particular, we evaluated the differential splicing of DDX5 and PKM transcripts. The aberrant splicing detected by the computational tool MATS was validated in leukemic cells, showing that different DDX5 splice variants are expressed in the parental vs. resistant cells. In these cells, we also observed a higher PKM2/PKM1 ratio, which was not detected in the Panc-1 gemcitabine-resistant counterpart compared to parental Panc-1 cells, suggesting a different mechanism of drug-resistance induced by gemcitabine exposure.

The cytotoxicity assays described in the protocol provide a reliable and robust method to assess the resistance of cancer cells to chemotherapeutic agents in vitro. By means of the MTT assay, sensitivity to Dex was determined in four T-ALL cell lines, including Dex-sensitive parental CEM-WT cells, and three Dex-resistant sublines: CEM/R30dm, CEM-R5 and CEM-C3. Two different concentration ranges had to be used due to the large difference in sensitivity between CEM-WT (2 &#181;M - ...

Watch this Scientific Journal Video about Using RNA-sequencing to Detect Novel Splice Variants Related to Drug Resistance in In Vitro Cancer Models at JoVE.com

Using RNA-sequencing to Detect Novel Splice Variants Related to Drug Resistance in In Vitro Cancer Models

Here we describe a protocol aimed at investigating the impact of aberrant splicing on drug resistance in solid tumors and hematological malignancies. To this goal, we analyzed the transcriptomic profiles of parental and resistant in vitro models through RNA-seq and established a qRT-PCR based method to validate candidate genes.

Using RNA-sequencing to Detect Novel Splice Variants Related to Drug Resistance in In Vitro Cancer Models

Drug resistance remains a major problem in the treatment of cancer for both hematological malignancies and solid tumors. Intrinsic or acquired resistance ...

Research

JoVE Journal

Cancer Research

Utilizing Functional Genomics Screening to Identify Potentially Novel Drug Targets in Cancer Cell Spheroid Cultures

The identification of functional driver events in cancer is central to furthering our understanding of cancer biology and indispensable for the discovery ...

Establishing Dual Resistance to EGFR-TKI and MET-TKI in Lung Adenocarcinoma Cells In Vitro with a 2-step Dose-escalation Procedure

Establishing Dual Resistance to EGFR-TKI and MET-TKI in Lung Adenocarcinoma Cells In Vitro with a 2-step Dose-escalation Procedure

Drug resistance is a major challenge in cancer therapy. The generation of resistant sublines in vitro is necessary for discovering novel mechanisms to ...

Generation of Prostate Cancer Cell Models of Resistance to the Anti-mitotic Agent Docetaxel

Microtubule targeting agents (MTAs) are a mainstay in the treatment of a wide range of tumors. However, acquired resistance to chemotherapeutic drugs is a ...

Evaluating In Vitro DNA Damage Using Comet Assay

Evaluating In Vitro DNA Damage Using Comet Assay

DNA damage is a common phenomenon for each cell during its lifespan, and is defined as an alteration of the chemical structure of genomic DNA. Cancer ...

Establishment of Gastric Cancer Patient-derived Xenograft Models and Primary Cell Lines

The use of preclinical models to advance our understanding of tumor biology and investigate the efficacy of therapeutic agents is key to cancer research. ...

In Vitro Assay to Study Tumor-macrophage Interaction

Tumor associated macrophages (TAMs) account for a large percentage of cells in the tumor mass for different types of cancers. Glioblastoma (GBM), a ...

In Vitro Evaluation of Oncogenic Transformation in Human Mammary Epithelial Cells

Tumorigenesis is a multi-step process in which cells acquire capabilities&#160;that allow their growth, survival, and dissemination under hostile ...

Patient-Derived Tumor Explants As a "Live" Preclinical Platform for Predicting Drug Resistance in Patients

An understanding of drug resistance and the development of novel strategies to sensitize highly resistant cancers rely on the availability of suitable ...

Pooled shRNA Library Screening to Identify Factors that Modulate a Drug Resistance Phenotype

Understanding clinically relevant driver mechanisms of acquired chemo-resistance is crucial for elucidating ways to circumvent resistance and improve ...

Quantifying Telomere Length Using Non-Radioactive Chemiluminescence Detection

Optimization of Performance Parameters of the TAGGG Telomere Length Assay

Optimization of Performance Parameters of the TAGGG Telomere Length Assay (Video) | JoVE 

Telomeres are repetitive sequences which are present at chromosomal ends; their shortening is a characteristic feature of human somatic cells. Shortening ...