Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 12 December 2023
Sec. Tissue Engineering and Regenerative Medicine

Spheroid size influences cellular senescence and angiogenic potential of mesenchymal stromal cell-derived soluble factors and extracellular vesicles

  • 1Department of Experimental Medicine, University of Genoa, Genoa, Italy
  • 2IRCCS Ospedale Policlinico San Martino, Genoa, Italy

Introduction: The secretome of mesenchymal stromal cells (MSCs) serves as an innovative tool employed in the regenerative medicine approach. In this particular context, three-dimensional (3D) culture systems are widely utilized to better replicate in vivo conditions and facilitate prolonged cell maintenance during culture. The use of spheroids enables the preservation of the classical phenotypical characteristics of MSCs. However, the distinct microenvironment within the spheroid may impact the secretome, thereby enhancing the angiogenic properties of adult MSCs that typically possess a reduced angiogenic potential compared to MSCs derived from perinatal tissues due to the hypoxia created in the internal region of the spheroid.

Methods: In this study, large spheroids (2,600 cells, ∼300 μm diameter) and small spheroids (1,000 cells, ∼200 μm diameter) were used to examine the role of spheroid diameter in the generation of nutrients and oxygen gradients, cellular senescence, and the angiogenic potential of secreted factors and extracellular vesicles (EVs).

Results: In this study, we demonstrate that large spheroids showed increased senescence and a secretome enriched in pro-angiogenic factors, as well as pro-inflammatory and anti-angiogenic cytokines, while small spheroids exhibited decreased senescence and a secretome enriched in pro-angiogenic molecules. We also demonstrated that 3D culture led to a higher secretion of EVs with classical phenotypic characteristics. Soluble factors and EVs from small spheroids exhibited higher angiogenic potential in a human umbilical vein endothelial cell (HUVEC) angiogenic assay.

Discussion: These findings highlighted the necessity of choosing the appropriate culture system for obtaining soluble factors and EVs for specific therapeutic applications.

1 Introduction

The secretion of paracrine factors and extracellular vesicles (EVs) by mesenchymal stromal cells (MSCs) has gained a lot of attention in the field of regenerative medicine (Rani et al., 2015; Ferreira et al., 2018; Praveen Kumar et al., 2019; Rezaie et al., 2019). MSCs are largely known for their immunomodulatory, anti-inflammatory, and anti-microbial effects (Kyurkchiev et al., 2014; Fierabracci et al., 2016; Marrazzo et al., 2019; Markov et al., 2021). However, their ability to stimulate angiogenesis has also been widely reported (Watt et al., 2013). MSCs indeed produce and secrete a plethora of pro-angiogenic factors, such as angiogenin, vascular endothelial growth factor-A (VEGF-A), interleukin-8 (IL-8), and hepatocyte growth factor (HGF), able to promote new vessel formation (Burlacu et al., 2013; Liu et al., 2015). Interestingly, MSCs isolated from different sources possess different angiogenic potentials. More deeply, MSCs obtained from the amniotic fluid and umbilical cord showed better support for the angiogenic process than bone marrow-derived MSCs, suggesting that MSCs from early developmental stages are more supportive than adult MSCs (Roura et al., 2012; Roubelakis et al., 2013; Pinto et al., 2020).

The limited angiogenic potential surely presents a limitation for the use of adult MSCs, but the most relevant constraint is related to the onset of cellular senescence during the long-term expansion required to obtain a clinically relevant cell amount. Cellular senescence is characterized by cell proliferation blockage, altered cell morphology (Neurohr et al., 2019), increased DNA damage, and a change in gene expression and a secretory profile known as the senescence-associated secretory phenotype (SASP). This is characterized by the high secretion of pro-inflammatory cytokines, reduced secretion of cytokines with immunomodulatory and anti-inflammatory effects, and reduced production of EV (Coppé et al., 2010; Turinetto et al., 2016; Jeske et al., 2021). Specifically, SASP is characterized by the upregulation of different pro-inflammatory cytokines, metalloproteases, and growth factors (Acosta et al., 2013; Lunyak et al., 2017). The activation of cellular senescence in in vitro cultures is usually related to the onset of replicative senescence, which can be caused by a plethora of factors, such as telomere shortening, DNA damage, and oxidative stress (Wagner et al., 2008; Pizzuti et al., 2023). Analysis of the cellular senescence of MSCs in vitro is extremely important for controlling the quality of cells used in cellular therapy approaches. It has been widely observed that the onset of replicative senescence may lead to a complete alteration of MSC characteristics and a reduction of their regenerative potential. In addition, MSC senescence has been proved to impair the immunomodulatory and anti-inflammatory abilities of MSCs (Meng et al., 2017).

Three-dimensional (3D) culture methods represent an alternative strategy that could potentially overcome the major limitations of classical two-dimensional (2D) cultures. These systems could better recapitulate the in vivo conditions, allowing them to recreate a more physiological environment with proper cell–cell and cell–matrix interactions (Soares et al., 2012; Costa et al., 2016; Langhans, 2018). Among 3D culture methods, spheroids are most widely used in cancer research, but they also represent a promising approach in other research fields, such as stem cells and regenerative medicine. Spheroid cultures of MSCs have indeed already been reported to allow more efficient cell organization and enhance the therapeutic potential of MSCs (Chen et al., 2021; Kouroupis and Correa, 2021). For instance, it has been demonstrated that MSCs cultured as spheroids possess enhanced immunomodulatory characteristics as well as increased anti-inflammatory effects (Regmi et al., 2021; Xie et al., 2021). 3D architecture also supports a more in vivo-like cell behavior, leading to a reorganization of the cell junction (i.e., with the expression of E-cadherin) and cytoskeleton and a drastic change in cell morphology (Cesarz and Tamama, 2016). Furthermore, some authors demonstrated that MSC spheroid cultures can stimulate an increase in cell viability (Regmi et al., 2021), differentiation capability, and regenerative potential (Ohori-Morita et al., 2022), sustained by a delay of replicative senescence (Cesarz and Tamama, 2016).

However, spheroid culture also presents some limitations due to the generation of different regions inside the aggregates. More specifically, the outer region, being exposed to external nutrients, is mostly populated by proliferating cells, the intermediate region hosts senescent cells, and the internal region becomes a necrotic and apoptotic core of cells due to the establishment of nutrients and oxygen gradients, causing cell suffering (Mueller-Klieser, 1984; Alvarez-Pérez et al., 2005; Curcio et al., 2007; Kouroupis and Correa, 2021). Several studies demonstrated that the control of spheroid diameter is crucial in the generation of these gradients, and a diameter of 200 μm was set as the threshold for the generation of the internal hypoxic core (Alvarez-Pérez et al., 2005; Curcio et al., 2007; Bhang et al., 2011). However, the formation of mild nutrients and oxygen gradients always occurs even in small-diameter spheroids. This process could work as internal priming for the cells, leading to an enhancement of the angiogenic properties of MSC-derived soluble factors and EVs (Cheng et al., 2013; Costa et al., 2017; Redondo-Castro et al., 2018; Ezquerra et al., 2021). The standardized production of spheroids with similar dimensions and characteristics is of extreme importance for their use in regenerative medicine. Moreover, the possible use of label-free methods for characterization, such as mass density analysis, could be an extremely useful system to determine the density and mass of the spheroid since alteration may suggest different secretory characteristics (Marrazzo et al., 2021).

The present study aims to demonstrate the hypothesis that spheroid diameter is crucial for the regulation of cell senescence and secretome composition, as well as for EV secretion. For this purpose, bone marrow-derived human MSCs were used to produce large- and small-sized spheroids. Cell properties, such as stemness and senescence, were analyzed to evaluate whether spheroid size might influence cell behavior. Moreover, the pro-angiogenic activity of secretome and EVs was evaluated on human umbilical vein endothelial cells (HUVECs).

2 Materials and methods

2.1 Mesenchymal stromal cell isolation and culture

MSCs were isolated from bone marrow aspirates derived from patients undergoing total or partial arthroscopy after obtaining informed consent (CER Liguria: 372/2019). In brief, bone marrow was washed with phosphate-buffered saline (PBS) and centrifuged at 300 g for 10 min. Nucleated cells were seeded at a density of 1.5 × 105 cells/cm2 and cultured in α-MEM GlutaMAX medium (Gibco, Waltham, MA, United States) supplemented with 10% fetal bovine serum (FBS) (Gibco), 100 U/mL penicillin/streptomycin (Euroclone, Milan, Italy), and 1 ng/mL fibroblast growth factor 2 (FGF-2) (PrepoTech) (complete medium). Cells were maintained in a humidified incubator at 37°C and 5% CO2. After 5 days, they were washed with PBS to remove unattached cells and a fresh, complete medium was added. At 90% confluence, cells were detached using trypsin/EDTA (Euroclone) and used for further experiments. Only cells at early passages (passage 2) were used for the subsequent analyses. For all the experiments, MSCs were cultured as spheroids and monolayer culture was used as a control.

2.2 Spheroid formation and culture

Spheroids were formed using a non-adherent MSC culture that stimulates MSC aggregation. 3D Petri Dish® molds (MicroTissues Inc.) were used following the manufacturer’s instructions. 3D Petri dish molds, containing 96 circular wells with a depth of 800 µm and a diameter of 400 μm, were placed on a 24-well plate and conditioned for 24 h in a complete medium. After conditioning, either 2.5 × 105 or 1 × 105 cells re-suspended in 75 μL of the complete medium were seeded into the mold, allowing spheroid formation of approximately 2,600 cells (large spheroids) and 1,000 cells (small spheroids); 10 min after the seeding of the cells, 1 mL of the additional complete medium was added outside the 3D Petri dish. Spheroids were maintained in culture for 72 h prior to being used for further analyses. Spheroids were monitored daily under the microscope, and images were acquired using a Leica DMi1 microscope (Leica Microsystem, Wetzlar, Germany). Changes in spheroid diameter over time were measured using ImageJ software.

2.3 MSC phenotypic characterization

The surface MSC phenotype was analyzed using flow cytometry analysis. 3D Petri dishes were inverted and centrifuged at 300 g for 10 min to remove spheroids from the agar. Spheroids and 2D MSCs were then trypsinized using 2.5% Trypsin (Gibco) for 15 min, and 50,000 cells/tube were incubated with specific primary antibodies conjugated with fluorescein isothiocyanate (FITC) or phycoerythrin (PE): FITC-anti-human CD31, FITC-anti-human CD34 (BD Bioscience, San Diego, CA, United States), PE-anti-human CD90, and PE-anti-human CD105 (BD Pharmingen, San Diego, CA, United States). Cells were also stained with the correlate IgG isotype controls conjugated with FITC (BD Bioscience) or PE (BD Pharmingen). Stainings were performed for 30 min at room temperature (RT) in the dark. Samples were run on a CytoFLEX Flow Cytometer (Beckman Coulter, Brea, CA, United States), and data were analyzed using FlowJo software (BD Biosciences).

2.4 Colony-forming unit (CFU) assay

After trypsinization, MSCs derived from both monolayer and spheroid cultures were seeded on a 10-cm-diameter Petri dish at a density of 10 cells/cm2 (Banfi et al., 2000). Cells were maintained in culture for 15 days in the complete medium, with refreshment every 3 days. Cells were then fixed in 10% buffered formalin (Bio Optica, Milan, Italy) for 10 min and stained with 1% methylene blue for 45 min. Stained dishes were scanned using the Epson Perfection 1260 scanner, and colony-forming units (CFUs) were analyzed using the ImageJ ColonyArea plugin (Guzmán et al., 2014). Colonies’ cell morphology was evaluated using a Leica DMi1 inverted phase microscope (Leica).

2.5 Viability analysis

A total of 20,000 cells were seeded on a glass coverslip for 2D MSC culture, and MSC spheroids were obtained as described above. Cell viability was assessed using a LIVE/DEAD® imaging assay (Thermo Fisher Scientific), following the manufacturer’s instructions. Images were acquired using an AxioPhot microscope (Carl Zeiss, Jena, Germany, for 2D cultured MSCs) and an EVOS FL microscope (Thermo Fisher Scientific, for 3D spheroids) and analyzed using ImageJ software.

2.6 Immunofluorescence staining

Expression of HIF-1α, VEGF-A, Ki-67, cleaved caspase-3, and p21 in large and small spheroids was assessed via immunofluorescence. For 2D MSC culture, 20,000 cells were seeded on a glass coverslip, and MSC spheroids were obtained as described above. For cleaved caspase-3, a positive control was performed by treating monolayer MSCs with 10 µM etoposide for 72 h to verify antibody specificity. After 72 h of culture, microtissues were fixed in 10% buffered formalin for 10 min, embedded in OCT (Bio Optica), and stored at −80°C. Cryosections of 8 μm thickness were obtained using a Leica Cryostat. Sections were permeabilized in 0.5% Triton X-100 in PBS for 15 min. Unspecific binding was prevented by blocking the samples with 3% bovine serum albumin (BSA)/0.5% Triton X-100 in PBS for 1 h at room temperature. Sections were incubated overnight at 4°C with specific primary antibodies directed against HIF-1α (dilution 1:500, Santa Cruz Biotechnology), VEGF-A (dilution 1:500, Santa Cruz Biotechnology), cleaved caspase-3 (dilution 1:200, Cell Signaling Technology), Ki-67 (dilution 1:6400, Cell Signaling Technology), and p21 (dilution 1:800, Cell Signaling Technology); all antibodies were diluted in 1.5% BSA/0.5% Triton X-100 in PBS. After several washes in PBS/0.5% Triton X-100, samples were incubated with specific secondary goat anti-rabbit/mouse antibodies, conjugated with Alexa Fluor 488 (dilution 1:300, Life Technologies), and diluted in 1.5% BSA/0.5% Triton X-100 for 75 min at room temperature. The cytoskeleton was visualized with 5 U/mL phalloidin-594 staining (Life-Tech). Cell nuclei were stained with DAPI (dilution 1:10,000, Thermo Fisher Scientific) for 5 min at room temperature, and slides were cover-slipped with an aqueous mounting medium. Negative controls were run in parallel for all the staining by incubating the sections only with secondary antibodies to exclude any unspecific binding. Images were acquired using an Apotome microscope (Carl Zeiss) and analyzed using ImageJ software. The mean diameter of cryosections was also calculated, and only cryosections of the inner region of the spheroids were considered for analysis, while those of the outer region were excluded.

2.7 RNA extraction and quantitative real-time polymerase chain reaction

Total RNA from both MSC monolayer and spheroid cultures was collected after 72 h of culture using TRIzol™ reagent (Thermo Fisher Scientific), following the manufacturer’s instructions. The RNA concentration was quantified by measuring absorbance at 260 nm using a spectrophotometer (BioSpectrometer, Eppendorf, Milan, Italy). RNA purity was checked considering the 260 nm/280 nm ratio, with values included from 1.5 to 2.1. Complementary DNA (cDNA) was synthesized starting from 2 μg of total RNA using the SuperScript™ VILO™ Master Mix (Thermo Fisher Scientific), following the manufacturer’s instructions. Transcript levels of target genes (CDK1, HIF-1α, LAMB1, and p21) were measured using quantitative real-time polymerase chain reaction (qRT-PCR) performed on a 7500 Fast Real-Time PCR System (Applied Biosystems, Waltham, MA, United States) using the SYBR Green Master Mix (BrightGreen 2X qPCR, ABM®, Richmond, Canada). Sequence of primer used to evaluate the expression of the target genes are listed in Table 1. The housekeeping gene GAPDH was used as an endogenous control for normalization. Data were analyzed using the ΔΔCt method, and the results were shown as fold changes of the target gene with respect to a 2D condition normalized to the fold change of GAPDH.

TABLE 1
www.frontiersin.org

TABLE 1. Primers used for quantitative real-time PCR.

2.8 Cell senescence evaluation: Sudan black B staining and senescence-associated β-galactosidase activity

Senescence levels were evaluated using Sudan black B (SBB) staining, a dye capable of labeling lipofuscin (Evangelou and Gorgoulis, 2017). For 2D MSC cultures, 20,000 cells were seeded on a glass coverslip, and MSC spheroids were obtained as described above. Positive control was performed by treating MSCs with 20 μM H2O2 for 2 h, and treated cells were then maintained in the complete medium for 24 h before staining (Kiyoshima et al., 2012). Both coverslips and MSC spheroids were fixed in 10% buffered formalin for 10 min, and spheroid cryosections were obtained, as previously explained. Sections and coverslips were dehydrated in 70% ethanol and stained with a 0.7% SBB solution (Evangelou and Gorgoulis, 2017). Cell nuclei were counterstained using a 0.1% nuclear fast red solution. Slides and coverslips were mounted with a glycerol mounting medium (Dako, Santa Clara, CA, United States). β-Galactosidase activity was assessed using a Senescence β-Galactosidase Staining Kit (Cell Signaling Technology). In particular, for spheroid staining, microtissues were fixed for 10 min using the solution provided with the kit (glutaraldehyde, methanol, and formaldehyde), and spheroid cryosections were obtained as previously explained. Final staining was performed on both cryosections and coverslips, following the manufacturer’s instructions. Coverslips were mounted using an aqueous mounting medium (Dako). For spheroid cryosections, cell nuclei were stained with Mayer’s hematoxylin (Sigma-Aldrich) for 5 min at room temperature, and slides were cover-slipped with an aqueous mounting medium. Images were acquired using an AxioPhot microscope (Carl Zeiss) and analyzed using ImageJ software. The mean diameter of cryosections was also calculated, and only cryosections of the inner region of the spheroids were considered for analysis, while cryosections of the outer region were excluded.

2.9 EV separation

After 72 h of seeding, MSC cultures were rinsed twice with PBS and maintained for 30 min in α-MEM GlutaMAX medium supplemented with only 100 U/mL penicillin/streptomycin (serum-free medium, SF). Medium was replaced with fresh SF medium, and cells were maintained in an incubator at 37°C with 5% CO2 for 24 h. After cell pre-conditioning, conditioned medium (CM) were collected and centrifuged at 300 g for 10 min at 4°C to remove cells and cell debris and at 2000 g for 20 min at 4°C to remove apoptotic bodies. CM were concentrated using Amicon Ultra Filters with a molecular weight cut-off of 100 KDa (Millipore, Burlington, MA, United States). Concentrated CM were centrifuged at 10,000 g for 40 min at 4°C to pellet the fraction of EVs enriched in medium-sized vesicles (mEVs). The supernatant was centrifuged at 100,000 g for 2 h at 4°C to pellet the fraction of EVs enriched in small-sized vesicles (sEVs). Supernatants derived from 100,000 g centrifugation were collected as depleted medium (DM). EV pellets were washed with PBS and centrifuged with the same g acceleration and time used for their separation at 4°C. After the washing step, EV pellets were re-suspended in 50 μL of 0.22-μm-filtered PBS for the subsequent analysis. A Beckman Coulter ultracentrifuge (Beckman Coulter Optima XPN-100 ultracentrifuge, Beckman Coulter) was used with swinging bucket rotors of type SW55Ti. The concentration of membrane-bound proteins of isolated EVs, used to indirectly quantify EV concentration and protein amount of CM and DM, was evaluated using a bicinchoninic acid assay (BCA), following the manufacturer’s instructions (Thermo Fisher Scientific). The total protein content was normalized to the initial number of seeded cells.

2.10 Protein composition analysis of MSC-derived conditioned medium

Analysis of protein composition was performed using the Proteome Profiler Human XL Cytokine Array Kit (R&D Systems), following the manufacturer’s instructions. In brief, 15 μg of proteins of CM were incubated with a cocktail of antibodies conjugated to the surface of a membrane overnight at 4°C. After several washing steps, the membrane was incubated with a cocktail of biotin-conjugated antibodies for 1 h at room temperature, then with HRP-conjugated streptavidin for 30 min at room temperature. After treatment of the membrane with substrate necessary for the HRP reaction (luminol and hydrogen peroxide) (LiteAblot® TURBO, Euroclone), chemiluminescence was impressed on autoradiographic film (Fujifilm). Membranes were scanned using the Epson Perfection 1260 scanner, and analysis was performed using ImageJ to detect the pixel intensity. Four independent experiments were performed. Intensity values were normalized to the control spots of the membrane. Gene Ontology (GO) was performed on proteins with a normalized pixel intensity ≥0.3 using DAVID, focusing on the “Biological process (BP)” term. Graphical representations of cytokine array analysis and GO were carried out using R libraries: pheatmap, ggplot2, and circlize (Gu et al., 2014), setting the adjusted p-value < 0.05 as significant.

2.11 VEGF quantification on conditioned medium

Conditioned medium for MSCs derived from the three different culture conditions were collected, and VEGF was quantified via enzyme-linked immunosorbent assay (ELISA) using the Human VEGF ELISA Kit (Thermo Fisher Scientific), following the manufacturer’s instructions. Data were normalized to the initial number of seeded cells.

2.12 Western blot analysis

For Western blot analysis, after separation, mEVs and sEVs were re-suspended in 50 μL of RIPA buffer (1% Nonidet P-40, 0.1% sodium deoxycholate, and 0.1% sodium dodecyl sulfate in PBS, pH 7.5) supplemented with protease inhibitors (Cell Signaling Technology, Danvers, MA, United States). Cell lysates, used as controls, were obtained by lysing 60,000 MSCs for monolayer culture and one microtissue for each condition in 100 μL of RIPA buffer supplemented with protease inhibitors. After protein quantification with the BCA assay, electrophoresis was run by loading 1 μg of proteins onto 4%–12% Tris-Glycine gels (Invitrogen), and proteins were finally transferred on a PVDF membrane (Millipore). After blocking with 5% non-fat dry milk in TTBS (500 mM NaCl, 0.1% v/v Tween, and 20 mM Tris, pH 7.5) for 1 h at room temperature, membranes were incubated overnight at 4°C with specific primary antibodies diluted in 2.5% non-fat dry milk in TTBS and directed against syntenin-1 (dilution 1:1000, ab133267, Abcam), flotillin-1 (dilution 1:10,000, ab133497, Abcam, Cambridge, United Kingdom), CD63 (dilution 1:500, 10628D, Invitrogen, Waltham, MA, United States), CD81 (dilution 1:1000, 555675, BD Biosciences), Grp94 (dilution 1:1000, ab3674, Abcam), and calnexin (dilution 1:1000, Abcam). After several washes in TTBS, membranes were incubated with specific HRP-conjugated goat anti-rabbit/mouse antibodies diluted in 2.5% non-fat dry milk in TTBS (dilution 1:2000, Cell Signaling Technology, Danvers, MA, United States) for 1 h at room temperature. The detection of proteins was performed by providing the substrate for the HRP reaction (LiteAblot® TURBO, Euroclone), and chemiluminescence was impressed on autoradiographic film (Fujifilm, Tokyo, Japan). Images were scanned using the Epson Perfection 1260 scanner.

2.13 Nanoparticle tracking analysis

EV samples were characterized and quantified using ZetaView (Particle Metrix GmbH, Germany) Nanoparticle Tracking Analysis (NTA), equipped with a sample cell and two lasers (488 nm and 640 nm). ZetaView 8.05.14_SP7 software was used. After calibration with 100 nm polystyrene beads, samples were diluted in 0.22-μm-filtered PBS 1X and injected using a 1-mL syringe. Size distribution analyses of 11 different positions were performed for each sample on at least three different sEV preparations. The particle number was then normalized to the initial number of seeded cells.

2.14 Non-conventional flow cytometry analysis of EVs

For flow cytometry characterization, EVs were suspended in 0.22-μm-filtered PBS/EDTA, distributed in flow cytometry tubes (1 × 108 EVs/tube), and stained with 1 μM CFDA-SE (VybrantTM CFDA SE Cell Tracer Kit, Thermo Fisher Scientific) at RT to visualize intact vesicles. Since CFDA-SE diffuses within vesicles only at RT, a control at 4°C was also included. A mixture of fluorescent beads with a diameter ranging from 100 to 900 nm (Megamix-Plus FSC and Megamix-Plus SSC, Biocytex, France) was used to discriminate EV size. Expression of typical EV markers, such as CD9 (APC Mouse Anti-Human CD9, 312108, BioLegend, San Diego, CA, United States), CD63 (PE-CF594 Mouse Anti-Human CD63, 565403, BD Biosciences), and CD81 (BV421 Mouse Anti-Human CD81, 740079, BD Biosciences), was evaluated within the CFDA-SE positive events. All the evaluations were performed using a CytoFLEX flow cytometer (Beckman Coulter), and data were analyzed using FlowJo. Dimensional analysis was performed by creating three dimensional gates (EVs <100 nm, 100 nm < EVs <160 nm, and 160 nm < EVs <900 nm).

2.15 Transmission electron microscopy

The EV preparations from 2D cultures and from 3D spheroids of MSCs were re-suspended in 20 μL PBS (pH 7.4) and fixed by adding an equal volume of 2% paraformaldehyde in 0.1 mol/L phosphate buffer (pH 7.4). EVs were then adsorbed for 10 min onto formvar-carbon-coated copper grids by floating the grids on 4 μL drops on parafilm. Subsequently, grids with adhered vesicles were rinsed in PBS and negatively stained with aqueous 2% uranyl acetate for 5 min at room temperature. Stained grids were then embedded in ice-cold 2.5% methylcellulose for improved preservation and air-dried before examination. Electron micrographs were obtained using a Hitachi 120 kV transmission electron microscope (HT7800 series, Tokyo, Japan) equipped with a MegaView III digital camera and RADIUS 2.0 software (EMSIS, Muenster, Germany).

2.16 HUVEC culture

HUVECs (ATCC® CRL-1730™) were plated on a T75-cm2 flask and cultured in F12-K medium (ATCC®), supplemented with 30 mg/mL of endothelial cell growth supplement (Sigma-Aldrich, St. Louis, MI, United States), 10% FBS, 100 U/mL of penicillin/streptomycin, and 0.1 mg/mL of heparin (Pharmatex) (complete medium). For the following test, HUVECs were used at passage number 5. For each test, positive and negative controls were carried out by incubating cells in complete and F12-K media, supplemented with 100 U/mL of penicillin/streptomycin (serum-free medium, SF medium), respectively. The effect of different fractions of the MSC-derived secretome on HUVEC proliferation, migration, and tube formation was evaluated.

2.17 HUVEC proliferation assay

HUVECs were seeded in a 96-well plate (10,000 cells/well) and cultured in the complete medium for 24 h to allow cell attachment. To synchronize the cell cycle, HUVECs were incubated with the SF medium for 2 h. HUVEC treatments were performed for 24 h using different secretome fractions (CM, DM, mEVs, and sEVs) derived from 30,000 MSCs. Cell proliferation was then measured using the BrdU proliferation assay, following the manufacturer’s instructions (Cell Proliferation ELISA, BrdU colorimetric, Roche, Basel, Switzerland). Data were normalized to a negative control.

2.18 HUVEC migration assay

To evaluate MSC-derived CM, DM, and EV ability to induce HUVEC migration, a Transwell system (Corning, New York, United States) was used. Transwells, with a pore size of 8 μm, were placed in a 24-well plate, and 50,000 cells re-suspended in the SF medium were seeded inside the inserts. Secretome fractions (CM, DM, mEVs, and sEVs) derived from 150,000 MSCs of the different culture conditions were put in the lower chamber. After 9 h of incubation at 37°C and 5% CO2, inserts were washed with PBS, and the non-migrated cells in the upper part of the membrane were gently removed using a cotton swab. Migrated cells were fixed in 10% buffered formalin for 10 min and stained with 1% methylene blue for 30 min. Images were acquired using a Leica DMi1 inverted phase microscope, and analysis was performed on at least five different areas using the ImageJ Cell Counter plugin. Data were normalized to a positive control.

2.19 Tube formation assay

Wells of a 96-well plate were coated with 50 μL of Matrigel (Corning) and incubated at 37°C for 30 min to allow Matrigel polymerization. Later, 20,000 HUVECs, re-suspended in the SF medium, were seeded on the top of the Matrigel and incubated at 37°C and 5% CO2 for 6 h in the presence of different secretome fractions (CM, DM, mEVs, and sEVs) derived from 60,000 MSCs. Images were acquired using a Leica DMi1 inverted phase microscope, and analysis was performed on at least three different areas using the ImageJ Angiogenesis Analyzer plugin (Carpentier et al., 2020). Data were normalized to a positive control.

2.20 Statistical analysis

Data were analyzed using GraphPad Prism 8.0 (GraphPad Software, Inc.), and they are shown as mean ± standard deviation (SD). Statistical analysis of differences between multiple groups was performed using a two-way or one-way ANOVA, followed by Tukey’s multiple comparison test. The level of significance was set at p < 0.05.

3 Results

3.1 MSCs in spheroid cultures maintain stemness characteristics but exhibit morphological alterations

Compared to conventional 2D culture, spheroids represent a promising alternative (Langan et al., 2016). On this basis, we evaluated whether spheroid size might influence the behavior and paracrine activity of human MSCs.

Two different sizes of spheroids were generated using human MSCs: large (diameter ∼300 μm) and small spheroids (diameter ∼200 μm), which originated from approximately 2,600 and 1,000 cells, respectively (Figure 1A). 2D cultured MSC show a classical fibroblast like spindle-shaped cells morphology (Figure 1B). Spheroid diameter was measured over time, and interestingly, both large and small spheroids showed a significant size decrease during culture (Figure 1C; Supplementary Table S2). In addition, we highlighted size stabilization in small spheroids after 72 h of culture, while the diameter of large spheroids showed a constant decreasing trend until 120 h of culture, when the diameter became stable.

FIGURE 1
www.frontiersin.org

FIGURE 1. Spheroid dimensional analysis. (A) Representative images of large and small spheroids at different culture time points (24, 48, 120, and 192 h); scale bar = 500 µm. (B) Representative image of 2D MSC culture used as the control; scale bar = 500 µm. (C) Dimensional analysis of spheroid diameter variation along culture time. Diameter was calculated on 100 spheroids for each time point and replicate using ImageJ. Data are represented as the mean ± SD. **** p-value <0.0001 (N = 4, two-way ANOVA and Tukey’s multiple comparison).

To verify that stemness properties were not affected by the culture method, surface phenotype, clonogenic potential, and trilineage differentiation potential (Supplementary Figure S1) were further investigated. The surface phenotype was evaluated using flow cytometry, assessing putative negative markers CD31 and CD34 and positive markers CD90 and CD105. No statistically significant alteration was observed between different experimental groups (Figures 2A, B), demonstrating that cells retain their typical phenotype independently from the culture condition.

FIGURE 2
www.frontiersin.org

FIGURE 2. MSC phenotype and self-renewal ability evaluation. (A) Representative histograms reporting flow cytometry analysis of MSCs cultured in 2D (purple), large (green), or small (orange) spheroids. Histogram curves identify cells reacting with CD31-, CD34-, CD90-, and CD105-specific antibodies. Area under the gray curve identifies the reaction of the cells with the correspondent non-reactive immunoglobulin of the same isotype (isotype control). Data are representative of three independent experiments. (B) Flow cytometry analysis of MSCs cultured in 2D, large, or small spheroids. Histogram shows the average percentage of positive events for each marker. Data are represented as the mean ± SD (N = 3). (C) Representative images of CFU-f (left) and a magnification showing the colony cell morphology (right); scale bar = 200 µm. (D) Number of colonies and (E) area of the Petri dish occupied by colonies of MSCs derived from 2D, large, and small spheroid cultures. Data are represented as the mean ± SD (N = 3).

The clonogenic potential, defined as the ability to form colonies starting from a single cell (CFU-f: colony-forming unit), was investigated. MSCs from 2D, large, and small spheroid culture conditions were seeded at low density, and their macro- and microscopic morphology was evaluated after 2 weeks. Microscopic morphology showed significant differences between different culture conditions (Figure 2C). MSC-derived colonies cultured in 2D showed fibroblast-like, spindle-shaped cell morphology with larger cells compared to spheroid-derived CFUs, which have a smaller size and elongated spindle-shaped morphology with a more compact organization. Interestingly, this peculiarity was enhanced in CFU derived from small spheroids. Furthermore, the total number of colonies was analyzed, and no significant difference was observed (Figure 2D). Finally, the surface area covered by colonies was analyzed (Figure 2E), and no significant difference was observed.

3.2 Spheroid size influences the establishment of a hypoxic core

Cell response to hypoxic stress is mainly mediated by hypoxia-inducible factors (HIFs). To study the influence of diameter on the onset of hypoxia in spheroids, mRNA expression levels of the transcriptional factor hypoxia-induced factor-1α (HIF-1α) were investigated using qRT-PCR (Figure 3A). HIF-1α was strongly overexpressed in large spheroids, compared to both 2D-cultured MSCs and small spheroids, where no differences were highlighted. To assess the localization of HIF-1α inside the 3D structure, its expression was further analyzed using immunofluorescence. Figure 3B shows the quantification of HIF-1α-positive cells, and Figure 3F shows a localization of HIF-1α in the central core of large spheroids, as highlighted by white arrows. No HIF-1α expression was found in small spheroids, suggesting that a small diameter may support the generation of a less strong gradient, allowing nutrients and oxygen to reach the central region.

FIGURE 3
www.frontiersin.org

FIGURE 3. Spheroid size influences the generation of a hypoxic core. (A) Quantitative real-time PCR for the expression of HIF-1α (hypoxia-inducible factor alpha). HIF-1α expression of the 2D culture was compared to the expression in the 3D conditions. Data are represented as the mean ± SD. *** p-value <0.001 (N = 3, one-way ANOVA). Quantification of (B) HIF-1α-, (C) ki67-, and (D) cleaved caspase-3-positive cells on immunostaining. Positive cells have been counted considering 10 different regions of interest (ROI) for each replicate, with a mean diameter of 342.29 ± 20.90 µm for large spheroids and 194.43 ± 9.20 µm for small spheroids. Data are represented as the mean ± SD. *** p-value <0.001 and **** p-value <0.0001 (N = 3, one-way ANOVA). (E) Quantification of VEGF in the conditioned medium using the ELISA assay. Data are represented as the mean ± SD. ** p-value <0.01 and *** p-value <0.001 (N = 3, one-way ANOVA). Representative images of immunostaining for (F) HIF-1α, (G) VEGF-A, (H) ki67, and (I) cleaved caspase-3 on different culture conditions. For all the images from the top, pictures show the nuclei (DAPI), phalloidin, and the investigated markers mentioned before, along with a merged view; scale bar = 25 µm. In all the images, white arrows indicate positive cells.

Based on the different expression levels of HIF-1α, the expression of VEGF-A, a marker deeply interconnected with HIF-1α and hypoxia, was also investigated. VEGF secretion in the conditioned medium was also quantified using ELISA (Figure 3E). A significantly higher concentration of VEGF was observed in both large and small spheroids compared to 2D-derived CM, suggesting that spheroid culture can stimulate a higher secretion of VEGF compared to the classical monolayer culture. Expression of VEGF-A was also investigated via immunofluorescence staining (Figure 3G), which highlighted a strong signal in the internal regions of large spheroids. On the contrary, a weak but still detectable signal was observed in small spheroids and totally absent in 2D cultures. All these data together suggest that VEGF-A secretion is partially dependent on the culture condition.

In addition, the expression of the proliferative marker, ki67, was also investigated (Figures 3C, H). A decreased expression was found in large spheroids, while enhanced ki67 expression was observed in small spheroids, particularly in the external regions (white arrows in Figure 3H). However, both large and small spheroids exhibited lower ki67 expression compared to the 2D cultured control, suggesting a reduced proliferative capability in 3D culture compared to monolayer culture.

On the other hand, the expression of cleaved caspase-3 was also investigated (Figures 3D, I) to assess whether 3D culture influences levels of apoptosis. Staining specificity was confirmed by treating MSCs for 72 h with 10 μM etoposide, a molecule that can induce apoptosis, as a positive control. Interestingly, similar levels of apoptosis were observed between the different conditions, suggesting that this alternative culture method does not enhance the apoptotic process. Cell death was also assessed via the LIVE/DEAD assay (Supplementary Figure S2), showing that no PI-positive cells are present in the structure.

3.3 Spheroid culture influences cellular senescence in MSCs

Culture conditions may have an important influence on cellular senescence onset. Senescence was investigated by SBB staining, β-galactosidase activity, and qRT-PCR on senescence-related genes. SBB is a dye used to stain lipofuscin, a lysosomal product accumulated inside the cells following lysosomal dysfunction, usually occurring in senescent cells (Figure 4A). A positive control (Ctr pos), obtained by treating MSCs with 20 μM H2O2, was also carried out to confirm dye specificity. Quantitative analysis (Figure 4B) showed that SBB-positive cell number significantly increased in large spheroids compared to 2D cultures, suggesting that spheroid size may influence cellular senescence. On the other hand, levels of cellular senescence are decreased in small spheroids compared to both 2D and large spheroid cultures.

FIGURE 4
www.frontiersin.org

FIGURE 4. Spheroid culture influences cellular senescence in MSCs. (A) Representative images of Sudan black B and β-galactosidase staining of the different culture conditions; scale bar = 50 µm. Ctr pos, positive control. (B) Quantification of SBB- and β-galactosidase-positive cells. SBB- and β-galactosidase-positive cells have been counted considering 10 different regions of interest for each replicate with a mean diameter of 338.14 ± 19.62 µm for large spheroids and 187.43 ± 13.26 µm for small spheroids. Data are represented as the mean ± SD. * p-value <0.05, ** p-value <0.01, and **** p-value <0.0001 (N = 4, one-way ANOVA). (C) Representative images of immunostaining for p21 on different culture conditions. From the top, pictures show the nuclei (DAPI), phalloidin, p21, and their merge view; scale bar = 25 µm. (D) Quantification of p21-positive cells on immunostaining. Positive cells have been counted considering 10 different regions of interest for each replicate, with a mean diameter of 332.12 ± 18.40 µm for large spheroids and 193.35 ± 7.22 µm for small spheroids. (E) Quantitative real-time PCR for the expression of p21, LAMB1 (lamin subunit beta 1), and CDK1 (cyclin-dependent kinase 1); data are represented as the mean ± SD. * p-value <0.05, ** p-value <0.01, and *** p-value <0.001 (N = 3, one-way ANOVA).

β-Galactosidase is an enzyme usually expressed by senescent cells. The evaluation of β-galactosidase activity is one of the gold standard methods used to evaluate the presence of senescent cells. Positive control, as also observed from SBB staining, shows the presence of a high percentage of senescent cells (Figure 4A), confirming staining specificity. Quantitative analysis (Figure 4B) revealed similar results to SBB staining. A higher number of β-galactosidase-positive cells were identified in large spheroids compared to 2D, while a comparable number of β-galactosidase-positive cells were found in small spheroids and the 2D condition.

Similar results were obtained by analyzing the expression of p21 using immunostaining (Figures 4C, D), a protein usually upregulated in senescent cells, where it inhibits the formation of the CDK2–cyclin E complex involved in the G1-S cell cycle checkpoint (González-Gualda et al., 2021; Schmitt et al., 2022). A higher number of p21-positive cells were observed in large spheroids compared to the 2D culture, while a similar number of p21-positive cells were found in small spheroids and the 2D culture. In addition, we highlighted that a percentage of senescent cells can also be detected in the 2D culture in all the performed analyses for SBB, β-galactosidase, and p21 staining (Figures 4A–D). This may be due to the old age (average 75 years) of MSC donors, which could explain the co-isolation of a senescent fraction of cells (Gnani et al., 2019; Kaminska et al., 2021).

In addition, cellular senescence was further investigated by assessing mRNA levels of genes normally dysregulated during the senescent process, such as p21, lamin subunit B1 (LAMB1), and cyclin-dependent kinase-1 (CDK1), using real-time PCR (Figure 4E). mRNA levels of p21 were upregulated in large spheroids compared to both 2D culture and small spheroids, where similar expression levels of p21 were found. LAMB1, a key component of the nuclear lamina, is another key hallmark of cellular senescence since it is usually downregulated in senescent cells, leading to the loss of nuclear integrity and DNA leakage in the cytoplasm (González-Gualda et al., 2021; Kumari and Jat, 2021). A slight, not-significant downregulation of this marker was observed in large spheroids compared to the 2D culture. Intriguingly, small spheroids showed a strong upregulation of LAMB1 compared to both 2D culture and large spheroids. Finally, no significant variations were observed in the expression of CDK-1.

3.4 MSC spheroid-derived secretome showed upregulation of pro-angiogenic and inflammatory cytokines

In order to evaluate the secretome composition, CM were analyzed using the XL Proteome Profiler to evaluate differently secreted molecules between the three culture conditions, and important differences were highlighted (Figure 5A; Supplementary Figure S3). Compared to 2D, large spheroid-derived secretome showed upregulation of different cytokines, such as fibroblast growth factor-19 (FGF-19), growth differentiation factor-15 (GDF-15), interleukin-8 (IL-8), monocyte chemoattractant protein-1 (MCP-1), macrophage migration inhibitory factor (MIF), matrix metallopeptidase-9 (MMP-9), osteopontin (OPN), stromal derived factor-1-alpha (SDF-1α), and VEGF-A, and the downregulation of pentraxin-3 (PTX3).

FIGURE 5
www.frontiersin.org

FIGURE 5. Secretome profiling. (A) Quantification of normalized pixel intensity for each modulated protein on cytokine arrays of CM derived from 2D, large, and small spheroids MSC culture. Data are represented as the mean ± SD; * p-value <0.05, ** p-value <0.01, *** p-value <0.001, and **** p-value <0.0001 (N = 4, one-way ANOVA). ANG, angiogenin; FGF-19, fibroblast growth factor 19; GDF-15, growth differentiation factor 15; IL-6, interleukin 6; IL-8, interleukin 8; MCP-1, monocyte chemoattractant protein 1; MIF, macrophage migration inhibitory factor; MMP-9, matrix metallopeptidase 9; OPN, osteopontin; PTX3, pentraxin 3; SDF1a, stromal cell-derived factor 1; THBS-1: thrombospondin 1; VEGF-A, vascular endothelial growth factor-A; HGF, hepatocyte growth factor. (B) Bubble plot of Gene Ontology analysis of the BP term on modulated cytokines identified with cytokine arrays. (C) Chord plot showing the cytokine involvement in the different biological processes identified with Gene Ontology analysis.

In addition, compared to small spheroid-derived secretome, large-spheroid-derived secretome showed upregulation of GDF-15, OPN, PTX3, THBS-1, and VEGF-A and downregulation of interleukin-6 (IL-6) and HGF (p = 0.08).

On the other hand, compared to 2D, small spheroid-derived secretome exhibited an upregulation of ANG, IL-6, VEGF-A, and MIF and a downregulation of PTX3.

Furthermore, as reported in Supplementary Figure S4, the different biological replicates analyzed have similar behavior underlying data consistency. In addition, the dendrogram located on top of the heatmap highlights that all the biological replicates derived from a specific condition are clustered in the same branch.

In Figure 5B, the main biological processes identified by GO analysis in different comparisons are presented, and in Figure 5C, the main factors involved in these different biological processes are shown.

3.5 MSC spheroids secrete a higher quantity of EVs compared to 2D

EVs were isolated from the conditioned medium (CM) obtained by maintaining MSCs in the SF medium for 24 h at 37°C and 5% CO2. The BCA assay revealed a higher protein concentration in sEVs derived from large and small spheroids, compared to 2D-derived sEVs, with no significant difference between large and small spheroid-derived sEVs (Figure 6A). No significant difference was highlighted in mEV protein concentration.

FIGURE 6
www.frontiersin.org

FIGURE 6. MSC spheroids secrete higher quantity of EVs compared to 2D. (A) Quantification of protein content on mEVs and sEVs using the BCA assay; data are normalized to the cell number (µg/106 cells). Data are represented as the mean ± SD; * p-value <0.05 and ** p-value <0.01 (N = 4, two-way ANOVA and Tukey’s multiple comparison). (B) NTA of particle concentration, expressed as particle/106 cells. Data are represented as the mean ± SD; * p-value <0.05, *** p-value <0.001, and **** p-value <0.0001 (N = 4, two-way ANOVA and Tukey’s multiple comparison). (C) Representative NTA size distribution analysis of mEVs and sEVs (N = 4). (D) Western blot analysis on 2D, large, and small spheroid MSC culture-derived mEVs and sEVs; cell lysates were also loaded as controls. Specific expression of CD63, CD81, flotillin-1, syntenin-1, and GRP94 was investigated. (E) Quantification on flow cytometry analysis of CD9-, CD63-, and CD81-positive events in the CFDA-SE gate. Data are shown as ratio between mean fluorescence intensity (MFI) of EVs stained with specific antibody and the correspondent isotype control (relative MFI) (N = 3). (F) TEM micrographs of isolated mEVs and sEVs. Scale bar = 200 nm. White arrows indicate EVs.

NTA of mEVs and sEVs (Figures 6B, C) showed an increased concentration of mEVs derived from small spheroids and a higher concentration of sEVs derived both from large and small spheroid MSCs compared to 2D. Furthermore, a higher concentration of sEVs was also observed in small spheroids compared to large spheroids. NTA of EV size showed a consistent reduction of dimension in sEVs derived from all conditions (Table 2) with no significant alteration between different culture conditions.

TABLE 2
www.frontiersin.org

TABLE 2. NTA of size of mEVs and sEVs isolated from 2D, large, and small spheroids. Data are reported as the mean ± SD (N = 4). mEVs: medium-sized EVs, sEVs: small-sized EVs.

Western blot analysis of mEVs, sEVs, and cell lysate (Figure 6D) was used to investigate the expression of classical vesicular markers (flotillin-1, syntenin-1, CD63, and CD81). As expected, the negative control GRP94, a chaperon needed to assist the folding and assembly of secreted and membrane proteins in ER (Eletto et al., 2010), was expressed only in cell lysate. Syntenin-1, a classical marker of sEVs, was expressed only in this population in all the culture conditions. Flotillin-1 was expressed in both EV populations and cell lysate. Interestingly, the expression of tetraspanins CD63 and CD81 was detected in sEVs derived from all culture conditions. CD81 was not observed in any population of mEVs, while CD63 was poorly observable only in mEVs derived from 2D and large spheroids and not observed in mEVs derived from small spheroids.

EVs were also analyzed using a non-conventional flow cytometry approach to evaluate the expression of tetraspanins (CD9, CD63, and CD81). This approach is based on the use of fluorescent dimensional beads of known diameter (Supplementary Figure S6A) and the lipophilic tracer, CFDA-SE, that passively diffuses inside EVs and interacts with intraventricular enzymes at RT (Supplementary Figure S6B). An appropriate dimensional gate was designed using EVs stained with CFDA-SE at 4°C, and the expression levels of CD9, CD63, and CD81 were assessed in the CFDA-SE-positive events. As reported in Figure 6E, CD9 is the least expressed marker, while expression levels of CD63 and CD81 were similar in different culture conditions. Moreover, similar levels of tetraspanin expression were observed between mEVs and sEVs. Dimensional analysis was performed by creating three different gates (EVs<100 nm, 100 nm<EVs<160 nm, and 160 nm<EVs<900 nm). As shown in Table 3 and Supplementary Figure S6D, an enrichment of the EVs <100 nm was observed in all the conditions, particularly in sEVs. Similar levels of EVs in the gate 100 nm<EVs<160 nm were present in the different EV preparations, and a small percentage of EVs with dimensions higher than 160 nm was observed. EVs were further characterized by transmission electron microscopy (Figure 6F). TEM micrographs confirmed that all the EV preparations are enriched in intact vesicles surrounded by a lipid bilayer (white arrows in Figure 6F).

TABLE 3
www.frontiersin.org

TABLE 3. Percentage of CDFA-SE-positive events in the considered three-dimensional gates (EVs <100 nm, 100 nm<EVs<160 nm, and 160 nm<EVs<900 nm). Data are reported as the mean ± SD (N = 3). mEVs: medium-sized EVs, sEVs: small-sized EVs.

3.6 Spheroid culture increased the angiogenic potential of MSC-derived soluble factors and EVs

The angiogenic potential of the different secretome fractions (CM, DM, mEVs, and sEVs) derived from 2D, large, and small spheroids was evaluated on HUVECs by assessing their proliferation, migration, and ability to create tubes. Cell treatments were performed using a known amount of the different secretome fractions based on the cell number, and HUVECs were treated in all tests, maintaining a cell ratio of 1:3 (HUVEC:MSC). First, the effect of the different secretome fractions on HUVEC proliferation was investigated (Figure 7A). Interestingly, only small spheroid-derived secretome showed a strong effect on proliferation induction. More deeply, both CM and DM were able to induce significantly higher HUVEC proliferation compared to 2D and large spheroids. In addition, small spheroid-derived sEVs also showed a proliferative effect significantly more robust compared to 2D-derived sEVs, suggesting a synergistic effect influenced by both vesicular and soluble components. On the contrary, mEVs derived from all culture conditions did not show any effect on HUVEC proliferation.

FIGURE 7
www.frontiersin.org

FIGURE 7. Effects of MSC secretome fractions on HUVEC proliferation, migration, and tube formation. (A) HUVEC proliferation was assessed via the BrdU incorporation assay; (B) migration was assessed via the Transwell migration assay. (C) Quantitative analysis of number of branches in the tube formation assay. Data are represented as the mean ± SD. * p-value <0.05, ** p-value <0.01, *** p-value <0.001, and **** p-value <0.0001 (N = 3, two-way ANOVA and Tukey’s multiple comparison). Ctr pos, positive control; Ctr neg, negative control; CM, conditioned medium; DM, depleted medium; mEVs, medium-sized EVs; sEVs, small-sized EVs. (D) Representative images of the HUVEC tube formation assay; scale bar = 500 µm.

Endothelial cells are able to migrate toward chemoattractant signals. A Transwell system was used to evaluate the migratory ability of HUVECs by inserting the different secretome fractions in the lower chamber and the cells in the upper chamber. Small spheroid-derived CM showed a stronger effect on HUVEC migration compared to large spheroid-derived CM, while the 2D-derived CM had low activity (Figure 7B). Small spheroid-derived DM showed a reduced effect compared to the CM. As observed before, small spheroid-derived sEVs induced a higher migratory effect on HUVECs compared to sEVs derived from the other conditions.

A tube formation assay was used to evaluate this specific aspect of the angiogenic process that starts only in the later stages, when cells migrate in the region of interest, proliferate to increase their number, and then start to create vessels (Adair and Montani, 2010). Different from the other tests, similar effects of CM and DM derived from the two spheroid cultures were observed (Figures 7C, D). In this assay, a significantly higher effect of both small spheroids-derived mEVs and sEVs was also observed.

Collectively, these results indicated that a possible synergic action of soluble factors and EVs can enhance the angiogenic potential of the small spheroid-derived secretome.

4 Discussion

The role of MSC-derived soluble factors and EVs in angiogenesis is crucial for the development of new therapeutic approaches, which may be used for the treatment of vasculogenesis-related disorders. In this study, we investigated how a different culture system could affect MSC properties and their angiogenic potential without the use of external priming such as hypoxia or cytokine treatments. Although these approaches are often successfully used to enhance specific MSC therapeutic potential, they present some limitations, including induction of immunogenicity, high costs, MSC source-related variability, and a lack of standardized protocols (Noronha et al., 2019; Noronha et al., 2021).

In this study, MSCs were cultured as spheroids of different sizes. Cellular senescence, secretome composition, EV secretion, and angiogenic potential were investigated to evaluate whether these aspects could be affected by the use of spheroid culture and the size of the spheroid.

First, the culture conditions do not affect the stemness of MSCs, which maintains their superficial phenotype and clonogenic potential, while cell morphology was influenced in terms of size reduction during culture time, which may be due to a cytoskeleton reorganization (Tsai et al., 2015; Cesarz and Tamama, 2016).

The cell morphology difference was mirrored in CFU-forming MSCs, where cells derived from 2D standard conditions showed enlarged morphology, in contrast to 3D conditions, especially in small spheroids, where smaller and extremely elongated cells were found. MSCs usually possess a classical fibroblastic morphology with spindle-shaped cells, which are lost when senescence occurs (Turinetto et al., 2016). As observed in our data, MSCs forming CFU under the 3D condition show a change in morphology compared to those under the 2D condition. This morphology alteration has already been described as a normally occurring event in MSC-cultured spheroids. The size of MSCs in spheroids is drastically smaller than cells in the 2D monolayer, resulting in a 75% reduction in individual cell volume (Cesarz and Tamama, 2016). It has been reported that the process of spheroid formation is characterized by an active cell rearrangement, where MSCs first form a loose network, followed by a gradual cell condensation over a central core. Spheroids compress progressively in the first few days due to a reduction in the amount of cytoplasm and cell volume (Bartosh et al., 2010).

The process of spheroid formation is responsible for the creation of a peculiar architecture, which generates a heterogeneous microenvironment characterized by gradients of different soluble molecules, including gases and nutrients, more accurately mimicking the in vivo milieu (Sart et al., 2014). Among these molecules, oxygen has been considered a key factor in regulating cell functions, including cell proliferation (Cheema et al., 2008). In more detail, an oxygen gradient is formed in the spheroids, with (i) a surface with a high concentration of oxygen, populated by cells with a high proliferation rate, (ii) a middle zone with a low concentration of oxygen, characterized by viable quiescent cells, and (iii) a hypoxic core in the central area, with necrotic and apoptotic cells. Oxygen availability is also known to regulate HIFs, which orchestrate the expression of hypoxia-associated genes involved in the regulation of stem cell properties, development, and metabolism (Grayson et al., 2007). We here demonstrated that HIF-1α is downregulated in small spheroids compared to large spheroids, suggesting that nutrients can diffuse more efficiently and reach the core. Therefore, the use of small spheroids could overcome this critical constraint, preventing this gradient formation (Mueller-Klieser, 1984).

The specific 3D organization also influences the diffusion of metabolic wastes, growth factors, and cytokines produced by MSC that might be translated into their increased internal levels (Sart et al., 2014). This could explain the upregulated expression of VEGF-A in large spheroids compared to small spheroids. Furthermore, substantial evidence has also reported that the hypoxic environment may be responsible for the upregulation of some growth factors, such as VEGF, FGF-2, and HGF (Bhang et al., 2011).

The effects of spheroid size on cell proliferation and apoptosis were also investigated. Interestingly, apoptotic levels were not enhanced in both spheroid conditions compared to the 2D culture. As reported in the literature, the short-term culture used in this study (96 h) and the dimension of the spheroids may not be sufficient to enhance the apoptotic process (Jauković et al., 2020). On the other hand, we observed that ki67 was strongly downregulated in large spheroids compared to small spheroids, suggesting impairment in cell proliferation. The literature indicates that less than 1% of 3D-cultured MSCs are actively dividing, with 90% arrested in the G0/G1 phase of the cell cycle (Bartosh et al., 2010). In addition, the block of cell proliferation in large spheroids may be related to different signals coming from suffering and senescent cells in the inner region that became stressed due to the lack of nutrients and oxygen, arresting the proliferative ability of more external cells.

Senescence is usually correlated with metabolic dysfunctions, leading to an alteration in MSC secretory ability, characterized by higher secretion of pro-inflammatory molecules capable of stimulating cellular senescence in neighboring cells (Marques et al., 2017). The presence of a higher number of senescent cells in large spheroids was confirmed by an increased activity of β-galactosidase (Hernandez-Segura et al., 2018) and accumulation of lipofuscin, linked with lysosomal dysfunction (Kumari and Jat, 2021).

Low proliferation levels in large spheroids are in accordance with the overexpression of p21 observed, a protein involved in cell cycle arrest and usually upregulated at the onset of replicative senescence (González-Gualda et al., 2021). In addition, we observed a downregulation of LAMB1, a protein needed to ensure nuclear integrity and avoid the leakage of DNA in the cytosol, which can enhance cellular damage and fuel SASP secretion (Kumari and Jat, 2021). These findings confirmed that the establishment of a hypoxic region in the spheroids could increase cell damage and lead to the onset of replicative senescence.

Conversely, the absence of a strong hypoxic core in small spheroids is accompanied by decreased cellular senescence, as evidenced by reduced β-galactosidase activity and lipofuscin accumulation, in addition to p21 downregulation. The upregulation of LAMB1 in small spheroids, usually downregulated in senescent cells (Freund et al., 2012), may minimize DNA loss in the cytoplasm, where it can fuel SASP secretion and enhance cellular damage.

Metabolic and redox alterations were also recently linked to the onset of replicative senescence. From the observed data, one hypothesis could be that in large spheroids, the hypoxic core led to an increase in ROS production, thereby stimulating the onset of replicative senescence (Yuan et al., 2020; Jeske et al., 2021). Given the significant difference highlighted between the two spheroid cultures, we suggest the adoption of small spheroids for the maintenance of MSCs, capable of reducing the onset of cellular senescence.

Different levels of cellular senescence are also mirrored by secretome composition, which is evaluated using an XL Proteome Profiler (Fernandez-Rebollo et al., 2020). Despite the limitation of the low number of cytokines analyzed in this method, we have observed an upregulation of cytokines known to be part of the SASPs (Basisty et al., 2020; Birch and Gil, 2020; Kumari and Jat, 2021) in large spheroids, such as GDF-15, IL-8, MCP-1, MMP-9, and THBS-1. As observed from GO analysis, the activation of the specific pathways ERK1 and ERK2, normally regulating cellular proliferation, might also be involved in the premature onset of cellular senescence in our large spheroid model (Zou et al., 2019). On the contrary, small spheroid-derived secretome exhibited expression of these cytokines comparable to 2D-cultured MSCs, confirming the hypothesis that the establishment of a hypoxic core is detrimental in terms of cellular senescence induction, while size reduction may help overcome this issue. As a confirmation of this hypothesis, GO analysis on large spheroids showed the presence of terms such as “inflammatory response,” confirming the upregulation of cytokines with pro-inflammatory activity, including the SASP factors previously cited, similarly downregulated under 2D and small spheroid conditions.

Remarkably, both 3D secretomes showed upregulation of MIF, a cytokine with proved rejuvenating effects on aged MSCs (Xia et al., 2015), and positive effects on cell proliferation and survival (Y. Zhang et al., 2019). The effect on cell proliferation was indeed detected only in small spheroids, where ki67 expression was maintained, contrary to large spheroids. This divergent cell response may be due to the high level of senescence-related cytokines, cell stress, and damage in the inner region of large spheroids. This hypothesis was confirmed by GO analysis, in which BPs related to negative regulation of apoptosis and proliferation, such as the MAPK pathway and protein kinase B signaling, were identified only in small spheroid-derived secretome.

Moreover, cytokine arrays performed on large spheroid secretome also revealed that large spheroid secretome is enriched in molecules showing pro-angiogenic, anti-angiogenic, and pro-inflammatory activity. More specifically, the identified cytokines included IL-8, known for its role in the process of tube formation in angiogenesis (Li et al., 2003) and possessing strong pro-inflammatory activity; MMP-9, a protease involved in ECM degradation (Almalki and Agrawal, 2016); SDF-1α, capable of stimulating HUVEC migration and tube formation (M. Zhang et al., 2017); and VEGF-A. On the other hand, a subset of anti-angiogenic cytokines was identified, including MCP-1, capable of inducing endothelial cell apoptosis (X. Zhang et al., 2011); PTX-3, an inflammation-related long pentraxin that blocks the angiogenic process by inhibiting FGFR activation (O’Neill et al., 2016); and THBS-1, a protein involved in the regulation of inflammatory process, also known as an endogenous inhibitor of angiogenesis, affecting endothelial cell proliferation, migration, and stimulating apoptosis (Lawler and Jack, 2012; Feng et al., 2015). GO analysis showed that the BP “inflammatory response” is activated by some cytokines with anti-angiogenic activity and is also able to induce stress in spheroid-cultured MSCs and stimulate the onset of cellular senescence.

On the contrary, small spheroid-derived secretome showed an upregulation of potent pro-angiogenic molecules, such as ANG, IL-6, able to stimulate HUVEC proliferation and migration (Fan et al., 2008), VEGF-A, and HGF, a cytokine with pro-angiogenic and anti-inflammatory properties (Rong et al., 2018). As also observed in GO analysis, terms related to “positive regulation of angiogenesis” as well as “negative regulation of apoptosis” (activated by molecules with strong pro-angiogenic activity such as HGF, VEGF-A, and IL-6) were found, suggesting the involvement of these cytokines in protecting endothelial cells from apoptosis. In addition, this culture condition also evidenced a downregulation of all the anti-angiogenic and pro-inflammatory cytokines (IL-8, GDF-15, MCP-1, MMP-9, PTX-3, and THBS-1). Due to the low number of cytokines evaluated by this method, further analysis (such as miRNA content evaluation and proteomic studies) should be performed to better understand the influence of a different culture system on secretome composition.

EVs derived from different culture conditions were also isolated and characterized. Surprisingly, culture conditions also affected the release of EVs, a cell product extremely important in the creation of cell-free therapeutic approaches. In this study, we demonstrated that MSC spheroid culture stimulates a higher release of sEVs compared to the 2D culture, showing the double advantage of delaying cellular senescence onset and easily obtaining a higher amount of EVs. sEVs maintain classical phenotypical characteristics, particularly in terms of tetraspanin expression, which are necessary molecules for sEV biogenesis (Raposo and Stoorvogel, 2013).

Studies on MSC-derived EVs are particularly focused on the creation of therapeutic tools for the treatment of inflammatory-related disorders, such as osteoarthritis and autoimmune diseases (Palamà et al., 2020). However, the use of MSC-derived EVs in the treatment of vasculogenesis-related disorders is not so exploited. In this context, the biological effect of different secretome and EV fractions derived from different culture conditions was investigated using an in vitro HUVEC angiogenic assay. According to secretome analysis, large spheroids-derived soluble factors (CM and DM) showed a stronger effect on HUVEC migration and tube formation compared to the 2D control. However, the potential was limited compared to the soluble factors (CM and DM) derived from small spheroids, which also exhibited a strong effect on HUVEC proliferation induction. The weaker pro-angiogenic effect of large spheroid-derived soluble factors could be explained by the presence of strong anti-angiogenic factors, highlighted by secretome analyses. These inhibitory factors were indeed downregulated in small-sized spheroids, allowing them to exert a higher angiogenic potential.

We also observed that small spheroid-derived sEVs possess an enhanced angiogenic potential capable of stimulating all the key parts of angiogenesis, while no significant effects were observed for large spheroid-derived sEVs. Regarding the different effects of secretome fractions on HUVECs, we suggest a synergic action of soluble factors and EVs derived from 3D culture in the stimulation of the different phases of the angiogenic process (Gorgun et al., 2021). This hypothesized synergism was confirmed by the evidence that the fraction enriched in only soluble factors (DM) and EV population alone have a minor effect compared to the complete conditioned medium.

Considering all the data together, spheroids represent a possible alternative culture system that may be used for the maintenance of MSCs. This study demonstrated the relevance of the size of the generated spheroids to avoid the generation of a hypoxic core that will lead to an increase in cellular senescence and an alteration in the MSC secretory profile. In particular, MSCs present in small spheroids maintain the classical phenotype and show reduced cellular senescence. Moreover, MSC spheroid secretome is enriched in pro-angiogenic cytokines, suggesting a potential use of small spheroid-derived soluble factors and sEVs as a cell-free therapeutic tool for the treatment of vasculogenesis-related disorders (myocardial infarction or stroke), in which the creation of new vessels after the traumatic event is needed (Rani et al., 2015; Todorova et al., 2017; Ferreira et al., 2018; Praveen Kumar et al., 2019; Rezaie et al., 2019).

Overall, MSC spheroids are known to display increased paracrine secretion, proliferation, stemness, and anti-inflammatory and pro-angiogenic properties compared to the 2D condition. Furthermore, several studies indicate that MSC spheroid bioprinted in a biomaterial significantly supports MSC functionality and improves MSC engraftment in damaged tissues, increasing MSC survival in vivo. On this basis, the development of a reproducible large-scale manufacturing protocol to produce MSC spheroids would promote their use in clinical and therapeutic applications (Murphy et al., 2017; Allen et al., 2019; Almeida et al., 2020; Yuan et al., 2022).

Data availability statement

All relevant data supporting the findings of this study are available within the paper, its supplementary information and from the corresponding authors upon reasonable request.

Ethics statement

The studies involving humans were approved by the CER Liguria: 372/2019 Comitato etico regionale della Liguria. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

MR: conceptualization, data curation, formal analysis, investigation, methodology, resources, visualization, and writing–original draft. DR: investigation, methodology, and writing–original draft. PA: Methodology, Investigation, writing-review and editing. MP: conceptualization, data curation, supervision, validation, and writing–original draft. CG: conceptualization, funding acquisition, project administration, supervision, and writing–review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors thank the University of Genoa for funding the acquisition of the HITACHI 120 kV TEM microscope HT7800 (Grant DR3404, 2018, Heavy Equipment).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2023.1297644/full#supplementary-material

References

Acosta, J. C., Banito, A., Wuestefeld, T., Georgilis, A., Janich, P., Morton, J. P., et al. (2013). A complex secretory program orchestrated by the inflammasome controls paracrine senescence. Nat. Cell Biol. 15 (8), 978–990. doi:10.1038/NCB2784

PubMed Abstract | CrossRef Full Text | Google Scholar

Adair, T. H., and Montani, J.-P. (2010). Angiogenesis. Available at: https://www.ncbi.nlm.nih.gov/books/NBK53242/.

Google Scholar

Allen, L. M., Matyas, J., Ungrin, M., Hart, D. A., Sen, A., and Rodríguez-Lozano, F. J. (2019). Serum-free culture of human mesenchymal stem cell aggregates in suspension bioreactors for tissue engineering applications. Available at: https://ClinicalTrials.gov.

CrossRef Full Text | Google Scholar

Almalki, S. G., and Agrawal, D. K. (2016). Effects of matrix metalloproteinases on the fate of mesenchymal stem cells. Stem Cell Res. Ther. 7 (1), 129–212. doi:10.1186/s13287-016-0393-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Almeida, F., Bernardes, N., Oliveira, F. D., Costa, A. C., Fernandes-Platzgummer, A., Farinha, J. P., et al. (2020). Scalable production of human mesenchymal stromal cell-derived extracellular vesicles under serum-/xeno-free conditions in a microcarrier-based bioreactor culture system. Front. Cell Dev. Biol. 8, 553444. doi:10.3389/fcell.2020.553444

PubMed Abstract | CrossRef Full Text | Google Scholar

Alvarez-Pérez, J., Ballesteros, P., and Cerdán, S. (2005). Microscopic images of intraspheroidal PH by 1H magnetic resonance chemical shift imaging of PH sensitive indicators. Magma (New York, N.Y.) 18 (6), 293–301. doi:10.1007/S10334-005-0013-Z

PubMed Abstract | CrossRef Full Text | Google Scholar

Banfi, A., Muraglia, A., Dozin, B., Mastrogiacomo, M., Cancedda, R., and Quarto, R. (2000). Proliferation kinetics and differentiation potential of ex vivo expanded human bone marrow stromal cells: implications for their use in cell therapy. Exp. Hematol. 28 (6), 707–715. doi:10.1016/S0301-472X(00)00160-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartosh, T. J., Ylöstalo, J. H., Mohammadipoor, A., Bazhanov, N., Coble, K., Claypool, K., et al. (2010). Aggregation of human mesenchymal stromal cells (MSCs) into 3D spheroids enhances their antiinflammatory properties. Proc. Natl. Acad. Sci. U. S. A. 107 (31), 13724–13729. doi:10.1073/PNAS.1008117107

PubMed Abstract | CrossRef Full Text | Google Scholar

Basisty, N., Kale, A., Jeon, O. H., Kuehnemann, C., Payne, T., Rao, C., et al. (2020). A proteomic atlas of senescence-associated secretomes for aging biomarker development. PLoS Biol. 18 (1), e3000599. doi:10.1371/JOURNAL.PBIO.3000599

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhang, S. H., Cho, S. W., La, W. G., Lee, T. J., Yang, H. S., Sun, A. Y., et al. (2011). Angiogenesis in ischemic tissue produced by spheroid grafting of human adipose-derived stromal cells. Biomaterials 32 (11), 2734–2747. doi:10.1016/J.BIOMATERIALS.2010.12.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Birch, J., and Gil, J. (2020). Senescence and the SASP: many therapeutic avenues. Genes and Dev. 34 (23–24), 1565–1576. doi:10.1101/GAD.343129.120

PubMed Abstract | CrossRef Full Text | Google Scholar

Burlacu, A., Grigorescu, G., Rosca, A. M., Preda, M. B., and Simionescu, M. (2013). Factors secreted by mesenchymal stem cells and endothelial progenitor cells have complementary effects on angiogenesis in vitro. Stem Cells Dev. 22 (4), 643–653. doi:10.1089/SCD.2012.0273

PubMed Abstract | CrossRef Full Text | Google Scholar

Carpentier, G., Berndt, S., Ferratge, S., Rasband, W., Cuendet, M., Uzan, G., et al. (2020). Angiogenesis analyzer for ImageJ — a comparative morphometric analysis of ‘endothelial tube formation assay’ and ‘fibrin bead assay. Sci. Rep. 10 (1), 11568. doi:10.1038/S41598-020-67289-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Cesarz, Z., and Tamama, K. (2016). Spheroid culture of mesenchymal stem cells. Stem Cells Int. 2016, 1–11. doi:10.1155/2016/9176357

CrossRef Full Text | Google Scholar

Cheema, U., Brown, R. A., Alp, B., and MacRobert, A. J. (2008). Spatially defined oxygen gradients and vascular endothelial growth factor expression in an engineered 3D cell model. Cell. Mol. Life Sci. 65 (1), 177–186. doi:10.1007/s00018-007-7356-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L. C., Wang, H. W., and Huang, C. C. (2021). Modulation of inherent niches in 3d multicellular msc spheroids reconfigures metabolism and enhances therapeutic potential. Cells 10 (10), 2747. doi:10.3390/cells10102747

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, N.-C., Chen, S. Y., Li, J. R., and Young, T. H. (2013). Short-term spheroid formation enhances the regenerative capacity of adipose-derived stem cells by promoting stemness, angiogenesis, and chemotaxis. Stem Cells Transl. Med. 2, 584–594. doi:10.5966/sctm.2013-0007

PubMed Abstract | CrossRef Full Text | Google Scholar

Coppé, J. P., Desprez, P. Y., Krtolica, A., and Campisi, J. (2010). The senescence-associated secretory phenotype: the dark side of tumor suppression. Annu. Rev. Pathology 5, 99–118. doi:10.1146/ANNUREV-PATHOL-121808-102144

CrossRef Full Text | Google Scholar

Costa, E. C., Moreira, A. F., de Melo-Diogo, D., Gaspar, V. M., Carvalho, M. P., and Correia, I. J. (2016). 3D tumor spheroids: an overview on the tools and techniques used for their analysis. Biotechnol. Adv. 34 (8), 1427–1441. doi:10.1016/J.BIOTECHADV.2016.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, M. H. G., McDevitt, T. C., Cabral, J. M., da Silva, C. L., and Ferreira, F. C. (2017). Tridimensional configurations of human mesenchymal stem/stromal cells to enhance cell paracrine potential towards wound healing processes. J. Biotechnol. 262, 28–39. doi:10.1016/J.JBIOTEC.2017.09.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Curcio, E., Salerno, S., Barbieri, G., De Bartolo, L., Drioli, E., and Bader, A. (2007). Mass transfer and metabolic reactions in hepatocyte spheroids cultured in rotating wall gas-permeable membrane system. Biomaterials 28 (36), 5487–5497. doi:10.1016/J.BIOMATERIALS.2007.08.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Eletto, D., Dersh, D., and Argon, Y. (2010). GRP94 in ER quality control and stress responses. Seminars Cell Dev. Biol. 21, 479–485. doi:10.1016/j.semcdb.2010.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Evangelou, K., and Gorgoulis, V. G. (2017). Sudan black B, the specific histochemical stain for lipofuscin: a novel method to detect senescent cells. Methods Mol. Biol. 1534, 111–119. doi:10.1007/978-1-4939-6670-7_10

PubMed Abstract | CrossRef Full Text | Google Scholar

Ezquerra, S., Zuleta, A., Arancibia, R., Estay, J., Aulestia, F., and Carrion, F. (2021). Functional properties of human-derived mesenchymal stem cell spheroids: a meta-analysis and systematic review. Stem Cells Int. 2021, 1–12. doi:10.1155/2021/8825332

CrossRef Full Text | Google Scholar

Fan, Y., Ye, J., Shen, F., Zhu, Y., Yeghiazarians, Y., Zhu, W., et al. (2008). Interleukin-6 stimulates circulating blood-derived endothelial progenitor cell angiogenesis in vitro. J. Cereb. Blood Flow Metabolism 28 (1), 90–98. doi:10.1038/SJ.JCBFM.9600509

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, N., Wang, Z., Zhang, Z., He, X., Wang, C., and Zhang, L. (2015). MiR-487b promotes human umbilical vein endothelial cell proliferation, migration, invasion and tube formation through regulating THBS1. Neurosci. Lett. 591 (March), 1–7. doi:10.1016/j.neulet.2015.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandez-Rebollo, E., Franzen, J., Goetzke, R., Hollmann, J., Ostrowska, A., Oliverio, M., et al. (2020). Senescence-associated metabolomic phenotype in primary and IPSC-derived mesenchymal stromal cells. Stem Cell Rep. 14 (2), 201–209. doi:10.1016/j.stemcr.2019.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, J. R., Teixeira, G. Q., Santos, S. G., Barbosa, M. A., Almeida-Porada, G., and Gonçalves, R. M. (2018). Mesenchymal stromal cell secretome: influencing therapeutic potential by cellular pre-conditioning. Front. Immunol. 9, 2837. doi:10.3389/FIMMU.2018.02837

PubMed Abstract | CrossRef Full Text | Google Scholar

Fierabracci, A., Del Fattore, A., Muraca, M., Vittorio Delfino, D., and Muraca, M. (2016). The use of mesenchymal stem cells for the treatment of autoimmunity: from animals models to human disease. Curr. Drug Targets 17 (2), 229–238. doi:10.2174/1389450116666150722140633

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, A., Laberge, R. M., Demaria, M., and Campisi, J. (2012). Lamin B1 loss is a senescence-associated biomarker. Mol. Biol. Cell 23 (11), 2066–2075. doi:10.1091/mbc.E11-10-0884

PubMed Abstract | CrossRef Full Text | Google Scholar

Gnani, D., Crippa, S., della Volpe, L., Rossella, V., Conti, A., Lettera, E., et al. (2019). An early-senescence state in aged mesenchymal stromal cells contributes to hematopoietic stem and progenitor cell clonogenic impairment through the activation of a pro-inflammatory program. Aging Cell 18 (3), 12933. doi:10.1111/ACEL.12933

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Gualda, E., Baker, A. G., Fruk, L., and Muñoz-Espín, D. (2021). A guide to assessing cellular senescence in vitro and in vivo. FEBS J. 288 (1), 56–80. doi:10.1111/FEBS.15570

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorgun, C., Ceresa, D., Lesage, R., Villa, F., Reverberi, D., Balbi, C., et al. (2021). Dissecting the effects of preconditioning with inflammatory cytokines and hypoxia on the angiogenic potential of mesenchymal stromal cell (MSC)-Derived soluble proteins and extracellular vesicles (EVs). Biomaterials 269 (February), 120633. doi:10.1016/j.biomaterials.2020.120633

PubMed Abstract | CrossRef Full Text | Google Scholar

Grayson, W. L., Zhao, F., and Bunnell, B. (2007). Hypoxia enhances proliferation and tissue formation of human mesenchymal stem cells. Biochem. Biophysical Res. Commun. 358 (3), 948–953. doi:10.1016/J.BBRC.2007.05.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, Z., Gu, L., Eils, R., Schlesner, M., and Brors, B. (2014). Circlize implements and enhances circular visualization in R. Bioinforma. Oxf. Engl. 30 (19), 2811–2812. doi:10.1093/BIOINFORMATICS/BTU393

PubMed Abstract | CrossRef Full Text | Google Scholar

Guzmán, C., Bagga, M., Kaur, A., Westermarck, J., and Abankwa, D. (2014). ColonyArea: an ImageJ plugin to automatically quantify colony formation in clonogenic assays. PLoS ONE 9 (3), e92444. doi:10.1371/journal.pone.0092444

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez-Segura, A., Nehme, J., and Demaria, M. (2018). Hallmarks of cellular senescence. Trends Cell Biol. 28 (6), 436–453. doi:10.1016/J.TCB.2018.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Jauković, A., Abadjieva, D., Trivanović, D., Stoyanova, E., Kostadinova, M., Pashova, S., et al. (2020). Specificity of 3D MSC spheroids microenvironment: impact on MSC behavior and properties. Stem Cell Rev. Rep. 16 (5), 853–875. doi:10.1007/s12015-020-10006-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeske, R., Yuan, X., Fu, Q., Bunnell, B. A., Logan, T. M., and Li, Y. (2021). In vitro culture expansion shifts the immune phenotype of human adipose-derived mesenchymal stem cells. Front. Immunol. 12 (March), 621744. doi:10.3389/fimmu.2021.621744

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaminska, A., Wedzinska, A., Kot, M., and Sarnowska, A. (2021). Effect of long-term 3D spheroid culture on WJ-MSC. Cells 10 (4), 719. doi:10.3390/CELLS10040719

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiyoshima, T., Enoki, N., Kobayashi, I., Sakai, T., Nagata, K., Wada, H., et al. (2012). Oxidative stress caused by a low concentration of hydrogen peroxide induces senescence-like changes in mouse gingival fibroblasts. Int. J. Mol. Med. 30 (5), 1007–1012. doi:10.3892/ijmm.2012.1102

PubMed Abstract | CrossRef Full Text | Google Scholar

Kouroupis, D., and Correa, D. (2021). Increased mesenchymal stem cell functionalization in three-dimensional manufacturing settings for enhanced therapeutic applications. Front. Bioeng. Biotechnol. 9 (February), 621748. doi:10.3389/FBIOE.2021.621748

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, R., and Jat, P. (2021). Mechanisms of cellular senescence: cell cycle arrest and senescence associated secretory phenotype. Front. Cell Dev. Biol. 9, 645593. doi:10.3389/fcell.2021.645593

PubMed Abstract | CrossRef Full Text | Google Scholar

Kyurkchiev, D., Bochev, I., Ivanova-Todorova, E., Mourdjeva, M., Oreshkova, T., Belemezova, K., et al. (2014). Secretion of immunoregulatory cytokines by mesenchymal stem cells. World J. Stem Cells 6 (5), 552. doi:10.4252/WJSC.V6.I5.552

PubMed Abstract | CrossRef Full Text | Google Scholar

Langan, L. M., Dodd, N. J. F., Owen, S. F., Purcell, W. M., Jackson, S. K., and Jha, A. N. (2016). Direct measurements of oxygen gradients in spheroid culture system using electron parametric resonance oximetry. PLoS ONE 11 (2), e0149492. doi:10.1371/JOURNAL.PONE.0149492

PubMed Abstract | CrossRef Full Text | Google Scholar

Langhans, S. A. (2018). Three-dimensional in vitro cell culture models in drug discovery and drug repositioning. Front. Pharmacol. 9, 6. doi:10.3389/FPHAR.2018.00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawler, P. R., and Jack, L. (2012). Molecular basis for the regulation of angiogenesis by thrombospondin-1 and -2. Cold Spring Harb. Perspect. Med. 2 (5), a006627. doi:10.1101/cshperspect.a006627

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, A., Dubey, S., Varney, M. L., Dave, B. J., and Singh, R. K. (2003). IL-8 directly enhanced endothelial cell survival, proliferation, and matrix metalloproteinases production and regulated angiogenesis. J. Immunol. 170 (6), 3369–3376. doi:10.4049/JIMMUNOL.170.6.3369

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Hao, H., Xia, L., Ti, D., Huang, H., Dong, L., et al. (2015). Hypoxia pretreatment of bone marrow mesenchymal stem cells facilitates angiogenesis by improving the function of endothelial cells in diabetic rats with lower ischemia. PloS One 10 (5), e0126715. doi:10.1371/JOURNAL.PONE.0126715

PubMed Abstract | CrossRef Full Text | Google Scholar

Lunyak, V. V., Amaro-Ortiz, A., and Gaur, M. (2017). Mesenchymal stem cells secretory responses: senescence messaging secretome and immunomodulation perspective. Front. Genet. 8 (DEC), 220. doi:10.3389/FGENE.2017.00220

PubMed Abstract | CrossRef Full Text | Google Scholar

Markov, A., Thangavelu, L., Aravindhan, S., Zekiy, A. O., Jarahian, M., Chartrand, M. S., et al. (2021). Mesenchymal stem/stromal cells as a valuable source for the treatment of immune-mediated disorders. Stem Cell Res. Ther. 12 (1), 192. doi:10.1186/S13287-021-02265-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Marques, J. C. M. T., Cornélio, D. A., and Nogueira Silbiger, V. (2017). Identification of new genes associated to senescent and tumorigenic phenotypes in mesenchymal stem cells. Sci. Rep. 7, 1–13. doi:10.1038/s41598-017-16224-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Marrazzo, P., Crupi, A. N., Alviano, F., Teodori, L., and Bonsi, L. (2019). Exploring the roles of MSCs in infections: focus on bacterial diseases. J. Mol. Med. (Berlin, Ger. 97 (4), 437–450. doi:10.1007/S00109-019-01752-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Marrazzo, P., Pizzuti, V., Zia, S., Sargenti, A., Gazzola, D., Roda, B., et al. (2021). Microfluidic tools for enhanced characterization of therapeutic stem cells and prediction of their potential antimicrobial secretome. Antibiotics 10 (7), 750. doi:10.3390/ANTIBIOTICS10070750

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, X., Xue, M., Xu, P., Hu, F., Sun, B., and Xiao, Z. (2017). MicroRNA profiling analysis revealed different cellular senescence mechanisms in human mesenchymal stem cells derived from different origin. Genomics 109 (3–4), 147–157. doi:10.1016/J.YGENO.2017.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mueller-Klieser, W. (1984). Method for the determination of oxygen consumption rates and diffusion coefficients in multicellular spheroids. Biophysical J. 46 (3), 343–348. doi:10.1016/S0006-3495(84)84030-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, K. C., Whitehead, J., Zhou, D., Ho, S. S., and Leach, J. K. (2017). Engineering fibrin hydrogels to promote the wound healing potential of mesenchymal stem cell spheroids. Acta Biomater. 64 (December), 176–186. doi:10.1016/J.ACTBIO.2017.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Neurohr, G. E., Terry, R. L., Lengefeld, J., Bonney, M., Brittingham, G. P., Moretto, F., et al. (2019). Excessive cell growth causes cytoplasm dilution and contributes to senescence. Cell 176 (5), 1083–1097.e18. doi:10.1016/J.CELL.2019.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Noronha, N., Mizukami, A., Caliári-Oliveira, C., Cominal, J. G., Rocha, J. L. M., Covas, D. T., et al. (2019). Priming approaches to improve the efficacy of mesenchymal stromal cell-based therapies. Stem Cell Res. Ther. 10, 131. doi:10.1186/s13287-019-1224-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Noronha, N. C., Mizukami, A., Orellana, M. D., Oliveira, M. C., Covas, D. T., Swiech, K., et al. (2021). Hypoxia priming improves in vitro angiogenic properties of umbilical cord derived-mesenchymal stromal cells expanded in stirred-tank bioreactor. Biochem. Eng. J. 168 (April), 107949. doi:10.1016/j.bej.2021.107949

CrossRef Full Text | Google Scholar

Ohori-Morita, Y., Niibe, K., Limraksasin, P., Nattasit, P., Miao, X., Yamada, M., et al. (2022). Novel mesenchymal stem cell spheroids with enhanced stem cell characteristics and bone regeneration ability. Stem Cells Transl. Med. 11 (4), 434–449. doi:10.1093/STCLTM/SZAB030

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Neill, C. L., Guduric-Fuchs, J., Chambers, S. E. J., O’Doherty, M., Bottazzi, B., Stitt, A. W., et al. (2016). Endothelial cell-derived pentraxin 3 limits the vasoreparative therapeutic potential of circulating angiogenic cells. Cardiovasc. Res. 112 (3), 677–688. doi:10.1093/cvr/cvw209

PubMed Abstract | CrossRef Full Text | Google Scholar

Palamà, M. E. F., Shaw, G. M., Carluccio, S., Reverberi, D., Sercia, L., Persano, L., et al. (2020). The secretome derived from mesenchymal stromal cells cultured in a xeno-free medium promotes human cartilage recovery in vitro. Front. Bioeng. Biotechnol. 8, 90. doi:10.3389/fbioe.2020.00090

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinto, D. S., Ahsan, T., Serra, J., Fernandes-Platzgummer, A., Cabral, J. M. S., and da Silva, C. L. (2020). Modulation of the in vitro angiogenic potential of human mesenchymal stromal cells from different tissue sources. J. Cell. Physiology 235 (10), 7224–7238. doi:10.1002/jcp.29622

PubMed Abstract | CrossRef Full Text | Google Scholar

Pizzuti, V., Paris, F., Marrazzo, P., Bonsi, L., and Alviano, F. (2023). Mitigating oxidative stress in perinatal cells: a critical step toward an optimal therapeutic use in regenerative medicine. Biomolecules 13, 971. doi:10.3390/BIOM13060971

PubMed Abstract | CrossRef Full Text | Google Scholar

Praveen Kumar, L., Kandoi, S., Misra, R., S., V., K., R., and Verma, R. S. (2019). The mesenchymal stem cell secretome: a new paradigm towards cell-free therapeutic mode in regenerative medicine. Cytokine and Growth Factor Rev. 46, 1–9. doi:10.1016/J.CYTOGFR.2019.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Rani, S., Ryan, A. E., Griffin, M. D., and Ritter, T. (2015). Mesenchymal stem cell-derived extracellular vesicles: toward cell-free therapeutic applications. Mol. Ther. J. Am. Soc. Gene Ther. 23 (5), 812–823. doi:10.1038/MT.2015.44

CrossRef Full Text | Google Scholar

Raposo, G., and Stoorvogel, W. (2013). Extracellular vesicles: exosomes, microvesicles, and friends. J. Cell Biol. 200, 373–383. doi:10.1083/jcb.201211138

PubMed Abstract | CrossRef Full Text | Google Scholar

Redondo-Castro, E., Cunningham, C. J., Miller, J., Brown, H., Allan, S. M., and Pinteaux, E. (2018). Changes in the secretome of tri-dimensional spheroid-cultured human mesenchymal stem cells in vitro by interleukin-1 priming. Stem Cell Res. Ther. 9 (1), 11. doi:10.1186/S13287-017-0753-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Regmi, S., Kumar Raut, P., Pathak, S., Shrestha, P., Park, P. H., and Jeong, J. H. (2021). Enhanced viability and function of mesenchymal stromal cell spheroids is mediated via autophagy induction. Autophagy 17 (10), 2991–3010. doi:10.1080/15548627.2020.1850608

PubMed Abstract | CrossRef Full Text | Google Scholar

Rezaie, J., Heidarzadeh, M., Hassanpour, M., Hassan, A., Shokrollahi, E., Ahmadi, M., et al. (2019). The angiogenic paracrine potential of mesenchymal stem cells. Update Mesenchymal Induc. Pluripotent Stem Cells 2019, 84433. doi:10.5772/INTECHOPEN.84433

CrossRef Full Text | Google Scholar

Rong, S., Wang, X. l., Wang, Y. c., Wu, H., Zhou, X. d., Wang, Z. k., et al. (2018). Anti-inflammatory activities of hepatocyte growth factor in post-ischemic heart failure. Acta Pharmacol. Sin. 39 (10), 1613–1621. doi:10.1038/APS.2018.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Roubelakis, M. G., Tsaknakis, G., Pappa, K. I., Anagnou, N. P., and Watt, S. M. (2013). Spindle shaped human mesenchymal stem/stromal cells from amniotic fluid promote neovascularization. PloS One 8 (1), e54747. doi:10.1371/JOURNAL.PONE.0054747

PubMed Abstract | CrossRef Full Text | Google Scholar

Roura, S., Juli, R., Soler-Botija, C., Pujal, J. M., Gálvez-Montón, C., Prat-Vidal, C., et al. (2012). Human umbilical cord blood-derived mesenchymal stem cells promote vascular growth in vivo. PloS One 7 (11), e49447. doi:10.1371/JOURNAL.PONE.0049447

PubMed Abstract | CrossRef Full Text | Google Scholar

Sart, S., Tsai, A. C., Li, Y., and Ma, T. (2014). Three-dimensional aggregates of mesenchymal stem cells: cellular mechanisms, biological properties, and applications. Tissue Eng. Part B, Rev. 20 (5), 365–380. doi:10.1089/TEN.TEB.2013.0537

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitt, C. A., Wang, B., and Demaria, M. (2022). Senescence and cancer - role and therapeutic opportunities. Nat. Rev. Clin. Oncol. 19 (10), 619–636. doi:10.1038/S41571-022-00668-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Soares, C. P., Midlej, V., Oliveira, M. E. W. d., Benchimol, M., Costa, M. L., and Mermelstein, C. (2012). 2D and 3D-organized cardiac cells shows differences in cellular morphology, adhesion junctions, presence of myofibrils and protein expression. PloS One 7 (5), e38147. doi:10.1371/JOURNAL.PONE.0038147

PubMed Abstract | CrossRef Full Text | Google Scholar

Todorova, D., Simoncini, S., Lacroix, R., Sabatier, F., and Dignat-George, F. (2017). Extracellular vesicles in angiogenesis. Circulation Res. 120, 1658–1673. doi:10.1161/CIRCRESAHA.117.309681

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, A. C., Liu, Y., Yuan, X., and Ma, T. (2015). Compaction, fusion, and functional activation of three-dimensional human mesenchymal stem cell aggregate. Tissue Eng. Part A 21 (9–10), 1705–1719. doi:10.1089/TEN.TEA.2014.0314

PubMed Abstract | CrossRef Full Text | Google Scholar

Turinetto, V., Vitale, E., and Giachino, C. (2016). Senescence in human mesenchymal stem cells: functional changes and implications in stem cell-based therapy. Int. J. Mol. Sci. 17 (7), 1164. doi:10.3390/IJMS17071164

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, W., Horn, P., Castoldi, M., Diehlmann, A., Bork, S., Saffrich, R., et al. (2008). Replicative senescence of mesenchymal stem cells: a continuous and organized process. PLOS ONE 3 (5), e2213. doi:10.1371/JOURNAL.PONE.0002213

PubMed Abstract | CrossRef Full Text | Google Scholar

Watt, S. M., Gullo, F., van der Garde, M., Markeson, D., Camicia, R., Khoo, C. P., et al. (2013). The angiogenic properties of mesenchymal stem/stromal cells and their therapeutic potential. Br. Med. Bull. 108 (1), 25–53. doi:10.1093/BMB/LDT031

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, W., Zhang, F., Xie, C., Jiang, M., and Hou, M. (2015). Macrophage migration inhibitory factor confers resistance to senescence through CD74-dependent AMPK-FOXO3a signaling in mesenchymal stem cells. Stem Cell Res. Ther. 6 (1), 82. doi:10.1186/S13287-015-0076-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, A. W., Zacharias, N. A., Binder, B. Y. K., and Murphy, W. L. (2021). Controlled aggregation enhances immunomodulatory potential of mesenchymal stromal cell aggregates. Stem Cells Transl. Med. 10 (8), 1184–1201. doi:10.1002/SCTM.19-0414

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, X., Liu, Y., Bijonowski, B. M., Tsai, A. C., Fu, Q., Logan, T. M., et al. (2020). NAD+/NADH redox alterations reconfigure metabolism and rejuvenate senescent human mesenchymal stem cells in vitro. Commun. Biol. 3 (1), 774–815. doi:10.1038/s42003-020-01514-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, X., Sun, L., Jeske, R., Nkosi, D., York, S. B., Liu, Y., et al. (2022). Engineering extracellular vesicles by three-dimensional dynamic culture of human mesenchymal stem cells. J. Extracell. Vesicles 11 (6), e12235. doi:10.1002/JEV2.12235

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, M., Qiu, L., Zhang, Y., Xu, D., Zheng, J. C., and Jiang, L. (2017). CXCL12 enhances angiogenesis through CXCR7 activation in human umbilical vein endothelial cells. Sci. Rep. 7 (1), 8289. doi:10.1038/s41598-017-08840-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Liu, X., Shang, H., Xu, Y., and Qian, M. (2011). Monocyte chemoattractant protein-1 induces endothelial cell apoptosis <italic&amp;gt;in vitro&amp;lt;/italic&amp;gt; through a p53-dependent mitochondrial pathway. Acta Biochimica Biophysica Sinica 43 (10), 787–795. doi:10.1093/abbs/gmr072

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Zhu, W., He, H., Fan, B., Deng, R., Hong, Y., et al. (2019). Macrophage migration inhibitory factor rejuvenates aged human mesenchymal stem cells and improves myocardial repair. Aging 11 (24), 12641–12660. doi:10.18632/AGING.102592

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, J., Lei, T., Guo, P., Yu, J., Xu, Q., Luo, Y., et al. (2019). Mechanisms shaping the role of ERK1/2 in cellular senescence (review). Mol. Med. Rep. 19 (2), 759–770. doi:10.3892/mmr.2018.9712

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: mesenchymal stromal cells, spheroids, senescence, secretome, extracellular vesicles, angiogenesis

Citation: Rovere M, Reverberi D, Arnaldi P, Palamà MEF and Gentili C (2023) Spheroid size influences cellular senescence and angiogenic potential of mesenchymal stromal cell-derived soluble factors and extracellular vesicles. Front. Bioeng. Biotechnol. 11:1297644. doi: 10.3389/fbioe.2023.1297644

Received: 20 September 2023; Accepted: 23 November 2023;
Published: 12 December 2023.

Edited by:

Bruce Alan Bunnell, University of North Texas Health Science Center, United States

Reviewed by:

Sara Al-Ghadban, University of North Texas Health Science Center, United States
Pasquale Marrazzo, University of Bologna, Italy
Xuegang Yuan, University of California, Los Angeles, United States
Carl Austin Gregory, Texas A and M Health Science Center, United States

Copyright © 2023 Rovere, Reverberi, Arnaldi, Palamà and Gentili. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chiara Gentili, Y2hpYXJhLmdlbnRpbGlAdW5pZ2UuaXQ=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.