Skip to main content

ORIGINAL RESEARCH article

Front. Immunol. , 17 March 2025

Sec. Cancer Immunity and Immunotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1560724

This article is part of the Research Topic The crosstalk between emerging cell death and immune microenvironment remodeling in cancer progression and treatment View all 10 articles

Integrated multi-omics landscape of non-small cell lung cancer with distant metastasis

Teng He&#x;Teng He1†Ting Jiang&#x;Ting Jiang2†Xiaoyuan Sun&#x;Xiaoyuan Sun1†Fang Yang&#x;Fang Yang2†Dan ZhangDan Zhang2Shan YaoShan Yao2Jiangrong Liao*Jiangrong Liao2*Xueling Wu*Xueling Wu1*
  • 1Department of Respiratory and Critical Care Medicine, Ren Ji Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
  • 2Department of Respiratory and Critical Care Medicine, Guizhou Aerospace Hospital, Guizhou, China

Background: Distant metastasis is one of the important factors affecting the prognosis of lung cancer patients. Extracellular vesicles (EVs) play an important role in the occurrence, development, and metastasis of cancer. However, it is currently unclear whether EVs in BALF are involved in distant tumor metastasis.

Methods: we collected bronchoalveolar lavage fluid (BALF) from patients with metastatic and non-metastatic non-small cell lung cancer (NSCLC) to isolate exosomes, which were then characterized by nanoparticle tracking analysis (NTA) and transmission electron microscopy (TEM), followed by comprehensive metabolomic and proteomic analysis to ultimately construct a distant metastasis prediction model for non-small cell lung cancer.

Results: Our research has found that the BALF of NSCLC patients is rich in EVs, which have typical morphology and size. There are significant differences in protein expression and metabolite types between patients with distant metastasis and those without distant metastasis. Sphingolipid metabolism pathways may be a key factor influencing distant metastasis in NSCLC. Subsequently, we constructed a predictive model for distant metastasis in NSCLC based on differentially expressed proteins identified by proteomics. This model has been proven to have high predictive value.

Conclusion: The multi-omic analysis generated in this study provided a global overview of the molecular changes, which may provide useful insight into the therapy and prognosis of NSCLC metastasis

Introduction

Regardless of the great advances made in lung cancer treatment such as molecular targeted therapies and immunologic control over the past decade, it remains one of the most deadliest forms of cancer globally, accounting for 24% of all cancer deaths (1). Approximately 57% of lung cancer patients have already developed distant metastasis at the time of initial diagnosis. Furthermore, tumor metastasis accounts for 67%-90% of all cancer mortality, with the figure exceeding 70% in lung cancer (2). Brain metastasis is one of the most common metastatic sites. The incidence rate is about 7.3 cases per 100000 people (3), of which about half are NSCLC patients. Furthermore, the one-year survival rate of NSCLC patients with distant metastasis is only 15% -19% (4). Therefore, the study of mechanisms of lung cancer metastasis will be of great significance.

Exosomes are a subtype of EVs, which are lipid bilayer-enclosed vesicles, spherical or cup-shaped, with diameters ranging from 30 to 150nm (5). Whether under physiological or pathological conditions, almost all cells secrete exosomes, widely distributed in physiological fluids, including blood, saliva, urine, and bronchoalveolar lavage fluid (6). They play a critical role in cell-cell communication and can regulate various biological and pathological processes such as inflammation and immune responses, tumor growth, and metastasis (7).

Lung tumor-derived exosomes has received increasing research attention in metastatic activity integrated with immune system evasion, epithelia-mesenchymal transition and angiogenesis. It has been demonstrated that tumor exosome integrins (ITGs) determine organotropic metastasis with exosome proteomics. Exosome-based ITGα6β4 and ITGα6β1 were associated with lung metastasis, while ITGαvβ5 was linked to liver metastasis. In their observations, exosome integrin isolated from circulating plasma of cancer patients uptake by resident cells activates Src phosphorylation and pro-inflammatory S100 gene expression (8). A recent study with plasma exosomes from lung cancer cases having liver metastasis (LM) or bone metastasis (BM) was subjected to proteomic profiles. SELL and MUC5B could be used as diagnostic markers of BM, while APOH, CD81, and CCT5 could help diagnose LM in local advanced lung cancer patients (9). Hypoxic bone marrow-derived mesenchymal stem cells (BMSC)-derived exosomal miRNAs promote metastasis of lung cancer with a microRNA array (10). Further study demonstrated that plasma exosome-mediated transfer of miR-193a-3p, miR-210-3p and miR-5100 could promote invasion of lung cancer cells by activating STAT3 signaling-induced EMT. Using lncRNA-seq, a novel exosomal long noncoding RNA (lncRNA), lnc-MLETA1 was found to promote tumor metastasis by regulating the miR-186-5p/EGFR and miR-497-5p/IGF1R axes in NSCLC (11).

BALF has been considered to be a sample that best represents the condition of pulmonary diseases, except for lung biopsy and surgical specimen. Exosomes in BALF has been shown to be a key factor involved in the growth and progression of lung cancer in vitro and in vivo (12). Certainly, BALF also has limitations, such as low protein content, the need for invasive procedures to obtain it, and the risk of contamination. However, knowledge of the role of BALF exosomes highlighting for being as biomarkers in diagnosis of cancer metastasis remain acquaint scarcely. In this study, exosomes extracted from BALF of patients with NSCLC were confirmed through NTA. Proteomic and metabolic analyses were subsequently conducted to explore the role of exosomes in NSCLC.

Methods

Ethical approval and clinical specimen collection

The study involving human participants was reviewed and approved by the Renji Hospital, Shanghai Jiaotong University School of Medicine, China (No. LY2024-179-A) and Guizhou Aerospace Hospital, China (No. (2024)1-050), and providing written informed consent in accordance with the regulations of the relevant Hospital Organizational Review Board and the Declaration of Helsinki. NSCLC patients were collected from June 2023 to June 2024 at the respiratory department of the aforementioned hospital. All patients were over 18 years old, confirmed by pathology as NSCLC with or without distant metastasis, and were able to obtain complete and comprehensive clinical data. Exclude patients with infectious diseases such as active tuberculosis, HIV, and viral hepatitis, other types of tumors, and those who refuse or are unable to cooperate with bronchial examinations. All patients received preliminary diagnosis and obtained BALF from the aforementioned hospital.

EV purification and characterization

BALF-derived EVs were isolated with a combined chromatography kit, following the manufacturer’s instructions. Briefly, 1 ml of BALF was concentrated using a 30K molecular weight cut-off (MWCO) spin filter (Amicon, Merck Millipore, MA, USA). Then 100 μL of concentrated BALF were loaded onto columns pre-packed with 500μL of Sepharose CL-6B resins (Cytiva). 200 μL filtrate were collected and added to columns containing a specific mixture of Fractogel EMD SO3- (M) (Millipore Sigma) and Capto Core 700 (Cytiva). After incubation, the columns were centrifuged at 500 x g for 30 seconds to separate and isolate the EVs, which were then prepared for further downstream analysis.

EVs were characterized by NTA and TEM. NTA was carried out using a ZetaView instrument (Particle Metrix, Germany) to assess the particle size distribution and concentration. For NTA measurement, according to the ZetaView manual, routine calibration and instrument performance verification were performed using NanoStandard™ Series 100 nm standards. Before each sample injection, the instrument was flushed with PBS to ensure the absence of interfering particles within the lens. Samples were diluted with PBS to achieve 50-400 particles per frame. Subsequently, three cycles were performed by scanning 11 cell positions each and capturing 30 frames per position.

TEM was utilized to examine the morphological features of the EVs. A 30 μl aliquot of the sample suspension was placed on a carbon-coated copper grid for 10 minutes. Excess liquid was blotted away, followed by staining with uranyl acetate solution. After removing the residual staining solution, the grid was air-dried on filter paper for 3 hours. Subsequently, the samples were then observed under the microscope (JEOL, Japan).

Proteomics sequencing

All samples were analyzed for proteomics by EVbio, China. Extract proteins from the above samples using protease inhibitors and ultrasonic lysis. Subsequently, reduction alkylation, SP3 enzyme digestion, and desalination treatment were performed, and the Vanquish Neo/Astral mass spectrometer was used to generate mass spectrometry detection raw data in data-independent acquisition (DIA) mode. The raw MS files were analyzed and searched against target protein database based on the species of the samples using DIA-NN. The parameters were set as follows: the protein modifications were carbamidomethylation (C) (fixed), oxidation (M) (variable), Acetyl (Protein N-term) (variable), the enzyme specificity was set to trypsin, the maximum missed cleavages were set to 2, Precursor Qvalue Cutoff was set to 20 ppm, and Protein Qvalue Cutoff was set to 20ppm.

Metabolomics sequencing

All samples were analyzed for metabolomics by EVbio, China. Analyze exosome samples using the Thermo Fisher Scientific U3000 liquid chromatograph, and perform substance annotation identification with Progenesis QI software. After peak alignment, baseline correction, and noise reduction, proceed with feature alignment as well as quantitative and qualitative analysis. Principles for quantitative screening of metabolites: RSD value is less than 30, more than 50% of the QC samples have expression (including 50%), and more than 30% of the experimental samples have expression (including 30%). Principles for qualitative screening of metabolites: eliminate metabolites with a fragmentation score of less than 30, and deduplicate the annotation results from the HMDB database.

Construction of a biomarker prediction model

We have constructed a predictive model for distant metastasis of NSCLC using selected proteins. This model is based on the Least Absolute Shrinkage and Selection Operator (LASSO) analysis, modeled with neural networks, and finally evaluated using the receiver operating characteristic (ROC) curve with an area under the curve (AUC) to calculate the diagnostic performance of the model.

Statistical analysis

Mann-Whitney U tests were used for comparisons of continuous variables, whereas Chi-square tests were used for categorical variables. Correlation analysis was assessed by Spearman correlation. P value of < 0.05 (two-tailed) was considered statistically significant.

Results

Participants and exosome characteristics

We collected BALF samples from 30 NSCLC patients, including 14 patients with distant metastasis (Group A) and 16 patients with non-distant metastasis (Group B). The clinical characteristics of patient are shown in Table 1. The workflow of proteomics and omics experiments for BALF samples is shown in Figure 1A. The average age of Group A is 67.6 (SD, 9.2) with males accounting for 78.6%. The average age of Group B is 68.8 (SD, 9.5) with males accounting for 81.3%. There was no significant difference in laboratory tests between Group A and Group B (Table 2). We conducted an analysis of the collected patient BALF using TEM and NTA identified exosomes with typical morphological and size characteristics (Figure 1B), with an average diameter of 126.7 nm in Group A and 113.7 nm in Group B (Figure 1C).

Table 1
www.frontiersin.org

Table 1. Basic characteristics of patients.

Figure 1
www.frontiersin.org

Figure 1. Overview of research design and characteristics of EVs. (A) The workflow of proteomics and omics experiments for BALF samples. (B) Typical TEM images of EVs in BALF (white arrows). Scale bar: 200 nm. (C) Using NAT to obtain typical images of the particle size distribution and concentration of EVs in BALF.

Table 2
www.frontiersin.org

Table 2. Hematological examination characteristics of patients.

Proteomic changes and biological pathways of EVs in BALF of NSCLC patients

After filtering the low abundant proteins, 4255 high-quality proteins were collected for the data analysis. The median protein number for 30 samples was 3361 (Figures 2A, B). Our preliminary analysis employed PCA (Figure 2C) and PLS-DA (Figure 2D) method to explore the clustering patterns among the groups, which demonstrated a distinct separation between Group A and Group B. A comprehensive differential expression analysis revealed that 209 proteins exhibited unique changes in NSCLC patients with distant metastasis compared to those without, with 109 proteins showing continuous downregulation and 100 proteins showing upregulation (Figure 2E). These findings demonstrated the characteristics of serum protein alterations correlating with the metastasis of the disease. Among the upregulated proteins, Myotubularin-related protein 3 (MTMR3) (Supplementary Figure S1A) and Interleukin-1 Receptor Associated Kinase 1(IRAK4) (Supplementary Figure S1B) are the two most significantly altered proteins. On the contrary, transmembrane conductance regulator (CFTR) is one of the most significantly downregulated proteins (Supplementary Figure S1C). Subsequently, KEGG analysis showed significant changes in the following pathways, Human immunodeficiency virus 1 infection, Autophagy, Rap1 signaling pathway, Estrogen signaling pathway, Sphingolipid signaling pathway, Phospholipase D signaling pathway, Neurotrophin signaling pathway (Figure 2F).

Figure 2
www.frontiersin.org

Figure 2. Proteomic changes and biological pathways of EVs in BALF of NSCLC patients. (A) Quantity of EVs proteins in BALF of each patient. (B) Differences in protein abundance between two groups of EVs. (C) PCA analysis of EVs proteins in distant metastasis group and non-distant metastasis group of lung cancer. (D) PLS-DA analysis of EVs proteins in group A and group B (E) Volcano plot shows the differential proteins between the group A and group B, each dot represents a protein, with red dots for proteins significantly upregulated in brain metastasis group and blue dots for proteins significantly downregulated in brain metastasis group. Fold change ≥ 2.0 and P < 0.05 is considered statistically significant. (F) KEGG pathway analysis of differential proteins.

Metabolomic changes of EVs in BALF of lung cancer patients

In our metabolic investigation, we cataloged a comprehensive array of 651 metabolites, encompassing amino acids, lipids, and other critical metabolites. The PLS and OPLS plots revealed a pronounced metabolic differentiation in B VS A groups (Figures 3A, B). The most significant changes in metabolites include 2-Keto-3-Deoxy-D-Manooctanoic Acid,3-Methyl-2-thioxoimidazolidin-4-one, L-Phenylalanine, Metyrosine, N-Palmitoyltaurine, PE(P-18:0_18:1(12Z)-2OH(9,10)), PG(i-24:0_PGF1alpha), Isoleucyl-Arginine, N-Acetylcysteinamide, Trimetaphosphoric acid (Figures 3C, D). Next, we further analyzed the signaling pathways involved through KEGG analysis, including Biosynthesis of alkaloids derived from, Aminoacyl-tRNA biosynthesis, Protein digestion and absorption, Biosynthesis of plant hormones, Biosynthesis of various secondary, 2-Oxocarboxylic acid metabolism, Phenylalanine metabolism, Central carbon metabolism in cancer, Cyanoamino acid metabolism, Mineral absorption, Glucosinolate biosynthesis, Biosynthesis of amino acids, Biosynthesis of plant secondary metabolites, Phenylalanine, tyrosine and tryptophan biosynthesis, Biosynthesis of phenylpropanoids, Phenylpropanoid biosynthesis, Biosynthesis of antibiotics, ABC transporters (Figure 3E).

Figure 3
www.frontiersin.org

Figure 3. Metabolomic changes of EVs in BALF of NSCLC patients. (A) PLS analysis of EVs metabolites in group A and group B (B) OPLS analysis of EVs metabolites in group A and group B (C) The heatmap showed the 22 metabolites with significant changes. (D) The lollipop chart shows the 10 down regulated proteins and 1 up regulated protein with the most significant changes in the distant metastasis group of lung cancer. (E) KEGG pathway analysis of differential metabolites.

Biomarker panel for early prediction of distant metastasis in NSCLC patients

We have constructed a biomarker model for the early prediction of lung cancer distant metastasis using the aforementioned proteomics data. 662 candidate variables obtained proteomics were screened by LASSO regression analysis (Figure 4A). Ultimately, a prediction model for lung cancer distant metastasis composed of 15 proteins was obtained (Figure 4B). Subsequently, we analyzed the Area Under the Curve (AUC) for each of the 15 proteins including AGAL, DCIL2, TPD54, CBL, PML, PCP, ASGL1, LEMD2, UBA5, PLCL2, CK052, H1.5, CEAM6, A1AG2 and FAK1 and found that their diagnostic performance ranged from 0.648 to 0.879 (Figures 4C–Q). The protein model composed of 15 types of proteins has shown excellent predictive power, with an AUC of 0.95 (95% CI 0.917-0.976), which proves this point (Figure 4R).

Figure 4
www.frontiersin.org

Figure 4. Biomarker panel for early prediction of distant metastasis in NSCLC patients. (A) Constructing LASSO regression analysis model using differential proteins. (B) Prediction model of biomarkers for distant metastasis of lung cancer containing 15 proteins. (C–Q) Receiver operating characteristic curve (ROC) of AGAL, DCIL2, TPD54, CBL, PML, PCP, ASGL1, LEMD2, UBA5, PLCL2, CK052, H1.5, CEAM6, A1AG2 and FAK1. (R) The overall predictive performance of protein models. 95% confidence interval displayed on the graph.

Discussion

This study is the first to use exosomes from patients’ BALF as the research object, and systematically elucidates the differences between NSCLC patients with metastasis and those without metastasis by combining proteomics and metabolomics. Furthermore, we utilized the aforementioned results to construct a predictive model containing 15 proteins, providing reliable biomarkers for the prediction of distant metastasis in NSCLC patients.

BALF is a common method of injecting saline into the alveoli and then retrieving alveolar cell components and secretions through negative pressure aspiration, which can well reflect the changes within the lung tissue and is therefore referred to as a liquid biopsy (13, 14). TME is shaped and conditioned by cancer cells, and it also serves as a critical target for therapy (15). EVs in BALF have become a new means of intercellular communication (16). Similarly, our research found that BALF of NSCLC patients contains a rich array of exosomes, which was confirmed by NTA.

Subsequently, we performed proteomic analysis on the collected exosome samples. The results indicated that there were abundant differences in protein expression between the distant metastasis group and the non-distant metastasis group. Moreover, KEGG analysis suggested the potential role of the sphingolipid signaling pathway particularly the key protein MTMR3. MTMR3 is a member of the myotubularin protein family (17). Essentially, it is a phosphatidylinositol 3-phosphatase that can hydrolyze phosphatidylinositol-3-phosphate, affecting phosphatidylinositol metabolism and thus participating in the sphingolipid metabolism pathway. The MTMR3 protein has been confirmed to be associated with the proliferation, migration, and autophagy of tumor cells in breast cancer (18, 19). Additionally, the MTMR3 protein plays an important role in phosphoinositide metabolism and is a crucial component of the sphingolipid metabolic pathway (17). Research indicates that the high expression level of MTMR3 may be correlated with a lower overall survival period and a shorter recurrence-free survival period in breast cancer patients (19). The underlying mechanism might be caused by MTMR3 influencing the transition of tumor cells from G1 to S phase. In oral cancer, miR-99a can inhibit the proliferation, migration and invasion of tumor cells by reducing the expression of MTMR3 (20). IRAK-4-NF-kB is an important component of the Toll-like receptor (TLR) signaling pathway (21). Immunohistochemical studies on melanoma have shown that phosphorylated IRAK-4 is widely expressed in diseased tissues (22). Moreover, silencing the gene expression of IRAK-4 by shRNA or the application of small molecule inhibitors can induce the proliferation of tumor cells in patients with T-acute lymphoblastic leukemia, suggesting that the signaling mediated by IRAK-4 is an important factor in disease progression (23). Interestingly, the inflammatory response can promote the expression of various inflammatory factors by activating IRAK-4- NF-kB, thereby participating in sphingolipid metabolism (24).

Furthermore, metabolomics analysis indicates that multiple significantly different metabolites are associated with the sphingolipid signaling pathway. Phosphatidylethanolamine (PE) itself is the final product of sphingolipid degradation. Specifically, sphingosine 1-phosphate (S1P) is decomposed into PE and hexadecenal under the action of sphingosine-1-phosphate lyase (SPL) (25). Additionally, PE serves as a primary supplier of phosphate groups in the synthesis of sphingomyelin (26). N-Palmitoyltaurine is a derivative of taurine, which is an important component of sphingosine in the brain. In addition, L-Phenylalanine is one of the precursors for the synthesis of ceramide, and it can also serve as a substrate to participate in the synthesis of sphingosine, thereby affecting sphingolipid metabolism. Research indicates that a novel probiotic, Christensenella minuta, can regulate cytochrome P450 metabolism, sphingolipid metabolism, PI3K-AKT pathway, and integrin pathway by affecting phenylalanine metabolism, thereby alleviating acetaminophen mediated liver injury (27). All of these evidences support the potential role of sphingolipid metabolism in the distant metastasis of NSCLC.

In this study, based on the results of exosomal proteomics, we constructed a model containing 15 proteins for predicting distant metastasis in NSCLC patients. This model provides an example for early identification of distant metastases in NSCLC. Interestingly, TPD54 has been shown to interact with intracellular nanovesicles and play a role in regulating cell migration and invasion. A previous study revealed that in prostate cancer patients, the mRNA and protein expression levels of TPD54 were significantly upregulated and closely associated with poor prognosis (28). In vitro experiments demonstrated that overexpression of the TPD54 gene in oral squamous cell carcinoma cells significantly increased the size of cell colonies (29). The Casitas B-lineage lymphoma (CBL) family, as a class of ubiquitin ligases comprising c-CBL, CBL-b, and CBL-c, can regulate signal transduction in multiple tyrosine kinase-dependent pathways to mediate tumor growth (30). Studies have shown that Cbl-b is upregulated in exhausted (PD1+Tim3+) CD8+ T cells, and knocking out the Cbl-b gene in CAR-T cells renders them resistant to exhaustion (31). Furthermore, in NSCLC xenograft models, reducing CBL-c expression diminishes tumor cell viability and clonogenicity, mediating the inhibition of tumor growth. Interestingly, silencing CBL-c expression significantly enhances the sensitivity of EGFR-mutant NSCLC cells to tyrosine kinase inhibitors. Therefore, targeting CBL-c represents a promising strategy to suppress tumor growth and improve the therapeutic response of EGFR-mutant NSCLC to tyrosine kinase inhibitor treatment (32). Similarly, by reducing the autophosphorylation kinase activity of focal adhesion kinase 1 (FAK1), the cell proliferation of FC-IBC02, SUM190, and KPL4 can be decreased (33). In summary, these biomarkers are closely related to the proliferation, migration, and invasion of tumor cells. In summary, the expression of these 15 proteins is closely associated with tumor cell proliferation, migration, and invasion. This model achieved an AUC of 0.95 in predicting distant metastasis of NSCLC. Additionally, cancer cells undergo profound alterations in energy metabolism during metastasis. Our study revealed significant changes in the sphingolipid metabolic pathway in NSCLC patients with distant metastasis. Therefore, combining proteomic models with metabolite detection may be one of the approaches to further improve the accuracy of this predictive model. Interestingly, in our study, the non-distant metastasis group included patients with stages I, II, and III, who also have the potential to develop distant metastasis, and there are differences among these stages. Therefore, the heterogeneity in staging between the groups may introduce confounding factors that could influence the results. Consequently, further subgroup analyses are needed to validate the independent effects of different stages versus metastasis itself.

Our research still has limitations. Firstly, further in vivo experiments are needed to confirm the role of the sphingolipid metabolism pathway in distant metastasis of NSCLC. Secondly, the distant metastasis prediction model contains a large number of proteins, which limits its widespread clinical application. Finally, the combined analysis of proteomics and metabolomics, along with errors in data processing, may affect the accuracy of the results.

In conclusion, this comprehensive study highlights the importance of lipid metabolism, in elucidating the pathogenesis of NSCLC metastasis. These findings hold significant implications for potential risk assessment and therapeutic strategies in lung cancer. Furthermore, following further detailed validation, the results of this multi-omics analysis are expected to advance toward broader clinical applications.

Data availability statement

The sequencing data for this article has been uploaded to the iProX database with the Project ID: IPX0011284000, and the specific access link is: https://www.iprox.cn/page/project.html?id=IPX0011284000.

Ethics statement

The studies involving humans were approved by Renji Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Guizhou Aerospace Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

TH: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing, Resources, Software, Supervision, Validation. TJ: Conceptualization, Investigation, Methodology, Writing – review & editing. XS: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. FY: Investigation, Methodology, Project administration, Writing – review & editing. DZ: Investigation, Methodology, Project administration, Writing – review & editing. SY: Investigation, Methodology, Project administration, Writing – review & editing. JL: Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing. XW: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by funding from Shanghai Science and Technology Innovation Action Plan Medical Innovation Research Special Project (23Y31900103), Shanghai Science and Technology Innovation Action Plan Western Medicine Guidance Research Special Project (20Z11901003) and Innovation Project of Renji Hospital Affiliated to Shanghai Jiao Tong University School of Medicine (2019NYLYCP0103).

Acknowledgments

We thank all the participants.

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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/fimmu.2025.1560724/full#supplementary-material

References

1. Yin L, Liu X, Shao X, Feng T, Xu J, Wang Q, et al. The role of exosomes in lung cancer metastasis and clinical applications: an updated review. J Trans Med. (2021) 19:312. doi: 10.1186/s12967-021-02985-1

PubMed Abstract | Crossref Full Text | Google Scholar

2. Sauer S, Reed DR, Ihnat M, Hurst RE, Warshawsky D, Barkan D. Innovative approaches in the battle against cancer recurrence: novel strategies to combat dormant disseminated tumor cells. Front Oncol. (2021) 11:659963. doi: 10.3389/fonc.2021.659963

PubMed Abstract | Crossref Full Text | Google Scholar

3. Singh R, Stoltzfus KC, Chen H, Louie AV, Lehrer EJ, Horn SR, et al. Epidemiology of synchronous brain metastases. Neurooncol Adv. (2020) 2:vdaa041. doi: 10.1093/noajnl/vdaa041

PubMed Abstract | Crossref Full Text | Google Scholar

4. Blandin Knight S, Crosbie PA, Balata H, Chudziak J, Hussell T, Dive C. Progress and prospects of early detection in lung cancer. Open Biol. (2017) 7:170070. doi: 10.1098/rsob.170070

PubMed Abstract | Crossref Full Text | Google Scholar

5. Kalluri R, LeBleu VS. The biology, function, and biomedical applications of exosomes. Sci (New York NY). (2020) 367:eaau6977. doi: 10.1126/science.aau6977

PubMed Abstract | Crossref Full Text | Google Scholar

6. Zhou B, Xu K, Zheng X, Chen T, Wang J, Song Y, et al. Application of exosomes as liquid biopsy in clinical diagnosis. Signal Transduction Targeted Ther. (2020) 5:144. doi: 10.1038/s41392-020-00258-9

PubMed Abstract | Crossref Full Text | Google Scholar

7. Steinbichler TB, Dudás J, Riechelmann H, Skvortsova I-I. The role of exosomes in cancer metastasis. Semin Cancer Biol. (2017) 44:170–81. doi: 10.1016/j.semcancer.2017.02.006

PubMed Abstract | Crossref Full Text | Google Scholar

8. Hoshino A, Costa-Silva B, Shen TL, Rodrigues G, Hashimoto A, Tesic Mark M, et al. Tumour exosome integrins determine organotropic metastasis. Nature. (2015) 527:1476–4687. doi: 10.1038/nature15756

PubMed Abstract | Crossref Full Text | Google Scholar

9. Li S, Qu Y, Liu L, Zhang X, He Y, Wang C, et al. Comparative proteomic profiling of plasma exosomes in lung cancer cases of liver and brain metastasis. Cell Bioscience. (2023) 13:180. doi: 10.1186/s13578-023-01112-5

PubMed Abstract | Crossref Full Text | Google Scholar

10. Zhang X, Sai B, Wang F, Wang L, Wang Y, Zheng L, et al. Hypoxic BMSC-derived exosomal miRNAs promote metastasis of lung cancer cells via STAT3-induced EMT. Molecular Cancer. (2019) 18:40. doi: 10.1186/s12943-019-0959-5

PubMed Abstract | Crossref Full Text | Google Scholar

11. Hsu XR, Wu JE, Wu YY, Hsiao SY, Liang JL, Wu YJ, et al. Exosomal long noncoding RNA MLETA1 promotes tumor progression and metastasis by regulating the miR-186-5p/EGFR and miR-497-5p/IGF1R axes in non-small cell lung cancer. J Experimental Clin Cancer Res. (2023) 42:283. doi: 10.1186/s13046-023-02859-y

PubMed Abstract | Crossref Full Text | Google Scholar

12. Yang Y, Ji P, Wang X, Zhou H, Wu J, Quan W, et al. Bronchoalveolar lavage fluid-derived exosomes: A novel role contributing to lung cancer growth. Front Oncol. (2019) 9:197. doi: 10.3389/fonc.2019.00197

PubMed Abstract | Crossref Full Text | Google Scholar

13. Liu Z, Yan J, Tong L, Liu S, Zhang Y. The role of exosomes from BALF in lung disease. J Cell Physiol. (2022) 237:161–8. doi: 10.1002/jcp.v237.1

Crossref Full Text | Google Scholar

14. Zhang H, Deng D, Li S, Ren J, Huang W, Liu D, et al. Bronchoalveolar lavage fluid assessment facilitates precision medicine for lung cancer. Cancer Biol Med. (2023) 21:230–51. doi: 10.20892/j.issn.2095-3941.2023.0381

PubMed Abstract | Crossref Full Text | Google Scholar

15. Jin Y, Cai S, Zhou Y, Guo D, Zeng Y, Xu W, et al. Targeting SLC7A11/xCT improves radiofrequency ablation efficacy of HCC by dendritic cells mediated anti-tumor immune response. Imeta. (2024) 3:e248. doi: 10.1002/imt2.v3.6

PubMed Abstract | Crossref Full Text | Google Scholar

16. Huang J, Ding M, Lu Y, Xu L, Zhang Y, Han S, et al. MiR-1246b, a novel miRNA molecule of extracellular vesicles in bronchoalveolar lavage fluid, promotes nodule growth through FGF14 in patients with lung cancer. Cell Death Dis. (2023) 14:789. doi: 10.1038/s41419-023-06218-9

PubMed Abstract | Crossref Full Text | Google Scholar

17. Walker DM, Urbé S, Dove SK, Tenza D, Raposo G, Clague MJ. Characterization of MTMR3. an inositol lipid 3-phosphatase with novel substrate specificity. Curr Biol. (2001) 11:1600–5. doi: 10.1016/S0960-9822(01)00501-2

PubMed Abstract | Crossref Full Text | Google Scholar

18. Wang B, Mao JH, Wang BY, Wang LX, Wen HY, Xu LJ, et al. Exosomal miR-1910-3p promotes proliferation, metastasis, and autophagy of breast cancer cells by targeting MTMR3 and activating the NF-κB signaling pathway. Cancer Lett. (2020) 489:87–99. doi: 10.1016/j.canlet.2020.05.038

PubMed Abstract | Crossref Full Text | Google Scholar

19. Wang Z, Zhang M, Shan R, Wang YJ, Chen J, Huang J, et al. MTMR3 is upregulated in patients with breast cancer and regulates proliferation, cell cycle progression and autophagy in breast cancer cells. Oncology Reports. (2019) 42:1915–23. doi: 10.3892/or.2019.7292

PubMed Abstract | Crossref Full Text | Google Scholar

20. Kuo YZ, Tai YH, Lo HI, Chen YL, Cheng HC, Fang WY, et al. MiR-99a exerts anti-metastasis through inhibiting myotubularin-related protein 3 expression in oral cancer. Oral Dis. (2014) 20:e65–75. doi: 10.1111/odi.2014.20.issue-3

PubMed Abstract | Crossref Full Text | Google Scholar

21. Bai YR, Yang WG, Hou XH, Shen DD, Zhang SN, Li Y, et al. The recent advance of Interleukin-1 receptor associated kinase 4 inhibitors for the treatment of inflammation and related diseases. Eur J Med Chem. (2023) 258:115606. doi: 10.1016/j.ejmech.2023.115606

PubMed Abstract | Crossref Full Text | Google Scholar

22. Srivastava R, Geng D, Liu Y, Zheng L, Li Z, Joseph MA, et al. Augmentation of therapeutic responses in melanoma by inhibition of IRAK-1,-4. Cancer Res. (2012) 72:6209–16. doi: 10.1158/0008-5472.CAN-12-0337

PubMed Abstract | Crossref Full Text | Google Scholar

23. Yadav H, Shirumalla RK. Emerging trends in IRAK-4 kinase research. Mol Biol Rep. (2023) 50:7825–37. doi: 10.1007/s11033-023-08438-w

PubMed Abstract | Crossref Full Text | Google Scholar

24. Sukocheva OA, Neganova ME, Aleksandrova Y, Burcher JT, Chugunova E, Fan R, et al. Signaling controversy and future therapeutical perspectives of targeting sphingolipid network in cancer immune editing and resistance to tumor necrosis factor-α immunotherapy. Cell Commun Signal. (2024) 22:251. doi: 10.1186/s12964-024-01626-6

PubMed Abstract | Crossref Full Text | Google Scholar

25. Atreya KB, Saba JD. Neurological consequences of sphingosine phosphate lyase insufficiency. Front Cell Neurosci. (2022) 16:938693. doi: 10.3389/fncel.2022.938693

PubMed Abstract | Crossref Full Text | Google Scholar

26. Hu K, Zhang Q, Chen Y, Yang J, Xia Y, Rao B, et al. Cryo-EM structure of human sphingomyelin synthase and its mechanistic implications for sphingomyelin synthesis. Nat Struct Mol Biol. (2024) 31:884–95. doi: 10.1038/s41594-024-01237-2

PubMed Abstract | Crossref Full Text | Google Scholar

27. Yao T, Fu L, Wu Y, Li L. Christensenella minuta alleviates acetaminophen-induced hepatotoxicity by regulating phenylalanine metabolism. Nutrients. (2024) 16(14):2314. doi: 10.3390/nu16142314

PubMed Abstract | Crossref Full Text | Google Scholar

28. Ren L, Chen J, Zhang X. Increased expression of tumor protein D54 is associated with clinical progression and poor prognosis in patients with prostate cancer. Oncol Lett. (2017) 14:7739–44. doi: 10.3892/ol.2017.7214

PubMed Abstract | Crossref Full Text | Google Scholar

29. Kato K, Mukudai Y, Motohashi H, Ito C, Kamoshida S, Shimane T, et al. Opposite effects of tumor protein D (TPD) 52 and TPD54 on oral squamous cell carcinoma cells. Int J Oncol. (2017) 50:1634–46. doi: 10.3892/ijo.2017.3929

PubMed Abstract | Crossref Full Text | Google Scholar

30. Ren J, Lv L, Tao X, Zhai X, Chen X, Yu H, et al. The role of CBL family ubiquitin ligases in cancer progression and therapeutic strategies. Front Pharmacol. (2024) 26:1432545. doi: 10.3389/fphar.2024.1432545

PubMed Abstract | Crossref Full Text | Google Scholar

31. Kumar J, Kumar R, Kumar Singh A, Tsakem EL, Kathania M, Riese MJ, et al. Deletion of Cbl-b inhibits CD8(+) T-cell exhaustion and promotes CAR T-cell function. J immunotherapy cancer. (2021) 9:e001688. doi: 10.1136/jitc-2020-001688

PubMed Abstract | Crossref Full Text | Google Scholar

32. Hong SY, Kao YR, Lee TC, Wu CW. Upregulation of E3 ubiquitin ligase CBLC enhances EGFR dysregulation and signaling in lung adenocarcinoma. Cancer Res. (2018) 78:4984–96. doi: 10.1158/0008-5472.CAN-17-3858

PubMed Abstract | Crossref Full Text | Google Scholar

33. Salem I, Alsalahi M, Chervoneva I, Aburto LD, Addya S, Ott GR, et al. The effects of CEP-37440, an inhibitor of focal adhesion kinase, in vitro and in vivo on inflammatory breast cancer cells. Breast Cancer Res. (2016) 18:37. doi: 10.1186/s13058-016-0694-4

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: lung cancer, metastasis, exosome, proteomics, metabolomic

Citation: He T, Jiang T, Sun X, Yang F, Zhang D, Yao S, Liao J and Wu X (2025) Integrated multi-omics landscape of non-small cell lung cancer with distant metastasis. Front. Immunol. 16:1560724. doi: 10.3389/fimmu.2025.1560724

Received: 14 January 2025; Accepted: 24 February 2025;
Published: 17 March 2025.

Edited by:

Bufu Tang, Fudan University, China

Reviewed by:

Lei Lei Liang, Chinese Academy of Medical Sciences and Peking Union Medical College, China
Fang Zhang, Jingling Hospital, China
Ziqiang Zhang, Shanghai Pudong Hospital, China

Copyright © 2025 He, Jiang, Sun, Yang, Zhang, Yao, Liao and Wu. 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: Jiangrong Liao, TGpyMzQwOUBzaW5hLmNvbQ==; Xueling Wu, d3V4dWVsaW5nNzZAMTI2LmNvbQ==

These authors have contributed equally to this work

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.

Research integrity at Frontiers

Man ultramarathon runner in the mountains he trains at sunset

94% of researchers rate our articles as excellent or good

Learn more about the work of our research integrity team to safeguard the quality of each article we publish.


Find out more