Low muscle mass in lung cancer is associated with an inflammatory and immunosuppressive tumor microenvironment
Journal of Translational Medicine volume 21, Article number: 116 (2023)
Computed tomographies (CT) are useful for identifying muscle loss in non-small lung cancer (NSCLC) cachectic patients. However, we lack consensus on the best cutoff point for pectoralis muscle loss. We aimed to characterize NSCLC patients based on muscularity, clinical data, and the transcriptional profile from the tumor microenvironment to build a cachexia classification model.
We used machine learning to generate a muscle loss prediction model, and the tumor's cellular and transcriptional profile was characterized in patients with low muscularity. First, we measured the pectoralis muscle area (PMA) of 211 treatment-naive NSCLC patients using CT available in The Cancer Imaging Archive. The cutoffs were established using machine learning algorithms (CART and Cutoff Finder) on PMA, clinical, and survival data. We evaluated the prediction model in a validation set (36 NSCLC). Tumor RNA-Seq (GSE103584) was used to profile the transcriptome and cellular composition based on digital cytometry.
CART demonstrated that a lower PMA was associated with a high risk of death (HR = 1.99). Cutoff Finder selected PMA cutoffs separating low-muscularity (LM) patients based on the risk of death (P-value = 0.003; discovery set). The cutoff presented 84% of success in classifying low muscle mass. The high risk of LM patients was also found in the validation set. Tumor RNA-Seq revealed 90 upregulated secretory genes in LM that potentially interact with muscle cell receptors. The LM upregulated genes enriched inflammatory biological processes. Digital cytometry revealed that LM patients presented high proportions of cytotoxic and exhausted CD8+ T cells.
Our prediction model identified cutoffs that distinguished patients with lower PMA and survival with an inflammatory and immunosuppressive TME enriched with inflammatory factors and CD8+ T cells.
Cachexia is a metabolic syndrome found in 50% of patients with lung cancer . It is associated with the risk of death in non-small cell lung cancer (NSCLC) —the most prevalent histological type . Low muscle mass in cachectic patients with advanced NSCLC is associated with pain and poor quality of life . Therefore, early detection of cachexia is of utmost importance in predicting NSCLC patient outcomes and guiding better treatment decisions.
Multiple screening tools are used to measure muscle loss . One feasible approach is computed tomography (CT), which is routinely performed for cancer diagnosis and follow-up . CT-based muscle quantification at the third lumbar vertebra (L3) level is commonly used to assess muscle mass in cancer patients . However, CT scans of NSCLC patients usually do not include L3 . Studies overcame this limitation by analyzing the pectoralis muscle area (PMA) [8, 9]. Our group demonstrated that PMA alone is a cachexia predictor in treatment-naive NSCLC patients . PMA quantification in lung cancer patients is associated with clinical outcomes; however, each study used different calculations and cutoffs to select cachectic patients [8, 9]. Although informative, the best PMA cutoff that classifies cachectic patients based on CTs remains to be defined.
Molecular markers also contribute to the clinical classification of cachexia; however, none are yet approved for clinical use. The synergic expression of tumor-derived cachexia-inducing factors (CIFs) is correlated with the prevalence of cachexia in different tumor types . For example, lung cancer cachexia is characterized by increased expression of CIFs that ultimately induce muscle wasting . NSCLC patients with low PMA present increased expression of IL8, IL6, and other critical CIFs . Cell types in the tumor microenvironment (TME) contribute distinctively by secreting cytokines or other CIFs that act on cell surface receptors of target tissues . Cancer samples can be divided into clinically relevant immune subtypes depending on TME composition and gene expression profile . Thus, computational deconvolution of the transcriptome from tumor bulk using digital cytometry and single-cell RNA-Seq (scRNA-Seq) can help infer the TME cellular proportions. TME analysis also has the potential to identify cells responsible for secreting CIF in cachectic patients.
The aim of this study was to build a machine learning (ML) model based on PMA and clinical data to identify NSCLC patients with low muscularity and poor prognosis. We additionally characterized the global secretome and TME immune fraction of these patients by utilizing tumor transcriptomes. Overall, the predictive model determined the best PMA cutoffs to select potentially cachectic patients. We demonstrated that these patients have a pro-inflammatory state characterized by high tumor gene expression of CIFs. Based on computational deconvolution of bulk transcriptomes, we also identified that patients with lower PMA and worse survival show a TME enriched with cytotoxic and exhausted CD8+ T cells.
Software, tools, and databases are described in Additional file 1.
CTs and clinical data were downloaded from The Cancer Imaging Archive (TCIA) database. The NSCLC-Radiogenomics collection contains 211 NSCLC patients [14, 15] (discovery set). CTs were taken at diagnosis, and patients were treated with surgery. We included 107 patients with primary tumors, presenting whole-body PET-CT before surgery, not subjected to chemo or radiotherapy, with arms raised during CT to avoid positioning bias, and overall survival data. An independent validation set with 36 primary NSCLC patients (all gave their informed consent prior to their inclusion in the study) from the Faculty of Medicine, São Paulo State University (UNESP)—Botucatu, Brazil (2012–2019) was used to confirm the prediction model (approved by the Faculty of Medicine Research Ethics Committee REB#45723921.0.0000.5411). Third-two patients from the validation set had survival data.
We analyzed the PMA using a single axial slice of the CT as previously described . We analyzed the total skeletal muscle cross-sectional area at the level of the L3 to classify patients with low muscle mass based on pre-established skeletal muscle indexes (SMI; men < 55 cm2/m2; women < 39 cm2/m2) . This classification was compared with the one generated by our model. We generated PMA-SMI (cm2/height2) to test whether height normalization affects our prediction model.
Cachexia classification model using machine learning
Predictive analysis based on a non-parametric risk prediction model was performed by the Classification And Regression Tree (CART) method using the regression analysis available in "rpart" using RStudio. We performed pruning to generate the tree with minimum × error (prune function; cp = 0.053). We used ten-fold cross-validation to assess the possibility of overfitting. We analyzed the age, tumor stage, PMA, and survival as variables. The imbalance between the two groups was addressed by a balanced resampling using random down-sampling. The Cutoff Finder package was independently applied to select the PMA cutoff that distinguishes patients' survival using R. Area Under the Curve-Receiver Operating Characteristics (AUC-ROC) analysis with a 95% confidence interval (CI) was used via the easyROC web tool to determine the diagnostic value of our PMA cutoffs in predicting low muscle mass. For this comparison, muscularity was defined using L3 SMI cutoffs. EasyROC was also used to test PMA-indexes in determining the LM and HM patients.
Tumor transcriptomic analysis
The tumor RNA-sequencing data (Illumina HiSeq 2500) of 130 subjects from NSCLC Radiogenomics collection were publicly available at GEO datasets (GSE103584). Of the 107 included patients, 46 had tumor transcriptomic data. We compared count data from LM patients (N = 9) with HM patients (N = 37) to perform differential gene expression analysis using the default setting of the BioJupies platform. The differentially expressed genes (DEGs) were selected according to fold change (FC) >|1.5| and P-value < 0.05. DEGs encoding secreted proteins were filtered based on the Human Protein Atlas database using the majority decision-based method for secreted proteins list of 2943 genes (Additional file 3: Table S1). Principal component analysis (PCA) was performed using the ClustVis webtool.
Enrichment and protein–protein interaction (PPI)
We used EnrichR to perform a functional enrichment analysis of DEGs. We selected the top 10 most enriched biological processes using gene ontology (GO) - 2018. PPI networks were conducted using the STRING tool. We considered experiments, databases, co-expression, neighborhood, and co-occurrence as active interaction sources. Only interactions with the highest confidence (interaction score > 0.9) were included, and disconnected nodes were omitted from the final network. A P-value < 0.01 was considered statistically significant. PPI network data was visualized and annotated using Cytoscape. The network was analyzed using the CytoNCA plugin, which calculated the node's betweenness centrality.
Tumor-muscle crosstalk prediction
We used the consensus list of ligand-receptors generated by Ramilowski et al.  to test if the skeletal muscles express receptors for the tumor-derived factors. We validated the gene expression of each predicted muscle receptor in human primary muscle cells from the FANTOM5 project, 491 muscle tissues from the GTEx portal, muscle tissues from the BioGPS website, muscle tissues from GeneAtlas U133A, and the gcrma human dataset. The alluvial diagram connecting the ligands to their receptors was generated using the SankeyMATIC tool.
Expression profile of secretome genes in malignant cells
To evaluate the abundance of DEGs found in the tumor of LM in malignant lung cells, we used transcriptomic data from the Cancer Cell Line Encyclopedia (CCLE). We selected eight cell lines with higher correlation with tumor tissues using the "TCGA-110CL Cell Line Panel". Considering that lung tumor organoids recapitulate the histology, gene expression, and genomic profile of the original tumor , we included lung tumor-derived organoids transcriptomic data available in the HCMI (Human Cancer Model Initiative) catalog (duplicates of HCM-CSHL-0058-C34). The selected cell lines and organoids are described in Additional file 3: Table S2.
Digital cytometry and scRNA-Seq
To evaluate whether immune infiltrating cells express the DEGs found in the LM, we performed digital cytometry analysis using the CIBERSORTx tool to impute the immune cell fractions of 22 cell types (LM22 matrix signature) from the discovery bulk RNA-Seq data. We applied the default settings and batch-correction. We also performed a gene expression imputation analysis to verify whether the over-expressed genes in LM patients are expressed in the six major leukocyte subsets (B cells, CD4+ T cells, CD8+ T cells, NKT cells, NK cells, and monocytes) of patients with NSCLC derived from scRNA-Seq data . The gene expression profile (GEP) was estimated by CIBERSORTx group mode analysis. In addition, the specific state of T cells (naive, cytotoxic, and exhausted) was evaluated using Immune Cell Abundance Identifier (ImmuCellAI).
Ligand–receptor interactions from gene expression
We used the computational framework CellPhoneDB to predict cell–cell communication and better understand the interplay of tumor immune cells using its repository of curated ligand–receptor interactions for scRNA-Seq. We used the default setting to select the statistically relevant interaction (P-value < 0.05) between the NSCLC T cells using the expression profile enriched in LM patients. This analysis was performed by re-analyzing the publicly available scRNA-Seq data from Guo et al.  (GSE99254).
Statistical analysis not previously reported was performed using GraphPad Prism® (GraphPad Software, v5.0, 2008 USA). The Mann–Whitney U-test was applied to compare the body composition variables between LM and HM patients (applied for non-parametric data). Chi-square analysis was used for categorical variables. Overall survival analysis was performed using the Kaplan–Meier method, and the log-rank test was used to compare the curves. Differences with a P-value < 0.05 were considered significant.
Predictive models identify the best PMA cutoff that classifies low muscle mass
We included 107 patients with NSCLC with a mean age of 68.7 years old (42 to 87). Sixty-eight were male, and 39 were female. Sex differences were considered when analyzing the PMA. We also identified differences according to the age range, especially in women; the older ones (> 65 years old) presented lower PMA than younger ones (≤ 65 years) (Additional file 2: Fig. S1A). PMA did not differ between young and older men (Additional file 2: Fig. S1A).
We normalized PMA with the z-score considering sex and age differences to make the groups comparable (Additional file 2: Fig. S1B). The CART analysis showed that patients with lower PMA are associated with a high risk of death followed by tumor stage (Fig. 1A). Using the Cutoff Finder, the same muscularity cutoff separated patients with poor survival and low PMA (low-muscularity patients, LM; n = 22) from patients with better survival and high PMA (high-muscularity patients, HM; n = 85) (Fig. 1B). Down-sampling for balanced comparisons using CART demonstrated models with similar cutoffs able to separate high-risk patients (Additional file 2: Fig. S1C). Using the L3 SMI as a muscularity parameter, we found that PMA cutoffs classify patients with low muscle mass with an AUC value of 0.84 (95% CI: 0.7–0.98; P < 0.001; Fig. 1C). Considering the relevance of height normalization, we also generated PMA-SMI. The LM group also presented a significant decrease in PMA-SMI (Table 1), and these values were consistent with the PMA classification for men and women (Additional file 2: Fig. S1D). Therefore, we conclude that both PMA and PMA height-normalized are good predictors.
By transforming the z-scored PMA values into absolute values, we established the PMA cutoffs to select LM patients according to age and sex: young male = PMA < 34.4 cm2; old male = PMA < 34.3 cm2; young female = PMA < 25.4 cm2; and old female = PMA < 18.6 cm2. LM patients had lower body mass index (BMI), SMI, and PMA-SMI. We observed increased smokers' frequency, high tumor grade, pleural invasion, and tumor recurrence in the LM patients (Table 1). An independent group of LM patients with NSCLC analyzed according to the PMA cutoffs generated during the discovery analysis confirmed the associations with poor prognosis (worse survival, recurrence, and smoking habits) (Table 1 and Fig. 1D). By combining the survival data of patients from the discovery and validation sets, we found a more significant log-rank P-value < 0.001 (Fig. 1E). Although not statistically significant, more patients with LM (31.3%) had KRAS mutations compared to HM (16.9%).
Tumors from patients with LM present a pro-inflammatory profile
We included nine LM patients and 37 HM patients with RNA-Seq data availability for the expression analysis and observed a similar expression profile using the PCA (Fig. 2A). These similarities resulted in considering the P-values < 0.05 to select DEGs for further analysis. We took into consideration the biological context of these DEGs to avoid misinterpretation due to false positives [20,21,22]. We identified 682 DEGs (332 were up and 350 were downregulated) in LM (Fig. 2B and Additional file 3: Table S3). Enrichment analysis of upregulated genes demonstrated leukocyte chemotaxis, cytokine/chemokine activity, and T-cell activation as the most enriched terms. The downregulated genes enriched terms associated with drug transport, cellular adhesion, and glucuronidation (Fig. 2C). We found 90 and 41 up and downregulated secretome genes in LM (Additional file 3: Table S3). We failed to identify significant interactions among the downregulated secretory genes in the PPI. Therefore, we focused on upregulated secretome genes that included several known cachexia mediators described in pan-cancer  and NSCLC  (Additional file 3: Table S4), such as IL6, IFNG, LIF, CCL2, and CSF3 (Fig. 2D). IL6 presented a central position in the network. CCL2 (MCP‐1) was the CIF with the highest abundance in LM tumors (Fig. 2E).
Tumor-derived factors are predicted to interact with skeletal muscle
To verify how the secretome profile from tumors of LM is potentially altering the skeletal muscle morphology, we constructed a tumor-muscle connectome by identifying the gene expression of the ligands produced by the tumors and their corresponding receptors in skeletal muscle cells and tissues (Fig. 3). We observed that several ligands have receptors expressed in abundance both in skeletal muscle cells and tissues. Interestingly, LIF, IFNG, and IL6 were the cytokines presenting receptors with high expression in muscle cells compared to other cytokines, such as CCLs and CXCLs, in which receptor expression is absent in muscle cell lines. Many receptors for CCLs and CXCLs are expressed in skeletal tissues, as demonstrated in the GTEx and BioGPS data, suggesting that these cytokines target muscle tissue cells other than myoblasts, myotubes, and satellite cells. This analysis ratifies that tumor-released factors may target skeletal muscle via specific receptors.
Malignant cells express low levels of inflammatory genes
Once we identified a cachexia-like expression profile in LM tumors, we next investigated the tumor cell types responsible for producing CIFs. First, to evaluate whether the malignant cells express the 332 upregulated genes found in the LM, we investigated the abundance of these DEGs using the expression profile of NSCLC cells (lines and organoids). We applied k-means clusterization and identified five clusters. We detected two clusters in secretome genes: cluster A represents abundant secretory genes in malignant cells, and cluster B represents low-level of secretory genes in malignant cells. The genes in cluster A enriched angiogenesis and leukocyte chemotaxis (THBS1 and VEGFC), and the genes present in cluster B are associated with inflammation. Among the genes not predicted to encode secreted proteins, we identified three clusters: clusters C and E, representing genes abundant in malignant cells associated with cell proliferation and mitosis, and cluster D, containing genes expressed at low levels in malignant cells and associated with T cell biology. Noteworthy, cluster E enriched the term "positive regulation of T cell apoptotic process" associated with IDO1 and CD274 (PDL1) (Fig. 4).
LM patients present high proportions of CD8+ T cells
Since we observed that many inflammation-associated genes are expressed at low levels by the malignant cells, we used CIBERSORTx to predict immune cells that contribute to the expression of these genes. LM patients presented high proportions of CD8+ T cells (Fig. 5A). To understand the molecular profile of CD8+ T cells in the TME, we used an independent dataset of NSCLC scRNA-Seq to analyze the expression profile of the 90 secretome genes identified in LM patients. We found CD8+ T cells expressing several of these genes in higher proportions than other immune cells (Fig. 5B), suggesting that CD8+ T cells may express secretome genes.
When plotting PCA using the expression profile enriched by CD8+ T cells, we better separated HM patients from LM patients (Fig. 5C). Cytotoxic and exhausted T cells were prevalent in LM compared to HM (Fig. 5D). Investigating the 332 upregulated genes in LM in another scRNA-Seq data from the T cell population of NSCLC using CellphoneDB, we predicted ligands (CCL3L3, CCL5, CCL3, and HLA-G) and receptors (CCR1 and LILRB2) communications among effectors, exhausted, pre-exhausted, and naive CD8+ T cells (Fig. 5E). We observed that exhausted and pre-exhausted CD8+ T cells presented increased interactions, specifically via CCR1 (exhausted)—CCL3L3 (pre-exhausted) and CCL3 (exhausted)—CCR1 (pre-exhausted) (Additional file 2: Figs. S2A and 5E). Moreover, genes expressed in exhausted tumor-infiltrating CD8+ T cells  were upregulated in LM patients, such as GZMB, VCAM1, IFNG, CCL3, FASLG, MYO7A, STMN1, SAMSN1, and ITGAE (Additional file 3: Table S5). CD8 and exhaustion scores are negatively correlated with PMA (Additional file 2: Fig. S2B).
We showed that a lower PMA evaluated in CTs used for NSCLC diagnosis before tumor resection predicts worse outcomes in patients with NSCLC. We also proposed PMA cutoffs for women and men at young and older ages based on ML analysis. Tumor transcriptome demonstrated that LM patients present a cachexia-like profile highly pro-inflammatory and express CIFs. Digital cytometry based on transcriptomic data revealed that TME is enriched with a suppressive immune microenvironment, mainly composed of CD8+ T cells. The prediction of ligand-receptor revealed candidates that may ultimately lead to muscle wasting.
Cancer cachexia classification remains an important clinical challenge. The milestones for the classification of cachexia are weight loss > 5% over the past 6 months or BMI < 20 and any degree of weight loss > 2% or lumbar SMI determined by CT imaging (men < 55 cm2/m2; women < 39 cm2/m2) . However, CTs are often taken in the thoracic region of lung cancer patients . To overcome this limitation, we constructed a cachexia classification model using ML to evaluate clinical parameters routinely available for NSCLC patients. We found that lower PMA can predict worse survival. The PMA (z-scored) cutoff identified in CART was also validated by independent analysis. These results demonstrated that PMA cutoffs are strong predictors of the poor prognosis of NSCLC patients with low muscle mass. Despite not demonstrating an association with worse overall survival, our research group previously demonstrated PMA cutoffs as a potential cachexia predictor in NSCLC patients. The cutoffs found in this study were PMA < 32.2 cm2 and < 21 cm2 for men and women, respectively, similar to those found herein using robust ML. CART was recently used to predict cachexia risk levels in incurable cancer patients with more than 88% of precision and sensitivity . These results reinforce the importance of using ML approaches applied to clinical features of cancer patients to classify cachexia.
We found a high frequency of smokers in the LM group, which corroborates the literature that indicates an increase in muscle loss in tobacco smokers  and confirms the significant role of nicotine in altering the tumor microenvironment and prognosis, as demonstrated for cachectic patients with pancreatic cancer . A high risk of tumor recurrence was found in LM patients and agreed with other cachexia studies . KRAS mutation was frequently found in our LM patients, and evidence suggests KRAS-mutated tumors are commonly associated with lung cancer cachexia . Since adenocarcinomas—the tumor with the highest prevalence of the KRAS mutation —occur at similar rates in LM (72.7%) and HM (85.7%) groups, we confirmed that KRAS mutations were associated with muscularity rather than tumor histopathology. Thus, studies investigating differences in the cachexia phenotype of preclinical models inoculated with KRAS-mutated NSCLC cell lines may help to understand the impact of KRAS mutation in the syndrome. These essential features can help guide cancer cachexia definition towards precision medicine in the future.
Our cachexia classification was validated by the increase of CIFs genes in the tumor of LM patients, such as LIF, IL6, IFNG, and CCL2 . We also identified increased expression of genes described in NSCLC patients with lower PMA, such as NCAM1, SCG2, CSF3, and CCL8 . To understand how the secretome influences skeletal muscle phenotype and function, we characterized the landscape of cell–cell interaction between tumor and skeletal muscle cells using transcriptomic data from both tissues. This approach has a higher resolution than proteomics and provides insights into cell and tissue communication . We found receptors for CIFs expressed at high levels in skeletal muscle cells, such as IFNG receptor (IFNGR2) and LIF receptors (LIFR and IL6ST), which can directly trigger atrophy . The gp130 family members, IL-11, IL-6, OSM, and LIF (upregulated in the LM) potentially impact skeletal muscle wasting by inducing the STAT3 signaling pathway in the preclinical model of cancer cachexia . These results confirm the importance of these tumor-derived factors in NSCLC cachexia. Previous studies have individually inhibited the classical CIFs described herein . Although preclinical studies inhibiting CIFs have shown some efficiency, these factors have not provided the same benefit in clinical trials . Our tumor-muscle interactome revealed a combinatorial action of secreted mediators from cells within the TME that directly acts on muscle wasting. More studies are needed to describe the causality of these ligand–receptors interactions and their combinatorial action in inducing muscle atrophy in NSCLC cachexia.
To investigate the tumor cells responsible for CIFs secretion, we further evaluated the expression profile of neoplastic lung cells. We found that inflammation-related genes are lowly expressed, while genes associated with vessel permeability and T-cell motility, recruitment, and activation are increased . We identified enrichment in genes associated with T cell apoptosis, specifically the genes CD274 (encoding PDL1) and IDO1. The processes involved in tumor immune evasion through IFN-γ (enriched in LM) have been reviewed . Anti-tumor cells secrete IFN-γ that induces genomic instability and/or immune evasive gene expression signature in malignant cells (including PDL1 and IDO1) . PDL1 is the ligand of the inhibitory receptor PD-1, responsible for CD8 dysfunction in NSCLC . IDO1 produces local depletion of amino acid tryptophan, contributing to the immunosuppressive microenvironment . These results corroborate our digital cytometry that demonstrates an increase in CD8+ T cells in TME of LM patients enriched with exhausted cells but also enriched with cytotoxic T cells.
Although cytotoxic CD8+ T cells infiltrated in TME are frequently associated with intense anticancer activity , the CD8+ T cells' signal in tumor tissues is not always related to better outcomes. Using TCGA data, six immune subtypes were identified. Subtype 2 (C2) is characterized by IFN-γ predominance, the highest content of M1 macrophages and CD8+ T cells. Cancer patients presenting subtype C2 showed less favorable survival despite having the highest lymphocyte infiltration . The authors proposed that these tumors are more aggressive or escape immune recognition by pre-existing remodeling. However, our results raise the possibility that patients with C2 tumors have less favorable survival due to cachexia. In fact, CD8+ T cells induce cachexia and adipose tissue wasting in experimental models of chronic infection via type I IFN .
Considering the immunosuppressive profile of LM tumors and immunotherapy promises to restore T cell cytotoxicity and its anticancer properties, immunotherapy would probably benefit NSCLC cachectic patients. Cachectic patients with NSCLC are resistant to immunotherapy ; thus, restoring the cytotoxicity of CD8+ T cells may be crucial for treatment response. Indeed, clinical trials testing an agent to reverse cachexia combined with an immune checkpoint inhibitor to treat advanced lung cancer patients is relevant since the treatment with ghrelin analog affects T cell development and proliferation . Therefore, more studies are needed to verify whether immunotherapy combined with CIF-targeted drugs has the potential to reverse cachexia and increase the immunotherapy response.
Strengths and limitations
This study is innovative and valuable for NSCLC cachexia. Here, we integrated imaging, clinical, and transcriptomic data to capture the complex interplay between tumor biology and host response associated with cachexia by applying bioinformatics and ML analyses. The predictive model generated herein was validated in an independent set of NSCLC, demonstrating the great potential of our methodology. However, it has limitations. Patients lack information on other comorbidities that may influence body composition and clinical outcomes. Considering that LM patients had a higher frequency of smokers, it is rational to consider that the skeletal muscle alterations identified are triggered by the combination of smoking and CIFs. However, it is difficult to determine the individual impact of each of these variables. We also failed to apply the golden standard definitions of cachexia, such as weight loss, in those patients. Moreover, even though transcriptomic data have high resolution and produce more robust information than proteomic data, the results should be interpreted carefully considering the complex transcriptional regulation in human cells. Therefore, the results of the secretome and ligand-receptor predictions presented in this study need experimental validation in cachectic patients. Furthermore, the immune TME profile of patients with LM patients must be validated by histological examination to check the correlation between CD8 exhaustion and cachexia in a larger number of samples.
We used systems biology to classify cachexia in treatment-naive NSCLC patients and characterized the TME immune fraction from the LM patients to understand molecular mechanisms by which the tumor induces skeletal muscle atrophy and influences poor outcomes. These results may contribute to the development of therapeutic strategies aiming to decrease muscle atrophy and deleterious inflammation, thus increasing the survival and quality of life of patients with NSCLC cachexia. Overall, the predictive model identified patients with lower muscularity and worse survival, who showed a TME with an inflammatory and immunosuppressive profile, expressing CIFs and potentially enriched with CD8+ T cells. These findings open new perspectives on the early cachexia diagnosis for NSCLC patients. In addition, this study paves the way for the identification of tumor cell-specific cachexia mediators and biomarkers. In the future, we hope that these results can serve as a basis for generating effective treatments for cachexia in patients with lung cancer.
Availability of data and materials
The datasets presented in this study can be found in online public repositories available under the accession numbers GSE103584 and GSE99254. The results shown here are in part based on data generated by the Genotype-Tissue Expression project (GTEx) (https://gtexportal.org/).
Non-small cell lung cancer
Pectoralis muscle area
The Cancer Imaging Archive
Differentially expressed genes
Body mass index
Skeletal muscle index
Single-cell RNA sequencing
Baracos VE, Martin L, Korc M, Guttridge DC, Fearon KCH. Cancer-associated cachexia. Nat Rev Dis Primer. 2018;4:17105.
Kimura M, Naito T, Kenmotsu H, Taira T, Wakuda K, Oyakawa T, et al. Prognostic impact of cancer cachexia in patients with advanced non-small cell lung cancer. Support Care Cancer. 2015;23:1699–708.
Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature. 2018;553:446–54.
Fearon K, Strasser F, Anker SD, Bosaeus I, Bruera E, Fainsinger RL, et al. Definition and classification of cancer cachexia: an international consensus. Lancet Oncol. 2011;12:489–95.
Buckinx F, Landi F, Cesari M, Fielding RA, Visser M, Engelke K, et al. Pitfalls in the measurement of muscle mass: a need for a reference standard: measurement of muscle mass. J Cachexia Sarcopenia Muscle. 2018;9:269–78.
Prado CM, Birdsell LA, Baracos VE. The emerging role of computerized tomography in assessing cancer cachexia. Curr Opin Support Palliat Care. 2009;3:269–75.
Swensen SJ, Jett JR, Hartman TE, Midthun DE, Sloan JA, Sykes A-M, et al. Lung Cancer Screening with CT: Mayo Clinic Experience. Radiology. 2003;226:756–61.
Kinsey CM, San José Estépar R, van der Velden J, Cole BF, Christiani DC, Washko GR. Lower pectoralis muscle area is associated with a worse overall survival in non-small cell lung cancer. Cancer Epidemiol Biomarkers Prev. 2017;26:38–43.
Cury SS, de Moraes D, Freire PP, de Oliveira G, Marques DVP, Fernandez GJ, et al. Tumor transcriptome reveals high expression of IL-8 in non-small cell lung cancer patients with low pectoralis muscle area and reduced survival. Cancers. 2019;11:1251.
Freire PP, Fernandez GJ, Moraes D, Cury SS, Dal Pai-Silva M, Reis PP, et al. The expression landscape of cachexia-inducing factors in human cancers. J Cachexia Sarcopenia Muscle. 2020;11:947–61.
Sørensen J. Lung Cancer Cachexia: Can Molecular Understanding Guide Clinical Management? Integr Cancer Ther. 2018;17:1000–8.
Argilés JM, Stemmler B, López-Soriano FJ, Busquets S. Inter-tissue communication in cancer cachexia. Nat Rev Endocrinol. 2019;15:9–20.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang T-H, et al. The immune landscape of cancer. Immunity. 2018;48:812-830.e14.
Bakr S, Gevaert O, Echegaray S, Ayers K, Zhou M, Shafiq M, et al. A radiogenomic dataset of non-small cell lung cancer. Sci Data. 2018;5: 180202.
Gevaert O, Xu J, Hoang CD, Leung AN, Xu Y, Quon A, et al. Non-small cell lung cancer: identifying prognostic imaging biomarkers by leveraging public gene expression microarray data—methods and preliminary results. Radiology. 2012;264:387–96.
Ramilowski JA, Goldberg T, Harshbarger J, Kloppmann E, Lizio M, Satagopam VP, et al. A draft network of ligand–receptor-mediated multicellular signalling in human. Nat Commun. 2015;6:7866.
Kim M, Mun H, Sung CO, Cho EJ, Jeon H-J, Chun S-M, et al. Patient-derived lung cancer organoids as in vitro cancer models for therapeutic screening. Nat Commun. 2019;10:3991.
Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. 2019;37:773–82.
Guo X, Zhang Y, Zheng L, Zheng C, Song J, Zhang Q, et al. Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing. Nat Med. 2018;24:978–85.
Zhao R, Pan Z, Li B, Zhao S, Zhang S, Qi Y, et al. Comprehensive analysis of the tumor immune microenvironment landscape in glioblastoma reveals tumor heterogeneity and implications for prognosis and immunotherapy. Front Immunol. 2022;13: 820673.
Cui J, Li F, Wang G, Fang X, Puett JD, Xu Y. Gene-expression signatures can distinguish gastric cancer grades and stages. PLoS ONE. 2011;6: e17819.
Lecante LL, Gaye B, Delbes G. Impact of in utero rat exposure to 17Alpha-ethinylestradiol or genistein on testicular development and germ cell gene expression. Front Toxicol. 2022;4: 893050.
Vagnildhaug OM, Brunelli C, Hjermstad MJ, Strasser F, Baracos V, Wilcock A, et al. A prospective study examining cachexia predictors in patients with incurable cancer. BMC Palliat Care. 2019;18:46.
Derksen JWG, Kurk SA, Oskam MJ, Peeters PHM, Punt CJA, Koopman M, et al. Factors contributing to cancer-related muscle wasting during first-line systemic treatment for metastatic colorectal cancer. JNCI Cancer Spectr. 2019;3: pkz014.
Underwood PW, Zhang DY, Cameron ME, Gerber MH, Delitto D, Maduka MU, et al. Nicotine induces IL-8 secretion from pancreatic cancer stroma and worsens cancer-induced cachexia. Cancers. 2020;12:329.
Marhold M, Topakian T, Unseld M. Sarcopenia in cancer—a focus on elderly cancer patients. Memo - Mag Eur Med Oncol. 2020. https://doi.org/10.1007/s12254-020-00637-6.
Shiono M, Huang K, Downey RJ, Consul N, Villanueva N, Beck K, et al. An analysis of the relationship between metastases and cachexia in lung cancer patients. Cancer Med. 2016;5:2641–8.
Garrido P, Olmedo ME, Gómez A, Paz Ares L, López-Ríos F, Rosa-Rosa JM, et al. Treating KRAS-mutant NSCLC: latest evidence and clinical consequences. Ther Adv Med Oncol. 2017;9:589–97.
Armingol E, Officer A, Harismendy O, Lewis NE. Deciphering cell–cell interactions and communication from gene expression. Nat Rev Genet. 2020;71:88.
Seto DN, Kandarian SC, Jackman RW. A key role for leukemia inhibitory factor in C26 cancer cachexia. J Biol Chem. 2015;290:19976–86.
Cao Z, Scott AM, Hoogenraad NJ, Osellame LD. Mediators and clinical treatment for cancer cachexia: a systematic review. JCSM Rapid Commun. 2021;4:166–86.
Kadambi A, Mouta Carreira C, Yun CO, Padera TP, Dolmans DE, Carmeliet P, et al. Vascular endothelial growth factor (VEGF)-C differentially affects tumor vascular function and leukocyte recruitment: role of VEGF-receptor 2 and host VEGF-A. Cancer Res. 2001;61:2404–8.
Mojic M, Takeda K, Hayakawa Y. The dark side of IFN-γ: its role in promoting cancer immunoevasion. Int J Mol Sci. 2017;19:89.
Zhang Y, Huang S, Gong D, Qin Y, Shen Q. Programmed death-1 upregulation is correlated with dysfunction of tumor-infiltrating CD8+ T lymphocytes in human non-small cell lung cancer. Cell Mol Immunol. 2010;7:389–95.
Raskov H, Orhan A, Christensen JP, Gögenur I. Cytotoxic CD8+ T cells in cancer and cancer immunotherapy. Br J Cancer. 2020;124(2):359–67.
Baazim H, Schweiger M, Moschinger M, Xu H, Scherer T, Popa A, et al. CD8+ T cells induce cachexia during chronic viral infection. Nat Immunol. 2019;20:701–10.
Philip M, Fairchild L, Sun L, Horste EL, Camara S, Shakiba M, et al. Chromatin states define tumour-specific T cell dysfunction and reprogramming. Nature. 2017;545:452–6.
Bonomi P, Fidler MJ, Shah P, Borgia J. Theoretical and practical implications of treating cachexia in advanced lung cancer patients. Cancers. 2019;11:1619.
We thank Dr. Ivan Rodrigo Wolf for the help with machine learning analysis.
This study was supported by the São Paulo Research Foundation—FAPESP (grant 2020/03854-4 to RFC) and by the National Council for Scientific and Technological Development, CNPq (#141601/2019-1 to SSC and #311530/2019-2 to RFC).
Ethics approval and consent to participate
Human studies have been performed in accordance with the national laws and ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments. The study has been performed in accordance with the national laws and ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments. The use of computed tomographies from 36 primary NSCLC patients from the Faculty of Medicine, São Paulo State University (UNESP) - Botucatu, Brazil (2012–2019) was approved by the Faculty of Medicine Research Ethics Committee (REB#45723921.0.0000.5411). All patients gave their informed consent prior to being included in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Information regarding software, tools, and databases used in the methodology.
Pectoralis muscle area (PMA) cut-offs predict cachexia. A) Boxplot of PMA (mm2) from 38 young males, 97 old males, 34 young females, and 42 old females. B) Boxplot of the same measurements from panel A z-score normalized for gender and age. C) Decision tree generated using Classification And Regression Trees (CART) and Kaplan–Meier survival analyses for down-sampling models. The bottom boxes indicate hazard ratios, the number of patients at risk in each leaf, and the percentage of patients in each leaf. Red-curve: PMA > cut-off and low-risk group; Black-curve: PMA < cut-off and high-risk group. D) Area under the ROC curve (AUC-ROC) that demonstrates the specificity and sensitivity of the PMA skeletal muscle index (PMA cm2/height2) in indicating sarcopenia (determined by the PMA cut-offs generated in Fig. 1A). Male (n = 66) and female (n = 23) groups included both young and older patients that had height information available. Figure S2. Exhausted CD8+ T cells interact with pre-exhausted CD8+ T cells through tumor-derived factors identified in low-muscularity patients and negatively correlate with pectoralis muscle areas (PMA). A) Heatmap demonstrating the number of interactions between each CD8 cell subset. CD8_C6-LAYN is the exhausted subtype, while CD8_C4-GZMK corresponds to the pre-exhausted subtype. Analysis was performed using CellPhoneDB using the single-cell RNAseq data publicly available at GSE99254. B) Pearson correlation plot showing a negative correlation between CD8 and exhaustion scores (identified in the digital cytometry analysis) and PMA. BMI: body mass index.
Human Protein Atlas (HPA) list of 2943 secretory genes predicted using majority decision-based method for secreted proteins (MDSEC). Table S2. Lung cancer cell lines and lung tumor 3D organoids. Table S3. List of differentially expressed genes in low-muscularity patients compared to high-muscularity patients. Table S4. List of cachexia-inducing factors described in the literature. Table S5. List of genes specifically expressed in exhausted tumor CD8+ T cells according to Guo et al., 2018.
About this article
Cite this article
Cury, S.S., de Moraes, D., Oliveira, J.S. et al. Low muscle mass in lung cancer is associated with an inflammatory and immunosuppressive tumor microenvironment. J Transl Med 21, 116 (2023). https://doi.org/10.1186/s12967-023-03901-5