Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 01 April 2022
Sec. Pharmacology of Anti-Cancer Drugs
This article is part of the Research Topic Multi-omics in Studying the Mechanisms of Anti-cancer Drugs Resistance and Toxicity View all 13 articles

MMP9 Expression Correlates With Cisplatin Resistance in Small Cell Lung Cancer Patients

  • 1Department of Oncology, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
  • 2Department of Neurology, Ganzhou People’s Hospital, Ganzhou, China
  • 3Department of Neurology, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China

Background: Cisplatin is the basis of the primary treatment for SCLC chemotherapy. However, the limited objective response rate and definite drug resistance greatly restrict the clinical potential and therapeutic benefits of cisplatin use. Therefore, it is essential to identify biomarkers that can discern the sensitivity of SCLC patients to cisplatin treatment.

Methods: We collected two SCLC cohorts treated with cisplatin that included mutation data, prognosis data and expression data. The sensitivity of cisplatin was evaluated by the pRRophetic algorithm. MCPcounter, quanTIseq, and xCell algorithms were used to evaluate immune cell score. GSEA and ssGSEA algorithms were used to calculate immune-related pathway scores. Univariate and multivariate Cox regression models were employed, and survival analysis was used to evaluate the prognostic value of the candidate genes.

Results: MMP9-High is related to improved clinical prognoses of patients with SCLC (HR = 0.425, p = 0.0085; HR = 0.365, p = 0.0219). Multivariate results showed that MMP-High could be used as an independent predictor of the prognosis of SCLC after cisplatin treatment (HR = 0.216, p = 0.00153; HR = 0.352; p = 0.0199). In addition, MMP9-High displayed a significantly lower IC50 value of cisplatin and higher immunogenicity than MMP9-Low SCLC. Compared with MMP9-Low SCLC, MMP9-High included significantly increased levels of T-cells, cytoxic lymphocytes, B-cells, NK-cells, and dense cells (DCS). Similarly, the activity of cytokine binding, B-cell, NK-cell mediated immune response chemokine binding, and antigen presentation pathways in MMP9-High was significantly higher than that in MMP9-Low.

Conclusion: In this study, we identified that MMP9-High could be potentially considered a novel biomarker used to ascertain the improved prognosis of SCLC patients after cisplatin treatment. Furthermore, we indicated that the tumor immune microenvironment of MMP9-High SCLC is mainly characterized by a large number of infiltrated activated immune cells as well as activated immune-related pathways.

Introduction

Small cell lung cancer (SCLC), which accounts for 13–15% of lung cancer, is a subtype of lung cancer known to exhibit high malignancy and poor prognosis (Sabari et al., 2017; Qiu et al., 2019; Luo et al., 2019). The median survival time of SCLC patients is only 8–13 months, and the 5-year survival rate is 1–5% (de Hoyos and DeCamp, 2014; Carter et al., 2014). Chemotherapy for SCLC is comprised of a combination approach that includes platinum-based antineoplastic drugs commonly used in the treatment of various cancers. Cisplatin a well-known, effective, and widely used first-line drug with an objective response rate (ORR) of 50–60% (Horn et al., 2018). However, almost all patients diagnosed with SCLC will inevitably present with drug resistance and tumor recurrence. Studies have shown that the functional mechanism of cisplatin resistance occurs as a result of structural changes to DNA or cytoplasm, the abnormality of DNA damage repair, the change of signaling triggered by molecular damage caused by cisplatin, and the change of compensatory survival signal (Galluzzi et al., 2014; Inoue et al., 2014; Song et al., 2015). However, the lack of biomarkers used to identify cisplatin sensitivity in SCLC population is a large clinical detriment. Therefore, it is of great importance to find biomarkers that pertain to the sensitivity of SCLC patients undergoing cisplatin treatment.

Evidence has recognized the close relationship between tumor immune microenvironment (TIME) and chemotherapeutic drugs. Tumor-associated macrophages (TAMs), tumor-associated neutrophils (TANs), myeloid-derived suppressor cells (MDSCs), regulatory T-cells, Immunosuppressive cells such as T-regulatory cells (Tregs), and regulatory B-cells (Bregs) can not only directly inhibit killer cells such as CTL and NK, but also interact with each other, enhancing the effect of inhibitory factors and facilitating the recruitment of more immunosuppressive cells. This process enables tumor cells to achieve immune escape (Lin et al., 2019). Additionally, TAMs are known to secrete TGF-β1, which results in the up-regulation of Gfi-1 expression in tumor cells. Gfi-1 expression in the promoter region effectively inhibits the expression of CTGF and HMGB1, which accordingly reduce the sensitivity of tumor cells to gemcitabine (Xian et al., 2017). TAMs expressed IGF act on the IGF1 receptor to promote chemotherapy resistance of gemcitabine and albumin-bound paclitaxel (Ireland et al., 2016). In addition, the study showed that during the chemotherapy of gemcitabine, TAMs, G-MDSCs, Tregs and T-cells decreased, and CTL increased (Mitchem et al., 2013; Eriksson et al., 2016). Gemcitabine promoted the up regulation of HLA-DR, CD40, CCR7 and the down regulation of CD163 and cd206, and induced M1 polarization (Di Caro et al., 2016). In conclusion, chemotherapy drugs can act on TIME, and TIME can also affect the effect of chemotherapy drugs.

Matrix metalloproteinases (MMPs) play an important role in tumorigenesis, development, invasion and metastasis. MMP9 is the largest molecular weighted enzyme in the MMP family. The function of MMP9 is closely related to tumor invasion and metastasis through its functional degradation of type IV and V collagen and gelatin (Mondal et al., 2020). Studies have shown that MMP9 expression is regulated in ovarian cancer, cervical cancer, non-small cell lung cancer and breast cancer, all of which demonstrate a close relationship to cisplatin sensitivity (Rauvala et al., 2006; Braicu et al., 2014; Qiao et al., 2020; Makhoul et al., 2016). In addition, Li et al. found evidence that MMP9 can regulate the biological function of monocytes (Zhou et al., 2012; Li et al., 2021; Xu et al., 2021). Furthermore, Xu et al. found that the high MMP9 group was significantly enriched in the immune response pathway and cytokine production pathway (Xu et al., 2021). However, at present, the expression of MMP9 and cisplatin sensitivity in SCLC patients has not been clarified and the associated relationship with the tumor immune microenvironment is not clear. Therefore, in this study, we aimed to explore the relationship between MMP9 expression and cisplatin sensitivity as well as elucidate the TIME in SCLC, so as to provide a theoretical basis for the precise treatment of SCLC and improve the clinical benefits to patients.

Methods

SCLC Cohort

We collected two published SCLC cohorts from the gene expression omnibus (GEO) database (Clough and Barrett, 2016), namely SCLC (George et al., 2015) and SCLC (Jiang et al., 2016). In this study, the inclusion criteria for SCLC patients were that the included SCLC patients must have survival data, expression data, and mutation data. According to the above inclusion criteria, a total of 68 SCLC patients were recorded from paper 1 (Supplementary Table S1) and 48 SCLC patients from paper 2 (Supplementary Table S2). We filtered the mutation data of the two SCLC cohorts according to the definition and type of non-synonymous mutation in a maftools R package (Mayakonda et al., 2018). The non-synonymous mutation data obtained after filtering was then used for subsequent analyses. The analysis process of this study is detailed in Figure 1.

FIGURE 1
www.frontiersin.org

FIGURE 1. The comprehensive study design.

Prediction of Cisplatin Sensitivity

We used the pRRophetic algorithm (Geeleher et al., 2014) to predict the IC50 value of cisplatin by constructing a relevant ridge regression model with a GDSC cell line expression profile as a training set, and the SCLC cohort as validation set. According to the median value of MMP9 expression, SCLC patients were divided into MMP9-High and MMP9-Low categories. The IC50 values of cisplatin between the two groups were analyzed using the Mann Whitney U test. Further, a Spearman correlation test was used to analyze the expression of MMP9 and the related IC50 value of cisplatin.

Analysis of Tumor Immune Microenvironment Infiltration

MCPcounter, quanTIseq and xCell algorithms were used to analyze the expression profile data of SCLC and obtain the overall score of immune cells (Becht et al., 2016; Aran et al., 2017; Plattner et al., 2020). GSEA was used to analyze the difference of signal pathway activity between MMP9-High and MMP9-Low classes (Reimand et al., 2019). ssGSEA (Hänzelmann et al., 2013) was used to evaluate the signal pathway activity of each SCLC patient in the c2 and c5 pathway sets according to the MsigDB database (Liberzon et al., 2011).

The Predictive Value of MMP9 in an Immunotherapy Cohort

We verified the prognostic value of MMP9 expression in the NSCLC cohort receiving ICIs treatment using a CAMOIP webpage tool (Lin et al., 2021a).

Immunogenicity Analysis

The mutation data and expression data of the TCGA cohort were downloaded using the TCGAbiolinks R package. The data of tumor mutation burden (TMB) and neoantigen loads (NALs) of the TCGA cohort were obtained from published literature (Thorsson et al., 2018).

Immunohistochemistry

Tissue samples were deparaffinized and rehydrated. After treatment with endogenous peroxidase blocking solution, they were treated with specific antibodies against MMP9 (ab76003, Abcam), overnight at 4°C. After they were washed with PBS, the samples were treated with horseradish peroxidase-conjugated anti-rabbit IgG (SV0002, Boster, Wuhan, China) and then stained with diaminobenzidine (DAB). All results were assessed by two pathologists. Expression levels were scored by multiplying the percentage of positive cells by the staining intensity. The positivity percentage was scored as 0 if <5% (negative), 1 if 5–30% (sporadic), 2 if 30–70% (focal) and 3 if 70% (diffuse) of the cells were stained; and staining intensity was scored as 0 for no staining, 1 for weak to moderate staining and 2 for strong staining. A score of ≥2 was regarded as ‘high’ and the score of <2 is regarded as “low” in immunohistochemical staining.

Cell Counting Kit-8 Assay

Cells were cultured at 5 × 103 cells per well in a 96-well plate with cytotoxic drugs for 24 h. Cytotoxic drugs (cisplatin and etoposide) were diluted to obtain different concentration gradients. Absorbance was detected at 450 nm after treatment with 10 μl CCK-8 reagent (Dojindo, Kumamoto) for 4 h. The experiments were performed in five replicate wells per sample and the assays were conducted in triplicate.

Quantitative Real-Time PCR

Total RNA was isolated using Trizol reagent (Invitrogen, United States), according to the manufacturer’s instructions. The quantity and purity of the total RNA was measured using the Nanodrop® ND1000 (Thermo Fisher) and the Agilent Bioanalyzer. Reverse transcription was performed with 2 μg of total RNA using M-MLV reverse transcriptase (Accurate Biology, AG11728) according to the manufacturer’s recommendations. Quantitative PCR was performed using CFX96 Touch Real-Time PCR Detection Instrument (BioRad, United States). Reactions were performed using SYBR® Green Premix Pro Taq HS qPCR Kit (ROX Plus) (AG11718). Values were normalized to GAPDH via the 2−ΔΔCt method.

Statistical Analysis

The statistical results of KM survival analysis were obtained by a log rank test, and the visual results of KM survival analyses were obtained by survivor and survminer. Univariate COX and multivariate COX models were used to verify whether MMP9 can be used as an independent predictor of the prognosis of SCLC patients treated with cisplatin. The different analysis of continuity variables between MMP9-High and MMP9-Low groups was completed by the Mann Whitney U test. In this study, all the analysis is based on R software. The p value is bilateral, and p < 0.05 is regarded as statistically significant.

Results

MMP9 is an Independent Predictor of the Prognosis of SCLC Treated With Cisplatin

In order to explore the influence of MMP9 on the prognosis of patients with SCLC treated with cisplatin, we used the univariate COX regression model and multivariate COX regression model to evaluate the SCLC cohort, including SCLC (George et al.) and SCLC (Jiang et al.). In SCLC (George et al.), we found that only MMP9-High was related to an improved clinical prognosis of patients with SCLC, while common clinical factors were not relevant in the prognosis of patients (Figure 2A; HR = 0.425; p = 0.0085). Multivariate results showed that MMP-High could be used as an independent predictor of prognosis of SCLC after cisplatin treatment (Figure 2A; HR = 0.216; p = 0.00153). Then, univariate COX and multivariate COX regression models were also applied to SCLC (Jiang et al.), and the results showed that MMP9-High was not only related to significantly prolonged clinical prognosis time (Figure 2B; HR = 0.365; p = 0.0219) but it can also be used as an independent predictor (Figure 2B; HR = 0.352; p = 0.0199). The KM curve shows that MMP9-High is related to significantly prolonged OS in patients with SCLC (Figure 2C: log-rank p = 0.007; HR = 0.44; Figure 2D: log-rank p = 0.017; HR = 0.39). We used a Sankey diagram to visualize the clinical features of two SCLC cohorts one by one (Figures 2E,F).

FIGURE 2
www.frontiersin.org

FIGURE 2. The prognostic value of MMP9. The univariate and multivariate COX regression model of the SCLC (George et al.) (A) and the SCLC (Jiang et al.) (B). Kaplan-Meier curves comparing the progress free survival (PFS) of patients with MMP9-High and patients with MMP9-Low in the SCLC (George et al.) (C) and the SCLC (Jiang et al.) (D). A Sankey diagram visualizing the clinical characteristics between MMP9-High and MMP9-Low patients in the SCLC (George et al.) (E) and the SCLC (Jiang et al.) (F).

MMP9 is Related to Cisplatin Sensitivity and Immunogenicity

We employed the pRRophetic algorithm to predict the cisplatin sensitivity of each SCLC patient to obtain an IC50 value. We found that MMP9-High had significantly lower IC50 values of cisplatin than MMP9-Low SCLC (Figure 3A; p < 0.05). Similarly, the expression of MMP9 was negatively correlated with the IC50 value of cisplatin (Figure 3B, p = 0.001; R = −0.38; method: spearman). We found that MMP9-High had significantly higher TMB than MMP9-Low (Figure 3C; All p < 0.05) and NALs (Figure 3D; all p < 0.05). Figure 3E shows that for SCLC (George et al.), the type and mutation frequency of the driver genes in the top 20 mutation frequencies in the cohort. The results showed that there was no significant difference between the driver genes in the top 20 mutation frequencies of MMP9-High and MMP9-Low. In SCLC (Jiang et al.), the mutation frequency of MMP9-High is significantly higher than that of MMP9-Low in PTPRB and NTRK (Figure 3F). Supplementary Figures S1A,B shows the mutual exclusion and co-occurrence of the top 20 driving mutations in MMP9-High and MMP9-Low, respectively.

FIGURE 3
www.frontiersin.org

FIGURE 3. The association between MMP9, immunogenicity, and cisplatin. (A) Comparison of IC50 values of cisplatin between MMP9-High and MMP9-Low tumors. (B) The association between the IC50 value of cisplatin and the expression of MMP9. (C) Comparison of the tumor mutation burden between MMP9-High and MMP9-Low tumors. (D) Comparison of neoantigen loads (NALs) between MMP9-High and MMP9-Low tumors. The top 20 mutated driver genes in the SCLC (George et al.) (E) and SCLC (Jiang et al.) (F). (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; Mann-Whitney U test).

MMP9 is Related to Activated Immune Cells

Under the MCPcounter algorithm, whether in SCLC (George et al.) or SCLC (Jiang et al.), we found that the TIME when MMP9-High was higher than MMP9-Low was significantly infiltrated with t-cells, cytoxic lymphocytes, B-cells, NK cells, and dense cells (DCS) (Figures 4A,B; all p < 0.05). Under the quanTIseq algorithm, we found that MMP9-High has significantly increased B cells and m1 macrophages when compared with MMP9-Low (Figure 4C: George et al.; Figure 4D: Jiang et al.). With the xCell algorithm, we found that, compared with MMP9-Low, MMP9-High has significant DCs and NKT (Figure 4E: George et al.; Figure 4F: Jiang et al.).

FIGURE 4
www.frontiersin.org

FIGURE 4. The association between MMP9 and immune cells. Comparison of immune cells estimated according to MCPcounter between MMP9-High and MMP9-Low tumors in the SCLC (George et al.) (A) and the SCLC (Jiang et al.) (B). Comparison of immune cells estimated by quanTIseq between MMP9-High and MMP9-Low tumors in the SCLC (George et al.) (C) and the SCLC (Jiang et al.) (D). Comparison of immune cells estimated by xCell between MMP9-High and MMP9-Low tumors in the SCLC (George et al.) (E) and the SCLC (Jiang et al.) (F). (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; Mann-Whitney U test).

MMP9 is Related to the Up-Regulation of the Immune-Related Signaling Pathway

In order to further explore the difference of signal pathway activity between MMP9-High and MMP9-Low, we used GSEA and ssGSEA to evaluate and calculate the signal pathway activity. Figures 5A,B shows that in SCLC (George et al.) and SCLC (Jiang et al.), MMP9-High has significantly increased activity of immune activation related pathways, such as cytokine binding and immune response mediated by B cells and NK cells when compared with MMP9-Low. In addition, the results of ssGSEA analysis showed that the activity of MMP9-High in B-cell, T-cell and NK cell activation, cytokine secretion, chemokine binding, and antigen presentation was significantly higher than that of MMP9-Low (Figures 5C,D).

FIGURE 5
www.frontiersin.org

FIGURE 5. The association between MMP9 and immune-related signaling. (A) The results of GSEA in the SCLC (George et al.) relating to immune signaling. (B) The results of GSEA in the SCLC (Jiang et al.) relating to immune signaling. The GSEA of hallmark gene sets was downloaded from the MSigDB, and each run was performed with 1000 permutations. Differences in pathway activities scored per cell by MMP9-High and MMP9-Low tumors in the SCLC (George et al.) (C) and the SCLC (Jiang et al.) (D). Shown are t values from a linear model.

MMP9 is Related to the Prognosis of Immunotherapy

In order to explore the role of MMP9 in the prognosis of patients receiving immune checkpoint inhibitors, we used CAMOIP as a web tool to verify the relationship between MMP9 and the prognosis of immunotherapy. We found that in the NSCLC (Kim et al.) cohort, the PFS time of MMP9-High was significantly longer than that of MMP9-Low (Figure 6A; log-rank p = 0.026; HR = 0.4). In another NSCLC (Hwang et al.) cohort, the PFS time of MMP9-High was significantly longer than that of MMP9-Low (Figure 6B; log-rank p = 0.029; HR = 0.37). These results suggest that MMP9 may be an important biomarker for the prognosis of immunotherapy.

FIGURE 6
www.frontiersin.org

FIGURE 6. The prognostic value of the MMP9 in immunotherapy. Kaplan-Meier curves comparing the progress free survival (PFS) of patients with MMP9-High and patients with MMP9-Low in the NSCLC (Kim et al.) (A) and the NSCLC (Hwang et al.) (B).

Small Cell Lung Cancer Model Proves that MMP9 is Related to Cisplatin Sensitivity

In this study, the IC50 values of two pairs of chemotherapy-sensitive and drug-resistant small cell lung cancer cells to the first-line chemotherapy drug (cisplatin) were detected using the CCK8 method. The results showed that the IC50 of chemotherapy-sensitive H69/H446 small cell lung cancer cells to cisplatin was significantly lower than that of the corresponding chemotherapy-resistant cells (H69AR/H446DDP) (Figures 7A,B). To detect the effect of MMP9 on the chemotherapy resistance of small cell lung cancer cells, this study used the CCK8 method to detect the IC50 value of chemotherapy-resistant small cell lung cancer cells with upregulated expression of MMP9 to first-line chemotherapy drugs (cisplatin and etoposide). The results showed that upregulating the expression level of MMP9 (H69-MMP9, H446-MMP9) in chemotherapy-resistant small cell lung cancer cells significantly reduced the IC50 value of SCLC cells to chemotherapy drugs (Figures 7C,D), and the subsequent chemotherapy resistance of cells decreased.

FIGURE 7
www.frontiersin.org

FIGURE 7. (A) The IC50 of H69/H69AR cells to cisplatin. (B) The IC50 of H446/H446DDP to cisplatin. (C) The IC50 of H69AR-NC/H69AR-MMP9 to cisplatin. (D) The IC50 of H446DDP-NC/H446DDP-MMP9 to cisplatin. (E) The expression level of MMP9 mRNA in H69, H69AR, H446 and H446DDP. (F) The expression level of CDYL protein in H69, H69AR, H446 and H446DDP. (G) The expression level of MMP9 mRNA in H69AR-NC, H69AR-MMP9, H446DDP-NC and H446DDP-MMP9. (H) The expression level of CDYL protein in H69AR-NC, H69AR-MMP9, H446DDP-NC and H446DDP-MMP9 (***p < 0.001).

Firstly, the basic expression levels of MMP9 mRNA and MMP9 protein in small cell lung cancer cells were detected by real-time fluorescence quantitative PCR and immunohistochemistry. The results of the real-time quantitative PCR showed that the expression level of MMP9 mRNA in chemotherapy-resistant H69AR cells and H446DDP cells was significantly lower than that in chemotherapy-sensitive H69 cells and H446 cells (Figure 7E) (p < 0.001). Immunohistochemically, the results also showed that the expression level of CDYL protein in chemotherapy-resistant H69AR cells and H446DDP cells was lower than that of their parents’ chemotherapy-sensitive H69 cells and H446 cells (Figure 7F). Based on the basic expression level of MMP9 in the above two pairs of chemotherapy-sensitive and chemotherapy-resistant small cell lung cancer cells, this study further uses lentivirus-mediated LV5-MMP9 to up-regulate the expression level of MMP9 in chemotherapy-resistant H69AR and H446DDP small cell lung cancer cells. The verification results of quantitative PCR, and immunohistochemistry showed that we successfully constructed the small cell lung cancer cell model with an up-regulated expression of CDYL (H69AR-MMP9, H446DDP-MMP9). (Figures 7G,H). (p < 0.001).

Discussion

The results of the univariate COX regression, multivariate COX regression, and KM analysis demonstrated that MMP9-High may be an independent predictor of improved prognosis in SCLC patients after receiving cisplatin. In addition, the expression of MMP9 was negatively correlated with the IC50 value of cisplatin. Based on the observations pertaining to MMP9, we analyzed the TIME of patients with SCLC. Compared with MMP9-Low, MMP9-High displayed significantly increased activated immune cells and an amplified active immune activation pathway. In addition, according to the immunotherapy cohort, MMP9 may represent a suitable novel biomarker for screening patients undergoing immunotherapy.

Remodeling the immunogenicity of tumor cells may be one of the reasons underlying the improved prognosis of patients with MMP9-High SCLC after cisplatin treatment. After cisplatin induces the death of tumor cells, it will release immunogenic substances originally located in tumor cells, thus activating the APC-mediated antigen presentation process, resulting in an anti-tumor immune response. Evidence from a study investigating the anti-tumor mechanism of cisplatin employed protein omics based on mass spectrometry in order to detect the content of protein in the supernatant of tumor cell culture before and after cisplatin administration. A total of 2,239 varieties of protein were identified, of which 526 types were up-regulated more than 3 times after cisplatin treatment, including tumor-driving genes such as NRAS, heat shock proteins, metabolic enzymes, and other proteins. Furthermore, APC stimulated by antigenic substances in these supernatants can significantly enhance the proliferation and functional level of CD8+T cells, suggesting that antigen release induced by chemotherapeutic drugs can significantly activate anti-tumor immune response through the antigen presentation system (Beyranvand Nejad et al., 2016; Lin et al., 2020).

In this study, we found that MMP9-High has significantly higher TMB and NALs than MMP9-Low. Up-regulation of immunocompetent ligand on the surface of tumor cells may contribute to the possible mechanisms of improved prognosis in patients with MMP9-High SCLC after cisplatin treatment. In addition to the above-mentioned release of immunogenic substances related to cell death, chemotherapeutic drugs can also affect the protein expression of tumor cells. Immunogenic tumor cells remodeling tumor cells can down-regulate the expression of MHCI to avoid the killing effect of CTL cells. It has been found that cisplatin can up-regulate MHCI on the surface of head and neck cancer cells and enhance the presentation of tumor antigens, thus promoting the recognition of tumor cells by CTL and the activation of CD8+T cells (Gameiro et al., 2012; Tran et al., 2017; Lin et al., 2021b). In addition, cisplatin can also lead to the up-regulation of MHCI expression in ovarian cancer cells (Grabosch et al., 2019). After subcutaneous inoculation of tumor cells pretreated with cisplatin in vitro, it was found that, when compared with the control group, the tumor cells treated with cisplatin were not likely to form tumors, and the expression of MHCI in grown tumors was higher (Nio et al., 2000). Similar results were obtained with regard to colon cancer cells (Ohtsukasa et al., 2003). We can infer that cisplatin can up-regulate the expression of MHCI in vitro and in vivo. In addition to tumor cells, chemotherapeutics can up-regulate the expression level of MHCI in antigen presenting cells (APCs) (Jackaman et al., 2012). In addition, studies reveal that cisplatin can improve the antigen presenting ability of APC such as DCs(Shurin et al., 2009; Zitvogel et al., 2013; Lin et al., 2021c). In this study, we found that MMP9-High has significantly infiltrated DCs, higher MHC, and increased antigen presentation activities compared to that of MMP9-Low.

Inflammatory TIME may be one of the reasons for the improved prognosis of patients with MMP9-High SCLC after cisplatin treatment. Chemotherapy drugs can enhance the sensitivity of tumor cells to immune killing. For example, cisplatin can enhance the sensitivity of tumor cells to the CTL specific killing effect (Ramakrishnan et al., 2012). In addition, active immune effector cells, such as NKs, serve to mediate cytotoxicity (Lichtenstein and Pende, 1986). Furthermore, cisplatin also causes effector cells to produce more cytokines that regulate and promote various immune responses (Kepp et al., 2013; Huang et al., 2021). Moreover, chemotherapy can eliminate immunosuppressive cells such as MDSC and Tregs by inducing apoptosis and through other mechanisms, so that immunotherapy can achieve the maximum efficacy. For example, the use of cisplatin before the injection of DNA vaccine encoding CRT can reduce the level of MDSC in tumor-bearing mice (Tseng et al., 2008; Li et al., 2020). Cisplatin is beneficial in inducing the formation of self-reactive T-cells and anti-tumor immune responses (Tseng et al., 2008). The CTLs and NKs with significant infiltration in MMP9-High, respectively, verified the above results.

However, this research has some limitations. First, the cohorts of SCLC are very limited, and this study only includes two of the SCLC cohorts recorded in the current public data. Secondly, the SCLC cohort lacks data on cisplatin drug sensitivity, and the data on cisplatin drug sensitivity in this study is based on the pRRophetic algorithm.

Conclusion

Based on the results obtained in this study, we identified that MMP9-High may be a potential new biomarker that facilitates the improved prognosis of SCLC patients after cisplatin treatment. In addition, the TIME of MMP9-High SCLC is primarily characterized by a large number of infiltrated activated immune cells and activated immune-related pathways.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

Conceptualization, RL, LW; Formal analysis, ZQ; Visualization, ZQ; Writing–original draft, XH, XW, QH; Writing–review and editing, ZQ, XH, XW, QH, QL, RL, and LW. All authors read and approved the final manuscript.

Funding

This work was supported by the Science and Technology Project of Education Department of Jiangxi Province (GJJ211513, 180812,180805), the Science and Technology Project of Health and Family Planning Commission of Jiangxi Province (Number 20195363, 20195361), the Haosen Oncology Research Fund (Y-HS2019/2-015), and the Science and Technology Project of Ganzhou City (GZ2021SF002).

Conflict of Interest

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.868203/full#supplementary-material

References

Aran, D., Hu, Z., and Butte, A. J. (2017). xCell: Digitally Portraying the Tissue Cellular Heterogeneity Landscape. Genome Biol. 18, 220. doi:10.1186/s13059-017-1349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Becht, E., Giraldo, N. A., Lacroix, L., Buttard, B., Elarouci, N., Petitprez, F., et al. (2016). Estimating the Population Abundance of Tissue-Infiltrating Immune and Stromal Cell Populations Using Gene Expression. Genome Biol. 17, 218. doi:10.1186/s13059-016-1070-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Beyranvand Nejad, E., van der Sluis, T. C., van Duikeren, S., Yagita, H., Janssen, G. M., van Veelen, P. A., et al. (2016). Tumor Eradication by Cisplatin Is Sustained by CD80/86-Mediated Costimulation of CD8+ T Cells. Cancer Res. 76, 6017–6029. doi:10.1158/0008-5472.CAN-16-0881

PubMed Abstract | CrossRef Full Text | Google Scholar

Braicu, E. I., Gasimli, K., Richter, R., Nassir, M., Kümmel, S., Blohmer, J. U., et al. (2014). Role of Serum VEGFA, TIMP2, MMP2 and MMP9 in Monitoring Response to Adjuvant Radiochemotherapy in Patients with Primary Cervical Cancer-Rresults of a Companion Protocol of the Randomized NOGGO-AGO Phase III Clinical Trial. Anticancer Res. 34, 385–391.

PubMed Abstract | Google Scholar

Carter, B. W., Glisson, B. S., Truong, M. T., and Erasmus, J. J. (2014). Small Cell Lung Carcinoma: Staging, Imaging, and Treatment Considerations. RadioGraphics 34, 1707–1721. doi:10.1148/rg.346140178

PubMed Abstract | CrossRef Full Text | Google Scholar

Clough, E., and Barrett, T. (2016). The Gene Expression Omnibus Database. Methods Mol. Biol. 1418, 93–110. doi:10.1007/978-1-4939-3578-9_5

PubMed Abstract | CrossRef Full Text | Google Scholar

de Hoyos, A., and DeCamp, M. M. (2014). Surgery for Small Cell Lung Cancer. Thorac. Surg. Clin. 24, 399–409. doi:10.1016/j.thorsurg.2014.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Caro, G., Cortese, N., Castino, G. F., Grizzi, F., Gavazzi, F., Ridolfi, C., et al. (2016). Dual Prognostic Significance of Tumour-Associated Macrophages in Human Pancreatic Adenocarcinoma Treated or Untreated with Chemotherapy. Gut 65, 1710–1720. doi:10.1136/gutjnl-2015-309193

PubMed Abstract | CrossRef Full Text | Google Scholar

Eriksson, E., Wenthe, J., Irenaeus, S., Loskog, A., and Ullenhag, G. (2016). Gemcitabine Reduces MDSCs, Tregs and TGFβ-1 while Restoring the Teff/treg Ratio in Patients with Pancreatic Cancer. J. Transl Med. 14, 282. doi:10.1186/s12967-016-1037-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Galluzzi, L., Vitale, I., Michels, J., Brenner, C., Szabadkai, G., Harel-Bellan, A., et al. (2014). Systems Biology of Cisplatin Resistance: Past, Present and Future. Cell Death Dis 5, e1257. doi:10.1038/cddis.2013.428

PubMed Abstract | CrossRef Full Text | Google Scholar

Gameiro, S. R., Caballero, J. A., and Hodge, J. W. (2012). Defining the Molecular Signature of Chemotherapy-Mediated Lung Tumor Phenotype Modulation and Increased Susceptibility to T-Cell Killing. Cancer Biother. Radiopharm. 27, 23–35. doi:10.1089/cbr.2012.1203

PubMed Abstract | CrossRef Full Text | Google Scholar

Geeleher, P., Cox, N., and Huang, R. S. (2014). pRRophetic: an R Package for Prediction of Clinical Chemotherapeutic Response from Tumor Gene Expression Levels. PLoS One 9, e107468. doi:10.1371/journal.pone.0107468

PubMed Abstract | CrossRef Full Text | Google Scholar

George, J., Lim, J. S., Jang, S. J., Cun, Y., Ozretić, L., Kong, G., et al. (2015). Comprehensive Genomic Profiles of Small Cell Lung Cancer. Nature 524, 47–53. doi:10.1038/nature14664

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabosch, S., Bulatovic, M., Zeng, F., Ma, T., Zhang, L., Ross, M., et al. (2019). Cisplatin-induced Immune Modulation in Ovarian Cancer Mouse Models with Distinct Inflammation Profiles. Oncogene 38, 2380–2393. doi:10.1038/s41388-018-0581-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hänzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: Gene Set Variation Analysis for Microarray and RNA-Seq Data. BMC Bioinformatics 14, 7. doi:10.1186/1471-2105-14-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Horn, L., Mansfield, A. S., Szczęsna, A., Havel, L., Krzakowski, M., Hochmair, M. J., et al. (2018). First-Line Atezolizumab Plus Chemotherapy in Extensive-Stage Small-Cell Lung Cancer. N. Engl. J. Med. 379, 2220–2229. doi:10.1056/NEJMoa1809064

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, W., Lin, A., Luo, P., Liu, Y., Xu, W., Zhu, W., et al. (2021). EPHA5 Mutation Predicts the Durable Clinical Benefit of Immune Checkpoint Inhibitors in Patients with Lung Adenocarcinoma. Cancer Gene Ther. 28, 864–874. doi:10.1038/s41417-020-0207-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Inoue, M., Koga, F., Yoshida, S., Tamura, T., Fujii, Y., Ito, E., et al. (2014). Significance of ERBB2 Overexpression in Therapeutic Resistance and Cancer-specific Survival in Muscle-Invasive Bladder Cancer Patients Treated with Chemoradiation-Based Selective Bladder-Sparing Approach. Int. J. Radiat. Oncol. Biol. Phys. 90, 303–311. doi:10.1016/j.ijrobp.2014.05.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Ireland, L., Santos, A., Ahmed, M. S., Rainer, C., Nielsen, S. R., Quaranta, V., et al. (2016). Chemoresistance in Pancreatic Cancer Is Driven by Stroma-Derived Insulin-like Growth Factors. Cancer Res. 76, 6851–6863. doi:10.1158/0008-5472.CAN-16-1201

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackaman, C., Majewski, D., Fox, S. A., Nowak, A. K., and Nelson, D. J. (2012). Chemotherapy Broadens the Range of Tumor Antigens Seen by Cytotoxic CD8(+) T Cells In Vivo. Cancer Immunol. Immunother. 61, 2343–2356. doi:10.1007/s00262-012-1307-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Huang, J., Higgs, B. W., Hu, Z., Xiao, Z., Yao, X., et al. (2016). Genomic Landscape Survey Identifies SRSF1 as a Key Oncodriver in Small Cell Lung Cancer. Plos Genet. 12, e1005895. doi:10.1371/journal.pgen.1005895

PubMed Abstract | CrossRef Full Text | Google Scholar

Kepp, O., Menger, L., Vacchelli, E., Locher, C., Adjemian, S., Yamazaki, T., et al. (2013). Crosstalk between ER Stress and Immunogenic Cell Death. Cytokine Growth Factor. Rev. 24, 311–318. doi:10.1016/j.cytogfr.2013.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Lin, A., Luo, P., Shen, W., Xiao, D., Gou, L., et al. (2020). DNAH10 Mutation Correlates with Cisplatin Sensitivity and Tumor Mutation burden in Small-Cell Lung Cancer. Aging (Albany NY) 12, 1285–1303. doi:10.18632/aging.102683

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Wang, H., Zhang, Y., Qiao, R., Xia, P., Kong, Z., et al. (2021). COL3A1 and MMP9 Serve as Potential Diagnostic Biomarkers of Osteoarthritis and Are Associated with Immune Cell Infiltration. Front. Genet. 12, 721258. doi:10.3389/fgene.2021.721258

PubMed Abstract | CrossRef Full Text | Google Scholar

Liberzon, A., Subramanian, A., Pinchback, R., Thorvaldsdóttir, H., Tamayo, P., and Mesirov, J. P. (2011). Molecular Signatures Database (MSigDB) 3.0. Bioinformatics 27, 1739–1740. doi:10.1093/bioinformatics/btr260

PubMed Abstract | CrossRef Full Text | Google Scholar

Lichtenstein, A. K., and Pende, D. (1986). Enhancement of Natural Killer Cytotoxicity by Cis-Diamminedichloroplatinum (II) In Vivo and In Vitro. Cancer Res. 46, 639–644.

PubMed Abstract | Google Scholar

Lin, A., Wei, T., Meng, H., Luo, P., and Zhang, J. (2019). Role of the Dynamic Tumor Microenvironment in Controversies Regarding Immune Checkpoint Inhibitors for the Treatment of Non-small Cell Lung Cancer (NSCLC) with EGFR Mutations. Mol. Cancer 18, 139. doi:10.1186/s12943-019-1062-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, A., Zhang, J., and Luo, P. (2020). Crosstalk between the MSI Status and Tumor Microenvironment in Colorectal Cancer. Front. Immunol. 11, 2039. doi:10.3389/fimmu.2020.02039

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, A., Wei, T., Liang, J., Qi, C., Li, M., Luo, P., et al. (2021). CAMOIP: A Web Server for Comprehensive Analysis on Multi-Omics of Immunotherapy in Pan-Cancer. China: bioRxiv.

Google Scholar

Lin, A., Qiu, Z., Zhang, J., and Luo, P. (2021). Effect of NCOR1 Mutations on Immune Microenvironment and Efficacy of Immune Checkpoint Inhibitors in Patient with Bladder Cancer. Front. Immunol. 12, 630773. doi:10.3389/fimmu.2021.630773

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, A., Xu, W., Luo, P., and Zhang, J. (2021). Mutations Status of Chemokine Signaling Pathway Predict Prognosis of Immune Checkpoint Inhibitors in Colon Adenocarcinoma. Front. Pharmacol. 12, 721181. doi:10.3389/fphar.2021.721181

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, P., Lin, A., Li, K., Wei, T., and Zhang, J. (2019). DDR Pathway Alteration, Tumor Mutation Burden, and Cisplatin Sensitivity in Small Cell Lung Cancer: Difference Detected by Whole Exome and Targeted Gene Sequencing. J. Thorac. Oncol. 14, e276–e279. doi:10.1016/j.jtho.2019.08.2509

PubMed Abstract | CrossRef Full Text | Google Scholar

Makhoul, I., Griffin, R. J., Siegel, E., Lee, J., Dhakal, I., Raj, V., et al. (2016). High-circulating Tie2 Is Associated with Pathologic Complete Response to Chemotherapy and Antiangiogenic Therapy in Breast Cancer. Am. J. Clin. Oncol. 39, 248–254. doi:10.1097/COC.0000000000000046

PubMed Abstract | CrossRef Full Text | Google Scholar

Mayakonda, A., Lin, D. C., Assenov, Y., Plass, C., and Koeffler, H. P. (2018). Maftools: Efficient and Comprehensive Analysis of Somatic Variants in Cancer. Genome Res. 28, 1747–1756. doi:10.1101/gr.239244.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchem, J. B., Brennan, D. J., Knolhoff, B. L., Belt, B. A., Zhu, Y., Sanford, D. E., et al. (2013). Targeting Tumor-Infiltrating Macrophages Decreases Tumor-Initiating Cells, Relieves Immunosuppression, and Improves Chemotherapeutic Responses. Cancer Res. 73, 1128–1141. doi:10.1158/0008-5472.CAN-12-2731

PubMed Abstract | CrossRef Full Text | Google Scholar

Mondal, S., Adhikari, N., Banerjee, S., Amin, S. A., and Jha, T. (2020). Matrix Metalloproteinase-9 (MMP-9) and its Inhibitors in Cancer: A Minireview. Eur. J. Med. Chem. 194, 112260. doi:10.1016/j.ejmech.2020.112260

PubMed Abstract | CrossRef Full Text | Google Scholar

Nio, Y., Hirahara, N., Minari, Y., Iguchi, C., Yamasawa, K., Toga, T., et al. (2000). Induction of Tumor-specific Antitumor Immunity after Chemotherapy with Cisplatin in Mice Bearing MOPC-104E Plasmacytoma by Modulation of MHC Expression on Tumor Surface. Anticancer Res. 20, 3293–3299.

PubMed Abstract | Google Scholar

Ohtsukasa, S., Okabe, S., Yamashita, H., Iwai, T., and Sugihara, K. (2003). Increased Expression of CEA and MHC Class I in Colorectal Cancer Cell Lines Exposed to Chemotherapy Drugs. J. Cancer Res. Clin. Oncol. 129, 719–726. doi:10.1007/s00432-003-0492-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Plattner, C., Finotello, F., and Rieder, D. (2020). Deconvoluting Tumor-Infiltrating Immune Cells From RNA-seq Data Using quanTIseq. Meth. Enzymol. 636, 261–285. doi:10.1016/bs.mie.2019.05.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiao, X., Gu, Y., Yu, J., Wang, J., Liu, X., Gu, M., et al. (2020). The Combination of CD147 and MMP-9 Serum Levels Is Identified as Novel Chemotherapy Response Markers of Advanced Non-small-cell Lung Cancer. Dis. Markers 2020, 8085053. doi:10.1155/2020/8085053

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Z., Lin, A., Li, K., Lin, W., Wang, Q., Wei, T., et al. (2019). A Novel Mutation Panel for Predicting Etoposide Resistance in Small-Cell Lung Cancer. Drug Des. Devel Ther. 13, 2021–2041. doi:10.2147/DDDT.S205633

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramakrishnan, R., Huang, C., Cho, H. I., Lloyd, M., Johnson, J., Ren, X., et al. (2012). Autophagy Induced by Conventional Chemotherapy Mediates Tumor Cell Sensitivity to Immunotherapy. Cancer Res. 72, 5483–5493. doi:10.1158/0008-5472.CAN-12-2236

PubMed Abstract | CrossRef Full Text | Google Scholar

Rauvala, M., Turpeenniemi-Hujanen, T., and Puistola, U. (2006). The Value of Sequential Serum Measurements of Gelatinases and Tissue Inhibitors during Chemotherapy in Ovarian Cancer. Anticancer Res. 26, 4779–4784.

PubMed Abstract | Google Scholar

Reimand, J., Isserlin, R., Voisin, V., Kucera, M., Tannus-Lopes, C., Rostamianfar, A., et al. (2019). Pathway Enrichment Analysis and Visualization of Omics Data Using g:Profiler, GSEA, Cytoscape and EnrichmentMap. Nat. Protoc. 14, 482–517. doi:10.1038/s41596-018-0103-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabari, J. K., Lok, B. H., Laird, J. H., Poirier, J. T., and Rudin, C. M. (2017). Unravelling the Biology of SCLC: Implications for Therapy. Nat. Rev. Clin. Oncol. 14, 549–561. doi:10.1038/nrclinonc.2017.71

PubMed Abstract | CrossRef Full Text | Google Scholar

Shurin, G. V., Tourkova, I. L., Kaneno, R., and Shurin, M. R. (2009). Chemotherapeutic Agents in Noncytotoxic Concentrations Increase Antigen Presentation by Dendritic Cells via an IL-12-dependent Mechanism. J. Immunol. 183, 137–144. doi:10.4049/jimmunol.0900734

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Y., Zhou, X., Bai, W., and Ma, X. (2015). FBW7 Increases Drug Sensitivity to Cisplatin in Human Nasopharyngeal Carcinoma by Downregulating the Expression of Multidrug Resistance-Associated Protein. Tumour Biol. 36, 4197–4202. doi:10.1007/s13277-015-3056-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorsson, V., Gibbs, D. L., Brown, S. D., Wolf, D., Bortone, D. S., Ou Yang, T. H., et al. (2018). The Immune Landscape of Cancer. Immunity 48, 812. e14. doi:10.1016/j.immuni.2018.03.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, L., Allen, C. T., Xiao, R., Moore, E., Davis, R., Park, S. J., et al. (2017). Cisplatin Alters Antitumor Immunity and Synergizes with PD-1/pd-L1 Inhibition in Head and Neck Squamous Cell Carcinoma. Cancer Immunol. Res. 5, 1141–1151. doi:10.1158/2326-6066.CIR-17-0235

PubMed Abstract | CrossRef Full Text | Google Scholar

Tseng, C. W., Hung, C. F., Alvarez, R. D., Trimble, C., Huh, W. K., Kim, D., et al. (2008). Pretreatment with Cisplatin Enhances E7-specific CD8+ T-Cell-Mediated Antitumor Immunity Induced by DNA Vaccination. Clin. Cancer Res. 14, 3185–3192. doi:10.1158/1078-0432.CCR-08-0037

PubMed Abstract | CrossRef Full Text | Google Scholar

Xian, G., Zhao, J., Qin, C., Zhang, Z., Lin, Y., and Su, Z. (2017). Simvastatin Attenuates Macrophage-Mediated Gemcitabine Resistance of Pancreatic Ductal Adenocarcinoma by Regulating the TGF-β1/Gfi-1 axis. Cancer Lett. 385, 65–74. doi:10.1016/j.canlet.2016.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, T., Gao, S., Liu, J., Huang, Y., Chen, K., and Zhang, X. (2021). MMP9 and IGFBP1 Regulate Tumor Immune and Drive Tumor Progression in Clear Cell Renal Cell Carcinoma. J. Cancer 12, 2243–2257. doi:10.7150/jca.48664

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., Zhang, J., and Chao, J. (2012). Porphyromonas Gingivalis Promotes Monocyte Migration by Activating MMP-9. J. Periodontal Res. 47, 236–242. doi:10.1111/j.1600-0765.2011.01427.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zitvogel, L., Galluzzi, L., Smyth, M. J., and Kroemer, G. (2013). Mechanism of Action of Conventional and Targeted Anticancer Therapies: Reinstating Immunosurveillance. Immunity 39, 74–88. doi:10.1016/j.immuni.2013.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: small lung cell cancer, cisplatin, MMP9, survival, resistance

Citation: Wu L, Wang X, He X, Li Q, Hua Q, Liu R and Qiu Z (2022) MMP9 Expression Correlates With Cisplatin Resistance in Small Cell Lung Cancer Patients. Front. Pharmacol. 13:868203. doi: 10.3389/fphar.2022.868203

Received: 02 February 2022; Accepted: 24 February 2022;
Published: 01 April 2022.

Edited by:

Clare Y. Slaney, Peter MacCallum Cancer Centre, Australia

Reviewed by:

Yonglin Yi, Southern Medical University, China
Xiaohua Li, Sixth People’s Hospital of Chengdu, China

Copyright © 2022 Wu, Wang, He, Li, Hua, Liu and Qiu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rongrong Liu, LRR1026@163.com; Zhengang Qiu, qiuzhengang@gmu.edu.cn

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.