Joint oxygen-glucose deprivation as the cause of necrosis in a tumor analog

The sandwich system was recently developed as an in vitro tumor analog. Like spheroids, sandwiches are organized, multicellular systems in which the interplay between diffusion and consumption leads to the formation of spatial gradients; a necrotic center and a viable cell border subsequently develop. Using sandwiches of the 9L and V79 cell lines, the effects of oxygen and glucose deprivation on the onset and formation of necrosis were investigated. The data indicate that in sandwiches necrosis is a result of a shortage of both substances. Complementary cell monolayer experiments to determine a number of consumption parameters were performed. On the basis of the data, we propose a joint oxygen-glucose deprivation model for V79 cell necrosis. It is assumed a cell dies when oxygen deprivation in conjunction with glucose deprivation lowers the cell's ATP production rate below a critical value. Interactions of the concentrations and consumptions of oxygen and glucose are analyzed theoretically; concentration profiles are obtained by numerically solving coupled non-linear integral equations arising from the diffusion equation. The predicted viable border widths are in good agreement with the observed values.     

Differences in the X-ray sensitivity of cells in different regions of the sandwich, a diffusion-limited system for cell growth

The sandwich system was recently developed as a tumor analog; like spheroids, sandwiches are diffusion-limited multicellular systems which exhibit a necrotic center and a viable cell border. Using sandwiches of the 9L cell line, we compared the X-ray sensitivity of cells in the inner half of the viable border, adjacent to the necrotic center, with that of cells in the outer half of the viable border, adjacent ot the medium. No cells were hypoxic at the time of irradiation. The cells in the inner half of the viable border exhibited an increased radioresistance over cells in the outer half. The effect was dose multiplying with a multiplying factor of 1.5. Besides the sandwich studies, the X-ray sensitivity of 9L plateau monolayer cultures (induced by starvation) was compared to exponentially growing monolayer cultures. The plateau cultures exhibited an increased radioresistance over the exponentially growing cultures. The effect was also dose multiplying with a multiplying factor of 1.5.     

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.     

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.     

ATP concentrations in multicellular tumor spheroids assessed by single photon imaging and quantitative bioluminescence

To evaluate the interrelationship among the cellular energy status and the development of necrosis in tumor microregions, local ATP concentrations and the extent of necrosis were determined in multicellular tumor spheroids, i.e., in spherical tumor cell aggregates. The spheroids were grown in rotated suspension cultures using EMT6 cells that were derived from a murine mammary sarcoma. The distribution of viable and necrotic cell areas was assessed by histological investigations. The regional distribution of ATP concentrations was measured with a novel technique using quantitative bioluminescence and single photon imaging. This method makes it possible to determine ATP concentrations in absolute terms with a spatial resolution at the level of a single cell. The results show that ATP concentrations in the center of EMT6 spheroids decrease from values of 1.0 to 1.5 mM in small spheroids with 300 microns in diameter to values close to or at the background level in 750 microns spheroids. Necrosis was detectable in spheroids larger than 300 microns, and virtually no spheroid without necrosis was found at sizes larger than 600 microns. Since the emergence of central necrosis precedes the drop in ATP to undetectably low values, the data suggest that energy metabolism is not or not directly involved in the development of necrosis in tumor spheroids under the growth conditions investigated.     

Cytophotometric measurement of the cellular DNA content of [3H]thymidine-labelled spheroids. Demonstration that some non-labelled cells have S and G2 DNA content

Spheroids from the V279-171b and MCa-11 cell lines were incubated continuously for 24 hr in [3H]thymidine for labelling of the outer cells of the viable rim. The spheroids were dispersed into single cells, and the DNA content of photomapped cells was measured by absorption cytophotometry. Autoradiographs were then prepared from which we ascertained cellular labelling. For spheroids of both cell lines, we found a larger proportion of cells with a G0/G1 DNA content among the non-labelled inner spheroid cells than among the labelled outer cells (P less than 0.001). This block of non-labelled spheroid cells in G0/G1 was not a cell cycle perturbation caused by the isotope for the MCa-11 spheroids. Approximately 8% of non labelled MCa-11 spheroid cells had S/G2 DNA content, suggesting that non-cycling cells in spheroids may be blocked in S and G2 as well as in the G0/G1 phase of the cell cycle.     

Modeling mass transfer in hepatocyte spheroids via cell viability, spheroid size, and hepatocellular functions

Hepatocyte aggregation into spheroids attributes to their increased activity, but in the absence of a vascular network the cells in large spheroids experience mass transfer limitations. Thus, there is a need to define the spheroid size which enables maximal cell viability and productivity. We developed a combined theoretical and experimental approach to define this optimal spheroid size. Hepatocyte spheroids were formed in alginate scaffolds having a pore diameter of 100 microm, in rotating T-flasks or spinners, to yield a maximal size of 100, 200, and 600 microm, respectively. Cell viability was found to decrease with increasing spheroid size. A mathematical model was constructed to describe the relationship between spheroid size and cell viability via the oxygen mass balance equation. This enabled the prediction of oxygen distribution profiles and distribution of viable cells in spheroids with varying size. The model describes that no oxygen limitation will take place in spheroids up to 100 microm in diameter. Spheroid size affected the specific rate of albumin secretion as well; it reached a maximal level, i.e., 60 microg/million cells/day in 100-microm diameter spheroids. This behavior was depicted in an equation relating the specific albumin secretion rate to spheroid size. The calculated results fitted with the experimental data, predicting the need for a critical number of viable hepatocytes to gain a maximal albumin secretion. Taken together, the results on mass transport in spheroids and its effects on cell viability and productivity provide a useful tool for the design of 3D scaffolds with pore diameters of 100 microm.     

In situ oxygen consumption rates of cells in V-79 multicellular spheroids during growth

The rate of consumption of oxygen by V-79 cells in multicellular spheroids was measured as a function of the spheroid diameter. In situ consumption was equal to that of exponentially growing cells for spheroids less than 200 micron in diameter. The rate of oxygen consumption decreased for cells in spheroids between 200 and 400 micron diameter to a value one-fourth the initial, then remained constant with further spheroid growth. Comparison of consumption rates for spheroid-derived cells before and after dissociation from the spheroid structure indicated that the spheroid microenvironment accounted for only 20% of the change in oxygen consumption rate. Cell-cell contact, cell packing, and cell volume were not critical parameters. Plateau-phase cells had a fivefold lower rate of oxygen consumption than exponential cells, and it is postulated that the spheroid quiescent cell population accounts for a large part of the intrinsic alteration in oxygen consumption of cells in spheroids. Some other mechanism must be involved in the regulation of cellular oxygen consumption in V-79 spheroids to account for the remainder of the reduction observed in this system.     

Further characterization of 4-bromomisonidazole as a potential detector of hypoxic cells

[14C]Bromomisonidazole was prepared by direct bromination of [ring-2] [14C]misonidazole in dioxane. The uptake and binding of the two labeled sensitizers were compared in vitro in 1-mm EMT-6 spheroids which contain a necrotic core. Using liquid scintillation counting it was shown that spheroids incubated with 50 microM [14C]bromomisonidazole concentrated drug above levels in the medium by 1 1/2 hr and achieved maximum concentration by 10 hr with no further increase at 23 hr. Spheroids incubated with 50 microM [14C]misonidazole may concentrate the sensitizer more slowly but ultimately reached the same fivefold increase over levels in the medium by 23 hr as was observed for bromomisonidazole. Autoradiographs prepared from spheroids after incubation with [14C]misonidazole or [14C]bromomisonidazole showed silver grains preferentially located over viable hypoxic cells in the inner half of the spheroid rim adjacent to the necrotic center, with lower grain density over nonviable necrotic areas and many fewer grains over oxic cells at the periphery of the spheroid. The results indicate that both severely and moderately hypoxic cells may preferentially bind [14C]bromomisondiazole. The data support the potential of radiolabeled bromomisonidazole for in vivo imaging pending additional studies of the metabolism of this agent.     

Modeling of self-organized avascular tumor growth with a hybrid cellular automaton

Pattern formation in multicellular spheroids is addressed with a hybrid lattice-gas cellular automaton model. Multicellular spheroids serve as experimental model system for the study of avascular tumor growth. Typically, multicellular spheroids consist of a necrotic core surrounded by rings of quiescent and proliferating tumor cells, respectively. Furthermore, after an initial exponential growth phase further spheroid growth is significantly slowed down even if further nutrient is supplied. The cellular automaton model explicitly takes into account mitosis, apoptosis and necrosis as well as nutrient consumption and a diffusible signal that is emitted by cells becoming necrotic. All cells follow identical interaction rules. The necrotic signal induces a chemotactic migration of tumor cells towards maximal signal concentrations. Starting from a small number of tumor cells automaton simulations exhibit the self-organized formation of a layered structure consisting of a necrotic core, a ring of quiescent tumor cells and a thin outer ring of proliferating tumor cells.     

A reduction in the in situ rates of oxygen and glucose consumption of cells in EMT6/Ro spheroids during growth

The rates of consumption of oxygen and glucose by EMT6/Ro cells in multicellular spheroids were measured at various times during normal growth. In situ spheroid cellular consumption rates were similar to those of exponentially growing single cells up to a spheroid diameter of 150 micron. Further growth resulted in decreases in the rates of both oxygen and glucose consumption which were correlated with the increase in spheroid diameter and cell number. At a diameter of 1300 micron, both rates of cellular consumption had decreased by a factor of 2.5. The rates of consumption per unit of nonnecrotic spheroid volume decreased in a similar manner. Measurements with single cells demonstrated that the rate of oxygen consumption was coupled with glucose concentration, and vice versa. The rates of consumption for cells dissociated from small spheroids indicated that there was some effect of the spheroid environment. As the spheroids grew, however, association in the spheroid structure accounted for a smaller proportion of the total observed reduction in the rates of nutrient consumption. The presence of central necrosis also appeared to have no effect on the rates of consumption of these nutrients. Spheroid-derived cells showed a decrease in cell volume with growth as the cells accumulated in a quiescent state. Measurements with single cells demonstrated that oxygen and glucose consumption were correlated with cell volume and with the development of nonproliferating cells. We conclude that the observed decrease in oxygen and glucose consumption with growth in spheroids is largely due to the progressive accumulation of cells in a quiescent state characterized by an inherently lower cellular rate of nutrient utilization.     

Necrotic core in EMT6/Ro tumour spheroids: Is it caused by an ATP deficit?

Although commonly related to nutrient deprivation, the cause of the formation of the necrotic core in the multicellular tumour spheroids is still a controversial issue. We propose a simple model for the cell ATP production that assumes glucose and lactate as the only fuel substrates, and describes the main reactions occurring in the glycolytic and the oxidative pathways. Under the key assumption that cell death occurs when ATP production falls to a critical level, we formulate a multiscale model that integrates the energy metabolism at the cellular level with the diffusive transport of the metabolites in the spheroid mass. The model has been tested by predicting the measurements of the necrotic radius obtained by Freyer and Sutherland (1986a) in EMT6/Ro spheroids under different concentrations of glucose and oxygen in the culture medium. The results appear to be in agreement with the hypothesis that necrosis is caused by ATP deficit.     

Evidence for a secreted chemorepellent that directs glioma cell invasion

Secreted chemotropic cues guide the migration of neuronal and glial cell precursors during neural development. It is not known if chemotropism contributes to directing the invasion of brain tissue by glioma cells. A model system has been developed that allows quantification of invasive behavior using gliomas spheroids embedded in collagen gels. Here we provide evidence that glioma spheroids secrete a chemorepellent factor(s) that directs cells away from the spheroid and into the collagen matrix. The relationship between total invasion, cell number, and implantation distance suggests that glioma cells respond to a gradient of the chemorepellent cue(s) that is well established at 48 h. C6 astrocytoma cells normally invade the collagen at an angle perpendicular to the spheroid edge. In contrast, an adjacent spheroid causes cells to turn away from their normal trajectory and slow their rate of invasion. Astrocytoma cells are repelled by an adjacent glioma spheroid but rapidly infiltrate astrocyte aggregates, indicating that astrocytes do not express the repellent cue. Uniform concentrations of repellent factor(s) in spheroid conditioned medium overwhelm endogenous gradients and render glioma cells less able to exhibit this chemotropic response. Concentration gradients of spheroid conditioned medium in cell migration assays also demonstrate the chemorepellent cue(s)'s tropic effect. Our findings indicate that glioma spheroids produce a secreted diffusible cue(s) that promotes glioma cell invasion. Identification of this factor(s) may advance current therapies that aim to limit tumor cell invasion.     

Cell migration in multicell spheroids: Swimming against the tide

Multicell spheroids, small spherical clusters of cancer cells, have become an importantin vitro model for studying tumour development given the diffusion limited geometry associated with many solid tumour growths. Spheroids expand until they reach a dormant state where they exhibit a grossly static three-layered structure. However, at a cellular level, the spheroid is demonstrably dynamic with constituent cells migrating from the outer well-nourished region of the spheroid toward the necrotic central core. The mechanism that drives the migrating cells in the spheroid is not well understood. In this paper we demonstrate that recent experiments on internationalization can be adequately described by implicating pressure gradients caused by differential cell proliferation and cell death as the primary mechanism. Although chemotaxis plays a role in cell movement, we argue that it acts against the passive movement caused by pressure differences.     

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.     

Formulation and numerical simulations of a continuum model of avascular tumor growth

In this paper we present a continuum mathematical model for a multicellular spheroid that mimics the micro-environment within avascular tumor growth. The model consists of a coupled system of non-linear convection-diffusion-reaction equations. This system is solved using a previously developed conservative Galerkin characteristics method. In the model considered, there are three cell types: the proliferative cells, the quiescent non-dividing cells which stay in the G₀ phase of the cell cycle and the necrotic cells. The model includes viable cell diffusion, diffusion of cellular material and the removal of necrotic cells. We assume that the nutrients diffuse passively and are consumed by the proliferative and quiescent tumor cells depending on the availability of resources (oxygen, glucose, etc.). The numerical simulations are performed using different sets of parameters, including biologically realistic ones, to explore the effects of each of these model parameters on reaching the steady state. The present results, taken together with those reported earlier, indicate that the removal of necrotic cells and the diffusion of cellular material have significant effects on the steady state, reflecting growth saturation, the number of viable cells, and the spheroid size.     

Imaging- and Flow Cytometry-based Analysis of Cell Position and the Cell Cycle in 3D Melanoma Spheroids

Three-dimensional (3D) tumor spheroids are utilized in cancer research as a more accurate model of the in vivo tumor microenvironment, compared to traditional two-dimensional (2D) cell culture. The spheroid model is able to mimic the effects of cell-cell interaction, hypoxia and nutrient deprivation, and drug penetration. One characteristic of this model is the development of a necrotic core, surrounded by a ring of G1 arrested cells, with proliferating cells on the outer layers of the spheroid. Of interest in the cancer field is how different regions of the spheroid respond to drug therapies as well as genetic or environmental manipulation. We describe here the use of the fluorescence ubiquitination cell cycle indicator (FUCCI) system along with cytometry and image analysis using commercial software to characterize the cell cycle status of cells with respect to their position inside melanoma spheroids. These methods may be used to track changes in cell cycle status, gene/protein expression or cell viability in different sub-regions of tumor spheroids over time and under different conditions.     

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.     

Automated selective dissociation of cells from different regions of multicellular spheroids

Summary In this report we describe a new apparatus which has been developed for the automated selective dissociation of multicellular spheroids into fractions of viable cells from different locations in the spheroid. This device is based on the exposure of spheroids to a 0.25% solution of trypsin under carefully controlled conditions, such that the cells are released from the outer spheroid surface in successive layers. Study of the spheroid size, number of cells per spheroid, and sections through the spheroid with increasing exposure to trypsin demonstrate the effectiveness of this technique. The technique has been successfully used on spheroids from five different cell lines over a wide range of spheroid diameters. We also present data detailing the effect of varying the dissociation temperature, the mixing speed, the trypsin concentration, and the number of spheroids being dissociated. The new apparatus has several advantages over previous selective dissociation methods and other techniques for isolating cells from different regions in spheroids, including: a) precise control over dissociation conditions, improving reproducibility; b) short time to recover cell fractions; c) ability to isolate large numbers of cells from many different spheroid locations; d) use of common, inexpensive laboratory equipment; and e) easy adaptability to new cell lines or various spheroid sizes. Applications of this method are demonstrated, including the measurement of nutrient consumption rates, regrowth kinetics, and radiation survivals of cells from different spheroid regions. This work was supported by grants CA-36535, CA-22585, and RR-02845 from the National Institutes of Health, Bethesda, MD, the National Flow Cytometry Resource (NIH grant RR-01315), and by the Department of Energy, Washington, DC.     

Optimization of the formation of embedded multicellular spheroids of MCF-7 cells: How to reliably produce a biomimetic 3D model

To obtain a multicellular MCF-7 spheroid model to mimic the three-dimensional (3D) of tumors, the microwell liquid overlay (A) and hanging-drop/agar (B) methods were first compared for their technical parameters. Then a method for embedding spheroids within collagen was optimized. For method A, centrifugation assisted cells form irregular aggregates but not spheroids. For method B, an extended sedimentation period of over 24 h for cell suspensions and increased viscosity of the culture medium using methylcellulose were necessary to harvest a dense and regular cell spheroid. When the number was less than 5000 cells/drop, embedded spheroids showed no tight cores and higher viability than the unembedded. However, above 5000 cells/drop, cellular viability of embedded spheroids was not significantly different from unembedded spheroids and cells invading through the collagen were in a sun-burst pattern with tight cores. Propidium Iodide staining indicated that spheroids had necrotic cores. The doxorubicin cytotoxicity demonstrated that spheroids were less susceptible to DOX than their monolayer cells. A reliable and reproducible method for embedding spheroids using the hanging-drop/agarose method within collagen is described herein. The cell culture model can be used to guide experimental manipulation of 3D cell cultures and to evaluate anticancer drug efficacy.