Transposon activation mutagenesis as a screening tool for identifying resistance to cancer therapeutics
© Chen et al; licensee BioMed Central Ltd. 2013
Received: 1 November 2012
Accepted: 19 February 2013
Published: 27 February 2013
The development of resistance to chemotherapies represents a significant barrier to successful cancer treatment. Resistance mechanisms are complex, can involve diverse and often unexpected cellular processes, and can vary with both the underlying genetic lesion and the origin or type of tumor. For these reasons developing experimental strategies that could be used to understand, identify and predict mechanisms of resistance in different malignant cells would be a major advance.
Here we describe a gain-of-function forward genetic approach for identifying mechanisms of resistance. This approach uses a modified piggyBac transposon to generate libraries of mutagenized cells, each containing transposon insertions that randomly activate nearby gene expression. Genes of interest are identified using next-gen high-throughput sequencing and barcode multiplexing is used to reduce experimental cost.
Using this approach we successfully identify genes involved in paclitaxel resistance in a variety of cancer cell lines, including the multidrug transporter ABCB1, a previously identified major paclitaxel resistance gene. Analysis of co-occurring transposons integration sites in single cell clone allows for the identification of genes that might act cooperatively to produce drug resistance a level of information not accessible using RNAi or ORF expression screening approaches.
We have developed a powerful pipeline to systematically discover drug resistance in mammalian cells in vitro. This cost-effective approach can be readily applied to different cell lines, to identify canonical or context specific resistance mechanisms. Its ability to probe complex genetic context and non-coding genomic elements as well as cooperative resistance events makes it a good complement to RNAi or ORF expression based screens.
KeywordsTransposon mutagenesis Chemotherapy Resistance Gene activation
The development of resistance to cancer therapeutics represents a major hindrance to the successful pharmacological treatment and eradication of tumors in patients. Although some progress has been made in combining or augmenting treatments to counteract resistance, a major obstacle is our limited understanding of the mechanisms of resistance to current or novel therapeutics. Drug resistance can be mediated by further genetic and/or epigenetic changes in the tumor and, with the advent of high throughput sequencing, it is now feasible to systematically survey mutations in tumor genomes from patients following resistance development. However, the identification of the relevant ‘driver’ mutations, and other potential targets in resistance pathways, remains challenging.
A complementary approach is to identify resistance pathways experimentally using in vitro culture or animal model systems. Findings from such studies can then be used to inform analysis of patient samples and develop therapies to counteract resistance. Direct experimental identification of resistance genes has focused largely on reverse genetic and chemical biology approaches, including cDNA and RNAi library screens [1, 2] or combined small molecule inhibitor and siRNA screens . Such approaches can require expensive reagents and specialized platforms, and the need to consistently deliver siRNAs limits their applicability. Perhaps more importantly, as reverse genetic approaches, they are biased toward previously characterized genetic elements.
Forward genetic approaches using mobile genetic elements provide a powerful alternative method for gene discovery that can overcome many of the limitations of reverse genetic approaches. Mutagenesis with mobile genetic elements that insert into the genome offers a great scope for screening as these provide readily detected tags to identify insertion sites, and can potentially either activate or disrupt gene expression. Retroviruses have been used for insertional mutagenesis to identify oncogenes and study therapeutic resistance in tumors [4–6], however they preferentially insert in regions of open chromatin and high gene expression, leading to potential bias in results from genome-wide screens. Furthermore, the requirements for viral long terminal repeats (LTRs) and other structural restrictions limit the use of complex DNA constructs, limiting its applications to loss-of-function mutagenesis  and specialized haploid cell lines .
Transposons, another class of mobile genetic elements , have increasingly been utilized as genetic tools in mammals after the discovery and engineering of two transposons, Sleeping Beauty (SB) and piggyBac (PB) [10–13]. A major advantage of transposons is the simplicity of their integration machinery, which permits the incorporation of long DNA sequences, including functional genetic elements such as promoters, transcriptional stops and splicing sequences. This flexibility has allowed development of a variety of powerful mutagenesis schemes [14, 15]. In their simplest application, transposons disrupt genes leading to loss of function, logically analogous to RNAi screens. With the incorporation of splice acceptors and reporter genes, transposons can also be used as an alternative to retroviral gene-traps [16, 17]. Such gene disruption approaches are the basis for genome-wide insertion libraries in mouse embryonic stem cells [14, 18]. Alternatively, inclusion of functional promoters within the transposon creates “activation tags” that cause expression of genes in which they land . Activation tagging has been used in mouse somatic models to identify oncogenes [20, 21]. This approach has great potential for gene discovery as it combines the strong phenotype of ‘gain-of-function’ approaches with the ability to probe the entire genome, including novel or uncharacterized genes and transcripts.
Here we report the development of transposon-based gene activation tagging for discovery of chemotherapeutic resistance genes. We constructed an activation PB transposon, generated mutagenesis libraries from several cancer cell lines, and characterized the mutations by sample barcoding and high-throughput sequencing. We validated this system by screening for genes involved in resistance to the microtubule targeting drug paclitaxel and identifying the multidrug resistance (MDR) gene ABCB1 as the primary gene target. Through further analysis of individual paclitaxel resistant clones, we also identify potential modifiers of ABCB1-mediated resistance. Hence, this study establishes a robust, flexible and adaptable system for identifying drug resistance.
Transposon plasmid PB-SB-PGK-neo-bpA and transposase plasmid pCMV-PBase were obtained from Pentao Liu of the Wellcome Trust Sanger Institute. This plasmid was designed as an insertion mutagen that disrupted the structure of the inserted host gene. Several changes were made in PB-SB-PGK-neo-bpA to convert it to an activating mutagen. The plasmid is first digested with HindIII restriction enzyme and calf intestinal phosphatase, and ligated with a PCR-amplified fragment containing the CMV enhancer and promoter sequence  and the splice donor from the rabbit beta-globin intron  to make pPB-SB-CMV-neo-SD. The pPB-SB-CMV-neo-SD plasmid was then digested with BglII and XmaI to remove the PGK-Neo-bpA cassette, and was ligated with a PCR-amplified SV40-driven puromycin cassette to provide a rapid selection marker to identify successful integrants. The final plasmid was sequence-verified and named pPB-SB-CMV-puro-SD.
Cell line transfection for library construction
To make a library, 1 × 107 cells were plated overnight in four T175 flasks at cell density of 1 × 105 cells per ml. HeLa and MCF7 were cultured in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with glutaMAX (Invitrogen) and 10% fetal bovine serum (FBS). T47D was cultured in RPMI with glutaMAX and 10% FBS. IMR32 was cultured in Eagle’s Minimum Essential Medium (EMEM) supplemented with 10% FBS. Cells were co-transfected with 36 μg pPB-SB-CMV-puro-SD and 36 μg pCMV-PBase plasmids using 216 μl Fugene 6 (Roche) and 4.5 ml serum-free OPTI-MEM. After three days, cells were treated with fresh media with 2 μg/ml puromycin and cultured for additional 7–10 days. Cells surviving antibiotics treatment were harvested and cryopreserved as transposon-tagged prescreened libraries. In total, eight independent libraries were constructed, two for each cell line. To measure the insertion numbers per cell, cells from the original prescreened HeLa library were diluted and plated in a 96-well plate at average one cell per well. Five single cell colonies were identified, expanded and harvested for analysis.
To determine transposition efficiency, cells were transfected as above. One day after transfection, one cell plate was trypsinized and re-plated to a 6-well plate at various dilution ratios. Cells were treated with puromycin three days after transfection until colonies could be stained with Methylene Blue for visual counting. Transposition efficiency was defined as the proportion of initially seeded cells that could form puromycin-selected colonies.
One million transposon-tagged cells from each library were plated in 100 mm tissue culture plates for drug treatment. Native untagged cells were similarly plated as study control. Paclitaxel dosages were 20 ng/ml for HeLa and MCF7, 15 ng/ml for T47D and 4 ng/ml for IMR32. Dosages were chosen as to sufficiently kill all parental cells within one week. Cells were treated until paclitaxel-resistant colonies were visible. Treatment time varied among cell lines depending on proliferation rates, and usually took ten days up to two weeks. Cells were then either harvested as resistant clones, or as resistant pools. To isolate resistant clones, colonies were picked from the drug-treated plates using 3 mm diameter cloning discs (Sigma), and expanded in 6-well plates in the presence of puromycin and paclitaxel. Cell clones exhibited stable resistance to both paclitaxel and puromycin, continuing to grow when retreated after 2 weeks culture in the absence of either agent. To harvest resistant pools, cells from the paclitaxel-treated plates were trypsinized and replated in the presence of puromycin and paclitaxel for one more week to remove any remaining non-resistant cells. These screens were performed on all eight libraries, including replicate screens for one library of each cell line.
Splinkerette PCR and nextgen sequencing for insertion site detection
Genomic DNA was harvested from samples using DNeasy Blood & tissue Kit (Qiagen). Insertion sites can be detected by splinkerette PCR, a modified version of ligation-mediated PCR . For the HeLa prescreened library, 3.3 μg genomic DNA was digested with 10 units of Csp6I (Fermentas) at 37°C for two hours, and ligated to 100 picomole double-stranded linker catalyzed by 2000 units of T4 DNA ligase at 16°C for overnight. The ligated sample was amplified with primers LP1 and PB51-IL in a 100 μl PCR reaction. Primer LP1 matches the linker sequences, and primer PB51-IL matches the transposon sequences. The thermo-cycling condition is the following: 3 min/94°C, 10 cycles of 15 sec/94°C; 30 sec/72°C with −1°C touchdown/cycle; 1 min/72°C, 20 cycles of 15 sec/94°C; 30 sec/62°C; 1 min/72°C, and 20 min/72°C. One microliter of the first PCR product was re-amplified in a 50 μl PCR using nested primers LP2a and PB52-ILa that contain Illumina single-end reaction adapter sequences for binding to the flowcell. Thermo-cycling condition was similar to that of the first PCR with 10 touchdown cycles and 10 regular cycles. Amplified products were purified using QIAquick PCR Purification Kit (Qiagen). For paclitaxel resistant pools and clones, 170ng genomic DNA was digested with 2 units of Csp6I and ligated to 10 picomole linkers. Up to 96 samples were processed with barcode linkers in a multi-well plate. Samples were pooled after PCR and purified. Sequencing was performed using Illumina HiSeq 50 Single Read following standard protocols except that sample loading density was reduced by 50% to avoid over-clustering due to the first 10 repetitive nucleotides. Each multiplexed cohort is loaded in one lane of the flow cell. Custom sequencing primer Seq-P1 matches the linker sequences prior to the barcodes and the read direction is opposite to CMV (Additional file 1: Table S1). Sequencing data were de-multiplexed and trimmed to remove the barcode plus 1 adjacent base remaining from ligation at the Csp6I half site, and any library adapter sequence present at the 3’ end of each read was removed. Reads of 7 bp or longer were retained and aligned to the hg19 reference genome using Bowtie alignment program , keeping only unique alignments placing the 5’ end of a trimmed read within 3 bp of a Csp6I site. All further analysis performed on the read counts at each Csp6I site.
TOPO cloning and sanger sequencing
Nested PCR products of resistant clones (7 from MCF7, 1 from HeLa, and 4 from T47D) were prepared as above and cloned into vector pCR2.1-TOPO (Invitrogen). Bacterial colonies were sequenced with primer PB5-ILseq (Additional file 1: Table S1) from the transposon side. Insertion sites were aligned using the BLAT function of the UCSC Genome Browser version hg19 (http://genome.ucsc.edu/cgi-bin/hgGateway).
Quantitation of mRNA expression
Total RNA was isolated using Qiagen RNeasy Mini kit. One microgram of total RNA was treated with RNase-free DNase to remove genomic DNA. First-strand cDNA was synthesized using Roche Transcriptor First Strand cDNA Synthesis kit, and quantitated by BIO-RAD SYBR Green. All reactions were normalized to actin.
Detection of chimeric mRNA
To detect the chimeric mRNA, mRNA was reverse-transcribed as above. 1 μl cDNA was PCR-amplified using a forward primer specific to the PB transposon sequence and the reverse primer matching the ABCB1 exon 3 sequence. The thermo-cycling conditions are: 3min/94°C, 30 cycles of 30 sec/94°C; 30 sec/55°C; 30 sec/72°C, and 5 min/72°C. PCR products were fractionated on 1.7% agarose gel. The control PCR used the primer pair provided in the cDNA synthesis kit to amplify the housekeeping gene hPBGD for 35 cycles with 50°C annealing temperature, and the PCR products were fractionated on a 3% gel.
Paclitaxel sensitivity assays
IMR32 Cells were reverse-transfected in a 96-well plate with either a control pCMV plasmid or pCMV6-ABCB1 plasmid (Origene). Each well contained 100 ng plasmid DNA, 0.3 μl Fugene 6 transfection reagent, and 10 μl OPTI-MEM, and 10,000 IMR32 cells in 100 μl antibiotics-free complete EMEM were seeded to each well. After two days, medium was replenished and cells were treated with serial-diluted paclitaxel for five days. Each sample was assayed with four replicate wells. Viability was measured by CellTiter-Glo (Promega) and data were processed using GraphPad Prism. Error bars represented standard error of means (SEM, n=4).
IMR32 cells were transfected with a control pCMV plasmid or pCMV6-ABCB1 plasmid respectively. Transfection was performed in a 6-well plate with each well containing 2 μg plasmid DNA, 6 μl Fugene 6 transfection reagent, 100 μl OPTI-MEM, and 200,000 IMR32 cells in 2 ml EMEM. After three days, cell were lysed with NP40 cell lysis buffer (Invitrogen) and sonicated to shear genomic DNA. Samples were diluted in SDS sample loading buffer, fractionated by SDS-PAGE (Bio-Rad), and transferred to polyvinylidene difluoride membrane. The membrane was blotted with MDR1/ABCB1 rabbit polyclonal antibody (Cell Signaling Technology #12273) diluted by 2,500-fold, and goat anti-rabbit IgG (Thermo Scientific #31460) diluted by 10,000-fold. For actin controls, the membrane was blotted with anti-actin rabbit monoclonal antibody diluted by 2,500-fold (Cell Signaling Technology #4970) and goat anti-rabbit IgG by 10,000-fold. Images were captured by G:Box (Syngene).
Statistical and bioinformatics methods
To identify potential insertion sites in analysis of resistant pools and clones, we first filtered sequencing data to exclude ‘background’ signal derived from contaminating non-resistant cells or the low incidence of PCR products from inappropriate linker reactions or PCR reactions. We assumed that such background signal would follow a Poisson distribution. This was supported by our observation that a frequency distribution of sequence analysis from resistant cells followed a bi-phasic distribution, with a large number of different sequences represented at low frequency (1–50 reads) which resembled a Poisson distribution, combined with a series distinct sequences present at high frequency (100 reads upwards). We selected sequences present at >100 reads for further analysis, which we estimate represents significant enrichment (p < 0.05) over background signal. For analysis of pools of resistant cells insertion sites and targeted genes were then compiled between all samples, removing any insertions seen twice in repeated analysis of the same sample. For analysis of sequences from resistant clones, samples were further filtered to identify the 1–10 sequences present at highest frequency in each clone, based on our previous analysis of the likely number of transposon insertions per cell. Clones were then clustered manually based on shared insertion sites, and any samples clearly derived from more than one clone excluded. Data were then visualized using Gene Pattern software (Broad Institute of MIT and Harvard).
We use the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool to perform functional analysis on genes enriched in the resistant pools (http://david.abcc.ncifcrf.gov). Only candidate genes identified above were used for analysis. Enriched genes were both listed as clusters and as an annotation chart.
To estimate the number of insertions needed to cover the genome, we assumed that only forward strand insertions within 64kb upstream could activate a gene based on our observation. We further postulated that the random event of integration within this 64kb region followed Poisson distribution. To achieve at least 1 insertion in 95% of total genes, the expected mean occurrence needed to be 3.0 [P (3.0, ≤0) = 0.05], which translated to 21.3kb gap between two insertions. Assuming genome size of 3 × 106 kb, 1× genome coverage would need 2.8 × 105 insertions. To achieve 2× coverage, the expected mean would be 4.75 [P (4.75, ≤1) = 0.05], equivalent to 4.4× 105 insertions.
Construction of gene activating transposons and generation of libraries of mutant cells
Characterization of insertion libraries
Although a previous smaller study reported a preference of PB for transcribed genes with 70 out of 104 insertions being intragenic , our study found that just 45.6% of total insertion sites were located within transcribed gene sequences. This observation was consistent with the fact that 40.8% of all TTAA sequences in the genome are intragenic, indicating that there was no major preference for the transposon to insert into transcribed sequences. In addition, particularly relevant for our gene activation strategy, given that our data (described below) indicate that the transposon can, at least in some instances, activate expression of genes at a range of up to 64kb, we found 63% of insertions were within 25kb of at least one gene, an arbitrary range we chose to assign genes to insertion sites. Furthermore, we found that the proportions of sense- versus antisense- strand insertions are equivalent, both for the 63% insertions and for all insertions, indicating that transcribed sequences did not affect insert orientations.
To gain comprehensive identification of insertions within individual cells, 5 clones from the Hela prescreened library were isolated and sequenced, using DNA barcoding (Figure 1D and Additional file 1: Table S1) to permit multiplexing of samples. We found that each colony contained between 1 and 11 insertions, with an average of 6 insertions (Figure 2B). Based on this result, we estimate that there may be up to 3.8 × 106 genomic insertion sites in our HeLa library of 6 × 105 independent clones. However, only a fraction of these insertions were revealed by our Illumina sequencing of the library, likely due to technical limitations of the amount of genomic DNA used as input or the efficiency of the PCR reactions.
The generation of cell clones also provided the opportunity to explore the effects of transposon insertion on gene expression, which is the key to our functional mutagenesis approach As illustrated by our analysis of ABCB1 in the next section, ‘sense’ insertions upstream of genes consistently resulted in increased expression as expected. In one clone in which the transposon inserted upstream of the gene in the reverse orientation, expression was also increased (Figure 2C ACADL). In contrast, intragenic insertion of the transposon caused decreased expression. Based on this characterization of individually targeted genes, we conclude that our ‘activation tagging’ approach will result in consistent strong stimulation of gene expression when inserted in the forward orientation upstream of genes, coupled with less predictable repression of expression for reverse direction and/or intragenic insertions.
Use of paclitaxel resistance screen to demonstrate the transposon functional mutagenesis approach
Identification of ABCB1as the benchmark resistant gene in all cell lines validated the mutagenesis screen
While several other genes were associated with multiple transposon insertions (Figure 3A), these genes were represented at significantly lower levels than ABCB1 (with none having greater than 6 independent insertions) and none were seen in all cell lines tested. To test whether the screen approach can enrich cell processes and pathways related to paclitaxel resistance, the 115 genes with ≥2 insertions were analyzed by functional and structural motifs using Database for Annotation, Visualization and Integrated Discovery tool (DAVID) (Additional file 4: Dataset S4) . There was strong enrichment for genes associated with microtubule components and cytoskeletal rearrangement, which are known paclitaxel targets (Figure 3C) [29–31]. Ion transport channels were likewise enriched, consistent with reports that ion channels utilize microfilaments for their function and are susceptible to paclitaxel, and that their expression can affect sensitivity to paclitaxel [32–38]. Thus, taken together, these results show that our transposon mutagenesis approach can readily identify major resistance mechanisms and provide potential insight into the biological processes targeted by the drug used to select resistant cells.
Use of clonal analysis to reveal gene interactions
Furthermore, the second most common hit in our pool analysis (Figure 3A), a transcription factor MEIS1, was only selected in IMR32 cells, and in clonal analysis was only seen in clones that also had insertions in ABCB1, implicating a possible role for MEIS1 in modifying ABCB1-mediated resistance, rather than inducing resistance alone. Insertion site orientation and positions of transposons suggested that enhanced resistance is associated with down regulation of MEIS1 expression (Figure 5D), although we were not able to directly verify this using siRNA-mediated gene knockdown (data not shown). Therefore, to look for independent evidence of MEIS1 function within ABCB1 context, we turned to our recently published database of drug sensitivity for a panel of cancer cell lines  and the publicly available Broad Institute Cancer Cell Line Encyclopedia (CCLE) microarray database . Our drug sensitivity database consists of 639 human cancer cell lines in combination with 130 targeted therapy or cytotoxic drugs, assayed in a 9-point 256-fold serial dilution setting. In total, 143 cell lines across diverse cancer types that have been assayed for paclitaxel sensitivity in our database overlapped with the CCLE cell line collection for gene expression, and were analyzed for correlation between paclitaxel sensitivity and expression of ABCB1 and MEIS1. As expected, there was a significant correlation between ABCB1 expression and paclitaxel sensitivity, but not between MEIS1 and paclitaxel sensitivity (Figure 5E). Instead, a negative correlation was observed between MEIS1 expression and paclitaxel sensitivity only in cell lines expressing high levels of ABCB1, but not in ABCB1-low cells, as predicted by clonal analysis. Although independent validation will be required to confirm the role of MEIS1 these data suggest that transposon activation mutagenesis and clonal analysis can be used to reveal interesting information such as primary resistant events and modifiers.
Acquired resistance to chemotherapeutic drugs remains a major hurdle to effective cancer treatment and eradication, and a better understanding of the genes and pathways that contribute to this is needed. Our data demonstrate that transposon mutagenesis provides a powerful, adaptable and cost-effective forward genetic approach for identifying resistance genes. The use of parallel screening in four separate cell lines with replicate samples to identify both common and cell-line restricted resistance gene candidates illustrates the potential of this system for gaining a deeper and more comprehensive view of resistance across the full spectrum of malignancy.
Transposon-based gene activation systems have a unique combination of properties that make them ideal for studies of tumor cell resistance. First, they are readily applied to gain-of-function genetic screens, unlike the vast majority of functional genetic interrogation approaches for mammalian genomes. This may be of particular relevance in tumors, where gene activation or amplification are common transforming events. The potential for long-range transcriptional activation and the need to hit only one allele of a gene in diploid cells provides higher effective genome coverage than gene inactivation strategies at comparable levels of mutagenesis. Furthermore, the presence of multiple effective mutations in single cells should allow for identification of ‘cooperating’ resistance genes, as suggested by our analysis of the ABCB1/MEIS1 interaction in resistant clones. A second advantage is the potential to identify resistance events occurring from changes in expression of uncharacterized or poorly understood genetic elements such as long intergenic non-coding (LINC) RNA or microRNAs. We have identified resistant cells bearing insertions in LINC-RNAs and unannotated transcripts and further investigation of these mutations could shed light on new genetic elements. Third, this system is readily transferable to new cell lines, including cells derived from patient samples, and libraries can be expanded and regenerated simply by re-transfection with the transposase. Hence transposon-based screens allow for the rapid generation of resistant clones to single drugs or combined therapies in specific tumor cells, providing insights into potential resistance mechanisms. These could then be used to guide design of new drug combinations and tailor treatments to particular tumor types.
Transposon-based screening has been used previously to identify potential mechanisms of resistance to the antibiotic puromycin and chemotherapeutic, vincristine . In those studies, transposon insertions were found primarily in Abcb1a/b (both drugs) and the closely related transporter Abcg2 (puromycin only), reinforcing our findings that overexpression of ABCB1 represents a major mechanism of drug resistance. However, additional candidate genes were not identified; this may be due to the absence of a splice donor in the transposon used, limiting the ability of the inserted promoter to activate gene expression, the presence of only one transposon per cell, or the relatively limited analysis of insertion events using capillary sequencing of isolated cell clones
Based on our finding that the transposon used here can exert a strong transcriptional activation effect at 64kb upstream from open reading frames, we estimate (using Poisson distribution) that libraries consisting of just 4.7 × 104 clones (2.8 × 105 insertions) or 7.3 × 104 clones (4.4 × 105 insertions) respectively could be potentially capable of activating 95% of genes by delivering at least one or multiple forward upstream insertions. Therefore in the case of HeLa cells we have approached meaningful close to genome-wide coverage for activation events. This is supported by our results with ABCB1, from which it is clear that all of our libraries had sufficient coverage to provide multiple insertions in a single strong resistance gene. We deduced that the low incidence of identification of other resistance genes could therefore reflect real differences in the ‘potency’ of individual genes to promote resistance, with only ABCB1 being sufficient while other events requiring additive effects of multiple genes to yield resistance in the high taxol concentration used here for selection. This is also supported by a prior transposon screen [41 described above] which identified only ABC-family transporters as potential resistance genes. Support for this also comes from our bioinformatics analysis, which revealed concordance of gene function or pathways between candidate resistance genes, with a strong enrichment in microtubule related biology previously linked to paclitaxel resistance.
However, even considering these potential limitations, our data identify new possible resistance genes and strengthen the evidence for previously identified candidates. As an example, clonal analysis of resistant cells strongly implicate MEIS1 as a modifier of ABCB1-mediated resistance, and this is further supported by our analysis of a large panel of tumor cells. MEIS1 is a class A homeodomain protein that acts as a cofactor for homeobox (HOX) proteins, and has been implicated as a critical downstream target of oncogenic fusion proteins in leukemia [42–45]. Although high expression promotes leukemia cell proliferation, silencing of MEIS1 increases resistance to the chemotherapeutic etoposide , in agreement with our findings in IMR32 cells. Our clonal analysis of mutations also reveal CXCR4, the receptor for the chemokine CXCL12 (stromal derived growth factor 1, SDF-1), as a potential resistance gene that functions independently of ABCB1 (Figure 5A). Up-regulation of CXCR4 is associated with increased metastasis and poor prognosis in various forms of cancer, in part due to effects on cellular phenotype, and is associated with chemotherapeutic resistance in numerous tumor models. For example, CXCR4 is upregulated in gefitinib-resistant non-small cell lung cancer cells and promotes epithelial-mesenchymal transition (EMT) and self-renewal activity . Likewise, CD133+ glioblastoma cancer stem cells with increased resistance to a range of chemotherapeutic agents, including paclitaxel, express high levels of CXCR4 , and high surface expression of this chemokine receptor is considered a marker of cancer stem cells . Of direct relevance to our results, increased CXCR4 expression and CXCL12/ CXCR4 signaling promote tumor cell resistance to the chemotherapeutic gemcitabine .
Other genes identified in our screen, such as ALK and PDE4D, increase tumor cell growth and protect from apoptosis, and may therefore promote resistance through these mechanisms.
Finally, further genes identified from independent hits in different cell lines, including the protocadherin PCDH15 and neuroblastoma breakpoint family member NBPF11 have not been implicated in tumor resistance, but based on the examples described above, may also have important roles.
We have developed a transposon-mediated activation mutagenesis and screen approach to systematically identify chemotherapy resistance. We demonstrated the feasibility of using this approach by identifying genes and pathways related to paclitaxel resistance. This system provided unbiased genome-wide coverage with sufficient depth to reliably capture the most well characterized mechanism of resistance in all cell lines processed and generate many believable additional hits based on available functional annotation. In addition to a dramatically lower cost and higher efficiency over RNAi and cDNA libraries, this approach does not require a priori knowledge of candidate genes, can survey untranscribed regions and can generate stable resistant clones pertinent to specific cell lines. Although further analysis will be required to dissect candidate genes, our findings highlight the potential for transposon-based functional genetics to aid in identifying both novel resistance genes and gene combinations. These may allow improved selection of chemotherapeutic drugs for particular classes of tumors, or the characterization of “resistance gene signatures” for new or ‘black box’ targeted therapeutics, allowing development of combination therapies to overcome potential resistance, and improve the efficacy and duration of new cancer therapies.
We thank Dr. Pentao Liu and Allan Bradley of Wellcome Trust Sanger Institute for gifting plasmids PB-SB-PGK-neo-bpA and pCMV-PBase.
This work was supported by National Institutes of Health grant (U01-AI070330). LC, AD, and CB received additional funding from the Wellcome Trust (086357).
- Meijer D, van Agthoven T, Bosma PT, Nooter K, Dorssers LC: Functional screen for genes responsible for tamoxifen resistance in human breast cancer cells. Mol Cancer Res. 2006, 4 (6): 379-386. 10.1158/1541-7786.MCR-05-0156.View ArticlePubMedGoogle Scholar
- Johannessen CM, Boehm JS, Kim SY, Thomas SR, Wardwell L, Johnson LA, Emery CM, Stransky N, Cogdill AP, Barretina J: COT drives resistance to RAF inhibition through MAP kinase pathway reactivation. Nature. 2010, 468 (7326): 968-972. 10.1038/nature09627.View ArticlePubMedPubMed CentralGoogle Scholar
- Gaither LA: Chemogenomics approaches to novel target discovery. Expert Rev Proteomics. 2007, 4 (3): 411-419. 10.1586/14789422.214.171.1241.View ArticlePubMedGoogle Scholar
- Kool J, Berns A: High-throughput insertional mutagenesis screens in mice to identify oncogenic networks. Nat Rev Cancer. 2009, 9 (6): 389-399.View ArticlePubMedGoogle Scholar
- Dorssers LC, van Agthoven T, Dekker A, van Agthoven TL, Kok EM: Induction of antiestrogen resistance in human breast cancer cells by random insertional mutagenesis using defective retroviruses: identification of bcar-1, a common integration site. Mol Endocrinol. 1993, 7 (7): 870-878. 10.1210/me.7.7.870.PubMedGoogle Scholar
- Miething C, Grundler R, Mugler C, Brero S, Hoepfl J, Geigl J, Speicher MR, Ottmann O, Peschel C, Duyster J: Retroviral insertional mutagenesis identifies RUNX genes involved in chronic myeloid leukemia disease persistence under imatinib treatment. Proc Natl Acad Sci U S A. 2007, 104 (11): 4594-4599. 10.1073/pnas.0604716104.View ArticlePubMedPubMed CentralGoogle Scholar
- Singhal R, Deng X, Chenchik AA, Kandel ES: Long-distance effects of insertional mutagenesis. PLoS One. 2011, 6 (1): e15832-10.1371/journal.pone.0015832.View ArticlePubMedPubMed CentralGoogle Scholar
- Carette JE, Raaben M, Wong AC, Herbert AS, Obernosterer G, Mulherkar N, Kuehne AI, Kranzusch PJ, Griffin AM, Ruthel G: Ebola virus entry requires the cholesterol transporter Niemann-Pick C1. Nature. 2011, 477 (7364): 340-343. 10.1038/nature10348.View ArticlePubMedPubMed CentralGoogle Scholar
- McClintock B: The origin and behavior of mutable loci in maize. Proc Natl Acad Sci U S A. 1950, 36 (6): 344-355. 10.1073/pnas.36.6.344.View ArticlePubMedPubMed CentralGoogle Scholar
- Ivics Z, Hackett PB, Plasterk RH, Izsvak Z: Molecular reconstruction of Sleeping Beauty, a Tc1-like transposon from fish, and its transposition in human cells. Cell. 1997, 91 (4): 501-510. 10.1016/S0092-8674(00)80436-5.View ArticlePubMedGoogle Scholar
- Elick TA, Bauser CA, Fraser MJ: Excision of the piggyBac transposable element in vitro is a precise event that is enhanced by the expression of its encoded transposase. Genetica. 1996, 98 (1): 33-41. 10.1007/BF00120216.View ArticlePubMedGoogle Scholar
- Fraser MJ, Ciszczon T, Elick T, Bauser C: Precise excision of TTAA-specific lepidopteran transposons piggyBac (IFP2) and tagalong (TFP3) from the baculovirus genome in cell lines from two species of Lepidoptera. Insect Mol Biol. 1996, 5 (2): 141-151. 10.1111/j.1365-2583.1996.tb00048.x.View ArticlePubMedGoogle Scholar
- Ding S, Wu X, Li G, Han M, Zhuang Y, Xu T: Efficient transposition of the piggyBac (PB) transposon in mammalian cells and mice. Cell. 2005, 122 (3): 473-483. 10.1016/j.cell.2005.07.013.View ArticlePubMedGoogle Scholar
- Chew SK, Rad R, Futreal PA, Bradley A, Liu P: Genetic screens using the piggyBac transposon. Methods. 2011, 53 (4): 366-371. 10.1016/j.ymeth.2010.12.022.View ArticlePubMedPubMed CentralGoogle Scholar
- Ivics Z, Li MA, Mates L, Boeke JD, Nagy A, Bradley A, Izsvak Z: Transposon-mediated genome manipulation in vertebrates. Nat Methods. 2009, 6 (6): 415-422. 10.1038/nmeth.1332.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang W, Bradley A, Huang Y: A piggyBac transposon-based genome-wide library of insertionally mutated Blm-deficient murine ES cells. Genome Res. 2009, 19 (4): 667-673. 10.1101/gr.085621.108.View ArticlePubMedPubMed CentralGoogle Scholar
- Horie K, Yusa K, Yae K, Odajima J, Fischer SE, Keng VW, Hayakawa T, Mizuno S, Kondoh G, Ijiri T: Characterization of Sleeping Beauty transposition and its application to genetic screening in mice. Mol Cell Biol. 2003, 23 (24): 9189-9207. 10.1128/MCB.23.24.9189-9207.2003.View ArticlePubMedPubMed CentralGoogle Scholar
- Mates L, Izsvak Z, Ivics Z: Technology transfer from worms and flies to vertebrates: transposition-based genome manipulations and their future perspectives. Genome Biol. 2007, 8 Suppl 1: S1-View ArticlePubMedGoogle Scholar
- Walden R, Fritze K, Hayashi H, Miklashevichs E, Harling H, Schell J: Activation tagging: a means of isolating genes implicated as playing a role in plant growth and development. Plant Mol Biol. 1994, 26 (5): 1521-1528. 10.1007/BF00016488.View ArticlePubMedGoogle Scholar
- Collier LS, Carlson CM, Ravimohan S, Dupuy AJ, Largaespada DA: Cancer gene discovery in solid tumours using transposon-based somatic mutagenesis in the mouse. Nature. 2005, 436 (7048): 272-276. 10.1038/nature03681.View ArticlePubMedGoogle Scholar
- Vassiliou GS, Cooper JL, Rad R, Li J, Rice S, Uren A, Rad L, Ellis P, Andrews R, Banerjee R: Mutant nucleophosmin and cooperating pathways drive leukemia initiation and progression in mice. Nat Genet. 2011, 43 (5): 470-475. 10.1038/ng.796.View ArticlePubMedPubMed CentralGoogle Scholar
- Boshart M, Weber F, Jahn G, Dorsch-Hasler K, Fleckenstein B, Schaffner W: A very strong enhancer is located upstream of an immediate early gene of human cytomegalovirus. Cell. 1985, 41 (2): 521-530. 10.1016/S0092-8674(85)80025-8.View ArticlePubMedGoogle Scholar
- Koh J, Enders GH, Dynlacht BD, Harlow E: Tumour-derived p16 alleles encoding proteins defective in cell-cycle inhibition. Nature. 1995, 375 (6531): 506-510. 10.1038/375506a0.View ArticlePubMedGoogle Scholar
- Jao LE, Maddison L, Chen W, Burgess SM: Using retroviruses as a mutagenesis tool to explore the zebrafish genome. Brief Funct Genomic Proteomic. 2008, 7 (6): 427-443. 10.1093/bfgp/eln038.View ArticlePubMedPubMed CentralGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang W, Lin C, Lu D, Ning Z, Cox T, Melvin D, Wang X, Bradley A, Liu P: Chromosomal transposition of PiggyBac in mouse embryonic stem cells. Proc Natl Acad Sci U S A. 2008, 105 (27): 9290-9295. 10.1073/pnas.0801017105.View ArticlePubMedPubMed CentralGoogle Scholar
- Gottesman MM, Fojo T, Bates SE: Multidrug resistance in cancer: role of ATP-dependent transporters. Nat Rev Cancer. 2002, 2 (1): 48-58. 10.1038/nrc706.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.View ArticlePubMedGoogle Scholar
- Zhang M, Wang W: Organization of signaling complexes by PDZ-domain scaffold proteins. Accounts Chem Res. 2003, 36 (7): 530-538. 10.1021/ar020210b.View ArticleGoogle Scholar
- Ranganathan R, Ross EM: PDZ domain proteins: scaffolds for signaling complexes. Curr Biol. 1997, 7 (12): R770-773. 10.1016/S0960-9822(06)00401-5.View ArticlePubMedGoogle Scholar
- Fanning AS, Anderson JM: Protein-protein interactions: PDZ domain networks. Curr Biol. 1996, 6 (11): 1385-1388. 10.1016/S0960-9822(96)00737-3.View ArticlePubMedGoogle Scholar
- Chen SZ, Jiang M, Zhen YS: HERG K+ channel expression-related chemosensitivity in cancer cells and its modulation by erythromycin. Cancer Chemother Pharmacol. 2005, 56 (2): 212-220. 10.1007/s00280-004-0960-5.View ArticlePubMedGoogle Scholar
- Kornau HC, Seeburg PH, Kennedy MB: Interaction of ion channels and receptors with PDZ domain proteins. Curr Opin Neurobiol. 1997, 7 (3): 368-373. 10.1016/S0959-4388(97)80064-5.View ArticlePubMedGoogle Scholar
- Stoeckel H, Takeda K: Plasmalemmal voltage-activated K(+) currents in protoplasts from tobacco BY-2 cells: possible regulation by actin microfilaments?. Protoplasma. 2002, 220 (1–2): 79-87.View ArticlePubMedGoogle Scholar
- Casini S, Tan HL, Demirayak I, Remme CA, Amin AS, Scicluna BP, Chatyan H, Ruijter JM, Bezzina CR, van Ginneken AC: Tubulin polymerization modifies cardiac sodium channel expression and gating. Cardiovasc Res. 2009, 85 (4): 691-700.View ArticlePubMedGoogle Scholar
- Miyano K, Tang HB, Nakamura Y, Morioka N, Inoue A, Nakata Y: Paclitaxel and vinorelbine, evoked the release of substance P from cultured rat dorsal root ganglion cells through different PKC isoform-sensitive ion channels. Neuropharmacology. 2009, 57 (1): 25-32. 10.1016/j.neuropharm.2009.04.001.View ArticlePubMedGoogle Scholar
- Montalbetti N, Li Q, Wu Y, Chen XZ, Cantiello HF: Polycystin-2 cation channel function in the human syncytiotrophoblast is regulated by microtubular structures. J Physiol. 2007, 579 (Pt 3): 717-728.View ArticlePubMedPubMed CentralGoogle Scholar
- Sonee M, Barron E, Yarber FA, Hamm-Alvarez SF: Taxol inhibits endosomal-lysosomal membrane trafficking at two distinct steps in CV-1 cells. Am J Physiol. 1998, 275 (6 Pt 1): C1630-1639.PubMedGoogle Scholar
- Garnett MJ, Edelman EJ, Heidorn SJ, Greenman CD, Dastur A, Lau KW, Greninger P, Thompson IR, Luo X, Soares J: Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature. 2012, 483 (7391): 570-575. 10.1038/nature11005.View ArticlePubMedPubMed CentralGoogle Scholar
- Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehar J, Kryukov GV, Sonkin D: The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012, 483 (7391): 603-607. 10.1038/nature11003.View ArticlePubMedPubMed CentralGoogle Scholar
- Kong J, Wang F, Brenton JD, Adams DJ: Slingshot: a PiggyBac based transposon system for tamoxifen-inducible 'self-inactivating' insertional mutagenesis. Nucleic Acids Res. 2010, 38 (18): e173-10.1093/nar/gkq658.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith LL, Yeung J, Zeisig BB, Popov N, Huijbers I, Barnes J, Wilson AJ, Taskesen E, Delwel R, Gil J: Functional crosstalk between Bmi1 and MLL/Hoxa9 axis in establishment of normal hematopoietic and leukemic stem cells. Cell Stem Cell. 2011, 8 (6): 649-662. 10.1016/j.stem.2011.05.004.View ArticlePubMedGoogle Scholar
- Zeisig BB, Milne T, Garcia-Cuellar MP, Schreiner S, Martin ME, Fuchs U, Borkhardt A, Chanda SK, Walker J, Soden R: Hoxa9 and Meis1 are key targets for MLL-ENL-mediated cellular immortalization. Mol Cell Biol. 2004, 24 (2): 617-628. 10.1128/MCB.24.2.617-628.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Lawrence HJ, Rozenfeld S, Cruz C, Matsukuma K, Kwong A, Komuves L, Buchberg AM, Largman C: Frequent co-expression of the HOXA9 and MEIS1 homeobox genes in human myeloid leukemias. Leukemia. 1999, 13 (12): 1993-1999. 10.1038/sj.leu.2401578.View ArticlePubMedGoogle Scholar
- Kroon E, Krosl J, Thorsteinsdottir U, Baban S, Buchberg AM, Sauvageau G: Hoxa9 transforms primary bone marrow cells through specific collaboration with Meis1a but not Pbx1b. EMBO J. 1998, 17 (13): 3714-3725. 10.1093/emboj/17.13.3714.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosales-Avina JA, Torres-Flores J, Aguilar-Lemarroy A, Gurrola-Diaz C, Hernandez-Flores G, Ortiz-Lazareno PC, Lerma-Diaz JM, de Celis R, Gonzalez-Ramella O, Barrera-Chaires E: MEIS1, PREP1, and PBX4 are differentially expressed in acute lymphoblastic leukemia: association of MEIS1 expression with higher proliferation and chemotherapy resistance. J Exp Clin Cancer Res. 2011, 30: 112-10.1186/1756-9966-30-112.View ArticlePubMedPubMed CentralGoogle Scholar
- Jung MJ, Rho JK, Kim YM, Jung JE, Jin YB, Ko YG, Lee JS, Lee SJ, Lee JC, Park MJ: Upregulation of CXCR4 is functionally crucial for maintenance of stemness in drug-resistant non-small cell lung cancer cells. Oncogene. 2013, 32 (2): 209-221. 10.1038/onc.2012.37.View ArticlePubMedGoogle Scholar
- Liu G, Yuan X, Zeng Z, Tunici P, Ng H, Abdulkadir IR, Lu L, Irvin D, Black KL, Yu JS: Analysis of gene expression and chemoresistance of CD133+ cancer stem cells in glioblastoma. Mol Cancer. 2006, 5: 67-10.1186/1476-4598-5-67.View ArticlePubMedPubMed CentralGoogle Scholar
- Miki J, Furusato B, Li H, Gu Y, Takahashi H, Egawa S, Sesterhenn IA, McLeod DG, Srivastava S, Rhim JS: Identification of putative stem cell markers, CD133 and CXCR4, in hTERT-immortalized primary nonmalignant and malignant tumor-derived human prostate epithelial cell lines and in prostate cancer specimens. Cancer Res. 2007, 67 (7): 3153-3161. 10.1158/0008-5472.CAN-06-4429.View ArticlePubMedGoogle Scholar
- Singh S, Srivastava SK, Bhardwaj A, Owen LB, Singh AP: CXCL12-CXCR4 signalling axis confers gemcitabine resistance to pancreatic cancer cells: a novel target for therapy. Br J Cancer. 2010, 103 (11): 1671-1679. 10.1038/sj.bjc.6605968.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/93/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.