Exosomes derived from M2 type tumor-associated macrophages promote osimertinib resistance in non-small cell lung cancer through MSTRG.292666.16-miR-6836-5p-MAPK8IP3 axis
Cancer Cell International volume 22, Article number: 83 (2022)
Osimertinib resistance limits the treatment of epidermal growth factor receptor-(EGFR)-mutated non-small-cell lung carcinoma (NSCLC). The mechanisms of osimertinib resistance need to be elucidated to determine alternative treatment strategies. This study explores the role of M2 type tumor-associated macrophage (TAM)-derived exosomal MSTRG.292666.16 in osimertinib resistance, and its related competing endogenous RNA (ceRNA) mechanism.
M2 type TAMs were induced with 200 ng/mL phorbol 12-myristate 13-acetate, 20 ng/mL IL-4 and IL-13, and M2 type macrophage markers were measured by RT-qPCR. Next, the exosomes were isolated and characterized. Tumor formation in nude mice was conducted using H1975 cells under different treatment conditions. Small RNA sequencing was performed on exosomes derived from sensitive and resistant plasma, and ceRNA networks were constructed. Fluorescence in situ hybridization was used to observe the localization of MSTRG.292666.16, and a ceRNA network (MSTRG.292666.16-miR-6836-5p-MAPK8IP3) was selected for further validation.
M2 type TAMs, and M2 type TAM-derived exosomes were successfully induced and isolated. Nude mice results showed that M2 type TAM-derived exosomes and MSTRG.292666.16 overexpression significantly increased tumor volume after administration of osimertinib for 4 weeks. M2 type TAMs were found in the resistant plasma, and MSTRG.292666.16 localized in the cytoplasm of H1975 cells. In addition, the genes in the ceRNA networks were significantly enriched in eight GO terms and seven KEGG pathways, including the MAPK signaling pathway. Subsequently, the levels of MSTRG.292666.16 and MAPK8IP3 significantly increased in both resistant plasma-derived exosomes and M2 type TAM-derived exosomes, while miR-6836-5p levels were significantly reduced. Finally, MSTRG.292666.16, miR-6836-5p, and MAPK8IP3 were part of the same network.
M2 type TAM-derived exosomes promoted osimertinib resistance in NSCLC by regulating the MSTRG.292666.16/miR-6386-5p/MAPK8IP3 axis.
Lung cancer is the most common cancer worldwide and has the highest mortality rate . Non-small cell lung cancer (NSCLC) accounts for 85% of all cases of lung cancer, with a 5-year survival rate of only 21% . Currently, surgery, chemotherapy, and radiation are the traditional treatments for lung cancer; however or patients with advanced disease, surgery is not sufficient to restrain cancer progression . With the continuous update of gene detection technology, epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs) have become an essential method for precision treatment of NSCLC . Osimertinib, a third-generation EGFR-TKI, selectively and irreversibly targets EGFR mutations (e.g., T790M [c.2369 C > T; p.Thr790Met], L858R [c.2572-2573 inseAG, c.2573T > G, c.2573delinsGA, and c.2574delinsGT; p.Leu858Arg], and C797S [c.2389T > A and c.2390G > C; p.Cys797Ser]), while sparing the wild-type EGFR tyrosine kinase . Recent studies have reported that the overall survival of patients with EGFR-mutated advanced NSCLC was up to 38.6 months after osimertinib treatment, and postoperative adjuvant therapy with osimertinib significantly increased disease-free survival [6, 7]. Therefore, osimertinib is a powerful tool in the first-line treatment of EGFR-mutated advanced NSCLC and in adjuvant therapy for NSCLC after surgery. However, osimertinib resistance is a major obstacle in the treatment of NSCLC, and the resistance mechanisms have not been fully elucidated.
Exosomes are nanoscale lipid inclusions that can carry RNA and proteins and are an important system for intercellular signal transduction and substance delivery . Previous studies have indicated that exosomes participate in various process of tumor drug resistance through the transmission of long non-coding RNA (lncRNA), microRNA (miRNA), and mRNA [9, 10]. Macrophages are one of the major populations of tumor-infiltrating immune cells, and tumor-associated macrophages (TAMs) are the primary actors and coordinators of tumor progression in the microenvironment . Macrophages are usually polarized to either M1 or M2 type; while M1 type macrophages have pro-inflammatory properties, M2 type macrophages contribute to anti-inflammatory response and immune homeostasis . M2 type TAMs secrete a variety of cytokines, chemokines, and exosomes, thus promoting invasion and chemoresistance of ovarian cancer cells , as well as transfer miR-21 to gastric cancer cells thus conferring cisplatin resistance . However, the effects of M2 type TAM-derived exosomes on osimertinib resistance in NSCLC remain unknown.
Competing endogenous RNA (ceRNA) can regulate the expression of downstream mRNAs by competing with shared miRNAs, and an imbalance of ceRNA networks influences cell processes and functions, thus leading to the occurrence and development of diseases [15, 16]. In ceRNA networks, lncRNAs, serving as ceRNAs, can participate in the initiation and progression of liver cancer by sponging miRNA response elements. A study by Xu et al. found that lncRNA MIR194-2HG was highly expressed in liver cancer samples and activated the Wnt/β-catenin signaling pathway through the miR-1207-5p/TCF19 axis, thereby playing important roles in the diagnosis and treatment of liver cancer . Our previous study showed that the level of lncRNA MSTRG.292666.16 was significantly higher in the osimertinib-resistant plasma-derived exosomes than in sensitive plasma-derived exosomes, and found that exosomal lncRNA MSTRG.292666.16 could promote acquired resistance to osimertinib and play a key role in drug resistance transmission . However, the ceRNA mechanisms of exosomal lncRNA MSTRG.292666.16 in osimertinib resistance in NSCLC are still unclear.
Therefore, in this study, we further confirmed the osimertinib resistance of M2 type TAM-derived exosomal lncRNA MSTRG.292666.16 to NSCLC in vivo. Moreover, small RNA sequencing of exosomes was performed and combined with the previous lncRNA sequencing data. Furthermore, the ceRNA mechanism of lncRNA MSTRG.292666.16 in NSCLC osimertinib resistance was investigated, and a ceRNA network of lncRNA MSTRG.292666.16-hsa-miR-6836-5p-MAPK8IP3 was selected for further validation. These results would provide new insights into the development of novel therapeutic targets and pathways to overcome osimertinib-acquired resistance in NSCLC.
Materials and methods
The patient samples were collected based on a previous study  with some minor modifications. Briefly, from October 2015 to September 2018, seven NSCLC patients were recruited from the Department of Oncology, Shanghai Pulmonary Hospital. The clinical characteristics of these patients are shown in Table 1. These patients were administrated first-generation EGFR-TKIs. However, after the development of acquired resistance, osimertinib was used to treat these patients for approximately nine months. Among them, four patients developed osimertinib acquired resistance, and three patients remained sensitive to osimertinib. After that, the patients were divided into sensitive (n = 3) and resistant (n = 4) groups, and the plasma of each patient was obtained for exosome isolation. The study was approved by the Ethics Committee of Shanghai Pulmonary Hospital, and informed consent was obtained from all participants.
Cell culture and induction of M2 type TAM
Human lung cancer H1975 cells and THP-1 cells were obtained from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). The H1975 cells were cultured in Dulbecco’s modified of Eagle’s medium (DMEM, Gibco, Grand Island, NY) supplemented with 10% fetal bovine serum (FBS, Gibco) and 1% penicillin/streptomycin (Gibco), and incubated in an incubator with 5% CO2 at 37 °C. THP-1 cells were maintained in RPMI1640 medium (Gibco) with 10% FBS, and then used to induce the formation of M2 type TAM as previously described . Briefly, THP-1 cells were divided into two groups: control and M2-TAM groups. The cells in the control group were cultured in RPMI1640 medium with 10% FBS. However, the cells in the M2-TAM group were incubated in medium containing 200 ng/mL phorbol 12-myristate 13-acetate (PMA, Sigma) for 24 h. Then, 20 ng/mL interleukin-4 (IL-4, PeproTech) and 20 ng/mL IL-13 (PeproTech) were added to the medium and cultured for another 72 h. Thereafter, the cells were harvested to determine the levels of M2 type macrophage markers (CD206, CD163, TGF-β, IL-10, and Arg-1) using reverse-transcription quantitative PCR (RT-qPCR). The successfully constructed cells were then used for exosome isolation.
Isolation and characterization of exosomes from the plasma and M2 type TAMs
The exosomes were isolated from the plasma and M2 type TAMs using differential centrifugation under temperature conditions at 4 °C, as previously reported . The plasma was sequentially centrifuged at 500×g for 10 min, 2000×g for 30 min, and 10,000×g for 30 min. The supernatant was filtered using 0.22 μm filters (Millipore, USA), followed by centrifugation at 120,000×g for 70 min. The sediment, including the exosomes, was resuspended in 200 μL of PBS. For the isolation from the M2 type TAMs, the cells were centrifuged at 300×g for 10 min, 10,000×g for 30 min, and 100,000×g for 70 min. The sediment was resuspended in 1 mL PBS, and then centrifuged at 10,000×g for 60 min. The exosomes were resuspended in 200 μL PBS and stored at −80 °C until use.
The concentrations of the isolated exosomes were measured using a BCA protein assay kit (BOSTER Biological Technology Co. Ltd, CA, USA). Then, a NanoSight NS300 particle size analyzer (NTA, Malvern Panalytical, Malvern, UK) was used to analyze the particle size of the exosomes based on the methods described by Soares et al. . According to the protocols of Zhu et al. , the morphology of the exosomes was visualized using transmission electron microscopy (TEM, JEOL LTD, Peabody, MA, USA). Finally, western blotting was used to determine the expression of exosome-specific markers (CD9, HSP70, and TSG101) with their corresponding antibodies (1:1000, Proteintech Group, Inc., Rosemont, IL, USA) as previously described .
Small RNA sequencing of plasma-derived exosomes
The exosomes isolated from the plasma prior to osimertinib treatment (sensitive group) and their corresponding drug-resistant plasma (resistant group) were delivered to Yanzai Biotechnology (Shanghai) Co. Ltd (Shanghai, China) for small RNA sequencing . The clean reads were aligned to the databases of human miRNA, genome, and Rfam, and Burrows–Wheeler Aligner (BWA) mapping software was used to analyze the gene expression profiles. Next, the DESeq algorithm was applied to screen the differentially expressed miRNAs (DEmiRNAs) with the thresholds of |log2Fold change (FC)|> 1 and false discovery rate (FDR) < 0.05. Then, the Miranda and RNAhybrid algorithms (http://www.microrna.org) were used to predict the target genes of the identified DEmiRNAs. Afterwards, these screened DEmiRNAs were submitted for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. P < 0.05 was considered as the criterion for significantly enriched GO terms and KEGG pathways.
Construction of ceRNA networks
A previous study by Deng et al.  reported 123 differentially expressed lncRNAs (DElncRNAs) between exosomes isolated from the plasma of sensitive and resistant groups through lncRNA sequencing. Here, we combined the sequencing results of lncRNAs and miRNAs to construct ceRNA networks. First, an online tool miRanda (http://cbio.mskcc.org/miRNA2003/miranda.html) was used to screen the lncRNA-miRNA pair. Subsequently, the RNAhybrid (http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/) and miRanda (http://www.microrna.org/microrna/home.do) databases were used to predict the target genes of the aforementioned DEmiRNAs screened from the lncRNA-miRNA co-expression pairs, and the co-expression pairs of miRNA-mRNA with the opposite expression direction were retained. After that, ceRNA networks were established using lncRNA-miRNA-mRNA co-expression pairs and visualized using Cytoscape (version 3.7.0). Because exosomal lncRNA MSTRG.292666.16 has been reported to be associated with osimertinib resistance in NSCLC , we chose a ceRNA relationship of lncRNA MSTRG.292666.16 -miRNA hsa-miR-6836-5p-mRNA MAPK8IP3 for further experiments.
Fluorescence in Situ Hybridization (FISH)
The localization of lncRNA MSTRG.292666.16 in H1975 cells was determined using a FISH assay kit (Guangzhou RiboBio Co., Ltd, Guangzhou, China) following the manufacturer’s protocols. The H1975 cells were seeded in a 24-well plate at a density of 5 × 104 cells/well and cultured to a confluence of approximately 70%. The cells were fixed with 4% paraformaldehyde for 10 min at room temperature (RT). After washing with PBS, 1 mL of pre-cooled transparent liquid (PBS containing 0.5% Triton X-100) was added to each well and incubated at 4 °C for 5 min. After removing the transparent liquid and a washing step, a prehybridization solution (200 μL) was added and blocked at 37 °C for 30 min. After removing the prehybridization solution, 2.5 μL of 20 μM lncRNA MSTRG.292666.16 FISH probe were mixed into 100 μL hybrid solution, the mixture was added to the cells and incubated in the dark at 37 °C overnight. In the dark at 42 °C, the cells were washed with hybridization solutions I, II, III, and PBS. Finally, the cells were sealed with DAPI, and a fluorescence microscope (IX70, Olympus Corporation, Japan) was used to observe the images.
Our previous in vitro study found that exosomal lncRNA MSTRG.292666.16 was closely associated with osimertinib resistance in NSCLC. Therefore, in this study, we mainly explored its roles in vivo. We first constructed H1975 cells overexpressing lncRNA MSTRG.292666.16, as previously described . Briefly, the pCDH vector and pCDH-MSTRG.292666.16 were prepared and provided by Yanzai Biotechnology (Shanghai) Co. Ltd. The H1975 cells were seeded into a 6-well plate at a density of 4 × 105 cells/well and cultured overnight. The next day, the medium was changed to serum-free medium, and the cells were transfected with 3 μg pCDH vector and 3 μg pCDH-MSTRG.292666.16 using Lipofectamine 3000 (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. After transfected for 6 h, the medium was replaced with complete medium and cultured for another 48. The level of MSTRG.292666.16 in the cells with different treatments was measured by RT-qPCR to evaluate the cell transfection efficiency. The sequences of MSTRG.292666.16 are shown in Table 2.
Dual-luciferase reporter gene assay
The sequences of WT/MUT MSTRG.292666.16 lncRNA inserts, miR-6836-5p mimics/negative control (NC) mimics, and WT/MUT MAPK8IP3 3′-untranslated regions (3′-UTRs) were synthesized and purchased from Yanzai Biotechnology (Shanghai) Co. Ltd. The psiCHECK2 vector (Yanzai Biotechnology (Shanghai) Co. Ltd) and pGL3 vector (Yanzai biotechnology (Shanghai) Co. Ltd) were used to generate the psiCHECK-MSTRG.292666.16 reporter plasmid (psiCHECK-MSTRG.292666.16 WT/MUT) and the 3′-UTR MAPK8IP3 reporter plasmid (pGL3-MAPK8IP3 WT/MUT), respectively. 293 T cells (Cell Bank, Chinese Academy of Sciences) were seeded into a 96-well plate and cultured overnight. The following day, the psiCHECK2-MSTRG.292666.16 WT/MUT (0.3 μg) or psiCHECK vector (0.3 μg), and pGL3-MAPK8IP3 WT/MUT (0.3 μg) or pGL3 vector (0.3 μg) were co-transfected into 293 T cells with NC mimics (100 nM) or miR-6836-5p mimics (100 nM) using Lipofectamine 3000 (Thermo Fisher Scientific) based on the manufacturer’s recommendations. After transfection for 6 h and culture for another 48 h, the cells were harvested to determine luciferase activity using a dual luciferase reporter system (Promega, Madison, WI, USA).
Total RNA was isolated from the plasma and M2 type TAMs using the mirVANA miRNA Isolation kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. The quality and concentration of total RNA was assessed using a microplate reader. Total RNA was reverse transcribed into cDNA using the PrimeScriptTM II 1st Strand cDNA synthesis Kit (Takara Biomedical Technology Co., Ltd., Beijing, China), according to the manufacturer’s protocols. For miRNA, the stem loop RT-qPCR method was used, and U6 was used as a reference gene . For lncRNA and mRNA, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) served as a housekeeping gene. The RT-qPCR reaction was initiated at 95 °C for 3 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. The levels of lncRNA MSTRG.292666.16, miR-6836-5p, MAPK8IP3, and TAM-related markers (CD206, CD163, TGF-β, IL-10, and Arg-1) were calculated using the 2−ΔΔCt method. The sequences of all primers are listed in Table 2.
Before animal experiments, H1975 cells were divided into four groups: control, H1975 + M2 type TAM-derived exosomes, H1975 + pCDH, and H1975 + pCDH-lncRNA groups. The cells in the H1975 + M2 type TAM-derived exosomes group were treated with 10 μg/mL M2 type TAM-derived exosomes, and the cells in the H1975 + pCDH and H1975 + pCDH-lncRNA groups were transfected with pCDH and pCDH-lncRNA, respectively. No treatment was applied for the cells in the control group. After that, the medium was removed, and the final concentration of cells in the different groups was adjusted to 2 × 107 cells/mL with PBS.
Twenty SPF BALB/C female nude mice aged 4–6 weeks were purchased from Shanghai SLAC Laboratory Animal Co., Ltd. (Shanghai, China). All mice were maintained under controlled temperature (24 ± 2 °C) and humidity (50 ± 5%) conditions, with a 12 h light/dark cycle. The rats were given ad libitum access to food and water during the experiment. After 7 days of acclimatization, the mice were randomly divided into four groups (n = 5 for each group): control, M2 TAM-derived Exos, pCDH, and pCDH-lncRNA groups. The mice in the M2 TAM-derived Exos, pCDH, and pCDH-lncRNA groups were subcutaneously injected with 0.1 mL H1975 cells pre-treated with M2 type TAM-derived exosomes (2 × 107 cells/mL), 0.1 mL H1975 cells transfected with pCDH (2 × 107 cells/mL) and 0.1 mL H1975 cells transfected with pCDH-lncRNA (2 × 107 cells/mL), respectively. The mice in the control group were injected with PBS. On the 10th day of modeling, osimertinib (2.5 mg/kg/d) was intraperitoneally injected into each mouse for four weeks. The long and short diameters of the tumors were measured weekly to calculate the tumor volume. After four weeks, the mice were sacrificed by cervical dislocation, and the tumor was removed, weighed, and photographed. All animal experiments were performed in accordance with the National Medical Advisory Committee (NMAC) guidelines using approved procedures of the Institutional Animal Care and Use Committee at Shanghai Pulmonary Hospital.
Data are expressed as mean ± standard deviation (SD), and GraphPad Prism 5 (GraphPad Software, CA, USA) was used for all statistical analyses. Student’s t-test was used for comparisons between the two groups. One-way analysis of variance (ANOVA) followed by Tukey’s test was used to compare the differences among more than two groups. Statistical significance was set at P < 0.05.
Induction of M2 phenotype in TAM
PMA, IL-4, and IL-13 were used to induce THP-1 cells to form M2 type TAMs, and the markers of M2 type macrophage (CD206, CD163, TGF-β, IL-10, and Arg-1) were measured. It was found that the expression levels of CD206, CD163, TGF-β, IL-10, and Arg-1 were significantly increased in the M2 type TAMs compared with those in THP-1 cells (P < 0.05, Fig. 1), which indicated that M2 type TAMs were successfully induced by PMA, IL-4, and IL-13.
Characterization of exosomes isolated from M2 type exosomes
TEM, NTA, and western blot were used to identify the exosomes from THP-1 cells and successfully induced M2 phenotype in TAMs. TEM results showed that the exosomes isolated from the TPH-1 cells and M2 type TAMs were nearly round with a diameter of approximately 100 nm (Fig. 2A). As shown in Fig. 2B, the major peaks of the exosomes isolated from the TPH-1 cells and M2 type TAMs were respectively 142 nm and 131 nm, and their average peaks were respectively 141.1 ± 48.3 nm and 128.6 ± 36.8 nm. These results were in accordance with the size distribution of exosomes, as previously reported . Finally, western blot showed that CD9, HSP70, and TSG101, which are exosome-specific markers, were all expressed in THP-1-derived exosomes and M2 type TAM-derived exosomes (Fig. 2C). These results indicate that exosomes were successfully isolated from THP-1 cells and M2 type TAMs.
The effects of M2 type TAM-derived exosomal MSTRG.292666.16 on osimertinib resistance in vivo
H1975 cells overexpressing lncRNA MSTRG.292666.16 were constructed to investigate the roles of exosomal lncRNA MSTRG.292666.16 in osimertinib-resistant NSCLC. No significant difference was found in the relative level of lncRNA MSTRG.292666.16 between the control and pCDH groups (P = 0.9560 > 0.05, Fig. 3A). The relative level of lncRNA MSTRG.292666.16 in the pCDH-lncRNA group was significantly was higher than that in the control group (P = 0.0018, Fig. 3A). These results suggest that H1975 cells overexpressing lncRNA MSTRG.292666.16 have been successfully established. Then, the effects of exosomal MSTRG.292666.16 on osimertinib resistance were determined by subcutaneously injecting the cells with different treatments into nude mice in vivo. After tumor formation, nude mice were administered osimertinib for four weeks. With the increase in time, there was no significant difference in tumor volume between the control and pCDH groups (P > 0.05, Fig. 3B). However, the tumor volume in the M2 TAM-derived Exos and pCDH-lncRNA groups significantly increased with time (P < 0.05, Fig. 3B). After four weeks, the tumor volumes in the different groups were measured. The tumor volume in the control and pCDH groups was respectively 277.53 ± 43.73 mm3 and 322.48 ± 85.21 mm3, which showed no significant difference between the two groups (P = 0.3642 > 0.05, Fig. 3C, D). Compared with the control group, the tumor volume was evidently increased after M2 type TAM-derived exosomes (P = 0.0003 < 0.05) and lncRNA MSTRG.292666.16 overexpression (P < 0.0001) pre-treatment (Fig. 3C, D). These in vivo results strongly suggest that M2 type TAM-derived exosomes and lncRNA MSTRG.292666.16 overexpression could enhance osimertinib resistance in NSCLC.
M2 type TAM-related markers in the plasma, and the localization of MSTRG.292666.16 in cells
To understand the mechanisms by which M2 type TAM-derived exosomes promote drug resistance through lncRNA MSTRG.292666.16, M2 type TAM-related markers were determined in the sensitive and resistant plasma. Compared to the sensitive plasma, the relative expression of CD206, CD163, TGF-β, and IL-10 in the resistant plasma was markedly upregulated (P < 0.05, Fig. 4A). These results showed the existence of M2 type TAM in the osimertinib-resistant plasma. In addition, FISH results showed that lncRNA MSTRG.292666.16 was located in the cytoplasm of H1975 cells (Fig. 4B). Therefore, we further investigated the ceRNA mechanism of M2 type TAM-derived exosomal lncRNA MSTRG.292666.16 in osimertinib resistance in NSCLC in the following study.
Screening DEmiRNAs in the plasma-derived exosomes and functional analyses
Exosomes isolated from the sensitive plasma and resistant plasma were identified previously  and then sent for miRNA sequencing. Upon sequencing, 680 DEmiRNAs were identified between the exosomes isolated from the sensitive and resistant plasma, including 393 upregulated miRNAs and 287 downregulated miRNAs in the resistant plasma-derived exosomes (Fig. 5A, Additional file 1: Table S1). Then, these DEmiRNAs were subjected to GO terms and KEGG pathway enrichment analysis. It was found that these DEmiRNAs were related to GO terms of “cell adhesion”, “cell migration”, “nervous system development”, “intracellular dephosphorylation” and “neurotrophin TRK receptor signaling pathway” (Fig. 5B). Besides, these DEmiRNAs were significantly enriched in pathways of “Wnt signaling pathway”, “Rap1 signaling pathway”, “Calcium signaling pathway”, “Cgmp-PKG signaling pathway”, “Ras signaling pathway”, “PI3K-Akt signaling pathway”, “HIF-1 signaling pathway” and “MAPK signaling pathway” (Fig. 5C).
Analysis of ceRNA networks
The analysis of ceRNA network was performed including the previous lncRNA sequencing data from the sensitive and resistant plasma . In total, 103 pairs of DElncRNA-DEmiRNA and 181 pairs of DEmiRNA-DEmRNA were obtained, and ceRNA networks were generated using DElncRNA-DEmiRNA-DEmRNA pairs (Fig. 6). Afterwards, the genes in the ceRNA networks were used for functional analyses. Based on the threshold of P < 0.05, 8 GO terms of biology process (“transcription, DNA-templated”, “positive/negative regulation of transcription from RNA polymerase II promoter”, “positive regulation of transcription, DNA-template”, “response to cytokine”, “hippocampus development”, “bone morphogenesis” and “negative regulation of endothelial cell proliferation”) and 7 KEGG pathways (“biosynthesis of antibiotics”, “terpenoid backbone biosynthesis”, “carbon metabolism”, “lysine degradation”, “MAPK signaling pathway”, “serotonergic synapse” and “lysosome”) were significantly enriched (Table 3).
Additionally, it has been reported that the MAPK signaling pathway plays an important role in cancer drug resistance mechanisms . Our study also found that MAPK signaling pathway was closely related to osimertinib resistance of NSCLC, and RPS6KA1, CACNA1B and MAPK8IP3 were significantly enriched in “MAPK signaling pathway” (Table 3). Based on the ceRNA networks, we found that lncRNA MSTRG.292666.16 was the hub lncRNA, and MAPK8IP3 was regulated by lncRNA MSTRG.292666.16/hsa-miR-6836-5p axis (Fig. 6). Therefore, we chose a ceRNA network of lncRNA MSTRG.292666.16-hsa-miR-6836-5p-MAPK8IP3 for further verification.
Validation of lncRNA MSTRG.292666.16-hsa-miR-6836-5p-mRNA-MAPK8IP3 in the exosomes
Compared with the sensitive plasma-derived exosomes, the levels of MSTRG.292666.16 and MAPK8IP3 were significantly increased in the resistant plasma-derived exosomes (P < 0.0001, Fig. 7A, C). The miR-6836-5p level in the exosomes isolated from sensitive and resistant plasma was opposite to that of MSTRG.292666.16 and MAPK8IP3 levels (Fig. 7B). For the exosomes isolated from the cells, the levels of MSTRG.292666.16 in the THP-1-derived exosomes and M2 type TAM-derived exosomes were 1.00 ± 0.043 and 29.26 ± 3.556, respectively, which showed a significant increase in the M2 type TAM-derived exosomes (P = 0.0002, Fig. 7D). A significant reduction in miR-6836-5p level in the M2 type TAM-derived exosomes compared with that in THP-1-derived exosomes was observed (P = 0.0063, Fig. 7E). The trend of MAPK8IP3 expression in the exosomes derived from different cells was similar to that of MSTRG.292666.16 (Fig. 7F). The level trends of MSTRG.292666.16, miR-6836-5p, and MAPK8IP3 in the exosomes determined by RT-qPCR were consistent with those measured by sequencing, which verified the relatively reliable sequencing results.
MSTRG.292666.16 interacted with miR-6836-5p and MAPK8IP3
Finally, a dual-luciferase reporter gene assay was performed to analyze the relationship between MSTRG.292666.16, miR-6836-5p, and MAPK8IP3. In the control group including psiCHECK2-control and pGL3-control plasmids, no significant differences were observed in luciferase activity between NC mimics and hsa-miR-6836-5p mimics (P < 0.05, Fig. 8G, H). In the psiCHECK2-lncRNA WT plasmid, the relative luciferase activity after co-transfected with miR-6836-5p mimics was signally decreased compared with NC mimics (P = 0.0022); while in the psiCHECK2-lncRNA MUT plasmid, the luciferase activity was significantly increased (P = 0.0105 < 0.05, Fig. 8G). These results indicate that MSTRG.292666.16 is the target of miR-6836-5p. Additionally, in the pGL3-MAPK8IP3 WT plasmid, the relative luciferase activity in the miR-6838-5p mimics group was significantly lower than that in the NC mimics group (P = 0.0001); however, after MAPK8IP3 mutation, the relative luciferase activity was significantly enhanced (P = 0.0037 < 0.05, Fig. 6H). These results indicate that MAPK8IP3 is the target of miR-6836-5p. Taken together, lncRNA MSTRG.292666.16, hsa-miR-6836-5p, and mRNA MAPK8IP3 could interact with each other.
Lung cancer, especially NSCLC, remains the leading cause of cancer death, and osimertinib is an important first-line drug for EGFR mutation-positive NSCLC . However, with the increasing time of use, the vast majority of patients develop acquired resistance to osimertinib, seriously affecting the survival of patients. Previous studies have reported that acquired EGFR p.T790M (c.2369 C > T; p.Thr790Met) and p.C797S (c.2389 T > A and c.2390G > C, p.Cys797Ser) mutations are closely associated with osimertinib resistance in NSCLC [28, 29]. The EGFR C797S mutation is the most common tertiary EGFR mutation, and cysteine is replaced by serine at codon 797 within the ATP-binding site in this mutation, resulting in the loss of covalent bonds between osimertinib and mutant EGFR . It is predictable that the C797S mutation can also develop cross-resistance to other irreversible EGFR-TKIs by preventing their binding to the EGFR active site . However, the other mechanisms related to osimertinib resistance in NSCLC remain unclear, and novel therapeutic strategies are urgent to be discovered.
Exosomes serve as tools for intercellular communication, and Wu et al. reported that exosomes derived from EGFR wild-type NSCLC could induce osimertinib resistance after being transferred to cancer cells . Our previous study found that osimertinib-resistant exosomes could induce H19575 cells resistant to osimertinib, and exosomal lncRNA MSTRG.292666.16 may be associated with osimertinib resistance in NSCLC . Macrophages can promote important homeostasis functions, such as endocytosis, phagocytosis, and signaling during inflammation . Therefore, this study proposed a hypothesis that lncRNA MSTRG.292666.16 may be transferred from M2 type TAMs of osimertinib-resistant lung cancer, and M2 type TAM-derived exosomes may influence osimertinib resistance in NSCLC through MSTRG.292666.16. CD163 and CD206 are both expressed by macrophages and are markers of macrophage activation . The expression of CD163 can be increased by IL-10 stimulation, and CD206 expression can be upregulated by IL-4 and IL-13 . TGF-β, IL-10, IL-4, and IL-13 are secreted by M2 type macrophages. Besides, M2 type TAMs have been reported to produce Arg-1 and promote wound healing by increasing the generation of Arg-1 and pro-chondrogenic cytokines (IL-10 and TGF-β) . Our study determined the expression of CD163, CD206, TGF-β, IL-10, and Arg-1 in the plasma of NSCLC patients and found that their expression was significantly upregulated in osimertinib-resistant plasma, which indicated the existence of M2 type TAMs in osimertinib-resistant plasma.
Thereafter, M2 type TAMs were successfully established, and exosomes were isolated from these cells. Tumor-forming experiments in nude mice showed that M2 type TAM-derived exosomes and MSTRG.292666.16 overexpression could enhance osimertinib resistance in NSCLC. A previous study demonstrated that lncRNA ZFAS1 was highly expressed in the tissues of esophageal squamous cell carcinoma, and exosomes with ZFAS1 overexpression significantly promoted cell-based xenograft tumor growth . Another study by Qu et al. showed that high levels of lncRNA ARSR were observed in sunitinib-resistant renal cell carcinoma cells, and locked nucleic acid (LNA)-mediated ARSR knockdown effectively restored the sensitivity of xenografts to synchronous sunitinib therapy . Furthermore, our results also showed that the MSTRG.292666.16 level was significantly higher in the M2 type TAM-derived exosomes than in the THP-1-derived exosomes. Together with our results, it can be inferred that M2 type TAM-derived exosomes may improve osimertinib resistance in NSCLC by upregulating MSTRG.292666.16.
Previous studies have shown that ceRNAs influence the progression of liver cancer, and cytoplasmic lncRNAs can function as ceRNAs regulating mRNA stabilization or translation by sponging miRNAs, and thus affecting signaling pathways [38, 39]. In this study, FISH results showed that lncRNA MSTRG.292666.16 localized in the cytoplasm of H1975 cells, which suggested its potential involvement in ceRNA mechanisms. Small sequencing identified 680 DEmiRNAs, including 393 upregulation and 287 downregulation in exosomes from resistant plasma. Then, by combining the miRNA sequencing data with our previous lncRNA sequencing data, ceRNA networks were constructed, and generated 103 pairs of DElncRNAs-DEmiRNAs and 181 pairs of DEmiRNAs-DEmRNAs. The genes in the ceRNA networks were significantly enriched in 8 GO terms and 7 KEGG pathways, such as “biosynthesis of antibiotics” (ICMT, OGDH, TKT, ZMPSTE24), “terpenoid backbone biosynthesis” (ICMT, ZMPSTE24), and “MAPK signaling pathway” (RPS6KA1, CACNA1B, MAPK8IP3). ICMT has been reported to facilitate the growth, migration, and doxorubicin resistance of hepatocellular carcinoma via various carcinogenic pathways . OGDH, a subunit of the OGDH complex, is involved in the TCA cycle, energy balance, and cellular signaling. Lu et al.  showed that OGDH, which serves as a positive regulator of gastric cancer progression, enhances mitochondrial function, and activates the Wnt/β-catenin signaling pathway. TKT was upregulated in the NSCLC cells, and its knockdown strengthened the action of gefitinib in NSCLC cells . Loss of ZMPSTE24 accelerates senescence associated with cancer . The MAPK signaling pathway has been shown to be activated in many cancers that are resistant to drugs , and our study found that RPS6KA1, CACNA1B, and MAPK8IP3 were significantly enriched in the MAPK signaling pathway. Fan et al.  demonstrated that the octreotide-paclitaxel conjugate could reverse paclitaxel resistance in A2780/Taxol human ovarian cancer cells by repressing the activity of the MAPK signaling pathway. These reports together with our results support that ICMT, OGDH, TKT, ZMPSTE24, RPS6KA1, CACNA1B, and MAPK8IP3, as well as pathways for biosynthesis of antibiotics, terpenoid backbone biosynthesis, and MAPK signaling pathway may play important roles in osimertinib resistance in NSCLC.
The MAPK signaling pathway has been reported to be activated in cancer drug resistant samples, and its abnormal activation could contribute to cell loss of differentiation and apoptosis , and could be a mechanism of resistance to osimertinib in EGFR-mutated NSCLC . Our established ceRNA networks demonstrated that RPS6KA1, CACNA1B and MAPK8IP3 were significantly enriched in “MAPK signaling pathway”, and lncRNA MSTRG.292666.16 was the hub node. Therefore, a ceRNA network of lncRNA MSTRG.292666.16- miR-6836-5p-MAPK8IP3 was chosen for further validation. The levels of MSTRG.292666.16 and MAPK8IP3 were significantly increased in resistant plasma-derived exosomes and M2 type TAM-derived exosomes, while miR-6836-5p levels were significantly evidently reduced, which was consistent with the expression trend of ceRNA mechanisms. The dual-luciferase reporter gene assay indicated that MSTRG.292666.16, miR-6836-5p, and MAPK8IP3 could interact with each other. A previous study reported that HOTAIR was the target of miR-331-3p, and HER2 was the target of miR-331-3p through a luciferase assay, which suggested that lncRNA HOTAIR could regulate HER2 expression by sponging miR-331-3p in gastric cancer . Another study found that lncRNA UCA1, miR-143, and FOSL2 interacted with each other in a luciferase reporter assay, and UCA1 packaged into exosomes could promote gefitinib resistance in NSCLC cells by mediating miR-143/FOSL2 . Therefore, we speculated that M2 type TAM-derived exosomes may enhance osimertinib resistance in NSCLC by regulating the MSTRG.292666.16/miR-6386-5p/MAPK8IP3 axis.
However, this study has some limitations. First, the specific ceRNA mechanism of MSTRG.292666.16/miR-6836-5p/MAPK8IP3 through M2 type TAM-derived exosomes in osimertinib-resistant NSCLC needs to be further investigated in vitro. Moreover, a larger sample size is needed to verify our conclusions, and further experiments should be conducted to investigate the effects of the identified genes and pathways on osimertinib resistance in NSCLC.
In conclusion, through sequencing, it was found that ICMT, OGDH, TKT, ZMPSTE24, RPS6KA1, CACNA1B, and MAPK8IP3, as well as pathways for biosynthesis of antibiotics, terpenoid backbone biosynthesis, and MAPK signaling pathway may play important roles in osimertinib resistance in NSCLC. Additionally, M2 type TAM-derived exosomes may promote osimertinib resistance in NSCLC via the MSTRG.292666.16/miR-6836-5p/MAPK8IP3 axis. Our work illuminated a novel mechanism of osimertinib resistance in NSCLC, and implied that M2 type TAM-derived exosomal MSTRG.292666.16/ miR-6836-5p/ MAPK8IP3 serves as a novel therapeutic target and may have potential therapeutic value in treating osimertinib resistance in NSCLC.
Availability of data and materials
The dataset used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Andrews Wright NM, Goss GD. Third-generation epidermal growth factor receptor tyrosine kinase inhibitors for the treatment of non-small cell lung cancer. Transl Lung Cancer Res. 2019;8(Suppl 3):S247–64.
Duma N, Santana-Davila R, Molina JR. Non-small cell lung cancer: epidemiology, screening, diagnosis, and treatment. Mayo Clin Proc. 2019;94(8):1623–40.
Pallis AG, Syrigos KN. Epidermal growth factor receptor tyrosine kinase inhibitors in the treatment of NSCLC. Lung Cancer. 2013;80(2):120–30.
Tan CS, Kumarakulasinghe NB, Huang YQ, Ang YLE, Choo JR, Goh BC, Soo RA. Third generation EGFR TKIs: current data and future directions. Mol Cancer. 2018;17(1):29.
Ramalingam SS, Vansteenkiste J, Planchard D, Cho BC, Gray JE, Ohe Y, Zhou C, Reungwetwattana T, Cheng Y, Chewaskulyong B, et al. Overall survival with osimertinib in untreated, EGFR-mutated advanced NSCLC. N Engl J Med. 2020;382(1):41–50.
Wu YL, Tsuboi M, He J, John T, Grohe C, Majem M, Goldman JW, Laktionov K, Kim SW, Kato T, et al. Osimertinib in resected EGFR-mutated non-small-cell lung cancer. N Engl J Med. 2020;383(18):1711–23.
Kalluri R, LeBleu VS. The biology, function, and biomedical applications of exosomes. Science. 2020;367(6478):eaau6977.
Li I, Nabet BY. Exosomes in the tumor microenvironment as mediators of cancer therapy resistance. Mol Cancer. 2019;18(1):32.
Wu S, Luo M, To KKW, Zhang J, Su C, Zhang H, An S, Wang F, Chen D, Fu L. Intercellular transfer of exosomal wild type EGFR triggers osimertinib resistance in non-small cell lung cancer. Mol Cancer. 2021;20(1):17.
Pathria P, Louis TL, Varner JA. Targeting tumor-associated macrophages in cancer. Trends Immunol. 2019;40(4):310–27.
Zhang Q, Xia J, Wang Y, Zhang J, Ji C, Cong R, Wang Y, Song N. Tumor infiltrating M2 macrophages could predict biochemical recurrence of localized prostate cancer after radical prostatectomy. Exp Cell Res. 2019;384(1):111588.
Nowak M, Klink M. The role of tumor-associated macrophages in the progression and chemoresistance of ovarian cancer. Cells. 2020;9(5):1299.
Zheng P, Chen L, Yuan X, Luo Q, Liu Y, Xie G, Ma Y, Shen L. Exosomal transfer of tumor-associated macrophage-derived miR-21 confers cisplatin resistance in gastric cancer cells. J Exp Clin Cancer Res. 2017;36(1):53.
Braga EA, Fridman MV, Moscovtsev AA, Filippova EA, Dmitriev AA, Kushlinskii NE. LncRNAs in ovarian cancer progression, metastasis, and main pathways: ceRNA and alternative mechanisms. Int J Mol Sci. 2020;21(22):8855.
Qi X, Zhang DH, Wu N, Xiao JH, Wang X, Ma W. ceRNA in cancer: possible functions and clinical implications. J Med Genet. 2015;52(10):710–8.
Xu G, Zhu Y, Liu H, Liu Y, Zhang X. LncRNA MIR194-2HG promotes cell proliferation and metastasis via regulation of miR-1207-5p/TCF19/Wnt/beta-catenin signaling in liver cancer. Onco Targets Ther. 2020;13:9887–99.
Deng Q, Fang Q, Xie B, Sun H, Bao Y, Zhou S. Exosomal long non-coding RNA MSTRG.292666.16 is associated with osimertinib (AZD9291) resistance in non-small cell lung cancer. Aging (Albany NY). 2020;12(9):8001–15.
Genin M, Clement F, Fattaccioli A, Raes M, Michiels C. M1 and M2 macrophages derived from THP-1 cells differentially modulate the response of cancer cells to etoposide. BMC Cancer. 2015;15:577.
Thery C, Amigorena S, Raposo G, Clayton A. Isolation and characterization of exosomes from cell culture supernatants and biological fluids. Curr Protoc Cell Biol. 2006; Chapter 3:Unit 3 22.
Soares Martins T, Catita J, Rosa MI, Silva OABCE, Henriques AG. Exosome isolation from distinct biofluids using precipitation and column-based approaches. PLoS One. 2018;13(6):e0198820.
Zhu Q, Li Q, Niu X, Zhang G, Ling X, Zhang J, Wang Y, Deng Z. Extracellular vesicles secreted by human urine-derived stem cells promote ischemia repair in a mouse model of hind-limb ischemia. Cell Physiol Biochem. 2018;47(3):1181–92.
Yin J, Zeng A, Zhang Z, Shi Z, Yan W, You Y. Exosomal transfer of miR-1238 contributes to temozolomide-resistance in glioblastoma. EBioMedicine. 2019;42:238–51.
Park ST, Kim J. Trends in next-generation sequencing and a new era for whole genome sequencing. Int Neurourol J. 2016;20(Suppl 2):S76-83.
Lin Y, Jin Y, Xu T, Zhou S, Cui M. MicroRNA-215 targets NOB1 and inhibits growth and invasion of epithelial ovarian cancer. Am J Transl Res. 2017;9(2):466–77.
Xu G, Ao R, Zhi Z, Jia J, Yu B. miR-21 and miR-19b delivered by hMSC-derived EVs regulate the apoptosis and differentiation of neurons in patients with spinal cord injury. J Cell Physiol. 2019;234(7):10205–17.
Lee S, Rauch J, Kolch W. Targeting MAPK signaling in cancer: mechanisms of drug resistance and sensitivity. Int J Mol Sci. 2020;21(3):1102.
Leonetti A, Sharma S, Minari R, Perego P, Giovannetti E, Tiseo M. Resistance mechanisms to osimertinib in EGFR-mutated non-small cell lung cancer. Br J Cancer. 2019;121(9):725–37.
Castellanos E, Feld E, Horn L. Driven by mutations: the predictive value of mutation subtype in EGFR-mutated non-small cell lung cancer. J Thorac Oncol. 2017;12(4):612–23.
Harrison PT, Vyse S, Huang PH. Rare epidermal growth factor receptor (EGFR) mutations in non-small cell lung cancer. Semin Cancer Biol. 2020;61:167–79.
Yoneda K, Imanishi N, Ichiki Y, Tanaka F. Treatment of non-small cell lung cancer with EGFR-mutations. J UOEH. 2019;41(2):153–63.
Veloso P, Fernandez A, Terraza-Aguirre C, Alvarez C, Vernal R, Escobar A, Hernandez M. Macrophages skew towards M1 profile through reduced CD163 expression in symptomatic apical periodontitis. Clin Oral Investig. 2020;24(12):4571–81.
Nielsen MC, Hvidbjerg Gantzel R, Claria J, Trebicka J, Moller HJ, Gronbaek H. Macrophage activation markers, CD163 and CD206, in acute-on-chronic liver failure. Cells. 2020;9(5):1175.
Mosser DM, Edwards JP. Exploring the full spectrum of macrophage activation. Nat Rev Immunol. 2008;8(12):958–69.
Fernandes TL, Gomoll AH, Lattermann C, Hernandez AJ, Bueno DF, Amano MT. Macrophage: a potential target on cartilage regeneration. Front Immunol. 2020;11:111.
Li Z, Qin X, Bian W, Li Y, Shan B, Yao Z, Li S. Exosomal lncRNA ZFAS1 regulates esophageal squamous cell carcinoma cell proliferation, invasion, migration and apoptosis via microRNA-124/STAT3 axis. J Exp Clin Cancer Res. 2019;38(1):477.
Qu L, Ding J, Chen C, Wu ZJ, Liu B, Gao Y, Chen W, Liu F, Sun W, Li XF, et al. Exosome-transmitted lncARSR promotes sunitinib resistance in renal cancer by acting as a competing endogenous RNA. Cancer Cell. 2016;29(5):653–68.
Chen DL, Lu YX, Zhang JX, Wei XL, Wang F, Zeng ZL, Pan ZZ, Yuan YF, Wang FH, Pelicano H, et al. Long non-coding RNA UICLM promotes colorectal cancer liver metastasis by acting as a ceRNA for microRNA-215 to regulate ZEB2 expression. Theranostics. 2017;7(19):4836–49.
Schmitt AM, Chang HY. Long noncoding RNAs in cancer pathways. Cancer Cell. 2016;29(4):452–63.
Xu J, Zhu Y, Wang F, Zhou Y, Xia G, Xu W. ICMT contributes to hepatocellular carcinoma growth, survival, migration and chemoresistance via multiple oncogenic pathways. Biochem Biophys Res Commun. 2019;518(3):584–9.
Lu X, Wu N, Yang W, Sun J, Yan K, Wu J. OGDH promotes the progression of gastric cancer by regulating mitochondrial bioenergetics and Wnt/beta-catenin signal pathway. Onco Targets Ther. 2019;12:7489–500.
Cao L, Hong W, Cai P, Xu C, Bai X, Zhao Z, Huang M, Jin J. Cryptotanshinone strengthens the effect of gefitinib against non-small cell lung cancer through inhibiting transketolase. Eur J Pharmacol. 2021;890:173647.
Lee JH, Yoo NJ, Kim MS, An CH, Lee SH. Mutational and expressional alterations of ZMPSTE24, DNA damage response-related gene, in gastric and colorectal cancers. Pathol Res Pract. 2016;212(12):1113–8.
Fan LL, Chen X, Zhang XY, Li ZM, Fan XM, Shen Y. Octreotide-paclitaxel conjugate reverses paclitaxel resistance by p38 mitogen-activated protein kinase (MAPK) signaling pathway in A2780/taxol human ovarian cancer cells. Med Sci Monit. 2020;26:e922612.
Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, Lu KH, Li JH, et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 2014;13:92.
Chen X, Wang Z, Tong F, Dong X, Wu G, Zhang R. lncRNA UCA1 promotes Gefitinib resistance as a ceRNA to target FOSL2 by sponging miR-143 in non-small cell lung cancer. Mol Ther Nucleic Acids. 2020;19:643–53.
This work was supported by Key projects of Shanghai Zhangjiang High Tech Park (2016-08).
Ethics approval and consent to participate
All subjects gave their informed consent for inclusion before they participated in the study. All the animal experiments were performed in line with the National Medical Advisory Committee (NMAC) guidelines using approved procedures of the Institutional Animal Care and Use committee at Shanghai Pulmonary Hospital.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
All differentially expressed miRNAs between the exosomes isolated from the sensitive and resistant plasma with the thresholds of |log2Fold change (FC)|> 1 and false discovery rate (FDR) < 0.05.
About this article
Cite this article
Wan, X., Xie, B., Sun, H. et al. Exosomes derived from M2 type tumor-associated macrophages promote osimertinib resistance in non-small cell lung cancer through MSTRG.292666.16-miR-6836-5p-MAPK8IP3 axis. Cancer Cell Int 22, 83 (2022). https://doi.org/10.1186/s12935-022-02509-x