High-throughput RNA interference screening using pooled shRNA libraries and next generation sequencing
- David Sims1,
- Ana M Mendes-Pereira1,
- Jessica Frankum1,
- Darren Burgess1,
- Maria-Antonietta Cerone1,
- Cristina Lombardelli1,
- Costas Mitsopoulos1,
- Jarle Hakas1,
- Nirupa Murugaesu1,
- Clare M Isacke1,
- Kerry Fenwick1,
- Ioannis Assiotis1,
- Iwanka Kozarewa1,
- Marketa Zvelebil1,
- Alan Ashworth1Email author and
- Christopher J Lord1Email author
© Sims et al.; licensee BioMed Central Ltd. 2011
Received: 13 June 2011
Accepted: 21 October 2011
Published: 21 October 2011
RNA interference (RNAi) screening is a state-of-the-art technology that enables the dissection of biological processes and disease-related phenotypes. The commercial availability of genome-wide, short hairpin RNA (shRNA) libraries has fueled interest in this area but the generation and analysis of these complex data remain a challenge. Here, we describe complete experimental protocols and novel open source computational methodologies, shALIGN and shRNAseq, that allow RNAi screens to be rapidly deconvoluted using next generation sequencing. Our computational pipeline offers efficient screen analysis and the flexibility and scalability to quickly incorporate future developments in shRNA library technology.
RNA interference (RNAi) facilitates the assessment of gene function by silencing gene expression using synthetic anti-sense oligonucleotides or plasmids. It exploits a physiological mechanism that represses gene expression, primarily by causing the degradation of mRNA transcripts. In mammalian cells, physiological RNAi is primarily mediated by non-protein-coding RNA transcripts, known as microRNAs (miRNAs). miRNAs are produced in a similar manner to mRNAs, but miRNAs are processed into shorter RNA species containing a hairpin structure, known as short-hairpin RNAs (shRNAs). shRNAs are in turn processed into short double-stranded pieces of RNA known as short interfering RNAs (siRNAs). Within the multi-protein RNA-induced silencing complex (RISC), one strand of a siRNA duplex binds a protein-coding mRNA transcript that bears a complementary nucleotide sequence. This interaction allows a nuclease in the RISC to cleave and destroy the protein-coding mRNA, therefore silencing the expression of the gene in a relatively sequence-specific manner.
The experimental use of synthetic siRNAs and shRNA-expressing plasmids has profoundly changed the way in which loss of function experiments can be performed. Previously, techniques that were either more time consuming (gene targeting), or capricious (antisense RNA), were used. Now libraries of RNAi reagents can be purchased and used to silence almost any gene at will. While siRNAs are typically used in multiwell plate-based screening, shRNAs are commonly used for pooled competitive screening approaches, often called barcode screening.
Barcode screening offers improvements in speed and scale compared to plate-based screening. In barcode screening, a large population of cells is infected or transfected with a pool of different shRNA vectors. Cells are then split into two groups and one group is treated differently from the other - for example, with a drug. After this selective pressure is applied, cells are harvested from both populations and integrated hairpins extracted from the genomic DNA of each population by PCR. The relative quantity of each hairpin in the two populations is then compared, to identify those genes that modulate the response to the perturbation in question. For example, in the case of drug screens, hairpins that are over- or under-represented in the drug treated sample compared to the control sample could be considered as targeting genes that modulate sensitivity or resistance to the drug, respectively.
Traditionally, Sanger sequencing has been used as a readout for positive selection screens. However, this approach is costly, time consuming and in general not scalable. In the case of negative selection screens, microarray hybridization is frequently used as a readout [1, 2]. This approach requires the production of custom microarray chips for each library, has a limited dynamic range and is restricted by the varying effectiveness of individual probes. Next generation sequencing (NGS) technologies have recently emerged as a cost-effective means of generating large quantities of sequence data in a short time. Using massively parallel sequencing in place of Sanger sequencing or microarray-based approaches offers several potential advantages in terms of flexibility of input library, scalability and dynamic range.
Already, a small number of laboratories have used shRNA screens coupled to NGS [1, 3–5]. One critical issue that limits the wider exploitation of this technology is the absence of a freely available and simple package for the analysis of shRNA NGS data. With this in mind, we describe here detailed protocols for pooled shRNA screening coupled to NGS screen deconvolution. As part of our optimization of this technology, we have also developed a computational pipeline to analyze NGS data from shRNA screens and describe two open source analysis packages, shALIGN and shRNAseq, designed to simplify barcode screen analysis. Using shRNA pools with engineered depletion, we also assess the sensitivity and reproducibility of this method. As the cost of both shRNA libraries and NGS is rapidly decreasing, these methods and analytical tools may aid the wider adoption of this powerful technology.
Results and discussion
Packaging of hairpin plasmid into lentiviral vectors requires large numbers of packaging cells and high transfection efficiency to ensure faithful representation of the plasmid pool in the viral supernatant. We have successfully employed two approaches to transfection of shRNA plasmids into packaging cell lines, calcium phosphate- and lipid-based transfection. Both methods were routinely used and returned viral supernatants of similar titer (data not shown). cDNA generated from viral supernatant was sequenced and compared to the plasmid DNA to ensure good representation of the library has been achieved (Figure 2a). Typically, cDNA from viral supernatant showed slightly greater variance in hairpin representation than the plasmid pool. Furthermore, hairpin representation at early time points post-viral integration demonstrated better correlation with plasmid representation than with virus (Figure 2a). This suggested that the viral cDNA preparation step was a considerable source of noise and thus plasmid shRNA sequence most likely represents a better reference for starting hairpin representation than virus. This analysis also demonstrated a high concordance between technical replicates, where the same DNA library was sequenced on different GAIIx runs.
Typically, lentiviral stocks were transduced using a multiplicity of infection (MOI) of 0.7 to reduce the likelihood of multiple integrations per cell and the emergence of combinatorial phenotypes. Accurate determination of viral titer in target cell lines allowed subsequent infection of screening cell lines at intended efficiencies. We tested a wide range of breast tumor cell line models and the majority infected at >60% using viral titers of 106 to 107 TU/ml (Figure 2b). Those that did not infect at high efficiency were puromycin selected to give a final green fluorescent protein (GFP)-positive cell population of >90% (Figure 2c).
Viral transduction and cell sampling
Regardless of the design of a particular screen, the manner in which the viral transduction and subsequent cell culture are performed is crucial to the success of the screen. The maintenance of hairpin representation (the number of cells infected with each shRNA) and logarithmic cell growth are of particular importance [3–5]. Throughout all shRNA barcode screens, we maintained an average representation of 1,000 cells per shRNA construct to maximize the potential for phenotypic effects from each shRNA being observed in the final analysis. Previous screens have shown similar levels of shRNA representation to be suitable for the detection of shRNA depletion and enrichment [3, 4]. Since barcode screening is a competitive growth screen, ensuring cells are in log growth at all times during the screen is critical to minimize changes in representation caused by localized restriction of cell growth due to over-confluence. Consequently, we recommend ensuring that cells are never allowed to achieve more than 70% confluence.
After viral integration and puromycin selection were complete, cultures were divided into two or more sets, depending on the experimental design. For example, in a typical drug sensitivity/resistance screen, cultures are divided into reference (vehicle treated) and test (drug treated) sets. Alternatively, in a simple viability screen, a sample of cells can be taken and stored for analysis at each passage, to generate a viability time-course. One of the strengths of the lentiviral system is the stable integration of hairpins; this allows the use of longer experimental time-courses than could generally be performed using siRNA screening. As a consequence, final screen results were typically assessed 2 to 3 weeks after dividing the cells into two arms. Every time the cultures were divided or sampled, aliquots were taken to assess the cell number (to construct growth curves) and the percentage of GFP-positive cells (to assess the number of cells required to maintain hairpin representation). To minimize screen variability, we use the same passage cells for each screen replicate. We also maintain consistent batches of media, serum, viral supernatant and tissue culture plasticware for all screen replicates, again to minimize experimental variation.
To enable sufficient representation of each shRNA in the screening pool, multiple PCR reactions were performed in parallel to generate the sequencing library from each shRNA pool. For example, to maintain a representation of 1,000 cells per shRNA in a pool of 10,000 shRNAs at an MOI of 0.7, PCR amplification from 1 × 107 cells is required. Since a diploid human cell contains approximately 6 pg genomic DNA, we performed PCR amplification from 60 μg total genomic DNA using 30 parallel PCR reactions, each with 2 μg of genomic DNA.
DNA quantification for next generation sequencing
First, we utilized the complement of the Illumina adapter sequences (p5 and p7, common to all Illumina sequencing libraries) as amplification primers in an intercalating dye-based quantitative PCR (qPCR) assay (SybrGreen). This is similar to an approach applied to the quantification of 454 Roche pyro-sequencing samples . Second, for the quantification of DNA containing hairpin inserts in shRNA-derived Illumina libraries, we established a second strategy that utilized a dual-labeled probe (DLP) hydrolysis qPCR assay (Taqman). Here we designed a DLP complementary to a constant internal region in shRNA-specific Solexa PCR products. We based this approach on a reported alternative Taqman assay defined for an established Illumina sequencing application, pair-end RNA sequencing (RNAseq) . Overall, the implementation of both of these innovative strategies allowed us to reliably quantify shRNA templates prior to massively parallel sequencing, leading to high numbers of mapped reads passing quality filters.
Sequencing data analysis
Detection of depleted hairpins in reference depletion screens
Although NGS performs better than microarray hybridization for screen deconvolution, the barcode screening format is subject to a significant degree of stochastic noise regardless of the hairpin frequency detection method, and biological replication of screens is almost certainly required to overcome this. We performed two biological replicates of the reference depletion screen. There was a good correlation between log normalized read counts from replicate screens (reference arm r2 = 0.89). Similarly, a high correlation of depletion-reference log ratios was observed between biological replicates (r2 = 0.92; Figure 6c), suggesting that screens are highly reproducible. As expected, higher levels of noise were observed between biological replicates than between technical replicates.
Using the reference-depletion experimental framework, we also titrated the number of PCR cycles required to give faithful representation of the starting material. This demonstrated that PCR amplification acted as a source of noise in the experiment and that reducing the number of PCR cycles resulted in a decreased false positive rate at a set false negative rate of 5% (Figure 6d). Consequently, we opted to use 26 PCR cycles as this generated a visible band on the gel but still resulted in a relatively low error rate.
Comparison of reference-depletion screens at different read depths
To test sensitivity in a genuine screen setting using different screening pool sizes, we performed a series of screens for viability in MCF7 cells. To assess screen performance, we used a number of genes where previous observations had shown that siRNA silencing inhibited MCF7 cells as measured using a viability assay based on cellular ATP levels . We validated the viability effects of a set of GIPZ library shRNAs targeting these genes using single hairpin GFP competition studies. We also validated a number of negative control non-targeting shRNAs using this method. This revealed a total of six shRNAs, which caused a >50% depletion in GFP-positive cells over 2 weeks, along with 11 non-targeting hairpins that showed no viability phenotype in the same period (Additional file 3).
These results suggested that the scoring system performs consistently regardless of pool size and the screen is as sensitive in a 10,000 shRNA pool as in a 2,000 shRNA pool. Therefore, to maximize throughput and minimize expense, we decided to use 10,000 shRNA pools for screening. Since the genome-wide library encompasses approximately 60,000 reagents, we reasoned that a pool size of approximately 10,000 shRNAs would be an ideal choice for this library. This would enable a genome-wide screen to be run on a single eight-lane Illumina flow-cell along with appropriate sequencing controls. Furthermore, a 10,000 shRNA pool size would require approximately 15,000,000 cells to be maintained in tissue culture over the course of the experiment, assuming an MOI of 0.7 and a representation of 1,000 cells per hairpin. This pool size enabled us to maintain the simplicity and reproducibility of the associated tissue culture work, which might become cumbersome in larger pool sizes.
To establish the reproducibility of the screening methodology, we repeated the screen for cell viability in MCF7 cells in a total of four biological replicates. The replicates showed high correlation despite being performed by different researchers at different times, suggesting that the screening method is robust (Figure 7d). Furthermore, >98% of library hairpins were identified in all screen replicates, compared to 49% of half-hairpin probes and 82% of barcode tiling probes having intensity above background in microarray-based approaches . This suggests a greatly improved sensitivity for the NGS profiling approach.
Here we describe detailed and optimized methods for high-throughput shRNA screening using NGS. Using massively parallel sequencing in place of microarray hybridization for deconvolution of shRNA barcode screen results offers several advantages. Firstly, we demonstrate that NGS profiling offers greater sensitivity, enabling more hairpins to be detected above background. Secondly, we show that NGS profiling has a better dynamic range when compared to literature examples of microarray-based approaches, since sequencing cannot be easily saturated in the same way that hybridization can, thus enabling a greater separation of true effects from background noise. Thirdly, sequencing of tags is both scalable and flexible, enabling new hairpins to be incorporated into the work schema without having to print a new batch of custom microarrays. Finally, our computational pipeline offers the ability to identify sequences where bases are not called due to sequencing errors, or mutated during PCR, since short reads can be matched inexactly to the reference library of shRNA barcode sequences.
We do note that one of the major limitations of any shRNA library is the efficacy of gene silencing, an issue we have not addressed here. However, the ability to rapidly assess silencing capacity and thus develop algorithms that would predict effective gene knockdown is improving . In addition, an increase in the redundancy of shRNA libraries (the number of different shRNA constructs per gene) will also improve the general effectiveness of pooled RNAi screens. Finally, improvements in sequencing technology will undoubtedly increase read number per lane; these developments will thus enable the use of larger shRNA pool sizes that could accommodate libraries with higher levels of redundancy. Nevertheless, even with the existing commercial shRNA libraries and also the ever-increasing availability of cost-effective NGS, the methods we describe here should enable the wider applicability of this powerful technology.
Materials and methods
Although the following methods are suitable for most viral shRNA libraries, the work described here used the Thermo Scientific Open Biosystems GIPZ Lentiviral human shRNAmir library (version 2). These methods could be used for other shRNA libraries, such as the RNAi consortium library , with altered PCR and sequencing primers. The GIPZ Lentiviral human shRNAmir library we used encompasses 61,416 distinct hairpin constructs targeting 15,739 human protein coding genes (based upon the Ensembl 56 build). In the GIPZ vector the 19-nucleotide siRNA sense sequence is inserted into a human mir-30 backbone . shRNA sequences are designed to have destabilized 5' ends in the antisense strand to encourage stand-specific incorporation into the RISC . The vector backbone includes a GFP-coding sequence that is transcribed as part of a bicistronic transcript with the shRNA sequence allowing the visualization of shRNAmir expressing cells, and a puromycin resistance marker for selecting infected cells.
LB media (1 ml) containing 50 μg/ml ampicillin was added to 96-deep-well microplates using a multidrop (Thermo Fisher, Waltham, MA, USA), and 1 μl of each bacterial inoculant (extracted from fully thawed 96-well glycerol stocks) was seeded per well using a Beckman FX robot. Culture plates were stacked evenly in shaking incubators (200 rpm) ensuring even air circulation. Following 16 hours of incubation at 37°C, cultures from a single batch of plates were pooled and plasmid DNA was isolated using a Plasmid Maxi kit (Qiagen, Crawley, UK) and normalized to a standard concentration (0.5 μg/ml). Plasmid DNA pools were combined in equal concentrations to create screening pools of different complexities.
For both calcium phosphate- and lipid-based transfection using lipofectamine 2000, 293T cells were seeded onto 10 cm dishes (50 to 80% confluent), and co-transfected with the desired shRNA library pool, along with the packaging plasmids psPAX2 and pMD2.G , in a Safety Category II facility. Thirty-six hours post-transfection, supernatant was collected, supplemented with 4 μg/ml polybrene and filtered through a 0.45 μm membrane. Viral supernatant from multiple 10 cm dishes was pooled, aliquoted and stored at -80°C for future use. Determination of viral titer allowed subsequent infection of screening cell lines at intended efficiencies. For this, both target and reference 293T cells were infected with 1:2 serial dilutions of the virus. Virus was removed 24 hours later and cells incubated for a further 48 hours at 37°C, after which the proportion of GFP-positive cells was determined by fluorescence-activated cell sorting (FACS) to provide an estimate of the fraction infected and viral titer.
Viral transduction and cell sampling
Cells were infected with viral pools using a MOI of 0.7 and media replenished after 16 hours. Seventy-two hours post-infection, when viral integration was presumed complete, cells were exposed to 1 mg/l puromycin for 2 days to select for cells with viral integration. Following puromycin selection, cultures were divided into two or more replica cultures, and continuously cultured for 2 to 3 weeks, dividing cultures regularly to ensure continued logarithmic growth and to maintain hairpin representation at approximately 1,000 cells per shRNA [4, 5]. Cells were harvested in suspension and either stored frozen at -20°C or processed immediately as described below.
Genomic DNA extraction and purification from cultured cells was carried out using a Gentra Puregene kit (Qiagen). shRNA sequences integrated into genomic DNA were recovered by PCR amplification using the following primers and a procedure based upon that described by Hannon, Elledge and Lowe [4, 5]: p5+mir30, 5'-AATGATACGGCGACCACCGACTAAAGTAGCCCCTTGAATTC-3'; p7+Loop, 5'-CAAGCAGAAGACGGCATACGATAGTGAAGCCACAGATGTA-3'. PCR was performed on PTC-225 3DNA Engine Tetrads (Bio-Rad, Hemel Hempstead, UK) using Amplitaq Gold polymerase and 20 to 33 cycles of denaturation (95°C, 30 s), annealing (52°C, 45 s) and extension (72°C, 60 s). In general, 2 μg genomic DNA was used per 100 μl PCR reaction and 48 PCR reactions were performed per 10,000 shRNA pool. PCR products from multiple parallel reactions were subsequently pooled, concentrated and purified using a QIAquick gel extraction kit (Qiagen).
Quantitative PCR DNA quantification
For absolute DNA quantification, real-time PCR was performed using two alternative strategies: SybrGreen and TaqMan. Both approaches used the appropriate universal 2X PCR mix from Applied Biosystems and the following oligonucleotides at 250 nM each: P5, AATGATACGGCGACCACCGA (20-mer); and P7, CAAGCAGAAGACGGCATACGA (21-mer). For TaqMan only, 250 nM of a 21-mer DLP was included (CCCTTGAATTCCGAGGCAGTA), with 5' reporter (6FAM fluorescein) and 3' quencher (tetramethylrhodamine (TAMRA)) as detectors, and ROX passive reference dye. Forty-cycle measurements of triplicate 25 μl reactions containing 10% template sample at a theoretical concentration of 10 pM (as defined by Agilent Bioanalyzer) were carried out in an ABI 7900 instrument. Solexa library concentrations were then inferred by comparing measurements to the standard curve using the Sequence Detection System (SDS) v2.2.1 software. To generate the standard curve, a 10-fold dilution series of a standard 100 nM sample (7.7 ng/μl for 117 bp) was used (final range 0.1 to 100 pM).
Following quantification, denatured shRNA-seq libraries (3 pM in NaOH) were pumped through eight-lane flowcell channels using the cluster station of an Illumina GAIIx sequencing platform. Bridge PCR was executed using the manufacturer's protocol (Amplification-linearization-blocking and multiple primer hybridization version 3). The sequencing primer (TAGCCCCTTGAATTCCGAGGCAGTAGGCA) was designed to sequence from two bases upstream of the 19 bp shRNA sense sequence. Twenty-six cycles of sequencing-by-synthesis (single read) were performed on an Illumina GAIIx according to the manufacturer's protocol (version 4).
Image analysis and base calling
Raw image data were analyzed using GA pipeline v1.4. PhiX was run as a control to ensure correct phasing. Base calling was performed by the Bustard package using the Chastity filter with a threshold of 0.6. The Chastity filter was applied to bases 3 to 21 of the 26 bases sequenced (the shRNA sense sequence) only. A maximum of two uncalled bases were allowed in these 19 bases. FASTQ files from this study are available at the ROCK web page  and at the European Nucleotide Archive .
shRNA library alignment
FASTq files generated by the GA pipeline were mapped to reference shRNA libraries using a bespoke Perl program, shALIGN. shALIGN trims short reads to the 19 bp sense sequence based on user-defined base positions, and groups identical reads disregarding sequence quality scores. shALIGN employs a hamming distance algorithm to align binned short read sequences to a user-defined reference shRNA library provided in a standard tab-delimited text format. The user is able to specify the maximum number of mismatches permissible between the short read and the reference. This alignment is equivalent to using Bowtie with the flowing flags: --best --strata -v 2 -m 1 -a. However, using shALIGN negates the need to construct and index a specific Bowtie reference library by inserting the sense sequence of each shRNA into a longer sequence. It also saves the need to write a script to parse the Bowtie output to count the number of reads mapped to each library construct at each edit distance. For the GIPZ library this was set at two mismatches as the vast majority of library sequences were separated by an edit distance of greater than two. All short reads that match to more than one library sequence at the same distance were excluded from further analysis, and logged as ambiguous. The shALIGN program outputs the total number of short reads matching to each library hairpin in each lane. We routinely load screen results into the ROCK breast cancer functional genomics database , where screen reagents are fully annotated to the latest genome build. ROCK provides a mechanism for the sharing and publication of raw and processed RNAi screen results. Source code and linux binaries for the shALIGN program are available from the ROCK web page .
To facilitate statistical analysis of screen results, we have developed a novel R package, shRNAseq. This package is based on the NChannelSet class in Bioconductor , originally designed to handle microarray gene expression data, and serves as a single user-friendly package encompassing all of the steps required. shRNAseq reads in matrices of short read counts per shRNA construct, generated by shALIGN, along with annotations of shRNA constructs and sequenced samples.
The package is designed to compare a pair of related screen conditions. For example, in a viability screen one would compare the hairpin representation in the starting population (which could be from the plasmid, virus or an early time point post-viral integration) to the hairpin representation at the end of the screen (for example, after 2 weeks). Alternatively, for a drug screen one would compare the drug and vehicle treated arms at the same time point. The package can analyze multiple screen replicates simultaneously. Each screen condition is loaded into the package separately, and annotated using files describing the shRNA constructs and the sample treatments. After data loading, read counts per hairpin are log2 transformed and then the ratio of the two screen conditions is calculated. This log ratio is normalized using a loess fit to the log mean read count for each screening pool. The distribution of normalized scores per pool is then rescaled by dividing by the pool MAD. If appropriate, screen biological replicates can be quantile (rank) normalized to ensure identical distributions, and scores can be summarized across replicates using a variety of methods (for example, median, regularized t-test). Finally, the package plots screen distributions before and after normalization and reports a table of normalized scores. The R package is accompanied by a detailed vignette describing the methodology and usage at the ROCK web page .
Hit detection was performed using three different methods. In the first method, replicate scores for each hairpin were summarized using the median and a hit threshold estimated from a quantile-quantile plot, identifying hairpin scores that significantly differed from the normal distribution. Second, the Gene Set Analysis R package  was used to look for enrichment or depletion of sets of hairpins targeting the same gene. This approach makes use of the hairpin redundancy within the library and works best in libraries containing multiple hairpins per gene. Finally, the RIGER algorithm  in the GENE-E java package  was also used to look for enrichment of shRNAs targeting the same gene. Here we used the log fold change metric and weighted average method. Source code and documentation for the shRNAseq R-package are available from the ROCK web page .
Full methods, including a detailed screening manual, are provided in Additional file 4.
Genome Analyzer IIx
green fluorescent protein
median absolute deviation
multiplicity of infection
next generation sequencing
RNA-induced silencing complex
short hairpin RNA
shRNA encoded within a miRNA
short interfering RNA
This work was supported by Breakthrough Breast Cancer, American Association of Cancer Research Stand Up to Cancer, The Breast Cancer Research Foundation and Cancer Research UK. We thank Drs. Greg Hannon, Joannas Zuber and Scott Lowe for sharing methods prior to publication. We acknowledge NHS funding to the Royal Marsden Hospital National Institute for Health Research Biomedical Research Centre.
- Burgess DJ, Doles J, Zender L, Xue W, Ma B, McCombie WR, Hannon GJ, Lowe SW, Hemann MT: Topoisomerase levels determine chemotherapy response in vitro and in vivo. Proc Natl Acad Sci USA. 2008, 105: 9053-9058. 10.1073/pnas.0803513105.PubMedPubMed CentralView ArticleGoogle Scholar
- Schlabach MR, Luo J, Solimini NL, Hu G, Xu Q, Li MZ, Zhao Z, Smogorzewska A, Sowa ME, Ang XL, Westbrook TF, Liang AC, Chang K, Hackett JA, Harper JW, Hannon GJ, Elledge SJ: Cancer proliferation gene discovery through functional genomics. Science. 2008, 319: 620-624. 10.1126/science.1149200.PubMedPubMed CentralView ArticleGoogle Scholar
- Bassik MC, Lebbink RJ, Churchman LS, Ingolia NT, Patena W, LeProust EM, Schuldiner M, Weissman JS, McManus MT: Rapid creation and quantitative monitoring of high coverage shRNA libraries. Nat Methods. 2009, 6: 443-445. 10.1038/nmeth.1330.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuber J, McJunkin K, Fellmann C, Dow LE, Taylor MJ, Hannon GJ, Lowe SW: Toolkit for evaluating genes required for proliferation and survival using tetracycline-regulated RNAi. Nat Biotechnol. 2011, 29: 79-83. 10.1038/nbt.1720.PubMedPubMed CentralView ArticleGoogle Scholar
- Silva JM, Marran K, Parker JS, Silva J, Golding M, Schlabach MR, Elledge SJ, Hannon GJ, Chang K: Profiling essential genes in human mammary cells by multiplex RNAi screening. Science. 2008, 319: 617-620. 10.1126/science.1149185.PubMedPubMed CentralView ArticleGoogle Scholar
- Quail MA, Kozarewa I, Smith F, Scally A, Stephens PJ, Durbin R, Swerdlow H, Turner DJ: A large genome center's improvements to the Illumina sequencing system. Nat Methods. 2008, 5: 1005-1010. 10.1038/nmeth.1270.PubMedPubMed CentralView ArticleGoogle Scholar
- Meyer M, Briggs AW, Maricic T, Hober B, Hoffner B, Krause J, Weihmann A, Paabo S, Hofreiter M: From micrograms to picograms: quantitative PCR reduces the material demands of high-throughput sequencing. Nucleic Acids Res. 2008, 36: e5-10.1093/nar/gkn201.PubMedPubMed CentralView ArticleGoogle 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: R25-10.1186/gb-2009-10-3-r25.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMedPubMed CentralView ArticleGoogle Scholar
- Boettcher M, Fredebohm J, Moghaddas Gholami A, Hachmo Y, Dotan I, Canaani D, Hoheisel JD: Decoding pooled RNAi screens by means of barcode tiling arrays. BMC Genomics. 2010, 11: 7-10.1186/1471-2164-11-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Iorns E, Turner NC, Elliott R, Syed N, Garrone O, Gasco M, Tutt AN, Crook T, Lord CJ, Ashworth A: Identification of CDK10 as an important determinant of resistance to endocrine therapy for breast cancer. Cancer Cell. 2008, 13: 91-104. 10.1016/j.ccr.2008.01.001.PubMedView ArticleGoogle Scholar
- Fellmann C, Zuber J, McJunkin K, Chang K, Malone CD, Dickins RA, Xu Q, Hengartner MO, Elledge SJ, Hannon GJ, Lowe SW: Functional identification of optimized RNAi triggers using a massively parallel sensor assay. Mol Cell. 2011, 41: 733-746. 10.1016/j.molcel.2011.02.008.PubMedPubMed CentralView ArticleGoogle Scholar
- Root DE, Hacohen N, Hahn WC, Lander ES, Sabatini DM: Genome-scale loss-of-function screening with a lentiviral RNAi library. Nat Methods. 2006, 3: 715-719. 10.1038/nmeth924.PubMedView ArticleGoogle Scholar
- Silva JM, Li MZ, Chang K, Ge W, Golding MC, Rickles RJ, Siolas D, Hu G, Paddison PJ, Schlabach MR, Sheth N, Bradshaw J, Burchard J, Kulkarni A, Cavet G, Sachidanandam R, McCombie WR, Cleary MA, Elledge SJ, Hannon GJ: Second-generation shRNA libraries covering the mouse and human genomes. Nat Genet. 2005, 37: 1281-1288.PubMedGoogle Scholar
- Klages N, Zufferey R, Trono D: A stable system for the high-titer production of multiply attenuated lentiviral vectors. Mol Ther. 2000, 2: 170-176. 10.1006/mthe.2000.0103.PubMedView ArticleGoogle Scholar
- ROCK - Breast Cancer Functional Genomics. [http://rock.icr.ac.uk/software/shrnaseq.jsp]
- The European Nucleotide Archive, Accession number ERP000908. [http://www.ebi.ac.uk/ena/data/view/ERP000908]
- Sims D, Bursteinas B, Gao Q, Jain E, MacKay A, Mitsopoulos C, Zvelebil M: ROCK: a breast cancer functional genomics resource. Breast Cancer Res Treat. 2010, 124: 567-572. 10.1007/s10549-010-0945-5.PubMedView ArticleGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView ArticleGoogle Scholar
- Gene Set Analysis. [http://www-stat.stanford.edu/~tibs/GSA/]
- Luo B, Cheung HW, Subramanian A, Sharifnia T, Okamoto M, Yang X, Hinkle G, Boehm JS, Beroukhim R, Weir BA, Mermel C, Barbie DA, Awad T, Zhou X, Nguyen T, Piqani B, Li C, Golub TR, Meyerson M, Hacohen N, Hahn WC, Lander ES, Sabatini DM, Root DE: Highly parallel identification of essential genes in cancer cells. Proc Natl Acad Sci USA. 2008, 105: 20380-20385. 10.1073/pnas.0810485105.PubMedPubMed CentralView ArticleGoogle Scholar
- GENE-E Java Package. [http://www.broadinstitute.org/cancer/software/GENE-E/]
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.