Identifying novel drug targets that transition from pre-clinical testing to human trials is a scientific priority. To that end, here we describe a functional genomics approach for examining the impact of gene depletion on cancer cell line spheroids, which more appropriately model human cancers in vivo.
The identification of functional driver events in cancer is central to furthering our understanding of cancer biology and indispensable for the discovery of the next generation of novel drug targets. It is becoming apparent that more complex models of cancer are required to fully appreciate the contributing factors that drive tumorigenesis in vivo and increase the efficacy of novel therapies that make the transition from pre-clinical models to clinical trials.
Here we present a methodology for generating uniform and reproducible tumor spheroids that can be subjected to siRNA functional screening. These spheroids display many characteristics that are found in solid tumors that are not present in traditional two-dimension culture. We show that several commonly used breast cancer cell lines are amenable to this protocol. Furthermore, we provide proof-of-principle data utilizing the breast cancer cell line BT474, confirming their dependency on amplification of the epidermal growth factor receptor HER2 and mutation of phosphatidylinositol-4,5-biphosphate 3-kinase (PIK3CA) when grown as tumor spheroids. Finally, we are able to further investigate and confirm the spatial impact of these dependencies using immunohistochemistry.
Solid tumors display significant histological, genetic and micro-environmental intra-tumor heterogeneity, which presents clinicians with a significant challenge in being able to treat patients successfully. The majority of models used to identify novel targeted therapies do not incorporate many of these features. Indeed, current targeted therapies utilized in the clinic have been developed in the last decade employing screening approaches that rely on cancer cell lines grown under two-dimensional (2D) culture conditions. Although this has brought about various successes, such as receptor tyrosine kinase inhibitors, it is becoming apparent that more complex models of cancer are required to fully appreciate the contributing factors that drive tumorigenesis in vivo and increase the number of new therapies that make the transition from pre-clinical models to clinical trials. Moreover, it is now well appreciated that 2D culture systems fail to reflect in vivo behavior 1,2. For example, in poorly vascularized tumors, as the demand for oxygen and nutrients within the microenvironment outstrip supply, regions of high and low deliveries develop. The presence of low oxygen (hypoxia) in tumors, as detected by the immunohistochemical staining of tumor sections for established hypoxic markers such as carbonic anhydrase IX (CAIX), correlates with a poorer clinical outcome in breast cancer 3,4 Thus incorporating features such as hypoxia into screening models may enhance our capacity to discover novel drug targets that will be more efficacious in vivo. Indeed, successfully targeting aggressive tumors that contain hypoxia is a clinical priority 5.
The altering of traditional 2D siRNA screening, in an attempt to more accurately recapitulate elements of the conditions encountered by cancer cells in the tumor microenvironment, has led to the identification of several genes that have been found to be important for tumor growth in vivo. These include functional genomic screens performed under low serum conditions 6, hypoxic conditions 7 and in combination 8. For example, silencing of 6-Phosphofructo-2-Kinase/Fructose-2,6-Biphosphatase 4 (PFKFB4), a protein responsible for regulating carbon entering glycolysis, only induced apoptosis in prostate cancer lines derived from metastasis when grown in low serum. Silencing of PFKFB4 in normal prostate cell lines under the same conditions had no effect, whereas, depletion of PFKFB4 completely ablated the growth of prostate cancer cell line xenografts 6.
In an extended panel of breast cancer cell lines, the silencing of mono-carboxylase transporter 4 (MCT4) preferentially led to the reduction of cell line growth under conditions of low oxygen. This vulnerability was validated in vivo in breast cancer cell line orthotropic xenografts. Perhaps most strikingly, silencing of Acetyl-CoA synthetase 2 (ACSS2), an enzyme responsible for converting acetate into acetyl-CoA, reduced cancer cell number under nutrient stress conditions (low oxygen and serum) but had little or no effect under normal culture conditions 8. Ablation of ACSS2 impacted the growth of breast and prostate cancer xenografts suggesting that nutrient gradients do not exist in isolation within the tumor microenvironment and that cells that reside in these regions are essential for tumor progression 8. Furthermore, ACSS2 was also found to be important in glioblastoma and hepatocellular carcinomas 9,10, suggesting that increased ACSS2 activity in tumors could be a fundamental mechanism that supports growth under unfavorable conditions.
Collectively, these studies prove that recapitulating the conditions encountered in vivo and performing siRNA screens allows for the identification of genes essential for cancer survival. As well as impacting 2D cancer cell growth under nutrient stress conditions, depletion of target genes in these studies inhibited cancer cell line spheroid growth, mirroring what was observed in tumor xenografts 6,8. Thus, cancer cell line spheroids contain several of the conditions encountered in the tumor microenvironment that impart sensitivity to ACSS2 silencing. Indeed, spheroids display nutrient gradients (Serum and oxygen), changes in pH, three-dimensional (3D) cell-cell contact, but also, alterations in proliferative celerity with cells undergoing cell cycle arrest and apoptosis. This is exemplified by the presence in cancer spheroids of necrotic regions, a feature not found in traditional 2D culture.
Cancer cell spheroids have already been used as more biologically relevant models to screen small molecule inhibitors but this only allows for the validation of compound efficacy or the repurposing of compounds originally design for other diseases 11. Current spheroid screening methods do not allow for the analysis of specific gene depletion in a high-throughput of high-content manner. Here, we describe for the first time, a functional genomics pipeline for uncovering specific gene dependencies utilizing small interfering RNA (siRNA) technology in cancer cell line spheroids. We designed a bespoke library with siRNAs targeting the 200 most frequently mutated genes in human breast cancers and assessed the impact of gene depletion on spheroid size and metabolic activity in BT474 breast cancer spheroids. We were able to robustly and reproducibly detect the impact of ERBB2 and PIK3CA silencing in 3D cultures. Moreover, we could assess the impact of gene depletion on the spatial architecture of BT474 spheroids using immunohistochemistry.
1. Preparation of 96-well siRNA Plates
NOTE: The outer edges of a 96-well plate are more prone to evaporation compared to other wells, thus limit the amount of siRNAs to 60 per 96-well plate. Fill in outer wells with plain medium or PBS to limit this. In addition, a validation screen of hits is recommended to alleviate plate biases.
2. Reverse Transfection of Cell Line
3. Automated Image Acquisition
4. Determining Cell Viability
Spheroid assays in ultra-low attachment 96-well plates provide a high-throughput phenotypic assessment for potential oncogenicity in a context that more readily recapitulates the physiological conditions found in tumors in vivo. Indeed, the cancer cell lines MCF10DCIS.com and BT474 form tight spheroid structures (Figure 1A) and immunohistological investigation of spheroid sections showed distinct spatial changes in cellular and nuclear morphology. Over time, some spheroids such as BT474 spheroids develop necrotic regions, a common feature of aggressive solid tumors (Figure 1B). Some spheroids do not develop necrotic cores, such as the MDA-MB-231 cell line, but do display marked variation in the proliferative marker Ki67, which inversely correlates with cleaved caspsase-3 expression, a marker of apoptosis (Figure 1C). To establish that multiple cell lines are indeed receptive to siRNA-mediated gene silencing BT474, MCF10DCIS.com, MDA-MB-231 and JIMT1 cells were reverse transfected with siRNA for seven days. The presence of transfection reagent (mock) or transfection of control siRNAs had no effect on spheroid viability, while silencing of the essential gene Ubiquitin B (UBB) significantly reduced spheroid viability in all the cell lines tested (Figure 1D).
We designed a human siRNA library that encompassed the most frequently mutated genes in unselected, ER+, HER2+ and triple negative breast cancers. This library consists of genes of known function, such as MYC, PIK3CA and TP53 and those whose contribution to carcinogenesis is not established. The screen also contains several non-targeting control siRNAs (Control #1, Control #2) and siRNAs targeting essential genes, such as PLK1 and UBB, which act as killing controls (Table 1). We chose to use the breast cancer cell lines BT474 as they readily form spheroids without the addition of reconstituted basement membrane, are established workhorse cell lines and have a known genomic architecture. For example, BT474 cells are positive for the estrogen receptor (ER+), overexpress the human epidermal growth factor receptor 2 (HER2+) and harbor mutations in TP53 (E285K) and PIK3CA (K111N)13.
Employing the protocol outlined above, we monitored the impact gene depletion had on spheroid size and viability after seven days of siRNA reverse transfection (Figure 1E). Interestingly, the majority of genes did not have a significant effect on spheroid area or viability (Figure 2A). Silencing of FOXO3, PIK3CA, ERBB2 and SF3B1 resulted in the most significant reproducible reduction in spheroid size. This reduction was also observed in spheroid viability after ERBB2 and SF3B1 silencing. Encouragingly, we confirmed the impact of PIK3CA, ERBB2 and SF3B1 silencing on spheroid size using bright-field microscopy (Figure 2B). We previously identified SF3B1 as an essential gene in numerous cell line models and thus siSF3B1 represents a good killing control in addition to UBB 14. Interestingly, of all 200 genes only the silencing of E-Cadherin resulted in a significant increase in spheroid viability (Figure 2A). Investigation of spheroid morphology showed that E-Cadherin silencing resulted in a complete breakdown of spheroid architecture, with viable cells resting on the bottom of the low attachment well (Figure 2B). Manual reinvestigation of the spheroid volume screen data showed that this had also been observed but had been eliminated from area quantification due to the object being above the set size restrictions. As previously highlighted, BT474 cells overexpress the receptor tyrosine kinase HER2 and harbor an oncogenic mutation in PIK3CA (K111N). We confirmed that silencing of ERBB2 and PIK3CA resulted in reduced spheroid viability, while transfection with non-targeting controls had no effect (Figure 2C).
Next we investigated the impact of siRNA depletion on spheroid histology. BT474 spheroids were reverse transfected with non-targeting control siRNA and siRNAs targeting PIK3CA, ERBB2 and UBB. Silencing of ERBB2 and UBB resulted in a reduction of the pro-proliferative marker Ki67 compared to control siRNA (Figure 2D). The activation of the pro-apoptotic marker cleaved caspase-3 was only observed after UBB silencing, suggesting that depletion of HER2 and PIK3CA did not result in apoptosis but were cytostatic rather than cytotoxic. Indeed, silencing of HER2 and PIK3CA did result in an increase in the protein expression of the cell cycle arrest protein p27 compared to control transfected spheroids.
Taken together, these results show that BT474 cells are driven by oncogenic HER2 and PIK3CA signaling when grown as 3D spheroids. More importantly, these results show that it is possible to design and implement a bespoke siRNA screening library of hundreds of genes in cancer cell line spheroids robustly and reproducibly.
Figure 1: Cell Line Optimization for 3-dimensional Growth. A. The breast cancer cell lines, MCF10DCIS.com and BT474 were cultured in low attachment plates for 7 days. Brightfield representative images were taken using an inverted microscope. Scale bars = 100 µm. B. MCF10DCIS.com and BT474 spheroids were cultured for 28 days. 100 μL of fresh media was replenished every 3 – 4 days. Spheroids were fixed in 3.8% formaldehyde, embedded, sectioned and stained with hematoxylin and eosin (H&E). Representative images are shown at low and high magnification. Scale bars represent 100 µm and 33 µm, respectively. C. MDA-MB-231 spheroids were cultured for 21 days. 100 μL of fresh media was replenished every 3 – 4 days. Spheroids were fixed in 3.8% formaldehyde, embedded, sectioned and stained with Ki67 and cleaved caspase-3. Representative images are shown. Scale bars = 100 µm. D. BT474, MCF10DCIS.com, MDA-MB-231, and JIMT1 cell lines were reverse transfected with mock (transfection reagent only) and the control siRNAs and UBB, ultra-low attachment plates were then spun to form spheroids. Fresh medium (100 µl) was added on days 1 and 4. After 7 days, cell viability was quantified. Data represent the mean ± SD of two independent biological replicates performed in triplicate normalized to control #1. Statistical significance was calculated using an unpaired Students t-test (p< 0.05). E. A flow diagram summarizing the reverse transfection protocol used to functionally interrogate gene dependency in cancer cell line spheroids. Please click here to view a larger version of this figure.
Figure 2: Functional Genomic Investigation of BT474 Spheroids Uncovers Oncogenic Dependencies. A. BT474 cells were reverse transfected with the 200-gene human siGENOME siRNA library in triplicate. Spheroid size and viability was observed. Raw data values were plate median normalized and z-scores were calculated to identify significant outliers greater than 1.7x the standard deviation of the plate median 15. Note outlying genes ERBB2, SF3B1, PLK1 and UBB and E-cad. siRNAs that significantly increased or reduced spheroid area and viability are shaded in blue and red, where red depicts the control siRNA's. B. Cells were reverse transfected with non-targeting control, E-Cadherin (E-Cad), PIK3CA, ERBB2 or UBB siRNA and then spun in an ultra-low attachment plate to form spheroids. After 7 days, brightfield representative spheroid images were taken using an inverted microscope. Scale bars = 100 µm. C. Cells were reverse transfected with non-targeting control, PIK3CA, ERBB2 or UBB siRNA and then spun in an ultra-low attachment plate to form spheroids. After 7 days, cell viability was quantified. Data represent the mean ± SD of two independent biological replicates performed in triplicate normalized to control #1. Statistical significance was calculated using an unpaired Students t-test (p < 0.05). D. After 7 days, spheroids were fixed, embedded, sectioned and stained for H&E, Ki67, Cleaved Caspase-3, and p27. Representative images are shown. Scale bars = 100 µm. Note the H&E of the siControl spheroids appear smaller due to an artefact of processing and individual intact spheroids chosen for staining. However this does not detract from the changes observed with the scanned images and cell viability results. Please click here to view a larger version of this figure.
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
Non-targeting 1 | TP53 | DST | KMT2C | Untreated 1 | FCGBP | ARID1B | FBXM7 | TTC40 | Non-targeting 2 | ||
GATA3 | TTN | MUC12 | MUC4 | AHNAK | HUWEI1 | DNAH11 | ITPR2 | ABCA13 | CREBBP | ||
MAP2K4 | PIK3CA | F5 | APOB | ANKRD30A | MUC17 | DNAH17 | LAMA2 | ACE | CSMD2 | ||
STARD9 | USH2A | FAT3 | LPR2 | CSMD3 | MYO18B | DNAH5 | MDN1 | ARHGAP5 | DNAH9 | ||
CXCR3 | MUC16 | RB1 | PKHD1L1 | DNAH2 | SYNE2 | DYNC1H1 | PCLO | CACNA1B | ERBB2 | ||
PLK1 | SYNE1 | LYST | PTEN | SPTA1 | Untreated 2 | VHL | RYR1 | COL6A3 | UBB | ||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
Non-targeting 1 | ENAM | RYR2 | BRCA2 | Untreated 1 | FMN2 | HECW1 | LAMB4 | SI | Non-targeting 2 | ||
FHOD3 | MACF1 | RYR3 | C2ORF16 | DMD | FRG1 | HERC2 | MYH11 | STAB1 | ZDBF2 | ||
GOLGA6L2 | NEB | SMG1 | CACNA1E | DNA14 | GCC2 | HIVEP2 | NIPBL | TANC1 | ZNF536 | ||
HMCN1 | NF1 | UBR5 | CACNA1F | DYNC2H1 | GON4L | HYDIN | PKD1L1 | TF | ANK3 | ||
HRNR | OBSCN | USP34 | CYMA5 | FAM208B | GPR112 | ITSN2 | RNF213 | TPR | ASPM | ||
PLK1 | PCDH15 | XIRP2 | COL7A1 | FLG2 | Untreated 2 | VHL | SAGE1 | UNC80 | UBB | ||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
Non-targeting 1 | DCHS2 | MUC5B | ZFHX4 | Untreated 1 | MYO9A | SPHKAP | CXORF22 | NCOR1 | VPS13D | ||
ATR | DMXL2 | MXRA5 | ANK2 | KIAA1210 | PRUNE2 | TCHH | DANH6 | NOTCH2 | ANKRD12 | ||
BIRC6 | DNAH10 | TENM1 | ATM | LRP1 | SCN10A | VPS13C | DNAH7 | SPEN | C5ORF42 | ||
CDH1 | DNAH3 | PEG3 | DIDO1 | MAP1A | SCN2A | IPTR3 | ERBB3 | SRRM2 | CCDC88A | ||
CUBN | NOCK11 | RELN | DNAH8 | MED12 | SHROOM2 | CEP350 | FAT4 | SZT2 | CHD4 | ||
PLK1 | EYS | SACS | KIAA1109 | MED13 | Untreated 2 | VHL | KMT2A | VPS13A | UBB | ||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
Non-targeting | QSER1 | ARID1A | WDFY3 | Untreated 1 | SDK1 | TEX15 | LAMA1 | Non-targeting 2 | |||
COL14A1 | SHROOM3 | ATRX | EFCAB5 | SF3B1 | CBFB | AHNAK2 | |||||
CSMD1 | TBX3 | KIAA0947 | FOXA1 | ITPR1 | DDX3X | KIF4A | |||||
MEFV | UBR4 | MYCBP2 | INPPL1 | FLG | HECTD4 | FAT2 | |||||
MGAM | VCAN | NBEAL1 | MAP3K1 | AKAP9 | GPR98 | FOXO3 | |||||
PLK1 | ZNF462 | SETX | NRP1 | HERC1 | Untreated 2 | VHL | UBB | ||||
Table 1: Plate Layout and Human siRNA Library De-convolution. The table contains the content of each of the siRNA pools and layout for the low attachment plates used in the screen.
Three-dimensional models of cancer are being increasingly employed to evaluate the efficacy of known and novel compounds that have been designed to selectively kill cancer cells. Cancer cell spheroids are structures that display conditions more similar to those encountered in tumors in vivo, thus compounds that have increased efficacy in 3D are more likely to have an effect in vivo. However, these modalities do not allow for the identification of potentially novel targets that have not been the subject of drug design that could have considerable efficacy in treating cancer.
We developed a siRNA functional genomics approach that allowed for durable gene silencing for up to seven days in cancer cell line spheroids. There are several critical steps to the protocol that require optimization before a siRNA screen can be performed. The ability to form reproducible viable spheroids on a large scale is essential. Moreover, appropriate transfection conditions should be rigorously optimized. We suggest trialing several different transfection reagents with appropriate non-targeting and killing controls prior to attempting the screen. We were able to show that several commonly used breast cancer cell lines, namely BT474, MCF10DCIS.com, MDA-MB-231 and JIMT1 were amenable to siRNA transfection. Furthermore, we provide proof-of-principle data screening the 200 most frequently mutated genes in breast cancer in BT474 spheroids, confirmed their dependency on amplification of HER2 and oncogenic mutation of PIK3CA. Interestingly, silencing of the transcription factor FOXO3 resulted in a reduction in spheroid size but no significant effect on viability. FOXO3 is known to regulate the response to hypoxia, altering the metabolic capacity of cancer cells allowing them to adapt more readily to their environment 16. This role could potentially interfere with the cell viability reading as it detects ATP abundance, one of the main products of cell metabolism.
In support of the observation of a reduction in spheroid size, it has been shown that silencing of FOXO3 in HeLa xenografts impaired tumor growth and induced apoptosis 17. It is important to note that some genes may impact the capacity of cancer cells to retain their 3D architecture, which could result in false positives. For example, knockdown of E-Cadherin resulted in dissolution of BT474 spheroid structure. This had been previously reported using E-Cadherin targeted antibodies 18. As with any screening platform, potential targets should be rescreened to assess the reproducibility of the effect observed. There are limitations to the technique, namely the transient nature of the siRNA-mediated gene knockdown. Sustained silencing longer than seven days was not achievable with siRNA.
The advantage of this approach is that it can be coupled with various other biometric dyes not just those that assess spheroid viability, for example, giving spatial information of spheroid hypoxia or monitoring cells undergoing apoptosis. Moreover, because the plate reader scans are relatively quick and non-invasive, the impact of siRNAs on spheroid size can be assessed over time rather than just at the experimental end point. Indeed, we are currently exploring several of these avenues within our screening pipeline. An alternative approach that utilizes 3D cultures to identify novel dependencies is the use of chemical libraries that inhibit either a broad range of targets or particular families of proteins. Indeed, Bitler et al. utilized this targeted approach to identify the synthetic lethality interaction between ARID1A status and EZH2 inhibitors in ovarian clear cell carcinomas 19. The discovery of CRISPR-Cas9 gene editing technology has also allowed for the development of genetic screens in organoid cultures and in vivo. However, this approach is dependent on having appropriate animal facilities and may be cost prohibitive 20.
In conclusion, we believe that we have outlined a protocol that more accurately models the oxygen and nutrient gradients, which are features of the tumor microenvironment in vivo, allowing for the identification of novel cancer targets or robust validation of established targets. Moreover, our protocol can be applied to any type of cell line that forms spheroids and hence can be routinely used in the cancer research community for high-throughput siRNA screens.
The authors have nothing to disclose.
The authors would like to thank Fredrik Walberg (FACS and imaging facility, ICR) for assistance with microscopy. This work was funded by Breast Cancer Now. RN is the recipient of a Breast Cancer Now Career Development Fellowship (2011MaySF01).
Lullaby | Oz Biosciences | LL70500 | lipid-based transfection reagent |
Viromer | Lipocalyx | VB-01LB-01 | virus-like polymer transfection reagent |
Ultra-low attachment plate | Corning | CLS7007 | 96 well plate |
Foil plate seals | ThermoFisher | AB-0626 | |
Luminescent cell viability dye | Promega | G7570 | CellTitre-Glo |
Pipette tips (200ul) | Starlab | S1111-0806 | |
Pipette tips (10ul) | Starlab | S1111-3800 | |
Pipette tips (1000ul) | Starlab | S1122-1830 | |
Serological pipettes (5ml) | Sarstedt | 86.1253.025 | |
Serological pipettes (10ml) | Sarstedt | 86.1254.025 | |
Serological pipettes (25ml) | Sarstedt | 86.1685.020 | |
RPMI Media | GIBCO | 11875-093 | |
DMEM Media | GIBCO | 11965-084 | |
Opti-MEM | GIBCO | 31985070 | |
Feta bovine serum | GIBCO | 16140063 | |
siRNA | Dharmacon | Cherry picked library | |
Countess Cell Counter | ThermoFisher Scientific | AMQAX1000 | |
Cell counting chamber slides | ThermoFisher Scientific | C10312 | |
Celigo S | Nexcelom | contact company | |
Victor X5 | Perkin Elmer | contact company | |
Benchtop centrifuge | Various | ||
Axiovert Inverted brightfield microscope | Zeiss | contact company | |
Tissue culture C02 Incubator | Various | ||
Mulitichannel pipette | Various |