EGFR/ErbB3 is significantly activated in gefitinib-resistant Kras mutant CRC cell lines, and ErbB3 is upregulated in Kras mutant CRC tissues
Initially, we constructed four groups of CRC-GR cell lines, Kras mutant cells (HCT116, LoVo, and SW620) and Kras wild type cells (HT-29), and the dose of gefitinib resistance reached 50 μm. After preliminary screening of the four groups of cells by phosphorylated RTK antibody array, 16 proteins were observed to exhibit altered phosphorylation levels in gefitinib-resistant cell lines. Statistical analysis revealed that EGFR and ErbB3 levels were changed in all three groups of Kras mutant cells (Fig. 1A and Fig. S1F). ErbB mRNA levels were analyzed in drug-resistant cell lines (Fig. 1B), and the data showed that the mRNA levels of both EGFR and ErbB3 were significantly increased in all four groups, while ErbB2/ErbB4 were not significantly changed. Furthermore, EGFR/ErbB3 protein levels were detected in four drug-resistant cell lines, and the data showed that total protein expression levels of EGFR and ErbB3 in GR Kras mutant cell lines were significantly increased, indicating that both EGFR and ErbB3 were significantly activated. However, the total and phosphorylated protein levels of EGFR/ErbB3 in HT-29 cell lines were not significantly changed (Fig. 1C). We also checked for the possible EGFR T790M mutation during TKI resistance, and no change in mutation level was detected in either the subcutaneous tumor model or in drug-resistant cells (Fig. S1B).
To investigate the significance of EGFR and ErbB3 in Kras-mutant CRC, we analyzed the expression of EGFR and ErbB3 in the TCGA database and the relationship with prognosis. The data showed that ErbB3 was the third most altered of 122 upregulated genes in Kras mutant tissues compared to other types CRC mutations (Fig. 1E). Further analysis showed that in CRC patients with high ErbB3 expression, high expression of EGFR was associated with a poorer prognosis than low expression of EGFR (Fig. 1D). Independent analysis of EGFR and ErbB3 revealed no correlation between EGFR expression and prognosis.
Overall, these data indicate that ErbB3 is highly expressed in tissues of patients with Kras mutations and in Kras mutant GR cell lines.
MiR-323a-3p directly targets EGFR/ErbB3 and inhibits CRC cell proliferation in vitro and tumor growth in vivo
To solve the limitations of Kras mutation in current clinical use of TKIs in CRC, we attempted to find microRNAs that can target both EGFR and ErbB3 and inhibit the expression of these two genes to improve the efficacy of targeted therapy. Combining four miRNA target molecular databases (TargetScan, HMDD, miRgator, and miRDB) for comprehensive analysis, miR-323a-3p was selected as our candidate microRNA (Fig. 2A). Next, a luciferase reporter assay was used to verify that miR-323a-3p directly binds to the 3’UTR of EGFR and ErbB3, and the inhibitory effect of miR-323a-3p on target genes was further verified in HCT116 cell lines (Fig. 2B and Fig. 2C). For Kras wild type CRC cells (HCT-8), the inhibitory effect of miR-323a-3p on target genes was also significant (Fig. S1D).
To determine the biological significance of miR-323a-3p in CRC, we analyzed miR-323a-3p expression and its relationship with prognosis in the TCGA database. The data showed that miR-323a-3p was expressed at low levels in CRC tissues (Fig. 2D), and this low expression was associated with poor prognosis in colorectal cancer (Fig. 1E). Then, we validated the miR-323a-3p expression levels by Quantitative real-time polymerase chain reaction (qRT-PCR) assay, and miR-323a-3p was significantly downregulated in CRC tissues compared to their adjacent tissues in our cohort (n=36, Fig. 2F).
Furthermore, we explored miR-323a-3p function in CRC cells. The results showed that the overexpression of miR-323a-3p in CRC cells significantly inhibited cell proliferation (Fig. 2H). Similarly, Kras wild type cells (HCT-8) also showed some inhibition (Fig. S2B). We also conducted a clone formation assay and observed that CRC cells transfected with miR-323a-3p formed significantly fewer clones than control cells (Fig. 2G). Moreover, miR-323a-3p also inhibited the migration and invasion of CRC cells (Fig. S2D and E). Significant differences in cell proliferation levels were observed after transfection of anti-miR-323a-3p into HCT116 and HCT-8 cells. After anti-miR-323a-3p transfection with HCT-8, the colony number and proliferation level increased, while HCT116 transfection caused no significant change (Fig. S2F).
Subsequently, we examined the effect of miR-323a-3p on the biological characteristics of CRC cells in vivo. We subcutaneously injected miR-323a-3p-overexpressing cell lines or control HCT116 cells into nude mice. As shown in Fig. 2I and Fig. S2G, 6 weeks after the injection of HCT116 cells overexpressing miR-323a-3p, compared to the control group, the tumor volume inhibition rate was 75%, and the tumor weight inhibition rate was 71%. In addition, the expression level of miR-323a-3p was detected in the tumor tissue of mice injected with HCT116 cells and was increased in the tumor tissue of mice injected with HCT116 cells transfected with miR-323a-3p (TUZHU).
In conclusion, these data indicate that low expression of miR-323a-3p is associated with poor prognosis in colorectal cancer. Overexpression of miR-323a-3p inhibits cell growth by targeting EGFR and ErbB3 directly in vivo and in vitro. Therefore, miR-323a-3p represents a tumor suppressor microRNA that inhibits EGFR/ErbB3 and improves the efficacy of targeted therapy for colorectal cancer.
MiR-323a-3p inhibits tumor growth by promoting apoptosis in CRC in vitro and in vivo
To further understand the function of miR-323a-3p, we performed hallmark gene enrichment analysis using RNA-seq data from miR-323a-3p-overexpressing cell lines. The data revealed that apoptosis was the second pathway in these enrichment pathways (Fig. 3a). Then, we detected changes in apoptosis-related proteins (caspase3/7/9) after administration of the miR-323a-3p agomir. The results showed that the levels of apoptotic markers were activated in HCT116 cells and subcutaneous tumors in response to overexpression miR-323a-3p (Fig. 3B and C). Similar results were obtained in HCT-8 cells (Fig. S3C).
To investigate the apoptotic function of miR-323a-3p in CRC, we selected 6 CRC cell lines (HCT116, LoVo, SW480, SW620, HCT-8 and HT-29) for flow cytometry analysis, and the results showed that the miR-323a-3p agomir had apoptotic effects on all 6 types of cells (Fig. 3D and E). Similar results were obtained for the changes in proliferation of the 6 types of cells by the detection of RTCA cell growth (Fig. 3F and Fig. S3a). Furthermore, we detected changes in target gene protein levels between cell lines after administration of agomir (Fig. 3G).
The above results indicated that miR-323a-3p inhibits cell proliferation by directly targeting EGFR and ErbB3 to induce apoptosis.
MiR-323a-3p inhibits EGFR/ErbB3 activation in GR CRC cells
ErbB3 is highly expressed in Kras mutant CRC tissues, enabling them to escape EGFR-TKI inhibition. EGFR/ErbB3 was activated in Kras mutant GR CRC cells and escaped gefitinib inhibition. MiR-323a-3p inhibited tumor growth in acquired GR Kras mutant cell lines and ErbB3 high expression Kras mutant CRC tissues by inhibiting the expression of EGFR and ErbB3.
First, the EGFR/ErbB3 inhibition rate of gefitinib was detected in the GR cell lines. Although p-EGFR in drug-resistant cell lines was reduced in response to treatment with 200 nM gefitinib for 48 h, protein levels of p-ErbB3 did not significantly change (Fig. 4B). However, the phosphorylated protein levels of EGFR/ErbB3 were significantly reduced after adding agomir for 72 h, and based on this, we found that the inhibitory effect of gefitinib and agomir was more significant (Fig. 4C). In addition, similar results were obtained in the Kras wild type cell line (HT-29) (Fig. S3D).
These results indicated that the upregulated phosphorylated EGFR/ErbB3 level of Kras mutant CRC cells after drug resistance could no longer be controlled at the effective level by TKIs, and the downstream PI3K/Akt pathway continued to be activated. However, the inhibitory effect of agomir on EGFR/ErbB3 in GR CRC cells was not affected, and the inhibitory effect was enhanced after simultaneous administration of agomir and gefitinib.
The molecular regulatory network of miR-323a-3p
MiR-323a-3p targets EGFR/ErbB3, promotes apoptosis, inhibits tumor growth in vivo and in vitro, and reverses GR due to its dual-targeting properties. However, its molecular mechanism has not yet been elucidated, so we continued to analyze the RNA-sequencing data of miR-323a-3p-overexpressing HCT116 cells and conducted a pathway analysis by combining the findings with the KEGG database (Fig. 5A). The results demonstrated that the differentially expressed genes were enriched in multiple EGFR/ErbB3-related pathways. Then, we detected the mRNA expression profile by RNA-seq, and the data showed that 88 genes were downregulated in miR-323a-3p HCT-116 cells compared to the control group (Fig. 5B). Then, we detected the expression of these genes in overexpressed HCT116 and HCT-8 cells by qRT-PCR (Fig. 5B), and the results showed that EGFR and ErbB3 were significantly downregulated. Then, we verified the protein level changes of the significantly different genes (PI3K, c-FOS, Akt, GSK-3β, ERK1/2, MMP9, PCNA and p21) through western blot and constructed a protein interaction network of the different genes in the above sequencing results using WEB-String (Fig. 5C and D). These results reveal the changes in the primary downstream phosphorylation levels of miR-323a-3p targeting EGFR/ErbB3 and the resulting molecular regulatory network. In addition, the rescue experiment we conducted with an EGFR agonist (NSC228155) is shown in Figure S4.
MiR-323a-3p and gefitinib synergistically inhibit CRC cell proliferation
Interestingly in the process of solving the drug resistance problem, we found that miR-323a-3p not only relieved CRC GR but also showed a stronger inhibitory effect on p-EGFR/p-ErbB3 protein levels in combination with gefitinib and miR-323a-3p than in the single administration group (Fig. 4C and Fig. 5A). Consequently, when exploring the apoptotic function of miR-323a-3p, we also attempted to combine gefitinib and miR-323a-3p and found that they had a strong effect on promoting apoptosis in CRC cells. Based on Bliss data analysis, the optimal concentration of the combination was finally determined (Fig. 5D). Reviewing the clinical trial data of gefitinib in the early stage of CRC, we found that, unlike the obvious effect of cetuximab on tumor size, previous studies found that the ErbB3/PI3K/Akt pathway also plays a negative regulatory role in trastuzumab-mediated p27Kip1-mediated G1 cell cycle arrest in breast cancer [14]. The lack of ErbB3 expression in poorly differentiated colorectal cancer cells enhances gefitinib sensitivity [15]. The antitumor mechanism of gefitinib is the promotion of cell cycle arrest and the reduction of the expression of tumor metastasis-related proteins [6] rather than the direct effect of miR-323a-3p on CRC cell apoptosis. This may also explain why gefitinib has long been a poor treatment for CRC. In this study, the synergistic effect of gefitinib and miR-323 in the combined treatment group was pro-apoptotic (Fig. 5B and C) based on the cell cycle arrest effect of gefitinib [7] and miR-323 promoting apoptosis through 3'‐ UTR binding to the EGFR/ErbB3 target gene. In addition, two inhibitors targeting individual ErbB family members have been shown to have significantly higher antitumor activity than single drug therapy alone [2].
MiR-323a-3p and gefitinib synergistically inhibit tumor growth, and miR-323a-3p blocks the formation acquired gefitinib resistance in a xenograft model
At present, clinical TKI treatment is ineffective for Kras mutant CRC, while agomir has a significant proapoptotic effect at the cellular level showing a synergy in combination with gefitinib on the proapoptotic effect. Therefore, we next evaluated the effect of miR-323a-3p combined with gefitinib on the growth of subcutaneous tumor HCT116 cells. Subcutaneous tumors grew rapidly from HCT116 cells, and gefitinib monotherapy did not inhibit their growth. Treatment with agomir after peritumor injection resulted in tumor shrinkage. These tumors rapidly regrew after discontinuation of agomir (Fig. 7A, B). After combined administration of miR-323a-3p and gefitinib, the subcutaneous tumors composed of HCT116 cells maintained a period of growth arrest, and the survival time was much longer than that of the monotherapy group.
Additionally, the HCT116 subcutaneous tumor model showed similar results in tumor weight and volume four weeks after administration (Fig. 7D, E), but the tumor weight and volume were slightly decreased after administration of miR-323a-3p or gefitinib, with tumor growth inhibition rates of 27% and 33%, respectively. After combination administration, tumor weight and volume were smaller than those of the single administration group, and tumor growth remained stagnant. The EGFR/ErbB3 protein expression results were also consistent with the cell experiments, indicating the synergistic effect of the combination drug treatment. More meaningfully, we extracted four groups of subcutaneous tumors for the generation of cell culture and the treatment of IC50 detection. The results showed that the single IC50 agomir group exhibited 47% inhibition compared to the control group. Furthermore, treatment for four weeks after lavage treatment, the primary cell IC50 increased by 52% compared to the control group, and the combination group of IC50 values were much lower than the single drug groups, up by 59% compared to the control group.
The above results indicate that the reversal of GR by agomir was not limited to the occurrence of drug resistance but also suggested that miR-323a-3p may prevent the formation of drug resistance by maintaining low levels of EGFR/ErbB3.