Oscillating growth patterns of multicellular tumour spheroids

The growth kinetics of 9L (rat glioblastoma cell line) and U118 (human glioblastoma cell line) multicellular tumour spheroids (MTS) have been investigated by non-linear least square fitting of individual growth curves with the Gompertz growth equation and power spectrum analysis of residuals. Residuals were not randomly distributed around calculated growth trajectories. At least one main frequency was found for all analysed MTS growth curves, demonstrating the existence of time-dependent periodic fluctuations of MTS volume dimensions. Similar periodic oscillations of MTS volume dimensions were also observed for MTS generated using cloned 9L cells. However, we found significant differences in the growth kinetics of MTS obtained with cloned cells if compared to the growth kinetics of MTS obtained with polyclonal cells. Our findings demonstrate that the growth patterns of three-dimensional tumour cell cultures are more complex than has been previously predicted using traditional continuous growth models.     

Oscillations in growth of multicellular tumour spheroids: a revisited quantitative analysis

Objectives: Multicellular tumour spheroids (MTS) provide an important tool for study of the microscopic properties of solid tumours and their responses to therapy. Thus, observation of large‐scale volume oscillations in MTS, reported several years ago by two independent groups ( 1 , 2 ), in our opinion represent a remarkable discovery, particularly if this could promote careful investigation of the possible occurrence of volume oscillations of tumours ‘in vivo’. Materials and methods: Because of high background noise, quantitative analysis of properties of observed oscillations has not been possible in previous studies. Such an analysis can be now performed, thanks to a recently proposed approach, based on formalism of phenomenological universalities (PUN). Results: Results have provided unambiguous confirmation of the existence of MTS volume oscillations, and quantitative evaluation of their properties, for two tumour cell lines. Proof is based not only on quality of fitting of the experimental datasets, but also on determination of well‐defined values of frequency and amplitude of the oscillations for each line investigated, which would not be consistent with random fluctuation. Conclusions: Biological mechanisms, which can be directly responsible for observed oscillations, are proposed, which relates also to recent work on related topics. Further investigations, both at experimental and at modelling levels, are also suggested. Finally, from a methodological point of view, results obtained represent further confirmation of applicability and usefulness of the PUN approach.     

Mathematical modelling of microenvironment and growth in EMT6/Ro multicellular tumour spheroids

In order to determine the role of micromilieu in tumour spheroid growth, a mathematical model was developed to predict EMT6/Ro spheroid growth and microenvironment based upon numerical solution of the diffusion/reaction equation for oxygen, glucose, lactate ion, carbon dioxide, bicarbonate ion, chlorine ion and hydrogen ion along with the equation of electroneutrality. This model takes into account the effects of oxygen concentration, glucose concentration and extracellular pH on cell growth and metabolism. Since independent measurements of EMT6/Ro single cell growth and metabolic rates, spheroid diffusion constants, and spinner flask mass transfer coefficients are available, model predictions using these parameters were compared with published data on EMT6/Ro spheroid growth and micro-environment. The model predictions of reduced spheroid growth due to reduced cell growth rates and cell shedding fit experimental spheroid growth data below 700 microns, but overestimated the spheroid growth rate at larger diameters. Predicted viable rim thicknesses based on predicted near zero glucose concentrations fit published viable rim thickness data for 1000 microns spheroids grown at medium glucose concentrations of 5.5 mM or less. However, the model did not accurately predict the onset of necrosis. Moreover, the model could not predict the observed decreases in oxygen and glucose metabolism seen in spheroids with time, nor could it predict the observed growth plateau. This suggests that other unknown factors, such as inhibitors or cell-cell contact effects, must also be important in affecting spheroid growth and cellular metabolism.     

Down-regulation of DNA mismatch repair enhances initiation and growth of neuroblastoma and brain tumour multicellular spheroids

Multicellular tumour spheroid (MCTS) cultures are excellent model systems for simulating the development and microenvironmental conditions of in vivo tumour growth. Many documented cell lines can generate differentiated MCTS when cultured in suspension or in a non-adhesive environment. While physiological and biochemical properties of MCTS have been extensively characterized, insight into the events and conditions responsible for initiation of these structures is lacking. MCTS are formed by only a small subpopulation of cells during surface-associated growth but the processes responsible for this differentiation are poorly understood and have not been previously studied experimentally. Analysis of gene expression within spheroids has provided clues but to date it is not known if the observed differences are a cause or consequence of MCTS growth. One mechanism linked to tumourigenesis in a number of cancers is genetic instability arising from impaired DNA mismatch repair (MMR). This study aimed to determine the role of MMR in MCTS initiation and development. Using surface-associated N2a and CHLA-02-ATRT culture systems we have investigated the impact of impaired MMR on MCTS growth. Analysis of the DNA MMR genes MLH1 and PMS2 revealed both to be significantly down-regulated at the mRNA level compared with non-spheroid-forming cells. By using small interfering RNA (siRNA) against these genes we show that silencing of MLH1 and PMS2 enhances both MCTS initiation and subsequent expansion. This effect was prolonged over several passages following siRNA transfection. Down-regulation of DNA MMR can contribute to tumour initiation and progression in N2a and CHLA-02-ATRT MCTS models. Studies of surface-associated MCTS differentiation may have broader applications in studying events in the initiation of cancer foci.     

Semiautomatic growth analysis of multicellular tumor spheroids

Multicellular tumor spheroids (MCTS) are routinely employed as three-dimensional in vitro models to study tumor biology. Cultivation of MCTS in spinner flasks provides better growing conditions, especially with regard to the availability of nutrients and oxygen, when compared with microtiter plates. The main endpoint of drug response experiments is spheroid size. It is common practice to analyze spheroid size manually with a microscope and an ocular micrometer. This requires removal of some spheroids from the flask, which entails major limitations such as loss of MCTS and the risk of contamination. With this new approach, the authors present an efficient and highly reproducible method to analyze the size of complete MCTS populations in culture containers with transparent, flat bottoms. MCTS sediments are digitally scanned and spheroid volumes are calculated by computerized image analysis. The equipment includes regular office hardware (personal computer, flatbed scanner) and software (Adobe Photoshop, Microsoft Excel, ImageJ). The accuracy and precision of the method were tested using industrial precision steel beads with known diameter. In summary, in comparison with other methods, this approach provides benefits in terms of semiautomation, noninvasiveness, and low costs.     

Growth characteristics of human melanoma multicellular spheroids in liquid-overlay culture: comparisons with the parent tumour xenografts

The growth characteristics of multicellular spheroids, derived from human melanoma xenografts and cultivated in liquid-overlay culture, were studied and compared with those of the parent tumours. Six of the seven melanomas investigated formed spheroids, which grew exponentially up to a volume of 1-2 X 10(7) microns 3 (a diameter of 270-340 microns) before the growth rate tapered off. The morphology of the spheroids varied considerably among the melanomas; some spheroids grew as densely packed, spherical structures of cells whereas others were loosely packed and showed an irregular shape. Central necrosis developed when the spheroids attained a diameter of 150-200 microns. The histological and cytological appearance of the spheroids was remarkably similar to that of the parent xenograft in five of the six cases. The sixth melanoma contained two subpopulations with distinctly different DNA content, one of which was predominant in the spheroids, the other in the tumours. This gave rise to clear histological and cytological differences. The volume-doubling time of the spheroids during the exponential growth phase ranged from 1.7 +/- 0.2 to 2.7 +/- 0.4 days and the fraction of cells in S from 13 +/- 1 to 28 +/- 2%. The volume-doubling time decreased with increasing fraction of cells in S, indicating that the differences in growth rate were due mainly to differences in the growth fraction or to differences in the duration of G1. The spheroid volume-doubling times did not correlate with those of the parent xenografts (Td = 4.2-22.5 days at V = 200 mm3), possibly because the cell loss factors of the xenografts were large and varied among the melanomas. The fractions of cells in G1/G0, S and G2 + M in the spheroids and the xenografts did not correlate either, but were found to be within the same narrow ranges in the spheroids and the xenografts--i.e. 50-80% (G1/G0), 10-30% (S) and 10-20% (G2 + M).     

Changes in P-glycoprotein activity are mediated by the growth of a tumour cell line as multicellular spheroids

Background  

Expression of P-glycoprotein (P-gp), the multidrug resistance (MDR) 1 gene product, can lead to multidrug resistance in tumours. However, the physiological role of P-gp in tumours growing as multicellular spheroids is not well understood. Recent evidence suggests that P-gp activity may be modulated by cellular components such as membrane proteins, membrane-anchoring proteins or membrane-lipid composition. Since, multicellular spheroids studies have evidenced alterations in numerous cellular components, including those related to the plasma membrane function, result plausible that some of these changes might modulate P-gp function and be responsible for the acquisition of multicellular drug resistance. In the present study, we asked if a human lung cancer cell line (INER-51) grown as multicellular spheroids can modify the P-gp activity to decrease the levels of doxorubicin (DXR) retained and increase their drug resistance.     

A model for the growth of multicellular spheroids

Abstract. Based on biological observations and the basic physical properties of tri-dimensional structures, a mathematical expression is derived to relate the growth rate of multicellular spheroids to some easily measurable parameters. This model involves properties both of the individual cells and of the spheroid structure, such as the cell doubling time in monolayer, the rate of cell shedding from the spheroid and the depth of the external rim of cycling cells. The derived growth equation predicts a linear expansion of the spheroid diameter with time. The calculated growth rate for a number of spheroid cell types is in good agreement with experimental data. The model provides a simple and practical view of growth control in spheroids, and is further adapted to include parameters presumably responsible for the growth saturation in large spheroids.     

Induction of apoptosis in the centre of multicellular tumour spheroids of colorectal adenocarcinomas--involvement of CD95 pathway and differentiation

Multicellular tumour spheroids (MCTS) are three-dimensional cell culture systems which are widely used in cancer research. They are characterized by an outer zone of proliferating cells, an inner region of differentiating quiescent cells and an area of so-called necrotic cell death in their centre. The exact cause of this cell death, a controversy for many years, was the aim of the present study. Our data show that cell death in the centre of MCTS of three colorectal adenocarcinoma cell lines (HRT-18, HT-29 and CX-2) was induced by apoptosis. Apoptotic cells were initially distributed at random but accumulated very quickly in the quiescent and central area at day 4-5, suggesting a time- rather than size-dependent synchronization of apoptosis parallel to the formation of the proliferation gradient in MCTS. To study mechanisms inducing apoptosis, the Fas-pathway was investigated. A cell-cell contact-dependent expression of CD95 was found in all MCTS. FasL was not detected in monolayer cultures, but was expressed in spheroids of HRT-18 and CX-2. We found that TNF alpha and TGF beta 1 activated the CD95 pathway in all three cell lines. Since both TNF-alpha and TGF-beta are known to be inducible by hypoxia in a variety of cell types, we suggest that these hypoxia-induced factors sensitize the CD95 pathway in the quiescent area of MCTS. Furthermore, a loss of the heat shock proteins 27, 32, 60, 73 and 90 was observed in the quiescent area of spheroids. This suggests that tumour cell differentiation in the inner region of MCTS may be an additional factor inducing apoptosis.     

Viable sorting of intact multicellular spheroids by flow cytometry

A flow cytometric method has been developed for sorting viable, intact multicellular spheroids in order to obtain uniformly-sized populations with diameters in the range of 50-100 microns. A FACS II instrument was modified for this purpose by installing a 200-microns-diameter exit orifice and by making adjustments in the sheath flow, oscillator frequency, and number of droplets sorted. Polystyrene microspheres (44 and 88 microns diameter) and 41-96-microns-diameter spheroids could be sorted and recovered with 70-100% efficiency, an improvement over previous reports. Unstained, viable spheroids were simultaneously analyzed for small-angle forward light scatter, 90 degree light scatter, and autofluorescence using a 488-nm laser operating at 100 mW. Analysis of the data demonstrated a considerable variation in both the 90 degrees light scatter and the autofluorescence signals for a given forward angle light scattering signal. By setting narrow sort windows on the forward angle light scattering signal and either the 90 degree light scatter or autofluorescence signals, uniformly spherical spheroid populations could be recovered. These sorted populations had coefficients of variation of the mean diameter in the range of 5-9%. This represents a variation of less than one cell diameter, and is a major improvement over any other technique. There was no significant difference in the subsequent growth rates of sorted spheroids compared to the unsorted spheroids. This technique will apply when uniform populations of small spheroids are required, such as investigations of the contact effect or in the initiation of growth curve studies.     

The Bcl-2 repertoire of mesothelioma spheroids underlies acquired apoptotic multicellular resistance

Three-dimensional (3D) cultures are a valuable platform to study acquired multicellular apoptotic resistance of cancer. We used spheroids of cell lines and actual tumor to study resistance to the proteasome inhibitor bortezomib in mesothelioma, a highly chemoresistant tumor. Spheroids from mesothelioma cell lines acquired resistance to bortezomib by failing to upregulate Noxa, a pro-apoptotic sensitizer BH3-only protein that acts by displacing Bim, a pro-apoptotic Bax/Bak-activator protein. Surprisingly, despite their resistance, spheroids also upregulated Bim and thereby acquired sensitivity to ABT-737, an inhibitor of anti-apoptotic Bcl-2 molecules. Analysis using BH3 profiling confirmed that spheroids acquired a dependence on anti-apoptotic Bcl-2 proteins and were ‘primed for death''. We then studied spheroids grown from actual mesothelioma. ABT-737 was active in spheroids grown from those tumors (5/7, ∼70%) with elevated levels of Bim. Using immunocytochemistry of tissue microarrays of 48 mesotheliomas, we found that most (33, 69%) expressed elevated Bim. In conclusion, mesothelioma cells in 3D alter the expression of Bcl-2 molecules, thereby acquiring both apoptotic resistance and sensitivity to Bcl-2 blockade. Mesothelioma tumors ex vivo also show sensitivity to Bcl-2 blockade that may depend on Bim, which is frequently elevated in mesothelioma. Therefore, mesothelioma, a highly resistant tumor, may have an intrinsic sensitivity to Bcl-2 blockade that can be exploited therapeutically.     

Increased thermoresistance developed during growth of small multicellular spheroids

Mammalian cells growing as multicell spheroids, an in vitro model of tumor microregions, have been shown previously to be more resistant than single cells from monolayer cultures to killing by ionizing radiation, hyperthermia, ultrasound, and chemotherapeutic drugs. Although the mechanisms by which cells in spheroids acquire these increased resistances are unknown, available evidence has indicated that intercellular contact mediates the process for ionizing radiation. This investigation was undertaken to evaluate the role of intercellular contact produced during growth of small spheroids on the sensitivity of EMT6/Ro mouse mammary tumor cells to moderate hyperthermia. Increased thermoresistance developed in small spheroids (approximately 70 micron diameter, 25 cells/spheroid), as measured by colony formation, after exposures to different temperatures in the range of 37 to 45 degrees C for periods less than or equal to 2 hr and at 42.5 degrees C for less than or equal to 8 hr. Experiments were performed to determine the relative contributions to this increased thermoresistance of 1) the extent of intercellular contact in spheroids of different cellular multiplicities, 2) differences in membrane damage influenced by trypsin heat treatment sequence, and 3) physiological changes associated with growth of cells as spheroids in suspension compared to monolayer culture. Treatment with trypsin prior to heating sensitized cells to killing by hyperthermia but did not account for the differential thermoresistance between cells from spheroids and monolayers. Spheroid multiplicity in the range of 1.16 to 76.2 cells/spheroid had no significant effect on cell survival after hyperthermia. However, cells grown in spinner suspension culture were more thermoresistant than cells from monolayer cultures and nearly as thermoresistant as cells in spheroids. From these data we conclude that the greater thermoresistance of EMT/Ro cells in spheroids is the result of cellular physiological changes associated with growth in suspension and is not mediated by intercellular contact.     

Spatio-temporal modeling of nanoparticle delivery to multicellular tumor spheroids

The inefficiency of nanoparticle penetration in tissues limits the therapeutic efficacy of such formulations for cancer applications. Recent work has indicated that modulation of tissue architecture with enzymes such as collagenase significantly increases macromolecule delivery. In this study we developed a mathematical model of nanoparticle penetration into multicellular spheroids that accounts for radially dependent changes in tumor architecture, as represented by the volume fraction of tissue accessible to nanoparticle diffusion. Parameters such as nanoparticle binding, internalization rate constants, and accessible volume fraction were determined experimentally. Unknown parameters of nanoparticle binding sites per cell in the spheroid and pore shape factor were determined by fitting to experimental data. The model was correlated with experimental studies of the penetration of 40 nm nanoparticles in SiHa multicellular spheroids with and without collagenase treatment and was able to accurately predict concentration profiles of nanoparticles within spheroids. The model was also used to investigate the effects of nanoparticle size. This model contributes toward the understanding of the role of tumor architecture on nanoparticle delivery efficiency.     

Metabolic imaging in multicellular spheroids of oncogene-transfected fibroblasts.

Four rat embryo fibroblast (REF) cell lines with defined oncogenic transformation were used to study the relationship between tumorigenic conversion, metabolism, and development of cell death in a 3D spheroid system. Rat1 (spontaneously immortalized) and M1 (myc-transfected) fibroblasts represent early nontumorigenic transformation stages, whereas Rat1-T1 (T24Ha-ras-transfected Rat1) and MR1 (myc/T24Ha-ras-co-transfected REF) cells express a highly tumorigenic phenotype. Localized ATP, glucose, and lactate concentrations in spheroid median sections were determined by imaging bioluminescence. ATP concentrations were low in the nonproliferating Rat1 aggregates despite sufficient oxygen and glucose availability and lack of lactate accumulation. In MR1 spheroids, a 50% decrease in central ATP preceded the development of central necrosis at a spheroid diameter of around 800 micrometer. In contrast, the histomorphological emergence of cell death at a diameter of around 500 micrometer in Rat1-T1 spheroids coincided with an initial steep drop in ATP. Concomitantly, reduction in central glucose and increase in lactate before cell death were recorded in MR1 but not in Rat1-T1 spheroids. As shown earlier, myc transfection confers a considerable resistance to hypoxia of MR1 cells in the center of spheroids, which is reflected by their capability to maintain cell integrity and ATP content in a hypoxic environment. The data obtained suggest that small alterations in the genotype of tumor cell lines, such as differences in the immortalization process, lead to substantial differences in morphological structure, metabolism, occurrence of cell death, and tolerance to hypoxia in spheroid culture.     

Visualizing the effect of tumor microenvironments on radiation-induced cell kinetics in multicellular spheroids consisting of HeLa cells

In this study, we visualized the effect of tumor microenvironments on radiation-induced tumor cell kinetics. For this purpose, we utilized a multicellular spheroid model, with a diameter of ∼500 μm, consisting of HeLa cells expressing the fluorescent ubiquitination-based cell-cycle indicator (Fucci). In live spheroids, a confocal laser scanning microscope allowed us to clearly monitor cell kinetics at depths of up to 60 μm. Surprisingly, a remarkable prolongation of G2 arrest was observed in the outer region of the spheroid relative to monolayer-cultured cells. Scale, an aqueous reagent that renders tissues optically transparent, allowed visualization deeper inside spheroids. About 16 h after irradiation, a red fluorescent cell fraction, presumably a quiescent G0 cell fraction, became distinct from the outer fraction consisting of proliferating cells, most of which exhibited green fluorescence indicative of G2 arrest. Thereafter, the red cell fraction began to emit green fluorescence and remained in prolonged G2 arrest. Thus, for the first time, we visualized the prolongation of radiation-induced G2 arrest in spheroids and the differences in cell kinetics between the outer and inner fractions.     

Redox-regulation of intrinsic prion expression in multicellular prostate tumor spheroids

The cellular function of the intrinsic prion protein (PrPc) remains largely unknown. In the present study PrPc expression was investigated in multicellular prostate tumor spheroids and was correlated to the intracellular redox state as evaluated using the fluorescent dye 2'7'-dichlorodihydrofluorescein diacetate (H2DCFDA). In small tumor spheroids (diameter 100 +/- 20 microm) reactive oxygen species (ROS) levels were increased as compared with large (diameter 250 +/- 50 microm) spheroids. ROS generation was mediated by the mitochondrial respiratory chain and a NADPH oxidaselike enzyme, because carbonylcyanide-m-chlorophenylhydrazone (CCCP), rotenone, and diphenylene iodonium chloride (DPI) significantly reduced ROS levels. The elevated ROS were correlated to an increased expression of PrPc, Cu/Zn superoxide dismutase (SOD-1), and catalase in small as compared with large spheroids. In large tumor spheroids, PrPc was predominantly expressed in the peripheral cell layers and colocalized with SOD-1 and catalase. Raising intracellular ROS in large tumor spheroids by hydrogen peroxide, menadione, buthionine sulfoximine (BSO), and incubation in glutamine-reduced medium increased PrPc expression. In small spheroids PrPc was downregulated after incubation with the radical scavengers dehydroascorbate (DHA) and vitamin E. Our data indicate that PrPc expression in tumor spheroids is related to the intracellular redox state and may participate in antioxidative defense.     

Growth and characterization of multicellular tumor spheroids of human bladder carcinoma origin

Summary We have examined the MGH-U1 human bladder carcinoma cell line and 12 primary bladder carcinoma biopsies for their ability to form spheroids in suspension culture and in multiwell dishes. MGH-U1 cells formed tightly packed spheroids with a necrotic center and viable rim whereas three sublines formed loose aggregates only. Spheroids formed from as few as 100 MGU-U1 cells placed into multiwells. MGH-U1 cells derived from spheroids formed new spheroids more rapidly and consistently than cells derived from monolayer culture. Spheroid diameter increased at a rapid rate of ∼100 μm/d in multiwell dishes, and necrosis occurred only in spheroids of diameter >1 mm. Spheroids placed in spinner culture at a higher concentration (∼1.5 spheroids/ml) grew more slowly and developed necrosis at smaller diameters. The width of the viable rim of spheroids grown in spinner culture was maintained at ∼190 μm over a wide range of spheroid diameters (400 to 1000 μm). Sequential trypsinization of spheroids, which stripped layers of cells from the spheroids, demonstrated no difference in the plating efficiency of cells derived from varying depths into the spheroid. Only one of the 12 primary bladder biopsy specimens demonstrated an ability to form spheroids. This biopsy, designated HB-10, formed spheroids that grew linearly over 40 d, formed colonies in methylcellulose culture and grew as xenografts in immune-deprived mice. These studies characterize the MGH-U1 spheroids that are useful in vitro models to study the effects of various treatments for solid tumors and demonstrate the limited capacity of cells from primary human bladder biopsies to form spheroids. Supported in part by a grant from the National Cancer Institute of Canada and by grant CA29526 NCI through the National Bladder Cancer Project, U.S.A.     

Penetration of anti-melanoma immunotoxin into multicellular tumor spheroids and cell kill effects

Summary In order to gain a better understanding of the interaction between immunotoxins and tumor cells at the level of three-dimensional tumor mass, we evaluated the cell kill effects of monoclonal antimelanoma-antibody/ricin-A-chain immunotoxin (ITN) on melanoma cells in multicellular tumor spheroids (MTS) as well as the penetration of ITN into MTS. For Minor melanoma cells in monolayer the ITN exerted cytotoxic effects after as little as 1 h of exposure. Increasing exposure time resulted in progressive increases in cytotoxic activity. In contrast, the cell kill effects of ITN were markedly delayed and reduced when Minor cells were in MTS. The ITN cytotoxic effects on the melanoma MTS were more than 100 fold less than those in monolayer. Patterns of ITN-induced cytotoxicities for Minor and for another melanoma cell line, DND-1A, were comparable. The native ricin A was more active against PC-10 squamous lung cancer cells than Minor cells, whereas the ITN was more cytotoxic against Minor cells than PC-10 cells, thus exhibiting selectivity. An autoradiographic study revealed time-dependent penetration of radiolabeled ITN from the surface of Minor MTS into the core. Incubation for 1 h resulted in the penetration of ITN into only the two or three outer layers of the Minor MTS, and low grain counts. Prolonged exposure resulted in inhomogeneous penetration of ITN into almost the entire melanoma MTS. Penetration of ITN into PC-10 MTS was extremely poor. The reduced cytotoxicity of ITN on melanoma cells in MTS as compared to cells grown in monolayer appears to correlate with its inhomogeneous distribution in the MTS. The delayed cytotoxicity of ITN is also consistent with its slow penetration into the core of the MTS.