Abstract
Nico Jacobi1, Rita Seeboeck2, Elisabeth Hofmann2, Helmut Schweiger1, Veronika Smolinska1, Thomas Mohr3, Alexandra Boyer1, Wolfgang Sommergruber4, Peter Lechner5, Corina Pichler-Huebschmann5, Kamil Önder6,7, Harald Hundsberger2, Christoph Wiesner2 and Andreas Eger1,2
1IMC University of Applied Sciences Krems, Department Life Sciences, Research Institute for Applied Bioanalytics and Drug Development, A-3500 Krems, Austria
2IMC University of Applied Sciences Krems, Department Life Sciences, Institute of Medical and Pharmaceutical Biotechnology, A-3500 Krems, Austria
3Medical University of Vienna, Institute for Cancer Research, A-1090 Vienna, Austria
4Boehringer Ingelheim, RCV GmbH and Co KG, A-1121 Vienna, Austria
5University Clinic Tulln, Department Surgery, A-3430 Tulln, Austria
6Research Program for Rational Drug Design in Dermatology and Rheumatology, Department of Dermatology, Paracelsus Medical University of Salzburg, A-5020 Salzburg, Austria
7ProComCure Biotech, A-5081 Anif, Austria
Correspondence to:
Andreas Eger, email: [email protected]
Keywords: drug discovery; 3D models; ERBB signaling; targeted drugs; oncogene addiction
Received: August 25, 2017 Accepted: October 27, 2017 Published: November 17, 2017
ABSTRACT
Complex three-dimensional (3D) in vitro models that recapitulate human tumor biology are essential to understand the pathophysiology of the disease and to aid in the discovery of novel anti-cancer therapies. 3D organotypic cultures exhibit intercellular communication, nutrient and oxygen gradients, and cell polarity that is lacking in two-dimensional (2D) monolayer cultures. In the present study, we demonstrate that 2D and 3D cancer models exhibit different drug sensitivities towards both targeted inhibitors of EGFR signaling and broad acting cytotoxic agents. Changes in the kinase activities of ErbB family members and differential expression of apoptosis- and survival-associated genes before and after drug treatment may account for the differential drug sensitivities. Importantly, EGFR oncoprotein addiction was evident only in the 3D cultures mirroring the effect of EGFR inhibition in the clinic. Furthermore, targeted drug efficacy was strongly increased when incorporating cancer-associated fibroblasts into the 3D cultures. Taken together, we provide conclusive evidence that complex 3D cultures are more predictive of the clinical outcome than their 2D counterparts. In the future, 3D cultures will be instrumental for understanding the mode of action of drugs, identifying genotype-drug response relationships and developing patient-specific and personalized cancer treatments.
INTRODUCTION
Therapeutic efforts are gradually moving away from the administration of cytostatic and cytotoxic agents towards the use of rationally designed drugs that specifically inhibit the activity of oncoproteins [1, 2]. In case the oncoproteins are vital for the maintenance of the malignant phenotype, their targeted inhibition can lead to cell cycle arrest, senescence, differentiation, or apoptosis [3–5]. Gain of function mutations or overexpression of oncogenes have been described to account for oncoprotein addiction [6–8]. One major goal of current drug discovery is to identify novel genotype-drug response relationships and to understand the molecular basis of oncoprotein addiction. This knowledge is instrumental for predicting the efficacy of targeted drugs in individual patients [9, 10]. The overall aim is the optimization of clinical outcomes through the effective personalization of treatment [11, 12].
The number of potential anti-cancer agents has increased steadily over the past decade. However, the number of drugs that successfully passed clinical development is still rather low [13]. Lack of clinical efficacy and unacceptable toxicity are two main reasons for drug failure [14]. Considering the high costs of clinical trials it is of prime importance that compounds, which eventually will turn out as poor performers in the clinic, are eliminated early in drug development, even before animal testing or clinical trials have started. To achieve that goal, several improvements in preclinical drug testing have been implemented in the last years [15–18]. When combined with genomics and proteomics, the detailed molecular characterization of the mode of action and the toxicity of drugs is feasible [19–21]. However, the ability of the novel methods to produce physiologically relevant information inevitably depends on the nature of the in vitro cancer models and how closely these cell-based models reflect the pathophysiology of the disease [22–25]. At present, two-dimensional (2D) cancer cell cultures are commonly used for the screening and functional characterization of anti-cancer drugs [5, 26–28]. The cells are mostly cultured on synthetic substrates such as glass or plastic. The artificial environment causes alterations in cancer cell morphology, cell-cell and cell-matrix interactions, and consequently physiological changes in proliferation, differentiation and metabolism [15, 29–31]. In addition, it has been questioned whether conventional 2D cultures retain their dependency on the very same activated oncoproteins as the tumor cells do in vivo. The changes in cell physiology in the 2D models might alter oncoprotein activity and hence influence or mask oncoprotein addiction. Therefore, 2D cancer models are often only poor predictors of clinical drug efficacy and toxicity [5, 15, 32–34].
In the last years, a large number of technologies have been developed that facilitate the organotypic cultivation of cells in three dimensions (3D) [35–38]. The 3D cell models have become attractive tools for investigating cancer cell proliferation, metabolism, differentiation, tumor-stroma crosstalk, invasion and metastasis [16, 27, 39]. Examples are spheroid cultures, the embedding of cells in synthetic or natural extracellular matrices or the cultivation of cells in microporous materials [38, 40, 41]. They consist of single cell types (homotypic cultures), various combinations of different cell types (heterotypic cultures) or represent in vitro cultures of cancer tissue explants [23, 42–46]. Many studies have demonstrated that 3D organotypic models much better recapitulate the anatomy, physiology and drug sensitivity of human cancers in vivo [27, 47–58]. However, so far none of these studies provided conclusive evidence of whether 3D cultures are more predictive of the clinical outcome than their 2D counterparts [33, 47, 59, 60].
In the present study, we used different 2D and 3D cancer cell cultures and evaluated drug efficacies, oncoprotein addiction and cell survival in response to targeted drugs interfering with epidermal growth factor receptor (EGFR) signaling. We could show that oncoprotein addiction and drug efficacies in the 3D but not in the 2D cancer models were comparable to drug responses in cancer patients. EGFR inhibition caused massive cancer cell death in 3D models expressing mutated EGFR, whereas all 2D cultures showed only highly attenuated responses. We could identify major changes in gene expression and HER family kinase activities in the 3D cultures that contribute to the differential drug sensitivities.
RESULTS
Over the last years, the technologies for cultivating cancer cells in 3D have improved significantly [22, 23, 61–70]. However, it still remains highly controversial if the 3D cancer models are superior to conventional 2D cell cultures in terms of predicting clinical drug efficacy [59, 60]. In the present study we have generated 2D and 3D lung cancer models and analyzed their sensitivities towards broad-spectrum cytotoxic agents or targeted inhibitors of the EGFR pathway. The cancer cells express either wild type EGFR or harbor distinct EGFR mutations in exons 19, 20 and 21 that are known to influence targeted drug sensitivity in vivo [71–73]. We quantitatively assessed drug efficacies in the genetically and phenotypically different 2D and 3D cancer models. The major goal was to determine whether genetic alterations of the EGFR influence drug efficacy differently in the 2D and 3D cultures and whether EGFR oncoprotein addiction in vivo can be recapitulated in any of the in vitro cultures.
EGFR status, cell morphology and invasive potential of lung cancer cells
For the initial drug testing, we have selected four different lung cancer cell lines derived from non-small cell lung carcinoma (NSCLC) patients (Figure 1). HCC827 cells harbor an in-frame deletion in exon 19 (E746_A750). These and similar mutations have been demonstrated to be essential for responding to anilinoquinazoline inhibitors such as gefitinib in the clinic [71, 73]. It is assumed that this therapeutic effect is based on the fact that the small deletions cause a repositioning of critical residues in the vicinity of the ATP-binding cleft. The conformational change stabilizes the interaction of these residues with both ATP and its competitive inhibitors [71, 73]. NCI-H1975 cells exhibit two different point mutations T790M and L858R in exons 20 and 21 respectively [74, 75]. The E746_A750 and the L858R mutations, often in combination with EGFR overexpression, are known to hyperactivate EGFR signaling upon EGF stimulation [71, 72]. The T790M mutation has been suggested to confer resistance to targeted EGFR kinase inhibitors [76, 77]. Both cell lines express E-cadherin and are able to form compact 3D microtumors (Figure 1). NCI-H1975 cells are moderately invasive when embedded in extracellular matrix (ECM) (Figure 1, organoids). On the other hand, the NCI-H1437 and Calu-1 cells express wild type EGFR and form compact or loose 3D microtumors, respectively. Calu-1 cells have lost E-cadherin expression and exhibit a mesenchymal cell morphology. Consequently, the cells are highly invasive when embedded in ECM (Figure 1, organoids).
Figure 1: Phenotype of lung cancer cells cultured in 2D and 3D. Lung cancer cells (HCC827, NCI-H1975, NCI-H1437 and Calu-1) were cultured for five days on conventional plastic dishes (monolayer) or for five days as three-dimensional (3D) cell aggregates (spheroids) using the hanging drop method. For embedding, spheroids were generated using the hanging drop method. After 24-48 hours the spheroids were embedded in Matrigel/collagen type I matrices and further cultivated for five days (organoids). Cells in 2D were stained with antibodies for the epithelial cell-cell adhesion molecule E-cadherin and the intermediate filament protein vimentin and processed for immunofluorescence confocal microscopy. Note that Calu-1 cells have lost the expression of E-cadherin and strongly expressed the mesenchymal marker protein vimentin. They exhibited an invasive behavior when the spheroids were embedded in Matrigel/collagen type I. Vimentin was also expressed in HCC827, NCI-H1975 and NCI-H1437 cells, in particular when seeded at low density. Scale bars: monolayer phase-contrast, 200 μm; monolayer E-cadherin, 75 μm; phase-contrast spheroid and organoid, 100 μm.
Efficacy of targeted drugs in 2D and 3D cancer cell models
First, we determined whether the 2D and 3D cultures display differences in their sensitivity towards targeted drugs interfering with EGFR kinase activity. For that purpose, we cultivated the four lung cancer cell lines on conventional plastic dishes in 2D or generated 3D microtumors (spheroids) that were either freely floating or embedded in ECM consisting of Matrigel and collagen type I (here referred to as organoids) (Figure 1). The cultures were treated with different concentrations of the EGFR kinase inhibitors gefitinib (Iressa®, AstraZeneca) or erlotinib (Tarceva®, Roche) for three days. Subsequently, cell proliferation and cytotoxicity were assesed using the alamarBlue® cell viability assay [78]. In HCC827 cells overexpressing the EGFR mutant E746_A750 we observed striking differences in cell viability between 2D and 3D. EGFR inhibition caused a strong decrease in cell viability in the 3D cultures (survival rate around five to ten percent) whereas cell survival in the 2D cultures was only marginally affected (survival rate around sixty to seventy percent) (Figure 2A and 2B). In additional control experiments, we generated 2D cultures with different cell densities that ranged from 5,000 to 15,000 cells in a 96-well plate. We could not detect a significant decrease in cell viability upon gefitinib or erlotinib treatment in the 2D cultures, irrespective of the cell density (Supplementary Figure 1).
Figure 2: Efficacy of gefitinib and erlotinib in 2D and 3D cultures. Lung cancer cells (HCC827, NCI-H1975, NCI-H1437 and Calu-1) were treated with different concentrations of gefitinib (A) or erlotinib (B) for 72 hours. Cell viability was assessed using the alamarBlue® assay according to the manufacturer’s instructions for 6 hours at 37°C. Each data point represents the mean ± s.d. of six independent measurements. Cell viability was significantly reduced in the HCC827 and NCI-H1975 cancer cells cultured in 3D. On the contrary, the 2D cultures showed only weak or no responses to treatment.
Next, we treated the HCC827 2D and 3D cultures with varying concentrations of gefitinib or erlotinib, but now used a cell permeable fluorogenic protease substrate (glycylphenylalanyl-aminofluorocoumarin, GF-AFC) that allows to selectively detect protease activity in viable cells. In line with the alamarBlue® assay, only the 3D cultures exhibited a significant reduction in cell viability after EGFR inhibition (Figure 3, AFC life stain). However, despite the large differences in cell surivival rates, the EC50 of the EGFR inhibitors was similar in the 2D and 3D cultures ranging from 2–96 nM (Figure 2A and 2B). Interestingly, rather low concentrations of gefitinib and erlotinib were sufficient to activate caspases in the HCC827 3D cultures, demonstrating that the decrease in cell viability ultimately resulted in apoptotic cell death (Figure 3, Caspase-Glo® 3/7).
Figure 3: EGFR inhibitors induce cell death in 3D cultures. HCC827 cells were treated with different concentrations of gefitinib (A) or erlotinib (B). Cell viability, relative to untreated controls, was measured after 72 hours using MultiTox-Fluor GF-AFC live stain multiplexed with the Caspase-Glo® 3/7 assay according to manufacturer’s instructions. Each data point represents the mean ± s.d. of at least four independent measurements. NLU: normalized luminescence units.
In contrast to the HCC827 cells, high concentrations of the EGFR inhibitors were necessary to decrease cell viability in the NCI-H1975 3D cultures (Figure 2A and 2B). In addition to the EGFR activating mutation L858R, the NCI-H1975 cells have acquired the T790M drug resistance mutation, and hence higher drug concentrations were necessary to overcome drug resistance in the 3D cultures [76, 79]. The NCI-H1975 2D cultures did not show a decrease in cell viability, even at the highest drug concentrations (Figure 2A and 2B). On the other hand, EGFR inhibition did not affect the viability of the NCI-H1437 and Calu-1 cells, irrespective of whether the cells were cultured in 2D or in 3D (Figure 2A and 2B). Despite the fact that both cell lines express significant levels of the wild type EGFR protein (Supplementary Figure 2) the drugs were non-effective (Figure 2A and 2B).
In addition, we also exposed the four cell lines to the MEK1/MEK2 inhibitor trametinib (Mekinist®, GlaxoSmithKline). The MEK1/MEK2 inactivation decreased the survival rate of the HCC827 cells in the 3D cultures, whereas no effects were observed in the 2D configuration (Supplementary Figure 3).
To get further evidence that lung cancer cells cultured in 3D react differently to EGFR kinase inhibitors, we treated three further lung cancer cell lines expressing mutated EGFR variants with gefitinib and erlotinib [71]. The cell line HCC4006 is harboring a similar exon 19 deletion (L747_A750) as the HCC827. Gefitinib and erlotinib significantly decreased cell viability in the 3D but not in the 2D cultures (Figure 4). The cell lines HCC2935 and NCI-H1650 exhibit the EGFR exon 19 mutations E746_S752del and E746_A750del respectively. However, both cell lines did not respond to drug treatment, neither in the 2D nor in the 3D setup (Figure 4). At present, the reason for drug failure in the two cell lines is unclear. HCC2935 cells also exhibit a mutation at S752I and NCI-H1650 cells contain an additional mutation at A750E [80]. It needs to be demonstrated whether these mutations are responsible for escape from treatment with gefitinib and erlotinib. Also gender specific tumor physiology might influence drug sensitivity [81]. In the future, 3D models will be instrumental for deciphering the precise molecular rationale for drug failure, and hence improve our understanding of personalized medicine. Taken together, we could demonstrate that the treatment of cancer cells with inhibitors targeting EGFR signaling can have different outcomes. Drug efficacy was dependent on both the genetic makeup and the cultivation conditions of the cancer cells. Only 3D cancer models exhibited a significant sensitivity towards inhibition of EGFR signaling, whereas the 2D cultures were only weakly responding to the drugs.
Figure 4: 3D cultures are instrumental for assessing oncogene addiction. HCC4006, HCC935 and NCI-H1650 were treated with different concentrations of gefitinib and erlotinib for 72 hours. Cell viability, relative to untreated controls, was determined using the alamarBlue® assay according to the manufacturer’s instructions. Each data point represents the mean ± s.d. of at least six independent measurements. HCC4006 cells exhibited strong sensitivity toward both EGFR inhibitors when cultivated in 3D.
Efficacy of broad-acting cytotoxic drugs in 2D and 3D cancer cell models
Next we tested the efficacy of chemotherapeutics in the different lung cancer models. Paclitaxel and cisplatin are broad-acting cytotoxic drugs that are commonly administered to NSCLC patients in the clinic [82–84]. All cell lines were rather insensitive to the genotoxic drug cisplatin and high concentrations of the drug were necessary to decrease the viability after three days of treatment (Supplementary Figure 4). We did not observe any differences in cisplatin sensitivity between the 2D and 3D cultures (Supplementary Figure 4). However, when treating the cells with the microtubule stabilizing agent paclitaxel, the cell lines reacted differently. High doses of paclitaxel were toxic for NCI-H1975 and Calu-1 cells cultured in 2D. The same cells were rather insensitive to high paclitaxel concentrations after emdedding the spheroids in ECM (Supplementary Figure 5, organoids). Interestingly, these two cell lines displayed a noticeable invasive behavior in ECM (see Figure 1). On the other hand, paclitaxel had weak or no effects in HCC827 and NCI-H1437 cells respectively, regardless of the culture condition.
Proliferation rate and cell size in 3D cultures are similar to cancer cells in tumor tissue
Cell proliferation and architecture are known predictors of drug sensitivity of cancer cells. Therefore, we assessed proliferation and cell size in 2D monolayers and in spheroids. Interestingly, the cells in spheroids exhibited an extraordinary low proliferation rate (Figure 5A). The mean doubling times (DT) of the four NSCLC cell lines grown as monolayers were around 1.5–3 days whereas the DT of the cells in the spheroids increased to 8–22 days. The low proliferation rate of cancer cells in spheroids is much closer to the tumor growth rate in breast, colon and lung tumors in vivo [85–88]. In addition, we observed a 2.6 to 3.2-fold reduction in cell size when cultivating the cancer cells as spheroids (Figure 5B), which is in line with reports published previously [89, 90]. Taken together, the high proliferation rate in 2D did not sensitize the cells to EGFR inhibitors.
Figure 5: Proliferation and cell size in 2D and 3D cultures. (A) 3000 cells were seeded to generate the 2D and 3D cultures. Every second day the cells were harvested by trypsinization for determining cell number and viability. Each data point represents the mean ± s.d. of three independent experiments. (B) Cell size distributions were assessed by measuring the diameter of >100 living adherent cells cultured on plastic dishes (2D) or >30 cells of H&E stained sections of paraformaldehyde fixed and paraffin embedded (FFPE) spheroids. Phase-contrast microscopy was performed using a Leica DMI6000B inverted microscope equipped with CTR6500 microscope drive control, DFC420C digital microscope camera with a 5 Megapixel CCD sensor and Leica Application Suite Version 3.8.0 software.
Expression and phosphorylation of ErbB proteins are altered in the 3D cultures
Differential drug sensitivities in the 3D models might be due to changes in the kinase activities of ErbB family members. Hence, we assessed ErbB protein and phosphorylation in 2D and 3D cultures via immunoblotting. We found that the total protein levels of EGFR and ErbB2 were significantly reduced in the 3D cultures (Figure 6A). Importantly, we also observed a marked downregulation of tyrosine phosphorylation when normalizing the phosphotyrosine levels to total intracellular EGFR or ErbB2 protein (Figure 6B). On the other hand, ErbB3 protein expression was significantly increased in the 3D cultures (Figure 6A). However, when normalizing the ErbB3 phosphotyrosine levels to total ErbB3 protein, we found that the phosphorylation of ErbB3 receptors in 3D was lower than in 2D (Figure 6B). In summary, protein expression and kinase activities of ErbB family members were substantially altered in the 3D cancer models which might at least partly explain the differential efficacies of the targeted inhibitors gefitinib, erlotinib and trametinib.
Figure 6: Expression and phosphorylation of ErbB family members in 2D and 3D cultures. (A) Immunoblots and quantification of total EGFR, ErbB2 and ErbB3 protein levels in HCC827 and HCC4006 cells cultured as monolayers or as spheroids. Total ErbB protein levels were normalized against actin. (B) Phosphorylation of the intracellular domain of EGFR, ErbB2 and ErbB3 on different tyrosines. Phosphorylation levels were normalized against the respective total ErbB protein, which in turn was not normalized to actin. Bars represent the mean ± s.d. of at least three independent experiments (except for p-ErbB3, Y1197).
2D and 3D cancer models exhibit differences in gene expression
To assess global gene expression in 2D versus 3D cultures we performed Affymetrix Gene Chip analysis using total RNA isolated from HCC827 cells. Microarray data can be accessed under Gene Expression Omnibus Accession Number (GSE102722). In the present study, we focused on the expression of apoptosis-associated genes, as they might have an impact on the differential drug efficacies in the 2D and 3D cultures. Affymetrix Gene Chip data were imported into R to be further processed using oligo and limma packages of R. As depicted in Figure 7A, the expression of a large array of genes involved in apoptosis was substantially increased in the 3D spheroids of HCC827 cells. A comprehensive list of genes annotated with the GO Term apoptotic process (GO:006915) can be found at http://amigo.geneontology.org/amigo/term/GO:0006915. Using the Ingenuity Pathway Analysis (IPA) software package we could identify major changes in gene expression in the TNF/FAS and growth factor signaling cascades (Figure 7B). The expression rate of MAP-kinases was mostly reduced, whereas the gene expression of signaling interfaces of the extrinsic and intrinsic apoptotic program was upregulated (Figure 7B). Major players in apoptosis that were differentially expressed in the 2D and 3D cultures include TNFR1, BFL-1, BAX, BAK, DIABLO, BCL2, NFkB, CASP3 and CASP7 (Figure 7B).
Figure 7: Expression of apoptosis- and survival-associated genes in 2D and 3D cultures. (A) Affymetrix Gene Chip data were imported into R and processed using oligo and limma packages of R. Quality control was done using the R package arrayQualityMetrics. Eigengene calculation was carried out using the WGCNA with modules defined by gene ontology term Apoptosis - GO:0006915. A gene was considered to be part of the module if it was annotated directly with the respective GO-Term. Note that the expression of genes involved in apoptosis was substantially altered in the 3D spheroid cultures. (B) Pathway analysis of apoptosis-associated genes using the Apoptosis Pathway of QIAGEN’s Ingenuity Pathway Analysis©. (C) 2D and 3D cultures were treated for 24 hours with 5 nM gefitinib. RNA was isolated und processed for RT-qPCR using the Human Apoptosis RT² Profiler PCR Array. The transcription of several genes associated with cell survival was downregulated in the 3D cultures upon gefitinib treatment. Each data point represents the mean of two independent experiments.
Next, we investigated if and to which extent the expression of apoptosis-associated genes is altered in the 2D and 3D cultures after treatment with tyrosine kinase inhibitors (TKI). 2D and 3D cultures of HCC827 cells were treated with 50 nM gefitinib for 24 hours. Total RNA was isolated and subjected to RT-qPCR using the Human Apoptosis RT² Profiler PCR Array (PAHS-012Z, QIAGEN) and Rotor-Gene Q platform. Quality control and normalization was performed with the RT2 Profiler PCR Array Data Analysis Webportal. We could show that in response to gefitinib treatment the expression of several genes that favor survival of cells (e.g. BCL-2 family members) was significantly downregulated in the 3D cultures (Figure 7C). Reduced expression of anti-apoptotic genes in response to TKI treatment might have a strong impact on the drug sensitivity of the 3D cultures. In summary, 3D cultivation caused major changes in the expression of apoptosis- and survival-associated genes, both before and after treatment with the EGFR inhibitor gefitinib.
Heterotypic 3D cancer models are hypersensitive to EGFR inhibition
The molecular crosstalk between tumor and stromal cells is known to influence cancer cell survival and proliferation, and consequently alter the efficacy of anti-cancer drugs [91]. In the present study, we generated cocultures consisting of spheroids (HCC827 cells) and cancer-associated fibroblasts (CAFs) that were embedded in a natural ECM made of Matrigel and collagen type I. Cocultures were grown for 10 days and then treated with 50 nM gefitinib for 24 hours. Cell death was assessed using antibodies that specifically recognize the cleaved caspase-3 fragment (CC3). Gefitinib triggered substantial apoptosis in the HCC827 homotypic cultures (Figure 8, homotypic culture). This is consistent with experiments in which the alamarBlue®, AFC life cell stain or Caspase-Glo® 3/7 assays were used to determine survival and cell death (Figures 2 and 3). Strikingly, in the cocultures the amount of apoptotic cancer cells increased significantly upon gefitinib application (Figure 8, heterotypic cultures). Hence, the CAFs increased the drug sensitivity of the HCC827 spheroids in the Matrigel/collagen matrix. Furthermore, we found that CAFs can trigger an invasive behavior in NCI-H1975 cells (Supplementary Figure 6). Cancer cell invasion was abrogated when treating the cells with the low molecular weight inhibitor crizotinib (Xalkori®, Pfizer), that targets the RTKs ALK, MET and ROS (Supplementary Figure 6). However, crizotinib did not induce cancer cell death but rather inhibited the invasive behavior of the cancer cells. In the CAF population themselves no apoptosis was evident after gefitinib or crizotinib treatment (Figure 8 and Supplementary Figure 6, heterotypic cultures). The molecular mechanisms of how CAFs increased drug sensitivity or triggered cancer cell invasion are unclear at present. Nevertheless, it shows that the integration of stromal cells into 3D cancer models can be a critical issue for properly evaluating the potency and the mode of action of anti-cancer drugs.
Figure 8: Cancer-associated fibroblasts alter drug responses in vitro. Organoids were generated from the lung adenocarcinoma cell line HCC827, that was cultured in the absence or presence of cancer-associated lung fibroblasts. Homo- and heterotypic cultures of HCC827 cancer cells without or with cancer-associated fibroblasts respectively, were treated with 50 nM gefitinib for 24 hours, processed for immunofluorescence microscopy using cleaved caspase-3 (CC3) and vimentin antibodies and mounted using Vectashield mounting medium containing DAPI. Scale bars: 75 μm.
DISCUSSION
In the past, the preclinical testing of anti-cancer drugs was mostly done with conventional 2D cell cultures [5, 92–94]. However, data obtained from 2D cultures often led to misinterpretations concerning the efficacy and toxicity of drugs [74, 95–100]. Drugs that showed impressive potencies against 2D monolayers often failed in subsequent cost-intensive clinical trials [13, 14]. To overcome these drawbacks, different technologies have been developed that allow the long-term cultivation of 3D organotypic cultures in vitro [62, 65, 101]. Evidence has been provided that 3D cultures recreate the physiological and mechanical cues that are typical for the cancer microenvironment [34, 68, 102]. They contain heterogeneous zones of proliferating, quiescent, and dying cells, which are likewise found in human tumor tissue [27, 103–106]. Thus, 3D models are likely to bring significant value to preclinical drug testing and bridge the gap between traditional in vitro monolayer cultures and expensive in vivo animal studies [26, 32, 35, 37, 45, 107]. However, to date the number of studies that quantitatively assessed drug efficacy in 3D models is still rather low. Further research is required to validate the applicability and the advantages of 3D models for drug testing.
Here we show that only 3D cultures (spheroids, organoids and heterotypic cultures) expressing mutated versions of the EGFR are sensitive to treatment with the low molecular weight drugs gefitinib, erlotinib and trametinib (Figures 2, 3, 4 and Supplementary Figure 3). On the contrary, the same cancer cells responded only weakly to the inhibitors when cultivated as 2D monolayers. The drug effects observed in 3D directly correlate with observations in the clinic. Lung cancer patients whose tumors express mutated EGFR often benefit from drug exposure. In response to treatment, the cancer cells undergo apoptosis and the tumors shrink extensively [108, 109]. Hence, our data provide conclusive evidence that 3D cancer models can predict clinical outcome with high accuracy. In line with our findings, recent reports also demonstrated differential drug efficacies of targeted inhibitors of ErbB proteins in 2D and 3D cultures [47, 110–113]. Eimer et al. demonstrated that three adherent cell lines derived from glioblastomas responded only moderately to erlotinib. Conversely, erlotinib induced a strong cell growth inhibition when the cells were cultivated as 3D neurospheres [111]. Treatment with trastuzumab induced a slight growth reduction in SKBR-3 cells cultured in 2D, whereas proliferation was strongly inhibited in the 3D spheroids mirroring the clinical benefit of trastuzumab for ErbB2-positive breast cancer patients [47]. Furthermore, in a panel of different breast cancer cell lines, the efficacies of the ErbB2-targeting drugs trastuzumab, pertuzumab and lapatinib were highly dependent on whether the cells were cultured as 2D monolayers or embdedded in 3D laminin-rich ECM gels [110]. Recently, H292 lung cancer cells grown in 2D and 3D were used to screen a compound library consisting of 41 anti-cancer agents [113]. The screening revealed pronounced differences in the potencies of the tested drug candidates in 2D and 3D. Interestingly, in some cases the 3D cultures appeared less and in other cases more sensitive to certain compounds [113].
Several studies provided evidence for a decreased potency of broad acting chemotherapeutics in 3D compared to 2D [113–115]. In our study, we found that two lung cancer cell lines were rather insensitive to high concentrations of paclitaxel when cultured as spheroids embedded in extracellular matrix (Supplementary Figure 5, organoids). Interestingly, the paclitaxel sensitive organoid cultures showed invasive behavior in long-term cultures (Figure 1). Therefore, cells cultivated in 3D might accurately predict the efficacies of drugs targeting the cytoskeleton. However, not all cancer cells grown in 3D are refractory to chemotherapeutics. They often show the same responses in 2D and 3D [113, 114, 116]. Also in our study, we did not observe any differences in the growth of 2D and 3D cultures that were treated with different cisplatin concentrations (Supplementary Figure 4).
Differences in gene expression, signal perception and integration, and intercellular communication might account for the pronounced differences in drug response between the 2D and 3D cancer models. We could observe major changes in the overall expression and phosphorylation of ErbB family members (Figure 6). In the 3D cultures, EGFR and ErbB2 expression was downregulated, whereas ErbB3 protein levels were increased. This was accompanied by a reduced phosphorylation of ErbB1, ErbB2, and ErbB3 proteins in the 3D cultures. Phosphorylation was reduced at multiple sites within the intracellular domain of the ErbB proteins (Figure 6). The low proliferation rate in spheroids (Figure 5A) might be a direct consequence of the reduced signaling activities of ErbB family members. A negative impact of TNF alpha signaling on EGFR activity has been reported recently [117]. We found a transcriptional upregulation of genes participating in TNF signaling (Figure 7B). Hence, reduced EGFR phosphorylation and signaling in 3D might be due to increased gene expression and activation of members of the tumor necrosis factor (TNF) alpha pathway (Figure 7B). In line with our findings, Ekert et al. found a downregulation of EGFR and c-MET proteins in lung cancer 3D cultures [112]. In addition, upon stimulation with ligands the extent of phosphorylation of EGFR and c-MET was lower in the 3D cultures [112]. Similarly, alterations in ErbB3 expression and phosphorylation were previously shown to modulate cell survival and growth in a 3D microenvironment [63].
Importantly, we detected major changes in gene expression when culturing the cancer cells in 3D. Differences in gene expression were particularly evident when focusing on critical regulators of apoptosis. The transcription of genes encoding for members of the TNF receptor family and different regulators of extrinsic and intrinsic apoptosis was upregulated in the 3D cultures (Figure 7B). In line with our findings, numerous previous studies have shown that 3D cultivation can significantly alter gene expression in cancer cells [118–120]. After treatment with gefitinib, the 3D cultures downregulated the transcription of genes critically involved in the regulation of cell survival (Figure 7C, BCL2 family members). Inhibition of EGFR activity might result in reduced AKT activity and thereby decreasing BCL2 expression. Such a link has recently been demonstrated in several studies [121–123]. Taken together, TKI treatment severely disturbs the balance between pro- and anti-apoptotic molecules in the 3D cultures, which eventually leads to apoptosis.
The cellular composition of the tumor stroma can significantly influence the outcome of drug treatments [91]. Gene expression and cellular phenotype are influenced by the surrounding microenvironment consisting of extracellular matrix (ECM) and different tumor stromal cells [102, 124–126]. Data from laser-captured microdissection and subsequent whole-genome Affymetrix GeneChip analyses revealed that 3D coculturing of colon cancer cells with cancer-associated fibroblasts alters gene expression in both stromal and cancer cells [102, 126]. Here, we found that the integration of cancer-associated fibroblasts into 3D cultures substantially increased the sensitivity of cancer cells towards the EGFR kinase inhibitor gefitinib (Figure 8). Interestingly, Wang et al. demonstrated recently that lung cancer cells became resistant to EGFR-TKI when cocultured in vitro with fibroblasts [127]. However, in that study coculturing was done in Transwell chambers where the cancer cells and fibroblasts were separated by the 8 μm pore filters [127]. In our opinion, the coculturing of spheroids and fibroblasts in extracellular matrix is more close to the in vivo setting. This has already been demonstrated in 3D cultures of breast, lung and colon cancer cells [39, 102, 124, 126]. Hence, the development of customized heterotypic 3D cancer models that incorporate stromal cells and paracrine signaling is instrumental for understanding the mode of action of drugs and for predicting clinical drug efficacies. Our findings are in line with other studies that demonstrated the impact of the tumor microenvironment on tumor cell physiology [91]. Also the mutational status of the tumor donor can influence the impact of CAFs on drug sensitivity. In specific cases CAFs can contribute to epithelial to mesenchymal transition (EMT) in tumor cells and induce resistance to low molecular weight TKIs [128, 129].
Taken together, we could show that 3D models facilitate a physiologically relevant evaluation of candidate drugs in vitro. A number of recent studies provided evidence that the generation of 3D cancer models can be automatized and used for high-throughput phenotype-based drug discovery [17, 18, 26, 27, 45, 116, 130, 131]. In the future, the 3D disease models will be beneficial for the drug approval process, increase cost-effectiveness and reduce the number of animal experiments. By means of the novel technologies different critical issues in drug discovery can be addressed including the mode of action of compounds, (non-) oncogene addiction, genotype-drug response relationships and synthetic lethality. The 3D cancer models will be instrumental for the identification of patient-specific and personalized cancer treatments and for understanding the molecular basis of drug resistance.
MATERIALS AND METHODS
Tumor cell lines and primary cells
HCC827 (female, 39 year), NCI-H1975 (female), NCI-H1437 (male, 60 years), HCC4006 (male, 50+ years), HCC2935 (male, 39 years), NCI-H1650 (male, 27 years) cells were originally derived from lung adenocarcinomas and Calu-1 (male, 47 years) from lung squamous cell carcinoma. All cell lines were obtained from the American Type Culture Collection (ATCC; Manassas, VA, USA). Cells were maintained in a humidified atmosphere at 37°C in 5% CO2 in RPMI-1640 or McCoy’s 5A-GlutaMAX supplemented with 10–15% fetal bovine serum, 2 mmol/L glutamine and antibiotics (Thermo Fisher Scientific, Madison, WI, USA). Fibroblasts were isolated from lung cancer tissue immediately after surgery (staging and grading: invasive, middle to low grade squamous cell carcinoma, G2-3 WHO pT2a G3 pL0 pN0 (0/22) pV0; mutational status of cancer cells unknown). The cancer samples were collected in accordance with the guidelines of the institutional ethics committee. Cylindrical tissue cores were generated using a motorized tissue coring press (Alabama Research and Development, Munford,TN, USA). Living tissue slices of 100–300 μm thickness and 5–12 mm diameter were cut from the tissue cores with the semi-automated precision-cut vibrating-blade microtome VT1200S (Leica, Wetzlar, Germany). The tissue slices were maintained at 37°C in 5% CO2 in RPMI-1640 medium supplemented with 20% fetal bovine serum (FBS), 2 mmol/L glutamine, antibiotics, 250 ng/mL amphotericin B, 10 μg/mL insulin, 10 μg/mL transferrin and 10 μg/mL selenious acid. After 5–15 days the tissue slices were carefully removed and the adherent cultures of fibroblasts were maintained at 37°C in 5% CO2 in serum-free fibroblast growth medium (PromoCell, Heidelberg, Germany).
The lentiviral system LentiORF™ was purchased from Thermo Fisher Scientific (Thermo Fisher Scientific, Madison, WI, USA). Primary cancer-associated fibroblasts (CAFs) and NCI-H1975 cancer cells were visualized by lentiviral-based stable expression of monomeric cyan fluorescent protein (mCFP) or Emerald green fluorescent protein (EGFP) according to manufacturer’s instructions by using the TLA-HEK293T packaging cell line, the transfection reagent Express-In™, control vectors and ORF plasmid DNA.
Spheroid generation and embedding
Spheroids were generated using the hanging drop method (GravityPLUS™ microtissue culture system; InSphero AG, Zurich, Switzerland) or using HydroCell™ 96-well flat-bottom plates (Nunc, Roskilde, Denmark) according to the manufacturer’s instructions. To generate organotypic cultures the spheroids were embedded in extracellular matrix (ECM) composed of equal amounts of growth medium, neutralized rat tail collagen type I and Matrigel® (Corning, NY, USA). For heterotypic cultures, spheroids consisting of 1.5 × 104 cancer cells were embedded together with 7.5 × 104 primary lung cancer fibroblasts (CAFs).
Immunofluorescence microscopy
Cells were fixed in 4% formaldehyde (Electron Microscopy Sciences, Hatfield, PA, USA), permeabilized in 0.2% Triton X-100 and incubated in 50 mM NH4Cl/0.1% glycine. Primary antibodies (E-cadherin, BD Biosciences, #610182; cleaved caspase-3, Cell Signaling Technology, #9661; vimentin, Dako, #M0725) were incubated over night at 4°C and secondary antibodies for 4 hours at room temperature. Samples were mounted using Vectashield mounting medium (H-1500, Vector Laboratories, Burlingame, CA, USA) containing 1.5 μg/mL 4′,6-diamidino-2-phenylindole (DAPI) and analyzed in a Leica DMIRE2 microscope equipped with a TCS SP2 confocal unit (Leica Microsystems, Heerbrugg, Switzerland).
SDS-PAGE and immunoblotting
Cells were lysed in situ in hot Laemmli sample buffer [132]. Cell lysates were further incubated for 5 minutes at 95°C. Primary antibodies were obtained from Abcam (Cambridge, MA, USA): phospho-ErbB3 Y1222 (#ab133445) and Y1289 (#ab76469) or Cell Signaling Technology (Danvers, MA, USA): ErbB1 (#4267), ErbB2 (#2165), ErbB3 (#12708), phospho-ErbB1 Y845 (#6963), Y992 (#2235), Y1068 (#3777), Y1086 (#2220), Y1148 (#4404), phospho-ErbB2 Y1196 (#6942), phospho-ErbB3 Y1197 (#4561), Y1328 (#8017) and actin (#3700). Secondary antibodies were derived from LI-COR Biosciences (Lincoln, NE, USA): IRDye 800CW goat anti-mouse IgG (#926-32210) and IRDye 800CW goat anti-rabbit IgG (#926-32211). Immunoblots were analyzed and quantified using the Odyssey CLx infrared imaging system (LI-COR Biosciences, Lincoln, NE, USA).
Cell viability and apoptosis assays
Gefitinib (#G-4408), erlotinib (#E-4007), trametinib (#T-8123), crizotinib (#C-7900) and paclitaxel (#P-9600) were purchased from LC Laboratories (Woburn, MA, USA). Cisplatin (#479306) was purchased from Sigma Aldrich (St. Louis, MO, USA). 2D and 3D cultures were incubated in the presence of the drugs for 72 hours. Cell viability was determined using the commercially available alamarBlue® assay kit (Thermo Fisher Scientific, Madison, WI, USA). Apoptosis was analyzed using the Caspase-Glo® 3/7 assay multiplexed with the MultiTox-Fluor GF-AFC life stain (Promega, Madison, WI, USA). Luminescence and fluorescence intensities were measured using the Paradigm detection platform (Molecular Devices, Sunnyvale, CA, USA). Statistical analysis was done using the GraphPad Prism Software 7.03 (GraphPad Software Inc., La Jolla, CA, USA). For the in situ analyses of apoptosis in cocultures the cells were incubated in the presence of 50 nM gefitinib or 20 nM crizotinib for 24 hours and processed for immunofluorescent microscopy using the cleaved caspase-3 antibody (Cell Signaling Technology, #9661) that specifically detects apoptotic cells.
Affymetrix GeneChips and RT-qPCR arrays
Total RNA was extracted using the miRNeasy Mini Kit (QIAGEN, Venlo, NED) according to the manufacturer’s instructions. RNA quality was assessed using the capillary gel electrophoresis platform Experion™ and Experion RNA StdSens Analysis Kit (#700-7104) from Bio-Rad according to manufacturer’s instructions (Software Version 3.1). The extracted RNA was amplified and labeled with the MassageAmpII-Biotin Enhanced Kit (Ambion, Austin, TX, USA). Fragmented antisense-RNA (15 μg) was used for hybridization of the Human Genome HuExon st 1.0 (Affymetrix). The arrays were hybridized and scanned using standard Affymetrix protocols. Microarray data were normalized using the Robust Multi-Array Analysis as implemented in Bioconductor [133, 134]. All analyses were performed with log2-transformed data. Hypothesis tests were performed using a modified t statistics with an empirical Bayes approach as implemented in Bioconductor LIMMA package [135]. P-values were adjusted by the false discovery rate method of Benjamini and Hochberg [136].
For RT-qPCR, cells were treated with 50 nM gefitinib for 24 hours and subjected to RT-qPCR using the Human Apoptosis RT² Profiler PCR Array (#PAHS-012Z, QIAGEN, Venlo, NED). Expression of genes was investigated in biological duplicates using the Rotor-Gene Q (QIAGEN, Venlo, NED). Quality control and normalization was performed with software provided by the RT2 Profiler PCR Array Data Analysis Webportal.
Bioinformatics
Affymetrix GeneChip data were imported into R and processed using oligo and limma packages of R [137, 138]. Quality Control was done using the R package arrayQualityMetrics [139]. To assess functionality with regard to the GO term apototic process (GO:006915) eigenexpression levels were used, as defined by the Singular Value Decomposition (SVD) of the expressionmatrices of genes annotated with the respective GO Term. Briefly, a gene was considered to be annotated with a specific GO-Term if it was either annotated by the term itself or one of its descendants (implicit annotation). The expression of genes implicitly annotated with the respective GO-Term was summarized using Single Value Decomposition as described by Alter et al. [140] using the WGCNA package of R [141]. Ingenuity pathways were generated using QIAGEN‘s Ingenuity Pathway Analysis (IPA®).
ACKNOWLEDGMENTS AND FUNDING
This work was supported with funds received from the Forschungs- und Bildungsges.m.b.H. (NFB) of Lower Austria, project numbers LS11-013 and LS13-012.
CONFLICTS OF INTEREST
The authors declare no conflicts of interest.
REFERENCES
1. Alaoui-Jamali MA, Morand GB, da Silva SD. ErbB polymorphisms: insights and implications for response to targeted cancer therapeutics. Front Genet. 2015; 6:17.
2. Hutchinson KE, Johnson DB, Johnson AS, Sanchez V, Kuba M, Lu P, Chen X, Kelley MC, Wang Q, Zhao Z, Kris M, Berger MF, Sosman JA, et al. ERBB activation modulates sensitivity to MEK1/2 inhibition in a subset of driver-negative melanoma. Oncotarget. 2015; 6:22348–22360. http://doi.org/10.18632/oncotarget.4255.
3. Weinstein IB. Cancer. Addiction to oncogenes—the Achilles heal of cancer. Science. 2002; 297:63–64.
4. Weinstein IB, Joe A. Oncogene addiction. Cancer research. 2008; 68:3077–3080; discussion 3080.
5. Sharma SV, Haber DA, Settleman J. Cell line-based platforms to evaluate the therapeutic efficacy of candidate anticancer agents. Nature reviews Cancer. 2010; 10:241–253.
6. Chen HY, Liu CH, Chang YH, Yu SL, Ho BC, Hsu CP, Yang TY, Chen KC, Hsu KH, Tseng JS, Hsia JY, Chuang CY, Chang CS, et al. EGFR-activating mutations, DNA copy number abundance of ErbB family, and prognosis in lung adenocarcinoma. Oncotarget. 2016; 7:9017–9025. http://doi.org/10.18632/oncotarget.7029.
7. Lee VH. Association of Exon 19 and 21 EGFR Mutation Patterns with Treatment Outcome after First-Line Tyrosine Kinase Inhibitor in Metastatic Non–Small-Cell Lung Cancer. Journal of Thoracic Oncology. 2013; 8:1148–1155.
8. Sharma SV, Settleman J. ErbBs in lung cancer. Experimental cell research. 2009; 315:557–571.
9. Crews KR, Hicks JK, Pui CH, Relling MV, Evans WE. Pharmacogenomics and individualized medicine: translating science into practice. Clinical pharmacology and therapeutics. 2012; 92:467–475.
10. Ni X, Zhang W, Huang RS. Pharmacogenomics discovery and implementation in genome-wide association studies era. Wiley interdisciplinary reviews Systems biology and medicine. 2013; 5:1–9.
11. Puri A, Saif MW. Pharmacogenomics update in pancreatic cancer. JOP : Journal of the pancreas. 2014; 15:114–117.
12. Karim NA, Bui H, Pathrose P, Starnes S, Patil N, Shehata M, Mostafa A, Rao M, Zarzour A, Anderson M. The use of pharmacogenomics for selection of therapy in non-small-cell lung cancer. Clinical Medicine Insights Oncology. 2014; 8:139–144.
13. Waring MJ, Arrowsmith J, Leach AR, Leeson PD, Mandrell S, Owen RM, Pairaudeau G, Pennie WD, Pickett SD, Wang J, Wallace O, Weir A. An analysis of the attrition of drug candidates from four major pharmaceutical companies. Nature reviews Drug discovery. 2015; 14:475–486.
14. Kola I, Landis J. Can the pharmaceutical industry reduce attrition rates? Nature reviews Drug discovery. 2004; 3:711–715.
15. Cukierman E, Pankov R, Yamada KM. Cell interactions with three-dimensional matrices. Current opinion in cell biology. 2002; 14:633–639.
16. Bissell MJ, Radisky DC, Rizki A, Weaver VM, Petersen OW. The organizing principle: microenvironmental influences in the normal and malignant breast. Differentiation; research in biological diversity. 2002; 70:537–546.
17. LaBarbera D, Reid B, Yoo B. The multicellular tumor spheroid model for high-throughput cancer drug discovery_LaBarbera_Expert Opin Drug Discov. 2012; 7:819–830.
18. Reid BG, Jerjian T, Patel P, Zhou Q, Yoo BH, Kabos P, Sartorius CA, Labarbera DV. Live multicellular tumor spheroid models for high-content imaging and screening in cancer drug discovery. Current chemical genomics and translational medicine. 2014; 8:27–35.
19. Gao Y, Holland RD, Yu LR. Quantitative proteomics for drug toxicity. Brief Funct Genomic Proteomic. 2009; 8:158–166.
20. Harrill AH, Rusyn I. Systems biology and functional genomics approaches for the identification of cellular responses to drug toxicity. Expert opinion on drug metabolism & toxicology. 2008; 4:1379–1389.
21. Stoehlmacher J. The impact of genomics and proteomics in the clinic: functional genetic polymorphisms and their value in response and toxicity prediction in solid tumours. Ann Oncol. 2006; 17:x263–268.
22. Unger C, Kramer N, Walzl A, Scherzer M, Hengstschlager M, Dolznig H. Modeling human carcinomas: physiologically relevant 3D models to improve anti-cancer drug development. Advanced drug delivery reviews. 2014; 79–80:50–67.
23. Thoma CR, Zimmermann M, Agarkova I, Kelm JM, Krek W. 3D cell culture systems modeling tumor growth determinants in cancer target discovery. Advanced drug delivery reviews. 2014; 69–70:29–41.
24. Smalley KS, Lioni M, Herlyn M. Life isn’t flat_taking cancer biology to the next dimension. In Vitro Cellular & Developmental Biology Animal. 2006; 42:242–247.
25. Boj SF, Hwang CI, Baker LA, Chio II, Engle DD, Corbo V, Jager M, Ponz-Sarvise M, Tiriac H, Spector MS, Gracanin A, Oni T, Yu KH, et al. Organoid models of human and mouse ductal pancreatic cancer. Cell. 2015; 160:324–338.
26. Krausz E, de Hoogt R, Gustin E, Cornelissen F, Grand-Perret T, Janssen L, Vloemans N, Wuyts D, Frans S, Axel A, Peeters PJ, Hall B, Cik M. Translation of a tumor microenvironment mimicking 3D tumor growth co-culture assay platform to high-content screening. Journal of biomolecular screening. 2013; 18:54–66.
27. Hirschhaeuser F, Menne H, Dittfeld C, West J, Mueller-Klieser W, Kunz-Schughart LA. Multicellular tumor spheroids: an underestimated tool is catching up again. Journal of biotechnology. 2010; 148:3–15.
28. Weinstein JN. Drug discovery: Cell lines battle cancer. Nature. 2012; 483:544–545.
29. Birgersdotter A, Sandberg R, Ernberg I. Gene expression perturbation in vitro--a growing case for three-dimensional (3D) culture systems. Semin Cancer Biol. 2005; 15:405–412.
30. Griffith LG, Swartz MA. Capturing complex 3D tissue physiology in vitro. Nat Rev Mol Cell Biol. 2006; 7:211–224.
31. Nelson CM, Bissell MJ. Of extracellular matrix, scaffolds, and signaling: tissue architecture regulates development, homeostasis, and cancer. Annu Rev Cell Dev Biol. 2006; 22:287–309.
32. Jacks T, Weinberg RA. Taking the Study of Cancer Cell Survival to a New Dimension. Cell. 2002; 111:923–925.
33. Abbott A. Cell culture: Biology’s new dimension. Nature. 2003; 424:870–872.
34. Yamada KM, Cukierman E. Modeling tissue morphogenesis and cancer in 3D. Cell. 2007; 130:601–610.
35. Breslin S, O’Driscoll L. Three-dimensional cell culture: the missing link in drug discovery. Drug discovery today. 2013; 18:240-249.
36. Jaganathan H, Gage J, Leonard F, Srinivasan S, Souza GR, Dave B, Godin B. Three-dimensional in vitro co-culture model of breast tumor using magnetic levitation. Sci Rep. 2014; 4:6468.
37. Pampaloni F, Reynaud EG, Stelzer EH. The third dimension bridges the gap between cell culture and live tissue. Nat Rev Mol Cell Biol. 2007; 8:839–845.
38. Kimlin LC, Casagrande G, Virador VM. In vitro three-dimensional (3D) models in cancer research: an update. Molecular carcinogenesis. 2013; 52:167–182.
39. Rudisch A, Dewhurst MR, Horga LG, Kramer N, Harrer N, Dong M, van der Kuip H, Wernitznig A, Bernthaler A, Dolznig H, Sommergruber W. High EMT Signature Score of Invasive Non-Small Cell Lung Cancer (NSCLC) Cells Correlates with NFkappaB Driven Colony-Stimulating Factor 2 (CSF2/GM-CSF) Secretion by Neighboring Stromal Fibroblasts. PloS one. 2015; 10:e0124283.
40. Zargar R, Nourmohammadi J, Amoabediny G. Preparation, characterization, and silanization of 3D microporous PDMS structure with properly sized pores for endothelial cell culture. Biotechnol Appl Biochem. 2016; 63:190–199.
41. Justice BA, Badr NA, Felder RA. 3D cell culture opens new dimensions in cell-based assays. Drug discovery today. 2009; 14:102–107.
42. Vickers AE, Fisher RL. Organ slices for the evaluation of human drug toxicity. Chem Biol Interact. 2004; 150:87–96.
43. Davies EJ, Dong M, Gutekunst M, Narhi K, van Zoggel HJ, Blom S, Nagaraj A, Metsalu T, Oswald E, Erkens-Schulze S, Delgado San Martin JA, Turkki R, Wedge SR, et al. Capturing complex tumour biology in vitro: histological and molecular characterisation of precision cut slices. Sci Rep. 2015; 5:17187.
44. de Graaf IA, Olinga P, de Jager MH, Merema MT, de Kanter R, van de Kerkhof EG, Groothuis GM. Preparation and incubation of precision-cut liver and intestinal slices for application in drug metabolism and toxicity studies. Nature protocols. 2010; 5:1540–1551.
45. Thoma CR, Stroebel S, Rosch N, Calpe B, Krek W, Kelm JM. A high-throughput-compatible 3D microtissue co-culture system for phenotypic RNAi screening applications. Journal of biomolecular screening. 2013; 18:1330–1337.
46. Kelm JM, Djonov V, Ittner LM, Fluri D, Born W, Hoerstrup SP, Fussenegger M. Design of custom-shaped vascularized tissues using microtissue spheroids as minimal building units. Tissue Eng. 2006; 12:2151–2160.
47. Pickl M, Ries CH. Comparison of 3D and 2D tumor models reveals enhanced HER2 activation in 3D associated with an increased response to trastuzumab. Oncogene. 2009; 28:461–468.
48. Zimmermann M, Box C, Eccles SA. Two-dimensional vs. three-dimensional in vitro tumor migration and invasion assays. Methods Mol Biol. 2013; 986:227–252.
49. Nelson CM, Bissell MJ. Modeling dynamic reciprocity: engineering three-dimensional culture models of breast architecture, function, and neoplastic transformation. Semin Cancer Biol. 2005; 15:342–352.
50. Edmondson R, Broglie JJ, Adcock AF, Yang L. Three-Dimensional Cell Culture Systems and Their Applications in Drug Discovery and Cell-Based Biosensors. Assay and Drug Development Technologies. 2014; 12:207–218.
51. Castello-Cros R, Khan DR, Simons J, Valianou M, Cukierman E. Staged stromal extracellular 3D matrices differentially regulate breast cancer cell responses through PI3K and beta1-integrins. BMC Cancer. 2009; 9:94.
52. Ghosh S, Spagnoli GC, Martin I, Ploegert S, Demougin P, Heberer M, Reschner A. Three-dimensional culture of melanoma cells profoundly affects gene expression profile: a high density oligonucleotide array study. J Cell Physiol. 2005; 204:522–531.
53. Rodriguez-Enriquez S, Gallardo-Perez JC, Aviles-Salas A, Marin-Hernandez A, Carreno-Fuentes L, Maldonado-Lagunas V, Moreno-Sanchez R. Energy metabolism transition in multi-cellular human tumor spheroids. J Cell Physiol. 2008; 216:189–197.
54. Shield K, Ackland ML, Ahmed N, Rice GE. Multicellular spheroids in ovarian cancer metastases: Biology and pathology. Gynecol Oncol. 2009; 113:143–148.
55. L’Esperance S, Bachvarova M, Tetu B, Mes-Masson AM, Bachvarov D. Global gene expression analysis of early response to chemotherapy treatment in ovarian cancer spheroids. BMC Genomics. 2008; 9:99.
56. Baharvand H, Hashemi SM, Kazemi Ashtiani S, Farrokhi A. Differentiation of human embryonic stem cells into hepatocytes in 2D and 3D culture systems in vitro. Int J Dev Biol. 2006; 50:645–652.
57. Desoize B, Jardillier J. Multicellular resistance: a paradigm for clinical resistance? Crit Rev Oncol Hematol. 2000; 36:193–207.
58. Friedrich J, Eder W, Castaneda J, Doss M, Huber E, Ebner R, Kunz-Schughart LA. A reliable tool to determine cell viability in complex 3-d culture: the acid phosphatase assay. Journal of biomolecular screening. 2007; 12:925–937.
59. Ravi M, Paramesh V, Kaviya SR, Anuradha E, Solomon FD. 3D cell culture systems: advantages and applications. J Cell Physiol. 2015; 230:16–26.
60. Amann A, Gamerith G, Huber JM, Zwierzina M, Hilbe W, Zwierzina H. Predicting drug sensitivity by 3D cell culture models. memo-Magazine of European Medical Oncology. 2015; 8:77–80.
61. Nath S, Devi GR. Three-dimensional culture systems in cancer research: Focus on tumor spheroid model. Pharmacol Ther. 2016; 163:94–108.
62. Jarde T, Lloyd-Lewis B, Thomas M, Kendrick H, Melchor L, Bougaret L, Watson PD, Ewan K, Smalley MJ, Dale TC. Wnt and Neuregulin1/ErbB signalling extends 3D culture of hormone responsive mammary organoids. Nat Commun. 2016; 7:13207.
63. Humtsoe JO, Pham E, Louie RJ, Chan DA, Kramer RH. ErbB3 upregulation by the HNSCC 3D microenvironment modulates cell survival and growth. Oncogene. 2016; 35:1554–1564.
64. Cisneros Castillo LR, Oancea AD, Stullein C, Regnier-Vigouroux A. Evaluation of Consistency in Spheroid Invasion Assays. Sci Rep. 2016; 6:28375.
65. Huch M, Gehart H, van Boxtel R, Hamer K, Blokzijl F, Verstegen MM, Ellis E, van Wenum M, Fuchs SA, de Ligt J, van de Wetering M, Sasaki N, Boers SJ, et al. Long-term culture of genome-stable bipotent stem cells from adult human liver. Cell. 2015; 160:299–312.
66. Weigelt B, Ghajar CM, Bissell MJ. The need for complex 3D culture models to unravel novel pathways and identify accurate biomarkers in breast cancer. Advanced drug delivery reviews. 2014; 69-70:42–51.
67. Stratmann AT, Fecher D, Wangorsch G, Gottlich C, Walles T, Walles H, Dandekar T, Dandekar G, Nietzer SL. Establishment of a human 3D lung cancer model based on a biological tissue matrix combined with a Boolean in silico model. Molecular oncology. 2014; 8:351–365.
68. Shamir ER, Ewald AJ. Three-dimensional organotypic culture: experimental models of mammalian biology and disease. Nat Rev Mol Cell Biol. 2014; 15:647–664.
69. Karthaus WR, Iaquinta PJ, Drost J, Gracanin A, van Boxtel R, Wongvipat J, Dowling CM, Gao D, Begthel H, Sachs N, Vries RG, Cuppen E, Chen Y, et al. Identification of multipotent luminal progenitor cells in human prostate organoid cultures. Cell. 2014; 159:163–175.
70. Gao D, Vela I, Sboner A, Iaquinta PJ, Karthaus WR, Gopalan A, Dowling C, Wanjala JN, Undvall EA, Arora VK, Wongvipat J, Kossai M, Ramazanoglu S, et al. Organoid cultures derived from patients with advanced prostate cancer. Cell. 2014; 159:176–187.
71. Lynch TJ, Bell DW, Sordella R, Gurubhagavatula S, Okimoto RA, Brannigan BW, Harris PL, Haserlat SM, Supko JG, Haluska FG, Louis DN, Christiani DC, Settleman J, et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004; 350:2129–2139.
72. Roskoski R Jr. The ErbB/HER family of protein-tyrosine kinases and cancer. Pharmacol Res. 2014; 79:34–74.
73. Kobayashi S, Boggon TJ, Dayaram T, Janne PA, Kocher O, Meyerson M, Johnson BE, Eck MJ, Tenen DG, Halmos B. EGFR mutation and resistance of non-small-cell lung cancer to gefitinib. N Engl J Med. 2005; 352:786–792.
74. Mukohara T, Engelman JA, Hanna NH, Yeap BY, Kobayashi S, Lindeman N, Halmos B, Pearlberg J, Tsuchihashi Z, Cantley LC, Tenen DG, Johnson BE, Janne PA. Differential effects of gefitinib and cetuximab on non-small-cell lung cancers bearing epidermal growth factor receptor mutations. Journal of the National Cancer Institute. 2005; 97:1185–1194.
75. Song JM, Anandharaj A, Upadhyaya P, Kirtane AR, Kim JH, Hong KH, Panyam J, Kassie F. Honokiol suppresses lung tumorigenesis by targeting EGFR and its downstream effectors. Oncotarget. 2016; 7:57752–57769. http://doi.org/10.18632/oncotarget.10759.
76. Bell DW, Gore I, Okimoto RA, Godin-Heymann N, Sordella R, Mulloy R, Sharma SV, Brannigan BW, Mohapatra G, Settleman J, Haber DA. Inherited susceptibility to lung cancer may be associated with the T790M drug resistance mutation in EGFR. Nat Genet. 2005; 37:1315–1316.
77. Shih JY, Gow CH, Yang PC. EGFR mutation conferring primary resistance to gefitinib in non-small-cell lung cancer. N Engl J Med. 2005; 353:207–208.
78. Riss TL, Moravec TA, Niles AL, Duellman S, Benink HA, Worzella TJ, Minor L. Cell Viability Assays BTI - Assay Guidance Manual. https://www.ncbi.nlm.nih.gov/books/NBK144065/.
79. Gow CH, Shih JY, Chang YL, Yu CJ. Acquired gefitinib-resistant mutation of EGFR in a chemonaive lung adenocarcinoma harboring gefitinib-sensitive mutation L858R. PLoS Med. 2005; 2:e269.
80. Phuchareon J, McCormick F, Eisele DW, Tetsu O. EGFR inhibition evokes innate drug resistance in lung cancer cells by preventing Akt activity and thus inactivating Ets-1 function. Proceedings of the National Academy of Sciences of the United States of America. 2015; 112:E3855–3863.
81. Ying H, Yang XD, Sun Z, Ning X, Wang Y, Bai C, Chen S, Wang Y. Lifestyle risks exposure and response predictor of gefitinib in patients with non-small cell lung cancer. Med Oncol. 2014; 31:220.
82. Rose MC, Kostyanovskaya E, Huang RS. Pharmacogenomics of cisplatin sensitivity in non-small cell lung cancer. Genomics, proteomics & bioinformatics. 2014; 12:198–209.
83. Fennell DA, Summers Y, Cadranel J, Benepal T, Christoph DC, Lal R, Das M, Maxwell F, Visseren-Grul C, Ferry D. Cisplatin in the modern era: The backbone of first-line chemotherapy for non-small cell lung cancer. Cancer Treat Rev. 2016; 44:42–50.
84. Joerger M, von Pawel J, Kraff S, Fischer JR, Eberhardt W, Gauler TC, Mueller L, Reinmuth N, Reck M, Kimmich M, Mayer F, Kopp HG, Behringer DM, et al. Open-label, randomized study of individualized, pharmacokinetically (PK)-guided dosing of paclitaxel combined with carboplatin or cisplatin in patients with advanced non-small-cell lung cancer (NSCLC). Ann Oncol. 2016; 27:1895–1902.
85. Friberg S, Mattson S. On the growth rates of human malignant tumors: implications for medical decision making. J Surg Oncol. 1997; 65:284–297.
86. Spratt JS, Meyer JS, Spratt JA. Rates of growth of human neoplasms: Part II. J Surg Oncol. 1996; 61:68–83.
87. Spratt JS Jr, Spratt TL. Rates of Growth of Pulmonary Metastases and Host Survival. Ann Surg. 1964; 159:161–171.
88. Arai T, Kuroishi T, Saito Y, Kurita Y, Naruke T, Kaneko M. Tumor doubling time and prognosis in lung cancer patients: evaluation from chest films and clinical follow-up study. Japanese Lung Cancer Screening Research Group. Jpn J Clin Oncol. 1994; 24:199–204.
89. Freyer JP, Schor PL. Regrowth kinetics of cells from different regions of multicellular spheroids of four cell lines. J Cell Physiol. 1989; 138:384–392.
90. Jiang X, Li H, Xie J, McKinley ET, Zhao P, Gore JC, Xu J. In vivo imaging of cancer cell size and cellularity using temporal diffusion spectroscopy. Magn Reson Med. 2017; 78:156–164.
91. McMillin DW, Negri JM, Mitsiades CS. The role of tumour-stromal interactions in modifying drug response: challenges and opportunities. Nature reviews Drug discovery. 2013; 12:217–228.
92. Garnett MJ, Edelman EJ, Heidorn SJ, Greenman CD, Dastur A, Lau KW, Greninger P, Thompson IR, Luo X, Soares J, Liu Q, Iorio F, Surdez D, et al. Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature. 2012; 483:570–575.
93. Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehar J, Kryukov GV, Sonkin D, Reddy A, Liu M, Murray L, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012; 483:603–607.
94. Shoemaker RH. The NCI60 human tumour cell line anticancer drug screen. Nature reviews Cancer. 2006; 6:813–823.
95. Jung SK, Lee MH, Lim DY, Lee SY, Jeong CH, Kim JE, Lim TG, Chen H, Bode AM, Lee HJ, Lee KW, Dong Z. Butein, a novel dual inhibitor of MET and EGFR, overcomes gefitinib-resistant lung cancer growth. Molecular carcinogenesis. 2014; 54:322–331.
96. Chan GK, Kleinheinz TL, Peterson D, Moffat JG. A simple high-content cell cycle assay reveals frequent discrepancies between cell number and ATP and MTS proliferation assays. PloS one. 2013; 8:e63583.
97. Jang WJ, Jung SK, Kang JS, Jeong JW, Bae MK, Joo SH, Park GH, Kundu JK, Hong YS, Jeong CH. Anti-tumor activity of WK88-1, a novel geldanamycin derivative, in gefitinib-resistant non-small cell lung cancers with Met amplification. Cancer science. 2014; 105:1245–1253.
98. Zou M, Xia S, Zhuang L, Han N, Chu Q, Chao T, Peng P, Chen Y, Gui Q, Yu S. Knockdown of the Bcl-2 gene increases sensitivity to EGFR tyrosine kinase inhibitors in the H1975 lung cancer cell line harboring T790M mutation. International journal of oncology. 2013; 42:2094–2102.
99. Song JY, Kim CS, Lee JH, Jang SJ, Lee SW, Hwang JJ, Lim C, Lee G, Seo J, Cho SY, Choi J. Dual inhibition of MEK1/2 and EGFR synergistically induces caspase-3-dependent apoptosis in EGFR inhibitor-resistant lung cancer cells via BIM upregulation. Investigational new drugs. 2013; 31:1458–1465.
100. Kwak EL, Sordella R, Bell DW, Godin-Heymann N, Okimoto RA, Brannigan BW, Harris PL, Driscoll DR, Fidias P, Lynch TJ, Rabindran SK, McGinnis JP, Wissner A, et al. Irreversible inhibitors of the EGF receptor may circumvent acquired resistance to gefitinib. Proceedings of the National Academy of Sciences of the United States of America. 2005; 102:7665–7670.
101. Merz F, Gaunitz F, Dehghani F, Renner C, Meixensberger J, Gutenberg A, Giese A, Schopow K, Hellwig C, Schafer M, Bauer M, Stocker H, Taucher-Scholz G, et al. Organotypic slice cultures of human glioblastoma reveal different susceptibilities to treatments. Neuro Oncol. 2013; 15:670–681.
102. Dolznig H, Rupp C, Puri C, Haslinger C, Schweifer N, Wieser E, Kerjaschki D, Garin-Chesa P. Modeling colon adenocarcinomas in vitro a 3D co-culture system induces cancer-relevant pathways upon tumor cell and stromal fibroblast interaction. The American journal of pathology. 2011; 179:487–501.
103. Frieboes HB, Zheng X, Sun CH, Tromberg B, Gatenby R, Cristini V. An integrated computational/experimental model of tumor invasion. Cancer research. 2006; 66:1597–1604.
104. Erlichman C, Tannock IF. Growth, characterization of multicellular tumor spheroids of human bladder carcinoma origin. In Vitro Cell Dev Biol. 1986; 22:449–456.
105. De Sousa EMF, Vermeulen L, Fessler E, Medema JP. Cancer heterogeneity--a multifaceted view. EMBO Rep. 2013; 14:686–695.
106. Mueller-Klieser W. Tumor biology and experimental therapeutics. Critical Reviews in Oncology/Hematology. 2000; 36:123–139.
107. Chignola R. Computational Challenges of tumor spheroid modeling. Journal of Bioinformatics and Computational Biology. 2010; 9:559–577.
108. Maemondo M, Inoue A, Kobayashi K, Sugawara S, Oizumi S, Isobe H, Gemma A, Harada M, Yoshizawa H, Kinoshita I, Fujita Y, Okinaga S, Hirano H, et al. Gefitinib or chemotherapy for non-small-cell lung cancer with mutated EGFR. N Engl J Med. 2010; 362:2380–2388.
109. Douillard JY, Ostoros G, Cobo M, Ciuleanu T, McCormack R, Webster A, Milenkova T. First-line gefitinib in Caucasian EGFR mutation-positive NSCLC patients: a phase-IV, open-label, single-arm study. Br J Cancer. 2014; 110:55–62.
110. Weigelt B, Lo AT, Park CC, Gray JW, Bissell MJ. HER2 signaling pathway activation and response of breast cancer cells to HER2-targeting agents is dependent strongly on the 3D microenvironment. Breast Cancer Res Treat. 2010; 122:35–43.
111. Eimer S, Dugay F, Airiau K, Avril T, Quillien V, Belaud-Rotureau MA, Belloc F. Cyclopamine cooperates with EGFR inhibition to deplete stem-like cancer cells in glioblastoma-derived spheroid cultures. Neuro Oncol. 2012; 14:1441–1451.
112. Ekert JE, Johnson K, Strake B, Pardinas J, Jarantow S, Perkinson R, Colter DC. Three-dimensional lung tumor microenvironment modulates therapeutic compound responsiveness in vitro—implication for drug development. PloS one. 2014; 9:e92248.
113. Lama R, Zhang L, Naim JM, Williams J, Zhou A, Su B. Development, validation and pilot screening of an in vitro multi-cellular three-dimensional cancer spheroid assay for anti-cancer drug testing. Bioorg Med Chem. 2013; 21:922–931.
114. Kim SH, Choi SJ, Kim YC, Kuh HJ. Anti-tumor activity of noble indirubin derivatives in human solid tumor models in vitro. Arch Pharm Res. 2009; 32:915–922.
115. Chambers KF, Mosaad EMO, Russell PJ, Clements JA, Doran MR. 3D Cultures of Prostate Cancer Cells Cultured in a Novel High-Throughput Culture Platform Are More Resistant to Chemotherapeutics Compared to Cells Cultured in Monolayer. PloS one. 2014; 9:e111029.
116. Friedrich J, Seidel C, Ebner R, Kunz-Schughart LA. Spheroid-based drug screen: considerations and practical approach. Nature protocols. 2009; 4:309–324.
117. Singhirunnusorn P, Ueno Y, Matsuo M, Suzuki S, Saiki I, Sakurai H. Transient suppression of ligand-mediated activation of epidermal growth factor receptor by tumor necrosis factor-alpha through the TAK1-p38 signaling pathway. J Biol Chem. 2007; 282:12698–12706.
118. Bissel MJ. The morphologies of breast cancer cell lines in three-dimensional assays correlate with their profiles of gene expression. Molecular oncology. 2007; 1:84–96.
119. Zschenker O, Streichert T, Hehlgans S, Cordes N. Genome-wide gene expression analysis in cancer cells reveals 3D growth to affect ECM and processes associated with cell adhesion but not DNA repair. PloS one. 2012; 7:e34279.
120. Luca AC, Mersch S, Deenen R, Schmidt S, Messner I, Schafer KL, Baldus SE, Huckenbeck W, Piekorz RP, Knoefel WT, Krieg A, Stoecklein NH. Impact of the 3D microenvironment on phenotype, gene expression, and EGFR inhibition of colorectal cancer cell lines. PloS one. 2013; 8:e59689.
121. Pugazhenthi S, Nesterova A, Sable C, Heidenreich KA, Boxer LM, Heasley LE, Reusch JE. Akt/protein kinase B up-regulates Bcl-2 expression through cAMP-response element-binding protein. J Biol Chem. 2000; 275:10761–10766.
122. Ricciardi MR, Mirabilii S, Licchetta R, Piedimonte M, Tafuri A. Targeting the Akt, GSK-3, Bcl-2 axis in acute myeloid leukemia. Adv Biol Regul. 2017;65:36–58.
123. Zhou C, Ding J, Wu Y. Resveratrol induces apoptosis of bladder cancer cells via miR21 regulation of the Akt/Bcl2 signaling pathway. Mol Med Rep. 2014; 9:1467–1473.
124. Kunz-Schughart LA, Heyder P, Schroeder J, Knuechel R. A heterologous 3-D coculture model of breast tumor cells and fibroblasts to study tumor-associated fibroblast differentiation. Experimental cell research. 2001; 266:74–86.
125. Bhadriraju K, Chen S. Engineering cellular microenvironments to improve cell-based drug testing. Drug discovery today. 2002; 7:612–620.
126. Rupp C, Scherzer M, Rudisch A, Unger C, Haslinger C, Schweifer N, Artaker M, Nivarthi H, Moriggl R, Hengstschlager M, Kerjaschki D, Sommergruber W, Dolznig H, et al. IGFBP7, a novel tumor stroma marker, with growth-promoting effects in colon cancer through a paracrine tumor-stroma interaction. Oncogene. 2015; 34:815–825.
127. Wang W, Li Q, Yamada T, Matsumoto K, Matsumoto I, Oda M, Watanabe G, Kayano Y, Nishioka Y, Sone S, Yano S. Crosstalk to stromal fibroblasts induces resistance of lung cancer to epidermal growth factor receptor tyrosine kinase inhibitors. Clin Cancer Res. 2009; 15:6630–6638.
128. Choe C, Shin YS, Kim C, Choi SJ, Lee J, Kim SY, Cho YB, Kim J. Crosstalk with cancer-associated fibroblasts induces resistance of non-small cell lung cancer cells to epidermal growth factor receptor tyrosine kinase inhibition. Onco Targets Ther. 2015; 8:3665–3678.
129. Mink SR, Vashistha S, Zhang W, Hodge A, Agus DB, Jain A. Cancer-associated fibroblasts derived from EGFR-TKI-resistant tumors reverse EGFR pathway inhibition by EGFR-TKIs. Molecular cancer research : MCR. 2010; 8:809–820.
130. Li L, Zhou Q, Voss TC, Quick KL, LaBarbera DV. High-throughput imaging: Focusing in on drug discovery in 3D. Methods. 2016; 96:97–102.
131. Howes AL, Richardson RD, Finlay D, Vuori K. 3-Dimensional culture systems for anti-cancer compound profiling and high-throughput screening reveal increases in EGFR inhibitor-mediated cytotoxicity compared to monolayer culture systems. PloS one. 2014; 9:e108283.
132. Laemmli UK. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970; 227:680–685.
133. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP. Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003; 31:e15.
134. 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, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004; 5:R80.
135. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004; 3:Article3.
136. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001; 125:279–284.
137. Smyth G. Limma: linear models for microarray data. New York: Springer. 2005.
138. Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010; 26:2363–2367.
139. Kauffmann A, Gentleman R, Huber W. arrayQualityMetrics—a bioconductor package for quality assessment of microarray data. Bioinformatics. 2009; 25:415–416.
140. Alter O, Brown PO, Botstein D. Singular value decomposition for genome-wide expression data processing and modeling. Proceedings of the National Academy of Sciences of the United States of America. 2000; 97:10101–10106.
141. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008; 9:559.
