Cetuximab-induced N6-methyladenosine modification of PUM1 mRNA enhanced stability in CRC cells
PUM1 mRNA content was markedly elevated in cetuximab-resistant CRC cells, relative to controls (Fig. 1A, DMSO vs. Cetuximab resistance, p < 0.05). Results of MeRIP-qPCR experiments showed thatm6A content of PUM1 was also enhanced in the cetuximab-resistant CRC cells (Fig. 1B, DMSO vs. Cetuximab resistance, p < 0.05). Knockdown experiments using siRNA showed that inhibition of METTL3 methyltransferase, but not KIAA1429 or METTL14 downregulation, was mediated by PUM1. Knockdown of the m6A demethylases, ALKBH5 and FTO, did not alter PUM1 expression in CRC cells (Figs. 1C and S1A, B). METTL3 knockdown was shown to greatly reduce the m6A modification of PUM1 mRNA in CRC cells (Figs. 1D and E, si-NC vs. si-METTL3, p < 0.05) and the METTL3-mediated m6A modification enhanced PUM1 levels. The YT521-B homology domain family (YTH) proteins have been shown to regulate mRNA stability via m6A [9]. YTHDF1 was physically associated with PUM1 mRNA in HCC cells (Figs. 1F, G, anti-IgG vs. anti-YTHDF1, p < 0.05) and YTHDF1 knockdown reduced PUM1 mRNA stability (Figs.S1C, 1H, I, si-NC vs. si-YTHDF1, p < 0.05). In addition, PUM1 transcript methylation was recognized by the reader “m6A” and YTHDF1 which maintained mRNA stability and enhanced expression through an m6A-YTHDF1-dependent mechanism.
PUM1 Promoted Monocyte-to-Macrophage Differentiation
Monocyte differentiation into macrophages or dendritic cells (DCs) is known to influence CRC resistance[18]. Low levels of the macrophage markers, CD14, CD71 and CD68, and DC markers, CD1a and CD83, in si-PUM1 CRC cells confirmed that PUM1 stimulated the differentiation of THP-1 cells into macrophages and not DCs (Figs. 2A, B, S2A: si-NC vs. si-PUM1, p < 0.05). Macrophages may be classified into M1 or M2 subtypes, depending on activation state and tumor-associated macrophages (TAMs) are known to be M2-like [19]. THP-1 cells exposed to si-PUM1 CRC cells had patterns of expression of the M1 marker genes, NOS2 and human leukocyte antigen DR alpha [HLA-DRA], which indicated an improved M1 phenotype (Figs. 2C and D, si-NC vs. si-PUM1, p < 0.05).
LV-PUM1 CRC cells showed upregulation of the macrophage markers, CD14, CD71 and CD68, and DC markers, CD1a and CD83 (Figs. 3A, B, S2B: LV-NC vs. LV-PUM1, p < 0.05). LV-PUM1 CRC cell lines prevented differentiation toward an M2 phenotype. TAMs showed enhanced M2 macrophage-specific marker genes, CD163 and MRC1, and genes associated with TAM polarization, CCL5 and interleukin-10 [IL-10] (Figs. 3C and D, LV-NC vs. LV-PUM1, p < 0.05). By contrast, when THP-1 cells were exposed to si-PUM1 CRC cells, NOS2 and HLA-DRA decreased (Figs. 3C and D, LV-NC vs. LV-PUM1, p < 0.05).
PUM1 activated the WNT pathway via modulation of DDX5 expression
DDX5 has been shown to promote CRC progression by stimulating the Wnt pathway [20]. Results of the TOP/FOP-flash luciferase assay showed enhanced Wnt signaling activation in PUM1 overexpressing cells whereas DDX5 knockdown suppressed Wnt signaling (Figs. 4A, S.2C, LV-PUM1 + si-NC vs.LV-PUM1 + si-DDX5, p < 0.05). PUM1 knockdown was demonstrated to suppress the Wnt pathway and DDX5 overexpression to prevent the suppression (Figs. 4B, S.2D, si-PUM1 + LV-Control vs. si-NC + LV-DDX5, p < 0.05). qRT-PCR results showed that elevation of PUM1 caused increasedAxin2, c-Myc and cyclin D1 content and DDX5 knockdown prevented these increases (Figs. 4C, E, G, LV-PUM1 + si-NC vs.LV-PUM1 + si-DDX5, p < 0.05). By contrast, reduction of PUM1 reduced Axin2, c-Myc and cyclin D1 content and DDX5 overexpression prevented this decrease (Figs. 4D, F, H, si-PUM1 + LV-Control vs. si-NC + LV-DDX5, p < 0.05). Hence, PUM1 modulated DDX5 expression to activate the WNT pathway.
Paclitaxel inhibited N6-demethyladenosine modification of PUM1 and suppressed the WNT pathway by diminishing PUM1 levels
Paclitaxel in combination with microRNA-34A has been shown to inhibit CRC progression [21]. Dot blot assays showed that paclitaxel treatment reduced the content of m6A in CRC cells (Fig. 5A) and PUM1 m6A content was shown to be reduced by MeRIP-PCR (Fig. 5B, DMSO vs. paclitaxel, p < 0.05). Paclitaxel also reduced levels of PUM1 and DDX5 mRNA transcripts in CRC cells (Figs. 5C and D, DMSO vs. paclitaxel, p < 0.05).
Paclitaxelreducedm6A modification of PUM1 in cetuximab-resistant CRC cells
A PDX model was developed from cetuximab-resistant CRC patients with varying levels of m6A modification of PUM1. Two CRC patients had a high level ofm6A modification and 2 had a low level (Fig. 6A, IgG vs. anti-m6A, p < 0.05). Paclitaxel treatment inhibited the m6A modification of PUM1in the PDX models and reduced PUM1 mRNA levels were seen (#1 and #2: DMSO vs. paclitaxel, p > 0.05; #3 and #4: DMSO vs. paclitaxel, p < 0.05). It is noteworthy that the PDXs with reduced m6A modification of PUM1 had marked growth inhibition on treatment with paclitaxel relative to controls (Fig. 6D, #1 and #2: DMSO vs. paclitaxel, p > 0.05; #3 and #4: DMSO vs. paclitaxel, p < 0.05). In summary, paclitaxel reduced the m6A modification of PUM1 in cetuximab-resistant CRC cells.
Increased PUM1 in human CRC was linked to elevated TAM levels and poor outcome
Levels of PUM1, CD86 and CD163 were measured in 30 pairs of CRC and adjacent non-cancer tissues (ANT) by RT-qPCR. PUM1 and CD163 mRNA levels were found to be elevated in CRC tissues relative to ANTs (Fig. 7A, Tumor vs. ANT, p < 0.05). Correlation analysis demonstrated that PUM1 transcript levels were inversely related to CD86 mRNA (R2 = 0.656, p < 0.001) and directly proportional to CD163 mRNA (R2 = 0.619, p < 0.001, Fig. 7B). PUM1 protein and macrophage invasion status were assessed in CRC tissues from 300 postoperative patients using IHC. PUM1 protein was elevated in CRC versus ANTs and normal tissues (Fig. 7C, Tumor Tissues vs. Adjacent normal, p < 0.05; Tumor Tissues vs. normal, p < 0.05). Numbers of CD206 + and CD163+ TAMs were also increased in tumor tissues with positive PUM1 expression (N = 150) relative to tumor tissues with negative PUM1 expression (N = 150) (Fig. 7D, p < 0.05). PUM1 levels were analyzed with respect to clinicopathological profile of 300 CRC patients (Table 1) and positive PUM1 expression correlated with tumor size > 5 cm (p = 0.005), presence of lymph node metastasis (p = 0.001) and advanced TNM stage (p < 0.001). K-M analysis showed that reduced PUM1 expression was associated with increased survival of CRC patients with TNM stage I + II disease (Figs. 7E, 7F, p < 0.05). Moreover, reduced PUM1 expression was associated with increased survival of CRC patients receiving cetuximab treatment (Figs. 7G and H, p < 0.05).