| Literature DB >> 31468733 |
Takashi Aizawa1, Hideaki Karasawa1, Ryo Funayama2, Matsuyuki Shirota3, Takashi Suzuki4, Shimpei Maeda1, Hideyuki Suzuki1, Akihiro Yamamura1, Takeshi Naitoh1, Keiko Nakayama2, Michiaki Unno1.
Abstract
Recent studies have shown that the tumor microenvironment plays a significant role in the progression of solid tumors. As an abundant component of the tumor microenvironment, cancer-associated fibroblasts (CAFs) have been shown to promote tumorigenesis and cancer aggressiveness, but their molecular characteristics remain poorly understood. In the present study, paired CAFs and normal fibroblasts (NFs) were isolated from five colorectal cancer (CRC) tissues from patients who underwent surgical resection. The gene expression profiles of CAFs and NFs identified by RNA sequencing were compared to understand the complex role of CAFs in cancer progression. Gene Set Enrichment Analysis revealed that the gene sets related to the Wnt signaling pathway were highly enriched in CAFs, as well as TGFβ signaling, which is considered to be a regulator of CAFs. Among the components of this pathway, Wnt2 was specifically expressed. The observations led us to speculate that Wnt2 is extremely involved in regulating CRC progression by CAFs. Thus, we performed immunohistochemical analysis on Wnt2 in 171 patients who underwent surgery for colorectal adenocarcinoma. Positive staining for Wnt2 was mainly observed in cancer stroma, although the immunoreactivity was weak in cancer cells. Wnt2 expression in CAFs was significantly correlated with depth of tumor (P < .001), lymph node metastasis (P = .044), TNM stage (P = .010), venous invasion (P < .001), and recurrence (P = .013). Subsequent in vitro analyses were conducted using conditioned medium (CM) from immortalized CAFs transfected with siRNA targeting Wnt2. As a result, cancer cell invasion and migration were significantly decreased in the CM from immortalized CAFs transfected with siRNA targeting Wnt2. Our findings indicated that Wnt2 protein released from CAFs enhances CRC cell invasion and migration. In conclusion, Wnt2 secreted by CAFs plays a key role in cancer progression and is a potential therapeutic target for CRC.Entities:
Keywords: RNA sequencing; Wnt2; cancer-associated fibroblast; colorectal cancer; gene set enrichment analysis
Mesh:
Substances:
Year: 2019 PMID: 31468733 PMCID: PMC6797671 DOI: 10.1002/cam4.2523
Source DB: PubMed Journal: Cancer Med ISSN: 2045-7634 Impact factor: 4.452
Characteristics of patients whose colorectal tissues were examined for isolating primary fibroblasts
| No. | Age | Gender | Location of cancer | Histological differentiation | TNM classification |
|---|---|---|---|---|---|
| 1 | 73 | Male | Rectal cancer | Well differentiated | T3N1bM1a StageIVA |
| 2 | 68 | Male | Rectal cancer | Moderate differentiated | T3N1bM1a StageIVA |
| 3 | 58 | Female | Rectal cancer | Moderate differentiated | T3N2bM0 StageIIIC |
| 4 | 66 | Male | Rectal cancer | Moderate differentiated | T4aN0M0 StageIIB |
| 5 | 61 | Male | Rectal cancer | Moderate differentiated | T4aN2aM1a StageIVA |
Figure 1Fibroblasts isolated from human colorectal cancer tissues. Isolated cells showed a spindle‐like shape, which makes no morphological distinction between CAFs and NFs (A). They were negative for the epithelial marker E‐cadherin (B). The expression of αSMA, vimentin, and FAP in these cells did not show a specific pattern between CAFs and NFs (C‐E)
Figure 2Analysis of gene expression profiling. The heat map demonstrated differentially expressed genes identified with RNA sequencing, including 283 upregulated genes and 301 downregulated genes in CAFs (A). Pathway analysis using Gene Set Enrichment Analysis revealed that the gene sets related to the Wnt signaling pathway were highly enriched in CAFs (B). CAFs had much higher WNT2 expression than NFs among the primary fibroblasts from each of the five cases (C). FPKM: fragments per kilobase of transcript per million mapped reads
KEGG Pathway analysis
| Name | SIZE | ES | NES | NOM | FDR q‐val |
|---|---|---|---|---|---|
| KEGG_PROTEASOME | 41 | 0.56163 | 1.9084287 | 0 | 0.025224 |
| KEGG_ARRHYTHMOGENIC_RIGHT_VENTRICULAR_CARDIOMYOPATHY_ARVC | 70 | 0.488159 | 1.8073922 | 0 | 0.042101 |
| KEGG_PATHOGENIC_ESCHERICHIA_COLI_INFECTION | 53 | 0.484691 | 1.7336141 | 0 | 0.062389 |
| KEGG_REGULATION_OF_ACTIN_CYTOSKELETON | 199 | 0.377994 | 1.669645 | 0 | 0.09594 |
| KEGG_FOCAL_ADHESION | 193 | 0.375368 | 1.6675874 | .002551 | 0.078578 |
| KEGG_HEDGEHOG_SIGNALING_PATHWAY | 51 | 0.471096 | 1.6389468 | .012245 | 0.085571 |
| KEGG_ECM_RECEPTOR_INTERACTION | 80 | 0.428931 | 1.6388332 | .008929 | 0.073641 |
| KEGG_VIRAL_MYOCARDITIS | 64 | 0.431181 | 1.6267112 | .002268 | 0.07074 |
| KEGG_CELL_ADHESION_MOLECULES_CAMS | 114 | 0.385292 | 1.6078509 | 0 | 0.072959 |
| KEGG_SMALL_CELL_LUNG_CANCER | 84 | 0.394938 | 1.5541589 | .011261 | 0.103223 |
| KEGG_BASAL_CELL_CARCINOMA | 50 | 0.436296 | 1.5480132 | .010246 | 0.098083 |
| KEGG_LEUKOCYTE_TRANSENDOTHELIAL_MIGRATION | 106 | 0.37584 | 1.5435936 | 0 | 0.093158 |
| KEGG_THYROID_CANCER | 29 | 0.480652 | 1.5274546 | .023504 | 0.097032 |
| KEGG_P53_SIGNALING_PATHWAY | 66 | 0.407492 | 1.5242089 | .016878 | 0.092584 |
| KEGG_HYPERTROPHIC_CARDIOMYOPATHY_HCM | 75 | 0.392986 | 1.5123453 | .002193 | 0.095073 |
| KEGG_RENAL_CELL_CARCINOMA | 69 | 0.3904 | 1.4776138 | .014925 | 0.118615 |
| KEGG_TIGHT_JUNCTION | 123 | 0.352359 | 1.4742628 | .006772 | 0.11432 |
| KEGG_WNT_SIGNALING_PATHWAY | 140 | 0.333361 | 1.4278333 | .008969 | 0.154247 |
| KEGG_TGF_BETA_SIGNALING_PATHWAY | 83 | 0.353121 | 1.4017543 | .031674 | 0.176525 |
| KEGG_ADHERENS_JUNCTION | 73 | 0.352914 | 1.3431549 | .059211 | 0.248549 |
| KEGG_GLYCOSAMINOGLYCAN_BIOSYNTHESIS_HEPARAN_SULFATE | 24 | 0.440585 | 1.3352597 | .112798 | 0.249587 |
| KEGG_MELANOGENESIS | 93 | 0.332118 | 1.3331608 | .051836 | 0.241666 |
SIZE: Number of genes in the gene set after filtering out those genes not in the expression dataset.
ES: Enrichment score for the gene set; that is, the degree to which this gene set is overrepresented at the top or bottom of the ranked list of genes in the expression dataset.
NES: Normalized enrichment score; that is, the enrichment score for the gene set after it has been normalized across analyzed gene sets.
NOM P‐val: Nominal P value; that is, the statistical significance of the enrichment score. The nominal P value is not adjusted for gene set size or multiple hypothesis testing; therefore, it is of limited use in comparing gene sets.
FDR q‐value: False discovery rate; that is, the estimated probability that the normalized enrichment score represents a false positive finding.
Top 20 genes in KEGG Wnt signaling pathway
| Name | Rank in gene list | Rank metric score | Running ES |
|---|---|---|---|
| WNT2 | 24 | 1.482293 | 0.03164 |
| SFRP2 | 32 | 1.371907 | 0.061677 |
| WNT5A | 138 | 1.013468 | 0.078918 |
| PRICKLE1 | 183 | 0.969715 | 0.098214 |
| VANGL2 | 256 | 0.896152 | 0.114493 |
| LEF1 | 448 | 0.784525 | 0.1224 |
| CCND1 | 452 | 0.781046 | 0.13955 |
| PLCB4 | 717 | 0.684703 | 0.141628 |
| LRP6 | 720 | 0.683883 | 0.156675 |
| FZD7 | 852 | 0.64638 | 0.164498 |
| PRKACG | 926 | 0.625185 | 0.174726 |
| WNT4 | 969 | 0.616328 | 0.186294 |
| WNT11 | 1159 | 0.572646 | 0.189608 |
| PORCN | 1171 | 0.570418 | 0.201696 |
| SFRP4 | 1177 | 0.567706 | 0.214022 |
| PPP2R5B | 1460 | 0.509164 | 0.21132 |
| FZD10 | 1696 | 0.471635 | 0.210117 |
| FZD8 | 1718 | 0.468426 | 0.21945 |
| CCND3 | 1841 | 0.450994 | 0.223392 |
| FZD4 | 1949 | 0.444956 | 0.227942 |
Rank in Gene List is the position of the gene in the ranked list of all genes included in our gene expression profile.
Rank Metric Score refers to the signal to noise ratio for each gene used to position the gene in the ranked list.
Running ES is the enrichment score at this point in the ranked list of genes.
Figure 3Wnt2 expression in colon tissue by immunohistochemistry. The images show representative immunohistochemical staining in normal colon tissue (A), cancer tissue with a score of 0 (B), cancer tissue with a score of 1 (C), and cancer tissue with a score of 2 (D). Immunoreactivity was found in cancer stroma, particularly in the fibromatous area, rather than cancer calls, although staining for Wnt2 were rarely observed in both epithelial and stromal cells in normal tissues
Correlations between Wnt2 expression and clinicopathological factors in patients with colorectal cancer
| Clinicopathological factors | Low expression (N = 85) | High expression (N = 86) |
|
|---|---|---|---|
| Age (y) | 69 (35‐92) | 69 (34‐89) | .461 |
| Gender | .826 | ||
| Male | 49 | 51 | |
| Female | 36 | 35 | |
| Location of tumor | .953 | ||
| Right side | 32 | 32 | |
| Left side | 53 | 54 | |
| Size (mm) | 31 (9‐130) | 39.5 (7‐125) | .110 |
| Histologic differentiation | .199 | ||
| Well or Moderate differentiated | 82 | 79 | |
| Others | 3 | 7 | |
| Depth of tumor |
| ||
| T1 + T2 | 42 | 14 | |
| T3 + T4 | 43 | 72 | |
| Lymph node metastasis |
| ||
| N0 | 61 | 49 | |
| N1 + N2 | 24 | 37 | |
| Distant metastasis | .080 | ||
| M0 | 74 | 66 | |
| M1 | 11 | 20 | |
| TNM Stage |
| ||
| I + II | 57 | 41 | |
| III + IV | 28 | 45 | |
| Lymphatic invasion | .303 | ||
| No | 31 | 25 | |
| Yes | 54 | 61 | |
| Venous invasion |
| ||
| No | 27 | 8 | |
| Yes | 58 | 78 | |
| Recurrence (Stage I‐III) |
| ||
| No | 69 | 52 | |
| Yes | 5 | 14 |
The values were presented as median (range) and analyzed using χ 2 test. All other values represented the number of cases and analyzed using Student's t test. P < .05 was considered significant.
P < .05 was defined as significant, and listed in bold.
Figure 4Kaplan‐Meier analysis of disease‐free survival (A), overall survival (B), and disease‐specific survival (C) of patients with colorectal cancer according to the expression of Wnt2 protein. There was no significant correlation between the high‐expression group and the low‐expression group in these variables
Figure 5Wnt2 protein secreted by CAFs promotes the invasion and migration of colorectal cancer cells in vitro. The qRT‐PCR analysis confirmed the downregulation of the Wnt2 mRNA expression levels in immortalized CAFs transfected with siRNA targeting Wnt2 compared with that in those transfected with non‐targeting siRNA (A). ELISA showed that conditioned medium (CM) from immortalized CAFs transfected with siRNA targeting Wnt2 contained lesser Wnt2 protein than that with non‐targeting siRNA (B). Although there was no difference in cancer cell progression among the three types of CM (C), cancer cell invasion and migration significantly decreased in CM from CAFs transfected with siRNA targeting Wnt2 (D, E). CM‐siWnt2‐1: conditioned medium derived from CAFs transfected with siWnt2‐1, CM‐siWnt2‐2: conditioned medium derived from CAFs transfected with siWnt2‐2, CM‐siNC: conditioned medium derived from CAFs transfected with non‐targeting siRNA. *P < .05 and **P < .01 versus control (left bar)