ABSTRACT
Blocking PD-1/PD-L1 signaling transforms cancer therapy and is assumed to unleash exhausted tumor-reactive CD8+ T cells in the tumor microenvironment (TME). However, recent studies have also indicated that the systemic tumor-reactive CD8+ T cells may respond to PD-1/PD-L1 immunotherapy. These discrepancies highlight the importance of further defining tumor-specific CD8+ T cell responders to PD-1/PD-L1 blockade. Here, using multiple preclinical tumor models, we revealed that a subset of tumor-specific CD8+ cells in the tumor draining lymph nodes (TdLNs) was not functionally exhausted but exhibited canonical memory characteristics. TdLN-derived tumor-specific memory (TTSM) cells established memory-associated epigenetic program early during tumorigenesis. More importantly, TdLN-TTSM cells exhibited superior anti-tumor therapeutic efficacy after adoptive transfer and were characterized as bona fide responders to PD-1/PD-L1 blockade. These findings highlight that TdLN-TTSM cells could be harnessed to potentiate anti-tumor immunotherapy.
Subject(s)
B7-H1 Antigen , Neoplasms , Humans , Programmed Cell Death 1 Receptor , CD8-Positive T-Lymphocytes , Immune Checkpoint Inhibitors , Tumor Microenvironment , Neoplasms/therapy , Neoplasms/pathology , Lymph Nodes/pathologyABSTRACT
In this article, we develop CausalEGM, a deep learning framework for nonlinear dimension reduction and generative modeling of the dependency among covariate features affecting treatment and response. CausalEGM can be used for estimating causal effects in both binary and continuous treatment settings. By learning a bidirectional transformation between the high-dimensional covariate space and a low-dimensional latent space and then modeling the dependencies of different subsets of the latent variables on the treatment and response, CausalEGM can extract the latent covariate features that affect both treatment and response. By conditioning on these features, one can mitigate the confounding effect of the high dimensional covariate on the estimation of the causal relation between treatment and response. In a series of experiments, the proposed method is shown to achieve superior performance over existing methods in both binary and continuous treatment settings. The improvement is substantial when the sample size is large and the covariate is of high dimension. Finally, we established excess risk bounds and consistency results for our method, and discuss how our approach is related to and improves upon other dimension reduction approaches in causal inference.
ABSTRACT
Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.
Subject(s)
Cell-Free Nucleic Acids , Deep Learning , Humans , Cell-Free Nucleic Acids/genetics , DNA Methylation , Biomarkers , Promoter Regions, Genetic , Biomarkers, Tumor/geneticsABSTRACT
Deciphering 3D genome conformation is important for understanding gene regulation and cellular function at a spatial level. The recent advances of single cell Hi-C technologies have enabled the profiling of the 3D architecture of DNA within individual cell, which allows us to study the cell-to-cell variability of 3D chromatin organization. Computational approaches are in urgent need to comprehensively analyze the sparse and heterogeneous single cell Hi-C data. Here, we proposed scDEC-Hi-C, a new framework for single cell Hi-C analysis with deep generative neural networks. scDEC-Hi-C outperforms existing methods in terms of single cell Hi-C data clustering and imputation. Moreover, the generative power of scDEC-Hi-C could help unveil the differences of chromatin architecture across cell types. We expect that scDEC-Hi-C could shed light on deepening our understanding of the complex mechanism underlying the formation of chromatin contacts.
Subject(s)
Chromatin , Chromosomes , Chromatin/genetics , Genome , DNA , Cluster AnalysisABSTRACT
BACKGROUND AND AIMS: RAD51 recombinase (RAD51) is a highly conserved DNA repair protein and is indispensable for embryonic viability. As a result, the role of RAD51 in liver development and function is unknown. Our aim was to characterize the function of RAD51 in postnatal liver development. APPROACH AND RESULTS: RAD51 is highly expressed during liver development and during regeneration following hepatectomy and hepatic injury, and is also elevated in chronic liver diseases. We generated a hepatocyte-specific Rad51 deletion mouse model using Alb -Cre ( Rad51 -conditional knockout (CKO)) and Adeno-associated virus 8-thyroxine-binding globulin-cyclization recombination enzyme to evaluate the function of RAD51 in liver development and regeneration. The phenotype in Rad51 -CKO mice is dependent on CRE dosage, with Rad51fl/fl ; Alb -Cre +/+ manifesting a more severe phenotype than the Rad51fl/fl ; Alb -Cre +/- mice. RAD51 deletion in postnatal hepatocytes results in aborted mitosis and early onset of pathological polyploidization that is associated with oxidative stress and cellular senescence. Remarkable liver fibrosis occurs spontaneously as early as in 3-month-old Rad51fl/fl ; Alb -Cre +/+ mice. While liver regeneration is compromised in Rad51 -CKO mice, they are more tolerant of carbon tetrachloride-induced hepatic injury and resistant to diethylnitrosamine/carbon tetrachloride-induced HCC. A chronic inflammatory microenvironment created by the senescent hepatocytes appears to activate ductular reaction the transdifferentiation of cholangiocytes to hepatocytes. The newly derived RAD51 functional immature hepatocytes proliferate vigorously, acquire increased malignancy, and eventually give rise to HCC. CONCLUSIONS: Our results demonstrate a novel function of RAD51 in liver development, homeostasis, and tumorigenesis. The Rad51 -CKO mice represent a unique genetic model for premature liver senescence, fibrosis, and hepatocellular carcinogenesis.
ABSTRACT
Elucidating the role of 3D architecture of DNA in gene regulation is crucial for understanding cell differentiation, tissue homeostasis and disease development. Among various chromatin conformation capture methods, HiChIP has received increasing attention for its significant improvement over other methods in profiling of regulatory (e.g. H3K27ac) and structural (e.g. cohesin) interactions. To facilitate the studies of 3D regulatory interactions, we developed a HiChIP interactions database, HiChIPdb (http://health.tsinghua.edu.cn/hichipdb/). The current version of HiChIPdb contains â¼262M annotated HiChIP interactions from 200 high-throughput HiChIP samples across 108 cell types. The functionalities of HiChIPdb include: (i) standardized categorization of HiChIP interactions in a hierarchical structure based on organ, tissue and cell line and (ii) comprehensive annotations of HiChIP interactions with regulatory genes and GWAS Catalog SNPs. To the best of our knowledge, HiChIPdb is the first comprehensive database that utilizes a unified pipeline to map the functional interactions across diverse cell types and tissues in different resolutions. We believe this database has the potential to advance cutting-edge research in regulatory mechanisms in development and disease by removing the barrier in data aggregation, preprocessing, and analysis.
Subject(s)
Chromatin , DNA , Cell Line , Chromatin/genetics , Gene Expression Regulation , Sequence Analysis, DNA/methods , Databases, GeneticABSTRACT
SignificanceThe present work might be significant for exploring advanced K-ion batteries with superb rate capability and cycle stability toward practical applications. The as-assembled K-ion half cell exhibits an excellent rate capability of 428 mA h g-1 at 100 mA g-1 and a high reversible specific capacity of 330 mA h g-1 with 120% specific capacity retention after 2,000 cycles at 2,000 mA g-1, which is the best among those based on carbon materials. The as-constructed full cell delivers 98% specific capacity retention over 750 cycles at 500 mA g-1, superior to most of those based on carbon materials that have been reported thus far.
ABSTRACT
Dendrobium officinale is a rare and precious medicinal plant. Southern blight is a destructive disease in the artificial cultivation of D. officinale, and one of its pathogens is Sclerotium delphinii. S. delphinii is a phytopathogenic fungus with a wide host range with extremely strong pathogenicity. In this study, S. delphinii was isolated from D. officinale with southern blight. Subsequently, this specific strain underwent thorough whole-genome sequencing using the PacBio Sequel II platform, which employed single-molecule real-time (SMRT) technology. Comprehensive annotations were obtained through functional annotation of protein sequences using various publicly available databases. The genome of S. delphinii measures 73.66 Mb, with an N90 contig size of 2,707,110 bp, and it contains 18,506 putative predictive genes. This study represents the first report on the genome size assembly and annotation of S. delphinii, making it the initial species within the Sclerotium genus to undergo whole-genome sequencing, which can provide solid data and a theoretical basis for further research on the pathogenesis, omics of S. delphinii.
ABSTRACT
MOTIVATION: The prediction of cancer drug response is a challenging subject in modern personalized cancer therapy due to the uncertainty of drug efficacy and the heterogeneity of patients. It has been shown that the characteristics of the drug itself and the genomic characteristics of the patient can greatly influence the results of cancer drug response. Therefore, accurate, efficient, and comprehensive methods for drug feature extraction and genomics integration are crucial to improve the prediction accuracy. RESULTS: Accurate prediction of cancer drug response is vital for guiding the design of anticancer drugs. In this study, we propose an end-to-end deep learning model named DeepAEG which is based on a complete-graph update mode to predict IC50. Specifically, we integrate an edge update mechanism on the basis of a hybrid graph convolutional network to comprehensively learn the potential high-dimensional representation of topological structures in drugs, including atomic characteristics and chemical bond information. Additionally, we present a novel approach for enhancing simplified molecular input line entry specification data by employing sequence recombination to eliminate the defect of single sequence representation of drug molecules. Our extensive experiments show that DeepAEG outperforms other existing methods across multiple evaluation parameters in multiple test sets. Furthermore, we identify several potential anticancer agents, including bortezomib, which has proven to be an effective clinical treatment option. Our results highlight the potential value of DeepAEG in guiding the design of specific cancer treatment regimens.
Subject(s)
Antineoplastic Agents , Neoplasms , Humans , Antineoplastic Agents/pharmacology , Antineoplastic Agents/therapeutic use , Neoplasms/drug therapy , Neoplasms/genetics , Bortezomib , Genomics , UncertaintyABSTRACT
It is well-recognized that blood flow at branches and bends of arteries generates disturbed shear stress, which plays a crucial in driving atherosclerosis. Flow-generated fluid shear stress (FSS), as one of the key hemodynamic factors, is appreciated for its critical involvement in regulating angiogenesis to facilitate wound healing and tissue repair. Endothelial cells can directly sense FSS but the mechanobiological mechanism by which they decode different patterns of FSS to trigger angiogenesis remains unclear. In the current study, laminar shear stress (LSS, 15 dyn/cm2) was employed to mimic physiological blood flow, while disturbed shear stress (DSS, ranging from 0.5 ± 4 dyn/cm2) was applied to simulate pathological conditions. The aim was to investigate how these distinct types of blood flow regulated endothelial angiogenesis. Initially, we observed that DSS impaired angiogenesis and downregulated endogenous vascular endothelial growth factor B (VEGFB) expression compared to LSS. We further found that the changes in membrane protein, migration and invasion enhancer 1 (MIEN1) play a role in regulating ERK/MAPK signaling, thereby contributing to endothelial angiogenesis in response to FSS. We also showed the involvement of MIEN1-directed cytoskeleton organization. These findings suggest the significance of shear stress in endothelial angiogenesis, thereby enhancing our understanding of the alterations in angiogenesis that occur during the transition from physiological to pathological blood flow.
Subject(s)
Angiogenesis , Endothelial Cells , Hemodynamics , Humans , Atherosclerosis/pathology , Cells, Cultured , Endothelial Cells/physiology , Intracellular Signaling Peptides and Proteins/metabolism , Neoplasm Proteins/metabolism , Stress, Mechanical , Vascular Endothelial Growth Factor B/metabolismABSTRACT
Rowell's syndrome is an autoimmune disease characterized by lupus erythematosus, erythema multiforme skin lesions, and speckled antinuclear antibody. We report the case of a woman who presented with erythema multiforme with target-type skin lesions and vulvar vegetation who fulfilled the criteria for Rowell's syndrome and condyloma acuminatum. The simultaneous occurrence of both conditions has rarely been reported in the literature.
Subject(s)
Erythema Multiforme , Lupus Erythematosus, Systemic , Female , Humans , Syndrome , Erythema Multiforme/diagnosis , Erythema Multiforme/pathology , Antibodies, AntinuclearABSTRACT
Small extracellular vesicles (sEVs) are heterogeneous biological nanoparticles (NPs) with wide biomedicine applications. Tracking individual nanoscale sEVs can reveal information that conventional microscopic methods may lack, especially in cellular microenvironments. This usually requires biolabeling to identify single sEVs. Here, we developed a light scattering imaging method based on dark-field technology for label-free nanoparticle diffusion analysis (NDA). Compared with nanoparticle tracking analysis (NTA), our method was shown to determine the diffusion probabilities of a single NP. It was demonstrated that accurate size determination of NPs of 41 and 120 nm in diameter is achieved by purified Brownian motion (pBM), without or within the cell microenvironments. Our pBM method was also shown to obtain a consistent size estimation of the normal and cancerous plasma-derived sEVs without and within cell microenvironments, while cancerous plasma-derived sEVs are statistically smaller than normal ones. Moreover, we showed that the velocity and diffusion coefficient are key parameters for determining the diffusion types of the NPs and sEVs in a cancerous cell microenvironment. Our light scattering-based NDA and pBM methods can be used for size determination of NPs, even in cell microenvironments, and also provide a tool that may be used to analyze sEVs for many biomedical applications.
Subject(s)
Extracellular Vesicles , Extracellular Vesicles/chemistry , Humans , Light , Nanoparticles/chemistry , Scattering, Radiation , Cellular Microenvironment , Particle Size , Diffusion , Tumor Microenvironment , Cell Line, Tumor , MotionABSTRACT
BACKGROUND: Generating elite rice varieties with high yield and superior quality is the main goal of rice breeding programs. Key agronomic traits, including grain size and seed germination characteristics, affect the final yield and quality of rice. The RGA1 gene, which encodes the α-subunit of rice G-protein, plays an important role in regulating rice architecture, seed size and abiotic stress responses. However, whether RGA1 is involved in the regulation of rice quality and seed germination traits is still unclear. RESULTS: In this study, a rice mutant small and round grain 5 (srg5), was identified in an EMS-induced rice mutant library. Systematic analysis of its major agronomic traits revealed that the srg5 mutant exhibited a semi-dwarf plant height with small and round grain and reduced panicle length. Analysis of the physicochemical properties of rice showed that the difference in rice eating and cooking quality (ECQ) between the srg5 mutant and its wild-type control was small, but the appearance quality was significantly improved. Interestingly, a significant suppression of rice seed germination and shoot growth was observed in the srg5 mutant, which was mainly related to the regulation of ABA metabolism. RGA1 was identified as the candidate gene for the srg5 mutant by BSA analysis. A SNP at the splice site of the first intron disrupted the normal splicing of the RGA1 transcript precursor, resulting in a premature stop codon. Additional linkage analysis confirmed that the target gene causing the srg5 mutant phenotype was RGA1. Finally, the introduction of the RGA1 mutant allele into two indica rice varieties also resulted in small and round rice grains with less chalkiness. CONCLUSIONS: These results indicate that RGA1 is not only involved in the control of rice architecture and grain size, but also in the regulation of rice quality and seed germination. This study sheds new light on the biological functions of RGA1, thereby providing valuable information for future systematic analysis of the G-protein pathway and its potential application in rice breeding programs.
Subject(s)
Oryza , Oryza/genetics , Seeds/genetics , Germination/genetics , Plant Breeding , Edible Grain/genetics , GTP-Binding ProteinsABSTRACT
Optimizing the active centers through reconstruction is recognized as the key to construct high-performance oxygen evolution reaction (OER) catalysts. Herein, a simple and rapid in situ leaching strategy to promote the self-reconstruction of NiFe-layered double hydroxides (LDHs) catalysts is reported. The trace Zn dopants are introduced in advance by a facile and one-step hydrothermal method, followed by leaching over the electrochemical activation process, which can remarkably reduce the formation potential of NiFeOOH active centers to enable the deeper self-reconstruction for the formation of abundant highly active centers. Moreover, the self-restructured NiFeOOH-VZn cannot only significantly lower the dehydrogenation energy barrier for the transformation from Ni(OH)2 to NiOOH, but also decrease the free energy barrier of rate determining step for the *OH converted to *O through a deprotonation process, thus significantly boosting the OER behaviors. As a proof of concept, the obtained NiFeOOH-VZn catalyst just requires a low overpotential of 240 mV at 10 mA cm-2, and delivers robust stability at 50 mA cm-2 over 120 h, which outperforms the benchmark of noble metal RuO2 and those of most non-noble metal catalysts ever reported.
ABSTRACT
BACKGROUND: The enhanced recovery after cardiac surgery is a bundle of measurements from preoperative to postoperative phases to improve patients' recovery. METHODS: This study is a multicenter, stepwise design, cluster randomized controlled trial. About 3,600 patients presenting during control and intervention periods are eligible if they are aged from 18 to 80 years old awaiting elective cardiac surgery with cardiopulmonary bypass (CPB). About 5 centers are randomly assigned to staggered start dates for one-way crossover from the control phase to the intervention phase. In the intervention periods, patients will receive ERAS strategy including preoperative, intraoperative, and postoperative approaches. During the control phase, patients receive usual care. The primary outcome consists of major adverse cardiac and cerebrovascular events (MACCEs), postoperative pulmonary complications (PPCs), and acute kidney injury (AKI). DISCUSSION: This study aims to compare the application of ERAS management protocol and traditional management protocol in adult cardiac surgery under extracorporeal circulation.
Subject(s)
Cardiac Surgical Procedures , Adolescent , Adult , Aged , Aged, 80 and over , Female , Humans , Male , Middle Aged , Young Adult , Acute Kidney Injury/etiology , Acute Kidney Injury/epidemiology , Cardiac Surgical Procedures/methods , Cardiac Surgical Procedures/adverse effects , Cardiopulmonary Bypass/methods , Enhanced Recovery After Surgery , Postoperative Complications/epidemiology , Postoperative Complications/prevention & controlABSTRACT
Renal fibrosis is a pathologic process that leads to irreversible renal failure without effective treatment. Epithelial-to-mesenchymal transition (EMT) plays a key role in this process. The current study found that aberrant expression of IL-11 is critically involved in tubular EMT. IL-11 and its receptor subunit alpha-1 (IL-11Rα1) were significantly induced in renal tubular epithelial cells (RTECs) in unilateral ureteral obstruction (UUO) kidneys, co-localized with transforming growth factor-ß1. IL-11 knockdown ameliorated UUO-induced renal fibrosis in vivo and transforming growth factor-ß1-induced EMT in vitro. IL-11 intervention directly induced the transdifferentiation of RTECs to the mesenchymal phenotype and increased the synthesis of profibrotic mediators. The EMT response induced by IL-11 was dependent on the sequential activation of STAT3 and extracellular signal-regulated kinase 1/2 signaling pathways and the up-regulation of metadherin in RTECs. Micheliolide (MCL) competitively inhibited the binding of IL-11 with IL-11Rα1, suppressing the activation of STAT3 and extracellular signal-regulated kinase 1/2-metadherin pathways, ultimately inhibiting renal tubular EMT and interstitial fibrosis induced by IL-11. In addition, treatment with dimethylaminomicheliolide, a pro-drug of MCL for in vivo use, significantly ameliorated renal fibrosis exacerbated by IL-11 in the UUO model. These findings suggest that IL-11 is a promising target in renal fibrosis and that MCL/dimethylaminomicheliolide exerts its antifibrotic effect by suppressing IL-11/IL-11Rα1 interaction and blocking its downstream effects.
Subject(s)
Epithelial-Mesenchymal Transition , Kidney Diseases , Ureteral Obstruction , Epithelial-Mesenchymal Transition/drug effects , Fibrosis , Interleukin-11/metabolism , Interleukin-11/pharmacology , Interleukin-11/therapeutic use , Kidney/pathology , Kidney Diseases/chemically induced , Kidney Diseases/prevention & control , Kidney Diseases/metabolism , Mitogen-Activated Protein Kinase 3/metabolism , Mitogen-Activated Protein Kinase 3/pharmacology , Transcription Factors/metabolism , Transforming Growth Factor beta1/metabolism , Ureteral Obstruction/drug therapy , Ureteral Obstruction/metabolism , Ureteral Obstruction/pathology , Animals , MiceABSTRACT
Imaging flow cytometry, which combines the advantages of flow cytometry and microscopy, has emerged as a powerful tool for cell analysis in various biomedical fields such as cancer detection. In this study, we develop multiplex imaging flow cytometry (mIFC) by employing a spatial wavelength division multiplexing technique. Our mIFC can simultaneously obtain brightfield and multi-color fluorescence images of individual cells in flow, which are excited by a metal halide lamp and measured by a single detector. Statistical analysis results of multiplex imaging experiments with resolution test lens, magnification test lens, and fluorescent microspheres validate the operation of the mIFC with good imaging channel consistency and micron-scale differentiation capabilities. A deep learning method is designed for multiplex image processing that consists of three deep learning networks (U-net, very deep super resolution, and visual geometry group 19). It is demonstrated that the cluster of differentiation 24 (CD24) imaging channel is more sensitive than the brightfield, nucleus, or cancer antigen 125 (CA125) imaging channel in classifying the three types of ovarian cell lines (IOSE80 normal cell, A2780, and OVCAR3 cancer cells). An average accuracy rate of 97.1% is achieved for the classification of these three types of cells by deep learning analysis when all four imaging channels are considered. Our single-detector mIFC is promising for the development of future imaging flow cytometers and for the automatic single-cell analysis with deep learning in various biomedical fields.
ABSTRACT
Ferroptosis is a novel type of cell death that is dependent on lipid peroxidation and iron accumulation, which distinguishes it from other types of programmed cell death. Current research indicates a significant association between ferroptosis and various pathological conditions, including cancer, neurological disorders, and cardiovascular diseases, albeit with a relatively unexplored role in dermatological afflictions. This paper elaborates on the mechanisms and signalling pathways of ferroptosis, summarizing the recent studies on ferroptosis and its related factors in dermatosis. Our objective is to shed light on novel perspectives and therapeutic strategies for dermatosis, enhancing the understanding of this under-researched area through this comprehensive review.
Subject(s)
Ferroptosis , Signal Transduction , Skin Diseases , Ferroptosis/physiology , Humans , Skin Diseases/metabolism , Iron/metabolism , Lipid Peroxidation , AnimalsABSTRACT
Objective: To assess global, regional and national trends in the impact of floods from 1990 to 2022 and determine factors influencing flood-related deaths. Methods: We used data on flood disasters from the International Disaster Database for 1990-2022 from 168 countries. We calculated the annual percentage change to estimate trends in the rates of people affected and killed by floods by study period, World Health Organization (WHO) region, country income level and flood type. We used multivariable logistic regression analysis to assess the factors associated with death from floods. Findings: From 1990 to 2022, 4713 floods were recorded in 168 countries, which affected > 3.2 billion people, caused 218 353 deaths and were responsible for more than 1.3 trillion United States dollars of economic losses. The WHO Western Pacific Region had the most people affected by floods (> 2.0 billion), accounting for 63.19% (2 024 599 380/3 203 944 965) of all affected populations. The South-East Asia Region had the most deaths (71 713, 32.84%). The African and Eastern Mediterranean Regions had the highest number of people affected and killed by floods per 100 000 population in 2022. The odds of floods causing more than 50 deaths were significantly higher in low-income countries (adjusted odds ratio: 14.34; 95% confidence interval: 7.46 to 30.04) compared with high-income countries. Numbers of people affected and mortality due to floods declined over time. Conclusion: Despite the decreases in populations affected and deaths, floods still have a serious impact on people and economies globally, particularly in lower-income countries. Action is needed to improve disaster risk management and flood mitigation.
Subject(s)
Floods , Humans , Global Health , Disasters , Developing Countries , Logistic Models , Natural DisastersABSTRACT
BACKGROUND: This study aimed to identify metabolic subtypes in ESCA, explore their relationship with immune landscapes, and establish a metabolic index for accurate prognosis assessment. METHODS: Clinical, SNP, and RNA-seq data were collected from 80 ESCA patients from the TCGA database and RNA-seq data from the GSE19417 dataset. Metabolic genes associated with overall survival (OS) and progression-free survival (PFS) were selected, and k-means clustering was performed. Immune-related pathways, immune infiltration, and response to immunotherapy were predicted using bioinformatic algorithms. Weighted gene co-expression network analysis (WGCNA) was conducted to identify metabolic genes associated with co-expression modules. Lastly, cell culture and functional analysis were performed using patient tissue samples and ESCA cell lines to verify the identified genes and their roles. RESULTS: Molecular subtypes were identified based on the expression profiles of metabolic genes, and univariate survival analysis revealed 163 metabolic genes associated with ESCA prognosis. Consensus clustering analysis classified ESCA samples into three distinct subtypes, with MC1 showing the poorest prognosis and MC3 having the best prognosis. The subtypes also exhibited significant differences in immune cell infiltration, with MC3 showing the highest scores. Additionally, the MC3 subtype demonstrated the poorest response to immunotherapy, while the MC1 subtype was the most sensitive. WGCNA analysis identified gene modules associated with the metabolic index, with SLC5A1, NT5DC4, and MTHFD2 emerging as prognostic markers. Gene and protein expression analysis validated the upregulation of MTHFD2 in ESCA. MTHFD2 promotes the progression of ESCA and may be a potential therapeutic target for ESCA. CONCLUSION: The established metabolic index and identified metabolic genes offer potential for prognostic assessment and personalized therapeutic interventions for ESCA, underscoring the importance of targeting metabolism-immune interactions in ESCA. MTHFD2 promotes the progression of ESCA and may be a potential therapeutic target for ESCA.