High-throughput RNA interference (RNAi) screening using a pool of lentiviral shRNAs can be a tool to detect therapeutically relevant synthetic lethal targets in malignancies. We provide a pooled shRNA screening approach to investigate the epigenetic effectors in acute myeloid leukemia (AML).
Understanding clinically relevant driver mechanisms of acquired chemo-resistance is crucial for elucidating ways to circumvent resistance and improve survival in patients with acute myeloid leukemia (AML). A small fraction of leukemic cells that survive chemotherapy have a poised epigenetic state to tolerate chemotherapeutic insult. Further exposure to chemotherapy allows these drug persister cells to attain a fixed epigenetic state, which leads to altered gene expression, resulting in the proliferation of these drug-resistant populations and eventually relapse or refractory disease. Therefore, identifying epigenetic modulations that necessitate the survival of drug-resistant leukemic cells is critical. We detail a protocol to identify epigenetic modulators that mediate resistance to the nucleoside analog cytarabine (AraC) using pooled shRNA library screening in an acquired cytarabine-resistant AML cell line. The library consists of 5,485 shRNA constructs targeting 407 human epigenetic factors, which allows high-throughput epigenetic factor screening.
Therapeutic options in acute myeloid leukemia (AML) have remained unchanged for nearly the past five decades, with cytarabine (AraC) and anthracyclines as the cornerstone for treating the disease. One of the challenges to the success of AML therapy is the resistance of leukemic stem cells to chemotherapy, leading to disease relapse1,2. Epigenetic regulation plays a vital role in cancer pathogenesis and drug resistance, and several epigenetic factors have emerged as promising therapeutic targets3,4,5. Epigenetic regulatory mechanisms affect proliferation and survival under continuous exposure to chemotherapeutic drugs. Studies in non-hematological malignancies have reported that a small fraction of cells that overcome the drug effect undergo various epigenetics modifications, resulting in those cells' survival6,7. However, the role of epigenetic factors in mediating acquired resistance to cytarabine in AML has not been explored.
High-throughput screening is an approach to drug discovery that has gained global importance over time and has become a standard method in different aspects to identify potential targets in cellular mechanisms, for pathway profiling, and at the molecular level8,9. The synthetic lethality concept involves the interaction between two genes where the perturbation of either gene alone is viable but of both genes simultaneously results in the loss of viability10. Exploiting synthetic lethality in cancer treatment could help identify and mechanistically characterize robust synthetic lethal genetic interactions11. We have adopted a combinatorial approach of high-throughput shRNA screening with synthetic lethality to identify the epigenetic factors responsible for acquired cytarabine resistance in AML.
Acute leukemias driven by chromosomal translocation of the mixed-lineage leukemia gene (MLL or KMT2A) are known to be associated with poor survival in patients. The resulting chimeric products of MLL gene rearrangements, i.e., MLL fusion proteins (MLL-FPs), can transform hematopoietic stem/progenitor cells (HSPCs) into leukemic blasts with the involvement of multiple epigenetic factors. These epigenetic regulators constitute a complicated network that dictates the maintenance of the leukemia program and, therefore, could form potential therapeutic targets. In this context, we used the MV4-11 cell line (harboring the MLL fusion gene MLL-AF4 with the FLT3-ITD mutation; termed as MV4-11 P) to develop the acquired cytarabine resistant cell line, termed as MV4-11 AraC R. The cell line was exposed to increasing doses of cytarabine with intermittent recovery from the drug treatment, known as a drug holiday. The half-maximal inhibitory concentration (IC50) was assessed by in vitro cytotoxicity assay.
We used the pooled epigenetic shRNA Library (see Table of Materials) driven by the hEF1a promoter with a pZIP lentiviral backbone. This library comprises shRNAs targeting 407 epigenetic factors. Each factor has 5-24 shRNAs, with a total of 5,485 shRNAs, including five non-targeting control shRNAs. The modified "UltrmiR" miR-30 scaffold has been optimized for efficient primary shRNA biogenesis and expression12,13.
The outline of this experiment is illustrated in Figure 1A. The current protocol focuses on RNAi screening using the epigenetic factor shRNA library in the MV4-11 AraC R cell line (Figure 1B), a suspension cell line. This protocol can be used to screen any targeted library in any drug-resistant cell line of one's choice. It should be noted that the transduction protocol will be different for adherent cells.
Follow the Institutional Biosafety Committee (IBSC) guidelines and use the proper facility to handle lentivirus (BSL-2). Personnel should be appropriately trained in the handling and disposal of lentivirus. This protocol follows the biosafety guidelines of Christian Medical College, Vellore.
1. Selection of the most potent promoter to obtain persistent and prolonged expression of the shRNAs
NOTE: It is essential to perform a transduction experiment using lentiviral vectors with different promoters that express fluorescence proteins to identify the promoter that provides a stable and long-term expression of shRNAs in the cell line selected for the experiment. The most commonly used promoters for this purpose are hEF1α (human elongation factor 1α), hCMV (human cytomegalovirus), and SSFV (spleen focus-forming virus) promoters that express green fluorescence protein (GFP) (Figure 2A).
2. Preparation of pooled lentiviral human epigenetic factor shRNA library
3. Estimation of the transduction efficiency of lentiviruses
4. Transduction of pooled epigenetic shRNA library in the drug-resistant cell line
5. Enrichment of GFP positive cells
NOTE: Expand the transduced cells by culturing them at a density of 0.5 x 106 cells/mL for 5-7 days. These cells are a mixed population of transduced and untransduced, selected based on GFP by sorting, as mentioned in the next step.
6. Dropout screening to identify epigenetic factors mediating drug resistance
7. Amplification of the integrated shRNAs by PCR
8. Next-generation sequencing and data analysis
The overall screening workflow is depicted in Figure 1A. In vitro cytotoxicity of the MV4-11 AraC R (48 h) revealed the IC50 to cytarabine in the MV4-11 AraC R to be higher than the MV4-11 P (Figure 1B). This cell line was used in the study as the model for screening the epigenetic factors responsible for cytarabine resistance.
Figure 2A shows the linearized pZIP vector maps with three different promoters: hEF1α (human elongation factor 1α), hCMV (human cytomegalovirus), and SSFV (spleen focus-forming virus), with the scrambled shRNA within the UltramiR sequence. Figure 2B illustrates the pZIP vector map used in the library experiment and the shRNA targeting the epigenetic factor with Zsgreen and puromycin as the selectable marker driven by the hEF1a promoter. These vectors were used in the selection of the most potent promoter experiment.
The selection of the most potent promoter to get a consistent and prolonged expression in the target cells is illustrated in Figure 3. The quantification of GFP measured by MFI showed more than 90% transduction efficiency at the end of 72 h in MV4-11 AraC R for all three promoters. The histogram for GFP expression shows a heterogeneous population for hCMV but a homogeneous one in the case of hEF1a and SFFV at day 3 of transduction (Figure 3A). The bar graph shows the GFP expression driven by three different promoters (hEF1α, hCMV, and SFFV) at day 3 of transduction in MV4-11 AraC R (Figure 3B). SFFV-driven GFP showed a reduction in MFI on day 5 of transduction (Figure 3C). No decrease in MFI was observed in the hEF1a-driven GFP transduced MV4-11 parental line post sorting. Thus, the hEF1a promoter was identified as a suitable promoter for the cell line (Figure 3D).
Figure 4A illustrates the transfection efficiency of the pooled shRNA library plasmids in 293T cells post 48 h of transfection. The GFP expression was bright, indicating good transfection efficiency. Post virus collection, the transduction efficiency of the prepared pooled virus was checked in 293T cells in three different concentrations (2 µL, 4 µL, and 8 µL). We observed that even 2 µL virus resulted in the same efficiency as that of 8 µL virus, thereby confirming the high viral titer (Figure 4B).
Figure 5 illustrates the pooled library transduction in the MV4-11 AraC R cell line and the determination of the lentiviral titer. The pooled shRNA library lentivirus was diluted 100x and then added to the MV4-11 AraC R cell line in different volumes (1 µL, 1.5 µL, 2 µL, and 2.5 µL). The efficiency was checked at the end of 72 h. The transduction efficiency was 30% for 2-2.5 µL of the virus, confirming one shRNA integration per cell (Figure 5A). Transduction with 2.5μl of pooled library virus in 1.1 x 107 MV4-11 AraC R cells resulted in 30% transduction efficiency. GFP-positive cells were sorted, representing the pooled shRNA library (Figure 5B).
After the transduction of pooled shRNA epigenetic lentivirus in the MV4-11 AraC R cell line, the GFP-positive cells were sorted on day 5 or day 7 (Figure 6). These sorted cells were subjected to prolonged exposure to cytarabine for 9 days. The viability was checked on the 9th day, and the cells were collected for DNA. (Figure 6A). The bar graph shows the reduction in the proliferation rate of the transduced cells in the presence of cytarabine (Figure 6B).
The extracted DNA from the samples was subjected to two rounds of PCR, and the final gel eluted product represented the enriched shRNAs quantified by NGS. Figure 7A shows the binding regions of the primer used in the 1st and the 2nd round of PCR, where the 1st round primers bind to the flanking regions of the integrated shRNA, and the 2nd forward primer binds to the loop sequence. The reverse primer binds to a region of the amplified vector in the 1st round of PCR. Figure 7B and Figure 7C illustrate the band size of the PCR product at the end of 1st (397 bp) and the 2nd round PCR (399 bp), which was gel-eluted, purified, and submitted for NGS analysis.
After subjecting the DNA to a standard quality control procedure, the samples were processed for NGS analysis. We used CRISPRCloud2, a user-friendly, cloud-based analysis platform for identifying enriched and depleted shRNAs in the pooled shRNA screening. Figure 8 illustrates the representation of enriched or depleted shRNAs targeting epigenetic factors that could mediate cytarabine resistance in AML.
Figure 1: Outline and model of the study. (A) Illustration of the overview of the workflow. (B) MV4-11 parental and AraC-resistant cells treated with increasing cytarabine concentrations (0.1 µM to 1000 µM) and viability-assessed by MTT assay. Please click here to view a larger version of this figure.
Figure 2. Illustration of the linearized vectors. (A) The three vector maps with three different promoters, hEF1α (human elongation factor 1α), hCMV (human cytomegalovirus), and SSFV (spleen focus-forming virus), with the scrambled shRNA within the UltramiR sequence. (B) Illustration of the vector map used in the library experiment with hEF1α promoter and the shRNA targeting the epigenetic factor with Zsgreen and puromycin as the selectable marker. Please click here to view a larger version of this figure.
Figure 3: Selection of the most potent promoter to obtain persistent and prolonged expression of the shRNAs. (A) Representative flow plots for GFP quantified by flow cytometry at the end of 72 h for hEF1a, hCMV, and SFFV promoters. hCMV shows a heterogeneous peak, while hEF1a and SFFV show a single homogenous peak. (B) Promoter efficiency in the MV4-11 AraC R cell line. (C) SFFV-driven GFP shows silencing of GFP in prolonged culture. (D) hEF1a-driven GFP cells show sustained expression post sorting of the cells. Please click here to view a larger version of this figure.
Figure 4: Checking the efficiency of the prepared pooled shRNA lentivirus. (A) Fluorescence microscopy images show the transfection efficiency of the pooled shRNA library transfected in 293T at the end of 48 h using 10x magnification. (B) The transduction efficiency of the pooled virus was assessed in 293T cells with varying virus volumes (2µL, 4µL, and 8µL) using 10x magnification. Please click here to view a larger version of this figure.
Figure 5: Pooled library transduction in the target cell line and determination of lentiviral titer. (A) MV4-11 AraC R cell line transduced with varying volumes of virus (1µL, 1.5µL, 2µL, and 2.5µL). The efficiency was checked at the end of 72 h. (B) Transduction of pooled library virus in 1 x 107 MV4-11 AraC R cells to achieve 30% transduction efficiency, then sorting of GFP-positive cells. Please click here to view a larger version of this figure.
Figure 6: Dropout screening to identify epigenetic factors mediating drug resistance. (A) Schematic illustration of drug treatment for the enrichment of resistant cells.Library infected cells were subjected to prolonged cytarabine exposure for 9 days, followed by checking the viability and collecting the cells for DNA extraction. (B) Reduction in viability to prolonged cytarabine exposure in resistant cell lines. Please click here to view a larger version of this figure.
Figure 7: Preparation of the amplicons representing the integrated shRNA. (A) The binding regions of the primer used in the 1st and the 2nd round of PCR, where the 1st round primers bind to the flanking regions of the integrated shRNA and the second forward primer binds to the loop sequence. The reverse binds to a region of the amplified vector region in the 1st PCR. (B) Illustration of the band size of the PCR product at the end of 1st round PCR (397 bp). (C) The 2nd round PCR product (399bp) was gel-eluted, purified, and given for NGS. Please click here to view a larger version of this figure.
Figure 8: Screening results in the AML cell line (MV-4-11 Ara-C R cell line). After subjecting the DNA to a standard quality control procedure, the samples were processed for NGS analysis. We used CRISPRCloud2, a user-friendly, cloud-based analysis platform, to identify enriched and depleted shRNAs in the pooled shRNA screening. The figure represents enriched or depleted shRNAs targeting epigenetic factors that could mediate cytarabine resistance in AML. Please click here to view a larger version of this figure.
Reagents | Amount to be added | Volume to be added |
OptiMEM Serum free medium | 1 ml | |
PP1- psPAX2 (1µg/µL) | 4 µg | 4 µg/1µg=4 µL |
PP2- pMD2.G (1µg/µL) | 4 µg | 4 µg/1µg=4 µL |
shRNA Library | 10,000 ng | 10000 ng/721ng= 13.8 µL |
Transfection Reagent TransIT-LTI | 30 µL | 30 µL |
Table 1: Preparation of Transfection Plasmid Mix (Calculation for 10cm plate) Please click here to download this Table.
Reagents | Amount for 1reaction |
10X KOD Buffer | 5µL |
DNTPs (2.5mM) | 5µL |
MgCl2 (25mM) | 4µL |
Betaine (5mM) | 5µL |
Forward Primer (10uM) | 1.7µL |
Reverse Primer (10uM) | 1.7µL |
KOD Hot Start Polymerase | 1.5µL |
gDNA (850ng) | |
H20 | |
Total | 50µL |
Table 2: PCR reaction mixture for 1st Round of PCR Please click here to download this Table.
Sample | Total Volume from Pooled PCR reaction | Amount of NT1 solution | Column used |
Sample | 1ml | 2ml | 1 |
PCR Purification: Mix 1 volume of the sample with 2 volumes of Buffer NT1 | |||
Column capacity- 700 µL | |||
Sample | Total Volume from Pooled PCR reaction | Amount of NT1 solution | Column used |
Gel containing DNA | For each 100 mg of agarose gel <2% | 200 µL | 1 |
Table 3: Calculations for purification of the products of 1st of PCR Please click here to download this Table.
Reagents | Calculations for one reaction |
10X KOD Buffer | 5 µL |
DNTP Mix | 5 µL |
MgSO4 (25 mM) | 4 µL |
Betaine (5 mM) | 2 µL |
Secondary Forward Primer (10 µM) | 7.5 µL |
Secondary Reverse Primer (10 µM) | 7.5 µL |
KOD Polymerase | 1.5 µL |
Pooled Primary PCR DNA (500 ng) | |
H20 | |
Total | 50 µL |
Table 4: PCR reaction mixture for 2nd Round of PCR Please click here to download this Table.
RNA interference (RNAi) is extensively used for functional genomics studies, which include siRNA and shRNA screening. The benefit of shRNA is that they can be incorporated into plasmid vectors and integrated into genomic DNA, resulting in stable expression and, thus, more prolonged knockdown of the target mRNA. A pooled shRNA library screening is robust and cost-effective compared to the conventional arrayed screens (siRNA). Identifying the essentiality of a specific class of proteins in a genome-wide screen can be cumbersome. Targeted screening approaches can reveal the essentiality of individual proteins for cancer cell survival or drug resistance within the specific class with reduced noise. This protocol highlights the RNAi screen of epigenetic factors to interrogate both essential and non-essential genes systematically, and it demonstrates the use of this approach as a functional platform allowing the identification of targets responsible for AML drug resistance.
This powerful technology demands certain precautions in the design and execution of the experiments. The first is the choice of promoters that do not undergo significant silencing during the cell culture for the expression of shRNAs as their robust expression is critical for the knockdown of target genes. Second, the maintenance of shRNA representation (the number of cells transduced with each shRNA) is of particular importance. Throughout the shRNA screen, an average representation of >500 cells per shRNA construct augments the potential for identifying the shRNAs that cause phenotypic effects.
The RNAi screening strategy is a competitive growth screen; therefore, it is critical to ensure the target cells are in the logarithmic growth phase and not over-confluent during the screen to minimize loss of representation caused by restricted growth. It is ideal to maintain the cells at 50%-70% confluence to minimize experimental variations. We recommend maintaining consistent lots of media, sera, viral supernatants, drugs, and tissue culture plasticware for all screen replicates to reduce experimental variation.
It is advisable to perform pilot experiments to optimize the transduction efficiency of the lentivirus to calculate the volume of viruses required to achieve one integration of shRNA per cell. To achieve this integration of one shRNA per cell, the transduction efficiency should be less than 30%. As the amount, of viruses to obtain 30% transduced cells may vary between the cell lines used for the experiments, a standardization experiment for each cell line is required. Transduction efficiencies can also be influenced by cell culture conditions and the rate of cell proliferation. Hence, functional titers should be determined in the pilot experiment using the same conditions used in the screening experiment. Maintaining 3-6 replicates is always advisable to get statistically significant hits.
The two major concerns of the shRNA screening approach are the off-target effects and the variable knockdown efficiencies of shRNAs16. These can be overcome by using multiple shRNAs of each gene. Adapting the prokaryotic immune system CRISPR-Cas9 in mammalian cells has allowed easy, efficient, and cost-effective genome editing in cultured human cells. This system has been exploited extensively to perform genome-wide genetic screens similar to those carried out using shRNA. The fate of cells harboring a particular single-guide RNA (sgRNA) sequence can be monitored using deep sequencing. The experimental setup outlined here can be used for CRISPR screening as well.
The authors have nothing to disclose.
This study is funded in part by a Department of Biotechnology grant BT/PR8742/AGR/36/773/2013 to SRV; and Department of Biotechnology India BT/COE/34/SP13432/2015 and Department of Science and Technology, India: EMR/2017/003880 to P.B. RVS and P.B. are supported by Wellcome DBT India Alliance IA/S/17/1/503118 and IA/S/15/1/501842, respectively. S.D. is supported by the CSIR-UGC fellowship, and S.I. is supported by an ICMR senior research fellowship. We thank Abhirup Bagchi, Sandya Rani, and the CSCR Flow Cytometry Core Facility staff for their help. We also thank MedGenome Inc. for helping with the high-throughput sequencing and data analysis.
Reagents | |||
100 bp Ladder Hyper Ladder | BIOLINE | BIO-33025 | |
1kb Ladder Hyper Ladder | BIOLINE | BIO-33056 | |
Agarose | Lonza Seachekm | 50004 | |
Betaine (5mM) | Sigma | B03001VL | |
Boric Acid | Qualigens | 12005 | |
Cell culture plasticware | Corning | as appicable | |
Cytosine β-D-arabinofuranoside hydrochloride | Sigma | C1768-500MG | |
DMEM | MP BIO | 91233354 | |
DMSO | Wak Chemie GMBH | Cryosure DMSO 10ml | |
EDTA | Sigma | E5134 | |
Ethidium Bromide | Sigma | E1510-10 mL | |
Fetal Bovine Serum | Thermo Fisher Scientific | 16000044 | |
Gel/PCR Purification Kit | MACHEREY-NAGEL | REF 740609.50 | |
Gibco- RPMI 1640 | Thermo Fisher Scientific | 23400021 | |
Glacial Acetic Acid | Thermo | Q11007 | |
hCMV GFP Plasmid | Transomics | TransOmics Promoter selection KIT | |
hEF1a GFPlasmid | Transomics | TransOmics Promoter selection KIT | |
HEK 293T | ATCC | CRL-11268 | |
HL60 cell line | ATCC | CCL-240 | |
KOD Hot Start Polymerase | Merck | 71086 | |
Molm13 cell line | Ellen Weisberg Lab, Dana Farber Cancer Institute, Boston, MA, USA | Dana Farber Cancer Institute, Boston, MA, USA | |
MV4-11 cell line | ATCC | CRL-9591 | |
Penicillin streptomycin | Thermo Fisher Scientific | 15140122 | |
psPAX2 and pMD2.G | Addgene | Addgene plasmid no.12260 & Addgene plasmid no. 12259 | |
Qubit dsDNA HS Assay Kit | Invitrogen | REF Q32854 | |
SFFV GFP Plasmid | Transomics | TransOmics Promoter selection KIT | |
shERWOOD-UltrmiR shRNA Library from Transomics | Transomics | Cat No. TLH UD7409; Lot No: A116.V 132.14 | |
Trans-IT-LTI Mirus | Mirus | Mirus Catalog Number: MIR2300 | |
Tris | MP Biomedicals | 0219485591 | |
Trypan Blue | Sigma-Aldrich | T8154-100ML | |
Ultra centrifuge Tubes | Beckman Coulter | 103319 | |
Equipments | |||
5% CO2 incubator | Thermo Fisher | ||
BD Aria III cell sorter | Becton Dickinson | ||
Beckman Coulter Optima L-100K- Ultracentrifugation | Beckman coulter | ||
Centrifuge | Thermo Multiguge 3SR+ | ||
ChemiDoc Imaging system (Fluro Chem M system) | Fluro Chem | ||
Leica AF600 | Leica | ||
Light Microscope | Zeiss Axiovert 40c | ||
Nanodrop | Thermo Scientific | ||
Qubit 3.0 Fluorometer | Invitrogen | ||
Thermal Cycler | BioRad |