Immune phenotypic linkage between colorectal cancer and liver metastasis

Unveiling major histocompatibility complex-mediated pan-cancer immune features by integrated single-cell and bulk RNA sequencing

Immune checkpoint inhibitors (ICIs) have transformed cancer therapy, yet persistent challenges such as low response rate and significant heterogeneity necessitate attention. The pivotal role of the major histocompatibility complex (MHC) in ICI efficacy, its intricate impacts and potentials as a prognostic marker, warrants comprehensive exploration. This study integrates single-cell RNA sequencing (scRNA-seq), bulk RNA-seq, and spatial transcriptomic analyses to unveil pan-cancer immune characteristics governed by the MHC transcriptional feature (MHC.sig). Developed through scRNA-seq analysis of 663,760 cells across diverse cohorts and validated in 30 solid cancer types, the MHC.sig demonstrates a robust correlation between immune-related genes and infiltrating immune cells, highlighting its potential as a universal pan-cancer marker for anti-tumor immunity. Screening the MHC.sig for therapeutic targets using CRISPR data identifies potential genes for immune therapy synergy and validates its predictive efficacy for ICIs responsiveness across diverse datasets and cancer types. Finally, analysis of cellular communication patterns reveals interactions between C1QC+macrophages and malignant cells, providing insights into potential therapeutic agents and their sensitivity characteristics. This comprehensive analysis positions the MHC.sig as a promising marker for predicting immune therapy outcomes and guiding combinatorial therapeutic strategies.

Single-cell RNA sequencing and cell-cell communication analysis reveal tumor microenvironment associated with chemotherapy responsiveness in ovarian cancer

BACKGROUND

To investigate the impact of the tumor microenvironment (TME) on the responsiveness to chemotherapy in ovarian cancer (OV).

METHODS

We integrated single cell RNA-seq datasets of OV containing chemo-response information, and characterize their clusters based on different TME sections. We focus on analyzing cell-cell communication to elaborate on the mechanisms by which different components of the TME directly influence the chemo-response of tumor cells.

RESULTS

scRNA-seq datasets were annotated according to specific markers for different cell types. Differential analysis of malignant epithelial cells revealed that chemoresistance was associated with the TME. Notably, distinct TME components exhibited varying effects on chemoresistance. Enriched SPP1+ tumor-associated macrophages in chemo-resistant patients could promote chemoresistance through SPP1 binding to CD44 on tumor cells. Additionally, the overexpression of THBS2 in stromal cells could promote chemoresistance through binding with CD47 on tumor cells. In contrast, GZMA in the lymphocytes could downregulate the expression of PARD3 through direct interaction with PARD3, thereby attenuating chemoresistance in tumor cells.

CONCLUSION

Our study indicates that the non-tumor cell components of the TME (e.g. SPP1+ TAMs, stromal cells and lymphocytes) can directly impact the chemo-response of OV and targeting the TME was potentially crucial in chemotherapy of OV.

Pan-cancer analysis and single-cell analysis reveals FAM110B as a potential target for survival and immunotherapy

Background: FAM110B belongs to the family that has a 110 sequence similarity (FAM110) and is located in the centrosome and mitotic spindle. FAM110B has been linked to tumor cell growth in earlier research. Uncertainty exists regarding FAM110B's function within the tumor microenvironment is unclear as well as pan-cancer. Methods: In order to assess the variation in FAM110B expression within normal and pan-cancer tissues, we combined the TCGA and GTEx databases. The cBioPortal database and the GSCALite platform were used to examine the variation in genome and methylation alteration of FAM110B. Cox regression, Kaplan-Meier, and SangerBox were employed to examine the clinical features and prognosis of FAM110B and pan-cancer. The purpose of the correlational research was to investigate the associations within immunerelated genes, tumor mutation burden, microsatellite instability, immune-related genes, and immunological checkpoints and FAM110B expression. ESTIMATE, EPIC, QUANTISEQ, and MCPCOUNTER methods were used to calculate the interaction among FAM110B expression as well as the tumor immune microenvironment. The immunoinfiltration and function of FAM110B were analyzed by single-cell databases (TISCH and CancerSEA). Finally, we evaluated the sensitivity of FAM110B to small-molecule medications through GDSC and CTRP databases. Results: The transcription and protein expression of FAM110B varies significantly throughout cancer types, and this has predictive value for the prognosis of some tumors; including brain lower grade glioma (LGG), stomach adenocarcinoma (STAD), pancreatic adenocarcinoma (PAAD), etc. In the tumor microenvironment, the expression level of FAM110B was associated with immune cell infiltration, immune checkpoint immune regulatory genes, tumor mutational burden, and microsatellite fragility to a certain extent. Conclusion: This work investigates the possibility of utility of FAM110B as a marker to forecast pan-cancer immunotherapy response, providing a theoretical basis for cancer therapy.

Single-Cell Analyses Reveal the Metabolic Heterogeneity and Plasticity of the Tumor Microenvironment during Head and Neck Squamous Cell Carcinoma Progression

Metabolic reprogramming is a hallmark of cancer. In addition to metabolic alterations in the tumor cells, multiple other metabolically active cell types in the tumor microenvironment (TME) contribute to the emergence of a tumor-specific metabolic milieu. Here, we defined the metabolic landscape of the TME during the progression of head and neck squamous cell carcinoma (HNSCC) by performing single-cell RNA sequencing on 26 human patient specimens, including normal tissue, precancerous lesions, early stage cancer, advanced-stage cancer, lymph node metastases, and recurrent tumors. The analysis revealed substantial heterogeneity at the transcriptional, developmental, metabolic, and functional levels in different cell types. SPP1+ macrophages were identified as a protumor and prometastatic macrophage subtype with high fructose and mannose metabolism, which was further substantiated by integrative analysis and validation experiments. An inhibitor of fructose metabolism reduced the proportion of SPP1+ macrophages, reshaped the immunosuppressive TME, and suppressed tumor growth. In conclusion, this work delineated the metabolic landscape of HNSCC at a single-cell resolution and identified fructose metabolism as a key metabolic feature of a protumor macrophage subpopulation. Significance: Fructose and mannose metabolism is a metabolic feature of a protumor and prometastasis macrophage subtype and can be targeted to reprogram macrophages and the microenvironment of head and neck squamous cell carcinoma.

The prognostic effect of tumor-associated macrophages in stage I-III colorectal cancer depends on T cell infiltration

BACKGROUND

Tumor-associated macrophages (TAMs) are associated with unfavorable patient prognosis in many cancer types. However, TAMs are a heterogeneous cell population and subsets have been shown to activate tumor-infiltrating T cells and confer a good patient prognosis. Data on the prognostic value of TAMs in colorectal cancer are conflicting. We investigated the prognostic effect of TAMs in relation to tumor-infiltrating T cells in colorectal cancers.

METHODS

The TAM markers CD68 and CD163 were analyzed by multiplex fluorescence immunohistochemistry and digital image analysis on tissue microarrays of 1720 primary colorectal cancers. TAM density in the tumor stroma was scored in relation to T cell density (stromal CD3+ and epithelial CD8+ cells) and analyzed in Cox proportional hazards models of 5-year relapse-free survival. Multivariable survival models included clinicopathological factors, MSI status and BRAFV600E mutation status.

RESULTS

High TAM density was associated with a favorable 5-year relapse-free survival in a multivariable model of patients with stage I-III tumors (p = 0.004, hazard ratio 0.94, 95% confidence interval 0.90-0.98). However, the prognostic effect was dependent on tumoral T-cell density. High TAM density was associated with a good prognosis in patients who also had high T-cell levels in their tumors, while high TAM density was associated with poorer prognosis in patients with low T-cell levels (pinteraction = 0.0006). This prognostic heterogeneity was found for microsatellite stable tumors separately.

CONCLUSIONS

This study supported a phenotypic heterogeneity of TAMs in colorectal cancer, and showed that combined tumor immunophenotyping of multiple immune cell types improved the prediction of patient prognosis.

Single-cell profiling reveals altered immune landscape and impaired NK cell function in gastric cancer liver metastasis

Gastric cancer (GC) is a substantial global health concern, and the development of liver metastasis (LM) in GC represents a critical stage linked to unfavorable patient prognoses. In this study, we employed single-cell RNA sequencing (scRNA-seq) to investigate the immune landscape of GC liver metastasis, revealing several immuno-suppressive components within the tumor immune microenvironment (TIM). Our findings unveiled an increased presence of cancer-associated fibroblasts (CAFs), myeloid-derived suppressor cell (MDSC)-like macrophages, tumor-associated macrophage (TAM)-like macrophages, and naive T cells, while conventional dendritic cells (cDCs) and effector CD8 T cells declined in LM. Additionally, we identified two distinct natural killer (NK) cell clusters exhibiting differential cytotoxicity-related gene expression, with cytotoxic NK cells notably reduced in LM. Strikingly, TGFβ was identified as an inducer of NK cell dysfunction, potentially contributing to immune evasion and tumor metastasis. In preclinical LM models, the combined approach of inhibiting TGFβ and transferring NK cells exhibited a synergistic impact, resulting in a significant reduction in liver metastasis. This work highlights the importance of understanding the complex immune dynamics within GC liver metastasis and presents a promising strategy combining TGFβ inhibition and NK-based immunotherapy to improve patient outcomes.

Advances in immune regulation of the G protein-coupled estrogen receptor

Estrogen and related receptors have been shown to have a significant impact on human development, reproduction, metabolism and immune regulation and to play a critical role in tumor development and treatment. Traditionally, the nuclear estrogen receptors (nERs) ERα and ERβ have been thought to be involved in mediating the estrogenic effects. However, our group and others have previously demonstrated that the G protein-coupled estrogen receptor (GPER) is the third independent ER, and estrogen signaling mediated by GPER is known to play an important role in normal physiology and a variety of abnormal diseases. Interestingly, recent studies have progressively revealed GPER involvement in the maintenance of the normal immune system, abnormal immune diseases, and inflammatory lesions, which may be of significant clinical value primarily in the immunotherapy of tumors. In this article, we review current advances in GPER-related immunomodulators and provide a theoretical basis and potential clinical targets to ameliorate immune-related diseases and immunotherapy for tumors.

Kombo knife combined with sorafenib in liver cancer treatment: Efficacy and safety under immune function influence

BACKGROUND

In the quest to manage hepatocellular carcinoma (HCC), the focus has shifted to a more holistic approach encompassing both data analytics and innovative treatments. Analyzing rich data resources, such as the cancer genome atlas (TCGA), and examining progressive therapies can potentially reshape the trajectory of HCC treatment.

AIM

To elucidate the immunological genes and the underlying mechanism of the combined Kombo knife and sorafenib regimen for HCC by analyzing data from TCGA and machine learning data.

METHODS

Immune attributes were evaluated via TCGA's postablation HCC RNA sequencing data. Using weighted gene coexpression network analysis and machine learning, we identified genes with high prognostic value. The therapeutic landscape and safety metrics of the integrated treatment were critically evaluated across cellular and animal models.

RESULTS

Immune genes-specifically, peptidylprolyl isomerase A and solute carrier family 29 member 3-emerged as significant prognostic markers. Enhanced therapeutic outcomes, such as prolonged progression-free survival and an elevated overall response rate, characterize the combined approach, with peripheral blood mononuclear cells displaying potent effects on HCC dynamics.

CONCLUSION

The combination of Kombo knife with sorafenib is an innovative HCC treatment modality anchored in immune-centric strategies.

Biomarkers to predict efficacy of immune checkpoint inhibitors in colorectal cancer patients: a systematic review and meta-analysis

Immune checkpoint inhibitors (ICIs) are approved to treat colorectal cancer (CRC) with mismatch-repair gene deficiency, but the response rate remains low. Value of current biomarkers to predict CRC patients' response to ICIs is unclear due to heterogeneous study designs and small sample sizes. Here, we aim to assess and quantify the magnitude of multiple biomarkers for predicting the efficacy of ICIs in CRC patients. We systematically searched MEDLINE, Embase, the Cochrane Library, and Web of Science databases (to June 2023) for clinical studies examining biomarkers for efficacy of ICIs in CRC patients. Random-effect models were performed for meta-analysis. We pooled odds ratio (OR) and hazard ratio (HR) with 95% confidence interval (CI) for biomarkers predicting response rate and survival. 36 studies with 1867 patients were included in systematic review. We found that a lower pre-treatment blood neutrophil-to-lymphocyte ratio (n=4, HR 0.37, 95%CI 0.21-0.67) predicts good prognosis, higher tumor mutation burden (n=10, OR 4.83, 95%CI 2.16-10.78) predicts response to ICIs, and liver metastasis (n=16, OR 0.32, 95%CI 0.16-0.63) indicates resistance to ICIs, especially when combined with VEGFR inhibitors. But the predictive value of tumor PD-L1 expression (n=9, OR 1.01, 95%CI 0.48-2.14) was insignificant in CRC. Blood neutrophil-to-lymphocyte ratio, tumor mutation burden, and liver metastasis, but not tumor PD-L1 expression, function as significant biomarkers to predict efficacy of ICIs in CRC patients. These findings help stratify CRC patients suitable for ICI treatments, improving efficacy of immunotherapy through precise patient management. (PROSPERO, CRD42022346716).

Distinctive multicellular immunosuppressive hubs confer different intervention strategies for left- and right-sided colon cancers

Primary colon cancers arising from the left and right sides exhibit distinct clinical and molecular characteristics. Sidedness-associated heterogeneity relies intricately on the oncogenic properties of cancer cells and multicellular interactions in tumor microenvironments. Here, combining transcriptomic profiling of 426,863 single cells from 105 colon cancer patients and validation with spatial transcriptomics and large-scale histological analysis, we capture common transcriptional heterogeneity patterns between left- and right-sided malignant epithelia through delineating two side-specific expression meta-programs. The proliferation stemness meta-program is notably enriched in left-sided malignant epithelia that colocalize with Mph-PLTP cells, activated regulatory T cells (Tregs), and exhausted CD8-LAYN cells, constituting the glucose metabolism reprogramming niche. The immune secretory (IS) meta-program exhibits specific enrichment in right-sided malignant epithelia, especially in smoking patients with right-sided colon cancer. The IShigh malignant epithelia spatially localize in hypoxic regions and facilitate immune evasion through attenuating Mph-SPP1 cell antigen presentation and recruiting innate-like cytotoxicity-reduced CD8-CD161 cells.

Heterogeneity of myeloid cells in common cancers: Single cell insights and targeting strategies

Tumor microenvironment (TME), is characterized by a complex and heterogenous composition involving a substantial population of immune cells. Myeloid cells comprising over half of the solid tumor mass, are undoubtedly one of the most prominent cell populations associated with tumors. Studies have unambiguously established that myeloid cells play a key role in tumor development, including immune suppression, pro-inflammation, promote tumor metastasis and angiogenesis, for example, tumor-associated macrophages promote tumor progression in a variety of common tumors, including lung cancer, through direct or indirect interactions with the TME. However, due to previous technological constraints, research on myeloid cells often tended to be conducted as studies with low throughput and limited resolution. For example, the conventional categorization of macrophages into M1-like and M2-like subsets based solely on their anti-tumor and pro-tumor roles has disregarded their continuum of states, resulting in an inadequate analysis of the high heterogeneity characterizing myeloid cells. The widespread adoption of single-cell RNA sequencing (scRNA-seq) in tumor immunology has propelled researchers into a new realm of understanding, leading to the establishment of novel subsets and targets. In this review, the origin of myeloid cells in high-incidence cancers, the functions of myeloid cell subsets examined through traditional and single-cell perspectives, as well as specific targeting strategies, are comprehensively outlined. As a result of this endeavor, we will gain a better understanding of myeloid cell heterogeneity, as well as contribute to the development of new therapeutic approaches.

Correlation analysis of interstitial maturity and prognosis of colorectal cancer: Meta-analysis

BACKGROUND

To investigate the relationship between interstitial maturity and prognosis of colorectal cancer.

AIM

To examine the correlation between interstitial maturity and the prognosis of colorectal cancer.

METHODS

The paper database PubMed, EMBASE, Cochranelibrary, Springerlink, CNKI, and Wanfang database were searched until December 2023. "tumor stroma maturity" "desmoplastic stroma reaction" "desmoplastic reaction" "stroma reaction" "degree of stroma reaction "" stroma classification" "stroma density" "colorectal cancer" "colon cancer" "rectal cancer" "prognosis" were searched for the search terms. Two system assessors independently screened the literature quality according to the inclusion exclusion criteria, Quality evaluation and data extraction were performed for the included literatures, and meta-analysis was performed for randomized control trials included at using Review Manager 5.2 software.

RESULTS

Finally, data of 9849 patients with colorectal cancer from 19 cosets in 15 literatures were included, including 4339 patients with mature type (control group), 3048 patients with intermediate type (intermediate group) and 2456 patients with immature type (immature group). The results of meta-analysis showed: Relapse-free survival [hazard ratio (HR) = 2.66, 95% confidence interval (CI): 2.30-3.08; P < 0.00001], disease-free survival (HR = 3.68, 95%CI: 2.33-5.81; P < 0.00001) and overall survival (HR = 1.70, 95%CI: 1.53-1.87; P < 0.00001) were significantly lower than those in mature group (control group); relapse-free survival (HR = 1.36, 95%CI: 1.17-1.59; P < 0.0001) and disease-free survival rate (HR = 1.85, 95%CI: 1.53-2.24; P < 0.0001) was significantly lower than the mature group (control group).

CONCLUSION

There is the correlation between tumor interstitial maturity and survival prognosis of colorectal cancer, and different degrees of tumor interstitial maturity have a certain impact on the quality of life of colorectal cancer patients.

Current understanding on TREM-2 molecular biology and physiopathological functions

Triggering receptor expressed on myeloid cells 2 (TREM-2), a glycosylated receptor belonging to the immunoglobin superfamily and especially expressed in the myeloid cell lineage, is frequently explained as a reminiscent receptor for both adaptive and innate immunity regulation. TREM-2 is also acknowledged to influence NK cell differentiation via the PI3K and PLCγ signaling pathways, as well as the partial activation or direct inhibition of T cells. Additionally, TREM-2 overexpression is substantially linked to cell-specific functions, such as enhanced phagocytosis, reduced toll-like receptor (TLR)-mediated inflammatory cytokine production, increased transcription of anti-inflammatory cytokines, and reshaped T cell function. Whereas TREM-2-deficient cells exhibit diminished phagocytic function and enhanced proinflammatory cytokines production, proceeding to inflammatory injuries and an immunosuppressive environment for disease progression. Despite the growing literature supporting TREM-2+ cells in various diseases, such as neurodegenerative disorders and cancer, substantial facets of TREM-2-mediated signaling remain inadequately understood relevant to pathophysiology conditions. In this direction, herein, we have summarized the current knowledge on TREM-2 biology and cell-specific TREM-2 expression, particularly in the modulation of pivotal TREM-2-dependent functions under physiopathological conditions. Furthermore, molecular regulation and generic biological relevance of TREM-2 are also discussed, which might provide an alternative approach for preventing or reducing TREM-2-associated deformities. At last, we discussed the TREM-2 function in supporting an immunosuppressive cancer environment and as a potential drug target for cancer immunotherapy. Hence, summarized knowledge of TREM-2 might provide a window to overcome challenges in clinically effective therapies for TREM-2-induced diseases in humans.

Identification of JUN gene and cellular microenvironment in response to PD-1 blockade treatment in lung cancer patients via single-cell RNA sequencing

Exploring the molecular mechanisms of PD-1/PDL-1 blockade for non-small cell lung cancer (NSCLC) would facilitate understanding for tumor microenvironment (TME) and development of individualized medicine. To date, biomarkers of response to PD-1 blockade therapy were still limited. In this study, we hypothesize that cell type in the tumor microenvironment can influence the effect of PD-1 blockade immunotherapy through specific genes. Therefore, we re-analyze the single-cell RNA sequencing data and validation in tissue from lung adenocarcinoma patients. Dynamic changes of cellular subpopulation were observed after anti-PD-1 immunotherapy among TMEs between primary/metastasis or good/poor response patients. Non-exhausted CD8 T cells and dysregulated genes were observed in responsing patients from PD-1 blockade therapy. Among all changed genes, JUN, involved in PD-1 blockade immunotherapy pathway, and could be considered as a PD-1 responsing biomarker.

Ginsenosides: an immunomodulator for the treatment of colorectal cancer

Ginsenosides, the primary bioactive ingredients derived from the root of Panax ginseng, are eagerly in demand for tumor patients as a complementary and alternative drug. Ginsenosides have increasingly become a "hot topic" in recent years due to their multifunctional role in treating colorectal cancer (CRC) and regulating tumor microenvironment (TME). Emerging experimental research on ginsenosides in the treatment and immune regulation of CRC has been published, while no review sums up its specific role in the CRC microenvironment. Therefore, this paper systematically introduces how ginsenosides affect the TME, specifically by enhancing immune response, inhibiting the activation of stromal cells, and altering the hallmarks of CRC cells. In addition, we discuss their impact on the physicochemical properties of the tumor microenvironment. Furthermore, we discuss the application of ginsenosides in clinical treatment as their efficacy in enhancing tumor patient immunity and prolonging survival. The future perspectives of ginsenoside as a complementary and alternative drug of CRC are also provided. This review hopes to open up a new horizon for the cancer treatment of Traditional Chinese Medicine monomers.

Exploring immune status in peripheral blood and tumor tissue in association with survival in patients with multi-organ metastatic colorectal cancer

Colorectal cancer (CRC) raises considerable clinical challenges, including a high mortality rate once the tumor spreads to distant sites. At this advanced stage, more accurate prediction of prognosis and treatment outcome is urgently needed. The role of cancer immunity in metastatic CRC (mCRC) is poorly understood. Here, we explore cellular immune cell status in patients with multi-organ mCRC. We analyzed T cell infiltration in primary tumor sections, surveyed the lymphocytic landscape of liver metastases, and assessed circulating mononuclear immune cells. Besides asking whether immune cells are associated with survival at this stage of the disease, we investigated correlations between the different tissue types; as this could indicate a dominant immune phenotype. Taken together, our analyses corroborate previous observations that higher levels of CD8+ T lymphocytes link to better survival outcomes. Our findings therefore extend evidence from earlier stages of CRC to indicate an important role for cancer immunity in disease control even after metastatic spreading to multiple organs. This finding may help to improve predicting outcome of patients with mCRC and suggests a future role for immunotherapeutic strategies.

Quantitative Proteomic Analysis Reveals Functional Alterations of the Peripheral Immune System in Colorectal Cancer

Colorectal cancer (CRC) is characterized by high morbidity, high mortality, and limited response to immunotherapies. The peripheral immune system is an important component of tumor immunity, and enhancements of peripheral immunity help to suppress tumor progression. However, the functional alterations of the peripheral immune system in CRC are unclear. Here, we used mass spectrometry-based quantitative proteomics to establish a protein expression atlas for the peripheral immune system in CRC, including plasma and five types of immune cells (CD4+ T cells, CD8+ T cells, monocytes, natural killer cells, and B cells). Synthesizing the results of the multidimensional analysis, we observed an enhanced inflammatory phenotype in CRC, including elevated expression of plasma inflammatory proteins, activation of the inflammatory pathway in monocytes, and increased inflammation-related ligand-receptor interactions. Notably, we observed tumor effects on peripheral T cells, including altered cell subpopulation ratios and suppression of cell function. Suppression of CD4+ T cell function is mainly mediated by high expression levels of protein tyrosine phosphatases. Among them, the expression of protein tyrosine phosphatase receptor type J (PTPRJ) gradually increased with CRC progression; knockdown of PTPRJ in vitro could promote T cell activation, thereby enhancing peripheral immunity. We also found that the combination of leucine-rich α-2 glycoprotein 1 (LRG1) and apolipoprotein A4 (APOA4) had the best predictive ability for colorectal cancer and has the potential to be a biomarker. Overall, this study provides a comprehensive understanding of the peripheral immune system in CRC. It also offers insights regarding the potential clinical utilities of these peripheral immune characteristics as diagnostic indicators and therapeutic targets.

Mechanism of bufalin inhibition of colon cancer liver metastasis by regulating M2-type polarization of Kupffer cells induced by highly metastatic colon cancer cells

Patients with metastatic colorectal cancer often have poor outcomes, primarily due to hepatic metastasis. Colorectal cancer (CRC) cells have the ability to secrete cytokines and other molecules that can remodel the tumor microenvironment, facilitating the spread of cancer to the liver. Kupffer cells (KCs), which are macrophages in the liver, can be polarized to M2 type, thereby promoting the expression of adhesion molecules that aid in tumor metastasis. Our research has shown that huachanshu (with bufalin as the main active monomer) can effectively inhibit CRC metastasis. However, the underlying mechanism still needs to be thoroughly investigated. We have observed that highly metastatic CRC cells have a greater ability to induce M2-type polarization of Kupffer cells, leading to enhanced metastasis. Interestingly, we have found that inhibiting the expression of IL-6, which is highly expressed in the serum, can reverse this phenomenon. Notably, bufalin has been shown to attenuate the M2-type polarization of Kupffer cells induced by highly metastatic Colorectal cancer (mCRC) cells and down-regulate IL-6 expression, ultimately inhibiting tumor metastasis. In this project, our aim is to study how high mCRC cells induce M2-type polarization and how bufalin, via the SRC-3/IL-6 pathway, can inhibit CRC metastasis. This research will provide a theoretical foundation for understanding the anti-CRC effect of bufalin.

Pan-Cancer Single-Cell and Spatial-Resolved Profiling Reveals the Immunosuppressive Role of APOE+ Macrophages in Immune Checkpoint Inhibitor Therapy

The heterogeneity of macrophages influences the response to immune checkpoint inhibitor (ICI) therapy. However, few studies explore the impact of APOE+ macrophages on ICI therapy using single-cell RNA sequencing (scRNA-seq) and machine learning methods. The scRNA-seq and bulk RNA-seq data are Integrated to construct an M.Sig model for predicting ICI response based on the distinct molecular signatures of macrophage and machine learning algorithms. Comprehensive single-cell analysis as well as in vivo and in vitro experiments are applied to explore the potential mechanisms of the APOE+ macrophage in affecting ICI response. The M.Sig model shows clear advantages in predicting the efficacy and prognosis of ICI therapy in pan-cancer patients. The proportion of APOE+ macrophages is higher in ICI non-responders of triple-negative breast cancer compared with responders, and the interaction and longer distance between APOE+ macrophages and CD8+ exhausted T (Tex) cells affecting ICI response is confirmed by multiplex immunohistochemistry. In a 4T1 tumor-bearing mice model, the APOE inhibitor combined with ICI treatment shows the best efficacy. The M.Sig model using real-world immunotherapy data accurately predicts the ICI response of pan-cancer, which may be associated with the interaction between APOE+ macrophages and CD8+ Tex cells.

HOXB9 promotes osteosarcoma cell survival and malignancy under glucose starvation via upregulating SPP1 expression

Homeobox B9 (HOXB9) has been shown to play a critical role in several tumors. However, the precise biological mechanisms and functions of HOXB9 in osteosarcoma remain largely unknown. In this study, we found that HOXB9 was increased upon glucose starvation. Elevated HOXB9 suppressed osteosarcoma cell death and supported cell growth and migration under glucose starvation. Further mechanistic studies demonstrated that HOXB9 directly bound to the promoter of secreted phosphoprotein 1 (SPP1) and transcriptionally upregulated SPP1 expression which then led cell death decrease and cell growth increase under glucose deprivation environment. Clinically, HOXB9 was significantly upregulated in osteosarcoma compared with normal tissues and increase of HOXB9 expression was positively associated with the elevation of SPP1 in osteosarcoma. Overall, our study illustrates that HOXB9 contributes to malignancy in osteosarcoma and inhibits cell death through transcriptional upregulating SPP1 under glucose starvation.

Multi-model analysis of gallbladder cancer reveals the role of OxLDL-absorbing neutrophils in promoting liver invasion

BACKGROUND

Gallbladder cancer (GBC) is the most common and lethal malignancy of the biliary tract that lacks effective therapy. In many GBC cases, infiltration into adjacent organs or distant metastasis happened long before the diagnosis, especially the direct liver invasion, which is the most common and unfavorable way of spreading.

METHODS

Single-cell RNA sequencing (scRNA-seq), spatial transcriptomics (ST), proteomics, and multiplexed immunohistochemistry (mIHC) were performed on GBC across multiple tumor stages to characterize the tumor microenvironment (TME), focusing specifically on the preferential enrichment of neutrophils in GBC liver invasion (GBC-LI).

RESULTS

Multi-model Analysis reveals the immunosuppressive TME of GBC-LI that was characterized by the enrichment of neutrophils at the invasive front. We identified the context-dependent transcriptional states of neutrophils, with the Tumor-Modifying state being associated with oxidized low-density lipoprotein (oxLDL) metabolism. In vitro assays showed that the direct cell-cell contact between GBC cells and neutrophils led to the drastic increase in oxLDL uptake of neutrophils, which was primarily mediated by the elevated OLR1 on neutrophils. The oxLDL-absorbing neutrophils displayed a higher potential to promote tumor invasion while demonstrating lower cancer cytotoxicity. Finally, we identified a neutrophil-promoting niche at the invasive front of GBC-LI that constituted of KRT17+ GBC cells, neutrophils, and surrounding fibroblasts, which may help cultivate the oxLDL-absorbing neutrophils.

CONCLUSIONS

Our study reveals the existence of a subset of pro-tumoral neutrophils with a unique ability to absorb oxLDL via OLR1, a phenomenon induced through cell-cell contact with KRT17+ GBC cells in GBC-LI.

Interactions between LAMP3+ dendritic cells and T-cell subpopulations promote immune evasion in papillary thyroid carcinoma

BACKGROUND

The incidence of papillary thyroid cancer (PTC) continues to rise all over the world, 10-15% of the patients have a poor prognosis. Although immunotherapy has been applied in clinical practice, its therapeutic efficacy remains far from satisfactory, necessitating further investigation of the mechanism of PTC immune remodeling and exploration of novel treatment targets.

METHODS

This study conducted a single-cell RNA sequencing (scRNA-seq) analysis using 18 surgical tissue specimens procured from 14 patients diagnosed with adjacent tissues, non-progressive PTC or progressive PTC. Key findings were authenticated through spatial transcriptomics RNA sequencing, immunohistochemistry, multiplex immunohistochemistry, and an independent bulk RNA-seq data set containing 502 samples.

RESULTS

A total of 151,238 individual cells derived from 18 adjacent tissues, non-progressive PTC and progressive PTC specimens underwent scRNA-seq analysis. We found that progressive PTC exhibits the following characteristics: a significant decrease in overall immune cells, enhanced immune evasion of tumor cells, and disrupted antigen presentation function. Moreover, we identified a subpopulation of lysosomal associated membrane protein 3 (LAMP3+) dendritic cells (DCs) exhibiting heightened infiltration in progressive PTC and associated with advanced T stage and poor prognosis of PTC. LAMP3+ DCs promote CD8+ T cells exhaustion (mediated by NECTIN2-TIGIT) and increase infiltration abundance of regulatory T cells (mediated by chemokine (C-C motif) ligand 17 (CCL17)-chemokine (C-C motif) receptor 4 (CCR4)) establishing an immune-suppressive microenvironment. Ultimately, we unveiled that progressive PTC tumor cells facilitate the retention of LAMP3+ DCs within the tumor microenvironment through NECTIN3-NECTIN2 interactions, thereby rendering tumor cells more susceptible to immune evasion.

CONCLUSION

Our findings expound valuable insights into the role of the interaction between LAMP3+ DCs and T-cell subpopulations and offer new and effective ideas and strategies for immunotherapy in patients with progressive PTC.

Integration of single-cell sequencing and bulk RNA-seq to identify and develop a prognostic signature related to colorectal cancer stem cells

The prognosis for patients with colorectal cancer (CRC) remains worse than expected due to metastasis, recurrence, and resistance to chemotherapy. Colorectal cancer stem cells (CRCSCs) play a vital role in tumor metastasis, recurrence, and chemotherapy resistance. However, there are currently no prognostic markers based on CRCSCs-related genes available for clinical use. In this study, single-cell transcriptome sequencing was employed to distinguish cancer stem cells (CSCs) in the CRC microenvironment and analyze their properties at the single-cell level. Subsequently, data from TCGA and GEO databases were utilized to develop a prognostic risk model for CRCSCs-related genes and validate its diagnostic performance. Additionally, functional enrichment, immune response, and chemotherapeutic drug sensitivity of the relevant genes in the risk model were investigated. Lastly, the key gene RPS17 in the risk model was identified as a potential prognostic marker and therapeutic target for further comprehensive studies. Our findings provide new insights into the prognostic treatment of CRC and offer novel perspectives for a systematic and comprehensive understanding of CRC development.

New Insight into Intestinal Mast Cells Revealed by Single-Cell RNA Sequencing

Mast cells (MCs) are tissue-resident immune cells distributed in all tissues and strategically located close to blood and lymphatic vessels and nerves. Thanks to the expression of a wide array of receptors, MCs act as tissue sentinels, able to detect the presence of bacteria and parasites and to respond to different environmental stimuli. MCs originate from bone marrow (BM) progenitors that enter the circulation and mature in peripheral organs under the influence of microenvironment factors, thus differentiating into heterogeneous tissue-specific subsets. Even though MC activation has been traditionally linked to IgE-mediated allergic reactions, a role for these cells in other pathological conditions including tumor progression has recently emerged. However, several aspects of MC biology remain to be clarified. The advent of single-cell RNA sequencing platforms has provided the opportunity to understand MCs' origin and differentiation as well as their phenotype and functions within different tissues, including the gut. This review recapitulates how single-cell transcriptomic studies provided insight into MC development as well as into the functional role of intestinal MC subsets in health and disease.

Mechanisms of myeloid-derived suppressor cell-mediated immunosuppression in colorectal cancer and related therapies

Severe immunosuppression is a hallmark of colorectal cancer (CRC). Myeloid-derived suppressor cells (MDSCs), one of the most abundant components of the tumor stroma, play an important role in the invasion, metastasis, and immune escape of CRC. MDSCs create an immunosuppressive microenvironment by inhibiting the proliferation and activation of immunoreactive cells, including T and natural killer cells, as well as by inducing the proliferation of immunosuppressive cells, such as regulatory T cells and tumor-associated macrophages, which, in turn, promote the growth of cancer cells. Thus, MDSCs are key contributors to the emergence of an immunosuppressive microenvironment in CRC and play an important role in the breakdown of antitumor immunity. In this narrative review, we explore the mechanisms through which MDSCs contribute to the immunosuppressive microenvironment, the current therapeutic approaches and technologies targeting MDSCs, and the therapeutic potential of modulating MDSCs in CRC treatment. This study provides ideas and methods to enhance survival rates in patients with CRC.

IL-2RG as a possible immunotherapeutic target in CRC predicting poor prognosis and regulated by miR-7-5p and miR-26b-5p

Despite advances in treatment strategies, colorectal cancer (CRC) continues to cause significant morbidity and mortality, with mounting evidence a close link between immune system dysfunctions issued. Interleukin-2 receptor gamma (IL-2RG) plays a pivotal role as a common subunit receptor in the IL-2 family cytokines and activates the JAK-STAT pathway. This study delves into the role of Interleukin-2 receptor gamma (IL-2RG) within the tumor microenvironment and investigates potential microRNAs (miRNAs) that directly inhibit IL-2RG, aiming to discern their impact on CRC clinical outcomes. Bioinformatics analysis revealed a significant upregulation of IL-2RG mRNA in TCGA-COAD samples and showed strong correlations with the infiltration of various lymphocytes. Single-cell analysis corroborated these findings, highlighting IL-2RG expression in critical immune cell subsets. To explore miRNA involvement in IL-2RG dysregulation, mRNA was isolated from the tumor tissues and lymphocytes of 258 CRC patients and 30 healthy controls, and IL-2RG was cloned into the pcDNA3.1/CT-GFP-TOPO vector. Human embryonic kidney cell lines (HEK-293T) were transfected with this construct. Our research involved a comprehensive analysis of miRPathDB, miRWalk, and Targetscan databases to identify the miRNAs associated with the 3' UTR of human IL-2RG. The human microRNA (miRNA) molecules, hsa-miR-7-5p and hsa-miR-26b-5p, have been identified as potent suppressors of IL-2RG expression in CRC patients. Specifically, the downregulation of hsa-miR-7-5p and hsa-miR-26b-5p has been shown to result in the upregulation of IL-2RG mRNA expression in these patients. Prognostic evaluation of IL-2RG, hsa-miR-7-5p, and hsa-miR-26b-5p, using TCGA-COAD data and patient samples, established that higher IL-2RG expression and lower expression of both miRNAs were associated with poorer outcomes. Additionally, this study identified several long non-coding RNAs (LncRNAs), such as ZFAS1, SOX21-AS1, SNHG11, SNHG16, SNHG1, DLX6-AS1, GAS5, SNHG6, and MALAT1, which may act as competing endogenous RNA molecules for IL2RG by sequestering shared hsa-miR-7-5p and hsa-miR-26b-5p. In summary, this investigation underscores the potential utility of IL-2RG, hsa-miR-7-5p, and hsa-miR-26b-5p as serum and tissue biomarkers for predicting CRC patient prognosis while also offering promise as targets for immunotherapy in CRC management.

Early Immune Remodeling Steers Clinical Response to First-Line Chemoimmunotherapy in Advanced Gastric Cancer

Adding anti-programmed cell death protein 1 (anti-PD-1) to 5-fluorouracil (5-FU)/platinum improves survival in some advanced gastroesophageal adenocarcinomas (GEA). To understand the effects of chemotherapy and immunotherapy, we conducted a phase II first-line trial (n = 47) sequentially adding pembrolizumab to 5-FU/platinum in advanced GEA. Using serial biopsy of the primary tumor at baseline, after one cycle of 5-FU/platinum, and after the addition of pembrolizumab, we transcriptionally profiled 358,067 single cells to identify evolving multicellular tumor microenvironment (TME) networks. Chemotherapy induced early on-treatment multicellular hubs with tumor-reactive T-cell and M1-like macrophage interactions in slow progressors. Faster progression featured increased MUC5A and MSLN containing treatment resistance programs in tumor cells and M2-like macrophages with immunosuppressive stromal interactions. After pembrolizumab, we observed increased CD8 T-cell infiltration and development of an immunity hub involving tumor-reactive CXCL13 T-cell program and epithelial interferon-stimulated gene programs. Strategies to drive increases in antitumor immune hub formation could expand the portion of patients benefiting from anti-PD-1 approaches.

SIGNIFICANCE

The benefit of 5-FU/platinum with anti-PD-1 in first-line advanced gastric cancer is limited to patient subgroups. Using a trial with sequential anti-PD-1, we show coordinated induction of multicellular TME hubs informs the ability of anti-PD-1 to potentiate T cell-driven responses. Differential TME hub development highlights features that underlie clinical outcomes. This article is featured in Selected Articles from This Issue, p. 695.

Association of preoperative aspartate aminotransferase to platelet ratio index with outcomes and tumour microenvironment among colorectal cancer with liver metastases

This study aims to investigate applicable robust biomarkers that can improve prognostic predictions for colorectal liver metastasis (CRLM) patients receiving simultaneous resection. A total of 1323 CRLM patients from multiple centres were included. The preoperative aspartate aminotransferase to platelet ratio index (APRI) level from blood of patients were obtained. Patients were stratified into a high APRI group and a low APRI group, and comparisons were conducted by analyzing progression-free survival (PFS), overall survival (OS) and postoperative early recurrence. Tumour samples of CRLM were collected to perform single-cell RNA sequencing and multiplex immunohistochemistry/immunofluorescence (mIHC/IF) to investigate the association of APRI levels and the tumour microenvironment of CRLM. Compared with APRI <0.33, PFS disadvantage (IPTW-adjusted HR = 1.240, P = 0.015) and OS disadvantage (IPTW- adjusted HR = 1.507, P = 0.002) of APRI ≥0.33 were preserved in the IPTW-adjusted Cox hazards regression analyses. An APRI ≥0.25 was associated with a significantly increased risk of postoperative early recurrence after adjustment (IPTW-adjusted OR = 1.486, P = 0.001). The external validation showed consistent results with the training cohort. In the high APRI group, the single-cell RNA sequencing results revealed a heightened malignancy of epithelial cells, the enrichment of inflammatory-like cancer-associated fibroblasts and SPP1+ macrophages associated with activation of malignant cells and fibrotic microenvironment, and a more suppressed-function T cells; mIHC/IF showed that PD1+ CD4+ T cells, FOXP3+ CD4+ T cells, PD1+ CD8+ T cells, FOXP3+ CD8+ T cells, SPP1+ macrophages and iCAFs were significantly increased in the intratumoral region and peritumoral region. This study contributed valuable evidence regarding preoperative APRI for predicting prognoses among CRLM patients receiving simultaneous resection and provided underlying clues supporting the association between APRI and clinical outcomes by single-cell sequencing bioinformatics analysis and mIHC/IF.

Identification of BGN positive fibroblasts as a driving factor for colorectal cancer and development of its related prognostic model combined with machine learning

BACKGROUND

Numerous studies have indicated that cancer-associated fibroblasts (CAFs) play a crucial role in the progression of colorectal cancer (CRC). However, there are still many unknowns regarding the exact role of CAF subtypes in CRC.

METHODS

The data for this study were obtained from bulk, single-cell, and spatial transcriptomic sequencing data. Bioinformatics analysis, in vitro experiments, and machine learning methods were employed to investigate the functional characteristics of CAF subtypes and construct prognostic models.

RESULTS

Our study demonstrates that Biglycan (BGN) positive cancer-associated fibroblasts (BGN + Fib) serve as a driver in colorectal cancer (CRC). The proportion of BGN + Fib increases gradually with the progression of CRC, and high infiltration of BGN + Fib is associated with poor prognosis in terms of overall survival (OS) and recurrence-free survival (RFS) in CRC. Downregulation of BGN expression in cancer-associated fibroblasts (CAFs) significantly reduces migration and proliferation of CRC cells. Among 101 combinations of 10 machine learning algorithms, the StepCox[both] + plsRcox combination was utilized to develop a BGN + Fib derived risk signature (BGNFRS). BGNFRS was identified as an independent adverse prognostic factor for CRC OS and RFS, outperforming 92 previously published risk signatures. A Nomogram model constructed based on BGNFRS and clinical-pathological features proved to be a valuable tool for predicting CRC prognosis.

CONCLUSION

In summary, our study identified BGN + Fib as drivers of CRC, and the derived BGNFRS was effective in predicting the OS and RFS of CRC patients.

Unravelling immune microenvironment features underlying tumor progression in the single-cell era

The relationship between the immune cell and tumor occurrence and progression remains unclear. Profiling alterations in the tumor immune microenvironment (TIME) at high resolution is crucial to identify factors influencing cancer progression and enhance the effectiveness of immunotherapy. However, traditional sequencing methods, including bulk RNA sequencing, exhibit varying degrees of masking the cellular heterogeneity and immunophenotypic changes observed in early and late-stage tumors. Single-cell RNA sequencing (scRNA-seq) has provided significant and precise TIME landscapes. Consequently, this review has highlighted TIME cellular and molecular changes in tumorigenesis and progression elucidated through recent scRNA-seq studies. Specifically, we have summarized the cellular heterogeneity of TIME at different stages, including early, late, and metastatic stages. Moreover, we have outlined the related variations that may promote tumor occurrence and metastasis in the single-cell era. The widespread applications of scRNA-seq in TIME will comprehensively redefine the understanding of tumor biology and furnish more effective immunotherapy strategies.

Investigating the Cell Origin and Liver Metastasis Factors of Colorectal Cancer by Single-Cell Transcriptome Analysis

Background

Colorectal cancer (CRC) is one of the deadliest causes of death by cancer worldwide. Liver metastasis (LM) is the main cause of death in patients with CRC. Therefore, identification of patients with the greatest risk of liver metastasis is critical for early treatment and reduces the mortality of patients with colorectal cancer liver metastases.

Methods

Initially, we characterized cell composition through single-cell transcriptome analysis. Subsequently, we employed copy number variation (CNV) and pseudotime analysis to delineate the cellular origins of LM and identify LM-related epithelial cells (LMECs). The LM-index was constructed using machine learning algorithms to forecast the relative abundance of LMECs, reflecting the risk of LM. Furthermore, we analyzed drug sensitivity and drug targeted gene expression in LMECs and patients with a high risk of LM. Finally, functional experiments were conducted to determine the biological roles of metastasis-related gene in vitro.

Results

Single-cell RNA sequencing analysis revealed different immune landscapes between primary CRC and LM tumor. LM originated from chromosomal variants with copy number loss of chr1 and chr6p and copy number gain of chr7 and chr20q. We identified the LMECs cluster and found LM-associated pathways such as Wnt/beta-catenin signaling and KRAS signaling. Subsequently, we identified ten metastasis-associated genes, including SOX4, and established the LM-index, which correlates with poorer prognosis, higher stage, and advanced age. Furthermore, we screened two drugs as potential candidates for treating LM, including Linsitinib_1510, Lapatinib_1558. Immunohistochemistry results demonstrated significantly elevated SOX4 expression in tumor samples compared to normal samples. Finally, in vitro experiments verified that silencing SOX4 significantly inhibited tumor cell migration and invasion.

Conclusion

This study reveals the possible cellular origin and driving factors of LM in CRC at the single cell level, and provides a reference for early detection of CRC patients with a high risk of LM.

Th17 Cells Secrete TWEAK to Trigger Epithelial-Mesenchymal Transition and Promote Colorectal Cancer Liver Metastasis

Liver metastasis is the leading cause of mortality in patients with colorectal cancer. Given the significance of both epithelial-mesenchymal transition (EMT) of tumor cells and the immune microenvironment in colorectal cancer liver metastasis (CRLM), the interplay between them could hold the key for developing improved treatment options. We employed multiomics analysis of 130 samples from 18 patients with synchronous CRLM integrated with external datasets to comprehensively evaluate the interaction between immune cells and EMT of tumor cells in liver metastasis. Single-cell RNA sequencing analysis revealed distinct distributions of nonmalignant cells between primary tumors from patients with metastatic colorectal cancer (mCRC) and non-metastatic colorectal cancer, showing that Th17 cells were predominantly enriched in the primary lesion of mCRC. TWEAK, a cytokine secreted by Th17 cells, promoted EMT by binding to receptor Fn14 on tumor cells, and the TWEAK-Fn14 interaction enhanced tumor migration and invasion. In mouse models, targeting Fn14 using CRISPR-induced knockout or lipid nanoparticle-encapsulated siRNA alleviated metastasis and prolonged survival. Mice lacking Il17a or Tnfsf12 (encoding TWEAK) exhibited fewer metastases compared with wild-type mice, while cotransfer of Th17 with tumor cells promoted liver metastasis. Higher TWEAK expression was associated with a worse prognosis in patients with colorectal cancer. In addition, CD163L1+ macrophages interacted with Th17 cells, recruiting Th17 via the CCL4-CCR5 axis. Collectively, this study unveils the role of immune cells in the EMT process and identifies TWEAK secreted by Th17 as a driver of CRLM.

SIGNIFICANCE

TWEAK secreted by Th17 cells promotes EMT by binding to Fn14 on colorectal cancer cells, suggesting that blocking the TWEAK-Fn14 interaction may be a promising therapeutic approach to inhibit liver metastasis.

IL-10 dampens antitumor immunity and promotes liver metastasis via PD-L1 induction

BACKGROUND & AIMS

The liver is one of the organs most commonly affected by metastasis. The presence of liver metastases has been reported to be responsible for an immunosuppressive microenvironment and diminished immunotherapy efficacy. Herein, we aimed to investigate the role of IL-10 in liver metastasis and to determine how its modulation could affect the efficacy of immunotherapy in vivo.

METHODS

To induce spontaneous or forced liver metastasis in mice, murine cancer cells (MC38) or colon tumor organoids were injected into the cecum or the spleen, respectively. Mice with complete and cell type-specific deletion of IL-10 and IL-10 receptor alpha were used to identify the source and the target of IL-10 during metastasis formation. Programmed death ligand 1 (PD-L1)-deficient mice were used to test the role of this checkpoint. Flow cytometry was applied to characterize the regulation of PD-L1 by IL-10.

RESULTS

We found that Il10-deficient mice and mice treated with IL-10 receptor alpha antibodies were protected against liver metastasis formation. Furthermore, by using IL-10 reporter mice, we demonstrated that Foxp3+ regulatory T cells (Tregs) were the major cellular source of IL-10 in liver metastatic sites. Accordingly, deletion of IL-10 in Tregs, but not in myeloid cells, led to reduced liver metastasis. Mechanistically, IL-10 acted on Tregs in an autocrine manner, thereby further amplifying IL-10 production. Furthermore, IL-10 acted on myeloid cells, i.e. monocytes, and induced the upregulation of the immune checkpoint protein PD-L1. Finally, the PD-L1/PD-1 axis attenuated CD8-dependent cytotoxicity against metastatic lesions.

CONCLUSIONS

Treg-derived IL-10 upregulates PD-L1 expression in monocytes, which in turn reduces CD8+ T-cell infiltration and related antitumor immunity in the context of colorectal cancer-derived liver metastases. These findings provide the basis for future monitoring and targeting of IL-10 in colorectal cancer-derived liver metastases.

IMPACT AND IMPLICATIONS

Liver metastasis diminishes the effectiveness of immunotherapy and increases the mortality rate in patients with colorectal cancer. We investigated the role of IL-10 in liver metastasis formation and assessed its impact on the effectiveness of immunotherapy. Our data show that IL-10 is a pro-metastatic factor involved in liver metastasis formation and that it acts as a regulator of PD-L1. This provides the basis for future monitoring and targeting of IL-10 in colorectal cancer-derived liver metastasis.

Comprehensive investigation into the influence of glycosylation on head and neck squamous cell carcinoma and development of a prognostic model for risk assessment and anticipating immunotherapy

Background

It has been well established that glycosylation plays a pivotal role in initiation, progression, and therapy resistance of several cancers. However, the correlations between glycosylation and head and neck squamous cell carcinoma (HNSCC) have not been elucidated in detail.

Methods

The paramount genes governing glycosylation were discerned via the utilization of the Protein-Protein Interaction (PPI) network and correlation analysis, coupled with single-cell RNA sequencing (scRNA-seq) analysis. To construct risk models exhibiting heightened predictive efficacy, cox- and lasso-regression methodologies were employed, and the veracity of these models was substantiated across both internal and external datasets. Subsequently, an exploration into the distinctions within the tumor microenvironment (TME), immunotherapy responses, and enriched pathways among disparate risk cohorts ensued. Ultimately, cell experiments were conducted to validate the consequential impact of SMS in Head and Neck Squamous Cell Carcinoma (HNSCC).

Results

A total of 184 genes orchestrating glycosylation were delineated for subsequent scrutiny. Employing cox- and lasso-regression methodologies, we fashioned a 3-gene signature, proficient in prognosticating the outcomes for patients afflicted with HNSCC. Noteworthy observations encompassed distinctions in the Tumor Microenvironment (TME), levels of immune cell infiltration, and the presence of immune checkpoint markers among divergent risk cohorts, holding potentially consequential implications for the clinical management of HNSCC patients.

Conclusion

The prognosis of HNSCC can be proficiently anticipated through risk signatures based on Glycosylation-related genes (GRGs). A thorough delineation of the GRGs signature in HNSCC holds the potential to facilitate the interpretation of HNSCC's responsiveness to immunotherapy and provide innovative strategies for cancer treatment.

Circulation immune cell landscape in canonical pathogenesis of colorectal adenocarcinoma by CyTOF analysis

Current studies on the immune microenvironment of colorectal cancer (CRC) were mostly limited to the tissue level, lacking relevant studies in the peripheral blood, and failed to describe its alterations in the whole process of adenocarcinoma formation, especially of adenoma carcinogenesis. Here, we constructed a large-scale population cohort and used the CyTOF to explore the changes of various immune cell subsets in peripheral blood of CRC. We found monocytes and basophils cells were significantly higher in adenocarcinoma patients. Compared with early-stage CRC, effector CD4+T cells and naive B cells were higher in patients with lymph node metastasis, whereas the basophils were lower. We also performed random forest algorithm and found monocytes play the key role in carcinogenesis. Our study draws a peripheral blood immune cell landscape of the occurrence and development of CRC at the single-cell level and provides a reference for other researchers.

Identification and prognostic analysis of candidate biomarkers for lung metastasis in colorectal cancer

Colorectal cancer (CRC) is one of the most prevalent types of malignant tumors. It's vital to explore new biomarkers and potential therapeutic targets in CRC lung metastasis through adopting integrated bioinformatics tools. Multiple cohort datasets and databases were integrated to clarify and verify potential key candidate biomarkers and signal transduction pathways in CRC lung metastasis. DAVID, STRING, UALCAN, GEPIA, TIMER, cBioPortal, THE HUMAN PROTEIN ATLAS, GSEA 4.3.2, FUNRICH 3.1.3, and R 4.2.3 were utilized in this study. The enriched biological processes and pathways modulated by the differentially expressed genes (DEGs) were determined with Gene Ontology, Kyoto Encyclopedia of Genes and Genomes. The search tool Retrieval of Interacting Genes and Cytoscape were used to construct a protein-protein interaction network among DEGs. Four hundred fifty-nine colorectal primary cancer and lung metastatic gene expression profiles were screened from 3 gene expression profiles (GSE41258, GSE68468, and GSE41568). Forty-one upregulated genes and 8 downregulated genes were identified from these 3 gene expression profiles and verified by the transcriptional levels of hub genes in other GEO datasets and The Cancer Genome Atlas database. Two pathways (immune responses and chemokine receptors bind chemokines), 13 key DEGs, 6 hub genes (MMP3, SFTPD, ABCA3, CLU, APOE, and SPP1), and 2 biomarkers (APOE, SPP1) with significantly prognostic values were screened. Forty-nine DEGs were identified as potential candidate diagnostic biomarkers for patients with CRC lung metastasis in present study. Enrichment analysis indicated that immune responses and chemokine receptors bind chemokines may play a leading role in lung metastasis of CRC, and further studies are needed to validate these findings.

Precursor exhausted CD8+T cells in colorectal cancer tissues associated with patient's survival and immunotherapy responsiveness

Exhausted CD8+T cells represent a distinct cellular lineage that emerges during both chronic infections and cancers. Recent studies have shown that persistent antigen exposure can drive the differentiation of precursor exhausted CD8+T cells, termed Tpex cells, which are characterized as TCF-1+PD-1+CD8+T cells. Elevated Tpex cell frequencies in the tumor microenvironment (TME) are associated with improved overall survival (OS) in cancer patients and heightened responsiveness to anti-PD-1 therapy. In our present study, we utilized multi-color immunohistochemistry (mIHC) to determine the localization and clinical implications of tumor-infiltrating Tpex cells within the TME of human colorectal cancer (CRC) tissues. We also conducted a multi-omics integrative analysis using single-cell RNA sequencing (scRNA-seq) data derived from both the murine MC38 tumor model and human CRC tissues. This analysis helped delineate the transcriptional and functional attributes of Tpex cells within the CRC TME. Furthermore, we employed spatial transcriptome sequencing data from CRC patients to investigate the interactions between Tpex cells and other immune cell subsets within the TME. In conclusion, our study not only established a method for Tpex cell detection using mIHC technology but also confirmed that assessing Tpex cells within the CRC TME could be indicative of patients' survival. We further uncovered the transcriptional and functional characteristics of Tpex cells in the TME and ascertained their pivotal role in the efficacy of immunotherapy against CRC.

Combination of IL-33 with PD-1 blockade augment mILC2s-mediated anti-tumor immunity

BACKGROUND

Group 2 innate lymphoid cells (ILC2s) represent one of the main tissue-specific innate lymphoid cell populations, which are key drivers of cytokine secretion in their occupational niche. However, the precise involvement of ILC2s in cancer immunity and their potential impact on immunotherapeutic approaches remain poorly understood.

METHODS

The proportion of ILC2s originating from various tissue sources were quantified through flow cytometry, along with the determination of CD4+ T cell and CD8+ T cell percentages. Flow cytometry was also employed to assess IFN-γ production and programmed cell death protein-1 (PD-1) expression in T cells. Immunohistochemistry was utilized to detect IL-33 expression in tumor tissues, while immunofluorescence was employed to confirm the infiltration of ILC2s in both murine and human tumor tissues.

RESULTS

In this study, we provide evidence that intra-tumoral ILC2s in lung adenocarcinoma (LUAD) exist in a quiescent state. However, the activation of intra-tumoral ILC2s is induced by IL-33 specifically in a natural ILC2s (nILC2, ST2+KLRG1-) phenotype. Considering the pivotal role of PD-1 in cancer immunotherapy and its immunoregulatory functions, we investigated the synergistic effects of IL-33 and anti-PD-1 and found that their combination enhances anti-tumor immunity and improves the efficacy of immunotherapy. Moreover, this combination leads to the upregulation of activated mature ILC2s (mILC2, ST2+KLRG1+) phenotype, thereby highlighting the activated ILC2s as a novel enhancer of the immunoregulatory properties of anti-PD-1.

CONCLUSIONS

Collectively, these findings underscore the significance of ILC2s and their contribution to the anti-tumor response in the context of cancer immunotherapy. Consequently, the simultaneous targeting of ILC2s and T cells represents a potentially promising and widely applicable strategy for immunotherapeutic interventions.

scRNA-seq of colorectal cancer shows regional immune atlas with the function of CD20+ B cells

Colorectal cancer (CRC) from different regions exhibits different histological, genetic characteristics, and molecular subtypes, even in response to conventional chemotherapies and immunotherapies. To characterize the immune landscape in different regions of CRC and search for potential therapeutic targets, we analyzed 39,484 single-cell transcription data from 19 samples of CRC and paired normal tissues from four regions to identify the immune characteristics of CRC among anatomic locations, especially in B cells. We discovered that immune cell infiltration in tumors significantly varied among different regions of CRC. B cells from right- and left-sided CRC had different development trajectories, but both had extensive interactions with myeloid cells and T cells. Survival analysis suggested that CD20+ B cells correlated with good prognosis in CRC patients, especially on the right side. Furthermore, the depletion of CD20+ B cells demonstrated that anti-CD20 promoted tumor growth progression and reversed the tumor-killing activity of anti-PD-1 treatment in vivo and in vitro. Our results highlight the characterization of the immune landscape of CRC in different regions. CD20+ B-cell infiltration has been associated with CRC patient prognosis and may promote the tumor-killing role of PD-1 antibodies.

IL-23 stabilizes an effector Treg cell program in the tumor microenvironment

Interleukin-23 (IL-23) is a proinflammatory cytokine mainly produced by myeloid cells that promotes tumor growth in various preclinical cancer models and correlates with adverse outcomes. However, as to how IL-23 fuels tumor growth is unclear. Here, we found tumor-associated macrophages to be the main source of IL-23 in mouse and human tumor microenvironments. Among IL-23-sensing cells, we identified a subset of tumor-infiltrating regulatory T (Treg) cells that display a highly suppressive phenotype across mouse and human tumors. The use of three preclinical models of solid cancer in combination with genetic ablation of Il23r in Treg cells revealed that they are responsible for the tumor-promoting effect of IL-23. Mechanistically, we found that IL-23 sensing represents a crucial signal driving the maintenance and stabilization of effector Treg cells involving the transcription factor Foxp3. Our data support that targeting the IL-23/IL-23R axis in cancer may represent a means of eliciting antitumor immunity.

SIMarker: Cellular similarity detection and its application to diagnosis and prognosis of liver cancer

BACKGROUND

The emergence of single-cell technology offers a unique opportunity to explore cellular similarity and heterogeneity between precancerous diseases and solid tumors. However, there is lacking a systematic study for identifying and characterizing similarities at single-cell resolution.

METHODS

We developed SIMarker, a computational framework to detect cellular similarities between precancerous diseases and solid tumors based on gene expression at single-cell resolution. Taking hepatocellular carcinoma (HCC) as a case study, we quantified the cellular and molecular connections between HCC and cirrhosis. Core analysis modules of SIMarker is publicly available at https://github.com/xmuhuanglab/SIMarker ("SIM" means "similarity" and "Marker" means "biomarkers).

RESULTS

We found PGA5+ hepatocytes in HCC showed cirrhosis-like characteristics, including similar transcriptional programs and gene regulatory networks. Consequently, the genes constituting the gene expression program of these cirrhosis-like subpopulations were designated as cirrhosis-like signatures (CLS). Strikingly, our utilization of CLS enabled the development of diagnosis and prognosis biomarkers based on within-sample relative expression orderings of gene pairs. These biomarkers achieved high precision and concordance compared with previous studies.

CONCLUSIONS

Our work provides a systematic method to investigate the clinical translational significance of cellular similarities between HCC and cirrhosis, which opens avenues for identifying similar paradigms in other categories of cancers and diseases.

Integrative single-cell transcriptomic analyses reveal the cellular ontological and functional heterogeneities of primary and metastatic liver tumors

BACKGROUND

The global cellular landscape of the tumor microenvironment (TME) combining primary and metastatic liver tumors has not been comprehensively characterized.

METHODS

Based on the scRNA-seq and spatial transcriptomic data of non-tumor liver tissues (NTs), primary liver tumors (PTs) and metastatic liver tumors (MTs), we performed the tissue preference, trajectory reconstruction, transcription factor activity inference, cell-cell interaction and cellular deconvolution analyses to construct a comprehensive cellular landscape of liver tumors.

RESULTS

Our analyses depicted the heterogeneous cellular ecosystems in NTs, PTs and MTs. The activated memory B cells and effector T cells were shown to gradually shift to inhibitory B cells, regulatory or exhausted T cells in liver tumors, especially in MTs. Among them, we characterized a unique group of TCF7+ CD8+ memory T cells specifically enriched in MTs that could differentiate into exhausted T cells likely driven by the p38 MAPK signaling. With regard to myeloid cells, the liver-resident macrophages and inflammatory monocyte/macrophages were markedly replaced by tumor-associated macrophages (TAMs), with TREM2+ and UBE2C+ TAMs enriched in PTs, while SPP1+ and WDR45B+ TAMs in MTs. We further showed that the newly identified WDR45B+ TAMs exhibit an M2-like polarization and are associated with adverse prognosis in patients with liver metastases. Additionally, we addressed that endothelial cells display higher immune tolerance and angiogenesis capacity, and provided evidence for the source of the mesenchymal transformation of fibroblasts in tumors. Finally, the malignant hepatocytes and fibroblasts were prioritized as the pivotal cell populations in shaping the microenvironments of PTs and MTs, respectively. Notably, validation analyses by using spatial or bulk transcriptomic data in clinical cohorts concordantly emphasized the clinical significance of these findings.

CONCLUSIONS

This study defines the ontological and functional heterogeneities in cellular ecosystems of primary and metastatic liver tumors, providing a foundation for future investigation of the underlying cellular mechanisms.

Molecular classification and tumor microenvironment characteristics in pheochromocytomas

Pheochromocytomas (PCCs) are rare neuroendocrine tumors that originate from chromaffin cells in the adrenal gland. However, the cellular molecular characteristics and immune microenvironment of PCCs are incompletely understood. Here, we performed single-cell RNA sequencing (scRNA-seq) on 16 tissues from 4 sporadic unclassified PCC patients and 1 hereditary PCC patient with Von Hippel-Lindau (VHL) syndrome. We found that intra-tumoral heterogeneity was less extensive than the inter-individual heterogeneity of PCCs. Further, the unclassified PCC patients were divided into two types, metabolism-type (marked by NDUFA4L2 and COX4I2) and kinase-type (marked by RET and PNMT), validated by immunohistochemical staining. Trajectory analysis of tumor evolution revealed that metabolism-type PCC cells display phenotype of consistently active metabolism and increased metastasis potential, while kinase-type PCC cells showed decreased epinephrine synthesis and neuron-like phenotypes. Cell-cell communication analysis showed activation of the annexin pathway and a strong inflammation reaction in metabolism-type PCCs and activation of FGF signaling in the kinase-type PCC. Although multispectral immunofluorescence staining showed a lack of CD8+ T cell infiltration in both metabolism-type and kinase-type PCCs, only the kinase-type PCC exhibited downregulation of HLA-I molecules that possibly regulated by RET, suggesting the potential of combined therapy with kinase inhibitors and immunotherapy for kinase-type PCCs; in contrast, the application of immunotherapy to metabolism-type PCCs (with antigen presentation ability) is likely unsuitable. Our study presents a single-cell transcriptomics-based molecular classification and microenvironment characterization of PCCs, providing clues for potential therapeutic strategies to treat PCCs.

Single-cell transcriptomic analysis deciphers heterogenous cancer stem-like cells in colorectal cancer and their organ-specific metastasis

OBJECTIVE

Metastasis is the major cause of cancer death. However, what types of heterogenous cancer cells in primary tumour and how they metastasise to the target organs remain largely undiscovered.

DESIGN

We performed single-cell RNA sequencing and spatial transcriptomic analysis in primary colorectal cancer (CRC) and metastases in the liver (lCRC) or ovary (oCRC). We also conducted immunofluorescence staining and functional experiments to examine the mechanism.

RESULTS

Integrative analyses of epithelial cells reveal a stem-like cell cluster with high protein tyrosine phosphatase receptor type O (PTPRO) and achaete scute-like 2 (ASCL2) expression as the metastatic culprit. This cell cluster comprising distinct subpopulations shows distinct liver or ovary metastatic preference. Population 1 (P1) cells with high delta-like ligand 4 (DLL4) and MAF bZIP transcription factor A (MAFA) expression are enriched in primary CRC and oCRC, thus may be associated with ovarian metastasis. P3 cells having a similar expression pattern as cholangiocytes are found mainly in primary CRC and lCRC, presuming to be likely the culprits that specifically metastasise to the liver. Stem-like cells interacted with cancer-associated fibroblasts and endothelial cells via the DLL4-NOTCH signalling pathway to metastasise from primary CRC to the ovary. In the oCRC microenvironment, myofibroblasts provide cancer cells with glutamine and perform a metabolic reprogramming, which may be essential for cancer cells to localise and develop in the ovary.

CONCLUSION

We uncover a mechanism for organ-specific CRC metastasis.

A multiomics analysis-assisted deep learning model identifies a macrophage-oriented module as a potential therapeutic target in colorectal cancer

Colorectal cancer (CRC) is a common malignancy involving multiple cellular components. The CRC tumor microenvironment (TME) has been characterized well at single-cell resolution. However, a spatial interaction map of the CRC TME is still elusive. Here, we integrate multiomics analyses and establish a spatial interaction map to improve the prognosis, prediction, and therapeutic development for CRC. We construct a CRC immune module (CCIM) that comprises FOLR2+ macrophages, exhausted CD8+ T cells, tolerant CD8+ T cells, exhausted CD4+ T cells, and regulatory T cells. Multiplex immunohistochemistry is performed to depict the CCIM. Based on this, we utilize advanced deep learning technology to establish a spatial interaction map and predict chemotherapy response. CCIM-Net is constructed, which demonstrates good predictive performance for chemotherapy response in both the training and testing cohorts. Lastly, targeting FOLR2+ macrophage therapeutics is used to disrupt the immunosuppressive CCIM and enhance the chemotherapy response in vivo.

T-cell infiltration and its regulatory mechanisms in cancers: insights at single-cell resolution

Tumor-infiltrating T cells recognize, attack, and clear tumor cells, playing a central role in antitumor immune response. However, certain immune cells can impair this response and help tumor immune escape. Therefore, exploring the factors that influence T-cell infiltration is crucial to understand tumor immunity and improve therapeutic effect of cancer immunotherapy. The use of single-cell RNA sequencing (scRNA-seq) allows the high-resolution analysis of the precise composition of immune cells with different phenotypes and other microenvironmental factors, including non-immune stromal cells and the related molecules in the tumor microenvironment of various cancer types. In this review, we summarized the research progress on T-cell infiltration and the crosstalk of other stromal cells and cytokines during T-cell infiltration using scRNA-seq to provide insights into the mechanisms regulating T-cell infiltration and contribute new perspectives on tumor immunotherapy.

Integrating single-cell and bulk RNA sequencing reveals CK19 + cancer stem cells and their specific SPP1 + tumor-associated macrophage niche in HBV-related hepatocellular carcinoma

PURPOSE

Cytokeratin 19-positive cancer stem cells (CK19 + CSCs) and their tumor-associated macrophages (TAMs) have not been fully explored yet in the hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC).

EXPERIMENTAL DESIGN

Single-cell RNA sequencing was performed on the viable cells obtained from 11 treatment-naïve HBV-associated HCC patients, including 8 CK19 + patients, to elucidate their transcriptomic landscape, CK19 + CSC heterogeneity, and immune microenvironment. Two in-house primary HCC cohorts (96 cases-related HBV and 89 cases with recurrence), TCGA external cohort, and in vitro and in vivo experiments were used to validate the results.

RESULTS

A total of 64,581 single cells derived from the human HCC and adjacent normal tissues were sequenced, and 11 cell types were identified. The result showed that CK19 + CSCs were phenotypically and transcriptionally heterogeneous, co-expressed multiple hepatics CSC markers, and were positively correlated with worse prognosis. Moreover, the SPP1 + TAMs (TAM_SPP1) with strong M2-like features and worse prognosis were specifically enriched in the CK19 + HCC and promoted tumor invasion and metastasis by activating angiogenesis. Importantly, matrix metalloproteinase 9 (MMP9) derived from TAM_SPP1, as the hub gene of CK19 + HCC, was activated by the VEGFA signal.

CONCLUSIONS

This study revealed the heterogeneity and stemness characteristics of CK19 + CSCs and specific immunosuppressive TAM_SPP1 in CK19 + HCC. The VEGFA signal can activate TAM_SPP1-derived MMP9 to promote the invasion and metastasis of CK19 + HCC tumors. This might provide novel insights into the clinical treatment of HCC patients.

Novel tumor-associated macrophage populations and subpopulations by single cell RNA sequencing

Tumor-associated macrophages (TAMs) are present in almost all solid tumor tissues. 16They play critical roles in immune regulation, tumor angiogenesis, tumor stem cell activation, tumor invasion and metastasis, and resistance to therapy. However, it is unclear how TAMs perform these functions. With the application of single-cell RNA sequencing (scRNA-seq), it has become possible to identify TAM subpopulations associated with distinct functions. In this review, we discuss four novel TAM subpopulations in distinct solid tumors based on core gene signatures by scRNA-seq, including FCN1 +, SPP1 +, C1Q + and CCL18 + TAMs. Functional enrichment and gene expression in scRNA-seq data from different solid tumor tissues found that FCN1 + TAMs may induce inflammation; SPP1 + TAMs are potentially involved in metastasis, angiogenesis, and cancer cell stem cell activation, whereas C1Q + TAMs participate in immune regulation and suppression; And CCL18 + cells are terminal immunosuppressive macrophages that not only have a stronger immunosuppressive function but also enhance tumor metastasis. SPP1 + and C1Q + TAM subpopulations can be further divided into distinct populations with different functions. Meanwhile, we will also present emerging evidence highlighting the separating macrophage subpopulations associated with distinct functions. However, there exist the potential disconnects between cell types and subpopulations identified by scRNA-seq and their actual function.

Single-cell analysis of immune and stroma cell remodeling in clear cell renal cell carcinoma primary tumors and bone metastatic lesions

BACKGROUND

Despite therapeutic advances, once a cancer has metastasized to the bone, it represents a highly morbid and lethal disease. One third of patients with advanced clear cell renal cell carcinoma (ccRCC) present with bone metastasis at the time of diagnosis. However, the bone metastatic niche in humans, including the immune and stromal microenvironments, has not been well-defined, hindering progress towards identification of therapeutic targets.

METHODS

We collected fresh patient samples and performed single-cell transcriptomic profiling of solid metastatic tissue (Bone Met), liquid bone marrow at the vertebral level of spinal cord compression (Involved), and liquid bone marrow from a different vertebral body distant from the tumor site but within the surgical field (Distal), as well as bone marrow from patients undergoing hip replacement surgery (Benign). In addition, we incorporated single-cell data from primary ccRCC tumors (ccRCC Primary) for comparative analysis.

RESULTS

The bone marrow of metastatic patients is immune-suppressive, featuring increased, exhausted CD8 + cytotoxic T cells, T regulatory cells, and tumor-associated macrophages (TAM) with distinct transcriptional states in metastatic lesions. Bone marrow stroma from tumor samples demonstrated a tumor-associated mesenchymal stromal cell population (TA-MSC) that appears to be supportive of epithelial-to mesenchymal transition (EMT), bone remodeling, and a cancer-associated fibroblast (CAFs) phenotype. This stromal subset is associated with poor progression-free and overall survival and also markedly upregulates bone remodeling through the dysregulation of RANK/RANKL/OPG signaling activity in bone cells, ultimately leading to bone resorption.

CONCLUSIONS

These results provide a comprehensive analysis of the bone marrow niche in the setting of human metastatic cancer and highlight potential therapeutic targets for both cell populations and communication channels.

A Comprehensive Comparison between Primary Liver Cancer and Liver Metastases through scRNA-Seq Data Analysis

Metastasis is one of the leading causes of cancer-related deaths. A comprehensive comparison of the differences between primary and metastatic cancers within the same organ can aid in understanding the growth mechanisms of cancer cells at metastatic sites, thereby helping to develop more effective targeted treatment strategies. Primary liver cancer is one of the most common types of cancer, and the liver is also one of the main metastatic sites. In this paper, we utilize single-cell RNA-Seq data to compare primary liver cancer and colorectal liver metastases from multiple perspectives, including cell types and proportions, activity of various cell types, cell-cell communication, mRNA expression differences within the same types of cells, key factors associated with cell proliferation, etc. Our analysis results show the following: (i) Compared to primary tissue, metastatic tissue contains more cytotoxic T cells and exhausted T cells, and it retains some specific characteristics of the primary site. (ii) Cells of the same type exhibit functional differences between primary and metastatic cancers, with metastatic cancer cells showing lower metabolism levels and immune cells exhibiting stronger immune activity. (iii) Interactions between monocytes and hepato-associated cells are strong in primary cancer, while depleted T cells frequently communicate with hepatocytes in metastatic cancer. (iv) Proliferation-related genes in primary and metastatic cancers are mainly involved in cell energy supply and basic metabolism activity, respectively.