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.     

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.     

Gene transfection of multicellular spheroid of hepatocytes on an artificial substrate

The handling of hepatocytes, a major cell population in the liver, is an important technique in both liver tissue engineering and hepatology. However, these cells are so fragile that it has been impossible to harvest hepatocytes with high viability from tissue culture dishes after a period of culture in vitro. In this study, we employed an artificial substrate for transfection of multilayer hepatocytes and harvested these cells with high viability after transfection. Hepatocytes cultured on an amphiphilic artificial substrate form multilayer aggregates (spheroids) in the presence of growth factors during gene transfection with cation liposomes. Compared to cells cultured on a collagen-coated plate, these spheroids are easily harvested with high viability by pipetting in EDTA solution. In addition, these spheroids rapidly spread on collagen after transfer from the artificial substrate, demonstrating that hepatocytes in the center of the spheroids were viable. Epidermal growth factor (EGF) increased the transfection efficiency into hepatocytes while hepatocyte growth factor (HGF) alone did not increase the efficiency. However, HGF synergestically increased the effect of EGF on transfection. Interestingly, this transfection required the process of spheroid formation because the gene was not transfected once the spheroid formation completed or under conditions where hepatocytes did not form spheroids. This method using spheroidal hepatocytes for in vitro transfection is promising for the development of ex vivo gene therapy.     

Formation of multicellular spheroids composed of adult rat hepatocytes in dishes with positively charged surfaces and under other nonadherent environments

Adult rat hepatocytes formed floating multicellular spheroids in primary culture in an uncoated plastic dish with a positively charged surface. Cells in the spheroids formed in such a simple way were similar to those formed in dishes coated with proteoglycan fraction isolated from rat liver reticulin fibers; in both cases, cells maintained high ability to produce albumin and poor ability to proliferate in response to epidermal growth factor. Coating dishes with albumin was also helpful in spheroid formation; coating with 2-hydroxymethyl methacrylate resulted in formation of incomplete spheroids. Elimination of serum factors was essential for the formation of spheroids; when cells were washed with serum-containing medium before seeding or if the medium was replaced with a serum-containing medium, spheroid formation was completely inhibited. Collagens, fibronectin, and laminin, all of which promote the adhesion and spreading of hepatocytes on substrates, inhibited spheroid formation. Furthermore, collagens disintegrated spheroids, and cells in the monolayer initiated proliferation. Thus, two distinct, mutually exclusive features of primary culture of adult hepatocytes apparently exist; monolayer culture with proliferative activity in an adherent environment and spheroid culture with poor proliferative activity and high albumin-producing ability in a nonadherent environment.     

Growth of cancer cell lines under stem cell-like conditions has the potential to unveil therapeutic targets

Malignant tumors comprise a small proportion of cancer-initiating cells (CIC), capable of sustaining tumor formation and growth. CIC are the main potential target for anticancer therapy. However, the identification of molecular therapeutic targets in CIC isolated from primary tumors is an extremely difficult task. Here, we show that after years of passaging under differentiating conditions, glioblastoma, mammary carcinoma, and melanoma cell lines contained a fraction of cells capable of forming spheroids upon in vitro growth under stem cell-like conditions. We found an increased expression of surface markers associated with the stem cell phenotype and of oncogenes in cell lines and clones cultured as spheroids vs. adherent cultures. Also, spheroid-forming cells displayed increased tumorigenicity and an altered pattern of chemosensitivity. Interestingly, also from single retrovirally marked clones, it was possible to isolate cells able to grow as spheroids and associated with increased tumorigenicity. Our findings indicate that short-term selection and propagation of CIC as spheroid cultures from established cancer cell lines, coupled with gene expression profiling, represents a suitable tool to study and therapeutically target CIC: the notion of which genes have been down-regulated during growth under differentiating conditions will help find CIC-associated therapeutic targets.     

Novel fluid shear‐based dissociation device for improved single cell dissociation of spheroids and cell aggregates

Biological industries commonly rely on bioreactor systems for the large‐scale production of cells. Cell aggregation, clumping, and spheroid morphology of certain suspension cells make their large‐scale culture challenging. Growing stem cells as spheroids is indispensable to retain their stemness, but large spheroids (>500 µm diameter) suffer from poor oxygen and nutrient diffusion, ultimately resulting in premature cell death in the centers of the spheroids. Despite this, most large‐scale bioprocesses do not have an efficient method for dissociating cells into single cells, but rely on costly enzymatic dissociation techniques. Therefore, we tested a proof‐of‐concept fluid shear‐based mechanical dissociator that was designed to dissociate stem cell spheroids and aggregates. Our prototype was able to dissociate cells while retaining high viability and low levels of apoptosis. The dissociator also did not impact long‐term cell growth or spheroid formation. Thus, the dissociator introduced here has the potential to replace traditional dissociation methods. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 34:293–298, 2018     

Fibroblast EXT1-Levels Influence Tumor Cell Proliferation and Migration in Composite Spheroids

Background

Stromal fibroblasts are important determinants of tumor cell behavior. They act to condition the tumor microenvironment, influence tumor growth, support tumor angiogenesis and affect tumor metastasis. Heparan sulfate proteoglycans, present both on tumor and stromal cells, interact with a large number of ligands including growth factors, their receptors, and structural components of the extracellular matrix. Being ubiquitously expressed in the tumor microenvironment heparan sulfate proteoglycans are candidates for playing central roles in tumor-stroma interactions. The objective of this work was to investigate the role of heparan sulfate expressed by stromal fibroblasts in modulating the growth of tumor cells and in controlling the interstitial fluid pressure in a 3-D model.

Methodology/Principal Findings

We generated spheroids composed of fibroblasts alone, or composite spheroids, composed of fibroblasts and tumor cells. Here we show that stromal fibroblasts with a mutation in the heparan sulfate elongating enzyme Ext1 and thus a low heparan sulfate content, formed composite fibroblast/tumor cell spheroids with a significant lower interstitial fluid pressure than corresponding wild-type fibroblast/tumor cell composite spheroids. Furthermore, immunohistochemistry of composite spheroids revealed that the cells segregated, so that after 6 days in culture, the wild-type fibroblasts formed an inner core and the tumor cells an outer layer of cells. For composite spheroids containing Ext1-mutated fibroblasts this segregation was less obvious, indicating impaired cell migration. Analysis of tumor cells expressing the firefly luciferase gene revealed that the changes in tumor cell migration in mutant fibroblast/tumor cell composite spheroids coincided with a lower proliferation rate.

Conclusions/Significance

This is the first demonstration that stromal Ext1-levels modulate tumor cell proliferation and affect the interstitial fluid pressure in a 3-D spheroid model. Learning how structural changes in stromal heparan sulfate influence tumor cells is essential for our understanding how non-malignant cells of the tumor microenvironment influence tumor cell progression.     

Partial reversion of the phenotype of a poorly differentiated hepatocellular carcinoma in a three-dimensional culture

In vivo, normal tissues and organs have a three-dimensional structure and function in a three-dimensional environment. The standard two-dimensional cell culture conditions drastically differ from those in vivo. For this reason, three-dimensional cultures based on different variants of the extracellular matrix are more adequate for analyzing normal and tumor cell growth. Culturing a poorly differentiated hepatocellular carcinoma in a collagen gel yielded spheroids whose growth pattern shifted towards the epithelial phenotype. The shift was expressed in changes in the cytoskeleton, enhanced formation of extracellular matrix fibrils between cells, and formation of fibronectin fibrils on the outer surface of spheroids. Analysis of 25 genes reflecting the level of morphological and functional hepatocyte differentiation showed that the expression of the gene encoding the transforming growth factor TGFβ2 was suppressed the most significantly.     

Biochemical and functional changes of rat liver spheroids during spheroid formation and maintenance in culture: I. morphological maturation and kinetic changes of energy metabolism, albumin synthesis, and activities of some enzymes

In the process of isolated single liver cells coming together to form three-dimensional spheroids, cells undergo dramatic environmental changes. How liver cells respond to these changes has not been well studied before. This study characterized the functional and biochemical changes during liver spheroid formation and maintenance. Spheroids were prepared in 6-well plates from freshly isolated liver cells from male Sprague rats by a gyrotatory-mediated method. Morphological formation, and functional and biochemical parameters of liver spheroids were evaluated over a period of 21 days in culture. Liver spheroid formation was divided into two stages, immature (1-5 days) and mature (>5 days), according to their size and shape, and changes in their functionality. Galactose and pyruvate consumption was maintained at a relatively stable level throughout the period of observation. However, glucose secretion and cellular GPT and GOT activities were higher in immature spheroids, decreased upto day 5 and remained stable thereafter. Cellular gamma-glutamyltransferase (gamma-GT) and lactate dehydrogenase (LDH) activities were initially undetectable or low and increased as spheroids matured. Albumin secretion decreased rapidly within the first 2 days and increased as spheroids matured. It is concluded that cells undergo functional and biochemical changes during spheroid formation following isolation of liver cells from intact tissue. Functionality and biochemical properties recovered and were maintained in mature spheroids. A relatively stable period (6-15 days) of functionality in mature spheroids was identified and is recommended for applications of the model.     

Development of intracytoplasmic lumens in a colon cancer cell line cultured on a non-adhesive surface

Cell-matrix interactions have important effects on phenotypic features, such as morphology, differentiation and cell growth. Several papers have suggested that when cell-matrix interactions are interrupted, cells grow as multicellular spheroids and eventually undergo apoptosis. We found that when ET(-), a laminin non-adherent colon cancer cell line, was cultured on poly-2-hydroxyethyl methacrylate (HEMA) coated plastic, the cells floated as cellular aggregates of spheroids or as single cells. Some of the single cells contained a very large intracytoplasmic lumen (ICL) and appeared similar to signet ring cells. These ICL were lined by a layer of short microvilli. The number of the cell did not increased cells when cultured on poly-HEMA. Another type of single cells, usually without ICL, demonstrated the characteristics of apoptotic cells by histologic examination. Acridine orange staining, flow cytometry and electron microscopy confirmed the apoptotic nature of those cells. In immunohistochemical staining for proliferating cell nuclear antigen, spheroids of cells and single cells with ICL were immunoreactive, while most of the single cells without ICL were negative. These results suggest that multicellular aggregation and formation of ICL were induced by the adaptation of ET(-) colon cancer cells in a harmful environment caused by reduced adhesiveness, and these changes might be related to cell survival.     

Formation and characterization of three dimensional human hepatocyte cell line spheroids on chitosan matrix for in vitro tissue engineering applications

Chitosan was used as a matrix to induce three-dimensional spheroids of HepG2 cells. Chitosan films were prepared and used for culturing Hep G2 cells. Attachment kinetics of the cells was studied on the chitosan films. The optimum seeding density of the Hep G2 cells, required for three-dimensional spheroid formation was determined and was found to be 5 × 10⁴/ml. The growth kinetics of Hep G2 cells was studied using (3-(4, 5-Dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide) (MTT) assay, and morphology of the cells was studied through optical photographs taken at various days of culture. The liver cell functions of the spheroids were determined by measuring albumin and urea secretions. The results obtained from these studies have shown that the culture of Hep G2 cells on chitosan matrix taking appropriate seeding density resulted in the formation of three-dimensional spheroids and exhibited higher amount of albumin and urea synthesis compared to monolayer culture. These miniature “liver tissue like” models can be used for in vitro tissue engineering applications like preliminary evaluation of the toxicity of drugs and chemicals.     

Characterization of a fetal urogenital sinus mesenchymal cell line U4F: Secretion of a negative growth regulatory activity

Summary Mesenchymal cell lines derived from fetal rat urogenital sinus organ cultures have been characterized to establish an in vitro system for addressing growth and differentiation regulatory factors involved in mesenchymal-epithelial interactions during prostate morphogenesis. A continuous cell line was developed and designated U4F. Immunocytochemical analysis showed vimentin intermediate filament content confirming a mesenchymal origin. Previous studies with urogenital sinus organ cultures have reported the expression of a negative growth activity, which is stimulatory to protein synthesis and secretion and alters phenotypic morphology of NBT-II bladder epithelial cells. Subconfluent and confluent U4F monolayers did not produce this growth inhibitory activity. Foci of stacked cells were observed 3 wk postconfluency, which evolved into multicellular spheroids. The negative growth activity was expressed in the conditioned medium coordinate with spheroid formation. Transplanted spheroids continued to express the growth inhibitory activity. Morphologic analysis of spheroids showed a cellular capsule and a core of extracellular matrix. A continuous cell strain (U4F1) with altered phenotypic properties, arose spontaneously from long-term U4F cultures. The U4F1 cell strain did not form spheroids, yet expressed the negative growth activity constitutively in monolayer culture. Analyses of physicochemical, immunological, and biological properties showed the activity is identical in conditioned media from urogenital sinus organ cultures, U4F spheroids, and U4F1 monolayers. Based on the combined properties, this activity cannot be ascribed to previously characterized negative growth factors. The establishment of this mesenchymal cell culture system will aid in the further identification of paracrine-acting growth and differentiation regulatory factors secreted by fetal mesenchyme.     

Aggregation of human dental pulp cells into 3D spheroids enhances their migration ability after reseeding

Multicellular three-dimensional (3D) spheroids allow intimate cell–cell communication and cell–extracellular matrix interaction. Thus, 3D cell spheroids better mimic microenvironment in vivo than two-dimensional (2D) monolayer cultures. The purpose of this study was to evaluate the behaviors of human dental pulp cells (DPCs) cultured on chitosan and polyvinyl alcohol (PVA) membranes. The protein expression of hypoxia-inducible factor 1-α (HIF-1α) and vascular endothelial growth factor (VEGF), and the migration ability of the DPCs from 2D versus 3D environments were investigated. The results showed that both chitosan and PVA membranes support DPCs aggregation to form multicellular spheroids. In comparison to 2D cultures on tissue culture polystyrene, DPC spheroids exhibited higher protein expression of HIF-1α and VEGF. The treatment with YC-1 (inhibitor to HIF-1α) blocked the upregulation of VEGF, indicating a downstream event to HIF-1α expression. When DPC spheroids were collected and subjected to the transwell assay, the cells growing outward from 3D spheroids showed greater migration ability than those from 2D cultures. Moreover, DPCs aggregation and spheroid formation on chitosan membrane were abolished by Y-27632 (inhibitor to Rho-associated kinases), whereas the inhibitory effect did not exist on PVA membrane. This suggests that the mechanism regulating DPCs aggregation and spheroid formation on chitosan membrane is involved with the Rho-associated kinase signaling pathway. In summary, the multicellular spheroid structure was beneficial to the protein expression of HIF-1α and VEGF in DPCs and enhanced the migration ability of the cells climbing from spheroids. This study showed a new perspective in exploring novel strategies for DPC-based research and application.     

Recent advances in three-dimensional multicellular spheroid culture for biomedical research

Many types of mammalian cells can aggregate and differentiate into 3-D multicellular spheroids when cultured in suspension or a nonadhesive environment. Compared to conventional monolayer cultures, multicellular spheroids resemble real tissues better in terms of structural and functional properties. Multicellular spheroids formed by transformed cells are widely used as avascular tumor models for metastasis and invasion research and for therapeutic screening. Many primary or progenitor cells on the other hand, show significantly enhanced viability and functional performance when grown as spheroids. Multicellular spheroids in this aspect are ideal building units for tissue reconstruction. Here we review the current understanding of multicellular spheroid formation mechanisms, their biomedical applications, and recent advances in spheroid culture, manipulation, and analysis techniques.     

Multicell spheroids as a model for cell kinetic studies

Cells growing in tissue culture as three-dimensional, multicellular aggregates called 'spheroids' typically show a decreasing growth fraction and development of quiescent subpopulations as the spheroids enlarge. Kinetic studies in a number of spheroid systems have indicated that the primary reason for the tumour-like growth is a progressive decrease in growth fraction, with only a modest elongation of cell cycle time in larger spheroids. In this paper, the cellular growth kinetics for spheroids of V79 Chinese hamster lung cells are reviewed, and the regrowth kinetics of cells resuming growth after recovery from quiescent regions of the spheroids are described. Further, the role of regrowth/repopulation in determining the spheroid response to anti-tumour cytotoxics is explored, with particular emphasis on treatment with cisplatin and etoposide. By separating the effects of cytotoxicity and regrowth in the overall spheroid response to anti-neoplastic drugs, it is suggested that 'drug resistance' in tumours can be a kinetic as well as a genetic problem.     

Serum, trypsin, and cell shape but not cell-to-cell contact influence the X-ray sensitivity of Chinese hamster V79 cells in monolayers and in spheroids

Nutrient concentration in the growth medium and trypsin affect cellular radiosensitivity in a manner that is related to cell shape (Reddy, Stevenson, and Lange, Int. J. Radiat. Biol. 55, 105-117 (1989); Reddy and Lange, Radiat. Res. 119, 338-347 (1989]. Hence we hypothesized that the concentration of serum in the medium could influence the X-ray sensitivity of cells and that the spread cells in monolayers and round cells in spheroids may differ in their response to the radiosensitizing effect of trypsin. We compared the X-ray sensitivity of monolayer and spheroid cells grown for 19 +/- 1 h in MEM supplemented with 5 or 15% serum. Cells were trypsinized and plated either immediately before, or 2.5 +/- 0.5 h after, irradiation and incubation for repair in situ. Survival of cells in monolayers and in spheroids was higher in MEM with 5% serum than with 15% serum. Trypsin treatment affected the shape and radiosensitivity of cells in monolayers but not in spheroids. When all cells were grown in the same serum concentration and a 2.5-h postirradiation incubation was allowed prior to trypsinization, the X-ray sensitivity of cells in spheroids was greater than that of cells in monolayers. The survival of cells in spheroids became equal to that of monolayer cells when cells in spheroids were converted to monolayers by placing them in 25-cm2 flasks and allowing them 3 h to attach and spread. Cell cycle distributions were nearly the same in monolayers and spheroids cultured in MEM with 5 or 15% serum. We conclude that: (1) serum concentration in the growth medium and trypsin do appear to contribute to the differences in the radiosensitivity of spheroids and monolayer V79 cells; (2) these differences are associated with changes in cell morphology.     

Invited review Multicell spheroids as a model for cell kinetic studies

Abstract. Cells growing in tissue culture as three-dimensional, multicellular aggregates called 'spheroids' typically show a decreasing growth fraction and development of quiescent subpopulations as the spheroids enlarge. Kinetic studies in a number of spheroid systems have indicated that the primary reason for the tumour-like growth is a progressive decrease in growth fraction, with only a modest elongation of cell cycle time in larger spheroids. In this paper, the cellular growth kinetics for spheroids of V79 Chinese hamster lung cells are reviewed, and the regrowth kinetics of cells resuming growth after recovery from quiescent regions of the spheroids are described. Further, the role of regrowth/repopulation in determining the spheroid response to anti-tumour cytotoxics is explored, with particular emphasis on treatment with cisplatin and etoposide. By separating the effects of cytotoxicity and regrowth in the overall spheroid response to anti-neoplastic drugs, it is suggested that 'drug resistance' in tumours can be a kinetic as well as a genetic problem.     

3D Spheroids’ Sensitivity to Electric Field Pulses Depends on Their Size

Dramatic differences of cells behavior exist between cells cultured under classical 2D monolayers and 3D models, the latter being closer to in vivo responses. Thus, many 3D cell culture models have been developed. Among them, multicellular tumor spheroid appears as a nice and easy-to-handle 3D model based on cell adhesion properties. It is composed of one or several cell types and is widely used to address carcinogenesis, or drugs screening. A few and recent publications report the use of spheroids to investigate electropermeabilization process. We studied the response of spheroids to electrical field pulses (EP) in terms of their age, diameter or formation technique. We found that small human HCT-116 colorectal spheroids are more sensitive to electric field pulses than larger ones. Indeed, the growth of spheroids with a diameter of 300 μm decreased by a factor 2 over 4 days when submitted to EP (8 pulses, lasting 100 μs at a 1,300 V/cm field intensity). Under those electrical conditions, 650 μm spheroids were not affected. These data were the same whatever the formation method (i.e. hanging drop and nonadherent techniques). These observations point out the fact that characteristics of 3D cell models have to be taken into account to avoid biased conclusions of experimental data.     

Biochemical and functional changes of rat liver spheroids during spheroid formation and maintenance in culture: II. nitric oxide synthesis and related changes

Liver cells isolated from intact tissue can reaggregate to form three-dimensional, multicellular spheroids in vitro. During this process, cells undergo a histological and environmental change. How cells respond biochemically to this change has not been studied in detail previously. We have investigated some biochemical changes in rat liver cells during the formation and maintenance of spheroids. Liver cells were isolated from male Sprague rats and spheroids cultured by a gyrotatory-mediated method. Liver cells were shown to respond to the isolation procedure and the formation of spheroids triggered histological environmental changes that increased arginine uptake, nitric oxide (NO) and urea syntheses, as well as raised levels of GSH, GSSG, glutamic acid and aspartic acid secretion within the first couple of days after cell isolation. Levels were maintained at a relatively stable level in the mature spheroids (>5 days) over the 3 week period of observation. P450 1A1 activity was lost in the first 2 days and gradually recovered thereafter. This study, for the first time, shows that liver cells after isolation and during spheroid formation actively uptake arginine and increase NO and urea syntheses. A high level of NO is likely to play an important role in modulating a series of biochemical changes in liver cells. It is considered that liver cells actively respond to the 'challenge' induced by the isolation procedure and subsequent histological environmental changes, and biochemical modulation and instability result. The stable cell-cell contacts and histological environment in mature spheroids permit and support functional recovery and maintenance in vitro. This period of stability permits the use of spheroids in toxicity studies to establish acute and chronic paradigms.     

Modeling mechanical inhomogeneities in small populations of proliferating monolayers and spheroids

Understanding the mechanical behavior of multicellular monolayers and spheroids is fundamental to tissue culture, organism development, and the early stages of tumor growth. Proliferating cells in monolayers and spheroids experience mechanical forces as they grow and divide and local inhomogeneities in the mechanical microenvironment can cause individual cells within the multicellular system to grow and divide at different rates. This differential growth, combined with cell division and reorganization, leads to residual stress. Multiple different modeling approaches have been taken to understand and predict the residual stresses that arise in growing multicellular systems, particularly tumor spheroids. Here, we show that by using a mechanically robust agent-based model constructed with the peridynamic framework, we gain a better understanding of residual stresses in multicellular systems as they grow from a single cell. In particular, we focus on small populations of cells (1–100 s) where population behavior is highly stochastic and prior investigation has been limited. We compare the average strain energy density of cells in monolayers and spheroids using different growth and division rules and find that, on average, cells in spheroids have a higher strain energy density than cells in monolayers. We also find that cells in the interior of a growing spheroid are, on average, in compression. Finally, we demonstrate the importance of accounting for stochastic fluctuations in the mechanical environment, particularly when the cellular response to mechanical cues is nonlinear. The results presented here serve as a starting point for both further investigation with agent-based models, and for the incorporation of major findings from agent-based models into continuum scale models when explicit representation of individual cells is not computationally feasible.