Next Article in Journal
Multi-Modal Multi-Spectral Intravital Macroscopic Imaging of Signaling Dynamics in Real Time during Tumor–Immune Interactions
Next Article in Special Issue
Human Blood Vessel Organoids Penetrate Human Cerebral Organoids and Form a Vessel-Like System
Previous Article in Journal
CLTA-4 Expression Is Associated with the Maintenance of Chronic Inflammation in Endometriosis and Infertility
Previous Article in Special Issue
Establishment of Human Leukocyte Antigen-Mismatched Immune Responses after Transplantation of Human Liver Bud in Humanized Mouse Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Organoids Are Limited in Modeling the Colon Adenoma–Carcinoma Sequence

1
Department of Surgical Oncology, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
2
Department of Surgical Oncology, Graduate School of Medicine, Gifu University, 1-1 Yanagido, Gifu 501-1194, Japan
3
Department of Gastroenterological Surgery, Yokohama City University Graduate School of Medicine, Yokohama 236-0004, Japan
4
Department of Biostatistics & Bioinformatics, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
5
Department of Surgery, Niigata University Graduate School of Medical and Dental Sciences, Niigata 951-8510, Japan
6
Department of Surgery, University at Buffalo Jacobs School of Medicine and Biomedical Sciences, The State University of New York, Buffalo, NY 14263, USA
7
Department of Breast Oncology and Surgery, Tokyo Medical University, 6-7-1 Nishishinjuku, Shinjuku, Tokyo 160-8402, Japan
8
Department of Breast Surgery, Fukushima Medical University School of Medicine, Fukushima 960-1295, Japan
*
Author to whom correspondence should be addressed.
Cells 2021, 10(3), 488; https://doi.org/10.3390/cells10030488
Submission received: 31 December 2020 / Revised: 19 February 2021 / Accepted: 21 February 2021 / Published: 25 February 2021
(This article belongs to the Special Issue Current Progress in Organ Regeneration: Cells, Organoids and Organs)

Abstract

:
The colon adenoma–carcinoma sequence is a multistep genomic-altering process that occurs during colorectal cancer (CRC) carcinogenesis. Organoids are now commonly used to model both non-cancerous and cancerous tissue. This study aims to investigate how well organoids mimic tissues in the adenoma–carcinoma sequence by comparing their transcriptomes. A total of 234 tissue samples (48 adenomas and 186 CRC) and 60 organoid samples (15 adenomas and 45 CRC) were analyzed. We found that cell-proliferation-related gene sets were consistently enriched in both CRC tissues and organoids compared to adenoma tissues and organoids by gene set enrichment analysis (GSEA). None of the known pathways in the colon adenoma–carcinoma sequence were consistently enriched in CRC organoids. There was no enrichment of the tumor microenvironment-related gene sets in CRC organoids. CRC tissues enriched immune-response-related gene sets, whereas CRC organoids did not. The proportions of infiltrating immune cells were different between tissues and organoids, whereas there was no difference between cancer and adenoma organoids. The amounts of cancer stem cells and progenitor cells were not different between CRC and adenoma organoids, whereas a difference was noted between CRC and adenoma tissues. In conclusion, we demonstrated that organoids model only part of the adenoma–carcinoma sequence and should be used with caution after considering their limitations.

1. Introduction

Colorectal cancer (CRC) is the third most diagnosed cancer worldwide [1]. Despite significant advancement in multimodal treatment options during recent years, CRC remains the second leading cause of cancer-related deaths, trailing only lung cancer [1,2]. To this end, there remains a need for novel therapeutic agents to treat this disease.
Pre-clinical models using cancer cell lines have been the main vehicles for drug development [3]. The advantage of this model is in the simplicity and reproducibility its offers in testing and screening drugs. However, the monoclonal cancer cells used in this model are unable to accurately reflect cancer biology or the tumor microenvironment (TME). In order to overcome this limitation, animal models are used [4,5,6]. However, syngeneic models only provide data on rodent cancer [5,6], while xenograft models, including patient-derived xenografts (PDX), ignore immune response [7]. The lack of an appropriate system to model human cancer is the reasons many phase III clinical trials fail [8]. Recently, cancer organoids have been used in pre-clinical studies for their superior efficiency in performing drug sensitivity tests as compared to PDX models [4]. PDX models can take 6–8 months before an investigation can be completed, whereas tumor organoids only require 1–3 months [4,9]. It is thought that cancer organoids mimic human cancer biology with higher fidelity than cancer cell lines, but it remains unclear how closely cancer organoids resemble human cancers and their TME [4].
The majority of CRCs are caused by progressive genomic alterations referred to as the adenoma–carcinoma sequence [10,11]. This sequence is usually initiated with a genetic alteration of the APC gene that allows formation on an adenoma, which is subsequently followed by the alterations of the KRAS and TP53 genes [11]. The dysregulation of the pathways related to cell proliferation and apoptosis, including the WNT, RAS-MAPK, PI3K, TGF-β, and TP53 pathways, progresses to the carcinogenesis of CRC [12]. Matano et al. utilized organoid models to recapitulate this sequence and analyzed the significance of each pathway involved [12]. Although the morphology, mutational status, and drug response of cancer organoids have been extensively studied [13,14,15], comprehensive transcriptomic analyses that compare the adenoma–carcinoma sequence in tissues and organoids has not been performed. To this end, we aimed to clarify whether organoid models mimic the adenoma–carcinoma sequence in fresh-frozen tissue samples obtained from two cohorts from the Gene Expression Omnibus (GEO) database—GSE41258 (adenoma and CRC tissue cohort) [16] and GSE57965 (adenoma and CRC organoid cohort) [12].
Recently, our group has performed in silico translational research to identify biomarkers [17,18,19], clinically relevant immune cells [20,21,22,23], predictive genes [24], as well as microRNAs in breast and gastrointestinal cancers [25]. We employed a computational algorithm, referred to as gene set enrichment analysis (GSEA), which enables us to analyze the differences in biological pathways between two distinct groups [26]. We utilized CIBERSORTx [27] and xCell [28], which uses transcriptomic data to estimate the composition of immune cells within human tumors. In this study, we hypothesized that mechanisms of the adenoma–carcinoma sequence are maintained in organoid models, and sought to clarify the utility of organoid models for future investigations.

2. Materials and Methods

2.1. Data Acquisition of Colorectal Cancer Cohorts

The transcriptomic and clinical data for colorectal cancer (CRC) and adenoma were obtained from three cohorts—one cohort of tissues (GSE41258) [16] and two cohorts of organoids (GSE57965 [12] and GSE74843) [29]. GSE41258 holds a total of 390 tissue samples, which include 186 primary CRCs, 48 adenomas (polyps), 1 high grade polyp, 47 liver metastases, 20 lung metastases, 2 microadenomas, 54 normal colons, 13 normal livers, 7 normal lungs, and 12 cell lines. Among these, we utilized a total of 234 (48 adenoma and 186 primary CRC) tissue samples. GSE57965 holds a total of 23 samples, which include 5 colon adenoma organoids, 3 colon epithelial organoids, 7 CRC organoids, 6 genome-engineered organoids, and 2 metastatic CRC organoids. Among these, we utilized a total of 12 (5 adenoma organoids and 7 CRC organoids) samples in the current study. GSE57965 holds a total of 58 samples, which include 10 adenoma organoids, 38 CRC organoids, 7 normal colorectal mucosa organoids, and 3 serrated organoids. Among these, we utilized a total of 48 samples (10 adenoma organoids and 38 CRC organoids). In an attempt to verify our results, we identified two other studies that analyzed the transcriptome of organoids. The study by van de Wetering et al. (GSE64392) [14] includes only CRC organoids and normal tissue organoids, but does not include adenoma organoids. Thus, it was unusable for our analyses. The study data in the study by Weeber et al. [13] was not deposited to the Gene Expression Omnibus, and therefore we were unable to access it. Given that all of the cohorts used in the current study are from a publicly accessible and deidentified database, ethics approval was waived by the Institutional Review Board.

2.2. Gene Set Enrichment Analysis (GSEA) and Single-Sample GSEA

The Broad Institute provided the publicly available software, gene set enrichment analysis (GSEA) (http://software.broadinstitute.org/gsea/index.jsp) [26]. The Broad Institute was founded by Massachusetts Institute of Technology and Harvard University (https://www.broadinstitute.org/) in 2004. In the current study, we utilized hallmark gene sets of Molecular Signature Database collections, as we have reported previously [17,18,25,30]. The statistical significance of the false discovery rate (FDR) was set to 0.25, as recommended by the GSEA developer Broad Institute. We also performed single-sample GSEA to demonstrate the absolute enrichment value of each sample by calculating an enrichment score for the gene set of interest [31].

2.3. CIBERSORTx

CIBERSORTx, the algorithm that enabled us to analyze the proportions of 22 immune cells within the tumor microenvironment of tissue samples, was developed and reported by Newman et al. [27].

2.4. Statistical Analysis

To perform the statistical analysis, we utilized the publicly available software R (version 4.02). The statistical significance was analyzed using Fisher’s exact test or one-way ANOVA. For the box plot figure, one-way ANOVA was performed. Median values are represented as black lines inside boxplots (Tukey type) and the spans of rectangles demonstrate interquartile ranges. The statistical significance was defined as a two-sided p value of less than 0.05.

3. Results

3.1. Cell-Proliferation-Related Gene Sets Were Enriched in Both Colorectal Cancer (CRC) Tissue and Organoid Cohorts

Cancer is more proliferative than benign adenoma tumors given the nature of malignancy. In order to investigate whether cancer organoids are more proliferative than adenoma organoids, we conducted gene set enrichment analysis (GSEA) of hallmark cell-proliferation-related gene sets on cohorts that included both adenoma and CRC transcriptomes from tissue (GSE41258, n = 48 and n = 186, respectively) and from organoids (GSE57965, n = 5 and n = 7, respectively). As expected, CRC tissues significantly enriched hallmark cell-proliferation-related gene sets, such as E2F target, G2M checkpoint, and mitotic spindle, but not Myc targets V1 or V2 (Figure 1A, Table S1; normalized enrichment score (NES) and false discovery rate (FDR); NES = −1.47, FDR = 0.133; NES = −1.56, FDR = 0.121; NES = −1.60, FDR = 0.167; NES = −0.94, FDR = 0.560; NES = −1.09, FDR = 0.439, respectively). MKI67, the most frequently used cell proliferation marker in the clinical setting, was significantly elevated in CRC tissue compared to adenoma tissue (Figure 1B, p = 0.008).
Interestingly, CRC organoids enriched hallmark cell-proliferation-related gene sets more than CRC tissue. CRC organoids enriched E2F targets, G2M checkpoint, Myc targets V1 and V2, but not mitotic spindle (Figure 1C, Table S2; NES = −1.58, FDR = 0.100; NES = −1.48, FDR = 0.132; NES = −1.52, FDR = 0.118; NES = −1.91, FDR = 0.002; NES = −1.09, FDR = 0.493, respectively). On the other hand, the expression levels of MKI67 were not statistically different between CRC and adenoma organoids, despite a difference in median value, likely because of the large variance and small sample size of CRC organoids (Figure 1D, p = 0.187). To overcome the sample size limitation, we performed a single-sample GSEA from the GSE57965 cohort. The results of single-sample GSEA were strikingly similar to that of larger cohort GSEA. The gene sets that were enriched in CRC tissue, such as E2F targets, G2M checkpoint, and Myc targets V1 and V2, were significantly enriched in CRC organoids as well. In addition, CRC organoids had a higher GSEA score in the single-sample GSEA than adenoma organoids (Figure 1E; p = 0.024, p = 0.037, p = 0.036, and p < 0.001 respectively). We additionally analyzed the GSE74843 cohort, which included adenoma organoids and CRC organoids. In this cohort, CRC organoids enriched the Myc target V2 gene set, which was further validated in the single-sample GSEA. CRC organoids demonstrated a higher single-sample GSEA score than adenoma organoids in the Myc target V2 gene set (Supplementary Figure S1).
The noted enrichment of cell-proliferation-related gene sets suggests that cancer organoids may be as proliferative as adenoma organoids and that CRC tissue may be as proliferative as adenoma tissue.

3.2. Among the Pathways in the Colon Adenoma–Carcinoma Sequence, Only the MTORC1 Gene Set Was Enriched in CRC Organoids

Several pathways are known to contribute to the development of CRC from a benign adenoma, also known as colon adenoma–carcinoma sequence [11,12]. It was of interest to investigate whether the same pathways were enriched in tumor tissues and organoids. GSEA of those pathways was conducted on tissue (GSE41258) and organoid (GSE57965) cohorts. In the tissue cohort, CRC tissue enriched all of the pathways involved in the adenoma–carcinoma sequence, including WNT β-catenin, KRAS signaling, MTORC1, and TGF-β pathways (Figure 2A, Table S1; NES = −1.36 FDR = 0.184; NES = −1.57 FDR = 0.136; NES = −1.58 FDR = 0.150; NES = −1.30 FDR = 0.234, respectively). On the other hand, only the MTORC1 pathway was enriched in CRC organoids in the GSE57965 cohort (Figure 2B, Table S2; NES = −1.38 FDR = 0.197). The TGF-β pathway was enriched in adenoma organoids (Figure 2B, NES = −1.68, FDR = 0.152), which was consistent with the GSE74843 cohort (Figure 2C; NES = −1.77 FDR = 0.182).
These results suggest that CRC organoids do not enrich as many of the adenoma–carcinoma sequence-related gene sets as CRC tissues.

3.3. Only CRC Tissues Enriched Tumor Immune Microenvironment (TME)-Related Gene Sets and Correlated with a Higher Infiltration of Stromal Cells When Compared to Adenoma Tissues

Given the nature of organoids, it was of interest to investigate whether they maintained the signaling pathways and contained the infiltrating cells of the tissue TME. CRC tissue enriched TME-related gene sets, including angiogenesis, but not adipogenesis (Figure 3A, Table S1; NES = −1.93 FDR = 0.024; NES = −0.86 FDR = 0.659, respectively). In contrast, CRC organoids did not enrich either of these gene sets (Figure 3B). CRC tissue analysis correlated with a higher infiltration of stromal cells in the TME as compared to adenoma tissues, including fibroblasts, lymphatic endothelial cells, microvascular endothelial cells, and pericytes (Figure 3C; p < 0.001, p < 0.001, p < 0.001, and p = 0.031). In addition, the stroma score was significantly higher in CRC tissues (Figure 3C; p < 0.001). In contrast, stromal cells were not highly infiltrated in CRC organoids and the corresponding stroma score was low (Figure 3D). In adenoma organoids from the GSE74843 cohort, only fibroblasts were highly infiltrated and the stroma score was significantly elevated (Supplementary Figure S2). These results imply that both adenoma and CRC organoids have poorly developed TMEs, meaning there is no difference in cell infiltration or TME-related gene expression profiles.

3.4. CRC Tissues Enriched Immune-Response-Related Gene Sets but CRC Organoids Did Not

It was previously demonstrated that highly proliferative cancer is associated with an enhanced immune response in the TME [32]. As expected, CRC tissue enriched immune-response-related gene sets, such as allograft rejection, interferon (IFN)-α response, IFN-γ response, inflammatory response, IL-6/JAK/STAT3, complement, and coagulation (Figure 4A, Table S1; NES = −1.43 FDR = 0.145; NES = −1.42 FDR = 0.150; NES = −1.48 FDR = 0.135; NES = −1.55 FDR = 0.119; NES = −1.44 FDR = 0.144; NES = −1.46 FDR = 0.138; NES = −1.52 FDR = 0.130, respectively). CRC organoids and adenoma organoids did not enrich any of the immune-response-related gene sets in either cohort (GSE57965 (Figure 4B) or GSE74843 (Supplementary Figure S3)). These results suggest that immune activity in CRC organoids is not different from immune activity in adenoma organoids.

3.5. The Proportions of Infiltrating Immune Cells Are Different between Tissue and Organoid Cohorts

Utilizing the CIBERSORT algorithm, the cell compositions of immune cell types were estimated from each transcriptomic profile. In adenoma tissues, nearly 60% of the cell population was composed of lymphocytes, half of which were B cells (Figure 5A). Among the myeloid cells, monocytes composed 29.8% of the population. In CRC tissues, nearly 75% of cells were lymphocytes. B cells were the most abundant overall cell type (Figure 5A) and granulocytes were the most abundant myeloid cell type in CRC tissues.
In adenoma organoids, the proportion of lymphocytes was similar to that in tissue samples, but the most abundant cell type was T cells (Figure 5B). Of the myeloid cells, granulocytes were the most abundant cell type in adenoma organoids. A similar trend was observed in CRC organoids, which contained a large proportion of lymphocytes, but the most abundant cell type was T cells. Regarding myeloid cells, the population of granulocytes was abundant in CRC organoids as well.

3.6. B Cells and T Cells Were Less Infiltrated and Macrophages Were More Infiltrated in CRC Tissues Compared to Adenoma Tissues, Whereas There Were No Differences in Organoids

The proportional differences of immune cells prompted us to analyze the cell composition differences of lymphocytes and myeloid cells between adenomas, CRC tissues, and organoids. Among T cells, infiltration of CD4 memory resting T cells and regulatory T cells (Tregs) was low, whereas follicular T cells were higher in CRC tissue (Figure 6A; p < 0.001, p = 0.0018, and p = 0.002, respectively). Only CD4 memory resting T cells were significantly lower in CRC organoids compared to adenoma organoids in both GSE57965 (Figure 6A; p = 0.012) and GSE74843 cohorts (Supplementary Figure S4A; p = 0.046). CD4-memory-activated cells were significantly higher in only the GSE57965 cohort (Figure 6A; p = 0.049). All types of macrophages (M0, M1, and M2) were highly infiltrated, whereas activated dendritic cells (DC) were less infiltrated in CRC tissue (Figure 6B; macrophages all p < 0.001, DC p = 0.046). CRC organoids demonstrated consistently higher infiltration of M0 cells in the GSE74843 cohort, whereas activated DCs (Supplementary Figure S4B; p = 0.015) were consistently infiltrated in CRC organoids from the GSE57965 cohort (Figure 6B; p = 0.005). The infiltration of memory B cells and plasma cells was low in CRC tissue compared to adenoma tissue, but there was no difference in B cell infiltration between CRC and adenoma organoids (Supplementary Figure S5A). There was no difference in infiltration of NK cells and eosinophils between CRC tissue, adenoma tissue, or organoids (Supplementary Figure S5B,C). Resting mast cells demonstrated higher infiltration in CRC organoids, but activated mast cells demonstrated higher infiltration in adenoma organoids from the GSE74843 cohort (Supplementary Figure S5C). Neutrophils were more prevalent in CRC tissues, but not in CRC organoids (Supplementary Figure S5).

3.7. There Was No Difference in the Infiltration of Stem Cells and Progenitor Cells between Adenoma and CRC Organoids

Given the nature of malignancy, we expected CRC tissues to have high infiltration of mesenchymal stem cells and progenitor cells. As expected, CRC tissues were associated with high infiltration of mesenchymal stem cells (Figure 7A; p < 0.001), but there was no difference in hematopoietic stem cell infiltration. In contrast, CRC organoids demonstrated no difference in either type of stem cell when compared to adenoma organoids in both GSE57965 and GSE74843 cohorts (Figure 7B, Supplementary Figure S6A). Unexpectedly, CRC tissues had less infiltration of multipotent progenitor cells and common myeloid progenitor cells, but had greater infiltration of megakaryocyte–erythroid progenitor cells (Figure 7C; p < 0.001, p < 0.001, and p < 0.001, respectively). Only megakaryocyte–erythroid progenitor cells were highly infiltrated in CRC organoids from the GSE74843 cohort (Supplementary Figure S6B). There was no difference in progenitor cells between CRC and adenoma organoids in the GSE57965 cohort (Figure 7D). These results imply that both CRC and adenoma organoids possess stem-cell-like characteristics, and thus there is no difference in stemness between the organoid types.

4. Discussion

In the current study, we found that cell-proliferation-related gene sets were enriched in both colorectal cancer (CRC) tissues and organoids. Only the MTORC1 pathway, among all of the pathways involved in the adenoma–carcinoma sequence (WNT β-catenin, KRAS signaling up, MTORC1, and TGF-β pathways), was enriched in CRC organoids of the GSE 57965 cohort. In comparison, all of the pathways were enriched in CRC tissue. CRC tissue enriched angiogenesis and infiltrated stromal cell gene sets, represented by fibroblast and vascular endothelial cell infiltration, whereas CRC organoids did not. CRC tissue enriched immune-response-related gene sets, which was not seen for CRC organoids. The ratio of immune cells was different between tissues and organoids. Lymphocytes were more abundant in CRC tissue compared to adenoma tissue. B cells and T cells were less infiltrated and macrophages were highly infiltrated in CRC tissue. There was no difference in lymphocyte or myeloid cell infiltration between CRC and adenoma organoids. There were no differences in cancer stem cells or progenitor cells between CRC and adenoma organoids. There was a difference in cancer stem and progenitor cells between CRC and adenoma tissue, which implies that organoids are inherently rich in these cells regardless of their derivation from cancer or adenoma.
Organoids are an important part of drug development and drug sensitivity screening, especially from the perspective of personalized medicine, given the ease with which they can be generated from patient tumors and the fidelity of their genomic profile to the original tumor. Organoids are more costly in comparison to cell line models, however they are reported to be more economical in comparison to PDX models [4]. As an advantage, organoids require less time to complete drug response assays compared to PDX models [9,14]. An expediated time frame could increase transability and support early clinical availability of drugs for clinical use.
CRC organoids enriched the cell-proliferation-related gene sets, which was concordant with CRC tissues. Multiple studies have investigated the gene expression alterations between tumors and their derivative organoids. Matano et al. reported that cancer organoids maintained the aggressive biology of cancer tissues by expressing CRC-specific genes [12]. Fujii et al. further demonstrated that CRC organoids preserved the same gene signatures and gene alterations of the original CRC tissue, as assessed by GSEA [29]. In agreement, Weeber et al. reported that cancer organoids demonstrated similar genetic profiles compared to the original tumor [13]. Furthermore, van de Wetering et al. reported the establishment of an organoid bank in which organoids possessed the genetic alterations representative of CRC, such as alterations in APC, TP53, PIK3CA, and KRAS [14]. In the current study, the transcriptome of the organoids was not compared to the tumor it was derived from. The original study by Matano et al. compared adenoma tissues and organoids to CRC tissues and organoids using GSEA to analyze the genomic fidelity of derived organoids to the original tumors [12]. The novelty of our approach is that we were able to study cancer biology as represented by the hallmarks of cancer progression (cell proliferation, inflammation, and metabolism) through the GSEA of hallmark gene sets of the Molecular Signature Database. We found that CRC organoids enriched the cell-proliferation-related gene sets more than adenoma organoids. The same relationship is present in the comparison between CRC tissues and adenoma tissues.
The mutational status and drug responses of organoids have been extensively studied [13,14,15]. However, this is the first study that has analyzed the biology of the adenoma–carcinoma sequence by comparing the transcriptomics of CRC and adenoma tissues and their respective organoids. We had expected that CRC organoids would enrich the gene sets related to the adenoma–carcinoma sequence. To our surprise, our analysis demonstrated the enrichment of only the MTORC1 pathway. These results may be due to the fact that GSEA is analyzing multiple genes instead of one certain gene.
A limitation of organoids is that they do not proportionally reflect the composition of the tumor microenvironment as represented by tissue-resident immune cells or vasculature promoting cells of the original cancer tissue [29,33]. The novelty of the current study is that it investigates the differences in stromal and vasculature formation between CRC organoids and adenoma organoids. The results indicated that there was lack of fidelity in the proportional loss of tissue-resident and vasculature-promoting cells between CRC and adenoma organoids versus tissues. This implies that immune and vascular cells are lost during the organoid development process, and thus no similarities were maintained.
In the current study, we found that the tumor immune microenvironment of the organoids is very different from that of tissues, and it does not reflect the changes of the adenoma–carcinoma sequence. This result suggests that interpretation of experimental results using organoids should be done with caution when analyzing immune cell functions. To overcome this issue, Neal et al. proposed a co-cultured system of cancer organoids and endogenous tumor infiltrating lymphocytes to test the drug response of anti-PD-1 or anti-PD-L1 immune check point inhibitors [33]. Outside of this artificial reconstitution of the TME, the complex interaction between organoids and immune cells cannot be reliably interpreted.
The infiltration of stem cells and progenitor cells, which are cells of multipotency and undifferentiation [34], did not differ between CRC and adenoma organoids, but differed between CRC and adenoma tissues. These results may reflect the ability of organoids to maintain stem cell characteristics, regardless of the tissue of origin (adenoma or cancer) [35,36]. Whether this stem-cell-rich environment of the organoid system interferes with its ability to model the original tumor tissue is to be determined with further investigation.
The current study has obvious limitations. One of the limitations stems from our inability to compare tissue samples and organoid samples in the same population. This was due to the lack of access to a cohort that has both tissue and organoid transcriptomic data. In addition, our results may not reflect the spatial difference within the bulk tumor and the heterogeneity among the patients due to the small number of organoid samples. Although we believe that the statistical significance is real when the difference persists despite the small sample size, we may not be encompassing many findings due to the small sample size of the organoids. Further investigation comparing organoid and tissue samples directly is warranted and would make our findings more intriguing and convincing.
In conclusion, we demonstrated that the organoids mimic only part of the adenoma–carcinoma sequence of colorectal cancer development. Thus, one needs to use organoids with caution due to their limitations.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/10/3/488/s1, Figure S1: Gene set enrichment analysis (GSEA) of cell-proliferation-related gene sets and analysis of MKI67 expression. Figure S2: GSEA of tumor immune microenvironment (TME)-related gene sets, the infiltration of stromal cells, and comparison of stroma scores. (A) GSEA of tissue sample (GSE74843). Figure S3: GSEA of gene sets associated with immune response. Analysis of tissue sample (GSE74843). Figure S4: Comparison of the infiltration of immune cells of GSE74843. Figure S5: Comparison of the infiltration of immune cells (B cells, NK cells, and granulocytes). Figure S6: Comparison of infiltration of stem cells and progenitor cells between adenomas and CRC tissues and organoids. Table S1: Gene sets enriched in colorectal cancer tissue (GSE41258). Table S2: Gene sets enriched in colorectal cancer organoids (GSE57965).

Author Contributions

Conceptualization, Y.T.; methodology, M.O.; formal analysis, Y.T. and W.T.; resources, L.Y., N.M., M.F., and K.Y.; writing—original draft preparation, Y.T.; writing—review and editing, A.P., and K.T.; supervision, K.T.; project administration, K.T. All authors have read and agreed to the published version of the manuscript.

Funding

Kazuaki Takabe, the corresponding author, is supported by Paul & Helen Ellis Charitable Trust US, as well as National Cancer Institute (NCI)/ National Institutes of Health (NIH) (R01CA160688, R01CA250412, R37CA248018, and Department of Defense - BCRP grant W81XWH-19-1-0674. Also, Roswell Park Comprehensive Cancer Center is supported by NCI/NIH (P30CA016056).

Institutional Review Board Statement

Given that all of the cohorts used in this study were from a publicly accessible and deidentified database, ethics approval was waived by the Institutional Review Board for the current study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  2. Joharatnam-Hogan, N.; Wilson, W.; Shiu, K.K.; Fusai, G.K.; Davidson, B.; Hochhauser, D.; Bridgewater, J.; Khan, K. Multimodal Treatment in Metastatic Colorectal Cancer (mCRC) Improves Outcomes-The University College London Hospital (UCLH) Experience. Cancers 2020, 12, 3545. [Google Scholar] [CrossRef]
  3. Lau, H.C.H.; Kranenburg, O.; Xiao, H.; Yu, J. Organoid models of gastrointestinal cancers in basic and translational research. Nat. Rev. Gastroenterol. Hepatol. 2020, 17, 203–222. [Google Scholar] [CrossRef]
  4. Weeber, F.; Ooft, S.N.; Dijkstra, K.K.; Voest, E.E. Tumor Organoids as a Pre-clinical Cancer Model for Drug Discovery. Cell Chem. Biol. 2017, 24, 1092–1100. [Google Scholar] [CrossRef] [PubMed]
  5. Rashid, O.M.; Nagahashi, M.; Ramachandran, S.; Dumur, C.; Schaum, J.; Yamada, A.; Terracina, K.P.; Milstien, S.; Spiegel, S.; Takabe, K. An improved syngeneic orthotopic murine model of human breast cancer progression. Breast Cancer Res. Treat. 2014, 147, 501–512. [Google Scholar] [CrossRef] [Green Version]
  6. Katsuta, E.; Rashid, O.M.; Takabe, K. Murine breast cancer mastectomy model that predicts patient outcomes for drug development. J. Surg. Res. 2017, 219, 310–318. [Google Scholar] [CrossRef] [PubMed]
  7. Kawaguchi, T.; Foster, B.A.; Young, J.; Takabe, K. Current Update of Patient-Derived Xenograft Model for Translational Breast Cancer Research. J. Mammary Gland Biol. Neoplasia 2017, 22, 131–139. [Google Scholar] [CrossRef] [PubMed]
  8. McMillin, D.W.; Negri, J.M.; Mitsiades, C.S. The role of tumour-stromal interactions in modifying drug response: Challenges and opportunities. Nat. Rev. Drug Discov. 2013, 12, 217–228. [Google Scholar] [CrossRef] [PubMed]
  9. Hidalgo, M.; Bruckheimer, E.; Rajeshkumar, N.V.; Garrido-Laguna, I.; De Oliveira, E.; Rubio-Viqueira, B.; Strawn, S.; Wick, M.J.; Martell, J.; Sidransky, D. A pilot clinical study of treatment guided by personalized tumorgrafts in patients with advanced cancer. Mol. Cancer Ther. 2011, 10, 1311–1316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Morson, B.C. Evolution of cancer of the colon and rectum. Cancer 1974, 34, 845–849. [Google Scholar] [CrossRef]
  11. Cross, W.; Kovac, M.; Mustonen, V.; Temko, D.; Davis, H.; Baker, A.M.; Biswas, S.; Arnold, R.; Chegwidden, L.; Gatenbee, C.; et al. The evolutionary landscape of colorectal tumorigenesis. Nat. Ecol. Evol. 2018, 2, 1661–1672. [Google Scholar] [CrossRef] [Green Version]
  12. Matano, M.; Date, S.; Shimokawa, M.; Takano, A.; Fujii, M.; Ohta, Y.; Watanabe, T.; Kanai, T.; Sato, T. Modeling colorectal cancer using CRISPR-Cas9-mediated engineering of human intestinal organoids. Nat. Med. 2015, 21, 256–262. [Google Scholar] [CrossRef] [PubMed]
  13. Weeber, F.; van de Wetering, M.; Hoogstraat, M.; Dijkstra, K.K.; Krijgsman, O.; Kuilman, T.; Gadellaa-van Hooijdonk, C.G.; van der Velden, D.L.; Peeper, D.S.; Cuppen, E.P.; et al. Preserved genetic diversity in organoids cultured from biopsies of human colorectal cancer metastases. Proc. Natl. Acad. Sci. USA 2015, 112, 13308–13311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Van de Wetering, M.; Francies, H.E.; Francis, J.M.; Bounova, G.; Iorio, F.; Pronk, A.; van Houdt, W.; van Gorp, J.; Taylor-Weiner, A.; Kester, L.; et al. Prospective derivation of a living organoid biobank of colorectal cancer patients. Cell 2015, 161, 933–945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Pauli, C.; Hopkins, B.D.; Prandi, D.; Shaw, R.; Fedrizzi, T.; Sboner, A.; Sailer, V.; Augello, M.; Puca, L.; Rosati, R.; et al. Personalized In Vitro and In Vivo Cancer Models to Guide Precision Medicine. Cancer Discov. 2017, 7, 462–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Sheffer, M.; Bacolod, M.D.; Zuk, O.; Giardina, S.F.; Pincas, H.; Barany, F.; Paty, P.B.; Gerald, W.L.; Notterman, D.A.; Domany, E. Association of survival and disease progression with chromosomal instability: A genomic exploration of colorectal cancer. Proc. Natl. Acad. Sci. USA 2009, 106, 7131–7136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Oshi, M.; Katsuta, E.; Yan, L.; Ebos, J.M.L.; Rashid, O.M.; Matsuyama, R.; Endo, I.; Takabe, K. A Novel 4-Gene Score to Predict Survival, Distant Metastasis and Response to Neoadjuvant Therapy in Breast Cancer. Cancers 2020, 12, 1148. [Google Scholar] [CrossRef]
  18. Oshi, M.; Tokumaru, Y.; Asaoka, M.; Yan, L.; Satyananda, V.; Matsuyama, R.; Matsuhashi, N.; Futamura, M.; Ishikawa, T.; Yoshida, K.; et al. M1 Macrophage and M1/M2 ratio defined by transcriptomic signatures resemble only part of their conventional clinical characteristics in breast cancer. Sci. Rep. 2020, 10, 16554. [Google Scholar] [CrossRef] [PubMed]
  19. Oshi, M.; Tokumaru, Y.; Angarita, F.A.; Yan, L.; Matsuyama, R.; Endo, I.; Takabe, K. Degree of Early Estrogen Response Predict Survival after Endocrine Therapy in Primary and Metastatic ER-Positive Breast Cancer. Cancers 2020, 12, 3557. [Google Scholar] [CrossRef]
  20. Oshi, M.; Asaoka, M.; Tokumaru, Y.; Yan, L.; Matsuyama, R.; Ishikawa, T.; Endo, I.; Takabe, K. CD8 T Cell Score as a Prognostic Biomarker for Triple Negative Breast Cancer. Int. J. Mol. Sci. 2020, 21, 6968. [Google Scholar] [CrossRef]
  21. Tokumaru, Y.; Oshi, M.; Katsuta, E.; Yan, L.; Huang, J.L.; Nagahashi, M.; Matsuhashi, N.; Futamura, M.; Yoshida, K.; Takabe, K. Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation. Int. J. Mol. Sci. 2020, 21, 5744. [Google Scholar] [CrossRef] [PubMed]
  22. Oshi, M.; Newman, S.; Tokumaru, Y.; Yan, L.; Matsuyama, R.; Kalinski, P.; Endo, I.; Takabe, K. Plasmacytoid Dendritic Cell (pDC) Infiltration Correlate with Tumor Infiltrating Lymphocytes, Cancer Immunity, and Better Survival in Triple Negative Breast Cancer (TNBC) More Strongly than Conventional Dendritic Cell (cDC). Cancers 2020, 12, 3342. [Google Scholar] [CrossRef]
  23. Oshi, M.; Asaoka, M.; Tokumaru, Y.; Angarita, F.A.; Yan, L.; Matsuyama, R.; Zsiros, E.; Ishikawa, T.; Endo, I.; Takabe, K. Abundance of Regulatory T Cell (Treg) as a Predictive Biomarker for Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer. Cancers 2020, 12, 3038. [Google Scholar] [CrossRef] [PubMed]
  24. Oshi, M.; Newman, S.; Murthy, V.; Tokumaru, Y.; Yan, L.; Matsuyama, R.; Endo, I.; Takabe, K. ITPKC as a Prognostic and Predictive Biomarker of Neoadjuvant Chemotherapy for Triple Negative Breast Cancer. Cancers 2020, 12, 2758. [Google Scholar] [CrossRef] [PubMed]
  25. Tokumaru, Y.; Katsuta, E.; Oshi, M.; Sporn, J.C.; Yan, L.; Le, L.; Matsuhashi, N.; Futamura, M.; Akao, Y.; Yoshida, K.; et al. High Expression of miR-34a Associated with Less Aggressive Cancer Biology but Not with Survival in Breast Cancer. Int. J. Mol. Sci. 2020, 21, 3045. [Google Scholar] [CrossRef]
  26. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Newman, A.M.; Steen, C.B.; Liu, C.L.; Gentles, A.J.; Chaudhuri, A.A.; Scherer, F.; Khodadoust, M.S.; Esfahani, M.S.; Luca, B.A.; Steiner, D.; et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 2019, 37, 773–782. [Google Scholar] [CrossRef] [PubMed]
  28. Aran, D.; Hu, Z.; Butte, A.J. xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017, 18, 220. [Google Scholar] [CrossRef] [Green Version]
  29. Fujii, M.; Shimokawa, M.; Date, S.; Takano, A.; Matano, M.; Nanki, K.; Ohta, Y.; Toshimitsu, K.; Nakazato, Y.; Kawasaki, K.; et al. A Colorectal Tumor Organoid Library Demonstrates Progressive Loss of Niche Factor Requirements during Tumorigenesis. Cell Stem Cell 2016, 18, 827–838. [Google Scholar] [CrossRef] [Green Version]
  30. Takeshita, T.; Yan, L.; Asaoka, M.; Rashid, O.; Takabe, K. Late recurrence of breast cancer is associated with pro-cancerous immune microenvironment in the primary tumor. Sci. Rep. 2019, 9, 16942. [Google Scholar] [CrossRef] [Green Version]
  31. Barbie, D.A.; Tamayo, P.; Boehm, J.S.; Kim, S.Y.; Moody, S.E.; Dunn, I.F.; Schinzel, A.C.; Sandy, P.; Meylan, E.; Scholl, C.; et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature 2009, 462, 108–112. [Google Scholar] [CrossRef] [PubMed]
  32. Oshi, M.; Takahashi, H.; Tokumaru, Y.; Yan, L.; Rashid, O.M.; Matsuyama, R.; Endo, I.; Takabe, K. G2M Cell Cycle Pathway Score as a Prognostic Biomarker of Metastasis in Estrogen Receptor (ER)-Positive Breast Cancer. Int. J. Mol. Sci. 2020, 21, 2921. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Neal, J.T.; Li, X.; Zhu, J.; Giangarra, V.; Grzeskowiak, C.L.; Ju, J.; Liu, I.H.; Chiou, S.H.; Salahudeen, A.A.; Smith, A.R.; et al. Organoid Modeling of the Tumor Immune Microenvironment. Cell 2018, 175, 1972–1988.e1916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Xu, H.; Lyu, X.; Yi, M.; Zhao, W.; Song, Y.; Wu, K. Organoid technology and applications in cancer research. J. Hematol. Oncol. 2018, 11, 116. [Google Scholar] [CrossRef] [PubMed]
  35. Bredenoord, A.L.; Clevers, H.; Knoblich, J.A. Human tissues in a dish: The research and ethical implications of organoid technology. Science 2017, 355. [Google Scholar] [CrossRef] [PubMed]
  36. Lancaster, M.A.; Knoblich, J.A. Organogenesis in a dish: Modeling development and disease using organoid technologies. Science 2014, 345, 1247125. [Google Scholar] [CrossRef]
Figure 1. Gene set enrichment analysis (GSEA) of cell-proliferation-related gene sets and analysis of MKI67 expression. (A) GSEA of adenoma vs. CRC tissue in GSE41258. (B) The MKI67 expression levels of adenoma and CRC in tissue. (C) GSEA of adenoma vs. CRC organoids in GSE57965. (D) The comparison of MKI67 expression levels between adenoma and CRC in organoids. Tukey-type boxplots demonstrate the median and interquartile level values. (E) Single-sample GSEA of GSE57965. Statistical significance was defined as a false discovery rate (FDR) < 0.25. Ade, adenoma; CRC, colorectal cancer; FDR, false discovery rate.
Figure 1. Gene set enrichment analysis (GSEA) of cell-proliferation-related gene sets and analysis of MKI67 expression. (A) GSEA of adenoma vs. CRC tissue in GSE41258. (B) The MKI67 expression levels of adenoma and CRC in tissue. (C) GSEA of adenoma vs. CRC organoids in GSE57965. (D) The comparison of MKI67 expression levels between adenoma and CRC in organoids. Tukey-type boxplots demonstrate the median and interquartile level values. (E) Single-sample GSEA of GSE57965. Statistical significance was defined as a false discovery rate (FDR) < 0.25. Ade, adenoma; CRC, colorectal cancer; FDR, false discovery rate.
Cells 10 00488 g001
Figure 2. GSEA of gene sets associated with an adenoma–carcinoma sequence. (A) Analysis of tissue sample (GSE41258). (B) Analysis of organoid sample (GSE57965). (C) Analysis of organoid sample (GSE74843). Statistical significance was defined as a false discovery rate (FDR) < 0.25. CRC, colorectal cancer; FDR, false discovery rate; NES, normalized enrichment score.
Figure 2. GSEA of gene sets associated with an adenoma–carcinoma sequence. (A) Analysis of tissue sample (GSE41258). (B) Analysis of organoid sample (GSE57965). (C) Analysis of organoid sample (GSE74843). Statistical significance was defined as a false discovery rate (FDR) < 0.25. CRC, colorectal cancer; FDR, false discovery rate; NES, normalized enrichment score.
Cells 10 00488 g002
Figure 3. GSEA of tumor immune microenvironment (TME)-related gene sets, the infiltration of stromal cells, and comparison of stroma scores. (A) GSEA of tissue sample (GSE41258). (B) GSEA of organoid sample (GSE57965). (C) Infiltration of stromal cells and comparison of stroma scores in tissue samples (GSE41258). (D) Infiltration of stromal cells and comparison of stroma scores in organoid samples (GSE57965). Tukey-type boxplots demonstrate the median and interquartile level values. Statistical significance was defined as a false discovery rate (FDR) < 0.25. Ade, adenoma; CRC, colorectal cancer. NES, normalized enrichment score.
Figure 3. GSEA of tumor immune microenvironment (TME)-related gene sets, the infiltration of stromal cells, and comparison of stroma scores. (A) GSEA of tissue sample (GSE41258). (B) GSEA of organoid sample (GSE57965). (C) Infiltration of stromal cells and comparison of stroma scores in tissue samples (GSE41258). (D) Infiltration of stromal cells and comparison of stroma scores in organoid samples (GSE57965). Tukey-type boxplots demonstrate the median and interquartile level values. Statistical significance was defined as a false discovery rate (FDR) < 0.25. Ade, adenoma; CRC, colorectal cancer. NES, normalized enrichment score.
Cells 10 00488 g003
Figure 4. GSEA of gene sets associated with immune response. (A) Analysis of tissue sample (GSE41258). (B) Analysis of organoid sample (GSE57965). Statistical significance was defined as a false discovery rate (FDR) < 0.25. CRC, colorectal cancer; FDR, false discovery rate.
Figure 4. GSEA of gene sets associated with immune response. (A) Analysis of tissue sample (GSE41258). (B) Analysis of organoid sample (GSE57965). Statistical significance was defined as a false discovery rate (FDR) < 0.25. CRC, colorectal cancer; FDR, false discovery rate.
Cells 10 00488 g004
Figure 5. Pie chart demonstrating the proportion of infiltrating immune cells. (A) Analysis of tissue sample (GSE41258). (B) Analysis of organoid sample (GSE57965). NK cell, natural killer cell.
Figure 5. Pie chart demonstrating the proportion of infiltrating immune cells. (A) Analysis of tissue sample (GSE41258). (B) Analysis of organoid sample (GSE57965). NK cell, natural killer cell.
Cells 10 00488 g005
Figure 6. Comparison of immune cell infiltration. (A) Analysis of T cells. (B) Analysis of monocytes. Tukey-type boxplots demonstrate the median and interquartile level values. Ade, adenoma; CRC, colorectal cancer; Tregs, regulatory T cells; M0, Macrophage M0; M1, Macrophage M1; M2, Macrophage M2; DC, Dendritic Cell.
Figure 6. Comparison of immune cell infiltration. (A) Analysis of T cells. (B) Analysis of monocytes. Tukey-type boxplots demonstrate the median and interquartile level values. Ade, adenoma; CRC, colorectal cancer; Tregs, regulatory T cells; M0, Macrophage M0; M1, Macrophage M1; M2, Macrophage M2; DC, Dendritic Cell.
Cells 10 00488 g006
Figure 7. Comparison of stem cell and progenitor cell infiltration between adenoma and CRC tissues and organoids. (A) Infiltration of stem cells in tissue samples (GSE41258). (B) Infiltration of stem cells in organoid samples (GSE57965). (C) Infiltration of progenitor cells in tissue samples (GSE41258). (D) Infiltration of progenitor cells in organoid samples (GSE57965). Tukey-type boxplots demonstrate the median and interquartile level values. Statistical significance was defined as a false discovery rate (FDR) < 0.25. Ade, adenoma; CRC, colorectal cancer.
Figure 7. Comparison of stem cell and progenitor cell infiltration between adenoma and CRC tissues and organoids. (A) Infiltration of stem cells in tissue samples (GSE41258). (B) Infiltration of stem cells in organoid samples (GSE57965). (C) Infiltration of progenitor cells in tissue samples (GSE41258). (D) Infiltration of progenitor cells in organoid samples (GSE57965). Tukey-type boxplots demonstrate the median and interquartile level values. Statistical significance was defined as a false discovery rate (FDR) < 0.25. Ade, adenoma; CRC, colorectal cancer.
Cells 10 00488 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tokumaru, Y.; Oshi, M.; Patel, A.; Tian, W.; Yan, L.; Matsuhashi, N.; Futamura, M.; Yoshida, K.; Takabe, K. Organoids Are Limited in Modeling the Colon Adenoma–Carcinoma Sequence. Cells 2021, 10, 488. https://doi.org/10.3390/cells10030488

AMA Style

Tokumaru Y, Oshi M, Patel A, Tian W, Yan L, Matsuhashi N, Futamura M, Yoshida K, Takabe K. Organoids Are Limited in Modeling the Colon Adenoma–Carcinoma Sequence. Cells. 2021; 10(3):488. https://doi.org/10.3390/cells10030488

Chicago/Turabian Style

Tokumaru, Yoshihisa, Masanori Oshi, Ankit Patel, Wanqing Tian, Li Yan, Nobuhisa Matsuhashi, Manabu Futamura, Kazuhiro Yoshida, and Kazuaki Takabe. 2021. "Organoids Are Limited in Modeling the Colon Adenoma–Carcinoma Sequence" Cells 10, no. 3: 488. https://doi.org/10.3390/cells10030488

APA Style

Tokumaru, Y., Oshi, M., Patel, A., Tian, W., Yan, L., Matsuhashi, N., Futamura, M., Yoshida, K., & Takabe, K. (2021). Organoids Are Limited in Modeling the Colon Adenoma–Carcinoma Sequence. Cells, 10(3), 488. https://doi.org/10.3390/cells10030488

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