A 3D ovarian cancer metastasis model using a decellularized peritoneal matrix to study therapy response

2.1 Materials

2.2 Generation of SIS from porcine intestines

Production of small intestinal submucosa matrix (SIS) by decellularization of porcine jejunum was performed as described previously (32). Briefly, segments of jejunum were cleaned with tap water and cut into 20 cm segments. Intestinal segments were inverted, and mucosa was scraped off using forceps, followed by incubation overnight in Ca²⁺- and Mg²⁺-free phosphate buffered saline (PBS^—) with 1% P/S. Prior to decellularization, specimens for DNA quantification, proteomic and histopathological analysis were collected. Samples for histopathological analysis were fixed in 4% paraformaldehyde (PFA) for 1 hour at room temperature and paraffin-embedded. Samples for DNA quantification and proteomics were frozen at −80 °C. For tissue decellularization, intestinal segments were first filled with 3.4% sodium deoxycholate (DOC; cat # D6750, Sigma Aldrich, St. Louis, USA) dissolved in Milli-Q (MQ) water and incubated at 4 °C for 1.5 hours with a magnetic stirrer. After washing 5 times during 1 hour with PBS^— containing 1% P/S, SIS pieces were stored in PBS^— containing 1% P/S overnight. The following day, SIS pieces were incubated with 0.167 mg/mL DNAse I (cat # 10104159001, Roche, Basel, Switzerland) in PBS with Ca²⁺ and Mg²⁺ (PBS⁺; cat # D8662, Sigma Aldrich, St. Louis, USA) for 2 hours at 37 °C under continuous agitation, then washed 3 times during 1 hour with PBS^—. Samples for histopathological analysis, DNA quantification, and proteomics were collected as described above. SIS pieces were stored in PBS^— containing 1% P/S at 4 °C. Generated SIS scaffolds were γ-sterilized with a dosage of 25 kGy (Gammatom srl, Guanzate, Italy) and stored in PBS^— with 1% P/S at 4 °C until further use.

2.3 Generation of PerMa from porcine and human peritoneum

Peritoneal matrix (PerMa) was generated by decellularization of parietal peritoneum from the anterior abdominal walls from three different sources: pigs, autopsy of women with healthy peritoneum, and women with ovarian cancer undergoing primary cytoreductive surgery. Before the decellularization, samples for DNA quantification, proteomic analysis, and histopathological analysis were collected, as described for SIS. After sampling, the peritoneum was washed for a minimum of 1 hour and up to 18 hours with PBS^— with 1% P/S. Submesothelial fat was removed by blunt dissection to decrease the thickness of the membranes and improve decellularization efficiency. For decellularization, the peritoneum was incubated in 3.4% DOC dissolved in MQ water for 4 hours at 4 °C on a tilting platform. After decellularization, the peritoneum was washed in PBS^—5 times for 1 hour at 4 °C, followed by treatment with 0.167 mg/mL DNAse I dissolved in PBS⁺ for 2 hours at 37 °C under continuous agitation. The peritoneum was then washed with PBS^—3 times for 1 hour at 4 °C on a tilting platform. Decellularized tissue samples for histopathological analysis, DNA quantification, and proteomics were collected as described above. The remaining decellularized scaffold was stored in PBS^— with 1% P/S at 4 °C. Generated PerMa scaffolds were γ-sterilized with a dosage of 25 kGy (Gammatom srl, Guanzate, Italy) and stored in PBS^— with 1% P/S at 4 °C until further use.

2.4 Histology and immunohistochemistry

Specimens from native organs, decellularized scaffolds, and scaffolds repopulated with cells were fixed with 4% PFA for 1 hour at room temperature and embedded in paraffin. Formalin-fixed paraffin-embedded (FFPE) samples were cut into 3 µm thick sections.

2.5 DNA quantification

Specimens of native organs and decellularized scaffolds were lyophilized overnight, weighed and chopped into 1×1 mm pieces, and digested with Proteinase K (Qiagen, Hilden, Germany) for 3 hours at 56 °C. DNA was extracted using the DNeasy blood and tissue kit (cat# 69504, Qiagen, Hilden, Germany) according to manufacturer’s instructions. DNA content was quantified with the Quant-iT^TM PicoGreen^TM dsDNA Assay Kit (ThermoFisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Measurements were performed at an excitation of 480 nm and an emission of 520 nm using a SpectraMax Gemini EM plate reader (Molecular Devices, San Jose, CA, USA). DNA values were determined from a Lambda-DNA standard curve.

2.6 3D printed crowns

Culture plate inserts, referred to as crowns, and staining rings, were custom-designed and 3D printed with a Form 2 3D printer (Formlabs, Somerville, MA, USA), using High Temp resin (Formlabs, Somerville, USA). The dimensions are given in Supplementary Fig. 2, and stl-files are given as Supplementary files. The post-processing included 15 minutes of soaking in isopropyl alcohol, followed by drying with compressed air to remove excess print material and curing in a UV chamber in two rounds of 20 minutes. The crowns and staining rings were autoclaved before use for cell cultures.

2.7 Multiphoton microscopy

2.8 Diffusion assay

Diffusion assay was performed as previously described (34). In short, FITC-dextran with an average molecular weight of 4 kD or 40 kD (cat # 46944 and FD40, Sigma Aldrich, St. Louis, MO, USA) was dissolved in PBS^— at a concentration of 10 µM, and low-molecular-weight fragments were removed with Amicon Ultra-4 Centrifugal Filters with a cut-off of 3 kD or 30 kD (cat# UFC8003 and UFC8030, Merck Millipore, Burlington, VT, USA). Scaffolds were mounted on 3D printed crowns to generate a Transwell-like system and incubated with PBS^— at 4 °C overnight. At the start of the experiment, 2200 µL fresh PBS^— was added to the outer compartment, and 1000 µL of 10 µM FITC-dextran solution was added to the inner chamber. The plate was placed in an incubator at 37 °C, 5% CO2 for the remainder of the experiment. At given time points (0, 5, 10, 20, 40, 60, 90, 120, 180, 240, 300, 360, 1440 minutes) 50 µL samples were taken from the outer compartment and the initial volume in the outer compartment was restored by addition of fresh PBS^—. The fluorescence signal was measured at an excitation of 490 nm and an emission of 525 nm using a SpectraMax Gemini EM plate reader (Molecular Devices, San Jose, CA, USA). The amount of FITC-dextran that diffused across the scaffold over time was determined from a FITC-dextran standard curve. The apparent permeability coefficient (Papp), giving the passive permeability of the substance per area over time, was calculated as described before (35). Briefly, Papp is calculated using the following formula: where dQ/dt is the rate of permeability (μmol/s), C0 is the initial concentration in the donor compartment, and A is the surface area of the filter (1.33 cm²).

2.9 Rheology

Rheological measurements were conducted using an Anton Paar Physica MCR302 (Anton Paar, Ostfildern, Germany) in a plate to plate set up equipped with a 25 mm diameter measuring probe, a solvent trap and a Peltier element. Samples were trimmed to fit the probe diameter and soaked with PBS^—. The gap width was determined by lowering the probe until a standard force of 1 N was reached. Measurements were performed at 37 °C. Initially, a frequency sweep from 0.1 to 100 rad/s was performed at 0.5% strain deformation, followed by an amplitude sweep with increasing strain deformation from 0.01 – 500% at an angular frequency of 10 rad/s. Mean values of storage and loss modulus as well as the damping factor were calculated from the linear region of amplitude sweep measurements. The flow point was determined as shear stress at the intersection point of storage and loss modulus (G’ = G’’).

2.10 Proteomics

2.11 SIS and PerMa cell cultures

2.12 Drug screen and viability assay

Carboplatin (cat # 507662, 10 mg/mL, Fresenius Kabi, Bad Homburg, Germany) was used within 48 hours after the seal was broken. Drug dilution series were prepared using cell culture medium. OV90 cells were used for all drug sensitivity assays.

To determine carboplatin EC50 values for the OV90 cell line cultured in a 2D monolayer, WST-1 assay was performed (cat # 11644807001, Roche, Basel, Switzerland). 2000 cells/well were seeded in 96 well cell culture plates and allowed to attach for 24 hours. The following day, the culture medium was removed, and carboplatin at different concentrations was added onto the plate. After 72 hours, WST-1 reagent was added (1:10) to each well and incubated for 2 hours at 37 °C, 5% CO2. The absorbance of formazan was measured at 440 nm with a Spectra Max Plus 384 microplate reader (Molecular Devices, San Jose, USA). The experiment was performed in triplicate and repeated three times.

Carboplatin sensitivity for OV90 in 2D culture was also assessed using apoptosis and necrosis as markers for cell viability. Five doses (10 µM, 31.6 µM, 45.6 µM (= EC50), 100 µM, and 316 µM) were chosen based on the WST-1 assay results. 90,000 cells/well were seeded in 6 well culture plates in 2 mL volume and allowed to attach for 24 hours. The following day, the culture medium was removed, and Carboplatin at different concentrations was added onto the plate at a volume of 3 mL per well. After 72 hours, the medium from each well was collected in individual tubes. The remaining attached cells were then trypsinized and combined with the medium containing detached cells. The cells were stained with Annexin V/PI (cat # V13245, Invitrogen, Waltham, MA, USA) and analyzed by flow cytometry for detection of apoptotic and necrotic cells. The experiment was performed with three technical replicates for each drug dose.

Based on the EC50 value for OV90 and the results from the apoptosis/necrosis assay, three drug doses (10 µM, 50 µM, and 100 µM) were chosen for further drug sensitivity testing in 3D culture. PerMa was mounted as a Transwell-like system and placed in 12 well plates in cell culture medium at 37 °C, 5% CO2 overnight. The following day, 25,000 cells/scaffold were seeded. On day three after seeding, the scaffolds were flipped, leaving the cells on the outer surface. On day seven after seeding, the PerMa scaffolds were imaged using a Dragonfly 505 confocal spinning disk system (Andor Technologies). Following imaging, the scaffolds were moved to a new 12-well plate containing carboplatin at different doses and cultured at 37 °C, 5% CO2. On day ten after seeding, the PerMa scaffolds were stained with an apoptosis/necrosis assay kit (cat # ab176749, Abcam, Cambridge, United Kingdom). Staining was performed according to the manufacturer’s recommendations using the 3D-printed staining ring. Briefly, the staining reagents and washing buffers were added in a volume of 300 µL to the staining ring, and the crowns with PerMa were placed on top of the staining ring. Staining was performed according to manufacturer’s recommendations. Scaffolds were washed in washing buffer for 1 minute, stained with 7-AAD, Apopxin Green, and CytoCalcein Violet for 30 minutes, followed by 3 x 1 minute washing in washing buffer. The scaffolds were imaged using the Dragonfly confocal microscope. Images were taken at 10x magnification, with a Z-stack depth of 150 to 300 µm with a 2 µm step size. Three representative imaging fields were selected for each scaffold. 7-AAD was imaged using a 561 nm laser line and a 575-625 nm filter. Apopxin Green was imaged using a 488 nm laser line and a 500-550 nm filter. CytoCalcein Violet was imaged using a 405 nm laser line and a 425-475 nm filter. Image analysis was performed in the IMARIS image analysis software using the surfaces and spots functions.

2.13 CAR T cytotoxicity testing

TAG72 CAR single chain variable fragment (scFv), huCC49, was derived from Kashmiri et al. (39). Muc16 CAR scFv (K101) was generated as described in Casey et al. (40). Both CARs are second-generation containing a CD8 hinge and CD8 transmembrane domain, a 4-1BB intracellular domain, and CD3ζ (41). Retroviral particles carrying the CAR vectors and transduction of T cells were done as described previously (42). Briefly, HEK cells were transfected with packaging vectors and a retroviral vector encoding the CAR constructs. After 2 days of harvesting, the supernatant, containing viral particles, was used in two rounds of transduction of PBMCs activated for 48 hours with anti-CD3 and anti-CD28 antibodies (coated on a 24-well plate at 1 µg/mL respectively). CAR expression on the T cells was assessed with an anti-G4S linker (Cell Signaling Technology, US) 3 and 8 days after transduction.

Cell cultures were prepared as described in section 2.11. Briefly, 25,000 OV90 cells were seeded per scaffold and treated with CAR T cells targeting TAG-72 or Mock T cells. The scaffolds were flipped and imaged with a Dragonfly 505 confocal microscope after three and seven days to confirm cell growth. On day eight, the scaffolds containing OV90 cells were flipped back so the cells were again on the inner surface and placed in new wells containing 1.5 mL medium. CAR T or mock T cells were labeled with CellTrace^TM Violet Cell Proliferation Kit (CTV; cat # C34557, Invitrogen, Waltham, USA) according to the manufactureŕs protocol and seeded in a volume of 0.5 mL to the inner chamber. The scaffolds were incubated at 37 °C, 5% CO2. After 24 hours, the scaffolds were flipped, so the cells were again on the outer surface and imaged with a Dragonfly 505 confocal microscope, as described above. CAR T/mock T cells were imaged using a 405 nm laser line and a 425-475 nm filter and GFP positive EOC cell lines using a 488 nm laser line and a 500-550 nm filter. Image analysis was performed in the IMARIS image analysis software using the surfaces function.

2.14 Statistical analyses

Statistical analyses were performed in PRISM Graphpad (version 10.1.0). Data is presented as mean ± standard deviation. For analyses comparing two datasets, unpaired t test was used. For comparison of more than two datasets, analysis of variance (ANOVA) with Tukey’s multiple comparison test was used. Significance was set to p < 0.05. Details about statistical analyses, including p values, are presented in the figure legends.

3.1 Chemical and enzymatic decellularization of peritoneum results in a cell-free ECM scaffold

3D in vitro organotypic models can be generated by culturing normal and malignant cells on decellularized tissue scaffolds, such as the porcine small intestinal submucosa (SIS) scaffold (30, 31). To establish a model that mimics peritoneal HGSOC metastases, we developed a decellularized tissue scaffold from porcine peritoneum, referred to as pPerMa. We compared this to the previously established porcine SIS (pSIS) and decellularized human peritoneum obtained from either autopsy of patients with no intra-abdominal disease (hPerMa) or surgery of EOC patients (ocPerMa; Fig. 1a).

• Download figure
• Open in new tab

Figure 1.

Decellularization of small intestine and peritoneum generates dECM scaffolds. a Schematic illustration of the generation of dECM scaffolds from porcine small intestine and peritoneum, and human peritoneum. b Photographs of porcine small intestine and porcine and human peritoneum before (left) and after (right) decellularization. Scale bars 30 mm. c Corresponding hematoxylin and eosin (H&E) staining of native and decellularized tissues. Scale bars 50 µm. d DNA content of native and decellularized porcine small intestine (n = 7 native and n = 7 decellularized), porcine peritoneum (n = 8 native and n = 8 decellularized), and human peritoneum (n = 4 native and n = 6 decellularized). Data represent mean ± SD, unpaired two-tailed t-test.

Porcine and human parietal peritoneum were collected from the anterior abdominal walls and decellularized using a chemical detergent followed by enzymatic DNA degradation (Supplementary Figure 1). We chose the detergent sodium deoxycholate (SDC) since it has been shown to effectively decellularize tissues of various thicknesses without compromising the structure of the ECM (43–45). The SIS scaffold was prepared as described previously (32).

To evaluate the decellularization efficacy, we assessed cell removal by histology and residual DNA content. Samples were taken from native intestine and peritoneum as well as from small intestinal submucosa (SIS) and peritoneal matrix (PerMa; Fig. 1b). The decellularization resulted in a cell-free scaffold as demonstrated by hematoxylin and eosin (H&E) staining (Fig. 1c). DNA quantification confirmed depletion of cell nuclei and genetic material to an average 20.03 ± 10.97 ng/mg in pSIS, 2.91 ± 4.62 ng/mg in pPerMa, and 12.48 ± 26.95 ng/mg in hPerMa (Fig. 1d). Altogether, the decellularization protocol results in an efficient removal of cells and a significant decrease in DNA content.

3.2 A custom-designed culture plate insert enables characterization and cell culture

To enable the characterization and the establishment of cell cultures with the dECM scaffolds, we engineered a 3D printed culture plate insert, referred to as a crown with attachment clip (Fig. 2a and Supplementary Fig. 2a,b). The crown serves three key functions: 1) elevating and securing the dECM scaffold above the surface of a multiwell culture plate, 2) establishing a Transwell®-like system where two chambers are separated solely by the scaffold, and 3) enabling confocal imaging, ensuring direct contact between the scaffold and the imaging surface. To enable the staining of scaffolds in a lower volume, we developed a 3D-printed staining ring (Fig. 2b, Supplementary Fig. 2c). The crowns and staining ring were developed to fit 12-well plates.

• Download figure
• Open in new tab

Figure 2.

Custom-designed 3D-printed culture plate inserts and staining rings enable cell culture, staining, and confocal imaging of 3D cultures using dECM scaffolds. a A 3D-printed culture plate insert, referred to as a crown, suspends and secures the dECM scaffold above the surface of a culture plate well, creating a Transwell®-like setup. b A staining ring enables staining of the dECM scaffold in a volume of 300 µL.

• Download figure
• Open in new tab

Figure 3.

The structure of major ECM proteins differs between the dECM scaffolds. a Representative images of Elastica v. Gieson and Alcian blue staining and immunohistochemical staining of laminin, fibronectin, collagen I, and collagen IV in pSIS, pPerMa, and hPerMa. Scale bars 50 µm. b Representative second harmonic generation (SHG) microscopy images of pSIS, pPerMa, hPerMa, and ocPerMa. 3D renderings (left) and single optical sections (right). Scale bars = 50 µm. c Schematic illustration of the setup for multiphoton imaging. The dECM scaffolds are mounted on crowns and placed bottom-up in a 6-well plate filled with PBS. The multiphoton objective is immersed in PBS. d Example images showing the masking of SHG images for The Workflow Of Matrix BioLogy Informatics (TWOMBLI) image analysis. e TWOMBLI analysis results. Data represent mean ± SD. *P < 0.05, **P < 0.01, one-way ANOVA with Tukey’s multiple comparison test. N = 3 for pSIS, n = 4 for pPerMa, n = 5 for hPerMa, n = 2 for ocPerMa. Three representative images were captured and analyzed for each sample. Images were cropped to 250 µm (x) * 250 µm (y) * 10 µm (z) stacks. Four z stacks were generated per image area to include the top 40 µm from the surface.

3.3 ThedECM scaffolds retain the structure of the native tissues and show similarities in structure

To assess the integrity of the ECM after tissue decellularization, we compared histological staining of major ECM components and the structure of the ECM with multiphoton microscopy (Fig. 3 and Supplementary Fig. 3). Histological staining showed similar patterns of laminin, fibronectin, collagen I, and collagen IV positivity in the native (Supplementary Fig. 3a) and decellularized (Fig. 3a) tissues from the same source. Second harmonic generation (SHG) and two-photon excitation fluorescence (TPEF) by multiphoton microscopy allow label-free visualization of non-centrosymmetric structures (typically collagen fibers) and tissue autofluorescence, respectively (48) (Fig. 3b). For multiphoton imaging, the scaffolds were mounted on crowns, placed in 6 well plates with the tissue facing up, and submerged in PBS (Fig. 3c). To assess the effect of decellularization on the matrix patterns, we used the ImageJ plugin TWOMBLI to quantitatively analyze the multiphoton images of our samples (33) (Fig. 3d). Fiber lacunarity, fractal dimension, % high-density matrix (HDM), branchpoints, and endpoints all remained unaltered between native and decellularized samples of the same tissue type (Fig. 3e). Only matrix alignment increased in the decellularized pPerMa (0.2931 ± 0.0531) compared to the native porcine peritoneum (0.1109 ± 0.0492), p = 0.0003 and curvature decreased in the decellularized pSIS (21.41 ± 2.011) compared to the native porcine intestine (25.44 ± 0.9989), p = 0.0358 (Supplementary Fig. 3c). In summary, the decellularization preserved the structure of major ECM components.

We also used histological staining and multiphoton microscopy to assess the differences between dECM scaffolds derived from different tissues. Alcian Blue stained positive only in pSIS and not pPerMa or hPerMa (Fig. 3a). Elastica/Van Gieson stain demonstrated elastic fibers interspersed between collagen fibers in the submesothelial stroma of human and porcine PerMa. The elastic fibers of pSIS were located deeper in the tissue. Immunohistochemical staining showed positivity for laminin, fibronectin, collagen I, and collagen IV in hPerMa (Fig. 3a). Collagen I and collagen IV stained positive in pSIS and pPerMa, while laminin stained only weakly positive in pSIS, and fibronectin did not stain in pSIS nor pPerMa. Of note, the content of high-density matrix, as assessed by TWOMBLI in the multiphoton images, was significantly lower in ocPerMa (95.82% ± 2.440) compared to hPerMa (99.47% ± 0.140, p = 0.0005) and pPerMa (98.67% ± 0.599, p = 0.0101) (Fig. 3e). Instead, the matrix had significantly more fiber alignment in pPerMa (0.2931 ± 0.0531) compared to pSIS (0.1163 ± 0.0353, p = 0.0047), hPerMa (0.1765 ± 0.0595, p = 0.0273) and ocPerMa (0.1388 ± 0.0211, p = 0.0238) (Fig. 3e). No significant differences were observed between the dECM scaffolds in lacunarity, fractal dimension, curvature, branchpoints, or endpoints (Fig. 3e). Overall, the structure of major ECM components was most comparable in the dECM scaffolds from different normal tissue sources. Of note, the ocPerMa displayed a reduced high-density matrix compared to the others, however it should be noted the small number of just two tissues.

3.4 PerMa is a viscoelastic, permeable scaffold

Next to structure, the physical properties of the ECM, such as elasticity or porosity are known to influence cellular behaviors, including spreading, growth, proliferation, migration, and differentiation (49). Therefore, we investigated whether the scaffolds differ in their physical characteristics. Oscillatory rheological measurements were conducted to assess the mechanical properties of the pPerMa in comparison to the pSIS and the hPerMa.

Frequency sweep measurements revealed a similar response of all three scaffolds to increasing angular frequencies (Fig. 4b). Storage modulus (G’), which is a measure of the elastic properties of a material, and loss modulus (G’’), describing the viscous material portion, exhibited parallel runs, confirming a viscoelastic solid behavior of all scaffolds. Shear moduli G’ and G’’ differed among the scaffolds (hPerMa < pSIS < pPerMa), indicating different viscoelastic properties. Strain sweep measurements confirmed this presumption (Fig. 4c). The pPerMa exhibited a significantly higher storage and loss modulus (846 ± 100 Pa and 122 ± 26 Pa) compared to the pSIS (464 ± 167 Pa and 60 ± 31 Pa) and the hPerMa (207 ± 156 Pa and 41 ± 27 Pa), demonstrating that the pPerMa is more rigid than hPerMa and the pSIS (Fig. 4d,e). Next, the damping factor (tan δ = G’’ / G’) was calculated to assess the elastic and viscous proportions of the scaffolds. The pPerMa (0.14 ± 0.02) and the pSIS (0.12 ± 0.02) exhibited a significantly lower damping factor compared to the hPerMa (0.23 ± 0.08) (Fig. 4f), suggesting an increased portion of elastic components in the pPerMa. This was corroborated by analysis of the flow point, defining the transition from solid to viscous behavior (G’ = G’’) of a viscoelastic material. For both, pPerMa (53 ± 14 Pa) and pSIS (46 ± 15 Pa), a significantly higher shear stress (τ) was required to reach the flow point compared to the hPerMa (18 ± 15 Pa) (Fig. 4g), indicating that both, pPerMa and pSIS contain more elastic components than the hPerMa. Summarized, the pPerMa is more rigid than the pSIS and the hPerMa while the viscoelastic behavior is more similar to the porcine SIS.

• Download figure
• Open in new tab

Figure 4.

Rheological analyses reveal differences in the stiffness and permeability of the dECM scaffolds. a Photographs of porcine SIS (blue), porcine PerMa (green), and human PerMa (yellow) mounted in the plate-to-plate rheometer. b,c Data plot of oscillatory measurements showing the storage (G’, square) and loss (G’’, circle) modulus of the pSIS, the pPerMa, and the hPerMa. Moduli are shown as a function of frequency (b, γ = 0.5%) and strain deformation (c, ω = 10 rad/s). d-g Rheological parameters of the dECM scaffolds calculated from the linear region of the amplitude sweep. Single data points represent measurements of individual biological samples. Data represent mean ± SD. *P < 0.05, **P < 0.01, ****P < 0.0001, one-way ANOVA with Tukey’s multiple comparison test. h Schematic illustration of the diffusion assay. The dECM scaffolds were mounted in a Transwell®-like setup. 4 or 40 kD FITC-Dextran was added to the inner chamber, and samples were taken from the outer chamber at different time points (0, 5, 10, 20, 40, 60, 90, 120, 180, 240, 300, 3600, and 1440 minutes). The fluorescence signal was analyzed with a plate reader. i,j Diffusion of 4 kD (i) and 40 kD (j) FITC-Dextran across the scaffolds over time. Data represent mean ± SD displayed as error shadows. The experiment was performed three times with three replicates for each scaffold. k,l The apparent permeability coefficient (Papp) for 4 kD (j) and 40 kD (k) FITC-Dextran calculated for different time frames. *P < 0.05, **P < 0.01, one-way ANOVA with Tukey’s multiple comparison test.

To evaluate scaffold permeability, we analyzed the diffusion of molecules with low (4 kDa) and high (40 kDa) molecular weights across the scaffolds. Decellularized scaffolds were mounted in a Transwell-like setup with an inner and an outer chamber separated only by the scaffolds (Fig. 4h). 4 kDa or 40 kDa FITC-Dextran were added to the inner chamber, and samples were taken from the outer chamber at different time points. Diffusion of FITC-Dextran across the scaffold was calculated as a percentage of total FITC-Dextran. The volume of the outer chamber was 68.75% of the total volume. All three matrices were permeable to low molecular weight FITC-Dextran (Fig. 4i). At 24 hours, significantly less 4 kDa FITC-Dextran had diffused across hPerMa, compared to pSIS and pPerMa (34.83% vs. 56.33% and 57.45%, respectively). The diffusion rate, given by the apparent permeability coefficient (Papp), was significantly lower for the hPerMa compared to pSIS and pPerMa for the complete period (p = 0.0078 and p = 0.0061, respectively; Fig. 4k). The matrices were only partially permeable to high molecular weight FITC-Dextran (Fig. 4j). At 24 hours, 4.00% for pSIS, 4.01% for pPerMa, and 2.72% for hPerMa of 40 kDa FITC-Dextran diffused, with no statistically significant difference in amount diffused FITC-Dextran or Papp value. For both low and high molecular weight FITC-Dextran, the diffusion rate was highest within the initial 60 minutes (Fig. 4 k,l). Interestingly, diffusion of 4 kD across the human PerMa was slower compared to the pPerMa and pSIS, indicating a lower permeability of the hPerMa.

3.5 The dECM scaffolds have unique matrisome compositions

To characterize the biomolecular composition of the native tissues and decellularized scaffolds, we performed label-free quantitative proteomics by liquid chromatography-tandem mass spectrometry (LC-MS/MS). The number of identified proteins varied between the groups ranging from approximately 2000 to 7000 proteins and was consistently lower in the dECM scaffolds compared to the corresponding native tissues (Supplementary Fig. 4a). This is expected after the removal of cellular content (50). The number of identified proteins was significantly higher in pSIS compared to pPerMa, hPerMa, and ocPerMa (Supplementary Fig. 4b).

Principal component analysis (PCA) showed that the biological replicates were similar (Fig. 5a,b). Both native and decellularized replicates clustered into separate groups, showing the reproducibility of the decellularization method. While the native human peritoneum showed variation between replicates, the biological replicates of decellularized hPerMa were more similar.

• Download figure
• Open in new tab

Figure 5.

Proteomic analysis by liquid-chromatography mass spectrometry (LC-MS/MS). a, b Principal component analysis (PCA) plots of the total proteome for porcine (a) and human (b) samples. c-f Composition of total matrisome proteins and the ten most abundant matrisome proteins in pSIS (n = 10), pPerMa (n = 11), hPerMa (n = 8), and ocPerMa (n = 2). The ratios of the top ten proteins are given as iBAQ intensity relative to the total iBAQ intensity (riBAQ) for each sample. Data represent mean ± SD. g Volcano plot of protein intensities in pSIS (n = 10) vs. pPerMa (n = 11). Matrisome proteins with at least a 3-fold change (log2) are labeled. h Profile bar charts for pSIS (n = 10), pPerMa (n = 11), hPerMa (n = 8), and ocPerMa (n = 2) of selected core matrisome proteins (n = 28). The proteins are sorted into collagens (top box), glycoproteins (middle box), and proteoglycans (bottom box). The ratios are given as relative iBAQ intensities. Data represent mean ± SD.

To evaluate whether key ECM proteins were retained after decellularization, we compared the abundance of core matrisome proteins and matrisome-associated factors, as defined by Naba et al. (37), in native and decellularized samples. The relative abundance of matrisome and matrisome-associated proteins vs. non-matrisome proteins was higher in all dECM scaffolds compared to their native counterparts (Supplementary Fig. 4c). The proportion of core matrisome and matrisome-associated proteins increased from 12% in pIntestine to 18% in pSIS, from 27% in pPeritoneum to 57% in pPerMa, from 28% in hPeritoneum to 66% in hPerMa, and from 13% in ocPeritoneum to 45% in ocPerMa. Concordance was observed between native and decellularized tissues when comparing the 10 most abundant matrisome proteins (Fig. 5c-f, Supplementary Fig. 4d). Volcano plots comparing native and decellularized tissues revealed enrichment of particularly core matrisome proteins in the dECM scaffolds. However, there was also a depletion of several matrisome proteins, especially matrisome-associated proteins (Supplementary Fig. 4 e-g). In summary, the decellularization resulted in the enrichment of matrisome proteins and retained the most abundant ECM proteins.

We then compared the different dECM scaffolds. The composition of matrisome and matrisome-associated proteins was similar across tissues, with some organ-specific differences (Fig. 5c-f). All dECM scaffolds were rich in glycoproteins and proteoglycans. The highest proportion was seen in hPerMa, where glycoproteins and proteoglycans together constituted 76% of the total matrisome proteins (Fig. 5e). This proportion was 59% for pSIS, 68% for pPerMa, and 66% for ocPerMa. The samples displayed great diversity in the composition of different categories of matrisome proteins. Collagen Type VI (COL6A1-3), glycoproteins such as dermatopontin (DPT), and proteoglycans such as lumican (LUM) and decorin (DCN) were found among the top 10 most abundant matrisome proteins. Several of these were identified in the top 10 lists across all groups. Despite the similarities in matrisome composition across all dECM scaffolds, the volcano plot comparing pSIS and pPerMa revealed differences in the abundance of non-matrisome and matrisome proteins (Fig. 5g). A panel of 28 core matrisome proteins revealed further differences in the relative abundance between the dECM scaffolds (Fig. 5h).

In summary, porcine PerMa represents a biologically relevant tissue scaffold with similarities to human PerMa in structure, physical properties, and biomolecular composition.

3.6 The pPerMa and the pSIS enable 3D in vitro cell cultures

Porcine tissues have several advantages over human tissues, including availability and less batch-to-batch variation. Therefore, porcine scaffolds were used to establish three-dimensional (3D) in vitro cell cultures. We first investigated the cytocompatibility of pSIS and pPerMa, and tested next, whether the ECM-specific differences in mechanical properties and protein composition of pSIS and pPerMa affected cellular behavior in vitro.

The matrices were mounted on the crowns, and cells were seeded on the surface of the inner chamber (Fig. 6a). Three EOC cell lines – OV90, Caov-3, and COV318 – were chosen based on their genetic resemblance to HGSOC (51). After three days, the scaffolds were flipped, leaving the cells on the outer surface where they were accessible for confocal imaging (Fig. 6a). The scaffolds were imaged on days 3, 7, 10, 14, and 21 after seeding (Fig. 6a). The surfaces function in the IMARIS image analysis software was used to estimate cell volume based on the fluorescent signal.

• Download figure
• Open in new tab

Figure 6.

3D cell cultures are established with the porcine dECM scaffolds. a Schematic illustration of the setup for cell culture. Scaffolds were mounted in a Transwell®-like setup and cells were seeded in the inner chamber. After three days, the scaffolds were flipped, allowing the cells to come in direct contact with the imaging surface, enabling confocal imaging. Scaffolds were kept in this reversed position in the culture plate between imaging. b Confocal images and representative images of histological staining of cell cultures on pSIS and pPerMa. Each row represents a cell line (COV318, OV90, and Caov3). COV318 and OV90 cells were seeded in a number of 25,000 cells/scaffold, while Caov3 cells were seeded at 400,000 cells/scaffold for monoculture and 200,000 cells/scaffold for coculture. For cocultures, BJ fibroblasts were seeded in a number of 12,500 cells/scaffold (COV318), 6,250 cells/scaffold (OV90), and 25,000 cells/scaffold (Caov3). Confocal images are 3D renderings of z-stacks, visualized in the x/y plane (left, quadratic images) and the x/z plane (right, rectangular images). Magenta represents tissue autofluorescence (ex = 561 nm, em = 575-625 nm), green represents GFP+ EOC cell lines (ex = 488 nm, em = 500-550 nm), and white represents miRFP670+ BJ fibroblasts (ex = 637 nm, em = 663-738 nm). Scale bars = 100 µm. Histological stainings were performed on a cross-section of the scaffolds, corresponding to the x/z plane in the confocal images. The monocultures are stained with EpCAM. Cocultures are stained with EpCAM (top) and hematoxylin & eosin (H&E; bottom). Scale bars = 50 µm. NA = not applicable, too few cells to stain. The experiment was performed two times with three replicates for each group. c Graphs showing the cell growth of COV318, OV90, and Caov3 cells on pSIS (monoculture) and pPerMa (monoculture and coculture with fibroblasts). The cell volume was estimated using the surfaces function in the IMARIS image analysis software. Data represent mean ± SD.

Both pSIS and pPerMa were cytocompatible and allowed cell growth. The three cell lines displayed distinct growth patterns. COV318 cells formed a monolayer on the surface of both pSIS and pPerMa, with similar growth kinetics (Fig. 6b, c, Supplementary Fig. 5a). There was no significant difference in cell volume between pSIS and pPerMa at any of the time points. OV90 cells formed superficial three-dimensional tumor-like structures, with an increased growth on pPerMa compared to pSIS (Fig. 6b, c, Supplementary Fig. 5b). The difference was statistically significant at all time points (p values 0.0001 day 3, 0.0032 day 7, 0.0101 day 10, 0.0004 day 14, and 0.0014 day 21 for unpaired, two-tailed t-test). Caov3 cells seeded in the same number as COV318 and OV90 did not grow, neither on the pSIS nor the pPerMa. When seeded at 400,000 cells/scaffold instead of 25,000, some cell growth was observed on pSIS but not on pPerMa (Fig. 6b, c, Supplementary Fig. 5c).

Several studies have found that cancer-associated fibroblasts (CAFs) promote peritoneal metastasis of EOC by contributing to adhesion and invasion (52–54). Therefore, we incorporated the human fibroblast cell line BJ in the 3D cell cultures on pPerMa. The BJ cell line was transduced with miRFP670, a far-red fluorophore, permitting simultaneous confocal imaging of GFP-positive tumor cells and fibroblasts. Interestingly, we found that Caov3 cells, which were unable to grow on pPerMa in monoculture, did grow in coculture with fibroblasts (Fig. 6b, c, Supplementary Fig. 6c). Moreover, the tumor cells formed clusters surrounded by fibroblasts, a pattern that is also seen in vivo (55). COV318 cells in coculture with fibroblasts grew as a monolayer separated from the fibroblasts, while OV90 cells grew interspersed with fibroblasts (Fig. 6b, Supplementary Fig. 6a, b).

3.7 The 3D cell cultures can be used for therapy sensitivity testing

To assess sensitivity to therapy in the 3D in vitro model, we developed a method to assess cell death using confocal microscopy. For the establishment of the technique, we treated the cells with carboplatin, a chemotherapeutical used in the treatment of EOC. To select carboplatin doses to be tested in the 3D system, we determined the EC50 value of carboplatin for the OV90 cell line in 2D culture by WST-1 assay (Supplementary Fig. 7a). The cells were treated for 72 hours with carboplatin at doses ranging from 1 µM to 3160 µM and an EC50 value of 45.60 µM for OV90 cells was determined. To validate the EC50 value, we performed a set of experiments to evaluate cell death by Annexin V/PI staining and flow cytometry in 2D cell culture (Supplementary Fig. 7b,c). The cells were treated for 72 hours with carboplatin doses ranging from 10 µM to 316 µM, including the EC50 dose of 45.6 µM. The cells displayed increasing cell death, defined as positive staining for Annexin V and/or PI (Q1-Q3) with increasing drug doses. The EC50 dose group had 54.5% cell death.

To assess cell death in 3D, we used confocal microscopy and stained for live, apoptotic, and necrotic markers (Supplementary Fig. 7d, e). OV90 cells cultured on pPerMa were treated with carboplatin at three different doses (10 µM, 50 µM, and 100 µM) for 72 hours. The scaffolds were stained with CytoCalcein Violet, labeling the cytoplasm of viable cells, Apopxin Green, a phosphatidylserine (PS) sensor labeling membranes of apoptotic cells, and 7-AAD, a DNA intercalator labeling cell nuclei in membrane-permeable necrotic cells. All dyes stained OV90 cells growing on pPerMa. The fluorescent signal could be quantified by the IMARIS image analysis software using the surfaces function for live cells and the spots function for apoptotic and necrotic cells (Supplementary Fig. 7e). However, the Apopxin Green dye also stained the scaffold itself and the staining of apoptotic cells was challenging to distinguish from background staining. For further therapy sensitivity testing, we chose to continue with only markers for live and necrotic cells.

To investigate the feasibility of our model system to assess sensitivity to novel cell-based therapies, we tested co-culture with CAR T cells. A CAR T cell product, directed towards tumor-associated glycoprotein 72 (TAG-72), was tested on TAG-72 positive OV90 cells. The CAR and Mock T cells were generated from PBMCs of three different healthy donors. The EOC cell lines were cultured on pPerMa for 8 days before treatment (Fig. 7a). Imaging was done on day seven to ensure growth on all scaffolds and assign them to treatment groups. On day eight, CellTrace Violet-stained T cells were seeded. After 24 hours, the scaffolds were flipped, stained with 7-AAD, and imaged with confocal microscopy. The volume of live tumor cells was estimated from the GFP+ signal using the surfaces function in the IMARIS image analysis software. The confocal images show the difference in live cells (green) and dead cells (magenta) in the treated and untreated groups (Fig. 7b). After 24 hours of treatment, there was a significant difference in live cell volume between OV90 cells treated with TAG-72 targeted CAR T cells (n = 9) compared to Mock T cells (n = 9, p = 0.0015) and untreated OV90 cells (n = 3, p = 0.0111; Fig. 7c). There was also a significant difference in live cell volume between OVCAR3 cells treated with MUC-16 targeted CAR T cells (n = 9) compared to Mock T cells (n = 9, p = 0.0149) and untreated (n = 3, p = 0.0054; Fig. 7d). In summary, our 3D in vitro model system on pPerMa can be used to assess sensitivity to traditional and novel therapies against EOC.

• Download figure
• Open in new tab

Figure 7.

CAR T cell sensitivity is assessed in the 3D cell cultures using confocal microscopy. a Schematic illustration of the setup for treatment with CAR T or Mock T cells. The EOC cell lines OV90 and OVCAR3 were seeded on PerMa. OV90 cells were seeded at 25,000 cells/scaffold and OVCAR3 cells at 100,000 cells/scaffold. After three days, the scaffolds were flipped and imaged with confocal microscopy. After four more days, the scaffolds were imaged again and flipped back to the original position with the cells in the inner chamber. The next day, CAR T or Mock T cells were seeded on the scaffolds. After 24 hours, the scaffolds were flipped, stained with 7-AAD, and imaged with confocal microscopy. b,d Representative confocal images of OV90 cells treated with CAR T cells targeting TAG72 or Mock T cells (b), and OVCAR3 cells treated with CAR T cells targeting Muc16 or Mock T cells (d). CAR T and Mock T cells were derived from three different donors. The experiment was performed with three replicates in each group. Magenta represents 7-AAD (ex = 561 nm, em = 575-625 nm), green represents GFP+ EOC cell lines (ex = 488 nm, em = 500-550 nm), and white represents CellTraceTM Violet stained T cells (ex = 405 nm, em = 425-475 nm). The T cells that were adhered to the scaffold and/or tumor cells, did not detach when the scaffolds were flipped and can be seen in white on the confocal images. Scale bars = 100 µm. NA = not applicable. c,e Graphs showing the estimated cell volume of OV90 cells (c) and OVCAR3 cells (e) untreated and treated with Mock T cells or CAR T cells (left). Graphs showing the change in estimated cell volume from the day before treatment to 24 hours after treatment for Mock T cells (middle) and CAR T cells (right). The cell volume was estimated using the surfaces function in the IMARIS image analysis software. Data represent mean ± SD, * P < 0.05, ** P < 0.01. Different symbols in the graphs represent donor 1 (circle), donor 2 (square), and donor 3 (triangle).