
==== Front
Cancer Med
Cancer Med
10.1002/(ISSN)2045-7634
CAM4
Cancer Medicine
2045-7634
John Wiley and Sons Inc. Hoboken

10.1002/cam4.70239
CAM470239
CAM4-2024-05-3006.R1
Research Article
Research Article
Extracellular matrix production and oxygen diffusion regulate chemotherapeutic response in osteosarcoma spheroids
Sagheb et al.
Sagheb Isabel S. 1
Coonan Thomas P. 1
Randall R. Lor 2
Griffin Katherine H. https://orcid.org/0000-0003-1207-4624
2 3 khgriffin@ucdavis.edu

Leach J. Kent https://orcid.org/0000-0002-1673-3335
1 2 jkleach@ucdavis.edu

1 Department of Biomedical Engineering University of California Davis California USA
2 Department of Orthopaedic Surgery UC Davis Health Sacramento California USA
3 School of Veterinary Medicine University of California Davis California USA
* Correspondence
J. Kent Leach and Katherine H. Griffin, Department of Orthopaedic Surgery, UC Davis Health, 4860 Y Street, Suite 3800, Sacramento, CA 95817, USA.
Email: jkleach@ucdavis.edu and khgriffin@ucdavis.edu

20 9 2024
9 2024
13 18 10.1002/cam4.v13.18 e7023925 8 2024
22 5 2024
03 9 2024
© 2024 The Author(s). Cancer Medicine published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Osteosarcoma (OS) survival rates and outcome have not improved in 50 years since the advent of modern chemotherapeutics. Thus, there is a critical need for an improved understanding of the tumor microenvironment to identify better therapies. Extracellular matrix (ECM) deposition and hypoxia are known to abrogate the efficacy of various chemical and cell‐based therapeutics. Here, we aim to mechanistically investigate the combinatorial effects of hypoxia and matrix deposition with the use of OS spheroids.

Methods

We use two murine OS cell lines with differential metastatic potential to form spheroids. We form spheroids of two sizes, use ascorbate‐2‐phosphate supplementation to enhance ECM deposition, and study cell response under standard (21% O2) and physiologic (5% O2) oxygen tensions. Finally, we examine chemotherapeutic responses to doxorubicin treatment.

Results

ECM production and oxygen tension are key determinants of spheroid size through cell organization based on nutrient and oxygen distribution. Interestingly, highly metastatic OS is more susceptible to chemotherapeutics compared to less metastatic OS when matrix production increases. Together, these data suggest that dynamic interactions between ECM production and oxygen diffusion may result in distinct chemotherapeutic responses despite inherent tumor aggressiveness.

Conclusion

This work establishes OS spheroids as a valuable tool for early OS tumor formation investigation and holds potential for novel therapeutic target and prognostic indicator discovery.

I describe an in vitro investigation to study the influence of extracellular matrix (ECM) content and spheroid diameter (as a surrogate for oxygen tension) on the efficacy of chemotherapy in osteosarcoma (OS). Interestingly, highly metastatic OS was more susceptible to chemotherapeutics compared to less metastatic OS when matrix production increased. These data suggest that dynamic interactions between ECM production and oxygen diffusion may result in distinct chemotherapeutic responses despite inherent tumor aggressiveness.

chemoresistance
chemotherapy
hypoxia
osteosarcoma
tumorigenesis
Foundation for the National Institutes of Health 10.13039/100000009 R01 AR079211 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:20.09.2024
Sagheb IS , Coonan TP , Randall RL , Griffin KH , Leach JK . Extracellular matrix production and oxygen diffusion regulate chemotherapeutic response in osteosarcoma spheroids. Cancer Med. 2024;13 :e70239. doi:10.1002/cam4.70239

Isabel S. Sagheb and Thomas P. Coonan contributed equally to this study.
==== Body
pmc1 INTRODUCTION

Osteosarcoma (OS) is the most common bone cancer in children and adolescents, resulting in the formation of malignant osteoid. 1 , 2 These tumors are highly aggressive with very high metastatic potential. Doxorubicin‐based chemotherapy remains the backbone of systemic treatment with no improvement in survival over the past 50 years. 3 This disappointing absence of clinical progress mandates further mechanistic investigation to improve our understanding of primary tumor development and the events leading to metastasis.

The majority of existing research on OS is conducted under standard culture conditions (21% O2) and on tissue culture plastic (TCP), 4 , 5 which fails to mimic the native environment of OS. As a sarcoma, OS is characterized by bulk stromal tissue formation that arises in the bone marrow niche and features an oxygen tension of 5% O2. 6 Furthermore, hallmarks of aggressive OS involve selective matrix‐related changes that are clinically recognized as rapid tumor growth through matrix deposition that outpaces the nutrient and vascular supply, resulting in necrotic, hypoxic lesions. 1 This, in turn, can reduce drug diffusion and targeting towards those regions. 1 , 2 Indeed, this close relationship between endogenous, interstitial matrix deposition and oxygen tension is widely studied in many solid tumors including OS. Many molecular and genetic studies show that master regulatory signaling pathways such as NOTCH and WNT, 7 which are known under homeostasis to regulate not only matrix cues but also proliferation and differentiation, are also disrupted by hypoxic conditions, largely through the HIF signaling pathway. 8 Indeed, multiple components of these pathways, from direct ligand targeting to microRNA biogenesis, are under investigation as potential therapeutics. 9 , 10 However, the vast majority of these studies have not modeled the hypoxic, three‐dimensional (3D) tumor microenvironment adequately, resulting in inconsequential conclusions and no appreciable change in clinical outcomes.

In the last 20 years, techniques to better model cancer environments have become a focus of the research field. 5 Going beyond genetic interrogations, these techniques largely aim to characterize the role of various extracellular matrix (ECM) cues, such as chemical and mechanical signals, through 3D culture with engineered scaffolds. 5 , 11 Though limited applications exist for OS, multiple groups report significant changes to both primary tumor and immortalized OS cell line behavior under 3D conditions. 11 , 12 For example, OS cultured on hydroxyapatite scaffolds, 13 in an engineered bone marrow environment, 14 and collagen hydrogels containing cancellous bone, laminin, and fibronectin, 15 have all shown stark differences in tumorigenic behavior compared to cells cultured in monolayer on TCP. Cellular aggregates formed from human OS cell lines were used to investigate responses to various clinically relevant chemotherapeutic regimes, yet these studies did not interrogate responses in conjunction with oxygen tension and matrix production. 16 Mechanistic studies also show that stiffness‐driven changes that arise with ECM remodeling alter OS proliferation, differentiation, and immune regulation. 17 , 18 Thus, the study of OS in 3D is necessary to accurately evaluate OS tumor formation and development, yet there are few reports that seek to combine this approach while also interrogating oxygen‐related changes.

Here, we examine the cell‐ECM interactions of primary OS tumorigenesis under physiologically relevant oxygen tensions with the use of OS aggregates, known as spheroids. Using two murine OS cell lines with differential metastatic potential, we formed spheroids of two different sizes, studied them under standard (21% O2) and physiologic (5% O2) oxygen tensions, and examined chemotherapeutic responses to doxorubicin treatment. We also sought to understand the role of matrix deposition in OS tumorigenesis, and thus enriched endogenous matrix production in spheroids by media supplementation with ascorbate‐2‐phosphate (A2P). We hypothesize that spheroids will exhibit increases in pro‐tumorigenic markers, including proliferation and chemoresistance, as a function of metastatic potential and oxygen tension. Furthermore, we hypothesize that ECM deposition within these spheroids regulates these behaviors.

2 METHODS

2.1 Cell culture

Highly metastatic OS (K7M2, ATCC, Manassas, VA) and less metastatic OS (K12, courtesy of Dr. Kurt Weiss, University of Pittsburgh Medical Center) were expanded and maintained under standard (21% O2) culture conditions. Cells were cultured in DMEM (Invitrogen, Carlsbad, CA) supplemented with 10% fetal bovine serum (FBS) (Atlanta Biologics, Flowery Branch, GA) and 1% penicillin–streptomycin (P/S) (Gemini Bio Products, West Sacramento, CA).

2.2 Spheroid formation and culture

OS cells were formed into spheroids composed of either 5,000 (5 K) cells or 10,000 (10 K) cells. Spheroids were formed using a previously described forced aggregation method. 19 , 20 Briefly, OS cells (either 1.45×105 cells/mL for 5,000 cell spheroids or 2.9×105 cells/mL for 10,000 cell spheroids) were pipetted into 1.5% agarose molds, each containing 29 microwells, in 24‐well plates, and the plates were centrifuged at 500 × g for 8 min. Plates were maintained statically in standard or physiologic culture conditions for 72 h to form spheroids, after which media was refreshed immediately following this incubation and every 2–3 days for the duration of the study. Spheroid formation after 72 h is defined as the “Day 0” timepoint. Due to their small size, unless otherwise stated, biochemical quantifications are presented where each well of 29 spheroids serves as a pooled replicate. For endogenous enrichment of matrix deposition, media was supplemented with 50 μg/mL ascorbate‐2‐phosphate (A2P, Sigma‐Aldrich, St. Louis, MO) during and after formation. For chemotherapeutic response studies, media was supplemented with 0.1 μM doxorubicin (Neta Scientific, Hainesport, NJ).

2.3 Cell proliferation and viability characterization

Metabolic activity was quantified with an alamarBlue assay (ThermoFisher Scientific, Waltham, MA) prior to DNA quantification. Samples were incubated with the alamarBlue reagent for 1.5 h, analyzed, washed twice with PBS, then collected in 300 μL passive lysis buffer (Promega Madison, WI). DNA content was determined using the Quant‐iT PicoGreen DNA Assay Kit (Invitrogen, Carlsbad, CA). Cell viability was evaluated by a Live/Dead assay with calceinAM and propidium iodide staining (Invitrogen). Samples were imaged with a confocal microscope (Leica Stellaris 5, Leica Microsystems, Deerfield, IL) and representative images were assessed with ImageJ.

2.4 Histology

Spheroids were fixed in 10% buffered formalin, embedded in HistoGel (ThermoFisher Scientific), then paraffin embedded and sectioned at 5 μm thickness. Sections were stained with hematoxylin and eosin (H&E) and imaged with 10X and 20X objectives using a Nikon Eclipse TE2000U microscope. To visualize regions with oxygen tensions lower than 1.3% 21 , 22 within the spheroids, pimonidazole hydrochloride (200 μM; Hypoxyprobe; Chemicon, Temecula, CA) was added to the spheroids and incubated under standard or physiologic culture conditions for 2 h. 23 Spheroids were then fixed and processed as described above. Samples were rehydrated, exposed to heat‐mediated antigen retrieval with sodium citrate buffer for 20 min, and incubated in blocking buffer composed of 10% goat serum and 10 mg/mL Bovine Serum Albumin (BSA) at room temperature for 30 min. Primary detection for pimonidazole hydrochloride was performed with a 1:100 dilution of the conjugated rat monoclonal anti‐pimonidazole IgG included in the kit, and cells were counterstained with DAPI. Other slides were rehydrated and exposed to enzyme mediated antigen retrieval with Proteinase K at 37°C for 30 min. Samples were blocked as described before, after which they were incubated with recombinant anti‐integrin α5 antibody (ab112183, Abcam, Cambridge MA) at a concentration of 1:100 overnight at 4°C. Slides were then treated with a secondary goat anti‐rabbit antibody conjugated with AF488 (ab150081, Abcam) at a concentration of 1:200 for 1 h at room temperature, and again counterstained with DAPI.

2.5 Collagen quantification

Collagen content was evaluated with a hydroxyproline kit (Chrondrex, Woodinville, WA), where hydroxyproline comprises an estimated 14% of total collagen amino acid content. 24 Samples were collected in 5 M HCl, hydrolyzed at 120°C for 26 h, and prepared per manufacturer instructions.

2.6 Mechanical characterization

Spheroid compressive storage moduli were examined using a MicroTester G2 (CellScale, Waterloo, ON). 25 Individual spheroids of approximately 500–600 μm were loaded onto an anvil in a water bath filled with PBS. The spheroids were compressed by 25% of their diameter by a stainless‐steel platen attached to a tungsten microbeam over 30 s. Microbeams with 0.2032 diameter were used for all tests. Displacement and force were tracked via MicroTester software. The linear region of the compressive modulus versus nominal strain graph was recorded as the calculated modulus.

2.7 Statistics

Data are presented as mean ± standard deviation. Shapiro–Wilk normality tests and F‐tests to compare variances were performed, and statistical analysis was conducted with Prism 10.2.0 (GraphPad, San Diego, CA) software utilizing either unpaired t‐test assuming Gaussian distribution and one‐way or two‐way analysis of variance (ANOVA) with post hoc Tukey's test depending on the number of groups and comparisons. Groups with different letters indicate statistically significant differences (p < 0.05), while groups with the same letters are not significant; ns denotes no significance among groups.

3 RESULTS

3.1 Spheroids formed with K7M2 and K12 cells require extended formation time

Spheroid formation, regardless of technique, is concluded when self‐assembled cell aggregates have deposited adequate extracellular matrix (ECM) for mechanical stability without disintegration. 26 For most adherent cell types, this process is complete after 48 h. However, spheroids formed with K12 cells were neither stable nor spherical even at 72 h (Figure S1). Surprisingly, we noted significant cell migration up the microwell wall beginning as early as 24 h after centrifugation. To encourage increased cell–cell interactions, we attempted spheroid formation with lower cell densities, higher agarose percentages, and increased centrifugation speeds (data not shown), none of which altered this behavior. Instead, we found that artifact from the 3D printed negative molds 19 created “steps” within the microwells. These steps facilitated K12 migration and egress from the focal center where spheroids normally form. This was remedied with the use of molds that were formed by casting in agarose 20 instead of 3D printing, thus eliminating the printed artifacts. However, it is important to note that this original migratory behavior is quite unusual and was an early indication that cell signaling and matrix deposition in K12s may contribute to cancerous processes. This was also validated when, even with the revised molds, we found that OS spheroids still required an extended formation time of 72 h before becoming mechanically stable without disintegration.

3.2 OS spheroids have similar diameters regardless of initial cell seeding density

We evaluated OS cell behavior during spheroid formation. To investigate how cell density and oxygen tension affected early OS spheroid formation, we formed spheroids of two densities, 5,000 cells/spheroid (5 K) and 10,000 cells/spheroid (10 K) and under two oxygen tensions, 21% and 5% O2. We measured the diameter of all spheroids (Figure 1A,B) from brightfield images (Figure 1C). Unlike other cells of mesenchymal origin, 26 OS cells did not compact throughout the course of spheroid formation. Rather, spheroids formed with K7M2 cells generally did not change size over time, and spheroids formed with K12 cells slightly increased in diameter. Surprisingly, we found no significant differences between spheroids made with 5,000 and 10,000 cells, suggesting an autonomous rearrangement to optimize cell contacts and nutrient distribution. This is further supported by the significant difference we detected in spheroid diameter as a function of oxygen tension (Figure S2) where spheroids, regardless of cell type and density, show a 10% increase in size when formed under 21% O2 compared to 5% O2 (p = 0.0069). Together, these data suggest that spheroids formed under 21% O2 are larger due to increased oxygen availability and subsequent diffusion.

FIGURE 1 OS spheroids have similar diameters regardless of initial cell density. OS spheroid diameter during formation of spheroids formed with 5,000 (5 K) (A) or 10,000 (10 K) (B) cells per spheroid. (C) Representative brightfield microscopy images of spheroids during formation. Scale bar = 500 μm. Data are mean ± SD (n = 10–12 single spheroids, distributed over at least three separate wells). Groups with statistically significant differences based on two‐way ANOVA do not share the same letters; ns denotes no significance among groups.

3.3 Spheroids with 10,000 OS cells have decreased viability

Due to the observed similarities in spheroid size, we next sought to interrogate spheroid viability 1 day after formation. We evaluated DNA content as an indicator of cell number and overall metabolic activity with an alamarBlue assay. We also directly examined spheroid viability with a Live/Dead assay. After 3 days of formation, DNA content for spheroids formed with 5,000 K12 cells at 21% O2 was significantly increased compared to all other 5 K groups (Figure 2A). We noted a similar increase in spheroids formed with 10,000 K12 cells at 21% O2 (Figure 2C). This suggests that K12 proliferation is more dependent on oxygen availability compared to K7M2 cells, which were consistently similar in DNA content regardless of cell density and oxygen tension. We also observed that metabolic output for spheroids formed with 5,000 cells is minimal (Figure 2B), whereas for spheroids formed with 10,000 cells, metabolic activity is elevated with an increased trend for those formed under 5% O2 compared to 21% O2 (Figure 2D). This indicates that spheroids formed with 10,000 cells, especially under hypoxic conditions, upregulate cellular and metabolic function, likely to maintain overall cell survival. This, in turn, suggests that 10,000 cell spheroids may approach the limit of nutrient and oxygen availability, typically considered 100–200 μm in length. 27 , 28 These quantitative data are supported by qualitative Live/Dead images, which reveal similar live cell distribution through the spheroid except in spheroids formed with 10,000 K7M2 cells under 21% O2 (Figure 2E). Here, we observe pronounced dead cells towards the center of the spheroid, but this is not seen under hypoxic conditions. Together with diameter discrepancies noted between oxygen tensions in Figure 1, this finding implies that during spheroid formation, oxygen sensing and regulation determines optimal spheroid size for nutrient and oxygen diffusion. Perhaps mimicking in vivo tumor formation, the larger K7M2 cells at 21% O2 may form then rapidly outgrow their oxygen supply, resulting in the presence of dead cells.

FIGURE 2 Spheroids formed with 10,000 OS cells have decreased viability. DNA content and metabolic activity of OS spheroids formed with 5,000 (5 K) (A, B) or 10,000 (10 K) (C, D) cells per spheroid. (E) Representative Live/Dead confocal microscopy images of spheroids as max projected z‐stacks after 72 h of formation. Live cells are green, and dead cells are red. Scale bar = 500 μm. Data are mean ± SD (n = 3 spheroid wells, 29 spheroids/well). Groups with statistically significant differences based on two‐way ANOVA do not share the same letters; ns denotes no significance among groups.

3.4 Spheroids of 5,000 OS cells contain more collagen than those of 10,000 OS cells

Spheroids are comprised solely of cell mass and endogenous secreted ECM. 26 Given that spheroids exhibited similar diameters but differences in DNA content, we next aimed to determine if ECM deposition counterbalanced cellular material to yield the observed similarities in spheroid size. To investigate this, we performed H&E staining for general morphology to evaluate cell and matrix distribution. We then directly interrogated ECM content with a hydroxyproline assay to measure collagen, the primary protein component of the ECM. 7 , 29 Spheroids stained with H&E exhibited uniform cell and matrix distribution (Figure 3A). Collagen content was often decreased in spheroids formed with 5,000 OS cells (Figure 3B,C), and the trends were roughly inversely proportional to DNA content presented in Figure 2. This confirms that the combined cellular and matrix content within OS spheroids of different cell numbers accounts for similarities in spheroid diameter.

FIGURE 3 Spheroids formed of 5,000 OS cells contain more collagen than those of 10,000 cells. (A) Representative H&E‐stained spheroids 1 day after formation. Scale bar = 100 μm. Collagen content as measured by hydroxyproline assay for spheroids formed with K7M2 (B) and K12 (C) cells. Data are mean ± SD (n = 3 spheroid wells, 29 spheroids/well). Groups with statistically significant differences based on two‐way ANOVA do not share the same letters; ns denotes no significance among groups.

3.5 Spheroids of 10,000 K12 cells have increased storage modulus compared to those with 5,000 cells

To further understand how ECM content may affect OS spheroid behavior, we examined individual spheroid compressive storage modulus as a functional output for ECM content. Using a MicroTester (Figure 4A), we loaded single spheroids onto an anvil and compressed them with a platen attached to a beam, which recorded force and displacement (Figure 4B, Video S1). Data were analyzed to yield compressive storage modulus. Compressive modulus exhibited a significant increase in 10 K compared to 5 K spheroids only for K12 spheroids (Figure 4E,F). We did not observe a similar trend in K7M2 spheroids formed under either oxygen tension (Figure 4C,D). Since collagen content in spheroids formed with 10,000 cells was consistently reduced compared to 5 K spheroids, these reported storage moduli were surprising, as collagen is largely thought to provide the structural integrity of most tissues. 29 , 30

FIGURE 4 Spheroids of 10,000 K12 cells have increased compressive storage modulus. A MicroTester (A, B) was used to measure the compressive storage modulus of spheroids formed under standard (21% O2) (C, E) and physiologic (5% O2) (D, F) culture conditions. Data are mean ± SD (n = 3–6 single spheroids, distributed over at least three separate wells). Groups with statistically significant differences based on unpaired t‐test do not share the same letters; ns denotes no significance among groups.

3.6 Spheroids formed with 10,000 K7M2 cells in standard conditions have increased integrin α5 expression

To elucidate the direct interactions between OS cells and ECM, we next investigated the role of integrins, which are cell‐ECM binding ligands. Integrins are known to play major roles throughout metastasis, from primary tumor cell egress, attachment to the vascular endothelium, and further migration to metastatic niches. 31 Capitalizing on the differential metastatic potentials of our cell lines, we chose an integrin subunit that could help delineate a mechanism for these differences and thus OS aggressiveness. Integrin subunit alpha5 (integrin α5) has been heavily investigated as a regulator of stem cell differentiation and migration, 32 bone formation, 33 and breast cancer in the context of osseus metastasis. 34 Though significantly less explored in OS, some studies show that inhibition of integrin α5 with other subunits, such as β1, results in decreased OS migration and proliferation. 35 Thus, we performed immunohistochemical staining on OS spheroids 1 day after formation for integrin α5 (Figure 5). We found that spheroids formed with 10 K K7M2 cells under standard (21% O2) culture conditions had a striking and distinct upregulation of integrin α5 near, but not on, the periphery of the spheroids. Interestingly, this localization corresponds with the Live/Dead staining pattern shown in Figure 2. This again harkens the idea that this particular group of spheroids may mimic in vivo tumor formation, with possible rapid migration away from a hypoxic and nutrient‐depleted core.

FIGURE 5 Spheroids formed with 10,000 K7M2 cells cultured under standard conditions have increased integrin α5 expression. Representative confocal microscopy images of sectioned spheroids with immunohistochemical staining for integrin α5 (green) counterstained with DAPI (blue). Scale bar = 200 μm.

We tested spheroids for regions with less than 1.3% O2 through pimonidazole treatment and staining to identify if there existed an anoxic region that would prompt the migration pattern we observed. Based on representative confocal microscopy images, we did not observe any areas of positive pimonidazole staining (Figure S3). Thus, none of these spheroids exhibited hypoxic regions with less than 1.3% O2, suggesting that an oxygen tension threshold between 1.3%–5% induces OS apoptosis, tumor necrosis, and cell migration away from the region.

3.7 Ascorbic acid supplementation does not change OS proliferation or spheroid compressive storage modulus

To interrogate the interplay between matrix deposition and OS behavior, we amplified ECM deposition by supplementing OS spheroids with 50 μg/mL A2P during formation. To characterize OS response, we again evaluated DNA content as an indicator of cell number with PicoGreen (Figure 6A,B), compressive storage modulus as a functional measure of ECM production with the MicroTester (Figure 6C,D), and collagen content with a hydroxyproline assay (Figure 6E,F). We also compared these A2P‐supplemented spheroids with non‐supplemented spheroid data previously shown. Broadly, we detected no statistically significant differences in spheroids formed with A2P across all outputs. However, as validation of the effects A2P supplementation, we did see consistent increasing trends in collagen content in spheroids formed with A2P across all conditions. Interestingly, Live/Dead staining showed increased dead cells across all groups (Figure 6G) compared to non‐A2P supplemented spheroids (Figure 2E), yet live and dead cells are evenly distributed within the spheroid. This perhaps suggests that increased ECM production interferes directly with oxygen and nutrient diffusion or indirectly through interruption of OS cell ability to sense oxygen and nutrient availability. Taken together, these data suggest that ECM production alone does not dictate OS proliferation and bulk tumor mechanical properties but may instead influence molecular signaling processes and subsequent cellular organization.

FIGURE 6 A2P does not change OS proliferation or spheroid compressive storage modulus. Spheroids formed with K7M2 and K12 cells were supplemented with 50 μg/mL A2P during formation to promote ECM deposition. One day after formation, OS spheroid DNA content (A, B), compressive storage modulus (C, D), and collagen content (E, F) were measured. (G) Live/Dead confocal microscopy images, as max projected z‐stacks, of A2P‐supplemented spheroids after 1 day after formation. Live cells are green, dead cells are red. Scale bar = 500 μm. Data are mean ± SD (n = 3, spheroid wells, 29 spheroids/well for A, B, E, F; n = 3–7, single spheroids, distributed over at least 3 separate wells for C, D). Groups with statistically significant differences based on two‐way ANOVA do not share the same letters; ns denotes no significance among groups.

3.8 Increased ECM content increases K7M2 susceptibility to doxorubicin but decreases K12 susceptibility

To enhance the clinical relevance of our findings, we tested OS spheroid response to doxorubicin (DOX). OS spheroids were supplemented with 0 or 50 μg/mL A2P and treated with 0.1 μM DOX for 7 days after formation. DNA content was measured as an indicator of cell number. Data were normalized to DNA content of spheroid groups 1 day after formation to yield normalized proliferation. DOX was not cytotoxic to OS spheroids without A2P, regardless of cell type and culture condition, suggesting that 3D culture as spheroids is inherently chemoprotective. Surprisingly, when spheroids formed with K7M2 cells and A2P were treated with DOX, there were stark decreases in proliferation (Figure 7A,B). However, K12 cells exhibited an entirely different behavior. When supplemented with A2P and treated with DOX, spheroids formed with K12 cells exhibited either the same (Figure 7C) or significantly increased proliferation (Figure 7D). This reveals that K12 cells with excess matrix production have enhanced chemoresistance, especially compared to K7M2s.

FIGURE 7 Increased ECM content influences OS spheroid susceptibility to doxorubicin as a function of metastatic potential. OS spheroids without A2P or supplemented with A2P were treated with 0.1 μM DOX for 7 days after formation. DNA content was measured for spheroids formed with K7M2 (A, B) and K12 (C, D) cells. Data are presented as Day 7 DNA content normalized to Day 1 DNA content. Data are mean ± SD (n = 3 spheroid wells, 29 spheroids/well). Groups with statistically significant differences based on two‐way ANOVA do not share the same letters; ns denotes no significance among groups.

4 DISCUSSION

It is broadly recognized that cellular processes involved in dedifferentiation, oxygen tension, and ECM deposition and remodeling are major drivers in the development of metastasis. 6 , 36 However, little is known about the mechanisms that determine these changes 37 and if early prognostic indicators specifically related to cellular functions, rather than genetic or familial changes which are historically unreliable for sarcomas, 2 , 38 , 39 may exist. Furthermore, it remains unclear whether aggressive tumors exhibit extensive ECM remodeling or if ECM changes drive tumor aggression. In this work, we studied the interactions between OS cells and ECM in early tumor formation with the use of an OS spheroid model. Specifically, we formed spheroids with highly metastatic murine K7M2 cells and less metastatic K12 cells under standard (21% O2) and physiologic (5% O2) conditions to delineate differences in early tumor formation and matrix changes as a function of OS cell aggressiveness and oxygen tension.

Our data confirm that spheroid size was largely determined by oxygen tension. Clinically, primary tumor necrosis is a valuable prognostic indicator. 40 Tumors that have outgrown their vascular and nutrient supply are characterized by low oxygen levels. 1 , 40 Indeed, we noted that large K7M2 spheroids formed under 21% O2 may show early signs of a necrotic core and cell migration away from that region. Though the clinical literature has no reported oxygen levels within necrotic OS tumors, multiple studies report the upregulation of hypoxia‐associated gene pathways, namely HIF‐1α. 41 , 42 Interestingly, HIF‐1α is most often upregulated under short periods of near anoxic levels (<0.1% O2), 41 but our data suggest that aggressive OS tolerance for hypoxia may decrease at levels above 1.3% O2. It is possible that the induced migration observed here is due to a well‐documented increase in K7M2 ezrin protein expression, 38 which has known correlations with hypoxia‐induced cell motility and tumor invasion in other cancers. 43 , 44

OS spheroids maximized matrix deposition and cellular content in relation to oxygen and nutrient availability. We evaluated matrix deposition primarily through collagen characterization, yet this oversimplifies in vivo tumor ECM. 36 Though collagen is the most abundant component of ECM, fibronectin plays an important role in the OS matrix and upregulation is associated with chemoresistance. 45 Indeed integrin α5 is a known ligand for fibronectin, 32 and studies with human MG‐63 OS spheroids show increased integrin α5 protein expression with hypoxic treatments. 46 Furthermore, a hallmark of osteoblastic OS is excessive, immature osteoid deposition with increased hydroxyapatite (HA). 15 This provides increased substrate stiffness from mineral content and additional cell adhesion and calcium nucleation sites. 47 Though unexplored here, other studies have shown that HA incorporation in human OS spheroids increases stemness markers and HIF‐1α. 12 Unfortunately, these cell lines are not osteogenic despite their osteoblastic origins, representing a limitation of this work.

The combined effects of hypoxia and matrix deposition reported here are also important to consider. As a sarcoma, clinical ECM deposition is so abundant that it effectively blockades the tumor and inhibits immune cell infiltration. 6 , 7 This process creates pockets of lower, but not anoxic, oxygen tension within tumors, where numerous hypoxia‐related pathways in addition to HIF‐1α become upregulated. 12 Multiple studies report that associated changes drive a more aggressive, though less differentiated phenotype. 48 Within cancer, tumor cell dedifferentiation is largely referred to as the epithelial to mesenchymal transition, where an upregulation in stem cell‐associated genes, such as NANOG, OCT‐4, and SOX‐2, induce a phenotype that is broadly known as a cancer stem cell (CSC). 6 , 48 CSCs are known to be more chemoresistant, yet often enter a quiescent state in hypoxic environments. 6 , 48 Specifically in OS, studies show that NOTCH signaling and associated increases in aldehyde dehydrogenase (ALDH) expression are reliable markers of stemness. 49 K7M2s, under standard culture conditions, have increased NOTCH signaling, yet hypoxia is known to upregulate NOTCH ligand expression. 49 As K12s show increased differentiation compared to K7M2s under baseline conditions, 38 it is possible that 5% O2 culture conditions may effectively increase K12 stemness. This would correspond with our observations that K12 spheroids supplemented with A2P at 5% O2 appear more chemoresistant, even compared to K7M2s under the same conditions due to relative changes in NOTCH expression. However, another possible mechanism to explain this finding is that ECM composition may also differ between cell lines, where K7M2s may produce more fibronectin. As discussed already, this is indirectly supported by the integrin α5 staining. These differences may promote subtle changes in diffusion that could impact chemotherapeutic resistance, and thus warrants further investigation with mass spectrometry. As such, this model could serve as an effective platform to investigate CSC‐related changes as a function of oxygen tension and differentiation.

In conjunction with future studies that further evaluate ECM changes, OS‐mediated matrix degradation and remodeling, and associated genetic changes, our findings provide a foundation for OS ECM research on prognostic indicators. However, there are notable limitations beyond those already mentioned. As cell lines, K7M2 and K12 cells have been selected for culture under 21% O2 during isolation. 50 This renders hypoxic studies less accurate. Additionally, the use of homotypic spheroids fails to capture important crosstalk between cell types. Future directions for this work may include co‐culture with fibroblasts to better understand the interplay between OS signaling and fibroblast‐driven changes to ECM production.

Overall, this study shows that ECM production and oxygen tension regulate OS behavior. Specifically, these two factors are influential determinants of spheroid size through cell organization based on nutrient and oxygen distribution, and we also show that these responses are generally not affected by metastatic potential. The exceptions to this are (1) spheroids of 10,000 K7M2 cells formed in 21% O2 feature unique viability localizations and cell migration patterns and (2) spheroids formed with 10,000 K12 cells and A2P under 5% O2, which, when treated with DOX, show increased chemoresistance, may be due to a hypoxia‐induced increase in stemness. These studies are then contextualized as they relate to key clinical markers, such as HIF‐1α, tumor ECM and mineral content, cancer stem cells, and chemo‐responsiveness. This work presents OS spheroids as an excellent platform for further research to mechanistically study ECM‐OS interactions. Further, after extensive validation, spheroids formed with primary tumor isolates have the potential to provide much needed, accurate personalized drug screening, as well as other mechanistic investigations for the discovery of novel clinical therapeutic targets and prognostic indicators.

AUTHOR CONTRIBUTIONS

Isabel S. Sagheb: Conceptualization (supporting); data curation (equal); investigation (equal); writing – review and editing (supporting). Thomas P. Coonan: Conceptualization (supporting); data curation (equal); investigation (equal); methodology (supporting). R. Lor Randall: Conceptualization (supporting); funding acquisition (lead); writing – review and editing (equal). Katherine H. Griffin: Conceptualization (equal); data curation (lead); formal analysis (equal); investigation (equal); methodology (equal); writing – original draft (lead); writing – review and editing (equal). J. Kent Leach: Conceptualization (equal); funding acquisition (lead); investigation (supporting); methodology (supporting); project administration (lead); resources (lead); supervision (lead); writing – review and editing (lead).

CONFLICT OF INTEREST STATEMENT

The authors declare no conflicts of interest.

Supporting information

Video S1. Representative video of mechanical testing of an osteosarcoma spheroid (K12 spheroid formed of 10,000 cells at 21% O2) on the MicroTester.

Data S1.

ACKNOWLEDGMENTS

This work was supported in part by the UC Davis School of Veterinary Medicine Endowment Funds and Graduate Student Support Program to KHG and the Shingle Springs Subaru “Share the Love” Campaign to RLR. This work was also partially supported by the National Institutes of Health (grant number R01 AR079211). JKL acknowledges financial support from the Lawrence J. Ellison Endowed Chair of Musculoskeletal Research.

DATA AVAILABILITY STATEMENT

The data that support the findings of this study are available from the corresponding author upon reasonable request.
==== Refs
REFERENCES

1 Geller DS , Gorlick R . Osteosarcoma: a review of diagnosis, management, and treatment strategies. Clin Adv Hematol Oncol. 2010;8 :705‐718.21317869
2 Eaton BR , Schwarz R , Vatner R , et al. Osteosarcoma. Pediatr Blood Cancer. 2021;68 (Suppl 2 ):e28352.32779875
3 Misaghi A , Goldin A , Awad M , Kulidjian AA . Osteosarcoma: a comprehensive review. SICOT J. 2018;4 :12.29629690
4 Peterse EFP , van Leeuwen TN , Cleton‐Jansen AM . In vitro studies of osteosarcoma: a researcher's perspective of quantity and quality. J Bone Oncol. 2017;7 :29‐31.28443232
5 Thai VL , Griffin KH , Thorpe SW , Randall RL , Leach JK . Tissue engineered platforms for studying primary and metastatic neoplasm behavior in bone. J Biomech. 2021;115 :110189.33385867
6 Yang C , Tian Y , Zhao F , et al. Bone microenvironment and osteosarcoma metastasis. Int J Mol Sci. 2020;21 (19 ):6985.32977425
7 Winkler J , Abisoye‐Ogunniyan A , Metcalf KJ , Werb Z . Concepts of extracellular matrix remodelling in tumour progression and metastasis. Nat Commun. 2020;11 :5120.33037194
8 Sahlgren C , Gustafsson MV , Jin S , Poellinger L , Lendahl U . Notch signaling mediates hypoxia‐induced tumor cell migration and invasion. Proc Natl Acad Sci USA. 2008;105 :6392‐6397.18427106
9 Yu F , Yu C , Li F , et al. Wnt/beta‐catenin signaling in cancers and targeted therapies. Signal Transduct Target Ther. 2021;6 :307.34456337
10 Rupaimoole R , Slack FJ . MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;16 :203‐222.28209991
11 De Luca A , Raimondi L , Salamanna F , et al. Relevance of 3d culture systems to study osteosarcoma environment. J Exp Clin Cancer Res. 2018;37 :2.29304852
12 Bassi G , Panseri S , Dozio SM , et al. Scaffold‐based 3D cellular models mimicking the heterogeneity of osteosarcoma stem cell niche. Sci Rep. 2020;10 :22294.33339857
13 Gonzalez Diaz EC , Lee AG , Sayles LC , et al. A 3D osteosarcoma model with bone‐mimicking cues reveals a critical role of bone mineral and informs drug discovery. Adv Healthc Mater. 2022;11 :e2200768.35767377
14 Griffin KH , Thorpe SW , Sebastian A , et al. Engineered bone marrow as a clinically relevant ex vivo model for primary bone cancer research and drug screening. Proc Natl Acad Sci USA. 2023;120 :e2302101120.37729195
15 Pavlou M , Shah M , Gikas P , Briggs T , Roberts SJ , Cheema U . Osteomimetic matrix components alter cell migration and drug response in a 3D tumour‐engineered osteosarcoma model. Acta Biomater. 2019;96 :247‐257.31302294
16 Collier CD , Wirtz EC , Knafler GJ , Morris WZ , Getty PJ , Greenfield EM . Micrometastatic drug screening platform shows heterogeneous response to MAP chemotherapy in osteosarcoma cell lines. Clin Orthop Relat Res. 2018;476 :1400‐1411.29481344
17 Chim LK , Williams IL , Bashor CJ , Mikos AG . Tumor‐associated macrophages induce inflammation and drug resistance in a mechanically tunable engineered model of osteosarcoma. Biomaterials. 2023;296 :122076.36931102
18 Luu AK , Viloria‐Petit AM . Targeting mechanotransduction in osteosarcoma: a comparative oncology perspective. Int J Mol Sci. 2020;21 (20 ):7595.33066583
19 Gonzalez‐Fernandez T , Tenorio AJ , Leach JK . Three‐dimensional printed stamps for the fabrication of patterned microwells and high‐throughput production of homogeneous cell spheroids. 3D Print Addit Manuf. 2020;7 :139‐147.32855996
20 Vorwald CE , Ho SS , Whitehead J , Leach JK . High‐throughput formation of mesenchymal stem cell spheroids and entrapment in alginate hydrogels. Methods Mol Biol. 2018;1758 :139‐149.29679328
21 Chapman JD , Franko AJ , Sharplin J . A marker for hypoxic cells in tumours with potential clinical applicability. Br J Cancer. 1981;43 :546‐550.7236492
22 Raleigh JA , Koch CJ . Importance of thiols in the reductive binding of 2‐nitroimidazoles to macromolecules. Biochem Pharmacol. 1990;40 :2457‐2464.2176499
23 Murphy KC , Hung BP , Browne‐Bourne S , et al. Measurement of oxygen tension within mesenchymal stem cell spheroids. J R Soc Interface. 2017;14 :20160851.28179546
24 Jimenez S , Harsch M , Rosenbloom J . Hydroxyproline stabilizes the triple helix of chick tendon collagen. Biochem Biophys Res Commun. 1973;52 :106‐114.4576266
25 Gonzalez‐Fernandez T , Tenorio AJ , Saiz AM Jr , et al. Engineered cell‐secreted extracellular matrix modulates cell spheroid mechanosensing and amplifies their response to inductive cues for the formation of mineralized tissues. Adv Healthc Mater. 2022;11 :e2102337.34968011
26 Gionet‐Gonzales MA , Leach JK . Engineering principles for guiding spheroid function in the regeneration of bone, cartilage, and skin. Biomed Mater. 2018;13 :034109.29460842
27 Wagner BA , Venkataraman S , Buettner GR . The rate of oxygen utilization by cells. Free Radic Biol Med. 2011;51 :700‐712.21664270
28 Grimes DR , Kelly C , Bloch K , Partridge M . A method for estimating the oxygen consumption rate in multicellular tumour spheroids. J R Soc Interface. 2014;11 :20131124.24430128
29 Frantz C , Stewart KM , Weaver VM . The extracellular matrix at a glance. J Cell Sci. 2010;123 :4195‐4200.21123617
30 Kular JK , Basu S , Sharma RI . The extracellular matrix: structure, composition, age‐related differences, tools for analysis and applications for tissue engineering. J Tissue Eng. 2014;5 :2041731414557112.25610589
31 Hamidi H , Ivaska J . Every step of the way: integrins in cancer progression and metastasis. Nat Rev Cancer. 2018;18 :533‐548.30002479
32 Neiman G , Scarafia MA , La Greca A , et al. Integrin alpha‐5 subunit is critical for the early stages of human pluripotent stem cell cardiac differentiation. Sci Rep. 2019;9 :18077.31792288
33 Mao L , Wang L , Xu J , Zou J . The role of integrin family in bone metabolism and tumor bone metastasis. Cell Death Dis. 2023;9 :119.
34 Pantano F , Croset M , Driouch K , et al. Integrin alpha5 in human breast cancer is a mediator of bone metastasis and a therapeutic target for the treatment of osteolytic lesions. Oncogene. 2021;40 :1284‐1299.33420367
35 Kawashima A , Kawahara E , Tokuda R , Nakanishi I . Tumour necrosis factor‐alpha provokes upregulation of alpha2beta1 and alpha5beta1 integrins, and cell migration in OST osteosarcoma cells. Cell Biol Int. 2001;25 :319‐329.11319838
36 Cui J , Dean D , Hornicek FJ , Chen Z , Duan Z . The role of extracelluar matrix in osteosarcoma progression and metastasis. J Exp Clin Cancer Res. 2020;39 :178.32887645
37 Griffin KH , Mizenko RR , Arun V , Carney RP , Leach JK . Extracellular vesicles from highly metastatic osteosarcoma cells induce pro‐tumorigenic macrophage phenotypes. Adv Biol. 2024;8 :e2300577.
38 Khanna C , Khan J , Nguyen P , et al. Metastasis‐associated differences in gene expression in a murine model of osteosarcoma. Cancer Res. 2001;61 :3750‐3759.11325848
39 Calvert GT , Randall RL , Jones KB , Cannon‐Albright L , Lessnick S , Schiffman JD . At‐risk populations for osteosarcoma: the syndromes and beyond. Sarcoma. 2012;2012 :152382.22550413
40 Beird HC , Bielack SS , Flanagan AM , et al. Osteosarcoma. Nat Rev Dis Primers. 2022;8 :77.36481668
41 Kumar H , Choi DK . Hypoxia inducible factor pathway and physiological adaptation: a cell survival pathway? Mediat Inflamm. 2015;2015 :584758.
42 Pierrevelcin M , Fuchs Q , Lhermitte B , et al. Focus on hypoxia‐related pathways in pediatric osteosarcomas and their druggability. Cells. 2020;9 :1998.32878021
43 Qureshi‐Baig K , Kuhn D , Viry E , et al. Hypoxia‐induced autophagy drives colorectal cancer initiation and progression by activating the PRKC/PKC‐EZR (ezrin) pathway. Autophagy. 2020;16 :1436‐1452.31775562
44 Maccio U , Mihic A , Lenggenhager D , et al. Hypoxia and ezrin expression in primary melanoma have high prognostic relevance. Int J Mol Sci. 2022;23 :10745.36142656
45 Kun‐Peng Z , Chun‐Lin Z , Xiao‐Long M , Lei Z . Fibronectin‐1 modulated by the long noncoding RNA OIP5‐AS1/miR‐200b‐3p axis contributes to doxorubicin resistance of osteosarcoma cells. J Cell Physiol. 2019;234 :6927‐6939.30204936
46 Indovina P , Rainaldi G , Santini MT . Hypoxia increases adhesion and spreading of MG‐63 three‐dimensional tumor spheroids. Anticancer Res. 2008;28 :1013‐1022.18507049
47 Qing F , Wang Z , Hong Y , et al. Selective effects of hydroxyapatite nanoparticles on osteosarcoma cells and osteoblasts. J Mater Sci Mater Med. 2012;23 :2245‐2251.22903597
48 Jubelin C , Munoz‐Garcia J , Cochonneau D , et al. Biological evidence of cancer stem‐like cells and recurrent disease in osteosarcoma. Cancer Drug Resist. 2022;5 :184‐198.35582537
49 Mu X , Isaac C , Greco N , et al. Notch signaling is associated with ALDH activity and an aggressive metastatic phenotype in murine osteosarcoma cells. Front Oncol. 2013;3 :143.23805413
50 Khanna C , Prehn J , Yeung C , Caylor J , Tsokos M , Helman L . An orthotopic model of murine osteosarcoma with clonally related variants differing in pulmonary metastatic potential. Clin Exp Metastasis. 2000;18 :261‐271.11315100
