Next Article in Journal
Electroconductive Nanofibrous Scaffolds Enable Neuronal Differentiation in Response to Electrical Stimulation without Exogenous Inducing Factors
Previous Article in Journal
A Virtual Inner Ear Model Selects Ramped Pulse Shapes for Vestibular Afferent Stimulation
Previous Article in Special Issue
Modeling and Validation of an Ultra-Compact Regenerative Liver Dialysis Device
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Role of Desmoplastic Physical Stroma in Pancreatic Cancer Progression Using a Three-Dimensional Collagen Matrix Model

1
Tsukuba Life Science Innovation Program (T-LSI), School of Comprehensive Human Sciences, University of Tsukuba, Tsukuba 305-8572, Japan
2
Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba 305-8565, Japan
3
Department of Research Promotion and Management, National Cerebral and Cardiovascular Center, Kishibe-Shimmachi, Suita 564-8565, Japan
4
Department of Clinical Chemistry and Laboratory Medicine, Kyushu University Hospital, Fukuoka 812-8582, Japan
5
Department of Clinical Chemistry and Laboratory Medicine, Graduate School of Medical Sciences, Kyushu University, Fukuoka 812-8582, Japan
6
School of Integrative & Global Majors, University of Tsukuba, Tsukuba 305-8572, Japan
*
Author to whom correspondence should be addressed.
Bioengineering 2023, 10(12), 1437; https://doi.org/10.3390/bioengineering10121437
Submission received: 27 November 2023 / Revised: 13 December 2023 / Accepted: 14 December 2023 / Published: 18 December 2023
(This article belongs to the Special Issue The New Frontiers of Artificial Organs Engineering)

Abstract

:
Pancreatic ductal adenocarcinoma (PDAC) is a refractory tumor with a poor prognosis, and its complex microenvironment is characterized by a fibrous interstitial matrix surrounding PDAC cells. Type I collagen is a major component of this interstitial matrix. Abundant type I collagen promotes its deposition and cross-linking to form a rigid and dense physical barrier, which limits drug penetration and immune cell infiltration and provides drug resistance and metabolic adaptations. In this study, to identify the physical effect of the stroma, type I collagen was used as a 3D matrix to culture Capan-1 cells and generate a 3D PDAC model. Using transcriptome analysis, a link between type I collagen-induced physical effects and the promotion of Capan-1 cell proliferation and migration was determined. Moreover, metabolomic analysis revealed that the physical effect caused a shift in metabolism toward a glycolytic phenotype. In particular, the high expression of proline in the metabolites suggests the ability to maintain Capan-1 cell proliferation under hypoxic and nutrient-depleted conditions. In conclusion, we identified type I collagen-induced physical effects in promoting Capan-1 cells, which cause PDAC progression, providing support for the role of dense stroma in the PDAC microenvironment and identifying a fundamental method for modeling the complex PDAC microenvironment.

1. Introduction

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer, accounting for >85% of all diagnoses [1]. PDAC is a refractory and aggressive malignancy with a dense fibrotic stroma that is associated with poor prognosis [2,3], reflected in a 5-year survival rate of <10% [4]. This dense stroma presents as a fibrous interstitial matrix surrounding PDAC cells, which limits drug penetration and immune cell infiltration and provides drug resistance and metabolic adaptations to PDAC [5].
Collagen, primarily types I, III, IV, V, and XV, serve as key components of the extracellular matrix (ECM) in PDAC [6]. Of these, type I collagen is the most abundant and is present three times higher than that in the normal pancreas [7,8,9,10]. It tends to deposit and crosslink excessively in the matrix, forming a stiff, physical barrier. Moreover, type I collagen is responsible for most desmoplastic reactions, disrupting the structure of the basement membrane and inducing PDAC cells to increase interstitial collagen levels, resulting in the deterioration of PDAC [11]. Type I collagen can interact with integrin α2β1 expressed on the surface of PDAC cells, promoting their proliferation and migration [6,10]. This process results in the loss of the E-cadherin complex and cell–cell adhesion, an important step in metastasis and epithelial-to-mesenchymal transition (EMT) [12]. Thus, the composition of the ECM plays a pivotal role in driving EMT, and a type I collagen-rich microenvironment favors a more aggressive PDAC phenotype.
Collagen-based 3D cultures have profound epigenetic effects on PDAC cells [11]. When PDAC cells are cultured in a 3D collagen matrix, the expression of histone acetyltransferases p300, P300/CBP-associated factor, and GCN5 are increased [13], which is accompanied by an elevated level of high-mobility group A2, an epigenetic regulator of proliferation, apoptosis, and DNA repair [14]. In addition, collagen-based 3D cultures promote angiogenesis and the secretion of growth factors, such as epidermal growth factor (EGF), matrix metalloproteinase 9 (MMP9), and vascular endothelial growth factor (VEGF) [15]. However, our understanding of epigenetic modifications remains incomplete, highlighting the need for further studies surrounding the impact of collagen-based 3D cultures on PDAC cells.
In this study, a type I collagen-based 3D PDAC model containing Capan-1 cells (a PDAC cell line) was constructed to better understand the impact of collagen on the progression of PDAC tumor microenvironments. This 3D culture model is more informative than 2D culture systems for exploring the spatial structure of tumors, growth patterns, and cell–cell interactions [16]. Remarkably, clinical-like PDAC glandular structures were formed inside the model, providing a platform for simple simulation of stromal physical effects on PDAC cells. To investigate the pathogenesis of collagen-based tumor tissues, we examined transcriptome changes in 2D and collagen-based 3D PDAC cultures. Furthermore, to comprehensively explore the type I collagen-induced physical effects at advanced stages of PDAC, we performed metabolomic analysis to identify key metabolites that influence biological pathways. Our findings suggest that type I collagen and its physical effects are crucial factors in the development of PDAC and that this 3D PDAC collagen model has significant implications for understanding the intricate mechanisms involving type I collagen and PDAC.

2. Materials and Methods

2.1. Cell and 2D Culturing Conditions

The human pancreatic cancer cell line Capan-1 (ATCC HTB-79) was cultured as per the previous protocol [17,18]. Briefly, Capan-1 cells were maintained in Dulbecco’s Modified Eagle’s medium (DMEM; FUJIFILM Wako Pure Chemical Corp., Osaka, Japan) supplemented with 20% fetal bovine serum (FBS), 1% nonessential amino acids, and 1% streptomycin–penicillin at 37 °C in a humidified atmosphere containing 5% CO2.

2.2. Preparation of PDAC-Embedded 3D Collagen Tissue Matrix

To establish the collagen-based 3D PDAC model, the Collagen–Capan-1 model, PDAC cells (1 × 105/mL) were suspended in neutralized type I collagen solution (IAC-30; Koken Co., Tokyo, Japan). The collagen mixture was poured into tissue culture plates pre-coated with 2-methacryloyloxyethyl phosphorylcholine (MPC) polymer (LIPIDURE-CM5206; NOF Co., Tokyo, Japan). The volume of collagen mixture was 800 μL/well when using 24-well plates and 200 μL/well when using 96-well plates. The corresponding volume of collagen solution was added to prepare collagen-only controls. Then, the plates were incubated at 37 °C for 30 min for gelation, and the gels were released from the wells using a vortex mixer. Additional fresh medium of 1000 μL and 100 μL was added to 24- and 96-well plates, respectively. The plates were placed in a 37 °C incubator with shaking at 95 rpm for extended culture. The medium was changed every 2–3 days and supplemented with 10% FBS, 1% nonessential amino acids, and 1% streptomycin–penicillin. Gel contraction was calculated as the percentage of lattice area relative to the initial gel area using ImageJ (1.53t Java 1.8.0_322, U.S. National Institutes of Health, Bethesda, MD, USA) software.

2.3. RNA Sequencing Analysis

RNA was extracted from Capan-1 cells that were normally cultured in dishes, referred to as the 2D Capan-1 cells, using NucleoSpin RNA (Macherey Nagel GmbH & Co., KG, Duren, Germany) and from the Collagen–Capan-1 model after 3 h of incubation using ISOGEN reagent (Nippon Gene; Tokyo, Japan), following the manufacturer’s instructions. RNA sequencing analysis of 2D Capan-1 and Collagen–Capan-1 was performed as previously described [19]. To elucidate the functions of differentially expressed genes, Gene Ontology (GO) enrichment analysis was performed using the online database for annotation, visualization, and integrated discovery (DAVID). Furthermore, Gene Set Enrichment Analysis (GSEA) was performed using GSEA software (GSEA software 4.2.3, Broad Institute, Cambridge, MA, USA). The raw sequences in FASTQ format are available from DDBJ (DRA017362).

2.4. Histological Analysis

To preserve the collagen matrix structure, the Collagen–Capan-1 model was encapsulated using iPGell (Genostaff, Tokyo, Japan). Paraffin-embedded sections with a thickness of 5 μm were then prepared. Hematoxylin and eosin (H&E) staining was performed for internal morphological analysis according to standard protocols.

2.5. Subcutaneous Transplantation In Vivo

Seven-week-old female nude mice (BALB/c-nu/nu; CLEA Japan, Tokyo, Japan) were maintained at the Animal Center of the National Institute of Advanced Industrial Science and Technology (AIST). The Collagen–Capan-1 model was prepared in 24-well plates and harvested on Day 7. Subsequently, Matrigel with three models was gelated to integrate one tumor fragment at room temperature. This fragment was subsequently transplanted into the left thigh of mice (n = 3). On Day 49, the mice were euthanized by cervical dislocation, and all subcutaneous tumors were isolated. All invasive procedures were conducted under isoflurane anesthesia, following guidelines approved by the Institutional Animal Care and Use Committee of the respective institutes of AIST.

2.6. Metabolomic Analysis

Primary metabolites were analyzed using liquid chromatography–mass spectrometry (LC-MS). Briefly, metabolites were extracted by vortexing the culture medium with ice-cold methanol. Subsequently, the metabolites were purified and analyzed by LC-MS according to the specific experimental methodology of a previous study [20].

2.7. Statistical Analysis and Visualization

The dataset containing metabolite information was converted to CSV format and then uploaded onto the MetaboAnalyst® platform (accessible at [https://www.metaboanalyst.ca] accessed on 2 November 2023). This platform is equipped to perform extensive processing and analysis of metabolomic data. To ensure data integrity, a default data integrity check was performed, following a filtration process based on the average intensity values. A significance level was set at p < 0.05, and the Pareto data scaling was used for normalization when comparing two group samples. Orthogonal-partial-least-squares discrimination (OPLS-DA) was conducted. The OPLS-DA results were visualized using principal component analysis (PCA) and heatmap cluster analyses. Additionally, the visualization of variable important in projection (VIP) plots was employed following the OPLS-DA results to pinpoint discriminative metabolite markers.

3. Results

3.1. Construction of PDAC Cell-Embedded 3D Collagen Tissue Matrix

Capan-1 cells were embedded in a type I collagen gel to construct a collagen-based 3D PDAC model (Figure 1A), referred to as the Collagen–Capan-1 model. This model was used to investigate the effects of type I collagen on the progression of PDAC. As shown in Figure 1B, the type I collagen gel with no PDAC cells acted as a control group (collagen-only); the Collagen–Capan-1 model presented a disk structure on Days 7 and 14, in contrast to the loose structure of the collagen-only group, suggesting that interactions between type I collagen and Capan-1 cells occurred within the model.

3.2. Effect of Type I Collagen-Derived Physical Stress on PDAC Cells

To understand the impact of type I collagen-derived physical stress on Capan-1 cells during the initial modeling stages, RNA sequencing analysis was used to evaluate alterations in global gene expression patterns of Capan-1 cells following co-culture with type I collagen for 3 h. Compared to Capan-1 cells that were cultured in dishes, referred to as the 2D Capan-1 cells, 1103 genes were differentially expressed in the transcriptome of the Collagen–Capan-1 model, of which 597 genes were upregulated, and 606 were downregulated (Figure 2A). Representative upregulated (FGF1, MYC, and IL6) and downregulated genes (SEMA4D, PCDH20, and ISLR) are shown in Figure 2A. In addition, detailed upregulated genes associated with tumor proliferation are shown in Figure 2B. On the other hand, downregulated genes associated with reduced cell adhesion, signaling increased metastasis and migration of PDAC cells, are shown in Figure 2C.
To further explore alterations involving biological functions, upregulated and downregulated genes with significant differences were subjected to Gene Ontology (GO) analysis, and the top five biological processes (BP) were listed. Notably, although the same BPs were listed in both up and down situations, the content gene sets were different. The representative GO terms related to upregulated genes were “positive regulation of cell proliferation “and “positive regulation of gene expression” (Figure 3A, Table S1), and that related to downregulated genes was “cell adhesion” (Figure 3B, Table S2), indicating the crucial role of type I collagen in proliferation and metastasis of PDAC cells. In addition, GSEA analysis was performed to identify pathways enriched in the ranked gene lists. The results demonstrated that actin-related pathways were enriched at the initial modeling stages, such as “Actin polymerization or depolymerization”, “Actin mediated cell contraction”, and “Actomyosin structure organization” (Figure 3C–E).

3.3. Internal Histological Analysis of Collagen–Capan-1 Model

As depicted in Figure 4A, B, H&E staining was used to identify the internal morphology of the Collagen–Capan-1 model. Compared with the cell-free model (collagen only), the Collagen–Capan-1 model displayed clinic-like glandular structures of PDAC on days 7 and 14. Notably, the Collagen–Capan-1 model exhibited a more mature glandular structure on Day 14, which potentially resembled an advanced stage of PDAC. Notably, these ductal structures were more abundant than those in 2D-cultured Capan-1 cells, which showed no such ducted glandular structures (Figure 4C) and recapitulated the glandular structures that were built in tumor subcutaneous transplantations in the in vivo model (Figure 4D). These findings highlight the significant role of collagen 3D matrix in supporting glandular structure formation in PDAC.

3.4. Identifying Secreted and Consumed Extracellular Metabolites in Collagen–Capan-1 Model

To investigate the role of type I collagen-induced physical effects on the metabolic mechanism of Capan-1 in the 3D matrix, metabolomic analysis of water-soluble metabolites collected from culture supernatants was performed using LC-MS. OPLS-DA statistical analysis was utilized to determine differential metabolic patterns in the Collagen–Capan-1 group compared to the 2D Capan-1 group (Figure 5A). In addition, a heatmap was generated to provide a detailed understanding of water-soluble metabolites, highlighting the top 25 metabolites with significant variation (Figure 5B). Visualization of the VIP plot also emphasized distinct metabolites between the two groups (Figure 5C). In detail, in the culture supernatant of the Collagen–Capan-1 model, a higher accumulation of glutamic acid and lactic acid was detected (Figure 5D), indicating a metabolic shift of Capan-1 cells toward glycolysis strongly in the 3D collagen matrix [21,22]. Interestingly, proline, known to support PDAC growth under hypoxic conditions as a nutrient, also exhibited elevated levels. This suggests that collagen has the potential to provide a hypoxic tumor microenvironment that is conducive to PDAC cell growth. Consumption of essential nutrients such as lysine and glutamine (Figure 5E) was also observed in Collagen–Capan-1 supernatants. Collectively, these findings suggest that PDAC cells can adapt to the metabolic mechanisms of glycolysis in a 3D collagen matrix, resembling clinical PDAC metabolism. These findings underscore the vital role of type I collagen in the context of PDAC.

4. Discussion

In this study, we established a type I collagen-based 3D PDAC model, offering insights into the physical effects of type I collagen on the pathogenesis and progression of PDAC. Histologically, a clinical-like glandular structure of PDAC was spatially developed in the 3D collagen matrix, facilitating a realistic understanding of interactions between PDAC cells and the interplay of collagen on PDAC cells. Transcriptomic analysis was employed to elucidate the potential of type I collagen-induced physical effects in promoting tumor proliferation and migration during the early modeling stage, mimicking the malignant role of elevated type I collagen in the pathogenic stage of PDAC. Furthermore, metabolomic analysis revealed that the type I collagen-based 3D matrix induced a metabolic shift toward a glycolytic phenotype in PDAC cells. Additionally, it provides proline to support the rapid proliferation of PDAC cells under hypoxic and nutrient-depleted conditions.
In GSEA analysis, actin-related biological properties were activated at the initial modeling stages, suggesting an interplay between collagen and PDAC cells [23,24]. Actin, an integral component of the cytoskeleton, is responsible for the maintenance and modification of cellular morphology [24,25]. Specifically, the dynamic reorganization of actin occurs under the influence of collagen, leading to a potential enhancement in cellular contraction, metastasis, and invasion by PDAC cells, which aligns with the reduction in “cell adhesion” in the GO analysis and the reduced expression of associated adhesion related molecule genes such as SEMA4D, PCDH20, and ISLR [26,27,28]. Moreover, Capan-1 cells may adjust their mechanical properties to adapt to the surrounding microenvironment, which is altered by collagen-induced mechanical stress via actin activation. This explains the positive relationship between mechanical stress in the tumor microenvironment and the capacity for proliferation and invasion by PDAC cells [29]. Collectively, these results provide a theoretical basis for collagen as a potential therapeutic target in PDAC.
PDAC cells, like other epithelial cancer cells, require substantial amounts of energy to support rapid proliferation. To meet this demand, PDAC cells can shift their metabolic pathways based on different tumor microenvironments. The adaptability of PDAC cells is a critical factor that contributes to the difficulty of PDAC treatment [30]. Alterations in the metabolic processes of Capan-1 cells between the 2D culture and the 3D collagen matrix also confirmed this hypothesis. In the Collagen–Capan-1 culture medium, a notable accumulation of lactic acid and glutamic acid was observed, suggesting that the high proliferation of Capan-1 cells relies on the glycolytic pathway, commonly known as the Warburg effect [31,32]. The excessive consumption of nutrients like lysine, valine, and histidine in the Collagen–Capan-1 medium, compared to the 2D Capan-1 culture, suggests a promotional effect of mechanical stress on the proliferation of Capan-1 cells. Moreover, increased proline expression implies that the physical effect and type I collagen can serve as a proline source to maintain PDAC cell proliferation under nutrient-depleted and hypoxic conditions [33,34]. Proline is catalyzed by proline dehydrogenase 1 (PEODH1) to generate glutamine, which is involved in the tricarboxylic acid (TCA) cycle for ATP production [20,35,36]. Furthermore, the notable reduction in tryptophan in Collagen–Capan-1 medium potentially indicates the activation of the tryptophan metabolic pathway in PDAC cells, implying that the mechanical stress of the 3D collagen matrix potentially supports the generation of immune evasion in this 3D PDAC model. Taken together, these findings demonstrate that PDAC cells employ a unique survival strategy in a 3D stromal environment, highlighting the matrix’s supportive role in PDAC progression.
Type I collagen is a primary component of ECM that contributes to drug resistance in PDAC [37]. For instance, it potentiates drug resistance to gemcitabine in pancreatic cancer cells [13,14]. First, type I collagen forms a physical barrier surrounding PDAC, impeding the delivery and efficacy of chemotherapeutic agents [38,39]. Additionally, type I collagen can lead to the development of PDAC resistance via interactions with cell-surface receptors and intercellular signaling pathways. Membrane-type 1-MMP (MT1-MMP) plays a central role in collagen-mediated drug resistance in PDAC [40]. Overexpression of MT1-MMP in the collagen microenvironment increases ERK1/2 phosphorylation and high-mobility group AT-hook 2 (HMGA2) expression, which further attenuates GEM-induced checkpoint arrest by acting as a means of drug resistance [41,42]. Consequently, type I collagen holds promise as a potential therapeutic target in PDAC. At the molecular level, a thorough understanding of the role of type I collagen in PDAC progression will be beneficial for future research focused on collagen-based drug resistance.
In this study, we conducted an in-depth exploration of type I collagen in PDAC, which is predominant in the ECM, to investigate its physical effects on PDAC cell proliferation. However, the absence of other cellular and extracellular matrix components in the Collagen–Capan-1 model makes it difficult to apply it to drug screening and research concerning disease mechanisms that require a higher demand for the simulation of PDAC microenvironments. Specifically, for instance, many cell types are presented in the PDAC stromal environment, including cancer-associated fibroblasts, immune cells, and epithelial cells. In addition, a diverse array of ECM components is included, comprising fibronectin, hyaluronic acid, laminin, matrix metalloproteinases [43], and multiple types of collagens, such as type III [44] and type IV [45]. Each of these components plays a significant role in PDAC progression and chemoresistance [46,47], and the incorporation of additional components into our 3D model offers the potential to simulate a more authentic PDAC stromal microenvironment. Therefore, to achieve a more realistic PDAC microenvironment, it is imperative to consider the inclusion of a wide range of PDAC components in our next research phase.

5. Conclusions

In the present study, we constructed a collagen-based 3D PDAC model with clinical glandular structures using Capan-1 cells. The physical role of type I collagen in promoting PDAC proliferation and migration was confirmed using molecular studies of global gene expression patterns and metabolomic analyses. These findings enhance our understanding of how collagen and mechanical stress, influenced by the stromal environment, affect the clinical PDAC microenvironment. This, in turn, facilitates further exploration into the remodeling of the PDAC microenvironment.

Supplementary Materials

The following supporting information can be downloaded from https://www.mdpi.com/article/10.3390/bioengineering10121437/s1. Table S1: Gene ontologies with the top five highest p-values for biological processes in the Collagen–Capan-1 model; Table S2: Gene ontologies with the top five lowest p-values of biological processes in the Collagen–Capan-1 model.

Author Contributions

Conceptualization, Y.S.K.; methodology, X.S., Y.N., M.Y., D.S., Y.K. and Y.S.K.; software, X.S., Y.N., M.Y., D.S., Y.K. and Y.S.K.; validation, X.S., Y.N. and Y.S.K.; formal analysis, X.S., Y.N. and Y.S.K.; data curation, X.S., Y.N. and Y.S.K.; writing—original draft preparation, X.S. and Y.S.K.; writing review and editing, X.S., Y.N., M.Y., D.S., Y.K. and Y.S.K.; visualization, X.S., Y.N. and Y.S.K.; supervision, Y.S.K.; project administration, Y.S.K.; funding acquisition, X.S. and Y.S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by JSPS KAKENHI (grant number: JP23K13919) and AMED (grant number: JP22ama121050j0001). X.S. is a recipient of the Japan Science and Technology Agency (JST)—SPRING Fellowship (grant number JPMJSP2124).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of the National Institute of Advanced Industrial Science and Technology (AIST) (protocol code: 2023-0310, date of approval: 1 July 2018).

Informed Consent Statement

Not acceptable.

Data Availability Statement

All other data supporting the findings of this study are available from the corresponding author upon request.

Acknowledgments

The authors thank the members of the Cellular and Molecular Biotechnology Research Institute. The authors are grateful to Yuka Akagi for the technical assistance and Yasuko Ozaki for secretarial assistance. Computations were partially performed on the NIG supercomputer at ROIS National Institute of Genetics.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rawla, P.; Sunkara, T.; Gaduputi, V. Epidemiology of Pancreatic Cancer: Global Trends, Etiology and Risk Factors. World J. Oncol. 2019, 10, 10–27. [Google Scholar] [CrossRef] [PubMed]
  2. Whatcott, C.J.; Diep, C.H.; Jiang, P.; Watanabe, A.; LoBello, J.; Sima, C.; Hostetter, G.; Shepard, H.M.; Von Hoff, D.D.; Han, H. Desmoplasia in Primary Tumors and Metastatic Lesions of Pancreatic Cancer. Clin. Cancer Res. 2015, 21, 3561–3568. [Google Scholar] [CrossRef] [PubMed]
  3. Shields, M.A.; Dangi-Garimella, S.; Redig, A.J.; Munshi, H.G. Biochemical Role of the Collagen-Rich Tumour Microenvironment in Pancreatic Cancer Progression. Biochem. J. 2012, 441, 541–552. [Google Scholar] [CrossRef] [PubMed]
  4. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer Statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef] [PubMed]
  5. Koenig, A.; Mueller, C.; Hasel, C.; Adler, G.; Menke, A. Collagen Type I Induces Disruption of E-Cadherin-Mediated Cell-Cell Contacts and Promotes Proliferation of Pancreatic Carcinoma Cells. Cancer Res. 2006, 66, 4662–4671. [Google Scholar] [CrossRef] [PubMed]
  6. Perez, V.M.; Kearney, J.F.; Yeh, J.J. The PDAC Extracellular Matrix: A Review of the ECM Protein Composition, Tumor Cell Interaction, and Therapeutic Strategies. Front. Oncol. 2021, 11, 751311. [Google Scholar] [CrossRef]
  7. Linder, S.; Castaños-Velez, E.; von Rosen, A.; Biberfeld, P. Immunohistochemical Expression of Extracellular Matrix Proteins and Adhesion Molecules in Pancreatic Carcinoma. Hepatogastroenterology 2001, 48, 1321–1327. [Google Scholar]
  8. Mollenhauer, J.; Roether, I.; Kern, H.F. Distribution of Extracellular Matrix Proteins in Pancreatic Ductal Adenocarcinoma and Its Influence on Tumor Cell Proliferation In Vitro. Pancreas 1987, 2, 14–24. [Google Scholar] [CrossRef]
  9. Imamura, T.; Iguchi, H.; Manabe, T.; Ohshio, G.; Yoshimura, T.; Wang, Z.H.; Suwa, H.; Ishigami, S.; Imamura, M. Quantitative Analysis of Collagen and Collagen Subtypes I, III, and V in Human Pancreatic Cancer, Tumor-Associated Chronic Pancreatitis, and Alcoholic Chronic Pancreatitis. Pancreas 1995, 11, 357–364. [Google Scholar] [CrossRef]
  10. Armstrong, T.; Packham, G.; Murphy, L.B.; Bateman, A.C.; Conti, J.A.; Fine, D.R.; Johnson, C.D.; Benyon, R.C.; Iredale, J.P. Type I Collagen Promotes the Malignant Phenotype of Pancreatic Ductal Adenocarcinoma. Clin. Cancer Res. 2004, 10, 7427–7437. [Google Scholar] [CrossRef]
  11. Weniger, M.; Honselmann, K.C.; Liss, A.S. The Extracellular Matrix and Pancreatic Cancer: A Complex Relationship. Cancers 2018, 10, 316. [Google Scholar] [CrossRef] [PubMed]
  12. Perl, A.-K.; Wilgenbus, P.; Dahl, U.; Semb, H.; Christofori, G. A Causal Role for E-Cadherin in the Transition from Adenoma to Carcinoma. Nature 1998, 392, 190–193. [Google Scholar] [CrossRef] [PubMed]
  13. Dangi-Garimella, S.; Sahai, V.; Ebine, K.; Kumar, K.; Munshi, H.G. Three-Dimensional Collagen I Promotes Gemcitabine Resistance In Vitro in Pancreatic Cancer Cells through HMGA2-Dependent Histone Acetyltransferase Expression. PLoS ONE 2013, 8, e64566. [Google Scholar] [CrossRef] [PubMed]
  14. Dangi-Garimella, S.; Krantz, S.B.; Barron, M.R.; Shields, M.A.; Heiferman, M.J.; Grippo, P.J.; Bentrem, D.J.; Munshi, H.G. Three Dimensional Collagen I Promotes Gemcitabine Resistance in Pancreatic Cancer through MT1-MMP-Mediated Expression of HMGA2. Cancer Res. 2011, 71, 1019–1028. [Google Scholar] [CrossRef] [PubMed]
  15. Biondani, G.; Zeeberg, K.; Greco, M.R.; Cannone, S.; Dando, I.; Dalla Pozza, E.; Mastrodonato, M.; Forciniti, S.; Casavola, V.; Palmieri, M.; et al. Extracellular Matrix Composition Modulates PDAC Parenchymal and Stem Cell Plasticity and Behavior through the Secretome. FEBS J. 2018, 285, 2104–2124. [Google Scholar] [CrossRef]
  16. Kapałczyńska, M.; Kolenda, T.; Przybyła, W.; Zajączkowska, M.; Teresiak, A.; Filas, V.; Ibbs, M.; Bliźniak, R.; Łuczewski, Ł.; Lamperska, K. 2D and 3D Cell Cultures—A Comparison of Different Types of Cancer Cell Cultures. Arch. Med. Sci. 2018, 14, 910–919. [Google Scholar] [CrossRef]
  17. Miyazaki, Y.; Oda, T.; Inagaki, Y.; Kushige, H.; Saito, Y.; Mori, N.; Takayama, Y.; Kumagai, Y.; Mitsuyama, T.; Kida, Y.S. Adipose-Derived Mesenchymal Stem Cells Differentiate into Heterogeneous Cancer-Associated Fibroblasts in a Stroma-Rich Xenograft Model. Sci. Rep. 2021, 11, 4690. [Google Scholar] [CrossRef] [PubMed]
  18. Miyazaki, Y.; Oda, T.; Mori, N.; Kida, Y.S. Adipose-Derived Mesenchymal Stem Cells Differentiate into Pancreatic Cancer-Associated Fibroblasts in Vitro. FEBS Open Bio 2020, 10, 2268–2281. [Google Scholar] [CrossRef]
  19. Imai, Y.; Mori, N.; Nihashi, Y.; Kumagai, Y.; Shibuya, Y.; Oshima, J.; Sasaki, M.; Sasaki, K.; Aihara, Y.; Sekido, M.; et al. Therapeutic Potential of Adipose Stem Cell-Derived Conditioned Medium on Scar Contraction Model. Biomedicines 2022, 10, 2388. [Google Scholar] [CrossRef]
  20. Nihashi, Y.; Song, X.; Yamamoto, M.; Setoyama, D.; Kida, Y.S. Decoding Metabolic Symbiosis between Pancreatic Cancer Cells and Cancer-Associated Fibroblasts Using Cultured Tumor Microenvironment. Int. J. Mol. Sci. 2023, 24, 11015. [Google Scholar] [CrossRef]
  21. Tao, J.; Yang, G.; Zhou, W.; Qiu, J.; Chen, G.; Luo, W.; Zhao, F.; You, L.; Zheng, L.; Zhang, T.; et al. Targeting Hypoxic Tumor Microenvironment in Pancreatic Cancer. J. Hematol. Oncol. 2021, 14, 14. [Google Scholar] [CrossRef] [PubMed]
  22. Guillaumond, F.; Leca, J.; Olivares, O.; Lavaut, M.-N.; Vidal, N.; Berthezène, P.; Dusetti, N.J.; Loncle, C.; Calvo, E.; Turrini, O.; et al. Strengthened Glycolysis under Hypoxia Supports Tumor Symbiosis and Hexosamine Biosynthesis in Pancreatic Adenocarcinoma. Proc. Natl. Acad. Sci. USA 2013, 110, 3919–3924. [Google Scholar] [CrossRef] [PubMed]
  23. Gardel, M.L.; Sabass, B.; Ji, L.; Danuser, G.; Schwarz, U.S.; Waterman, C.M. Traction Stress in Focal Adhesions Correlates Biphasically with Actin Retrograde Flow Speed. J. Cell Biol. 2008, 183, 999–1005. [Google Scholar] [CrossRef] [PubMed]
  24. Izdebska, M.; Zielińska, W.; Grzanka, D.; Gagat, M. The Role of Actin Dynamics and Actin-Binding Proteins Expression in Epithelial-to-Mesenchymal Transition and Its Association with Cancer Progression and Evaluation of Possible Therapeutic Targets. BioMed Res. Int. 2018, 2018, 4578373. [Google Scholar] [CrossRef] [PubMed]
  25. Matsuda, Y.; Ishiwata, T.; Kawamoto, Y.; Kawahara, K.; Peng, W.-X.; Yamamoto, T.; Naito, Z. Morphological and Cytoskeletal Changes of Pancreatic Cancer Cells in Three-Dimensional Spheroidal Culture. Med. Mol. Morphol. 2010, 43, 211–217. [Google Scholar] [CrossRef]
  26. Ch’ng, E.S.; Kumanogoh, A. Roles of Sema4D and Plexin-B1 in Tumor Progression. Mol. Cancer 2010, 9, 251. [Google Scholar] [CrossRef] [PubMed]
  27. Imoto, I.; Izumi, H.; Yokoi, S.; Hosoda, H.; Shibata, T.; Hosoda, F.; Ohki, M.; Hirohashi, S.; Inazawa, J. Frequent Silencing of the Candidate Tumor Suppressor PCDH20 by Epigenetic Mechanism in Non-Small-Cell Lung Cancers. Cancer Res. 2006, 66, 4617–4626. [Google Scholar] [CrossRef]
  28. Takahashi, M.; Kobayashi, H.; Mizutani, Y.; Hara, A.; Iida, T.; Miyai, Y.; Asai, N.; Enomoto, A. Roles of the Mesenchymal Stromal/Stem Cell Marker Meflin/Islr in Cancer Fibrosis. Front. Cell Dev. Biol. 2021, 9, 749924. [Google Scholar] [CrossRef]
  29. Hadden, M.; Mittal, A.; Samra, J.; Zreiqat, H.; Sahni, S.; Ramaswamy, Y. Mechanically Stressed Cancer Microenvironment: Role in Pancreatic Cancer Progression. Biochim. Biophys. Acta Rev. Cancer 2020, 1874, 188418. [Google Scholar] [CrossRef]
  30. Biancur, D.E.; Kimmelman, A.C. The Plasticity of Pancreatic Cancer Metabolism in Tumor Progression and Therapeutic Resistance. Biochim. Biophys. Acta Rev. Cancer 2018, 1870, 67–75. [Google Scholar] [CrossRef]
  31. Warburg, O. On the Origin of Cancer Cells. Science 1956, 123, 309–314. [Google Scholar] [CrossRef] [PubMed]
  32. Li, R.; Li, H.; Zhu, L.; Zhang, X.; Liu, D.; Li, Q.; Ni, B.; Hu, L.; Zhang, Z.; Zhang, Y.; et al. Reciprocal Regulation of LOXL2 and HIF1α Drives the Warburg Effect to Support Pancreatic Cancer Aggressiveness. Cell Death Dis. 2021, 12, 1106. [Google Scholar] [CrossRef] [PubMed]
  33. Kamphorst, J.J.; Nofal, M.; Commisso, C.; Hackett, S.R.; Lu, W.; Grabocka, E.; Vander Heiden, M.G.; Miller, G.; Drebin, J.A.; Bar-Sagi, D.; et al. Human Pancreatic Cancer Tumors Are Nutrient Poor and Tumor Cells Actively Scavenge Extracellular Protein. Cancer Res. 2015, 75, 544–553. [Google Scholar] [CrossRef]
  34. Olivares, O.; Mayers, J.R.; Gouirand, V.; Torrence, M.E.; Gicquel, T.; Borge, L.; Lac, S.; Roques, J.; Lavaut, M.-N.; Berthezène, P.; et al. Collagen-Derived Proline Promotes Pancreatic Ductal Adenocarcinoma Cell Survival under Nutrient Limited Conditions. Nat. Commun. 2017, 8, 16031. [Google Scholar] [CrossRef] [PubMed]
  35. Liu, W.; Le, A.; Hancock, C.; Lane, A.N.; Dang, C.V.; Fan, T.W.-M.; Phang, J.M. Reprogramming of Proline and Glutamine Metabolism Contributes to the Proliferative and Metabolic Responses Regulated by Oncogenic Transcription Factor C-MYC. Proc. Natl. Acad. Sci. USA 2012, 109, 8983–8988. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, W.; Phang, J.M. Proline Dehydrogenase (Oxidase), a Mitochondrial Tumor Suppressor, and Autophagy under the Hypoxia Microenvironment. Autophagy 2012, 8, 1407–1409. [Google Scholar] [CrossRef]
  37. Liang, C.; Shi, S.; Meng, Q.; Liang, D.; Ji, S.; Zhang, B.; Qin, Y.; Xu, J.; Ni, Q.; Yu, X. Complex Roles of the Stroma in the Intrinsic Resistance to Gemcitabine in Pancreatic Cancer: Where We Are and Where We Are Going. Exp. Mol. Med. 2017, 49, e406. [Google Scholar] [CrossRef]
  38. Yilmaz, M.; Christofori, G. EMT, the Cytoskeleton, and Cancer Cell Invasion. Cancer Metastasis Rev. 2009, 28, 15–33. [Google Scholar] [CrossRef]
  39. Makareeva, E.; Han, S.; Vera, J.C.; Sackett, D.L.; Holmbeck, K.; Phillips, C.L.; Visse, R.; Nagase, H.; Leikin, S. Carcinomas Contain a Matrix Metalloproteinase-Resistant Isoform of Type I Collagen Exerting Selective Support to Invasion. Cancer Res. 2010, 70, 4366–4374. [Google Scholar] [CrossRef]
  40. Yeon, S.-E.; No, D.Y.; Lee, S.-H.; Nam, S.W.; Oh, I.-H.; Lee, J.; Kuh, H.-J. Application of Concave Microwells to Pancreatic Tumor Spheroids Enabling Anticancer Drug Evaluation in a Clinically Relevant Drug Resistance Model. PLoS ONE 2013, 8, e73345. [Google Scholar] [CrossRef]
  41. Senthebane, D.A.; Jonker, T.; Rowe, A.; Thomford, N.E.; Munro, D.; Dandara, C.; Wonkam, A.; Govender, D.; Calder, B.; Soares, N.C.; et al. The Role of Tumor Microenvironment in Chemoresistance: 3D Extracellular Matrices as Accomplices. Int. J. Mol. Sci. 2018, 19, 2861. [Google Scholar] [CrossRef] [PubMed]
  42. Sterzyńska, K.; Klejewski, A.; Wojtowicz, K.; Świerczewska, M.; Nowacka, M.; Kaźmierczak, D.; Andrzejewska, M.; Rusek, D.; Brązert, M.; Brązert, J.; et al. Mutual Expression of ALDH1A1, LOX, and Collagens in Ovarian Cancer Cell Lines as Combined CSCs- and ECM-Related Models of Drug Resistance Development. Int. J. Mol. Sci. 2018, 20, 54. [Google Scholar] [CrossRef] [PubMed]
  43. Xu, S.; Xu, H.; Wang, W.; Li, S.; Li, H.; Li, T.; Zhang, W.; Yu, X.; Liu, L. The Role of Collagen in Cancer: From Bench to Bedside. J. Transl. Med. 2019, 17, 309. [Google Scholar] [CrossRef] [PubMed]
  44. Di Martino, J.S.; Nobre, A.R.; Mondal, C.; Taha, I.; Farias, E.F.; Fertig, E.J.; Naba, A.; Aguirre-Ghiso, J.A.; Bravo-Cordero, J.J. A Tumor-Derived Type III Collagen-Rich ECM Niche Regulates Tumor Cell Dormancy. Nat. Cancer 2022, 3, 90–107. [Google Scholar] [CrossRef] [PubMed]
  45. Öhlund, D.; Franklin, O.; Lundberg, E.; Lundin, C.; Sund, M. Type IV Collagen Stimulates Pancreatic Cancer Cell Proliferation, Migration, and Inhibits Apoptosis through an Autocrine Loop. BMC Cancer 2013, 13, 154. [Google Scholar] [CrossRef]
  46. Verderio, E.A.M.; Johnson, T.; Griffin, M. Tissue Transglutaminase in Normal and Abnormal Wound Healing: Review Article. Amino Acids 2004, 26, 387–404. [Google Scholar] [CrossRef]
  47. Lee, J.; Condello, S.; Yakubov, B.; Emerson, R.; Caperell-Grant, A.; Hitomi, K.; Xie, J.; Matei, D. Tissue Transglutaminase Mediated Tumor-Stroma Interaction Promotes Pancreatic Cancer Progression. Clin. Cancer Res. 2015, 21, 4482–4493. [Google Scholar] [CrossRef]
Figure 1. Construction of PDAC-embedded 3D collagen tissue matrix. (A) Graphical representation of PDAC collagen model establishment. (B) Physical characteristics of PDAC collagen model. The Collagen-only model was used as a control. Both groups were modeled in 24-well plates and observed on Days 1, 7, and 14.
Figure 1. Construction of PDAC-embedded 3D collagen tissue matrix. (A) Graphical representation of PDAC collagen model establishment. (B) Physical characteristics of PDAC collagen model. The Collagen-only model was used as a control. Both groups were modeled in 24-well plates and observed on Days 1, 7, and 14.
Bioengineering 10 01437 g001
Figure 2. Differential gene expression analysis of PDAC cells triggered by collagen 3D culturing at the initial modeling phase, as revealed by RNA-seq analysis (n = 1, respectively). (A) The moving average (MA) plot displays differentially expressed genes in PDAC cells following 3 h culture in collagen 3D matrix. Expressed genes with significant differences were adjusted by p < 0.05 and log2 [fold-change] ≥ 2. The gene expression of PDAC cells in 2D culture as a basis. Upregulated genes are represented in red (n = 597), while downregulated genes are represented in blue (n = 606). (B) Representative upregulated markers associated with cell proliferation in the Collagen–Capan-1 model. (C) Representative downregulated markers associated with cell adhesion in the Collagen–Capan-1 model. The trimmed mean of M values (TMM) expression level in Collagen-Capan-1 is represented by red bars, and the corresponding TMM expression in 2D Capan-1 is represented by blue bars.
Figure 2. Differential gene expression analysis of PDAC cells triggered by collagen 3D culturing at the initial modeling phase, as revealed by RNA-seq analysis (n = 1, respectively). (A) The moving average (MA) plot displays differentially expressed genes in PDAC cells following 3 h culture in collagen 3D matrix. Expressed genes with significant differences were adjusted by p < 0.05 and log2 [fold-change] ≥ 2. The gene expression of PDAC cells in 2D culture as a basis. Upregulated genes are represented in red (n = 597), while downregulated genes are represented in blue (n = 606). (B) Representative upregulated markers associated with cell proliferation in the Collagen–Capan-1 model. (C) Representative downregulated markers associated with cell adhesion in the Collagen–Capan-1 model. The trimmed mean of M values (TMM) expression level in Collagen-Capan-1 is represented by red bars, and the corresponding TMM expression in 2D Capan-1 is represented by blue bars.
Bioengineering 10 01437 g002
Figure 3. Alterations in global gene expression patterns of PDAC cells identified via RNA-seq analysis after 3 h incubation with type I collagen (n = 1, respectively). (A) GO terms (biological process) enriched in the Collagen–Capan-1 model. The red bar indicates GO terms related to promotional tumor functions, such as cell proliferation. (B) GO terms (biological process) reduced in the Collagen–Capan-1 model. The red bar indicates GO terms related to tumor migration and metastasis functions, such as cell adhesion. (C) GSEA of Collagen–Capan-1 model, compared to 2D Capan-1 culture. Enrichment plots of expression signatures of “Actin polymerization or depolymerization”, (D) “Action mediated cell contraction”, and (E) “Actomyosin structure organization”.
Figure 3. Alterations in global gene expression patterns of PDAC cells identified via RNA-seq analysis after 3 h incubation with type I collagen (n = 1, respectively). (A) GO terms (biological process) enriched in the Collagen–Capan-1 model. The red bar indicates GO terms related to promotional tumor functions, such as cell proliferation. (B) GO terms (biological process) reduced in the Collagen–Capan-1 model. The red bar indicates GO terms related to tumor migration and metastasis functions, such as cell adhesion. (C) GSEA of Collagen–Capan-1 model, compared to 2D Capan-1 culture. Enrichment plots of expression signatures of “Actin polymerization or depolymerization”, (D) “Action mediated cell contraction”, and (E) “Actomyosin structure organization”.
Bioengineering 10 01437 g003
Figure 4. Internal histological analysis of 3D collagen PDAC model. (A) Representative images of H&E staining in Collagen–Capan-1 model on Day 7. (B) Representative images of H&E staining in Collagen–Capan-1 model on Day 14. D indicates duct formation. Scale bars = 200 μm, 200 μm, and 50 μm, respectively. (C) Morphology of Capan-1 cells in 2D culture. Scale bar, 200 μm. (D) Histological analysis of H&E-stained tumor cells generated subcutaneously in vivo in the Collagen–Capan-1 model. D indicates duct formation. Scale bar = 100 μm.
Figure 4. Internal histological analysis of 3D collagen PDAC model. (A) Representative images of H&E staining in Collagen–Capan-1 model on Day 7. (B) Representative images of H&E staining in Collagen–Capan-1 model on Day 14. D indicates duct formation. Scale bars = 200 μm, 200 μm, and 50 μm, respectively. (C) Morphology of Capan-1 cells in 2D culture. Scale bar, 200 μm. (D) Histological analysis of H&E-stained tumor cells generated subcutaneously in vivo in the Collagen–Capan-1 model. D indicates duct formation. Scale bar = 100 μm.
Bioengineering 10 01437 g004
Figure 5. Metabolic alterations to Capan-1 cells in 3D collagen matrix. (A) Different metabolites detected in the supernatants of Capan-1 cells in 2D culture and collagen 3D matrix. Metabolomic datasets were subjected to OPLS-DA analysis using MetaboAnalyst 5.0. The OPLS-DA score plot visually distinguishes between two sample groups (n = 6 samples per group) based on metabolomic data: the supernatant of Capan-1 cell in 2D culture is represented in red, while the supernatant of Capan-1 cells in collagen 3D matrix is represented in green. The model consists of one predictive x-score component: component t [1] and one orthogonal x-score component to [1]. t [1] explains 29.6% of the predictive variation in x, and to [1] explains 15.5% of the orthogonal variation in x. (B) Heatmap comparing metabolites changed in the supernatants of Capan-1 cells in 2D culture and collagen 3D matrix. (C) VIP plot corresponding to the score plot of OPLS-DA to visualize metabolite markers that contributed to the discrimination between the supernatants of Capan-1 cells in 2D culture and collagen 3D matrix. (D) Representative upregulated metabolites (proline, glutamic acid, and lactic acid) in the supernatants of Capan-1 cells cultured in collagen 3D matrix. (E) Representative downregulated metabolites (lysine, glutamine, valine) in the supernatants of Capan-1 cells in collagen 3D matrix culture. The metabolite concentration of Capan-1 collagen is represented by the red bar, while the metabolite concentration of 2D Capan-1 is represented by the blue bar.
Figure 5. Metabolic alterations to Capan-1 cells in 3D collagen matrix. (A) Different metabolites detected in the supernatants of Capan-1 cells in 2D culture and collagen 3D matrix. Metabolomic datasets were subjected to OPLS-DA analysis using MetaboAnalyst 5.0. The OPLS-DA score plot visually distinguishes between two sample groups (n = 6 samples per group) based on metabolomic data: the supernatant of Capan-1 cell in 2D culture is represented in red, while the supernatant of Capan-1 cells in collagen 3D matrix is represented in green. The model consists of one predictive x-score component: component t [1] and one orthogonal x-score component to [1]. t [1] explains 29.6% of the predictive variation in x, and to [1] explains 15.5% of the orthogonal variation in x. (B) Heatmap comparing metabolites changed in the supernatants of Capan-1 cells in 2D culture and collagen 3D matrix. (C) VIP plot corresponding to the score plot of OPLS-DA to visualize metabolite markers that contributed to the discrimination between the supernatants of Capan-1 cells in 2D culture and collagen 3D matrix. (D) Representative upregulated metabolites (proline, glutamic acid, and lactic acid) in the supernatants of Capan-1 cells cultured in collagen 3D matrix. (E) Representative downregulated metabolites (lysine, glutamine, valine) in the supernatants of Capan-1 cells in collagen 3D matrix culture. The metabolite concentration of Capan-1 collagen is represented by the red bar, while the metabolite concentration of 2D Capan-1 is represented by the blue bar.
Bioengineering 10 01437 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Song, X.; Nihashi, Y.; Yamamoto, M.; Setoyama, D.; Kunisaki, Y.; Kida, Y.S. Exploring the Role of Desmoplastic Physical Stroma in Pancreatic Cancer Progression Using a Three-Dimensional Collagen Matrix Model. Bioengineering 2023, 10, 1437. https://doi.org/10.3390/bioengineering10121437

AMA Style

Song X, Nihashi Y, Yamamoto M, Setoyama D, Kunisaki Y, Kida YS. Exploring the Role of Desmoplastic Physical Stroma in Pancreatic Cancer Progression Using a Three-Dimensional Collagen Matrix Model. Bioengineering. 2023; 10(12):1437. https://doi.org/10.3390/bioengineering10121437

Chicago/Turabian Style

Song, Xiaoyu, Yuma Nihashi, Masamichi Yamamoto, Daiki Setoyama, Yuya Kunisaki, and Yasuyuki S. Kida. 2023. "Exploring the Role of Desmoplastic Physical Stroma in Pancreatic Cancer Progression Using a Three-Dimensional Collagen Matrix Model" Bioengineering 10, no. 12: 1437. https://doi.org/10.3390/bioengineering10121437

APA Style

Song, X., Nihashi, Y., Yamamoto, M., Setoyama, D., Kunisaki, Y., & Kida, Y. S. (2023). Exploring the Role of Desmoplastic Physical Stroma in Pancreatic Cancer Progression Using a Three-Dimensional Collagen Matrix Model. Bioengineering, 10(12), 1437. https://doi.org/10.3390/bioengineering10121437

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop