Graphical abstract graphic file with name fx1.jpg [61]Open in a new tab Highlights * • CRC shows a fetal/regenerative/revival state following YAP activation * • YAP transcription is promoted by activation of AP-1 and TEAD via mechanotransduction * • CRC cells in a fetal/regenerative/revival state shows chemoresistance * • Chemoresistance can be resolved by targeting key modules such as TEAD and TNF in CRC __________________________________________________________________ Classification Description: Cancer; Stem cell plasticity; Transcriptomics Introduction Cell plasticity can be displayed through a variety of processes such as epithelial-mesenchymal transition (EMT), dedifferentiation, and transdifferentiation.[62]^1^,[63]^2 In colorectal cancer (CRC) biology, cell plasticity has been widely recognized to be involved in the generation of distant metastasis and acquisition of chemoresistance, which leads to a poor prognosis.[64]^2 For instance, in vivo settings, CRC cells in the invasive front undergo EMT.[65]^3^,[66]^4 More recently, cell fate changes toward a fetal/regenerative/revival state are particularly focused on acquiring better understanding of the cell plasticity of CRC cells. These include a cell fate conversion between LGR5 negative and positive CRC cells[67]^5 or a presence of CRC cells in a fetal progenitor-like state.[68]^6 Importantly, CRC cells in a regenerative state have also been reported to show chemoresistance.[69]^7 However, the underlying mechanisms of cell fate conversion toward a fetal/regenerative/revival state in human CRC have not been fully elucidated yet. It has been reported that cell fate change is inhibited in cancer cells that harbor specific mutations,[70]^8 raising questions about the relationship between cancer mutation profiles and cell plasticity. We and others have previously reported that cell fate conversion toward a fetal/regenerative/revival state is induced in the normal intestinal epithelium [71]^9^,[72]^10^,[73]^11 upon inflammation and tissue regeneration. The process is orchestrated, at least partially, by collagen type I deposition in the mesenchyme following YAP/TAZ activation during tissue inflammation.[74]^9^,[75]^12 Accordingly, the organoid culture system of normal colonic epithelium in purified collagen type I[76]^13 faithfully recapitulates the process.[77]^9^,[78]^12 This enables the establishment of cystic organoids in a fetal/regenerative/revival state driven by the cooperative activation of several transcriptional factors such as YAP/TAZ-TEAD, AP-1, and RUNX2.[79]^12 In the current study, we hypothesize that the ECM niche plays a crucial role in understanding cell plasticity in human CRC. To test this hypothesis, we exposed human CRC organoids with naturally selected genetic mutation profiles to two conditions with different ECM components. Primary cancer cells cultured on a gel predominantly composed of collagen type I showed a transcriptomic signature closer to a fetal/regenerative/revival state with activated YAP transcription, at least within tested genetic mutation constellations. Mechanistically, YAP transcription was enhanced through the activation of AP-1 and TEAD-dependent transcription, coupled with the suppression of intestinal lineage-determining transcription factors such as CDX2, HNF4A, GATA6, and FOXA2.[80]^14 Interestingly, specific genomic mutations, such as RAS mutations, seemed to be associated with the enhancement of the YAP activation. These results suggest that CRC cells with genetic mutations maintain their epithelial pattern in response to the surrounding ECM niche. This, in turn, forms the basis for a cancer-specific reaction with a loss of colonic epithelial determinants and gain of a fetal/regenerative/revival signature. The present study also confirmed that an induced fetal/regenerative/revival state is correlated with chemoresistance to typical cytotoxic anticancer agents in CRC cells. However, CRC cells in an induced fetal/regenerative/revival state showed a clear sensitivity to the inhibition of TNF, which was identified as the largest hub gene in the gene network of these cells. Effective suppression of these cells was also achieved through the inhibition of crucial modules in YAP transcription by using a TEAD inhibitor. By using two in vitro models with different ECM components, it would be possible to set up a novel platform to formulate a unique approach to CRC research. Results Identifying conditions that induce a fetal/regenerative/revival state in established colorectal cancer cells from primary human tumors Invasion of CRC is known to be characterized by a transient loss of differentiation.[81]^3^,[82]^15^,[83]^16 A histopathological analysis of CRC surgical specimens validated the heterogeneous distribution of tumor cells in terms of their degree of differentiation ([84]Figure 1A). Differentiated and polarized tumor glands were observed in the tumoral center. In contrast, poorly differentiated tumor cells were observed occasionally in the invasive front ([85]Figure 1B). Subsequently, we analyzed the distribution of major ECM components in human CRC tissues to identify conditions that induce a fetal/regenerative/revival state in human CRC. Among various types of collagen, which is one of the most prevalent ECM scaffolding proteins in the mesenchyme,[86]^17 collagen type I is known to be especially expressed in the mesenchyme with an increase in the invasive front of CRC, where poorly differentiated cancer cells are present.[87]^3 IHC analysis of CRC surgical specimens confirmed that collagen type I expression was higher in the invasive front than in the central differentiated area, particularly with an abundant fibril formation in the mesenchyme ([88]Figure 1B). It was also revealed that hyaluronan was expressed with a marked increase in the mesenchymal compartment of the invasive front compared to the tumoral center ([89]Figure 1B). As for a component of basement membranes, laminin surrounded differentiated cancer cells in the tumoral center. Notably, its expression was lost in the invasive front ([90]Figure 1B), as previously reported.[91]^3 The distribution of collagen type IV was slightly upregulated in the mesenchymal compartment of the invasive front compared to the tumoral center ([92]Figure 1B). Figure 1. [93]Figure 1 [94]Open in a new tab Identifying conditions that induce a fetal/regenerative/revival state in established colorectal cancer cells from primary human tumors (A) A tiled image of a colorectal cancer surgical specimen is shown. Scale bar, 500 μm. Regions indicated with squares are enlarged in (B) for both the tumoral center and invasive front. (B) Differentiated and polarized tumor glands are observed in the tumoral center (top), whereas de-differentiated tumor cells are observed in the invasive front (bottom) after H.E. staining. Immunohistochemical images of collagen type I, hyaluronan, laminin, and collagen type IV (all green) are shown in the tumoral center (top) and invasive front (bottom). Images are counterstained with E-cadherin (red). For collagen type I and type IV, magnified images of the regions indicated with squares are shown in adjacent panels. Scale bars, 100 μm. Col type I; collagen type I, Col type IV; collagen type IV, E-cad; E-cadherin. (C) A diagram of different culture methods is presented. Matrigel cancer organoid (3D MG) is cultured in laminin-enriched Matrigel under L-WRN conditioned medium (upper panel). 2D HC is cultured on the mixed gel composed of collagen type I, collagen type IV, and hyaluronan under a serum-free ENRW medium (lower panel). (D) Both 3D MG (upper panel) and 2D HC (lower panel) show a propagation over time. Scale bars, 200 μm (lower panel) and 100 μm (lower panel). (E) Phalloidin staining (gray) reveals that the apical F-actin belt is well formed in 3D MG (upper panel), whereas a massive formation of actin stress fiber is observed in 2D HC (lower panel). Images are counterstained with DAPI (blue). Scale bars, 50 μm. (F) The mutational fingerprint of the nine CRC surgical specimens used for primary culture for the most common alteration, as determined by next-generation sequencing assay employing the Illumina AmpliSeq Cancer HotSpot Panel v2, is displayed. The color code indicates the type of genetic alteration, distinguishing between missense mutation, nonsense mutation, frameshift mutation, and splice mutation. See also [95]Figure S1, and [96]Videos S1, [97]S2, [98]S3, [99]S4, and [100]S5. We then explored whether combining collagen type I, hyaluronan, and collagen type IV allows CRC cells to expand in vitro on the synthetic gel (HC gel). Primary cancer cells showed a certain propagation as monolayer sheets, referred to as 2D HC, while standard budding CRC organoids were formed in laminin-enriched Matrigel (3D Matrigel; 3D MG) ([101]Figures 1C and 1D). The composition ratio of collagen type I, hyaluronan, and collagen type IV was 70%v/v, 15%v/v, and 15%v/v respectively. Notably, the HC gel did not allow for the 3D expansion of CRC cells ([102]Figure S1A). Ki-67 positive CRC cells were found in 2D HC, indicating that the proliferation of CRC cells was rigorously maintained in vitro ([103]Figure S1B). Time-lapse imaging showed active mitotic cell divisions in the central part of 2D HC ([104]Video S1 and [105]S2). Upon detailed inspection, multipolar cell division reflecting chromosomal instability[106]^18 ([107]Video S3) or a lamellipodia-like structure[108]^19 was well observed in 2D HC ([109]Video S4). In addition, the single-cell amoeboid-like cell displaying a rounded morphology[110]^20 was actively moving on the gel ([111]Video S5). Considering that these motilities are accompanied by the reorganization of the actin cytoskeleton,[112]^21 we investigated the intracellular actin distribution. In 3D MG, the apical F-actin belt was well formed, indicating that cellular polarity was maintained ([113]Figure 1E). In contrast, 2D HC showed a massive formation of actin stress fibers ([114]Figure 1E). Video S1. Time-lapse imaging of 2D HC, related to Figure 1 [115]Download video file^ (25.3MB, mp4) Video S2. Active mitotic cell divisions in 2D HC, related to Figure 1 [116]Download video file^ (10.6MB, mp4) Video S3. A Multipolar cell division in 2D HC, related to Figure 1 [117]Download video file^ (7MB, mp4) Video S4. A lamellipodia-like structure in 2D HC, related to Figure 1 [118]Download video file^ (24.2MB, mp4) Video S5. The single-cell amoeboid-like cell actively moving on the HC gel, related to Figure 1 [119]Download video file^ (24.1MB, mp4) Intriguingly, 2D HC converted efficiently into 3D MG with budding formation when transferred to Matrigel. Conversely, 3D MG could be passaged as 2D HC when seeded on the HC gel ([120]Figure S1C). The morphological reversibility observed between 3D MG and 2D HC suggests that cancer cell phenotype could flexibly shift in an environment-dependent manner, rather than different culture conditions leading to the selection of a specific subset of CRC cells. In the current study, the organoids were established from nine cases of left-sided differentiated CRC, including eight cases of rectal CRC (sample from #1–#8) and one case of descending CRC (sample #9) ([121]Table 1). The mutational fingerprint of these CRC surgical specimens used for primary culture was evaluated using a next-generation sequencing assay, employing the Illumina AmpliSeq Cancer HotSpot Panel v2. Regarding major significantly mutated genes in non-hypermutated CRC,[122]^22 the positivity for APC, TP53, KRAS, PIK3CA, and FBXW7 mutation was observed in 77.8% (7/9 cases), 88.9% (8/9 cases), 33.3% (3/9 cases), 11.1% (1/9 cases), and 22.2% (2/9 cases), respectively. These percentages suggest that the case selection process did not introduce any obvious bias ([123]Figure 1G). Table 1. Clinical information of colorectal cancer (CRC) surgical specimens for the establishment of 3D MG and 2D HC Cases Gender Location Pathological Diagnosis TNM classification Pathological stage #1 Male Rectum tub1>tub2>por pT3N0M0 ⅡA #2 Male Rectum tub1>tub2>por pT2N0M0 Ⅰ #3 Male Rectum tub2>tub1>por pT2N0M0 Ⅰ #4 Male Rectum tub1>tub2 pT3N0M0 ⅡA #5 Male Rectum tub1>tub2 pT2N0M0 Ⅰ #6 Male Rectum tub2>tub1 pT3N0M0 ⅡA #7 Female Rectum tub1>tub2>por pT2N0M0 ⅡA #8 Female Rectum tub1>tub2 pT4aN2a.m.1a IVA (Liver) #9 Male Descending tub2>por>tub1 pT3N0M0 ⅡA [124]Open in a new tab CRC staging in the current study is performed using the TNM classification system based on guidelines of the 8th edition of the Union for International Cancer Control (UICC). tub1, well differentiated tubular adenocarcinoma; tub2, moderately differentiated tubular adenocarcinoma; por, poorly differentiated adenocarcinoma. RNA sequencing revealed the acquisition of fetal/regenerative/revival signature in 2D HC To compare the transcriptional profile of 2D HC and 3D MG established from nine CRC surgical specimens mentioned above, the original 3D MG lines and their derived 2D HCs were subjected to bulk RNA sequencing analysis. Differentially expressed genes (DEG) (log[2]FC ≥ 1, p value < 0.05) were selected following DESeq2 analysis ([125]Table S1). A total of 865 genes were identified as upregulated in 2D HC (2D HC gene set) and 501 genes were upregulated in 3D MG (3D MG gene set). In the GO TERM (molecular function) analysis, signaling receptor activator activity, cytokine activity, extracellular matrix (ECM) structural constituent, and heparin binding were related to the 2D HC gene set ([126]Figure 2A). In the KEGG enrichment analysis, mechanotransduction-related pathways,[127]^23 such as ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathway, and inflammatory-related pathways, such as cytokine-cytokine receptor interaction and TNF signaling pathway, were related to the 2D HC gene set ([128]Figure 2B). The gene set enrichment analysis (GSEA) of the 2D HC gene set revealed significant enrichment in hallmarks such as “TNFA-SIGNALING-VIA-NFKB” (normalized enrichment score (NES) = 2.98, the false discovery rate (FDR) < 0.01), “EPITHELIAL-MESENCHYMAL-TRANSITION” (NES = 2.83, FDR<0.01), “INTERFERON-GAMMA-RESPONSE” (NSE = 2.61, FDR<0.01), “INTERFERON-ALPHA-RESPONSE” (NES = 2.48, FDR<0.01), “INFLAMMATORY-RESPONSE” (NES = 2.40, FDR<0.01) ([129]Figure 2C). Moreover, GSEA with a reported fetal signature[130]^24 showed a highly significant enrichment in 2D HC (NES = 2.09, FDR<0.01) ([131]Figure 2C). Furthermore, GSEA with a reported tumor budding signature, which consists of upregulated genes in tumor budding cells compared to those of the tumoral center[132]^16 showed a highly significant enrichment in 2D HC (NES = 1.89, FDR<0.01) ([133]Figure 2C). In short, 2D HC demonstrated significant enrichment in inflammation-related genes, epithelial mesenchymal transition, and fetal signature. These phenotypes, indicative of cellular plasticity, suggest that the ECM-mediated mechanotransduction could be one of the drivers of cell fate conversion in human CRC. Figure 2. [134]Figure 2 [135]Open in a new tab RNA sequencing analysis revealed the acquisition of a fetal/regenerative/revival signature in 2D HC (A) The top 10 GO-TERM (molecular function) of the upregulated gene set in 2D HC compared to 3D MG (2D HC gene set) is shown. (B) The top 20 KEGG pathways enriched in the 2D HC gene set are shown. (C) The GSEA of hallmark gene sets, fetal gene signature, and tumor budding signature enriched in 2D HC gene set are shown. NES; normalized enrichment score, FDR; the false discovery rate. (D) The heatmap displays Z score-transformed relative expression levels of genes for crypt base columnar (CBC) signature (LGR5, SMOC2, OLFM4, LRIG1 and EPHB3), regenerative stem cell (RSC) signature (TNF, TNFAPI3, CXCL2, ITGA2, RHOD, ICAM1, PLAU and PLAUR), and representative upregulated genes in normal COL sphere (TACSTD2, ANXA1, LY6D and CYR61) in nine cultures of 3D MG and 2D HC. COL; collagen. (E) A volcano plot of differentially expressed genes included in 3D MG and 2D HC gene sets is indicated. Representative genes are indicated in a plot. See also [136]Table S1. According to a previous report,[137]^7 the heterogeneous cell population in CRC tissues was broadly classified into the Lgr5 positive crypt base columnar stem cell signature (CBC signature)[138]^25 and the Lgr5 negative regenerative stem cell signature (RSC signature).[139]^9^,[140]^26 Although there were some variations in each sample, the expression of CBC signature genes such as LGR5, SMOC2, OLFM4, LRIG1, EPHB3, and IRS1 were generally upregulated in 3D MG. In contrast, the expression of RSC signature genes such as TNF, TNFAIP3, CXCL2, ITGA2, RHOD, ICAM1, PLAU, and PLAUR were generally upregulated in 2D HC ([141]Figure 2D). Representative genes enriched in the collagen sphere of normal intestines, such as TACSTD2, ANXA1, LY6D, and CYR61,[142]^9 were also generally upregulated in 2D HC ([143]Figure 2D). The results suggest that the phenotypic differences between 3D MG and 2D HC reflect the two significant characteristics of CRC tissues, and such phenotypical shifts are induced in the ECM niche-dependent manner. Considering that both SMOC2 and OLFM4 expression decreases in CRC cells at the invasive front,[144]^27^,[145]^28 the differences in transcriptional profiles between 3D MG and 2D HC may reflect the heterogeneity caused by the spatial distribution of ECM components. The 2D HC gene set includes ECM-related genes such as COL1A1, COL4A1, COL4A2, COL6A1, and FN1; collagen receptor genes such as ITGA2; mechanosensor-relative genes such as CD31; fetal marker genes such as TACSTD2, ANXA1, LY6D, and CYR61, which are enriched in the collagen sphere[146]^9; cytoskeleton-related genes such as KRT80, AHNAK2, and AKAP12; ECM degradation-related genes such as PLAU, PLAUR, SERPINE1, and TIMP2; mucin secretion-related genes such as MUC2, MUC16, and TFF2; TNF signaling pathway-related genes such as TNF, CCL2, CXCL2, TNFAIP3, ICAM1, IL1B, CSF2, and MMP14; and members of the AP-1 transcriptional factor family such as BATF, ATF3, and MAFK ([147]Figure 2E). The expression of some 2D HC markers in the invasive front of CRC tissues To verify the biological significance of upregulated genes in 2D HC through the histopathological assessment of CRC tissues, we focused on the overlapping gene signature between the tumor budding gene set, as the tumor budding signature is based on the histological selection of the cell population.[148]^16 Among the 29 genes highlighted ([149]Figure S2A), we identified several cytoskeleton-related genes. Notably, KRT80, a key component of cytoplasmic intermediate filaments,[150]^29 AHNAK2, linked to actin-dependent cytoskeletal formation,[151]^30^,[152]^31 and AKAP12, a critical regulator of actin cytoskeleton reorganization[153]^32 were discerned. We further validated the expression of these three molecules at the protein level through an IHC analysis. Consistent with the gene expression data ([154]Figure 3A), upregulated protein expression of these molecules in 2D HC was well verified in biological triplicates (#1–3) ([155]Figures 3B and 3C). Importantly, an IHC analysis of the three corresponding CRC surgical specimens in sample #1–3 also revealed the higher expression of these cytoskeleton-related molecules in the poorly differentiated cancer cells in the invasive front compared to the differentiated cancer cells in the tumoral center, as shown in the sequential sections including H&E staining ([156]Figures 3D–3F). These results imply that the upregulation of these three molecules in both 2D HC and the tumor invasive front might reflect some commonness between them, although the analyses are too limited to reach a definitive conclusion. Moreover, within the 29 overlapping genes, there might be additional candidate genes warranting further exploration. Figure 3. [157]Figure 3 [158]Open in a new tab The expression of some 2D HC markers in the invasive front of CRC tissues (A) A heatmap of Z score-transformed relative expression levels of KRT80, AHNAK2, and AKAP12, normalized by RLE in 9 cultures of 3D MG and 2D HC is shown. (B) Immunohistochemical images of cytoskeleton-related proteins: KRT80, AHNAK2, and AKAP12 (red) of 3D MG (upper panel) and 2D HC (lower panel) are shown. In MG panels, the region of organoids is indicated with a dashed white line in the left panel. In the right panels, merged images with E-cadherin (gray) are shown. Scale bars, 50 μm. (C) The protein expression levels of KRT80, AHNAK2, and AKAP12 are verified to be upregulated in 2D HC in each case (#1–3), as indicated by the bar graph. Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, unpaired t-test. (D) A tiled image of a colorectal cancer surgical specimen is shown. Scale bar, 500 μm. Regions indicated with squares are enlarged in [159]Figure 5E for both the tumoral center and the invasive front. (E) Differentiated and polarized tumor glands are observed in the tumoral center (top), whereas de-differentiated tumor cells are observed in the invasive front (bottom) after H.E. staining. Immunohistochemical images of KRT80, AHNAK2, and AKAP12 (red) of the tumoral center (upper panel) and the invasive front (lower panel) in the sequential sections are shown. In the H.E. staining map image, the region of the tumoral center is highlighted with a white square, while the invasive front is indicated by a yellow square. In the right panels, merged images with E-cadherin (gray) are shown. Scale bars. 100 μm. (F) The signal intensity of KRT80, AHNAK2, and AKAP12 in immunofluorescence images are shown by the bar graph in individual cases #1–3. Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, unpaired t-test. See also [160]Figure S2. We subsequently analyzed the relationship between the gene expression level and the clinicopathological characteristics of CRC. According to the data from the GEPIA website, KRT80 expression was closely associated with the pathological stage in CRC ([161]Figure S2B). AHNAK2 and AKAP12 are not only significantly upregulated in parallel with cancer progression in advanced histopathological stages but also associated with a poor overall survival and disease-free survival rate ([162]Figures S2C and S2D). These results suggest that cell fate change of CRC cells toward a fetal/regenerative/revival state is involved in clinical features, including a prognosis and progression of CRC. AP-1 and TEAD-mediated transcription and the loss of lineage determination synergize to control the fate of cancer cells To elucidate the ECM-dependent cell fate conversion in CRC, we investigated the transcriptional regulatory mechanisms in detail. We performed a KEGG pathway analysis of upregulated genes identified by RNA-seq (log2FC ≥ 1, p value < 0.05) in both 3D MG and 2D HC across all nine cases. When compared with normal colon crypts isolated from noncancerous tissues of the surgical specimens of cases #1–3 ([163]Table S2), the Wnt signaling pathway was enriched in both 3D MG and 2D HC ([164]Figure 4A). It was speculated to be consistent with CRC-specific gene mutations such as in APC and β-catenin. The Hippo signaling pathway was also enriched in both 3D MG and 2D HC, however ECM receptor interaction and focal adhesion were characteristically identified in 2D HC ([165]Figure 4A). These observations suggest the potential existence of additional amplification in YAP/TAZ/TEAD-dependent signal cascade in 2D HC, akin to the activation of the cascade via mechanotransduction in the normal colon epithelium. Among the typical YAP target genes,[166]^33 six genes (CYR61, CTGF, ANKRD1, F3, LATS2, and TGFB2) showed more than a 2-fold increase of the average TPM value in 2D HC compared to 3D MG ([167]Figure S3A). Moreover, in the analysis of the expression levels of individual genes after subclustering based on their RAS mutational status,[168]^34 expression markedly spiked in 2D HC with RAS mutations in cases #1 and #7 ([169]Figure 4B). The expression level of cytoskeleton-related genes such as KRT80, AHNAK2, and AKAP12, which were revealed to be enriched in the invasive front in the current study, showed a similar trend to the representative YAP signature. Notably, a more distinct distribution was observed upon further investigation through subclustering by RAS mutation status ([170]Figure 4C). This observation aligns with published results showing that the RAS oncogenic program modulates cellular mechanosensing by increasing actin stress fibers.[171]^35 Considering that the RAS mutational status is also reported to be an important factor in determining the inflammatory state,[172]^36 we visualized the expression levels of TNF in RNA sequencing tracks after subclustering by RAS mutational status. The results clearly revealed that TNF is upregulated in both RAS mutation-positive and negative cases, with a notable increase observed in cases #1 and #7 within the RAS mutation-positive group ([173]Figure 4D). Figure 4. [174]Figure 4 [175]Open in a new tab ECM niche-dependent AP-1 and TEAD-mediated transcription and the loss of lineage determination synergize to control the fate of cancer cells (A) The top 10 KEGG pathway enriched in upregulated genes in nine cultures (#1–9) of 2D HC (left panel) or 3D MG (right panel) compared to normal colonic epithelial crypts in three samples (#1–3) is shown. (B) Violin plot analyses comparing the distribution of total transcripts per million (TPM) values of indicated genes (CYR61, CTGF, ANKRD1, F3, LATS2, and TGFB2) in nine cultures (#1–9) of 3D MG and 2D HC and three normal colonic epithelial crypts are shown. Samples are subcategorized by RAS mutation status. WT; wild-type, MT; mutation. The color and shape of each point in the cancer cell lines corresponds to each case. (C) Violin plot analyses comparing the distribution of total transcripts per million (TPM) values of indicated genes (KRT80, AHNAK2, and AKAP12) in nine cultures of 3D MG and 2D HC and three normal colonic epithelial crypts are shown. Samples are subcategorized by RAS mutation status. WT; wild-type, MT; mutation. The color and shape of each point in cancer cell lines correspond to each case. (D) A representative RNA sequencing track visualized in Integrative Genomic Viewer (IGV) at TNF loci is shown. Samples are subclustered by RAS mutation status. TNF is generally upregulated in 2D HC compared to 3D MG in both RAS wild-type (WT) and RAS mutant (MT) cases. The data range of the RNA sequencing track was set as 400 and 2000 in RAS WT and RAS MT cases, respectively. (E) A principal component analysis (PCA) plot from the ATAC sequencing analysis (ATAC-seq) of three 3D MG and 2D HC cultures is shown. (F) Representative motifs among the top 25 most highly enriched transcriptional factor motifs in either 3D MG (2D HC^CLOSE) or 2D HC (2D HC^OPEN) are shown. See also [176]Figure S3, [177]Tables S2 and [178]S3. To dissect the involvement of epigenetic factors, we applied ATAC sequencing analysis (ATAC-seq) in 3D MG and 2D HC in cases #1–3. PCA revealed that the chromatin dynamics of 3D MG and 2D HC were divided into different clusters ([179]Figure 4E). Differentially accessible regions were analyzed using edgeR (log2FC ≥ 1, p value < 0.05) ([180]Table S3). 1464 and 540 enriched accessible peaks were identified in 2D HC (2D HC^open sites) and 3D MG (2D HC^close sites). By integrating ATAC-seq and RNA-seq, TNF, PLAU, F3, MUC16, CCL2, and ATF3 were identified in the overlap between 2D HC up and 2D HC^open genes ([181]Figure S3B). In contrast, several genes among the CBC signature, including OLFM4, SMOC2, BLNK, IRS1, and PHGDH were included in the overlap between 2D HC down and 2D HC^close genes ([182]Figure S3B). To identify the candidate TFs governing cellular reprogramming, we identified enriched motifs in 2D HC^open and 2D HC^close sites. The ATAC-based motif analysis revealed the robust and widespread induction of AP-1 family members such as FRA1, BATF, and ATF3, followed by gradual induction of TEAD family members in the top TFs binding 2D HC^open sites ([183]Figures 4F and [184]S3C). These results are consistent with the previous report on the cooperation between AP-1 and the YAP/TAZ/TEAD complex.[185]^37 Moreover, NF-κB-p65, which is known to be activated by TNF and required for the expression of TNF downstream genes,[186]^38 was also found in the top TFs ([187]Figure S3B). On the other hand, lineage-determining TFs of colonic epithelium such as CDX2, HNF4A, GATA6 and FOXA2[188]^14 occupied the top TFs binding 2D HC^close sites ([189]Figures 4F and [190]S3C). Consequently, it is shown that CRC cells with genetic mutations respond to the surrounding microenvironment, including the ECM niche via 1) the activation of AP-1 and TEAD-dependent transcription and 2) the repression of lineage-determining related transcription simultaneously, which would subsequently lead to a distinct YAP state. TNF serves as a major effector in fetal/regenerative/revival cells which exhibit the potential to convert back again into stem-like CRC cells We performed a gene network analysis of the 2D HC gene set. The protein-protein interaction network of the 2D HC gene set was constructed by STRING and visualized in Cytoscape. TNF was revealed to be the largest node (degree = 120). Furthermore, IL1B (degree = 95), CCL2 (degree = 60), ICAM1 (degree = 51), and CSF2 (degree = 46), which are all closely related to TNF signaling pathway according to the KEGG-database, were also found in the top 12 largest nodes ([191]Figure 5A). Other hub nodes ranked in the top 12 are FN1 (degree = 105), CXCL8 (degree = 78), COL1A1 (degree = 56), THBS1 (degree = 56), SERPINE1 (degree = 41), CD31 (degree = 42), and MYD88 (degree = 39) ([192]Figure 5A). Given that TNF was identified as a top-ranked pathway in the enrichment analysis of RNA-seq ([193]Figure 2C) and exhibited increased chromatin accessibility in 2D HC in the ATAC-seq ([194]Figure S3B; [195]Table S3), it was assumed to be one of the critical effectors of a fetal/regenerative/revival state. To investigate the effect of TNF on cancer cell phenotype, we applied TNF to 3D MG at a concentration of 50 ng/mL for 10 days and found an increase in representative YAP target genes (F3, CYR61, and CTGF) and a fetal marker gene (TACSTD2) as predicted. However, stem cell markers such as SMOC2 and LGR5 showed a decrease upon exposure to TNF ([196]Figure 5B). Interestingly, when medium supernatant from 2D HC was applied to 3D MG for five days, an increase of YAP target genes (F3, CYR61, and CTGF) and a fetal marker gene (TACSTD2), as well as a decrease of stem cell markers (OLFM4, SMOC2, and LGR5), were induced, suggesting that TNF is secreted from 2D HC cancer cells. This may act on secretary cancer cells themselves or neighboring cancer cells in an autocrine/paracrine manner ([197]Figure 5B). Importantly, the most prominent change in these genes was observed in 2D HC, suggesting that the importance lies not only in TNF, the end product induced by a series of cell fate conversion, but also in the upstream transcriptional dynamics of AP-1 and TEAD involved in mechanotransduction ([198]Figure 5B). Figure 5. [199]Figure 5 [200]Open in a new tab TNF serves as a major effector in fetal/regenerative/revival cells which exhibit the potential to convert back again into stem-like CRC cells (A) The protein-protein interaction network of the 2D HC gene set was constructed by STRING and visualized in Cytoscape. Each circle represents a hub node, and the diameter of the circle corresponds with the degree of each node. The largest node, TNF, and the other 11 largest nodes are indicated with gene names. (B) The quantitative polymerase chain reaction (qPCR) of indicated genes, including typical YAP target genes (F3, CYR61 and CTGF), a fetal marker gene (TACSTD2), and stem cell marker genes (OLFM4, SMOC2 and LGR5) in human colorectal cancer (CRC) organoids cultured under different conditions is shown. The average Ct value of each gene is normalized by the average Ct value of GAPDH, and the normalized value is plotted on the Y axis. The average expression level of seven cultures of 3D MG, 3D MG-TNF, 3D MG-supernatant and 2D HC is presented for each gene. Data are represented as mean ± SD. ∗p < 0.05, paired t-test. 3D MG-TNF; Matrigel CRC organoids were exposed to TNF (50 ng/mL) for 10 days 3D MG-Supernatant; Matrigel CRC organoids were cultured in the 2D HC medium supernatant for 5 days. (C) The qPCR results for the indicated genes evaluated in (B) are shown. The average expression level of seven cultures of 3D MG, 2D HC and 3D MG from 2D HC is presented for each gene. Data are represented as mean ± SD. ∗p < 0.05, paired t-test. 3D MG from 2D HC; Matrigel CRC organoids which is passaged from 2D HC. According to the previous reports [201]^5^,[202]^39 multiple cell types including LGR5^high stem-like CRC cells can be reestablished from LGR5^low CRC cells. Therefore, the ability of fetal/regenerative/revival cells to produce stem cells is one of the critical matters related to the observed cell fate conversion. We evaluated this issue by transferring 2D HC into Matrigel to assess whether a fetal/regenerative/revival state would be lost or if a stem cell phenotype could arise. In the qPCR analyses, YAP target genes (F3, CYR61, and CTGF) and a fetal marker gene (TACSTD2) were all repressed after replating from 2D HC to 3D MG. In contrast, stem cell markers (OLFM4, SMOC2, and LGR5) all exhibited an increase. These results illustrate that LGR5^high stem-like CRC can arise from fetal/regenerative/revival cells, supporting the previously reported fate-shift between LGR5^high and LGR5^low state ([203]Figure 5C). Chemoresistance in fetal/regenerative/revival cells can be overcome by inhibiting TEAD/TNF To explore the relation of cancer cell plasticity toward a fetal/regenerative/revival state and the acquisition of chemoresistance, we investigated whether 2D HC was resistant to chemotherapeutic agents. Therapeutic profiling was performed on seven cultures of 3D MG and their derived 2D HC, and these primary cancer cells were treated by three standard chemotherapeutic agents commonly used in CRC patients: 5-Fluorouracil (5-FU), Oxaliplatin, and Irinotecan (CPT-11),.[204]^40 First, 3D MG and 2D HC were treated with 5-FU at various concentrations ranging from 1 to 120 μM for 72 h. At the end of each experiment, cellular viability assay was performed. While 3D MG showed a significant reduction in cellular viability at any concentration, cellular viability remained significantly higher in 2D HC ([205]Figure 6A). The IC50 of 5-FU was significantly higher in 2D HC (IC50 = 91.1 μM) than 3D MG (IC50 = 12.0 μM) ([206]Figure 6A). 2D HC also showed a clear resistance to oxaliplatin and CPT-11 ([207]Figures 6B and 6C). IC50 of oxaliplatin in 2D HC and 3D MG was 45.4 μM and 19.0 μM respectively. IC50 of CPT-11 in 2D HC and 3D MG was 59.1 μM and 15.2 μM respectively ([208]Figures 6B and 6C). The pharmacotyping confirmed that CRC cells in a fetal/regenerative/revival state have a chemoresistance to cytotoxic chemotherapeutic agents in vitro. Figure 6. [209]Figure 6 [210]Open in a new tab Chemoresistance in fetal/regenerative/revival cells can be overcome by inhibiting TEAD/TNF (A) Cell proliferation assay after 72 h treatment with 5-Fluorouracil (5-FU) in seven 3D MG and 2D HC cultures are shown (left panel). The median inhibitory concentration (IC50) of 5-FU is shown (right panel). Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, paired t-test. (B) Cell proliferation assay after 72 h treatment with Oxaliplatin in seven 3D MG and 2D HC cultures is shown (left panel). IC50 of Oxaliplatin is shown (right panel). Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, paired t-test. (C) Cell proliferation assay after 72 h treatment with CPT-11 in seven 3D MG and 2D HC cultures are shown (left panel). IC50 of CPT-11 is shown (right panel). Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, paired t-test. (D) Cell proliferation assay after 72 h treatment with the TEAD inhibitor in seven 3D MG and 2D HC cultures is shown (left panel). IC50 of the TEAD inhibitor is shown (middle panel). Cell viability after 72 h treatment with 5-FU at 20 and 40 μM, TEAD inhibitor at 10 μM and a combination of 5-FU at 20 μM plus TEAD inhibitor at 10 μM is shown (right panel). Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, paired t-test. (E) Cell proliferation assay after 72 h treatment with the TNF inhibitor in seven 3D MG and 2D HC cultures is shown (left panel). IC50 of the TNF inhibitor is shown (middle panel). Cell viability after 72 h treatment with 5-FU at 20 and 40 μM, TNF inhibitor at 7.5 μM and a combination of 5-FU at 20 μM plus TNF inhibitor at 7.5 μM is shown (right panel). Data are represented as mean ± SD. Statistical significance: ∗p < 0.05, paired t-test. (F) Representative cell-cycle profiling of control and 5-FU treated 3D MG (upper panel) and 2D HC (lower panel) determined by fluorescence image cytometer NucleoCounter NC-250 system is shown. Cell populations in S and G2/M phases are indicated by a square. The >4N cell populations are indicated by a dashed square. (G) Representative images of cleaved caspase3 staining (red) of 2D HC after 72 h treatment of DMSO (control) (upper panel) or 5-FU at 40 μM (lower panel) are shown. Images are counterstained with DAPI (blue). The adjacent panels present merged images with Phalloidin (gray) for both control and 5-FU. Scale bars, 50 μm. (H) Structural changes in actin filaments induced by TNF inhibitor administration are depicted by phalloidin staining. Representative images of phalloidin staining (gray) of 2D HC after 72 h treatment of DMSO (control) (upper panel) or TNF inhibitor at 10μΜ are shown. In the adjacent panels, merged images with DAPI (blue) are presented for control and TNF inhibitor. Scale bars, 50 μm. See also [211]Figure S4 and [212]Videos S6 and [213]S7. Importantly, however, when 3D MG and 2D HC were exposed to a potent TEAD inhibitor (K-975), a small molecule that can inhibit YAP/TAZ-TEAD protein-protein interactions,[214]^41 2D HC was efficiently eliminated ([215]Figure 6D). IC50 of the TEAD inhibitor (K-975) in 2D HC and 3D MG was 13.0 μM and 72.1 μM respectively ([216]Figure 6D). Moreover, the combination of the TEAD inhibitor at 10 μM plus 5-FU at 20 μM showed more efficient eradication of 2D HC than both 5-FU monotherapy at 20 and 40 μM and the TEAD inhibitor monotherapy at 10 μM ([217]Figure 6D). Additionally, when 3D MG and 2D HC were exposed to a TNF inhibitor (CAS 1049741-03-8), a small molecule that directly binds to TNF and blocks the formation of trimers,[218]^42 the response showed the opposite pattern to that seen in the cell viability assay of standard chemotherapeutic agents, showing the significant sensitivity of 2D HC to a TNF inhibitor ([219]Figure 6E). IC50 of the TNF inhibitor in 2D HC and 3D MG was 4.6 μM and 20.2 μM respectively ([220]Figure 6E). The combination of the TNF inhibitor at 7.5 μM plus 5-FU at 20 μM showed more efficient eradication of 2D HC than both 5-FU monotherapy at 20 and 40 μM and the TNF inhibitor at 7.5 μM ([221]Figure 6E). To understand the mechanism by which 2D HC exhibits chemoresistance, we turn to cell cycle analysis. In the control state prior to 5-FU treatment, the population of the S phase, which was obviously detected in 3D MG, showed a clear decrease in 2D HC, suggesting that 2D HC has more quiescent characteristics ([222]Figures 6F and [223]S4A). Interestingly, the percentage of >4N cells was statistically higher in 2D HC, which reflects the existence of polyploid cancer cells associated with chromosomal instability ([224]Figures 6F and [225]S4A),[226]^43 supporting the presence of multipolar cell division occasionally seen in time-lapse imaging of 2D HC ([227]Video S3). In 3D MG, 5-FU treatment targeting cancer cells in the S phase[228]^44 caused a significant reduction of the cell population in the S phase and a significant accumulation of cells in the sub-G1 phase, reflecting the cell-cycle arrest in 3D MG ([229]Figures 6F and [230]S4B). However, 2D HC did not show significant changes in the cell cycle following 5-FU treatment ([231]Figures 6F and [232]S4B), which aligns with the established concepts of a slow-cycling state in drug-persistent tumor cells.[233]^45 Cleaved caspase-3 immunostaining revealed that 2D HC did not show apparent changes in cleaved caspase-3 expression after 72 h 5-FU treatment at 40μm M ([234]Figure 6G). The time-lapse imaging, which was performed immediately after application of 5-FU at 40 μM, confirmed that the active cell division was maintained, and the monolayer expansion was observed over time, suggesting that 2D HC acquired chemoresistance to 5-FU ([235]Video S6). However, administration of TNF inhibitor at 10μM induced progressive shrinkage of cancer cells, and time-lapse imaging ([236]Video S7) revealed a complete loss of autonomous cell motility. Phalloidin staining revealed that the formation of actin stress fibers was strongly hindered ([237]Figure 6H). Video S6. Time-lapse imaging of 2D HC in the presence of 5-Fluorouracil, related to Figure 6 [238]Download video file^ (37.1MB, mp4) Video S7. Time-lapse imaging of 2D HC in the presence of TNF inhibitor, related to Figure 6 [239]Download video file^ (116.5MB, mp4) To evaluate cell death caused by the cytotoxic regent functionally, we performed the replating assay and verified whether 3D MG and 2D HC could be passaged after the treatment of 5-FU and could regrow as 3D MG ([240]Figure S4C). The replating assay revealed that 3D MG could not be passaged following the treatment of 5-FU at 40 μM, suggesting the effective induction of cell death ([241]Figure S4D). However, 2D HC was successfully passaged as spherical growing organoids even after the treatment of 5-FU at the same concentration, indicating that 5-FU did not induce total cell death in 2D HC ([242]Figure S4D). In contrast, the replating assay revealed that 3D MG showed non-affected regrowth followed by passage after the administration of TNF inhibitor at 10 μM, whereas 2D HC did not show re-expansion at all, demonstrating that cell death was truly induced in 2D HC ([243]Figures S4E and S4F). These results indicate that a TEAD inhibitor or a TNF inhibitor could be effective in treating aggressive, drug-tolerant CRC cells in a fetal/regenerative/revival state within the limited genetic mutation profile analyzed in the current study. Discussion In the current study, we revealed that cellular plasticity toward a fetal/regenerative/revival state was induced following activation of YAP transcription via mechanotransduction by using two in-vitro culture systems with different ECM components. Interestingly, 2D HC showed a distinct signature from 3D MG through multiple differences in culture conditions, including the ECM component, physical property (2D or 3D), or the serum in the culture medium. However, 2D HC characteristically showed enrichment of ECM-receptor interaction and focal adhesion according to transcriptome profiling using RNA-seq. This implies that the ECM component plays a pivotal role in determining the cell fate. On the whole, the fetal signature[244]^24 as well as the EMT signature were both upregulated in 2D HC, implying that the ECM niche-dependent mechanotransduction synergizes with the genetic mutation-based oncogenic program to promote cell fate conversion and form the landscape of heterogeneity in CRC. The differences in transcriptome profile between 3D MG and 2D HC reflected the two major characteristics of CRC: CBC signature and RSC signature,[245]^7 illustrating the scheme that spatial distribution of the ECM components between the tumoral center and invasive front may influence the cancer phenotype. Expression of cytoskeleton-related genes such as KRT80, AKAP12, and AHNAK2 among upregulated genes in 2D HC were elevated in poorly differentiated CRC cells in the invasive front among limited samples. Additionally, these genes were correlated with the clinical staging and prognosis of CRC. These results suggest that differences in the ECM components of the tumor microenvironment might partially induce morphological or phenotypic heterogeneity. Intriguingly, 2D HC showed the upregulation of mucin secretion-related genes, including MUC16, MUC2, and TFF2. As MUC2 expression was downregulated when normal colonic epithelial cells were cultured in collagen type I,[246]^12 the acquisition of mucinous phenotype in CRC cells under collagen-enrich conditions serves as evidence of how the regenerative reaction equipped in the normal colonic epithelium is altered in CRC cells harboring genetic mutations. Given that the expression of MUC16, MUC2, and TFF2 are reported to be associated with advanced tumor stage or poor prognosis,[247]^46^,[248]^47^,[249]^48 2D HC is a valid tool in the study of CRC of an aggressive phenotype. The importance of the YAP/TAZ/TEAD complex as a regulator of CRC phenotype has been recently reported.[250]^49 In the current study, 2D HC showed a distinct YAP state, especially in the RAS mutation-positive cases. Epigenetic characterization of 2D HC by ATAC-seq revealed that the binding regions of colonic epithelium-specific lineage-determining transcriptional factors such as CDX2, HNF4A, GATA6, and FOXA2[251]^14 were all closed, whereas AP-1 and TEAD binding regions were open. This illustrates that CRC cells harboring mutations in the Wnt signaling pathway achieve further upregulation of the YAP/TAZ target genes. Direct verification of the interaction between AP-1 and YAP/TAZ/TEAD complex will be essential to identify specific YAP regulatory elements in human CRC. This is particularly important to formulate therapeutic strategies for CRC. As we presented in this study, 3D MG showed sensitivity to 5-FU, Oxaliplatin, and CPT-11, whereas 2D HC showed a distinct resistance to these cytotoxic agents, illustrating that cancer cells in a fetal/regenerative/revival state display chemoresistance. Importantly, this study identified that a TEAD inhibitor can be effectively used to eradicate these cancer cells, as the phenotype majorly relies on the activated YAP/TAZ-TEAD transcription mechanistically. TNF, which is identified as the largest hub gene and as a major effector in a fetal/regenerative/revival state might also be a promising therapeutic target. In fact, the administration of a TNF inhibitor, a small molecule directly binding to TNF, induced disruption of actin-stress fibers, leading to a progressive shrinkage of cancer cells and complete loss of autonomous cell motility. It is difficult to determine whether the mechanism of action of TNF is autocrine or paracrine. We hypothesize that both mechanisms are likely to be involved. In any case, an inflammatory network associated with TNF activation seems to play a central role inducing a fetal/regenerative/revival state of CRC cells. The entire network could, therefore, potentially be a druggable target in drug-tolerant persistent CRC cells. In a clinical setting, TNF inhibitors have been the mainstay of treatment for inflammatory bowel disease and autoimmune disease.[252]^50 In the oncology field, a TNF inhibitor was also demonstrated to enhance the anti-metastatic effect of cisplatin in a phase I trial in stage IV lung adenocarcinomas.[253]^51 Further mechanistic studies are warranted to elucidate how TNF inhibition suppresses the transcriptional dynamics via cooperative molecular function between AP-1 and YAP/TAZ/TEAD complex. The current study also clearly illustrates the need for further research on mesenchymal cells, including cancer-associated fibroblasts (CAFs), which produce both ECM proteins and cytokines.[254]^52 In CRC, trophoblasts were reported to produce hyaluronan, which induces mesenchymal phenotype with a poor prognosis.[255]^53 Furthermore, some soluble factors released by fibroblasts were also reported to drive the cell fate conversion toward a fetal/regenerative/revival state in murine CRC models.[256]^54 In addition, a previous report showed that a fetal/regenerative/revival state induced by CAFs was associated with chemoresistance.[257]^55 We believe that the scheme presented by our study is useful to understand the functions of CAFs in detail. Limitations of the study This study’s limitation lies in its focus on a limited number of cases, specifically those involving moderately differentiated CRC samples. However, our culture method could provide a platform for analyzing cellular reprogramming in other colorectal malignancies, such as mucinous adenocarcinoma and neuroendocrine carcinoma. Further analysis in a larger number of cases would be helpful to categorize the typical status of the ECM-mediated mechanotransduction and its resultant reprogramming based on the genetic mutation profile of an individual cancer. Although 3D MG is useful for analyzing the effect of several cytokines, such as the WNT ligand, on cellular plasticity by changing the medium composition, the collaborative application of 3D MG and 2D HC could broaden CRC research and lead to the development of effective therapeutic strategies targeting both genetic mutations and cancer cell plasticity. STAR★Methods Key resources table REAGENT or RESOURCE SOURCE IDENTIFIER Antibodies __________________________________________________________________ Mouse anti-E-Cadherin BD Biosciences Cat# 610181, RRID:[258]AB_397580 Rabbit anti-Collagen I Bioss Cat# bs-10423R Biotinylated hyaluronan binding protein Hokudo Cat# BC41 Rabbit anti-Collagen IV Proteintech Cat# 19674-1-AP, RRID:[259]AB_2878595 Rabbit anti-Laminin-1+2 Abcam Cat# ab7463, RRID:[260]AB_305933 Rabbit anti-KRT80 Proteintech Cat# 16835-1-AP, RRID:[261]AB_1851273 Rabbit anti-AHNAK2 Atlas Antibodies Cat# HPA004145, RRID:[262]AB_1859546 Rabbit anti-AKAP12 Proteintech Cat# 25199-1-AP, RRID:[263]AB_2879954 Rabbit anti-Ki67 Abcam Cat# ab15580, RRID:[264]AB_443209 Rabbit anti-Cleaved Caspase-3 (Asp175) Cell Signaling Technology Cat# 9661, RRID:[265]AB_2341188 Alexa Fluor 555 Goat anti Mouse IgG (H + L) Thermo Fisher Scientific Cat# A-21422, RRID:[266]AB_2535844 Alexa Fluor 555 Goat anti Mouse IgG2a (H + L) Thermo Fisher Scientific Cat# A-21137, RRID:[267]AB_2535776 Alexa Fluor 488 Donkey anti Rabbit IgG (H + L) Thermo Fisher Scientific Cat# [268]R37118, RRID:[269]AB_2556546 Alexa Fluor 555 Donkey Anti-Rat IgG (H + L) Abcam Cat# ab150154, RRID:[270]AB_2813834 DyLight 488 Streptavidin Vector Laboratories Cat# SA-5488, RRID:[271]AB_2336405 __________________________________________________________________ Biological samples __________________________________________________________________ Surgical specimens of human colorectal cancer Tokyo Medical and Dental University Hospital Yoshihito Kano __________________________________________________________________ Chemicals, peptides, and recombinant proteins __________________________________________________________________ Growth factor reduced Matrigel Corning Cat# 356231 Cellmatrix Type I-A Nitta Gelatin Cat# 631-00651 Cellmatrix Type IV Nitta Gelatin Cat# 638-05921 HyStem Sigma Cat# HYS020-1KT murine EGF Peprotech Cat# 315-09 murine Noggin R&D systems Cat# 1967-NG murine R-spondin 1 R&D systems Cat# 3474-RS Murine Wnt3a Abcam Cat# ab81484 N-2 supplement Gibco Cat# 17502-048 B-27 supplement Gibco Cat# 17504-044 Nicotinamide Sigma-Aldrich Cat# 72340 N-Acetyl-L cysteine Sigma-Aldrich Cat# A9165 Gastrin I Sigma-Aldrich Cat# G9145 A83-01 Tocris Cat# 2939 Prostaglandin E2 Nacalai Tesque Cat# 29334-21 Bovine Serum Albumin Sigma Cat# A9576 Y-27632 Tocris Cat# 1254 TNF-α Peprotech Cat# 300-01A 5-Fluorouracil Selleck Chemicals Cat# S1209 Oxaliplatin Selleck Chemicals Cat# S1224 Irinotecan (CPT-11) Selleck Chemicals Cat# S1198 TNF inhibitor Calbiochem Cat# 654256 TEAD inhibitor Selleck Chemicals Cat# E1329 __________________________________________________________________ Critical commercial assays __________________________________________________________________ DNeasy Blood & Tissue Kits Qiagen Cat# 69504 RNeasy Micro Kit Qiagen Cat# 74004 SuperScriptⅢ Reverse Transcriptase Invitrogen Cat# 18080044 SYBR green master mix Qiagen Cat# 20415 Cell Titer 96 Aqueous One Solution Cell Proliferation Assay Kit Promega Cat# G3581 VECTASHIELD Antifade Mounting Medium with DAPI Vector laboratories Cat# H-1200-10 Phalloidin FITC Reagent Abcam Cat# ab235137 __________________________________________________________________ Deposited data __________________________________________________________________ Bulk RNA-seq of human CRC patient tissue and organoids and cancer sheets This study SRA: PRJNA1044535 Bulk ATAC-seq of human CRC organoids and cancer sheets This study SRA: PRJNA1044535 Mendeley Data: [272]https://doi.org/10.17632/gbfgn8kfh8.1 __________________________________________________________________ Experimental models: Cell lines __________________________________________________________________ Human CRC Matrigel cancer organoid (3D MG) This study N/A Human CRC cancer sheet (2D HC) This study N/A L-WRN cells for L-WRN conditioned medium ATCC CRL-3276 __________________________________________________________________ Oligonucleotides __________________________________________________________________ Primer: AHNAK2 Fwd: CAGCAGGGGGTAGCTTTCAA This study N/A Primer: AHNAK2 Rev: CTTGCTGTTGCACCAAGTCC This study N/A Primer: KRT80 Fwd: TTCGACCTCGGGCATCTCTA This study N/A Primer: KRT80 Rev: AACTCCATGTCTGTGCGCTT This study N/A Primer: AKAP12 Fwd: AGAAGAGGTTGTGGCCGAAG This study N/A Primer: AKAP12 Rev: CCACCTCCGTTTTCTGCTCT This study N/A Primer: F3 Fwd: AGCTTTTGAGGGGCTGACTT This study N/A Primer: F3 Rev: GGAAGGTGCCCAGAATACCA This study N/A Primer: CYR61 Fwd: CGCCTTGTGAAAGAAACCCG This study N/A Primer: CYR61 Rev: GGTTCGGGGGATTTCTTGGT This study N/A Primer: CTGF Fwd: GGGAAATGCTGCGAGGAGT This study N/A Primer: CTGF Rev: TCCAGTCGGTAAGCCGC This study N/A Primer: TACSTD2 Fwd: CCCCTTTCGGTCCAACAACA This study N/A Primer: TACSTD2 Rev: AACTCCTCTTCTCCTCGGGT This study N/A Primer: OLFM4 Fwd: AGCTGGAGGTGGAGATAAGAAA This study N/A Primer: OLFM4 Rev: CTGAACCACAGACGGTTTGC This study N/A Primer: SMOC2 Fwd: GGCACATCCACAAGGTACGA This study N/A Primer: SMOC2 Rev: GTCCAGAACGCTGGTCAGAA This study N/A Primer: LGR5 Fwd: CCTTCCAACCTCAGCGTCTT This study N/A Primer: LGR5 Rev: AGGGATTGAAGGCTTCGCAA This study N/A __________________________________________________________________ Software and algorithms __________________________________________________________________ RIAS visualization tool Rhelixa [273]https://www.rhelixa.com/rias/ R package – ClusterProfiler Wu T et al.[274]^56 [275]https://bioconductor.org/packages/release/bioc/html/clusterProfile r.html ShinyGO v.0.77 Ge SX et al.[276]^57 [277]http://bioinformatics.sdstate.edu/go/ Gene set enrichment analysis tool software (version; 4.3.0) BROAD INSTITUTE [278]https://www.gsea-msigdb.org/gsea/index.jsp Cytoscape software Version 3.9.1 Doncheva NT et al.[279]^58 [280]http://cytoscape.org RRID:[281]SCR_003032 Fiji Schneider CA et al.[282]^59 [283]http://fiji.sc RRID:[284]SCR_002285 Adobe Photoshop (version 24.0.1 x64) Adobe [285]https://www.adobe.com/jp/products/photoshop.html RRID:[286]SCR_014199 GraphPad Prism software (version; 9.3.1) GraphPad Prism [287]https://www.graphpad.com/RRID:SCR_002798 NDP.view2 Hamamatsu Cat# U12388-01 IGV (version; 2.4.19) Robinson JT et al.[288]^60 [289]https://software.broadinstitute.org/software/igv/ __________________________________________________________________ Other __________________________________________________________________ Coster 24-well Clear TC-treated Multiple Well Plates Corning Cat# 3524 Nunc 48-well Cell-Culture Treated Multidishes Thermo Scientific Nunc Cat# 150687 Falcon 96-well Clear Flat Bottom TC-treated Culture Microplate Corning Cat# 353072 EZVIEW Glass Bottom Culture Plate LB 24well Iwaki Cat# 5826-024 Matsunami 8-well Chamber Plate Matsunami Cat# SCS-N28 Matsunami 2-well Chamber Plate Matsunami Cat# SCS-N22 Tissue-Tek O.C.T. Compound Sakura Finetek USA Cat# 4583 High Profile Microtome Blade Feather Cat# FHP-01 [290]Open in a new tab Resource availability Lead contact Further information and requests regarding resources and reagents should be directed to and will be fulfilled by the lead contact, Shiro Yui (yui.arm@tmd.ac.jp). Materials availability There are restrictions to the availability of unique materials generated in this study, such as human colorectal cancer (CRC) organoids, due to ethical regulations within Tokyo Medical and Dental University (TMDU). Data and code availability * • Bulk RNA-seq data and ATAC-seq data which have been deposited into Sequence Read Archive in National Center for Biotechnology Information and Mendeley Data are publicly available. Accession numbers are listed in the [291]key resources table. * • This paper does not report the original code. * • Any additional information required to reanalyze the data reported in this paper is available from the [292]lead contact upon request. Experimental model and study participant details Ethics approval and consent for the use of human surgical specimens Surgical specimens were obtained at Tokyo Medical and Dental University (TMDU) hospital. The Scientific Ethics Committee of TMDU approved the use of this material for research purposes (M2019-285), and informed consent was obtained adequately from all patients. Study participants In total, 9 patients of Japanese descent (females: n = 2: males: n = 7, between 31 and 83 years old) volunteered to participate in the current study. The clinicopathological features of the patients are shown in [293]Table 1. Method details Primary culture of human colorectal cancer cells Human CRC cells were harvested from surgical specimens. The 5mm–10 mm pieces of CRC tissue were dissected from a cancerous lesion of the CRC surgical specimen and cut into small pieces. The fragments were washed in 30 mL Phosphate Buffer Saline (PBS; Nacalai Tesque, #14249-24) twice and subsequently in 2% Mucofilin (Eisai) once, followed by 4°C incubations in 15 mM EDTA for 20 min. After being washed in 30 mL PBS, the tissues were incubated in 10 mL collagenase solution (Sigma-Aldrich; 6.25 mg/ml in PBS, #C7657) at 37°C for 20 min. After settling down, the supernatant was transferred into another 50 mL conical tube by passing through a 70 μm pore-size mesh filter (SPL Life Science, #93070). The total volume was adjusted to 10 mL with 0.1% v/v BSA/PBS into a 15 mL conical tube, and cells were pelleted at 500g for 3 min. When establishing the Matrigel CRC organoids (3D Matrigel; 3D MG), the pellets were suspended in Matrigel (Corning, #356231) and plated onto a 24-well plate (Flat bottom; Corning, #3524) or 48-well plate (Flat bottom; Thermo Scientific Nunc, #150687). Advanced DMEM/F12 (Gibco, #12634-010) with GlutaMAX (Gibco; 1% v/v, #35050-061) and Penicillin/Streptomycin (Nacalai Tesque; 1% v/v, #26253-84) was used as basal medium. The culture medium contains murine EGF (PeproTech; 50 ng/ml, #315-09), L-WRN conditioned medium (20%v/v), Nicotinamide (Sigma-Aldrich; 10 mM, #72340), N-2 (Gibco; 1%v/v, #17502-048)/B-27 supplement (Gibco; 2%v/v, #17504-044), N-Acetyl-L cysteine (Sigma-Aldrich; 1 mM, #A9165), Gastrin I (Sigma-Aldrich; 10 nM, #G9145) and A83-01 (Tocris; 500 nM, #2939). We also established in vitro culture system of primary CRC cells using Hyaluronan/Collagen gel (HC gel) composed of Cellmatrix Type I-A (collagen type I; 70% v/v, Nitta Gelatin, #631–00651), Cellmatrix Type IV (collagen type IV; 15% v/v, Nitta Gelatin, #638–05921), and HyStem (hyaluronan; 15% v/v, Sigma-Aldrich, #HYS020-1KT), and named it 2D HC. In this case, the pellets were suspended on HC gel and placed on the bottom of a 24-well plate. The culture medium for 2D HC contained the same basal medium above with supplementation of murine EGF (50 ng/ml), murine Noggin (R&D Systems; 100 ng/ml, #1967-NG), mouse R-spondin1 (R&D Systems; 500 ng/ml, #3474-RS), murine Wnt3a (Abcam; 300 ng/ml, #ab81484), Nicotinamide (10 mM), N-2 (1%v/v)/B-27 supplement (2%v/v), Prostaglandin E2 (Nacalai Tesque; 2.5 μM, #29334-21) and Bovine Serum Albumin (BSA; 1%, Sigma-Aldrich, #A9576). The medium was subsequently changed every 2–3 days. Y-27632 (Tocris; 10 μM, #1254) was supplemented for the first 2–3 days after the suspension of the pellets. 2D HC was expanded following the passage of 3D MG on HC gel, which was plated at the bottom of a 24-well plate for RNA and ATAC sequencing analysis, on the bottom of a chamber slide (MATSUNAMI, #SCS-N28, SCS-N22) or a 24-well glass plate (Flat bottom; IWAKI, #5826-024) for immunohistochemical analysis, and on the bottom of a 96-well plate (Flat bottom; Corning, #353072) for cell viability assays using chemotherapeutic agents. For the replating assays, 2D HC was sequentially passaged into 3D MG. For passaging, a monolayer cancer sheet was isolated from HC gel by collagenase solution (6.25 mg/ml in PBS) and was mechanically dissociated into fragments by pipetting 15 times using P1000 pipettes. These fragments were subsequently cultured in Matrigel under an L-WRN-conditioned medium as described above. Y-27632 was supplemented for the first 2–3 days after passage. DNA extraction The pellets of CRC cells were partially cryopreserved in CELLBANKER2 (ZENOGEN PHARMA) for DNA extraction. The total DNA of CRC cells was extracted using the DNeasy Blood & Tissue Kit (Qiagen, #69504) according to the manufacturer’s protocol. Target enrichment and sequencing For library construction, the quality of genomic DNA was assessed using the Genomic DNA ScreenTape on a TapeStation 4200 (Agilent Technologies, Santa Clara, CA, USA). The hotspot regions of approximately 2,800 COSMIC mutations from 50 oncogenes and tumor suppressor genes were enriched by multiplex PCR for 10 ng of genomic DNA obtained from each sample, according to the manufacturer’s instructions. The PCR-amplicon-based sequencing libraries were generated using AmpliSeq Library PLUS for Illumina (Illumina). The qualities of the libraries were assessed using High Sensitivity D1000 ScreenTape on a TapeStation 4200. The multiplex barcoded and equally-pooled libraries were sequenced using MiSeq (Illumina) in 151-base-pair (bp) paired-end reads. Alignment and variant calling were performed using the DNA Amplicon Workflow on an Illumina BaseSpace Sequence Hub ([294]https://jp.illumina.com/products/by-type/informatics-products/bas espace-sequence-hub/apps.html). Functional annotations and the possible effects of variants were added using Illumina Annotation Engine on the BaseSpace. RNA extraction The total RNA of 3DMG, 2DHC, and normal colonic epithelial crypts were extracted using RNeasy Micro Kit (Qiagen, #74004) according to the manufacturer’s protocol. Normal colonic epithelial crypts were isolated in the same way as described above after being harvested from a non-cancerous tissue of CRC surgical specimen. The isolated crypts were cryopreserved for RNA extraction. cDNA synthesis and qRT-PCR The synthesis of cDNA involved using superscript III reverse transcriptase (Invitrogen, #18080044) and random primers with 100-500 ng total RNA. Each reverse-transcription product was subjected to a PCR reaction using SYBR green master mix (QIAGEN, #20415), which was run on StepOnePlus Real-Time PCR system (Applied Biosystems). Each assay was performed in seven cases using specific primer sequences. Primer sequences are listed in the [295]key resources table. The expression level was determined by the ΔCt method, and mRNA levels were normalized to that of GAPDH, shown as the fold-change relative to control. RNA sequencing analysis RNA sequencing analysis of 3D MG and 2D HC established from human CRC specimens and normal colonic epithelial crypts was performed by TaKaRa Bio, Inc. (Kusatsu, Japan). Total RNA was reversed transcribed into cDNA with a Clontech SMART-Seq v4 Low Input RNA Kit according to the manufacturer’s instructions. RNA sequencing libraries were constructed from the amplified cDNA using Illumina Nextera XT DNA Library Prep Kits, validated using an Agilent 4200 TapeStation, and sequenced on the Illumina NovaSeq 6000 platform. Gene quantification results were normalized by TPM (transcripts per million). Differentially expressed gene (DEG) analysis was performed with DESeq2 (Version 1.24.0), following normalization of the raw read counts by relative log normalization (RLE). The heatmap of Z score-transformed relative expression levels of particular genes, and a volcano plot of DEGs were generated using the RIAS visualization tool ([296]https://www.rhelixa.com/rias/) with RLE values. Gene Ontology (GO) enrichment analysis was performed and visualized based on DEGs using the Bioconductor package, clusterProfiler4.0 [297]^56 in R software. A KEGG pathway enrichment analysis was performed and visualized based on DEGs using ShinyGO v.0.77 ([298]http://bioinformatics.sdstate.edu/go/).[299]^57 A Gene set enrichment analysis (GSEA) was performed using GSEA tool software provided by the Broad Institute ([300]https://www.gsea-msigdb.org/gsea/index.jsp). The hallmark gene set collection included in the Molecular Signatures Database,[301]^61 human fetal colon epithelium gene signature,[302]^24 and human tumor budding gene signature[303]^16 were used for the analysis. A protein-protein interaction network analysis of differentially expressed genes was constructed by STRING and visualized in Cytoscape software version 3.9.1..[304]^58 Database analyses We investigated the correlation between mRNA expression of KRT80, AHNAK2, and AKAP12 and clinical stages and subsequently analyzed the prognostic value of these genes in CRC patients using the Gene Expression Profiling Interactive Analysis (GEPIA) browser ([305]http://gepia.cancer-pku.cn/). Time-lapse imaging The monolayer expansion of 2D HC was recorded by time-lapse imaging every 5 min (BZ-X710, KEYENCE). Just before live imaging, 5-Fluorouracil and a TNF inhibitor were added to the culture medium. Throughout the imaging process, the culture condition was maintained under 5% CO2 at 37°C, in a stage-top incubator (Tokai Hit, Shizuoka, Japan). The images were analyzed using the BZ-X Analyzer 1.4.0.1 (KEYENCE). Serial images of 2D HC were converted to a time-lapse video in chronological order. Histology and immunohistochemistry For histology and immunohistochemistry, 3D MG was fixed with 4% Paraformaldehyde (Nacalai Tesque, #09154-85) overnight. Samples were embedded in Tissue-Tek O.C.T. Compound (Sakura) and frozen. 8μm-thick cryosections were prepared by Tissue-Tek Polar Microtome/Cryostat Polar (Sakura). 2D HC was cultured on the chamber plate for immunohistochemistry and was subsequently fixed with 4% Paraformaldehyde overnight. CRC tissues were fixed using 10% formaldehyde, embedded in paraffin, cut into 3μm-thick sections using a microtome, and then fixed on slides. Slides of CRC tissue were deparaffinized and subjected to hematoxylin and eosin (H.E.) staining. A tiled image of H.E. staining was obtained by using NanoZoomer S210 (Hamamatsu Photonics, Hamamatsu, Japan). Immunohistochemical staining of 3D MG, 2D HC and parental CRC tissue was performed with various antibodies, including mouse E-Cadherin (BD Biosciences, #610181), rabbit anti-collagen type I (Bioss, #bs-10423R), biotinylated hyaluronan binding protein (Hokudo, #BC41), rabbit anti-collagen type IV (Proteintech, #19674-1-AP), rabbit anti-laminin-1+2 (Abcam, #ab7463), rabbit anti-KRT80 (Proteintech, #16835-1-AP), rabbit anti-AHNAK2 (Atlas Antibodies, #HPA004145), rabbit anti-AKAP12 (Proteintech, #25199-1-AP), rabbit anti-Ki-67 (Abcam, #ab15580) and Cleaved Caspase-3 (Asp175) antibody (Cell Signaling Technology, #9661). Phalloidin FITC Reagent (Abcam, #ab235137) was used for IHC staining of both MG and HC. Nuclei were counterstained with VECTASHILED mounting medium (Vector, #H-1200). Images were captured with a laser scanning confocal microscope (Leica TCS SP8) and were analyzed in Fiji[306]^59 and Adobe Photoshop CS6. Quantification of immunohistochemical analysis For the in vitro CRC cells (#1–3) a series of images of 3D MG and 2D HC were acquired under the same confocal microscope settings. They were subsequently converted into a grayscale adapting the same pre-set. The signal intensity of KRT80, AHNAK2, and AKAP12 was determined by measuring around 200 independent points of one organoid in ImageJ using the multi-point tool and taking their average. As for 3D MG, the signal intensity was determined for 3–5 different organoids in the individual cases. As for 2D HC, the signal intensity was determined in 3–5 different colonies of the cancer sheet of each case (#1–3). In this assay, 2D HC samples were processed as 8μm-thick cryosections corresponding to target samples of 3D MG. In addition, background intensity was determined in each image, and this value was subtracted from the actual measured value to determine the correction value. These corrected values were subjected to the statistical analysis described below. For the parental CRC surgical specimens from #1–3, a series of images of the tumoral center and the invasive front were captured under the same setting of a confocal microscope and were converted into a grayscale adapting the same pre-set. The signal intensity of KRT80, AHNAK2, and AKAP12 was determined by measuring around 200 independent points of the tissues for each image in ImageJ using the multi-point tool and taking their average. The signal intensity was determined in 3–5 different lesions of the tumoral center and the invasive front, respectively. Background intensity was determined in each image, and this value was subtracted from the actual measured value to determine the correction value. These corrected values were subjected to the statistical analysis described below. ATAC sequencing analysis A quantity of 5,000 to 10,000 cryopreserved cells were thawed to room temperature, pelleted, and washed with cold PBS. The cell pellets were resuspended and tagmented using the enzyme and buffer provided in the ATAC-Seq Kit (Active Motif, #53150). Tagmented DNA was then purified using the DNA Purification Buffers, amplified with ten cycles of PCR, and subsequently purified by SPRI bead solution. The resulting material was quantified by the KAPA Library Quantification Kit for Illumina platforms (KAPA Biosystems) and sequenced with PE150 sequencing on the Illumina Novaseq 6000. The quality of the raw paired-end sequence reads was assessed with FastQC (0.11.7). Illumina Nextera adapter sequences were trimmed by Skewer (0.2.2). The trimmed reads were subsequently aligned to the hg38 version of the human genome using Bowtie2 (2.3.4.2), and non-uniquely mapping fragments were removed by Samtools (1.9). Subsequently, the Bowtie2-resultant. SAM files were converted into binary format (BAM) files with Samtools (1.9). The reads alignment was adjusted by using the ATAC-seq QC R package (1.8.5) and duplicate reads were removed using the MarkDuplicates tool by Picard (2.18.11). Peak calling was performed with MACS2 (2.1.2) with the following setting as default: --shift −75 --extsize 150 --nomodel -q 0.01. The peak regions from each sample were merged with Bedtools (2.27.1). Read summarization was performed using featureCounts (1.6.3) and normalized read counts were obtained by trimmed mean of M values (TMM). Principal component analysis (PCA) of the normalized counts was conducted, and each sample was projected onto the 2D plane of the first and second PCA axes using stats (3.6.1) and gplots (3.0.1.1) R packages. Identification of differentially accessible regions (DARs) was performed by edgeR. For track visualization, normalized bigwig files were uploaded to Integrative Genomics Viewer (IGV).[307]^60 Motif enrichment was identified using the hypergeometric test (HOMER). Cell cycle profiling Cell-cycle profiling of 3D MG and 2D HC was determined by fluorescence image cytometer NucleoCounter NC-250TM system (M&S TechnoSystems, Inc.) according to the manufacturer’s protocol. Chemotherapeutic profiling For the drug sensitivity test, 3D MG was seeded on a 48-well plate, and 2D HC was seeded on a 96-well plate. After confirming that 3D MG formed a budding structure and 2D HC widely spread throughout the well after a few days, chemotherapeutic agents such as 5-Fluorouracil (5-FU) (Selleck Chemicals, #S1209), Oxaliplatin (Selleck Chemicals, #S1224), Irinotecan (CPT-11; Sellek Chemicals, #S1198), a TNF inhibitor (Calbiochem, #654256), and a TEAD inhibitor (Selleck Chemicals, #E1329) were added to the culture medium at various concentrations. The concentrations of 5-FU, Oxaliplatin, and CPT-11 were 1, 10, 20, 40, 80, and 120 μM. The concentrations of the TNF inhibitor were 1, 5, 7.5, 10, 20, and 40 μM. The concentrations of the TEAD inhibitor were 0.1, 0.5, 1, 10, 20, 50, 100, 200, and 400 μM. 3D MG and 2D HC were exposed to the chemotherapeutic agents for 72 h, while DMSO-only treated cells were included as controls. At the end of treatment in each experiment, cellular viability was assessed using Cell Titer 96 Aqueous One Solution Cell Proliferation Assay Kit (Promega, #G3581). The median inhibitory concentration (IC50), which was defined as the concentration producing a 50% decrease in the viability of cells compared with untreated controls, was calculated using GraphPad Prism computer software (version; 9.3.1). Quantification and statistical analysis All error bars indicate the SD. Statistical assessment of overlap between different gene sets was done using Fishers’ exact t-test, assuming that the number of total genes is 25,000 in humans in [308]Figure S2A. Student t-tests were performed for determining statistical significance in [309]Figures 3C, 3F, [310]5B, 5C, [311]6A–6E, S4A, and S4B, using GraphPad Prism computer software (version; 9.3.1). All p values of <0.05 were considered significantly different. Acknowledgments