A genome scale overexpression screen to reveal drug activity in human cells
© Arnoldo et al.; licensee BioMed Central Ltd. 2014
Received: 3 September 2013
Accepted: 22 April 2014
Published: 29 April 2014
Target identification is a critical step in the lengthy and expensive process of drug development. Here, we describe a genome-wide screening platform that uses systematic overexpression of pooled human ORFs to understand drug mode-of-action and resistance mechanisms. We first calibrated our screen with the well-characterized drug methotrexate. We then identified new genes involved in the bioactivity of diverse drugs including antineoplastic agents and biologically active molecules. Finally, we focused on the transcription factor RHOXF2 whose overexpression conferred resistance to DNA damaging agents. This approach represents an orthogonal method for functional screening and, to our knowledge, has never been reported before.
Biological systems tend to remain phenotypically stable in the face of environmental challenges and genetic changes. As such, genetic perturbation has become an efficient technique to dissect cellular functions. In this study, we used the concept of modulating gene dosage in human cells to gain insight into drug mode of action. Understanding the primary mechanism of action as well as the potential polypharmacological effects of a drug can provide insight into how to reduce detrimental side effects, uncover new applications for novel indication and explain resistance mechanisms[2–4].
Drug resistance in malignant tissues can be categorized into three main mechanisms: (i) drug distribution/metabolism (pharmacokinetics), (ii) heterogeneity of cancer cells and (iii) tumor micro-environment[5, 6]. Among the other cellular mechanisms, gene overexpression (for example, by amplification of the drug target) can titrate a drug’s effect. This is exemplified by the classic case of methotrexate resistance through the amplification of the gene DHFR in neoplastic tissue from an individual with disseminated small-cell lung cancer that relapsed during methotrexate chemotherapy. Other forms of overexpression resistance include the up-regulation of 1) efflux pumps (for example, ATP-binding cassette transporters), 2) survival mechanisms (for example, anti-apoptotic proteins), 3) DNA damage repair, 4) pathways for drug inactivation and 5) the overexpression of target isotypes. Additionally, proteins downstream of the inhibited target can be modulated in such a manner as to bypass the toxic effect of a drug.
To dissect some of these mechanisms by which drugs act within the cell, we postulated that, when overexpressed, genes conferring resistance to a lethal chemical treatment can illuminate the drug mode of action. Overexpressing genes to confer resistance to an otherwise toxic compound is a well-established concept. Early studies showed that a streptomycin resistance gene cassette could function in a bacterial tetracycline resistant plasmid and it was shown in other studies that cloning of the mouse dihydrofolate reductase into a bacterial plasmid provided resistance to trimethoprim in Escherichia coli. A large-scale approach using the overexpression of a yeast genomic DNA library to identify genes conferring resistance to specific drugs has validated the concept of gene dosage for the discovery of drug targets in eukaryotes.
A logical extension of this concept is to employ newly available biological tools that enable the systematic perturbation of all protein function in human cultured cells using RNA interference (RNAi), CRISPR (clustered regularly interspaced short palindromic repeats), cDNA libraries, transposons or small molecule inhibitors in combination with gene overexpression. Although testing all important therapeutic drugs against all proteins from the approximately 20,000 genes of the human genome in all differentiated cell types is not currently feasible, we report the initial development of a gain-of-function approach to identify drug mode of action by the overexpression of 12,200 human ORFs (hORFs)[11, 12] in the human cell line HEK293.
In summary, we have developed a new experimental pipeline to identify genes whose up-regulation suppresses the toxic effect of chemicals in human cells and have successfully applied this strategy to seven pharmacological compounds.
Cell culture and plasmids
Human embryonic kidney HEK293 (ATCC), HEK293T (ATCC), non-small cell lung carcinoma NCI-H1299 (ATCC), chondrosarcoma SW1353 (gift from Johanne Martel-Pelletier, University of Montreal, Canada) and pancreatic adenocarcinoma HPAC (ATCC) cells were maintained in DMEM (Wisent Inc. Montreal, Quebec, Canada) supplemented with 10% fetal bovine serum (Gibco, Carlsbad, California, US) at 37°C and 5% CO2. The rtTA expressing HEK293_M2 (ATCC), human breast adenocarcinoma MCF7_M2 (ATCC) and adenocarcinomic human alveolar basal epithelial A549_M2 (ATCC) cells were maintained in DMEM (Wisent Inc.) supplemented with 10% fetal bovine serum (Gibco) at 37°C and 5% CO2. Blast phase chronic myelogenous leukemia K562 cells (gift from Reinhart Reithmeier, University of Toronto, Canada) were maintained in RPMI1640 (Sigma, St. Louis, Missouri, US) supplemented with 10% fetal bovine serum (Gibco) at 37°C and 5% CO2. Multiple myeloma U266B1 cells (gift from Aaron Schimmer, Princess Margaret Cancer Centre, Canada) were maintained in Iscove's modified Dulbecco's medium (IMDM; Gibco) supplemented with 10% fetal bovine serum (Gibco) at 37°C and 5% CO2. Thyroid gland medullary carcinoma TT cells (gift from Gilbert Cote, MD Anderson Cancer Center, USA) were maintained in Ham’s F-12 K (Gibco) supplemented with 10% fetal bovine serum (Gibco) at 37°C and 5% CO2.
The lentiviral destination vector pLD-T-IRES-Venus-WPRE-STOP was cloned using the following procedure: the Gateway cassette from pLS-Dest-EcF (gift of Dr Tony Pawson, University of Toronto, Canada) was subcloned into pLJM17 (J Moffat, unpublished) using the restriction enzymes MluI and XbaI, resulting in the pLD-puro-TRE (T) vector. IRES-Venus-WPRE were PCR amplified from the pSLIK-Venus using the following primers: XcmI_IRES_F (5′-GCGCCTTTTCCAAGGCAGCCCTGGAATTCCGCCCCTCTCCCTCC) and NsiI_WPRE_R (5′-AAACAATGCATGTCGACGCGGGGAGGCGGCCCAAAGGGAGATCC). The IRES-Venus-WPRE amplicon was cloned into the pLD-puro-T vector using the restriction enzymes XcmI and NsiI to replace the human phosphoglycerate kinase promoter and puromycin resistance gene. A stop codon was introduced directly downstream of the Gateway cassette by digesting pLD-T-IRES-Venus-WPRE with the restriction enzyme XbaI and blunted using DNA polymerase I, large (Klenow) fragment (New England Biolabs Inc., Beverly, Massachusetts, US).
The piggyBac PB-TGcMV-Neo plasmid used for the hit confirmation was derived from the PB-TET (AddGene, Cambridge, Massachusetts, US). The 'IRES-beta-Geo' fragment from PB-TET was replaced by the 'promoter PGK-neomycin resistance gene' fragment by homologous recombination in Saccharomyces cerevisiae. Briefly, PB-TET was digested with RsrII and ApaI to remove the beta-geo gene. The fragment coding for the promoter PGK and neomycin resistance gene, the URA3 cassette and the yeast 2 μ-origin of replication were PCR amplified and cloned into the digested PB-TET plasmid by homologous recombination in yeast. Yeast was then transformed with the p414-Cre plasmid in order to excise the URA3 cassette and the 2 μ-origin of replication after expression of the Cre recombinase.
All short hairpin RNAs (shRNAs) against the RHOXF2 and green fluorescent protein (GFP) mRNAs were derived from the RNAi Consortium (TRC) lentiviral libraries and obtained in the pLKO.1 vector. Five shRNAs against RHOXF2 were tested and were designed against the following regions: 5′-CGGGATGAGAGATGATTACTT (RHOXF2_shRNA_1); 5′-GACGAGAAAGAACTACAGGAT (RHOXF2_shRNA_2); 5′-AGAAGCATGAATGTGACTGAA (RHOXF2_shRNA_3); 5′-GAGGGCATTAATGGCAAGAAA (RHOXF2_shRNA_4); 5′-GTCGCTTACTGAAGAGGTCAA (RHOXF2_shRNA_5).
Library preparation, virus production and lentiviral infection
The 12,212 hORF collection (representing 10,214 distinct genes) represents version 3.1 of an ongoing effort to create a complete human set of protein-encoding genes[11, 12]. Initially based on the Mammalian Gene Collection cDNA collection, that work transfers full-length ORFs (excluding the 5′ and 3′ mRNA untranslated regions) into a Gateway system. Version 3.1 was obtained from Open Biosystems (Thermo Scientific, Waltham, Massachusetts, US) and regroups non-fully sequenced verified and non-clonal ORFs (limitations corrected in the 8.1 version). The collection was divided into 34 minipools of 376 hORFs; for each minipool, the hORFs were cloned en masse from the pDONR223 into the lentiviral expression vector pLD-T-IRES-Venus-WPRE-STOP by Gateway LR reaction. After electroporation, transformants were selected on Luria Broth (LB) plus ampicillin and the lentiviral vector was extracted.
Lentivirus was produced by normalizing the amount of DNA for the 34 hORF minipools and by co-transfection with the packaging plasmid psPAX2 and the envelope plasmid pMD2.G into the packaging HEK293T cells using FuGENE (Roche Mississauga, Ontario, Canada).
HEK293_M2 cells were then infected at a multiplicity of infection of 0.3. Doxycycline (2 μg/ml) was added to induce the expression of the Venus fluorescent protein and Venus-positive cells were sorted using a BD FACS Aria Cell Sorter (East Ruherford, New Jersey, US). Because the expression of the gene coding for Venus is linked to the expression of the hORFs, selection of the fluorescent cells allowed the selection for non-silent, stably integrated lentivirus as well as functional inducible hORFs.
Genomic DNA extraction, library preparation and data analysis for the hORFeome representation in HEK293_M2 cells
We thawed 6.3 × 106 HEK293_M2 cells harboring the virally integrated human ORFeome version 3.1 collection and these were cultured for 1 day to recover. Genomic DNA was recovered from 16 × 106 cells and 6 μg of genomic DNA (gDNA were used to amplify the collection (three PCR reactions with 2 μg of gDNA). Following PCR, the amplified ORFs were used directly to prepare a Nextera sequencing library (Illumina Nextera DNA Sample Preparation Kit, San Diego, California, US) and sequenced on an Illumina MiSeq according to the manufacturer’s protocol. Fifty-nucleotide paired-end reads were aligned against a reference database consisting of the ORF sequences in the hORFeome 3.1 collection using bwa (version 0.6.1). Coverage at each base position in each reference sequence was determined using the bedtools program genomeCoverageBed (v2.14.2). A custom Perl script was used to summarize the coverage levels to determine the mean coverage for each reference sequence and the proportion of each reference sequence with different levels of coverage. Raw data have been deposited in the ArrayExpress repository under accession number E-MTAB-2498.
hORFeome drug target screen
Aliquots of 500,000 HEK293_M2 cells were seeded in a T175 flask and grown for 15 hours before the induction of hORF expression by addition of doxycycline. After 18 hours of induction, cells formed micro-colonies that were cultured in the presence of drug or DMSO (as control) for 2 to 3 weeks until distinct colonies were detectable.
After each screen, surviving cells were harvested and gDNA was extracted as follows. Cells were lysed using SNET buffer (10 mM Tris pH 8.0, 0.1 M EDTA, 0.5% SDS, 0.1 mg/ml Proteinase K, 25 μg/ml RNAse A). gDNA was isolated using one volume of phenol/chloroform/isoamyl alcohol (25:24:1) pH 8.0, precipitated by adding two volumes of 95% ethanol (-20°C), washed with ice cold 70% ethanol and resuspended in Tris-HCl pH 8.0.
During the screen selection experiments, only a minority of hORFs protected HEK293_M2s against the cytotoxic effect of the drug. Compared to the initial cell population, HEK293_M2 cells harboring the most resistant genes should be over-represented. In order to further amplify the set of resistant genes from the resulting reduced population, hORFs were PCR amplified from only 80 ng of genomic DNA. The sequences of the forward and reverse primers specific for the inserted hORFs were 5′-CGGTACCCGGGGATCCTCTAGTCAGCTGAC and 5′-CCATTTGTCTCGAGGTCGAGAATTCTAGCTAGAATC, respectively. Each reaction was carried out in a 50 μl volume containing 25 μl of 2× Phusion Flash High-Fidelity Master Mix (Finnzymes, Espoo, Finland), 200 nM of each primer and 80 ng of genomic DNA. The PCR profile was: 1 minute at 98°C for one cycle; 10 s at 98°C, 20 s at 65°C, 4 minutes at 72°C for 35 cycles; 10 minutes at 72°C for one cycle. PCR products were purified using QIAquick PCR Purification Kit (QIAGEN Venlo, Limburg, Netherlands), 150 ng of purified PCR product was biotinylated using the BioPrime DNA Labeling Kit (Invitrogen, Carlsbad, California, US) and unincorporated biotin-14-dCTP was removed by passing the samples through Sephadex G-50 columns (GE Healthcare Fairfield, Connecticut, US). Sample (150 ng) was added to the Affymetrix GeneChip Human Gene 1.0 ST array and hybridized at 45°C for 17 hours with a rotation of 60 rpm. Chips were washed and stained with SAPE (2× MES staining buffer, 20 mg/mL bovine serum albumin (BSA) and 1 mg/ml streptavidin-phycoerythrin), washed on an Affymetrix fluidics station and scanned.
Representation of the complete hORFeome collection was assessed using the same experimental conditions with some modifications. gDNA was extracted from at least 16 × 106 cells (approximately 1,300 cells per hORF), 3 PCR reactions were performed using 2 μg of gDNA each (approximately 75 genomes per hORFs), two biotinylation reactions were performed using 500 ng of combined PCR products and 3.5 μg of sample was added to an Affymetrix GeneChip Human Gene 1.0 ST array.
Microarray raw data have been deposited in the ArrayExpress repository under accession number E-MTAB-2493.
Selection and individual validation of over-represented hORFs
The hORFs conferring resistance to the drug were identified in a step-wise procedure. Screen results were first visualized by plotting the log2 of the signal intensity for each hORF retrieved via PCR from cells cultured in the presence of drug (on the x-axis) and by plotting the log2 ratio of the signal intensity for each hORF of the cells cultured in the presence of drug divided by the signal intensity for each hORF of the cells grown in the presence of DMSO (on the y-axis). Genes with a log2(drug/DMSO) >3 and a log2(drug) >6 enrichment were then listed and ranked. Primary hits selected for validation were those previously enriched genes in common among the drug screen replicates (at least two independent experiments).
For validation, 10,000 HEK293-M2 cells were seeded in 300 μl of DMEM plus 10% iFBS plus penicillin/streptomycin in a 48-well plate pretreated with poly-L-lysine (Sigma, P4832). After 12 hours of incubation, gene expression was induced by the addition of doxycycline (2 μg/ml final concentration). After 18 hours of induction, drug was added and cell growth was assessed 3 days after drug addition by sulforhodamine B (SRB) assay as previously described.
The HEK293_M2 + RHOXF2 stable cell line (PB-TGcMV-Neo) were seeded in an eight-well chamber slide (BD Falcon, East Ruherford, New Jersey, US.) and cultivated overnight in the presence or absence of doxycyclin (2 μg/ml final). Cells were treated with 50 nM of mitomycin C for 2 days before fixation with 2% paraformaldehyde and permeabilization with 0.3% Triton X-100. Cells were blocked for 30 minutes with blocking/dilution buffer (10% goat serum, 0.5% NP-40, 0.5% saponin, 1× phosphate-buffered saline (PBS)) and incubated with the anti-phospho-histone H2A.X (serine 139, clone JBW301; Millipore, Billerica, Massachusetts, US) overnight at 4°C. Cells were then incubated with an Alexa 546 anti-mouse antibody for 1 hour at room temperature, incubated with DAPI (0.8 μg/ml; Sigma) for 10 minutes and mounted using ProLong Gold (Invitrogen). Images were captured using a 40× dry objective. Nuclei and γ-H2A.X foci were quantified using CellProfiler.
RNA extraction and library preparation for gene expression profiling of the RHOXF2 overexpressing cell line
The HEK293_M2 + RHOXF2 stable cell line (PB-TGcMV-Neo) was seeded in 10-cm dishes, cultivated 24 hours in the presence or absence of doxycycline (2 μg/ml final) before being treated with 40 nM of mitomycin C or DMSO (control) during a 2-day period. Total RNA was extracted (RNeasy Mini Kit, QIAGEN), mRNA-focused libraries were generated from 1 μg of total RNA (Illumina TruSeq RNA Sample Preparation Kit V.2) and sequenced using an Illumina HiSeq2000 according to the Illumina protocol. Raw data have been deposited in the ArrayExpress repository under accession number E-MTAB-2497.
RNA-sequencing data analysis
Single-end reads, 51 nucleotides in length, were generated on the Illumina HiSeq. Sequence data were aligned to the UCSC hg19 reference genome using TopHat (v2.0.0), provided with a RefSeq GTF file and instructed to align only across known junctions. Differential expression, as fragments per kilobase of exon per million fragments mapped (FPKM) differences, were generated by analysis with cuffdiff (v1.1.0) using the indicated comparisons. Based on their differential expression, potential activated and repressed isoforms were identified and related enriched biological processes were assessed using GOrilla.
RHOXF2 quantitative RT-PCR
Total RNA from 20 different normal human tissues (each pool from 3 donors) were obtained from the FirstChoice Human Total RNA Survey Panel (Ambion, Austin, Texas, US). Total RNA (5 μg) from each tissue was used for oligo(dT)12-18-primed reverse transcription using the SuperScript II reverse transcriptase (Invitrogen). Quantitative PCR was performed using the LightCycler 480 SYBR Green I Master (Roche) and the LightCycler 480 System with the following primers:
QhRHOXF1_1F (5′-ACCGTGTTCTACTGCCTGAGTGTA) and QhRHOXF1_1R (5′-TTCATGCCGTTCTCGTGGTTCACA) on RHOXF1 exon 1; QhRHOXF1_2F (5′-TGGAGGAGCTGGAAAGTGTT) and QhRHOXF1_2.2R (5′-GGCCCTTTTATTCTTAAACC) spanning RHOXF1 exon 1 and exon 3; QhRHOXF2_1F (5′-CCGGACCAGTGTAGCCAGTA) and QhRHOXF2_1R (5′-TCTTTTTCTTCTCCGCCTTG) spanning RHOXF2 exon 1 and exon 2; QhRHOXF2_2F (5′-ATGGTGCTGTCGCTTACTGA) and QhRHOXF2_2R (5′-TCGAGGTCTCCTTCCCATAG) on RHOXF2 exon 2; QhRHOXF2_4F (5′-CAGCGGGATGAGAGATGATT) and QhRHOXF2_4R (5′-TTGGGGAATGTGAAAGAAGG) on RHOXF2 exon 3.
Housekeeping/reference genes were: QhCYCG_F (5′-CTTGTCAATGGCCAACAGAGG) and QhCYCG_R (5′-GCCCATCTAAATGAGGAGTTGGT) on CYCG (cyclophilin G); QhGUSB_F (5′-ACGCAGAAAATATGTGGTTGGA) and QhGUSB_R (5′-GCACTCTCGTCGGTGACTGTT) on GUSB (beta-glucuronidase); QhActbF1 (5′-GAAGTCCCTTGCCATCCTAAAAG) and QhActbR1 (5′-AGGACTGGGCCATTCTCCTTA) on ACTB (beta-actin); QhEef1a1F2 (5′-CTGAACCATCCAGGCCAAAT) and QhEef1a1R2 (5′-AGCCGTGTGGCAATCCA) on EEF1A1 (eukaryotic translation elongation factor 1 alpha 1); QhGAPDH_F1 (5′-CATGAGAAGTATGACAACAGCCT) and QhGAPDH_R1 (5′-AGTCCTTCCACGATACCAAAGT) on GAPDH (glyceraldehyde-3-phosphate dehydrogenase).
Formalin-fixed, paraffin-embedded sections of human normal testis were deparaffinized for 12 hours at 58°C and rehydrated. Heat-induced epitope retrieval was in sodium citrate buffer pH 6.0 using a pressure cooker. All antibodies were diluted in Tris-buffered Saline (TBS) with 1% BSA. Sections were blocked in 10% goat serum (Gibco) for 2 hours at room temperature. The anti-RHOXF2 primary antibody (Sigma, HPA003314) was used at a 1:150 dilution and incubated overnight at 4°C in a humid chamber. Endogenous peroxidase activity was suppressed by incubating the section in 0.3% H2O2 for 15 minutes. Horse radish peroxidase (HRP)-conjugated secondary antibody (anti-rabbit IgG-HRP; GE Healthcare) was used at a 1:250 dilution and incubated for 1 hour at room temperature. 3,3′-Diaminobenzidine (Abcam) was used as chromogen. Sections were counterstained with Mayer’s hematoxylin (Sigma) before dehydration and mounting.
For the detection of endogenous RHOXF2, total cell lysates were prepared from 200,000 cells. Briefly, samples were lysed in 2× SB loading buffer (60 mM Tris-HCl pH6.8, 3% SDS, 10% glycerol, 0.05% bromophenol blue, 10% beta-mercaptoethanol) and subjected to SDS-PAGE (12%). The RHOXF2-specific mouse polyclonal antibody (Abcam, ab67811) and the tubulin specific rat antibody (Abcam, ab6160) were used at a dilution of 1:800 and 1:40,000 respectively. Immunological complexes were visualized by enhanced chemiluminescence using HRP-conjugated anti mouse or anti rat IgG.
Drug dose-response curves
K562 cells (100,000) were cultured in 600 μl of media in a 48-well plate for 24 hours before the addition of diluted drug. Two days after drug treatment, the percentage of drug inhibition was assessed by counting K562 using a Coulter particle count and cell analyser and compared to cells grown in the presence of DMSO only.
For the adherent HPAC, SW1353, and NCI-H1299 cell lines, 5,000 cells were seeded in 300 μl of media in 96-well plate and incubated for 24 hours prior to drug addition. Following two days of drug treatment, the percentage of drug inhibition was estimated by SRB assay and compared to cells grown in DMSO as control.
Development of the assay: DHFR and methotrexate
We implemented our approach on a larger scale by cloning, en masse, a pool of 376 random hORFs (including DHFR) into the lentiviral vector. Corresponding lentivirus were used to generate a population of HEK293_M2 cells with the stably integrated hORFs. Gene expression was induced as before by addition of doxycycline and cells were exposed to a lethal concentration of methotrexate. Ten days after treatment, surviving cells were harvested, genomic DNA was extracted, and hORFs were PCR amplified and hybridized on an Affymetrix GeneChip Human Gene 1.0 ST array to determine which hORFs conferred resistance to methotrexate. DHFR was found as the predominant gene in this pool of 376 hORFs whose overexpression provided resistance to methotrexate, validating our strategy at the scale of several hundred pooled clones (Additional file4).
The same approach was next applied to the complete collection of 12,212 hORFs. Here, we divided the collection into 34 minipools of 376 hORFs each, transferred them en masse into the pLD-T-IRES-Venus-WPRE-STOP lentivirus vector and then pooled 34 minipools together to produce sufficient quantities of lentivirus for the infection of HEK293_M2 cells at a multiplicity of infection of 0.3 (Figure 1). Because induction of the Venus marker is dependent on hORF expression, addition of doxycycline for 24 hours generated sufficient fluorescence to sort those cells with stable lentivirus integration (Additional file5). The presence and relative abundance of the virally integrated hORFs were assessed using next-generation sequencing. Based on the identification of at least 95% of the sequence with a minimum of 10× coverage, 87% of the mappable hORFs were considered as present in the HEK293_M2 cells. Further analysis of the sequencing results demonstrated a size-dependent relative abundance of the hORFs, with smaller clones being more represented in the population (Additional file6). An aliquot of 500,000 HEK293_M2 cells (an average of 41 cells per hORF) was seeded in a T175 flask, and grown for 15 hours before the induction of hORF expression by addition of doxycycline. After 18 hours of doxycycline induction, cells formed micro-colonies, which were cultured in the presence of methotrexate for 2 weeks. hORFs derived from the surviving cells were identified by PCR amplification and microarray hybridization. As we found for the minipool experiment, DHFR was the top gene conferring resistance to methotrexate when overexpressed (Figure 2c). These results encouraged us to test the same strategy for drugs with diverse modes of action.
Screens with additional drugs
Characteristics of screened compounds
Dermatologic agent (phase I)
Antineoplastic agent (phase II)
Antiobesity preparations, excluding diet products
Cytotoxic antibiotics and related substances
Other antineoplastic agents
Cytotoxic antibiotics and related substances
Experimental Hsp90 inhibitor
Antineoplastic agent (phases I and II)
Number of validated/tested candidate genes
Inner membrane protein, mitochondrial (mitofilin)
Mitogen-activated protein kinase kinase 1 interacting protein 1
Monocyte to macrophage differentiation-associated
Protein tyrosine phosphatase, non-receptor type 2
Rhox homeobox family, member 2
Chromosome 20 open reading frame 54
Coiled-coil domain containing 45
E74-like factor 5 (ets domain transcription factor)
Protein tyrosine phosphatase, non-receptor type 2
Rhox homeobox family, member 2
Transmembrane protein 150
Ubiquitin specific peptidase like 1
Zinc finger protein 64 homolog (mouse)
Eukaryotic translation initiation factor 4B
Screens for the antifolates methotrexate and aminopterin
For the antifolate screens, the single confirmed ORF was DHFR. This result highlights the fact that when a single hORF confers a strong growth advantage, it can dominate the other clones in the population and potentially mask the detection of other potential protective genes. One way to address this and potentially improve the dynamic range for this drug class would be to repeat these screens with a pool of hORFs lacking DHFR.
Among the potential expected hits present in our collection was the enzyme gamma-glutamyl hydrolase (GGH), which catalyzes the removal of polyglutamates from methotrexate. Therefore, overexpression of GGH would be expected to decrease the intracellular presence of methotrexate-polyglutamate, thereby increasing the efflux of methotrexate out of the cell and decreasing cytotoxicity. However, GGH overexpression was demonstrated to be insufficient to produce resistance to the drug in a human fibrosarcoma (HT-1080) and a human breast carcinoma (MCF-7) cell line.
DHFRL1 (dihydrofolate reductase like-1) consistently appeared as a top hit for screens with both methotrexate and aminopterin but failed to confer resistance when individually expressed in HEK293_M2 cells. Based on the strong similarity between DHFR and DHFRL1 (a total of 14 non-synonymous changes), we hypothesized that the observed unspecific signal was probably due to cross-hybridization. Interestingly, a similar problem occurred when the same strategy using the same hORFeome library was tested in the model organism S. cerevisiae. Both DHFR and DHFRL1 were enriched in the screens against methotrexate (and their expression confirmed at the protein level) but only DHFR conferred methotrexate resistance (data not shown).
Orlistat is a lipase inhibitor that reportedly targets LPL (lipoprotein lipase), PNLIP (pancreatic lipase) and FASN (fatty acid synthase). Although LPL is absent from our collection, PNLIP and FASN are present. However, genes for the latter two did not confer overexpression resistance.
Our screen identified the gene coding for monoglyceride lipase (MGLL) as a resistance hit. MGLL hydrolyzes intracellular triglyceride stores in adipocytes and other cells to fatty acids and glycerol. The enzyme might also be involved in the hydrolysis of monoglycerides. In vitro experiments suggest that orlistat can inhibit the MGLL-like activity in rat-cerebellar membranes and rat cerebellar homogenate, rendering the human MGLL a potential direct target for orlistat.
Another validated hit, PON3, is also involved in lipid metabolism. The gene encodes paraoxonase 3, which is secreted into the bloodstream and associates with high-density lipoprotein. The protein can also rapidly hydrolyze lactones and inhibit the oxidation of low-density lipoprotein.
Although the exact mechanism of action of bleomycin is unknown, available evidence indicates that its main mode of action is via the inhibition of DNA synthesis, with additional evidence for its inhibition of RNA and, to a lesser extent, protein synthesis.
We show that overexpression of IMMT (inner membrane protein, mitochondrial) partially rescues bleomycin toxicity. Bleomycin treatment induces damage of both nuclear and mitochondrial DNA. A recent study demonstrated that the mitochondrial localization of PARP-1 requires interaction with IMMT and that it is involved in the maintenance of mitochondrial DNA integrity. The authors also suggest that IMMT overexpression results in an increase of PARP-1 in the extracellular compartment. We therefore speculate that IMMT overexpression recruits sufficient PARP-1 into the mitochondrion to maintain mitochondrial DNA integrity and promote cell survival.
Cisplatin and mitomycin C screens
RHOXF2 and PTPN2 were found to suppress the toxicity of cisplatin and mitomycin C. Together with the antifolates, these DNA damaging agents were the only screens yielding identical hits. Also, as for methotrexate and aminopterin, both these alkylating agents have the same mechanism of action.
Another validated hit for mitomycin C, ZFP64 (zinc finger protein 64 homolog), is noteworthy because it has been shown to be potentially phosphorylated by ATM and ATR, the major signal transducing kinases of the DNA damage response in response to ionizing radiation. This observation is consistent with a potential role for ZFP64 in the repair of DNA damage generated by both alkylating agents and ionizing radiation.
RHOXF2 overexpression confers resistance to several DNA damaging agents
We next asked if RHOXF2′s protective effect was restricted to cisplatin and mitomycin C or if this gene, when overexpressed, could also suppress the toxicity of DNA damaging agents with different mechanisms of action. Accordingly, we assessed the cytotoxicity of antimetabolites (methotrexate, 5-fluorouracil), a radiomimetic (bleomycin), replication inhibitors (camptothecin, doxorubicin, etoposide), mitaplatin, 4-nitroquinoline-1-oxide and several experimental platinum anticancer agents (PT-ACRAMTU(EN), PT-ACRAMTU(PN), PT-ATUCA, PT-AMIDIN) on HEK293-M2 cells expressing or not the transcription factor (Figure 3c; Additional file8). With the exception of the topoisomerase I inhibitor camptothecin, all drugs whose toxicity was suppressed by RHOXF2 overexpression directly bind to DNA to exert their effects.
These observations motivated us to ask (i) if and how RHOXF2, an uncharacterized transcription factor, modulates gene expression in HEK293_M2 cells and (ii) if the nature of the regulated genes can illuminate the protective response to DNA damaging agents. We cultured HEK293_M2 cells with or without overexpression of RHOXF2 in the presence of 40 nM of mitomycin C or DMSO (control) and profiled their transcriptional response by RNA-seq. The sequencing data confirmed the RHOXF2 transcript was absent in the uninduced HEK293_M2 samples and verified its strong induction in the presence of doxycycline regardless of the drug treatment (Additional file9). Neither RHOXF1 nor RHOXF2B expression was affected by RHOXF2 overexpression (Additional file9). Interestingly, and despite their limited number, Gene Ontology analysis of the genes whose expression was induced during RHOXF2 overexpression showed that the most significantly enriched pathways are related to DNA damage processes (P < 7 × 10-5 and P < 3 × 10-4), stress response (P < 3 × 10-4), apoptosis (P < 4 × 10-4) and cell cycle (P < 7 × 10-4) (Figure 3d; Additional file9). As anticipated, in the absence of RHOXF2 and the presence of mitomycin C, we found a strong enrichment for genes related to the response to stress (P < 3 × 10-5) and regulation of protein metabolic process (P < 9.5 × 10-5). Finally, in the presence of RHOXF2, genes involved in the regulation of signal transduction (P < 4 × 10-4) were induced during drug treatment. Therefore, we suggest that when treated with DNA damaging agents, the protective effect conferred by overexpression of RHOXF2 is due, in part, to its activation of several key genes involved in the response to DNA damage.
Finally, we asked if the RHOXF2-dependent protective effect observed in HEK293_M2 cells was restricted to one cell type. We then tested if RHOXF2 overexpression could confer resistance to the rtTA-expressing human breast adenocarcinoma cell line MCF7_M2 and the adenocarcinomic human alveolar basal epithelial cell line A549_M2. As quantified by both a two-fold shift in IC50 (44.91 μM to 88.17 μM) and by the 21% decrease in toxicity at 3 μM, RHOXF2 overexpression protected MCF7_M2 cells against cisplatin-mediated toxicity (Additional file10a). In the absence of RHOXF2, addition of doxycycline resulted in a modest increase in IC50 (41.06 μM to 58.59 μM; empty vector). The effect of the vector alone, while detectable, is less than that observed in the presence of RHOXF2. In contrast, although RHOXF2 overexpression in A549_M2 cells provided some apparent resistance (shift in IC50 from 448 nM to 1,099 nM), the measured effect was, in fact, largely due to doxycycline (shift in IC50 from 718 nM to 1,454 nM with empty vector) (Additional file10b). In conclusion, our experiments suggest that the protective effect of RHOXF2 is not confined to a single genetic background and expression of the transcription factor is sufficient to protect both HEK293 and MCF7 cells against the DNA damaging agent cisplatin. As with all large-scale screens, the biological significance of this observation should be further confirmed, for example, by interrogating a larger panel of cell lines.
Based on its viability and transcriptional phenotypes in HEK293_M2 cells, we further investigated the potential role for RHOXF2 in neoplastic cells.
In normal tissue, RHOXF2 expression is confined to testis but RHOXF2 is expressed in several cancer cell lines
Computational analysis predicts RHOXF2 to be a testicular cancer candidate gene and, in fact, the transcript was detected in several types of human testicular cancers. Furthermore, in the present study, we show that RHOXF2 overexpression conferred resistance to cisplatin (Figure 3a) and that RHOXF2 is found in several other cancer cell lines, the non-small cell lung carcinoma NCI-H1299, the blast phase chronic myelogenous leukemia K562, the multiple myeloma U266B1, the thyroid gland medullary carcinoma TT (CRL-1803) and the chondrosarcoma SW1353 (Figure 4c; Additional file11).
In the chronic myelogenous leukemia K562 cell line, RHOXF2 depletion modulates the response to antineoplastic agents
Our study showed that RHOXF2 overexpression confers resistance to cisplatin and mitomycin C in HEK293_M2 cells, in which the protein is not endogenously expressed (Figure 4c; Additional file9). We first generalized this observation by positing that overexpression or ectopic expression of certain genes greatly reduces some detrimental drug effects on the growth of particular cancer cell lines. We then reasoned that if this is truly drug-gene specific, the opposite mechanism should experimentally verify it. More specifically, depletion of the endogenously expressed protein should confer sensitivity to the cell line during drug treatment. The most direct translation of this gene-drug specificity could then be exploited to understand drug resistance in cancer cell lines. Therefore, we asked if sensitivity to DNA damaging agents was observed when RHOXF2 is knocked down in cancer cells where the protein is normally expressed. We generated several stable cell lines to deplete RHOXF2 (NCI-H1299, K562, SW1353, and the human pancreatic HPAC as negative control) with five shRNAs that target the mRNA. Two shRNAs against GFP were used as a negative control. Two shRNAs showed a strong knock-down of the protein in all three RHOXF2-positive cell lines (Figure 4d; Additional file12e,h).
K562 showed a significant increase in sensitivity for both cisplatin (IC50: 3.231 μM to 2.052 μM) and mitomycin C (IC50: 377.8 nM to 250.4 nM) when the transcription factor was depleted, demonstrating the potential role of RHOXF2 in response to DNA damaging agents and particularly in response to DNA crosslinking agents (Figure 4e,f). Independent K562 transduction resulted in identical sensitivity to the drugs, showing the phenotype is independent of the lentiviral integration site. We ruled out the possibility that the calculated drug sensitivity simply reflected a growth difference among cell lines by demonstrating that the RHOXF2_shRNA and GFP_shRNA control cell lines divided at the same rate (Figure 4e,f, insets). Finally, we assessed the robustness of our knock-downs in the presence of the drugs by confirming that the modulation of RHOXF2 protein level is maintained even when the K562 cells were exposed to cisplatin or mitomycin C at IC70 (Additional file13a,b).
As expected, no change was measured for the RHOXF2-negative HPAC cells (Additional file12a,b). Interestingly, no increase in sensitivity to cisplatin and mitomycin C was observed for NCI-H1299 and SW1353 knock-down cell lines (Additional file12c,d,f,g).
The modest increase in sensitivity of K562 cells to DNA damaging agents can be explained by different factors. First, transcription factors are generally present in limited quantities in cells[38–40]. Here, although clearly depleted in the presence of the shRNAs, RHOXF2 is still present at a relatively high level in K562 cells (Figure 4d). Therefore, we speculate that despite the efficient knock-down, the limited shift in IC50 is in part due to the presence of enough protein to sustain a certain level of resistance during drug treatment. In order to test our first hypothesis, a complete knock-out of the transcription factor would be necessary (for example, using CRISPR technology). Secondly, a recent study in mouse embryonic stem cells reports that knock-down of a few transcription factors is associated with substantial transcriptome change. These data not only indicate the robustness of the pluripotency gene network but also suggest that perhaps RHOXF2 knock-down hardly perturbs the transcription factor network in the K562 cellular context.
We have reduced to practice an elegant method to simultaneously evaluate thousands of genes for their capacity to provide a survival advantage in the presence of a toxic dose of drug using cells expressing the entire hORF collection and a microarray or sequencing-based readout.
The availability of hORF collections represents great progress compared to traditional approaches (random mutagenesis, cDNA libraries, and so on) because it permits the systematic interrogation of a gene's biological function in a targeted manner. Such gain-of-function approaches would be bolstered when combined with the reciprocal loss-of-function approach. To date, several studies have already explored the availability of such collections. Most of the studies tested the hORFs function in an array (well-based) format[43–45] that simplifies the readout and analysis of the results. Here, we tested the function of thousands of ORFs simultaneously in a competitive manner. Our pool approach greatly reduced the tedious and expensive upstream steps of cloning, transduction and testing of the genes. The primary challenges to the pool approach are (i) potential loss of ORFs during cloning and infection (leading to representation variability), (ii) the potential masking effect of overrepresented clones in a population during selection (for example, DHFR in the case of methotrexate) and (iii) the necessity for dedicated analysis of the final results (hits prioritization). In addition to future enhancements of the human ORFeome collection, other promising gain-of-function approaches to assess drug mode of action are the development of functional variomics technology and the exploration of human isoform space in human cell lines.
In developing the assay, we first determined the optimal period of time between the beginning of the gene expression and the addition of the chemical compound. To address that question, we quantified the number of Venus-positive cells 12 hours, 24 hours and 48 hours after induction of expression of the hORF collection. Because the human gene and the Venus marker are co-transcribed, quantifying the level of Venus expression by flow cytometry serves as a proxy for the level of human gene expression. Because the number of Venus-positive cells peaks at 24 hours after doxycycline addition (Additional file5), we postulated that 18 hours post-induction would be sufficient time after which to add drug to test the human gene for any protective effects.
We then addressed the sensitivity of the technique by improving the ratio between the number of seeded cells, the growth area and the time of drug addition. Seeding 500,000 cells in a T175 flask 33 hours before drug addition provided enough time and density for isolated cells to form separate micro-colonies that should overexpress only one or a few genes. This 'clonogenic assay approach' was more sensitive to lower drug concentrations and provided a better selective environment for the surviving cells, presumably by minimizing signaling from a large population of dying cells that would arise at higher plating densities. Similarly, adding fresh medium and drug every 48 hours for the duration of the screen improved the selection process by washing away the dead and dying adherent cells that could have interfered with the subsequent screen readout.
During the course of the screen, the cell population followed a stereotypical pattern of growth, starting with a dramatic decrease in cell number and disappearance of the majority of cells, followed by a stabilization of the population and concluding with the development of medium to large size colonies originating from the initial surviving cells. Daily visual inspection of the cell population was crucial to determine when to stop the drug selection and harvest the cells, typically after 20 to 30 doubling times.
Several refinements may enhance this assay. First, improving the coverage and expression quality of the human genes using a more complete collection of fully sequenced ORFs as such collections become available. Moreover, as illustrated in Additional file6, relative abundance of the virally integrated hORFs is size-dependent. Therefore, to correct for that variation, preparation of the normalized amount of plasmid DNA for lentivirus production could be done as follows: i) reduce the pool size (less than 376 hORFs per pool as in our study); ii) hORFs should be grouped by gene length. One can also match the particular drug under interrogation to disease-appropriate cell types and environmental conditions. Whenever possible, experiments should be performed with an initial high number of cells to guarantee the full representation of the collection during drug target screening. Although DNA damaging resistance was observed in the RHOXF2-overexpressing HEK293_M2 cells, it would be informative to further study if that observation is dependent on RHOXF2 expression status. Finally, next-generation sequencing will facilitate the throughput of hit detection and reduce the number of false positives arising from microarray cross-hybridization.
Our study revealed several genes whose functions relate to the action of the drugs tested. RHOXF2 is of a particular interest because of (i) its strong rescue phenotype in HEK293_M2 cells, (ii) its capacity to suppress both cisplatin and mitomycin C toxicity, and (iii) the dearth of functional information on this transcription factor. Homeobox genes encode transcription factors that play a central role during embryogenesis. RHOXF2 belongs to the Rhox family of genes, which are expressed not only during embryogenesis but also in adult reproductive tissue. In normal tissue, RHOXF2 expression is likely restricted to testis but its function remains largely unknown[35, 48]. Here, we demonstrated that (i) the presence of RHOXF2 was sufficient to lower both cisplatin and mitomycin C toxicity in HEK293_M2 cells, and (ii) the transcription factor might exert its effect via the activation of genes involved in response to DNA damage. Furthermore, RHOXF2 protein expression was detected in several cancer cell lines. In one of these lines (chronic myelogenous leukemia/K562) knocking down the protein increased the toxicity for both DNA damaging agents. Interestingly, although no growth disadvantage was observed in K562 cells partially depleted for RHOXF2 (Figure 4e,f), its depletion via shRNA increased the sensitivity of this leukemia cell line to several pharmacological compounds. We therefore speculate that although not a driver of tumorigenesis, RHOXF2 could confer on the neoplastic cells a selective advantage during drug treatment.
We demonstrate the usefulness of our approach in order to understand drug mechanisms of action by combining genetic perturbation with drug treatment. Here we recapitulated examples of drug resistance and uncovered several unanticipated drug-target interactions by virtue of the suppression of chemical toxicity observed upon gene overexpression. We also provide a step-by-step protocol for performing such screens and for analyzing the data (Additional file1). Moreover, our innovative overexpression approach could also be applied to identify genes whose up-regulation could be toxic or enhance cell proliferation in different cell types. For example, the assay could be useful to uncover tumor-specific activity in the context of genome engineering as well as for the discovery of cell-type-specific oncogenes.
bovine serum albumin
Dulbecco's modified Eagle's medium
horse radish peroxidase
open reading frame
polymerase chain reaction
short hairpin RNA
CN and GG were supported by a grant from the Canadian Cancer Society (#20830) and the NHGRI.
- Jarosz DF, Taipale M, Lindquist S: Protein homeostasis and the phenotypic manifestation of genetic diversity: principles and mechanisms. Annu Rev Genet. 2010, 44: 189-216. 10.1146/annurev.genet.40.110405.090412View ArticlePubMedGoogle Scholar
- Nobeli I, Favia AD, Thornton JM: Protein promiscuity and its implications for biotechnology. Nat Biotechnol. 2009, 27: 157-167. 10.1038/nbt1519View ArticlePubMedGoogle Scholar
- Hopkins AL: Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008, 4: 682-690. 10.1038/nchembio.118View ArticlePubMedGoogle Scholar
- Strebhardt K, Ullrich A: Paul Ehrlich’s magic bullet concept: 100 years of progress. Nat Rev Cancer. 2008, 8: 473-480. 10.1038/nrc2394View ArticlePubMedGoogle Scholar
- Agarwal R, Kaye SB: Ovarian cancer: strategies for overcoming resistance to chemotherapy. Nat Rev Cancer. 2003, 3: 502-516. 10.1038/nrc1123View ArticlePubMedGoogle Scholar
- Seruga B, Ocana A, Tannock IF: Drug resistance in metastatic castration-resistant prostate cancer. Nat Rev Clin Oncol. 2011, 8: 12-23. 10.1038/nrclinonc.2010.136View ArticlePubMedGoogle Scholar
- Curt GA, Carney DN, Cowan KH, Jolivet J, Bailey BD, Drake JC, Chien Song KS, Minna JD, Chabner BA: Unstable methotrexate resistance in human small-cell carcinoma associated with double minute chromosomes. N Engl J Med. 1983, 308: 199-202. 10.1056/NEJM198301273080406View ArticlePubMedGoogle Scholar
- Cohen SN, Chang AC, Boyer HW, Helling RB: Construction of biologically functional bacterial plasmids in vitro. Proc Natl Acad Sci U S A. 1973, 70: 3240-3244. 10.1073/pnas.70.11.3240PubMed CentralView ArticlePubMedGoogle Scholar
- Chang AC, Nunberg JH, Kaufman RJ, Erlich HA, Schimke RT, Cohen SN: Phenotypic expression in E. coli of a DNA sequence coding for mouse dihydrofolate reductase. Nature. 1978, 275: 617-624. 10.1038/275617a0View ArticlePubMedGoogle Scholar
- Rine J, Hansen W, Hardeman E, Davis RW: Targeted selection of recombinant clones through gene dosage effects. Proc Natl Acad Sci U S A. 1983, 80: 6750-6754. 10.1073/pnas.80.22.6750PubMed CentralView ArticlePubMedGoogle Scholar
- Rual JF, Hirozane-Kishikawa T, Hao T, Bertin N, Li S, Dricot A, Li N, Rosenberg J, Lamesch P, Vidalain PO, Clingingsmith TR, Hartley JL, Esposito D, Cheo D, Moore T, Simmons B, Sequerra R, Bosak S, Doucette-Stamm L, Le Peuch C, Vandenhaute J, Cusick ME, Albala JS, Hill DE, Vidal M: Human ORFeome version 1.1: a platform for reverse proteomics. Genome Res. 2004, 14: 2128-2135. 10.1101/gr.2973604PubMed CentralView ArticlePubMedGoogle Scholar
- Lamesch P, Li N, Milstein S, Fan C, Hao T, Szabo G, Hu Z, Venkatesan K, Bethel G, Martin P, Rogers J, Lawlor S, McLaren S, Dricot A, Borick H, Cusick ME, Vandenhaute J, Dunham I, Hill DE, Vidal M: hORFeome v3.1: a resource of human open reading frames representing over 10, 000 human genes. Genomics. 2007, 89: 307-315. 10.1016/j.ygeno.2006.11.012View ArticlePubMedGoogle Scholar
- Shin KJ, Wall EA, Zavzavadjian JR, Santat LA, Liu J, Hwang JI, Rebres R, Roach T, Seaman W, Simon MI, Fraser ID: A single lentiviral vector platform for microRNA-based conditional RNA interference and coordinated transgene expression. Proc Natl Acad Sci U S A. 2006, 103: 13759-13764. 10.1073/pnas.0606179103PubMed CentralView ArticlePubMedGoogle Scholar
- Yang X, Boehm JS, Salehi-Ashtiani K, Hao T, Shen Y, Lubonja R, Thomas SR, Alkan O, Bhimdi T, Green TM, Johannessen CM, Silver SJ, Nguyen C, Murray RR, Hieronymus H, Balcha D, Fan C, Lin C, Ghamsari L, Vidal M, Hahn WC, Hill DE, Root DE: A public genome-scale lentiviral expression library of human ORFs. Nat Methods. 2011, 8: 659-661. 10.1038/nmeth.1638PubMed CentralView ArticlePubMedGoogle Scholar
- ArrayExpress. [http://www.ebi.ac.uk/arrayexpress/],
- Vichai V, Kirtikara K: Sulforhodamine B colorimetric assay for cytotoxicity screening. Nat Protoc. 2006, 1: 1112-1116. 10.1038/nprot.2006.179View ArticlePubMedGoogle Scholar
- Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang IH, Friman O, Guertin DA, Chang JH, Lindquist RA, Moffat J, Golland P, Sabatini DM: Cell Profiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 2006, 7: R100- 10.1186/gb-2006-7-10-r100PubMed CentralView ArticlePubMedGoogle Scholar
- Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009, 10: 48- 10.1186/1471-2105-10-48PubMed CentralView ArticlePubMedGoogle Scholar
- Abcam Immunohistochemistry Procedure. [http://docs.abcam.com/pdf/protocols/Short_IHC-P_protocol.pdf],
- Vander Heiden MG: Targeting cancer metabolism: a therapeutic window opens. Nat Rev Drug Discov. 2011, 10: 671-684. 10.1038/nrd3504View ArticlePubMedGoogle Scholar
- Farber S, Diamond LK, Mercer RD, Sylvester RF, Wolff JA: Temporary remissions in acute leukemia in children produced by folic acid antagonist, 4-aminopteroyl-glutamic acid (aminopterin). N Engl J Med. 1948, 238: 787-793. 10.1056/NEJM194806032382301View ArticlePubMedGoogle Scholar
- Shi Y, Burn P: Lipid metabolic enzymes: emerging drug targets for the treatment of obesity. Nat Rev Drug Discov. 2004, 3: 695-710. 10.1038/nrd1469View ArticlePubMedGoogle Scholar
- Hogan S, Fleury A, Hadvary P, Lengsfeld H, Meier MK, Triscari J, Sullivan AC: Studies on the antiobesity activity of tetrahydrolipstatin, a potent and selective inhibitor of pancreatic lipase. Int J Obes. 1987, 11: 35-42.PubMedGoogle Scholar
- Rose JL, Reeves KC, Likhotvorik RI, Hoyt DG: Base excision repair proteins are required for integrin-mediated suppression of bleomycin-induced DNA breakage in murine lung endothelial cells. J Pharmacol Exp Ther. 2007, 321: 318-326. 10.1124/jpet.106.113498View ArticlePubMedGoogle Scholar
- Kelland LR: Emerging drugs for ovarian cancer. Expert Opin Emerg Drugs. 2005, 10: 413-424. 10.1517/14728126.96.36.1993View ArticlePubMedGoogle Scholar
- Iyer VN, Szybalski W: A molecular mechanism of mitomycin action: linking of complementary DNA strands. Proc Natl Acad Sci U S A. 1963, 50: 355-362. 10.1073/pnas.50.2.355PubMed CentralView ArticlePubMedGoogle Scholar
- Egorin MJ, Lagattuta TF, Hamburger DR, Covey JM, White KD, Musser SM, Eiseman JL: Pharmacokinetics, tissue distribution, and metabolism of 17-(dimethylaminoethylamino)-17-demethoxygeldanamycin (NSC 707545) in CD2F1 mice and Fischer 344 rats. Cancer Chemother Pharmacol. 2002, 49: 7-19. 10.1007/s00280-001-0380-8View ArticlePubMedGoogle Scholar
- Cole PD, Kamen BA, Gorlick R, Banerjee D, Smith AK, Magill E, Bertino JR: Effects of overexpression of gamma-Glutamyl hydrolase on methotrexate metabolism and resistance. Cancer Res. 2001, 61: 4599-4604.PubMedGoogle Scholar
- Karlsson M, Reue K, Xia YR, Lusis AJ, Langin D, Tornqvist H, Holm C: Exon-intron organization and chromosomal localization of the mouse monoglyceride lipase gene. Gene. 2001, 272: 11-18. 10.1016/S0378-1119(01)00559-5View ArticlePubMedGoogle Scholar
- Saario SM, Laitinen JT: Monoglyceride lipase as an enzyme hydrolyzing 2-arachidonoylglycerol. Chem Biodivers. 2007, 4: 1903-1913. 10.1002/cbdv.200790158View ArticlePubMedGoogle Scholar
- Chen J, Stubbe J: Bleomycins: towards better therapeutics. Nat Rev Cancer. 2005, 5: 102-112. 10.1038/nrc1547View ArticlePubMedGoogle Scholar
- Shen CC, Wertelecki W, Driggers WJ, LeDoux SP, Wilson GL: Repair of mitochondrial DNA damage induced by bleomycin in human cells. Mutat Res. 1995, 337: 19-23. 10.1016/0921-8777(95)00008-8View ArticlePubMedGoogle Scholar
- Rossi MN, Carbone M, Mostocotto C, Mancone C, Tripodi M, Maione R, Amati P: Mitochondrial localization of PARP-1 requires interaction with mitofilin and is involved in the maintenance of mitochondrial DNA integrity. J Biol Chem. 2009, 284: 31616-31624. 10.1074/jbc.M109.025882PubMed CentralView ArticlePubMedGoogle Scholar
- Matsuoka S, Ballif BA, Smogorzewska A, McDonald ER, Hurov KE, Luo J, Bakalarski CE, Zhao Z, Solimini N, Lerenthal Y, Shiloh Y, Gygi SP, Elledge SJ: ATM and ATR substrate analysis reveals extensive protein networks responsive to DNA damage. Science. 2007, 316: 1160-1166. 10.1126/science.1140321View ArticlePubMedGoogle Scholar
- Wayne CM, MacLean JA, Cornwall G, Wilkinson MF: Two novel human X-linked homeobox genes, hPEPP1 and hPEPP2, selectively expressed in the testis. Gene. 2002, 301: 1-11. 10.1016/S0378-1119(02)01087-9View ArticlePubMedGoogle Scholar
- Hofmann O, Caballero OL, Stevenson BJ, Chen YT, Cohen T, Chua R, Maher CA, Panji S, Schaefer U, Kruger A, Lehvaslaiho M, Carninci P, Hayashizaki Y, Jongeneel CV, Simpson AJ, Old LJ, Hide W: Genome-wide analysis of cancer/testis gene expression. Proc Natl Acad Sci U S A. 2008, 105: 20422-20427. 10.1073/pnas.0810777105PubMed CentralView ArticlePubMedGoogle Scholar
- Korkola JE, Houldsworth J, Chadalavada RS, Olshen AB, Dobrzynski D, Reuter VE, Bosl GJ, Chaganti RS: Down-regulation of stem cell genes, including those in a 200-kb gene cluster at 12p13.31, is associated with in vivo differentiation of human male germ cell tumors. Cancer Res. 2006, 66: 820-827. 10.1158/0008-5472.CAN-05-2445View ArticlePubMedGoogle Scholar
- Vaquerizas JM, Kummerfeld SK, Teichmann SA, Luscombe NM: A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009, 10: 252-263. 10.1038/nrg2538View ArticlePubMedGoogle Scholar
- Ishihama Y, Schmidt T, Rappsilber J, Mann M, Hartl FU, Kerner MJ, Frishman D: Protein abundance profiling of the Escherichia coli cytosol. BMC Genomics. 2008, 9: 102- 10.1186/1471-2164-9-102PubMed CentralView ArticlePubMedGoogle Scholar
- Yang VW: Eukaryotic transcription factors: identification, characterization and functions. J Nutr. 1998, 128: 2045-2051.PubMedGoogle Scholar
- Nishiyama A, Sharov AA, Piao Y, Amano M, Amano T, Hoang HG, Binder BY, Tapnio R, Bassey U, Malinou JN, Correa-Cerro LS, Yu H, Xin L, Meyers E, Zalzman M, Nakatake Y, Stagg C, Sharova L, Qian Y, Dudekula D, Sheer S, Cadet JS, Hirata T, Yang HT, Goldberg I, Evans MK, Longo DL, Schlessinger D, Ko MS: Systematic repression of transcription factors reveals limited patterns of gene expression changes in ES cells. Sci Rep. 2013, 3: 1390-PubMed CentralView ArticlePubMedGoogle Scholar
- Hoon S, Smith AM, Wallace IM, Suresh S, Miranda M, Fung E, Fung E, Proctor M, Shokat KM, Zhang C, Davis RW, Giaever G, St Onge RP, Nislow C: An integrated platform of genomic assays reveals small-molecule bioactivities. Nat Chem Biol. 2008, 4: 498-506. 10.1038/nchembio.100View ArticlePubMedGoogle Scholar
- Skalamera D, Dahmer M, Purdon AS, Wilson BM, Ranall MV, Blumenthal A, Gabrielli B, Gonda TJ: Generation of a genome scale lentiviral vector library for EF1alpha promoter-driven expression of human ORFs and identification of human genes affecting viral titer. PLoS One. 2012, 7: e51733- 10.1371/journal.pone.0051733PubMed CentralView ArticlePubMedGoogle Scholar
- Skalamera D, Ranall MV, Wilson BM, Leo P, Purdon AS, Hyde C, Nourbakhsh E, Grimmond SM, Barry SC, Gabrielli B, Gonda TJ: A high-throughput platform for lentiviral overexpression screening of the human ORFeome. PLoS One. 2011, 6: e20057- 10.1371/journal.pone.0020057PubMed CentralView ArticlePubMedGoogle Scholar
- Reece-Hoyes JS, Barutcu AR, McCord RP, Jeong JS, Jiang L, MacWilliams A, Yang X, Salehi-Ashtiani K, Hill DE, Blackshaw S, Zhu H, Dekker J, Walhout AJ: Yeast one-hybrid assays for gene-centered human gene regulatory network mapping. Nat Methods. 2011, 8: 1050-1052. 10.1038/nmeth.1764PubMed CentralView ArticlePubMedGoogle Scholar
- Huang Z, Chen K, Zhang J, Li Y, Wang H, Cui D, Tang J, Liu Y, Shi X, Li W, Liu D, Chen R, Sucgang RS, Pan X: A functional variomics tool for discovering drug-resistance genes and drug targets. Cell Rep. 2013, 3: 577-585. 10.1016/j.celrep.2013.01.019PubMed CentralView ArticlePubMedGoogle Scholar
- Maclean JA, Chen MA, Wayne CM, Bruce SR, Rao M, Meistrich ML, Macleod C, Wilkinson MF: Rhox: a new homeobox gene cluster. Cell. 2005, 120: 369-382. 10.1016/j.cell.2004.12.022View ArticlePubMedGoogle Scholar
- MacLean JA, Wilkinson MF: The Rhox genes. Reproduction. 2010, 140: 195-213. 10.1530/REP-10-0100View ArticlePubMedGoogle Scholar
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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.