Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has high relapse and low 5-year survival rates. Single-cell profiling in relapsed HCC may aid in the design of effective anticancer therapies, including immunotherapies. We profiled the transcriptomes of ∼17,000 cells from 18 primary or early-relapse HCC cases. Early-relapse tumors have reduced levels of regulatory T cells, increased dendritic cells (DCs), and increased infiltrated CD8+ T cells, compared with primary tumors, in two independent cohorts. Remarkably, CD8+ T cells in recurrent tumors overexpressed KLRB1 (CD161) and displayed an innate-like low cytotoxic state, with low clonal expansion, unlike the classical exhausted state observed in primary HCC. The enrichment of these cells was associated with a worse prognosis. Differential gene expression and interaction analyses revealed potential immune evasion mechanisms in recurrent tumor cells that dampen DC antigen presentation and recruit innate-like CD8+ T cells. Our comprehensive picture of the HCC ecosystem provides deeper insights into immune evasion mechanisms associated with tumor relapse.

Toward Ultrastable Metal Anode/Li6PS5Cl Interface via an Interlayer as Li Reservoir.

All-solid-state sulfide-based Li metal batteries are promising candidates for energy storage systems. However, thorny issues associated with undesired reactions and contact failure at the anode interface hinder their commercialization. Herein, an indium foil was endowed with a formed interlayer whose surface film is enriched with LiF and LiIn phases via a feasible prelithiation route. The lithiated alloy of the interlayer can regulate Li+ flux and charge distribution as a Li reservoir, benefiting uniform Li deposition. Meanwhile, it can suppress the reductive decomposition of the Li6PS5Cl electrolyte and maintain sufficient solid-solid contact. In situ impedance spectra reveal that constant interface impedance and fast charge transfer are realized by the interlayer. Further, long-term Li stripping/plating over 2000 h at 2.55 mA cm-2 is demonstrated by this anode. All-solid-state cells employing a LiCoO2 cathode and the Pre In anode can work for over 700 cycles with a capacity retention of 96.15% at 0.5 C.

Investigations of sequencing data and sample type on HLA class Ia typing with different computational tools.

Human leukocyte antigen (HLA) can encode the human major histocompatibility complex (MHC) proteins and play a key role in adaptive and innate immunity. Emerging clinical evidences suggest that the presentation of tumor neoantigens and neoantigen-specific T cell response associated with MHC class I molecules are of key importance to activate the adaptive immune systemin cancer immunotherapy. Therefore, accurate HLA typing is very essential for the clinical application of immunotherapy. In this study, we conducted performance evaluations of 4 widely used HLA typing tools (OptiType, Phlat, Polysolver and seq2hla) for predicting HLA class Ia genes from WES and RNA-seq data of 28 cancer patients. HLA genotyping data using PCR-SBT method was firstly obtained as the golden standard and was subsequently compared with HLA typing data by using NGS techniques. For both WES data and RNA-seq data, OptiType showed the highest accuracy for HLA-Ia typing than the other 3 programs at 2-digit and 4-digit resolution. Additionally, HLA typing accuracy from WES data was higher than from RNA-seq data (99.11% for WES data versus 96.42% for RNA-seq data). The accuracy of HLA-Ia typing by OptiType can reach 100% with the average depth of HLA gene regions >20x. Besides, the accuracy of 2-digit and 4-digit HLA-Ia typing based on control samples was higher than tumor tissues. In conclusion, OptiType by using WES data from control samples with the high average depth (>20x) of HLA gene regions can present a probably superior performance for HLA-Ia typing, enabling its application in cancer immunotherapy.

RRS1 shapes robust root system to enhance drought resistance in rice.

A strong root system facilitates the absorption of water and nutrients from the soil, to improve the growth of crops. However, to date, there are still very few root development regulatory genes that can be used in crop breeding for agriculture. In this study, we cloned a negative regulator gene of root development, Robust Root System 1 (RRS1), which encodes an R2R3-type MYB family transcription factor. RRS1 knockout plants showed enhanced root growth, including longer root length, longer lateral root length, and larger lateral root density. RRS1 represses root development by directly activating the expression of OsIAA3 which is involved in the auxin signaling pathway. A natural variation in the coding region of RRS1 changes the transcriptional activity of its protein. RRS1T allele, originating from wild rice, possibly increases root length by means of weakening regulation of OsIAA3. Knockout of RRS1 enhances drought resistance by promoting water absorption and improving water use efficiency. This study provides a new gene resource for improving root systems and cultivating drought-resistant rice varieties with important values in agricultural applications.

Chemically Induced Activity Recovery of Isolated Lithium in Anode-free Lithium Metal Batteries.

The anode-free lithium metal battery is considered to be an excellent candidate for the new generation energy storage system because of its higher energy density and safety than the traditional lithium metal battery. However, the continuous generation of SEI or isolated Li hinders its practical application. In general, the isolated Li is considered electrochemically inactive because it loses electrical connection with the current collector. Here we show an abnormal phenomenon that the lost capacity appears to be recovered after cycles when the isolated Li reconnects with a deposited Li metal layer. The isolated Li reconnection is ascribed to the chemical induction of the block copolymer coating. The migration of Li+ is affected by the electron delocalization and the electron cloud density of the polymer, which determine the conversion direction of Li+. Based on the mechanism, we propose a strategy to slow down the capacity decay of the anode-free lithium metal battery.

Genomic correlates of programmed cell death ligand 1 (PD-L1) expression in Chinese lung adenocarcinoma patients.

BACKGROUND: Although PD-L1 expression is a crucial predictive biomarker for immunotherapy, it can be influenced by many factors. METHODS: A total of 248 Chinese patients with lung adenocarcinoma was retrospectively identified. Data for clinical features, gene alternations, signaling pathways and immune signatures was analyzed among negative expression group (TPS < 1%, n = 124), intermediate expression group (1% ≤ TPS < 50%, n = 93), and high expression group (TPS ≥ 50%, n = 38). Clinical outcomes among different expression groups were also evaluated from public database. RESULTS: Firstly, high tumor mutation burden was significantly associated with high PD-L1 expression in these Chinese patients with lung adenocarcinoma. In addition, gene alternations including TP53, PRKDC, KMT2D, TET1 and SETD2 apparently occurred in high PD-L1 expression group. Moreover, pathway analysis showed that mutations involving in DDR pathway, TP53 pathway, cell-cycle pathway and NOTCH pathway were obviously varied among three PD-L1 expression groups. Besides, most of patients in high PD-L1 expression group from TCGA database were determined as high-grade immune subtypes (C2-C4), showing significant higher proportions of IFN-gamma, CD8+ T-cells, NK cells, NK CD56 dim cells, Th1 cells, Th2 cells (P < 0.0001). Moreover, SETD2 mutation slightly correlated with overall survival from MSKCC cohort (HR 1.92 [95%CI 0.90-4.10], P = 0.085), and the percentage of IFN-gamma was significantly higher in SETD2 mutant group than in wild-type group (P < 0.01). CONCLUSIONS: This study illustrated in-depth genomic correlates of PD-L1 expression in Chinese lung adenocarcinoma patients and relevant immune signatures from public database, which might interpret more potential molecular mechanisms for immunotherapy in NSCLC.

Developing a 5-gene prognostic signature for cervical cancer by integrating mRNA and copy number variations.

BACKGROUND: Cervical cancer is frequently detected gynecological cancer all over the world. This study was designed to develop a prognostic signature for an effective prediction of cervical cancer prognosis. METHODS: Differentially expressed genes (DEGs) were identified based on copy number variation (CNV) data and expression profiles from different databases. A prognostic model was constructed and further optimized by stepwise Akaike information criterion (stepAIC). The model was then evaluated in three groups (training group, test group and validation group). Functional analysis and immune analysis were used to assess the difference between high-risk and low-risk groups. RESULTS: The study developed a 5-gene prognostic model that could accurately classify cervical cancer samples into high-risk and low-risk groups with distinctly different prognosis. Low-risk group exhibited more favorable prognosis and higher immune infiltration than high-risk group. Both univariate and multivariate Cox regression analysis showed that the risk score was an independent risk factor for cervical cancer. CONCLUSIONS: The 5-gene prognostic signature could serve as a predictor for identifying high-risk cervical cancer patients, and provided potential direction for studying the mechanism or drug targets of cervical cancer. The integrated analysis of CNV and mRNA expanded a new perspective for exploring prognostic signatures in cervical cancer.

Two Mn(II) and Zn(II) Coordination Polymers: Aqueous-Phase Detection of Nitroaromatic Explosives and Protective Effect on Atherosclerosis.

Via the solvothermal reaction between Zn(II) or Mn(II) salts and 5-(3,4-dicarboxylphenoxy)nicotinic acid (H3L) ligand, a trifunctional N,O-building block having three diverse kinds of functional groups (O-ether, N-pyridyl and COOH), two new coordination polymers (CPs) could be generated, and their chemical formulae respectively are {[Mn3(L)2(H2O)2]·4H2O} (1) and {[Zn(HL)]·NMP} (2). The complex 2 based on Zn(II) possesses high efficiency of fluorescence quenching for the nitrophenol (2,4,6-trinitrophenol, TNP; 4-nitrophenol, 4-NP; 3-nitrophenol, 3-NP; 2-nitrophenol, 2-NP) in the aqueous solution. Furthermore, the treatment activity of compounds on the atherosclerosis was assessed, and relevant mechanism was investigated. First of all, the ELISA assay was used to measure the content of the inflammatory cytokines released into the plasma. Besides, the levels of the NF-κb signaling pathway in the vascular endothelial cells were measured with real time RT-PCR. The hemolysis test was conducted in this research to measure the biocompatibility of the new compound.

Differential expression of miR-142-3p protects cardiomyocytes from myocardial ischemia-reperfusion via TLR4/NFkB axis.

Our research aims to explore the impact of miR-142 on myocardial apoptosis in the mouse ischemia and reperfusion (IR) model and investigate the underlying mechanisms at the molecular level. A considerable downregulation of miR-142 was observed in the cardiac area of mice post IR modeling. To understand the regulatory function of IR-induced miR-142 downregulation, the animals were categorized into four groups: IR model group; IR + agomir-142 group (IR mice treated with agomir-142); IR + antagomir-142 group (IR mice treated with antagomir-142); IR + agomir-142 + negative control (NC) group (IR mice processed with agomir-NC). The results indicated that agomir-142 upregulation was capable of shrinking IR damage-triggered infarction of the ventriculus sinister, strengthening myocardial function, and guarding against cardiomyocyte apoptosis, whereas further decreased miR-142 with antagomir-142 infection displayed negative influence of miR-142 against mice IR damage. In the cellular assay, miR-142 overexpression significantly improved proliferation and inhibited the apoptosis of neonatal rat cardiomyocytes (NRCs). Moreover, we found that miR-142 reduced the Bcl-2/Bax ratio and upregulated hydrogen peroxide (H2 O2 )-induced caspase-3 expression. Furthermore, transfection with an miR-142 mimic prevented the upregulation of TLR4/NFkB expression and activation in H2 O2 -treated NRCs. Our findings also revealed that miR-142 is linked to the 3'-untranslated area of the TLR4 gene. In addition, TLR4 overexpression considerably ablated the protective effects of miR-142 in terms of the cell viability of H2 O2 -treated NRCs. Taken together, miR-142 agomir injection in mice and miR-142 mimic transfection in NRCs plays a role in protecting the heart from IR damage and malfunction via the TLR4/NFkB axis both in vivo and in vitro.

Stress hyperglycemia ratio and the clinical outcome of patients with heart failure: a meta-analysis.

Background: Stress hyperglycemia ratio (SHR) is a newly suggested measure of stress-induced hyperglycemia that combines both short-term and long-term glycemic conditions. The study aimed to explore the association between SHR and the incidence of adverse clinical events with heart failure (HF) through a meta-analysis. Methods: Cohort studies relevant to the aim of the meta-analysis were retrieved by search of electronic databases including PubMed, Web of Science, Embase, Wanfang, and CNKI. A random-effects model was used to combine the data by incorporating the influence of between-study heterogeneity. Results: Ten studies involving 15250 patients with HF were included. Pooled results showed that compared to patients with lower SHR at baseline, those with a higher SHR were associated with an increased risk of all-cause mortality during follow-up (risk ratio [RR]: 1.61, 95% confidence interval [CI]: 1.17 to 2.21, p = 0.003; I2 = 82%). Further meta-regression analysis suggests that different in the cutoff of SHR significantly modify the results (coefficient = 1.22, p = 0.02), and the subgroup analysis suggested a more remarkable association between SHR and all-cause mortality in studies with cutoff of SHR ≥ 1.05 than those with cutoff of SHR < 1.05 (RR: 2.29 versus 1.08, p for subgroup difference < 0.001). Subsequent meta-analyses also showed that a high SHR at baseline was related to the incidence of cardiovascular death (RR: 2.19, 95% CI: 1.55 to 3.09, p < 0.001; I2 = 0%), HF-rehospitalization (RR: 1.83, 95% CI: 1.44 to 2.33, p < 0.001; I2 = 0%), and major adverse cardiovascular events (RR: 1.54, 95% CI: 1.15 to 2.06, p = 0.004; I2 = 74%) during follow-up. Conclusion: A high SHR at baseline is associated with a poor clinical prognosis of patients with HF. Systematic review registration: https://inplasy.com, identifier INPLASY202430080.

Haematococcus pluvialis polysaccharides improve microbiota-driven gut epithelial and vascular barrier and prevent alcoholic steatohepatitis development.

Algal polysaccharides possess many biological activities and health benefits, such as antioxidant, anti-tumor, anti-coagulant, and immunomodulatory potential. Gut microbiota has emerged as one of the major contributor in mediating the health benefits of algal polysaccharides. In this study we showed that Haematococcus pluvialis polysaccharides (HPP) decreased serum transaminase levels and hepatic triglyceride content, alleviated inflammation and oxidative stress in the liver of chronic and binge ethanol diet-fed mice. Furthermore, HPP reduced endotoxemia, improved gut microbiota dysbiosis, inhibited epithelial barrier disruption and gut vascular barrier (GVB) damage in ethanol diet-fed mice. Co-housing vehicle-fed mice with HPP-fed mice alleviated ethanol-induced liver damage and endotoxemia. Moreover, fecal microbiota transplantation from HPP-fed mice into antibiotic-induced microbiota-depleted recipients also alleviated ethanol-induced liver injury and improved gut epithelial and vascular barrier. Our study demonstrated that HPP ameliorated ethanol-induced gut epithelial and vascular barrier dysfunction through alteration of gut microbiota, therefore preventing alcoholic liver damage.

Impact of deoxynivalenol on rumen function, production, and health of dairy cows: Insights from metabolomics and microbiota analysis.

Deoxynivalenol contamination in feed and food, pervasive from growth, storage, and processing, poses a significant risk to dairy cows, particularly when exposed to a high-starch diet; however, whether a high-starch diet exacerbates these negative effects remains unclear. Therefore, we investigated the combined impact of deoxynivalenol and dietary starch on the production performance, rumen function, and health of dairy cows using metabolomics and 16 S rRNA sequencing. Our findings suggested that both high- and low-starch diets contaminated with deoxynivalenol significantly reduced the concentration of propionate, isobutyrate, valerate, total volatile fatty acids (TVFA), and microbial crude protein (MCP) concentrations, accompanied by a noteworthy increase in NH3-N concentration in vitro and in vivo (P < 0.05). Deoxynivalenol altered the abundance of microbial communities in vivo, notably affecting Oscillospiraceae, Lachnospiraceae, Desulfovibrionaceae, and Selenomonadaceae. Additionally, it significantly downregulated lecithin, arachidonic acid, valine, leucine, isoleucine, arginine, and proline metabolism (P < 0.05). Furthermore, deoxynivalenol triggered oxidative stress, inflammation, and dysregulation in immune system linkage, ultimately compromising the overall health of dairy cows. Collectively, both high- and low-starch diets contaminated with deoxynivalenol could have detrimental effects on rumen function, posing a potential threat to production performance and the overall health of cows. Notably, the negative effects of deoxynivalenol are more pronounced with a high-starch diet than a low-starch diet.

Metabolomic analysis of rumen-protected branched-chain amino acids in primiparous dairy cows.

Introduction: Peripartal cows are susceptible to a negative energy balance due to inadequate nutrient intake and high energy requirements for lactation. Improving the energy metabolism of perinatal dairy cows is crucial in increasing production in dairy cows. Methods: In this study, we investigated the impact of rumen-protected branched-chain amino acid (RPBCAA) on the production performance, energy and lipid metabolism, oxidative stress, and immune function of primiparous dairy cows using metabolomics through a single-factor experiment. Twenty healthy primiparous Holstein cows were selected based on body condition scores and expected calving date, and were randomly divided into RPBCAA (n = 10) and control (n = 10) groups. The control group received a basal diet from calving until 21 d in milk, and the RPBCAA group received the basal diet and 44.6 g/d RPLeu, 25.14 g/d RPIle, and 25.43 g/d RPVal. Results: In comparison to the control group, the supplementation of RPBCAA had no significant effect on milk yield and milk composition of the dairy cows. Supplementation with RPBCAA significantly increased the concentrations of insulin, insulin growth factor 1, glucagon, and growth hormones, which are indicators of energy metabolism in postpartum cows. The very low density lipoprotein, fatty acid synthase, acetyl coenzyme A carboxylase, and hormone-sensitive lipase contents of the RPBCAA group were significantly greater than that of the control group; these metrics are related to lipid metabolism. In addition, RPBCAA supplementation significantly increased serum glutathione peroxidase and immunoglobulin G concentrations and decreased malondialdehyde concentrations. Liquid chromatography-mass spectrometry (LC-MS) analysis revealed 414 serum and 430 milk metabolic features. Supplementation with RPBCAA primarily increased concentrations of amino acid and lipid metabolism pathways and upregulated the abundance of serotonin, glutamine, and phosphatidylcholines. Discussion: In summary, adding RPBCAA to the daily ration can influence endocrine function and improve energy metabolism, regulate amino acid and lipid metabolism, mitigate oxidative stress and maintain immune function on primiparous cows in early lactation.

Electrochemical Performances of Nickel-Rich Single-Crystal LiNi0.83 Co0.12 Mn0.05 O2 Cathode Material for Lithium-Ion Batteries Synthesized by Tuning Li+ /Ni2+ Mixing.

Single-crystal nickel-rich materials are promising alternatives to polycrystalline cathodes owing to their excellent structure stability and cycle performance while the cathode material usually appears high cation mixing, which may have a negative effect on its electrochemical performance. The study presents the structural evolution of single-crystal LiNi0.83 Co0.12 Mn0.05 O2 in the temperature-composition space using temperature-resolved in situ XRD and the cation mixing is tuned to improve electrochemical performances. The as-synthesized single-crystal sample shows high initial discharge specific capacity (195.5 mAh g-1 at 1 C), and excellent capacity retention (80.1 % after 400 cycles at 1 C), taking account of lower structure disorder (Ni2+ occupying Li sites is 1.56 %) and integrated grains with an average of 2-3 µm. In addition, the single-crystal material also displays a superior rate capability of 159.1 mAh g-1 at the rate of 5 C. This excellent performance is attributed to the rapid Li+ transportation within the crystal structure with fewer Ni2+ cations in Li layer as well as intactly single grains. In sum, the regulation of Li+ /Ni2+ mixing provides a feasible strategy for boosting single-crystal nickel-rich cathode material.

Molecular characteristics of multiple primary pulmonary nodules under a three-dimensional reconstruction model and relevant multi-omics analyses: a case report.

Background: In addition to CT images and pathological features, many other molecular characteristics remain unknown about multiple primary lung cancer (MPLC) from intrapulmonary metastatic lung cancer. Case presentation: In this study, we reported a patient with an early-stage MPLC with both adenocarcinoma in situ (AIS) subtype and minimally invasive adenocarcinoma (MIA) subtype. The patient was diagnosed with more than 10 nodules and underwent precise surgery assisted by three-dimensional (3D) reconstruction at the left upper lung lobe. Whole-exome sequencing (WES) and multiple immunohistochemistry (mIHC) were performed to reveal the genomic profiling and tumor microenvironments of multiple nodules in this patient with MPLC. Based on 3D reconstruction location information, we found that the genomic and pathological results of adjacent lymph nodes were quite different. On the other hand, PD-L1 expression and the proportion of infiltrating lymphocytes in tumor microenvironments were all at a low status and did not vary in adjacent lymph nodes. Additionally, maximum diameter and tumor mutational burden levels were found to be significantly associated with CD8+ T cell proportion (p<0.05). Besides, CD163+ macrophages and CD4+ T cell proportion were higher in MIA nodules than in AIS nodules (p<0.05). This patient reached a recurrence-free survival of 39 months. Conclusion: Generally, in addition to CT imaging and pathological results, genomic profiling and tumor microenvironments may facilitate identifying the potential molecular mechanisms and clinical outcomes in patients with early-stage MPLC.

Neoadjuvant therapy with immune checkpoint blockade, antiangiogenesis, and chemotherapy for locally advanced gastric cancer.

Despite neoadjuvant/conversion chemotherapy, the prognosis of cT4a/bN+ gastric cancer is poor. Immune checkpoint inhibitors (ICIs) and antiangiogenic agents have shown activity in late-stage gastric cancer, but their efficacy in the neoadjuvant/conversion setting is unclear. In this single-armed, phase II, exploratory trial (NCT03878472), we evaluate the efficacy of a combination of ICI (camrelizumab), antiangiogenesis (apatinib), and chemotherapy (S-1 ± oxaliplatin) for neoadjuvant/conversion treatment of cT4a/bN+ gastric cancer. The primary endpoints are pathological responses and their potential biomarkers. Secondary endpoints include safety, objective response, progression-free survival, and overall survival. Complete and major pathological response rates are 15.8% and 26.3%. Pathological responses correlate significantly with microsatellite instability status, PD-L1 expression, and tumor mutational burden. In addition, multi-omics examination reveals several putative biomarkers for pathological responses, including RREB1 and SSPO mutation, immune-related signatures, and a peripheral T cell expansion score. Multi-omics also demonstrates dynamic changes in dominant tumor subclones, immune microenvironments, and T cell receptor repertoires during neoadjuvant immunotherapy. The toxicity and post-surgery complications are limited. These data support further validation of ICI- and antiangiogenesis-based neoadjuvant/conversion therapy in large randomized trials and provide candidate biomarkers.

Erratum to 'Evolution and genotypic characteristics of small cell lung cancer transformation in non-small cell lung carcinomas' [Journal of the National Cancer Center, 1 (2021), 4: 153-162].

[This corrects the article DOI: 10.1016/j.jncc.2021.11.001.].

Comprehensive Analysis of Histone Modifications in Hepatocellular Carcinoma Reveals Different Subtypes and Key Prognostic Models.

Histone modification, an important epigenetic mechanism, is related to the carcinogenesis of hepatocellular carcinoma (HCC). In three datasets, we screened 88 epigenetic-dysregulated PCGs (epi-PCGs) , which were significantly associated with HCC survival and could cluster HCC into three molecular subtypes. These subtypes were associated with prognosis, immunomodulatory alterations, and response to different treatment strategies. Based on 88 epi-PCGs in the TCGA training set, a risk prediction model composed of 4 epi-PCGs was established. The model was closely related to the clinicopathological features and showed a strong predictive ability in different clinical subgroups. In addition, the risk prediction model was an independent prognostic factor for patients with HCC. The significance of epi-PCGs in HCC is revealed by our data analysis.

BUB1 Is Identified as a Potential Therapeutic Target for Pancreatic Cancer Treatment.

Pancreatic cancer is one of the most challenging cancer types in clinical treatment worldwide. This study aimed to understand the tumorigenesis mechanism and explore potential therapeutic targets for patients with pancreatic cancer. Single-cell data and expression profiles of pancreatic cancer samples and normal tissues from multiple databases were included. Comprehensive bioinformatics analyses were applied to clarify tumor microenvironment and identify key genes involved in cancer development. Immense difference of cell types was shown between tumor and normal samples. Four cell types (B cell_1, B cell_2, cancer cell_3, and CD1C+_B dendritic cell_3) were screened to be significantly associated with prognosis. Three ligand-receptor pairs, including CD74-MIF, CD74-COPA, and CD74-APP, greatly contributed to tumorigenesis. High expression of BUB1 (BUB1 Mitotic Checkpoint Serine/Threonine Kinase) was closely correlated with worse prognosis. CD1C+_B dendritic cell_3 played a key role in tumorigenesis and cancer progression possibly through CD74-MIF. BUB1 can serve as a prognostic biomarker and a therapeutic target for patients with pancreatic cancer. The study provided a novel insight into studying the molecular mechanism of pancreatic cancer development and proposed a potential strategy for exploiting new drugs.

Exploring the Interplay between Metabolism and Tumor Microenvironment Based on Four Major Metabolism Pathways in Colon Adenocarcinoma.

Tumor metabolism plays a critical role in tumor progression. However, the interaction between metabolism and tumor microenvironment (TME) has not been comprehensively revealed in colon adenocarcinoma (COAD). We used unsupervised consensus clustering to establish three molecular subtypes (clusters) based on the enrichment score of four major metabolism pathways in TCGA-COAD dataset. GSE17536 was used as a validation dataset. Single-cell RNA sequencing data (GSE161277) was employed to further verify the reliability of subtyping and characterize the correlation between metabolism and TME. Three clusters were identified and they performed distinct prognosis and molecular features. Clust3 had the worst overall survival and the highest enrichment score of glycolysis. 86 differentially expressed genes (DEGs) were identified, in which 11 DEGs were associated with favorable prognosis and 75 DEGs were associated with poor prognosis. Striking correlations were observed between hypoxia and glycolysis, clust3 and hypoxia, and clust3 and angiogenesis (P < 0.001).We constructed a molecular subtyping system which was effective and reliable for predicting COAD prognosis. The 86 identified key DEGs may be greatly involved in COAD progression, and they provide new perspectives and directions for further understanding the mechanism of metabolism in promoting aggressive phenotype by interacting with TME.