Integrated analysis of multimodal single-cell data with structural similarity

ARID5A orchestrates cardiac aging and inflammation through MAVS mRNA stabilization

Elucidating the regulatory mechanisms of human cardiac aging remains a great challenge. Here, using human heart tissues from 74 individuals ranging from young (≤35 years) to old (≥65 years), we provide an overview of the histological, cellular and molecular alterations underpinning the aging of human hearts. We decoded aging-related gene expression changes at single-cell resolution and identified increased inflammation as the key event, driven by upregulation of ARID5A, an RNA-binding protein. ARID5A epi-transcriptionally regulated Mitochondrial Antiviral Signaling Protein (MAVS) mRNA stability, leading to NF-κB and TBK1 activation, amplifying aging and inflammation phenotypes. The application of gene therapy using lentiviral vectors encoding shRNA targeting ARID5A into the myocardium not only mitigated the inflammatory and aging phenotypes but also bolstered cardiac function in aged mice. Altogether, our study provides a valuable resource and advances our understanding of cardiac aging mechanisms by deciphering the ARID5A-MAVS axis in post-transcriptional regulation.

Identification of novel metabolism-related biomarkers of Kawasaki disease by integrating single-cell RNA sequencing analysis and machine learning algorithms

Background

The bile acid metabolism (BAM) and fatty acid metabolism (FAM) have been implicated in Kawasaki disease (KD), but their precise mechanisms remain unclear. Identifying signature cells and genes related to BAM and FAM could offer a deeper understanding of their role in the pathogenesis of KD.

Method

We analyzed the public single-cell RNA sequencing (scRNA-seq) dataset GSE1687323 to characterize the immune cell-type landscape in KD. Gene sets related to BAM and FAM were collected from the Gene Set Enrichment Analysis (GSEA) database and previous literature. We analyzed the cellular heterogeneity of BAM and FAM at the single-cell level using R packages. Through differential expressed genes (DEG) analysis, high-dimensional Weighted Correlation Network Analysis (hdWGCNA) and machine learning algorithms, we identified signature genes associated with both BAM and FAM. The cellular expression patterns of signature genes were further validated using our own scRNA-seq dataset. Finally, quantitative real-time PCR (qRT-PCR) was performed to validate the expression levels of signature genes in KD, and Receiver Operating Characteristic (ROC) curve analysis was conducted to evaluate their diagnostic potential.

Results

Enhanced BAM and FAM were detected in monocytes and natural killer (NK) cells from KD in the public scRNA-seq dataset. Our scRNA-seq data confirmed the signature genes identified by machine learning algorithms: Vimentin (VIM) and chloride intracellular channel 1 (CLIC1) were upregulated in monocytes, while integrin subunit beta 2 (ITGB2) was elevated in NK cells of KD. qRT-PCR results also validated the bioinformatic analysis. Moreover, these genes demonstrated significant diagnostic potential. In the training dataset (GSE68004), the area under the curve (AUC) values and 95% CI were as follows: VIM: 0.914 (0.863-0.966), ITGB2: 0.958 (0.925-0.991), and CLIC1: 0.985 (0.969-1). The validation dataset (GSE73461) yielded similarly robust results, with AUC values and 95% CI: VIM: 0.872 (0.811-0.934), ITGB2: 0.861 (0.795-0.928), and CLIC1: 0.893 (0.837-0.948).

Conclusion

This study successfully identified and validated VIM and CLIC1 in monocytes, as well as ITGB2 in NK cells, as novel metabolism-related genes in KD. These findings suggest that BAM and FAM may play crucial roles in KD pathogenesis. Furthermore, these signature genes hold promising potential as diagnostic biomarkers for KD.

Cancer-associated fibroblasts gene signature: a novel approach to survival prediction and immunotherapy guidance in colon cancer

Background

Fibroblasts can regulate tumour development by secreting various factors. For COAD survival prediction and CAFs-based treatment recommendations, it is critical to comprehend the heterogeneity of CAFs and find biomarkers.

Methods

We identified fibroblast-associated specific marker genes in colon adenocarcinoma by single-cell sequencing analysis. A fibroblasts-related gene signature was developed, and colon adenocarcinoma patients were classified into high-risk and low-risk cohorts based on the median risk score. Additionally, the impact of these risk categories on the tumor microenvironment was evaluated. The ability of CAFGs signature to assess prognosis and guide treatment was validated using external cohorts. Ultimately, we verified MAN1B1 expression and function through in vitro assays.

Results

Relying on the bulk RNA-seq and scRNA-seq data study, we created a predictive profile with 11 CAFGs. The profile effectively differentiated survival differences among cohorts of colon adenocarcinoma patients. The nomogram further effectively predicted the prognosis of COAD patients, with low-risk patients having a better prognosis. A higher immune infiltration rate and lower IC50 values of anticancer drugs were significant in the high-risk group. In cellular experiments, Following MAN1B1 knockdown, in cell assays, the colony formation, migration, and invasion ability of HCT116 and HT29 cell lines decreased.

Conclusion

Our CAFG signature provides important insights into the role of CAF cells in influencing COAD prognosis. It may also serve as a guide for selecting immunotherapy options and predicting chemotherapy responses in COAD patients.

Integration of single-cell and bulk RNA sequencing data using machine learning identifies oxidative stress-related genes LUM and PCOLCE2 as potential biomarkers for heart failure

Oxidative stress (OS) is a pivotal mechanism driving the progression of cardiovascular diseases, particularly heart failure (HF). However, the comprehensive characterisation of OS-related genes in HF remains largely unexplored. In the present study, we analysed single-cell RNA sequencing datasets from the Gene Expression Omnibus and OS gene sets from GeneCards. We identified 167 OS-related genes potentially linked to HF by applying algorithms, such as AUCell, UCell, singscore, ssgsea, and AddModuleScore, combined with correlation analysis. Subsequently, we used feature selection algorithms, including least absolute shrinkage and selection operator, XGBoost, Boruta, random forest, gradient boosting machines, decision trees, and support vector machine recursive feature elimination, to identify lumican (LUM) and procollagen C-endopeptidase enhancer 2 (PCOLCE2) as key biomarker candidates with significant diagnostic potential. Bulk RNA-sequencing confirmed their elevated expression in patients with HF, highlighting their predictive utility. Single-cell analysis further revealed their upregulation primarily in fibroblasts, emphasising their cell-specific role in HF. To validate these findings, we developed a transverse aortic constriction-induced HF mouse model that showed enhanced cardiac OS activity and significant PCOLCE2 upregulation in the HF group. These results provide strong evidence of the involvement of OS-related mechanisms in HF. Herein, we propose a diagnostic strategy that provides novel insights into the molecular mechanisms underlying HF. However, further studies are required to validate its clinical utility and ensure its application in the diagnosis of HF.

Differential Phagocytosis induces Diverse Macrophage Activation States in Malignant Gliomas

Diffuse midline glioma (DMG) and Glioblastoma are malignant brain tumors in pediatric and adult patients. The current standard-of-care treatment for DMG is radiotherapy (RT), whereas GBM treatment includes surgery, followed by RT and chemotherapy. Although RT is known to modulate immune responses in cancer and enhance the effectiveness of myeloid checkpoint blockade, the downstream macrophage responses to differential phagocytosis induction remain poorly understood. This study examined macrophage-mediated phagocytosis caused by either RT, anti-CD47 checkpoint blockade, or their combination. We found that RT increased the expression of several damage-associated molecular patterns on the surface of glioma cell lines. Furthermore, RT enhanced anti-CD47-mediated macrophage phagocytosis of glioma cell lines in vitro . Single-cell RNA-sequencing revealed the diverse transcriptional and functional signatures of human macrophage subsets that either promoted or inhibited phagocytosis of glioma cells pretreated with RT, anti-CD47 therapy, or both. Consistent with these results, the combination therapy significantly reduced tumor growth, prolonged survival in glioma-bearing mice, and induced distinct macrophage activation states in vivo compared to either treatment alone. These findings highlight the plasticity and heterogeneity of macrophage responses during phagocytosis and provide compelling evidence for combining RT with anti-CD47 therapy as a promising therapeutic strategy for glioma treatment.

Novel post-translational modification learning signature reveals B4GALT2 as an immune exclusion regulator in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) presents significant challenges in prognosis and treatment efficacy evaluation. While post-translational modifications are known to influence tumor progression, their prognostic value in LUAD remains largely unexplored.

METHODS

We developed a post-translational modification learning signature (PTMLS) using machine learning techniques, analyzing data from 1231 LUAD patients across seven global cohorts. The signature's efficacy in predicting immunotherapy response was evaluated using 12 immunotherapy cohorts spanning multiple cancer types (n=1201). An in-house LUAD tissue cohort (n=171) was used to validate beta-1,4-galactosyltransferase 2's (B4GALT2's) prognostic significance. The role of B4GALT2 in immune exclusion was investigated through in vivo and in vitro experiments.

RESULTS

The established PTMLS exhibited exceptional predictive capabilities in LUAD patient outcomes, surpassing the efficacy of 98 existing LUAD prognostic indicators. The system's predictive value was validated across diverse malignancy categories for immunotherapeutic response assessment. From a biological perspective, significant correlations were observed between PTMLS and immunological parameters, whereby elevated PTMLS levels were characterized by attenuated immune responses and immunologically cold neoplastic features. Within the PTMLS framework, B4GALT2 was identified as a crucial molecular component (r=0.82, p<0.05), and its heightened expression was linked to unfavorable clinical outcomes in LUAD cases, particularly in specimens exhibiting CD8-depleted phenotypes. The spatial distribution patterns between B4GALT2 and immune cell populations, specifically CD8+ T lymphocytes and CD20+ B lymphocytes, were elucidated through multiplexed immunofluorescence analysis. Laboratory investigations subsequently established B4GALT2's regulatory influence on LUAD cellular expansion in both laboratory cultures and animal models. Significantly, suppression of B4GALT2 was found to enhance CD8+ T lymphocyte populations and their functional status, thereby potentiating anti-programmed cell death protein 1 immunotherapeutic efficacy in animal studies. This phenomenon was characterized by reduced CD62L+CD8 T lymphocyte levels alongside elevated GZMB+/CD44+/CD69+CD8 T cell populations.

CONCLUSION

The developed PTMLS system represents an effective instrument for individualized prognostic evaluation and immunotherapy stratification in both LUAD and diverse cancer populations. The identification of B4GALT2 as a previously unrecognized oncogenic factor involved in immune exclusion presents a novel therapeutic avenue for LUAD treatment and immunotherapy optimization.

PANoptosis-Related Optimal Model (PROM): A Novel Prognostic Tool Unveiling Immune Dynamics in Lung Adenocarcinoma

Background: PANoptosis, a recently characterized inflammatory programmed cell death modality orchestrated by the PANoptosome complex, integrates molecular mechanisms of pyroptosis, apoptosis, and necroptosis. Although this pathway potentially mediates tumor progression, its role in lung adenocarcinoma (LUAD) remains largely unexplored. Methods: Through comprehensive single-cell transcriptomic profiling, we systematically identified critical PANoptosis-associated gene signatures. Prognostic molecular determinants were subsequently delineated via univariate Cox proportional hazards regression analysis. We constructed a PANoptosis-related optimal model (PROM) through the integration of 10 machine learning algorithms. The model was initially developed using The Cancer Genome Atlas (TCGA)-LUAD cohort and subsequently validated across six independent LUAD cohorts. Model performance was evaluated using mean concordance index. Furthermore, we conducted extensive multiomics analyses to delineate differential pathway activation patterns and immune cell infiltration profiles between PROM-stratified risk subgroups. Results: Cellular populations exhibiting elevated PANoptosis signatures demonstrated enhanced intercellular signaling networks. PROM demonstrated superior prognostic capability across multiple validation cohorts. Receiver operating characteristic curve analyses revealed area under the curve values exceeding 0.7 across all seven cohorts, with several achieving values above 0.8, indicating robust discriminative performance. The model score exhibited significant correlation with immunological parameters. Notably, high PROM scores were associated with attenuated immune responses, suggesting an immunosuppressive tumor microenvironment. Multiomics investigations revealed significant alterations in critical oncogenic pathways and immune landscape between PROM-stratified subgroups. Conclusion: This investigation establishes PROM as a clinically applicable prognostic tool for LUAD risk stratification. Beyond its predictive utility, PROM elucidates PANoptosis-associated immunological and biological mechanisms underlying LUAD progression. These findings provide novel mechanistic insights into LUAD pathogenesis and may inform the development of targeted therapeutic interventions and personalized treatment strategies to optimize patient outcomes.

Single-cell analysis reveals cellular heterogeneity, gene expression profiles, and pathway dynamics in acne vulgaris

Acne vulgaris is a common dermatological condition, particularly in adolescents, and is associated with significant physical and psychological impacts. Its pathogenesis involves genetic factors, sebaceous gland dysfunction, inflammatory responses, and alterations in skin microbiota. Despite advancements in treatment, a comprehensive understanding of its cellular and molecular mechanisms remains limited, making therapeutic strategies challenging. Single-cell RNA sequencing (scRNA-seq) data from six acne vulgaris patients were obtained from the GEO database (GSE175817), filtered, and integrated. Differentially expressed genes (DEGs) were identified and subjected to functional enrichment analysis to explore the underlying pathogenesis. Gene set variation analysis (GSVA) was used to investigate metabolic pathway alterations, and ligand-receptor interactions were analyzed to examine cell-to-cell communication between lesional and non-lesional skin tissues. The integration of scRNA-seq data yielded the identification of eight distinct cell clusters, including endothelial cells, myeloid cells, lymphocytes, melanocytes, sebaceous gland cells, smooth muscle cells, keratinocytes and fibroblasts. The proportions of lymphoid and myeloid cells were found to be significantly different between the lesional and non-lesional sites. The differential expression of genes (DEGs) was found to be significantly specific to the different cell clusters. Abnormal intercellular communication was found to result in a substantial increase in the number and intensity of communications in the area of acne vulgaris lesions. Moreover, specific ligand-receptor pairs for SPP1 and IL6 associated with acne vulgaris were identified. Furthermore, the presence of specific alterations in metabolic pathways, including riboflavin metabolism, niacin metabolism, and lipoic acid metabolism, was observed. Our findings demonstrate the cellular heterogeneity and dysfunction of intercellular communication and metabolic signaling in the lesional skin tissues of patients with acne vulgaris. These findings have important implications for understanding the complex biological processes of acne vulgaris.

Comprehensive Analysis of m6A-Related Programmed Cell Death Genes Unveils a Novel Prognostic Model for Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) involves complex dysregulated cellular processes, including programmed cell death (PCD), influenced by N6-methyladenosine (m6A) RNA modification. This study integrates bulk RNA and single-cell sequencing data to identify 43 prognostically valuable m6A-related PCD genes, forming the basis of a 13-gene risk model (m6A-related PCD signature [mPCDS]) developed using machine-learning algorithms, including CoxBoost and SuperPC. The mPCDS demonstrated significant predictive performance across multiple validation datasets. In addition to its prognostic accuracy, mPCDS revealed distinct genomic profiles, pathway activations, associations with the tumour microenvironment and potential for predicting drug sensitivity. Experimental validation identified RCN1 as a potential oncogene driving LUAD progression and a promising therapeutic target. The mPCDS offers a new approach for LUAD risk stratification and personalised treatment strategies.

Exploring cellular diversity in lung adenocarcinoma epithelium: Advancing prognostic methods and immunotherapeutic strategies

Immunotherapy has brought significant advancements in the treatment of lung adenocarcinoma (LUAD), but identifying suitable candidates remains challenging. In this study, we investigated tumour cell heterogeneity using extensive single-cell data and explored the impact of different tumour cell cluster abundances on immunotherapy in the POPLAR and OAK immunotherapy cohorts. Notably, we found a significant correlation between CKS1B+ tumour cell abundance and treatment response, as well as stemness potential. Leveraging marker genes from the CKS1B+ tumour cell cluster, we employed machine learning algorithms to establish a prognostic and immunotherapeutic signature (PIS) for LUAD. In multiple cohorts, PIS outperformed 144 previously published signatures in predicting LUAD prognosis. Importantly, PIS reliably predicted genomic alterations, chemotherapy sensitivity and immunotherapy responses. Immunohistochemistry validated lower expression of immune markers in the low-PIS group, while in vitro experiments underscored the role of the key gene PSMB7 in LUAD progression. In conclusion, PIS represents a novel biomarker facilitating the selection of suitable LUAD patients for immunotherapy, ultimately improving prognosis and guiding clinical decisions.

Single-cell and bulk transcriptomic datasets enable the development of prognostic models based on dynamic changes in the tumor immune microenvironment in patients with hepatocellular carcinoma and portal vein tumor thrombus

Background

Hepatocellular carcinoma (HCC) patients exhibiting portal vein tumor thrombosis (PVTT) face a high risk of rapid malignant progression and poor outcomes, with this issue being compounded by a lack of effective treatment options. The integration of bulk RNA-sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq) datasets focused on samples from HCC patients with PVTT has the potential to yield unprecedented insight into the dynamic changes in the tumor microenvironment (TME) and associated immunological characteristics in these patients, providing an invaluable tool for the reliable prediction of disease progression and treatment responses.

Methods

scRNA-seq data from both primary tumor (PT) and PVTT cells were downloaded from the Gene Expression Omnibus (GEO) database, while the International Cancer Genome Consortium (ICGC) and Cancer Genome Atlas (TCGA) databases were used to access bulk RNA-seq datasets. scRNA-seq, clustering, GSVA enrichment, mutational profiling, and predictive immunotherapeutic treatment analyses were conducted using these data with the goal of systematically assessing the heterogeneity of PT and PVTT cells and establishing a model capable of predicting immunotherapeutic and prognostic outcomes in patients with HCC.

Results

These analyses revealed that PVTT cells exhibited patterns of tumor proliferation, stromal activation, and low levels of immune cell infiltration, presenting with immune desert and immune rejection-like phenotypes. PT cells, in contrast, were found to exhibit a pattern of immunoinflammatory activity. Core PVTT-associated genes were clustered into three patterns consistent with the tumor immune rejection and immune desert phenotypes. An established clustering model was capable of predicting tumor inflammatory stage, subtype, TME stromal activity, and patient outcomes. PVTT signature genes were further used to establish a risk model, with the risk scores derived from this model providing a tool to evaluate patient clinicopathological features including clinical stage, tumor differentiation, histological subtype, microsatellite instability status, and tumor mutational burden. These risk scores were also able to serve as an independent predictor of patient survival outcomes, responses to adjuvant chemotherapy, and responses to immunotherapy. In vitro experiments were used to partially validate the biological prediction results.

Conclusion

These results offer new insight into the biological and immunological landscape of PVTT in HCC patients, By utilizing individual patient risk scores, providing an opportunity to guide more effective immunotherapeutic interventional efforts.

Leveraging mitochondrial-programmed cell death dynamics to enhance prognostic accuracy and immunotherapy efficacy in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is a highly heterogeneous disease, posing significant challenges to accurate prognosis prediction. Mitochondria play a central role in the energy metabolism of eukaryotic cells and can influence programmed cell death (PCD) mechanisms, which are critical in tumorigenesis and cancer progression. However, the prognostic significance of the interplay between mitochondrial function and PCD in LUAD requires further investigation.

METHODS

We analyzed data from 1231 LUAD patients across seven global cohorts to develop a mitochondrial-related PCD signature (MPCDS) using machine learning. Validation was done using six immunotherapy cohorts (LUAD, melanoma, clear cell renal cell carcinoma; n=935) and a pan-cancer cohort of 21 tumor types. An in-house LUAD tissue cohort (n=100) confirmed the prognostic significance of nucleoside diphosphate kinase 4 (NME4). In vivo and in vitro experiments explored NME4's role in immune exclusion.

RESULTS

The MPCDS demonstrated strong predictive performance for prognosis in LUAD patients, surpassing 114 previously published LUAD signatures. Additionally, MPCDS effectively predicted outcomes in immunotherapy patients (including those with LUAD, melanoma, and clear cell renal cell carcinoma). Biologically, MPCDS was significantly associated with immune features, with the high MPCDS group exhibiting reduced immune activity and a tendency towards cold tumors. NME4, a key gene within the MPCDS (correlation=0.55, p<0.05), was associated with poorer prognosis in LUAD patients with high expression, particularly in CD8 desert phenotypes, as validated by our in-house cohort. Multiplex immunofluorescence confirmed the spatial colocalization and exclusion relationship between NME4 and immune cells such as CD3+ T cells and CD20+ B cells. Further experiments revealed that NME4 regulated the proliferation and invasion of LUAD cells both in vitro and in vivo. Importantly, inhibiting NME4 increased the abundance and activity of CD8+ T cells and enhanced the antitumor immunity of anti-programmed cell death protein-1 therapy in vivo.

CONCLUSION

The MPCDS provides personalized risk assessment and immunotherapy interventions for individual LUAD patients. NME4, a key gene within the MPCDS, has been identified as a novel oncogene associated with immune exclusion and may serve as a new target for LUAD intervention and immunotherapy.

Unraveling the role of low-density lipoprotein-related genes in lung adenocarcinoma: Insights into tumor microenvironment and clinical prognosis

BACKGROUND

The hypothesized link between low-density lipoprotein (LDL) and oncogenesis has garnered significant interest, yet its explicit impact on lung adenocarcinoma (LUAD) remains to be elucidated. This investigation aims to demystify the function of LDL-related genes (LRGs) within LUAD, endeavoring to shed light on the complex interplay between LDL and carcinogenesis.

METHODS

Leveraging single-cell transcriptomics, we examined the role of LRGs within the tumor microenvironment (TME). The expression patterns of LRGs across diverse cellular phenotypes were delineated using an array of computational methodologies, including AUCell, UCell, singscore, ssGSEA, and AddModuleScore. CellChat facilitated the exploration of distinct cellular interactions within LDL_low and LDL_high groups. The findmarker utility, coupled with Pearson correlation analysis, facilitated the identification of pivotal genes correlated with LDL indices. An integrative approach to transcriptomic data analysis was adopted, utilizing a machine learning framework to devise an LDL-associated signature (LAS). This enabled the delineation of genomic disparities, pathway enrichments, immune cell dynamics, and pharmacological sensitivities between LAS stratifications.

RESULTS

Enhanced cellular crosstalk was observed in the LDL_high group, with the CoxBoost+Ridge algorithm achieving the apex c-index for LAS formulation. Benchmarking against 144 extant LUAD models underscored the superior prognostic acuity of LAS. Elevated LAS indices were synonymous with adverse outcomes, diminished immune surveillance, and an upsurge in pathways conducive to neoplastic proliferation. Notably, a pronounced susceptibility to paclitaxel and gemcitabine was discerned within the high-LAS cohort, delineating prospective therapeutic corridors.

CONCLUSION

This study elucidates the significance of LRGs within the TME and introduces an LAS for prognostication in LUAD patients. Our findings accentuate putative therapeutic targets and elucidate the clinical ramifications of LAS deployment.

ScRNAbox: empowering single-cell RNA sequencing on high performance computing systems

BACKGROUND

Single-cell RNA sequencing (scRNAseq) offers powerful insights, but the surge in sample sizes demands more computational power than local workstations can provide. Consequently, high-performance computing (HPC) systems have become imperative. Existing web apps designed to analyze scRNAseq data lack scalability and integration capabilities, while analysis packages demand coding expertise, hindering accessibility.

RESULTS

In response, we introduce scRNAbox, an innovative scRNAseq analysis pipeline meticulously crafted for HPC systems. This end-to-end solution, executed via the SLURM workload manager, efficiently processes raw data from standard and Hashtag samples. It incorporates quality control filtering, sample integration, clustering, cluster annotation tools, and facilitates cell type-specific differential gene expression analysis between two groups. We demonstrate the application of scRNAbox by analyzing two publicly available datasets.

CONCLUSION

ScRNAbox is a comprehensive end-to-end pipeline designed to streamline the processing and analysis of scRNAseq data. By responding to the pressing demand for a user-friendly, HPC solution, scRNAbox bridges the gap between the growing computational demands of scRNAseq analysis and the coding expertise required to meet them.

scACT: Accurate Cross-modality Translation via Cycle-consistent Training from Unpaired Single-cell Data

Single-cell sequencing technologies have revolutionized genomics by enabling the simultaneous profiling of various molecular modalities within individual cells. Their integration, especially cross-modality translation, offers deep insights into cellular regulatory mechanisms. Many methods have been developed for cross-modality translation, but their reliance on scarce high-quality co-assay data limits their applicability. Addressing this, we introduce scACT, a deep generative model designed to extract cross-modality biological insights from unpaired single-cell data. scACT tackles three major challenges: aligning unpaired multi-modal data via adversarial training, facilitating cross-modality translation without prior knowledge via cycle-consistent training, and enabling interpretable regulatory interconnections explorations via in-silico perturbations. To test its performance, we applied scACT on diverse single-cell datasets and found it outperformed existing methods in all three tasks. Finally, we have developed scACT as an individual open-source software package to advance single-cell omics data processing and analysis within the research community.

Machine learning algorithms integrate bulk and single-cell RNA data to unveil oxidative stress following intracerebral hemorrhage

BACKGROUND

Increased oxidative stress (OS) activity following intracerebral hemorrhage (ICH) had significantly impacting patient prognosis. Identifying optimal genes associated with OS could enhance the understanding of OS after ICH.

METHODS

We employed single-cell RNA sequencing (scRNA-seq) to investigate the heterogeneity of OS across various cellular tiers following ICH, aiming to acquire biological insights into ICH. We utilized AUCell, Ucell, singscore, ssgsea, and AddModuleScore algorithms, along with correlation analysis, to identify hub genes influencing high OS post-ICH. Furthermore, we employed four machine learning algorithms, eXtreme Gradient Boosting, Boruta, Random Forest, and Least Absolute Shrinkage and Selection Operator, to identify the optimal feature genes. To validate the accuracy of our analysis, we conducted validation in ICH animal experiments.

RESULTS

After analyzing the scRNA-seq dataset using various algorithms, we found that OS activity exhibited heterogeneity across different cellular layers following ICH, with particularly heightened activity observed in monocytes. Further integration of bulk data and machine learning algorithms revealed that ANXA2 and COTL1 were closely associated with high OS after ICH. Our animal experiments demonstrated an increase in OS expression post-ICH. Additionally, the protein expression of ANXA2 and COTL1 was significantly elevated and co-localized with microglia. Pearson correlation coefficient analysis revealed a significant correlation between ANXA2 and OS, indicating strong consistency (r = 0.84, p < 0.05). Similar results were observed for COTL1 and OS (r = 0.69, p < 0.05).

CONCLUSIONS

Following ICH, ANXA2 and COTL1 might penetrate the brain via monocytes, localize within microglia, and enhance OS activity. This might help us better understand OS after ICH.

scCross: a deep generative model for unifying single-cell multi-omics with seamless integration, cross-modal generation, and in silico exploration

Single-cell multi-omics data reveal complex cellular states, providing significant insights into cellular dynamics and disease. Yet, integration of multi-omics data presents challenges. Some modalities have not reached the robustness or clarity of established transcriptomics. Coupled with data scarcity for less established modalities and integration intricacies, these challenges limit our ability to maximize single-cell omics benefits. We introduce scCross, a tool leveraging variational autoencoders, generative adversarial networks, and the mutual nearest neighbors (MNN) technique for modality alignment. By enabling single-cell cross-modal data generation, multi-omics data simulation, and in silico cellular perturbations, scCross enhances the utility of single-cell multi-omics studies.

Deciphering the prognostic role of endoplasmic reticulum stress in lung adenocarcinoma: integrating prognostic prediction and immunotherapy strategies

Endoplasmic reticulum stress (ERS) is a critical factor influencing lung adenocarcinoma (LUAD) progression and patient outcomes. In this study, we analyzed gene expression data from LUAD samples sourced from The Cancer Genomic Atlas and Gene Expression Omnibus databases. Utilizing advanced statistical methods including LASSO and Cox regression, we developed a ERS-associated signature (ERAS) based on ten ERS-related genes. This model stratified patients into high- and low-risk groups, with the high-risk group exhibiting decreased survival rates, elevated tumor mutational burden, and heightened chemotherapy sensitivity. Additionally, we observed lower immune and ESTIMATE scores in the high-ERAS group, indicating a potentially compromised immune response. Experimental validation through quantitative real-time polymerase chain reaction confirmed the utility of our model. Furthermore, we constructed a nomogram to predict 1-, 3-, and 5-year survival rates, providing clinicians with a valuable tool for personalized patient management. In conclusion, our study demonstrates the efficacy of the ERAS in identifying high-ERAS LUAD patients, offering promising implications for improved prognostication and treatment strategies.

Reparative effects after low-dose radiation exposure: Inhibition of atherosclerosis by reducing NETs release

OBJECTIVE

The cardiovascular system effects of environmental low-dose radiation exposure on radiation practitioners remain uncertain and require further investigation. The aim of this study was to initially investigate and explore the mechanisms by which low-dose radiation may contribute to atherosclerosis through a multi-omics joint comprehensive basic experiment.

METHODS

We used WGCNA and differential analyses to identify shared genes and potential pathways between radiation injury and atherosclerosis sequencing datasets, as well as tissue transcriptome immune infiltration level extrapolation and single-cell transcriptome data correction using the CIBERSORT deconvolution algorithm. Animal models were constructed by combining a high-fat diet with 5 Gy γ-ray whole-body low-dose ionizing radiation. The detection of NETs release was validated by enzyme-linked immunosorbent assay.

RESULTS

Analysis reveals shared genes in both datasets of post-irradiation and atherosclerosis, suggesting that immune system neutrophils may be a key node connecting radiation to atherosclerosis. NETs released by neutrophil death can influence the development of atherosclerosis. Animal experiments showed that the number of neutrophils decreased (P < 0.05) and the concentration of NETs reduced after low-dose radiation compared with the control group, and the concentration of NETs significantly increased (P < 0.05) in the HF group. Endothelial plaques were significantly increased in the high-fat feed group and significantly decreased in the low-dose radiation group compared with the control group.

CONCLUSIONS

Long-term low-dose ionizing radiation exposure stimulates neutrophils and inhibits their production of NETs, resulting in inhibition of atherosclerosis.

Unravelling tumour cell diversity and prognostic signatures in cutaneous melanoma through machine learning analysis

Melanoma, a highly malignant tumour, presents significant challenges due to its cellular heterogeneity, yet research on this aspect in cutaneous melanoma remains limited. In this study, we utilized single-cell data from 92,521 cells to explore the tumour cell landscape. Through clustering analysis, we identified six distinct cell clusters and investigated their differentiation and metabolic heterogeneity using multi-omics approaches. Notably, cytotrace analysis and pseudotime trajectories revealed distinct stages of tumour cell differentiation, which have implications for patient survival. By leveraging markers from these clusters, we developed a tumour cell-specific machine learning model (TCM). This model not only predicts patient outcomes and responses to immunotherapy, but also distinguishes between genomically stable and unstable tumours and identifies inflamed ('hot') versus non-inflamed ('cold') tumours. Intriguingly, the TCM score showed a strong association with TOMM40, which we experimentally validated as an oncogene promoting tumour proliferation, invasion and migration. Overall, our findings introduce a novel biomarker score that aids in selecting melanoma patients for improved prognoses and targeted immunotherapy, thereby guiding clinical treatment decisions.

Gene regulatory landscape of cerebral cortex folding

Folding of the cerebral cortex is a key aspect of mammalian brain development and evolution, and defects are linked to severe neurological disorders. Primary folding occurs in highly stereotyped patterns that are predefined in the cortical germinal zones by a transcriptomic protomap. The gene regulatory landscape governing the emergence of this folding protomap remains unknown. We characterized the spatiotemporal dynamics of gene expression and active epigenetic landscape (H3K27ac) across prospective folds and fissures in ferret. Our results show that the transcriptomic protomap begins to emerge at early embryonic stages, and it involves cell-fate signaling pathways. The H3K27ac landscape reveals developmental cell-fate restriction and engages known developmental regulators, including the transcription factor Cux2. Manipulating Cux2 expression in cortical progenitors changed their proliferation and the folding pattern in ferret, caused by selective transcriptional changes as revealed by single-cell RNA sequencing analyses. Our findings highlight the key relevance of epigenetic mechanisms in defining the patterns of cerebral cortex folding.

Deciphering lung adenocarcinoma evolution: Integrative single-cell genomics identifies the prognostic lung progression associated signature

We employed single-cell analysis techniques, specifically the inferCNV method, to dissect the complex progression of lung adenocarcinoma (LUAD) from adenocarcinoma in situ (AIS) through minimally invasive adenocarcinoma (MIA) to invasive adenocarcinoma (IAC). This approach enabled the identification of Cluster 6, which was significantly associated with LUAD progression. Our comprehensive analysis included intercellular interaction, transcription factor regulatory networks, trajectory analysis, and gene set variation analysis (GSVA), leading to the development of the lung progression associated signature (LPAS). Interestingly, we discovered that the LPAS not only accurately predicts the prognosis of LUAD patients but also forecasts genomic alterations, distinguishes between 'cold' and 'hot' tumours, and identifies potential candidates suitable for immunotherapy. PSMB1, identified within Cluster 6, was experimentally shown to significantly enhance cancer cell invasion and migration, highlighting the clinical relevance of LPAS in predicting LUAD progression and providing a potential target for therapeutic intervention. Our findings suggest that LPAS offers a novel biomarker for LUAD patient stratification, with significant implications for improving prognostic accuracy and guiding treatment decisions.

Spatiotemporal multi-omics: exploring molecular landscapes in aging and regenerative medicine

Aging and regeneration represent complex biological phenomena that have long captivated the scientific community. To fully comprehend these processes, it is essential to investigate molecular dynamics through a lens that encompasses both spatial and temporal dimensions. Conventional omics methodologies, such as genomics and transcriptomics, have been instrumental in identifying critical molecular facets of aging and regeneration. However, these methods are somewhat limited, constrained by their spatial resolution and their lack of capacity to dynamically represent tissue alterations. The advent of emerging spatiotemporal multi-omics approaches, encompassing transcriptomics, proteomics, metabolomics, and epigenomics, furnishes comprehensive insights into these intricate molecular dynamics. These sophisticated techniques facilitate accurate delineation of molecular patterns across an array of cells, tissues, and organs, thereby offering an in-depth understanding of the fundamental mechanisms at play. This review meticulously examines the significance of spatiotemporal multi-omics in the realms of aging and regeneration research. It underscores how these methodologies augment our comprehension of molecular dynamics, cellular interactions, and signaling pathways. Initially, the review delineates the foundational principles underpinning these methods, followed by an evaluation of their recent applications within the field. The review ultimately concludes by addressing the prevailing challenges and projecting future advancements in the field. Indubitably, spatiotemporal multi-omics are instrumental in deciphering the complexities inherent in aging and regeneration, thus charting a course toward potential therapeutic innovations.

Scalable integration of multiomic single-cell data using generative adversarial networks

MOTIVATION

Single-cell profiling has become a common practice to investigate the complexity of tissues, organs, and organisms. Recent technological advances are expanding our capabilities to profile various molecular layers beyond the transcriptome such as, but not limited to, the genome, the epigenome, and the proteome. Depending on the experimental procedure, these data can be obtained from separate assays or the very same cells. Yet, integration of more than two assays is currently not supported by the majority of the computational frameworks avaiable.

RESULTS

We here propose a Multi-Omic data integration framework based on Wasserstein Generative Adversarial Networks suitable for the analysis of paired or unpaired data with a high number of modalities (>2). At the core of our strategy is a single network trained on all modalities together, limiting the computational burden when many molecular layers are evaluated.

AVAILABILITY AND IMPLEMENTATION

Source code of our framework is available at https://github.com/vgiansanti/MOWGAN.

Clinical prognostication and immunotherapy response prediction in esophageal squamous cell carcinoma using the DNA damage repair-associated signature

BACKGROUND

The relationship between DNA damage repair (DDR) and cancer is intricately intertwined; however, its specific role in esophageal squamous cell carcinoma (ESCC) remains enigmatic.

METHODS

Employing single-cell analysis, we delineated the functionality of DDR-related genes within the tumor microenvironment (TME). A diverse array of scoring mechanisms, including AUCell, UCell, singscore, ssgsea, and AddModuleScore, were harnessed to scrutinize the activity of DDR-related genes across different cell types. Differential pathway alterations between high-and low-DDR activity cell clusters were compared. Furthermore, leveraging multiple RNA-seq datasets, we constructed a robust DDR-associated signature (DAS), and through integrative multiomics analysis, we explored differences in prognosis, pathways, mutational landscapes, and immunotherapy predictions among distinct DAS groups.

RESULTS

Notably, high-DDR activity cell subpopulations exhibited markedly enhanced cellular communication. The DAS demonstrated uniformity across multiple datasets. The low-DAS group exhibited improved prognoses, accompanied by heightened immune infiltration and elevated immune checkpoint expression. SubMap analysis of multiple immunotherapy datasets suggested that low-DAS group may experience enhanced immunotherapy responses. The "oncopredict" R package analyzed and screened sensitive drugs for different DAS groups.

CONCLUSION

Through the integration of single-cell and bulk RNA-seq data, we have developed a DAS associated with prognosis and immunotherapy response. This signature holds promise for the future stratification and personalized treatment of ESCC patients in clinical settings.

Droplet Hi-C for Fast and Scalable Profiling of Chromatin Architecture in Single Cells

Comprehensive analysis of chromatin architecture is crucial for understanding the gene regulatory programs during development and in disease pathogenesis, yet current methods often inadequately address the unique challenges presented by analysis of heterogeneous tissue samples. Here, we introduce Droplet Hi-C, which employs a commercial microfluidic device for high-throughput, single-cell chromatin conformation profiling in droplets. Using Droplet Hi-C, we mapped the chromatin architecture at single-cell resolution from the mouse cortex and analyzed gene regulatory programs in major cortical cell types. Additionally, we used this technique to detect copy number variation (CNV), structural variations (SVs) and extrachromosomal DNA (ecDNA) in cancer cells, revealing clonal dynamics and other oncogenic events during treatment. We further refined this technique to allow for joint profiling of chromatin architecture and transcriptome in single cells, facilitating a more comprehensive exploration of the links between chromatin architecture and gene expression in both normal tissues and tumors. Thus, Droplet Hi-C not only addresses critical gaps in chromatin analysis of heterogeneous tissues but also emerges as a versatile tool enhancing our understanding of gene regulation in health and disease.

Prognostic and immunotherapeutic potential of regulatory T cell-associated signature in ovarian cancer

Tumour-induced immunosuppressive microenvironments facilitate oncogenesis, with regulatory T cells (Tregs) serving as a crucial component. The significance of Treg-associated genes within the context of ovarian cancer (OC) remains elucidated insufficiently. Utilizing single-cell RNA sequencing (scRNA-Seq) for the identification of Treg-specific biomarkers, this investigation employed single-sample gene set enrichment analysis (ssGSEA) for the derivation of a Treg signature score. Weighted gene co-expression network analysis (WGCNA) facilitated the identification of Treg-correlated genes. Machine learning algorithms were employed to determine an optimal prognostic model, subsequently exploring disparities across risk strata in terms of survival outcomes, immunological infiltration, pathway activation and responsiveness to immunotherapy. Through WGCNA, a cohort of 365 Treg-associated genes was discerned, with 70 implicated in the prognostication of OC. A Tregs-associated signature (TAS), synthesized from random survival forest (RSF) and Least Absolute Shrinkage and Selection Operator (LASSO) algorithms, exhibited robust predictive validity across both internal and external cohorts. Low TAS OC patients demonstrated superior survival outcomes, augmented by increased immunological cell infiltration, upregulated immune checkpoint expression, distinct pathway enrichment and differential response to immunotherapeutic interventions. The devised TAS proficiently prognosticates patient outcomes and delineates the immunological milieu within OC, offering a strategic instrument for the clinical stratification and selection of patients.

Purine metabolism in lung adenocarcinoma: A single-cell analysis revealing prognostic and immunotherapeutic insights

Lung adenocarcinoma (LUAD) is a prevalent subtype of lung cancer, yet the contribution of purine metabolism (PM) to its pathogenesis remains poorly elucidated. PM, a critical component of intracellular nucleotide synthesis and energy metabolism, is hypothesized to exert a significant influence on LUAD development. Herein, we employed single-cell analysis to investigate the role of PM within the tumour microenvironment (TME) of LUAD. PM scoring (PMS) across distinct cell types was determined using AUCell, UCell, singscore and AddModuleScore algorithms. Subsequently, we explored communication networks among cells within high- and low-PMS groups, establishing a robust PM-associated signature (PAS) utilizing a comprehensive dataset comprising LUAD samples from TCGA and five GEO datasets. Our findings revealed that the high-PMS group exhibited intensified cell interactions, while the PAS, constructed using PM-related genes, demonstrated precise prognostic predictive capability. Notably, analysis across the TCGA dataset and five GEO datasets indicated that low-PAS patients exhibited a superior prognosis. Furthermore, the low-PAS group displayed increased immune cell infiltration and elevated CD8A expression, coupled with reduced PD-L1 expression. Moreover, data from eight publicly available immunotherapy cohorts suggested enhanced immunotherapy outcomes in the low-PAS group. These results underscore a close association between PAS and tumour immunity, offering predictive insights into genomic alterations, chemotherapy drug sensitivity and immunotherapy responses in LUAD. The newly established PAS holds promise as a valuable tool for selecting LUAD populations likely to benefit from future clinical stratification efforts.

HyGAnno: hybrid graph neural network-based cell type annotation for single-cell ATAC sequencing data

Reliable cell type annotations are crucial for investigating cellular heterogeneity in single-cell omics data. Although various computational approaches have been proposed for single-cell RNA sequencing (scRNA-seq) annotation, high-quality cell labels are still lacking in single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) data, because of extreme sparsity and inconsistent chromatin accessibility between datasets. Here, we present a novel automated cell annotation method that transfers cell type information from a well-labeled scRNA-seq reference to an unlabeled scATAC-seq target, via a parallel graph neural network, in a semi-supervised manner. Unlike existing methods that utilize only gene expression or gene activity features, HyGAnno leverages genome-wide accessibility peak features to facilitate the training process. In addition, HyGAnno reconstructs a reference-target cell graph to detect cells with low prediction reliability, according to their specific graph connectivity patterns. HyGAnno was assessed across various datasets, showcasing its strengths in precise cell annotation, generating interpretable cell embeddings, robustness to noisy reference data and adaptability to tumor tissues.

Ex vivo expansion potential of murine hematopoietic stem cells is a rare property only partially predicted by phenotype

The scarcity of hematopoietic stem cells (HSCs) restricts their use in both clinical settings and experimental research. Here, we examined a recently developed method for expanding rigorously purified murine HSCs ex vivo. After 3 weeks of culture, only 0.1% of cells exhibited the input HSC phenotype, but these accounted for almost all functional long-term HSC activity. Input HSCs displayed varying potential for ex vivo self-renewal, with alternative outcomes revealed by single-cell multimodal RNA and ATAC sequencing profiling. While most HSC progeny offered only transient in vivo reconstitution, these cells efficiently rescued mice from lethal myeloablation. The amplification of functional HSC activity allowed for long-term multilineage engraftment in unconditioned hosts that associated with a return of HSCs to quiescence. Thereby, our findings identify several key considerations for ex vivo HSC expansion, with major implications also for assessment of normal HSC activity.

Comprehensive scRNA-seq Model Reveals Artery Endothelial Cell Heterogeneity and Metabolic Preference in Human Vascular Disease

Vascular disease is one of the major causes of death worldwide. Endothelial cells are important components of the vascular structure. A better understanding of the endothelial cell changes in the development of vascular disease may provide new targets for clinical treatment strategies. Single-cell RNA sequencing can serve as a powerful tool to explore transcription patterns, as well as cell type identity. Our current study is based on comprehensive scRNA-seq data of several types of human vascular disease datasets with deep-learning-based algorithm. A gene set scoring system, created based on cell clustering, may help to identify the relative stage of the development of vascular disease. Metabolic preference patterns were estimated using a graphic neural network model. Overall, our study may provide potential treatment targets for retaining normal endothelial function under pathological situations.

Effective multi-modal clustering method via skip aggregation network for parallel scRNA-seq and scATAC-seq data

In recent years, there has been a growing trend in the realm of parallel clustering analysis for single-cell RNA-seq (scRNA) and single-cell Assay of Transposase Accessible Chromatin (scATAC) data. However, prevailing methods often treat these two data modalities as equals, neglecting the fact that the scRNA mode holds significantly richer information compared to the scATAC. This disregard hinders the model benefits from the insights derived from multiple modalities, compromising the overall clustering performance. To this end, we propose an effective multi-modal clustering model scEMC for parallel scRNA and Assay of Transposase Accessible Chromatin data. Concretely, we have devised a skip aggregation network to simultaneously learn global structural information among cells and integrate data from diverse modalities. To safeguard the quality of integrated cell representation against the influence stemming from sparse scATAC data, we connect the scRNA data with the aggregated representation via skip connection. Moreover, to effectively fit the real distribution of cells, we introduced a Zero Inflated Negative Binomial-based denoising autoencoder that accommodates corrupted data containing synthetic noise, concurrently integrating a joint optimization module that employs multiple losses. Extensive experiments serve to underscore the effectiveness of our model. This work contributes significantly to the ongoing exploration of cell subpopulations and tumor microenvironments, and the code of our work will be public at https://github.com/DayuHuu/scEMC.

Alzheimer DataLENS: An Open Data Analytics Portal for Alzheimer's Disease Research

Background

Recent Alzheimer's disease (AD) discoveries are increasingly based on studies from a variety of omics technologies on large cohorts. Currently, there is no easily accessible resource for neuroscientists to browse, query, and visualize these complex datasets in a harmonized manner.

Objective

Create an online portal of public omics datasets for AD research.

Methods

We developed Alzheimer DataLENS, a web-based portal, using the R Shiny platform to query and visualize publicly available transcriptomics and genetics studies of AD on human cohorts. To ensure consistent representation of AD findings, all datasets were processed through a uniform bioinformatics pipeline.

Results

Alzheimer DataLENS currently houses 2 single-nucleus RNA sequencing datasets, over 30 bulk RNA sequencing datasets from 19 brain regions and 3 cohorts, and 2 genome-wide association studies (GWAS). Available visualizations for single-nucleus data include bubble plots, heatmaps, and UMAP plots; for bulk expression data include box plots and heatmaps; for pathways include protein-protein interaction network plots; and for GWAS results include Manhattan plots. Alzheimer DataLENS also links to two other knowledge resources: the AD Progression Atlas and the Astrocyte Atlas.

Conclusions

Alzheimer DataLENS is a valuable resource for investigators to quickly and systematically explore omics datasets and is freely accessible at https://alzdatalens.partners.org.

Cross-species single-cell landscape of vertebrate pineal gland

The pineal gland has evolved from a photoreceptive organ in fish to a neuroendocrine organ in mammals. This study integrated multiple daytime single-cell RNA-seq datasets from the pineal glands of zebrafish, rats, and monkeys, providing a detailed examination of the evolutionary transition at single-cell resolution. We identified key factors responsible for the anatomical and functional transformation of the pineal gland. We retrieved and integrated daytime single-cell transcriptomic datasets from the pineal glands of zebrafish, rats, and monkeys, resulting in a total of 22 431 cells after rigorous quality filtering. Comparative analysis was then conducted to elucidate the evolution of pineal cells, their photosensitivity, their role in melatonin production, and the signaling processes within the glands of these species. Our analysis identified distinct cellular compositions of the pineal gland in zebrafish, rats, and monkeys. Zebrafish photoreceptors exhibited comprehensive phototransduction gene expression, while specific genes, including transducin (Gngt1, Gnb3, and Gngt2) and phosducin (Pdc), were consistently present in mammalian pinealocytes. We found transcriptional similarities between the pineal gland and retina, underscoring shared evolutionary and functional pathways. Zebrafish displayed unique light-responsive circadian gene activity compared to rats and monkeys. Key ligand-receptor interactions were identified, especially involving MDK and PTN, influencing melatonin synthesis across species. Furthermore, we observed species-specific GPCR (G protein-coupled receptors) expressions related to melatonin synthesis and their alignment with retinal expressions. Our findings also highlighted specific transcription factors (TFs) and regulatory networks associated with pineal gland evolution and function. Our study provides a detailed analysis of the pineal gland's evolution from fish to mammals. We identified key transcriptional changes and controls that highlight the gland's functional diversity. Notably, we found significant ligand-receptor interactions influencing melatonin synthesis and demonstrated parallels between pineal and retinal expressions. These insights enhance our understanding of the pineal gland's role in phototransduction, melatonin production, and circadian rhythms in vertebrates.

Characterization of the Prognosis and Tumor Microenvironment of Cellular Senescence-related Genes through scRNA-seq and Bulk RNA-seq Analysis in GC

BACKGROUND

Cellular senescence (CS) is thought to be the primary cause of cancer development and progression. This study aimed to investigate the prognostic role and molecular subtypes of CS-associated genes in gastric cancer (GC).

MATERIALS AND METHODS

The CellAge database was utilized to acquire CS-related genes. Expression data and clinical information of GC patients were obtained from The Cancer Genome Atlas (TCGA) database. Patients were then grouped into distinct subtypes using the "Consesus- ClusterPlus" R package based on CS-related genes. An in-depth analysis was conducted to assess the gene expression, molecular function, prognosis, gene mutation, immune infiltration, and drug resistance of each subtype. In addition, a CS-associated risk model was developed based on Cox regression analysis. The nomogram, constructed on the basis of the risk score and clinical factors, was formulated to improve the clinical application of GC patients. Finally, several candidate drugs were screened based on the Cancer Therapeutics Response Portal (CTRP) and PRISM Repurposing dataset.

RESULTS

According to the cluster result, patients were categorized into two molecular subtypes (C1 and C2). The two subtypes revealed distinct expression levels, overall survival (OS) and clinical presentations, mutation profiles, tumor microenvironment (TME), and drug resistance. A risk model was developed by selecting eight genes from the differential expression genes (DEGs) between two molecular subtypes. Patients with GC were categorized into two risk groups, with the high-risk group exhibiting a poor prognosis, a higher TME level, and increased expression of immune checkpoints. Function enrichment results suggested that genes were enriched in DNA repaired pathway in the low-risk group. Moreover, the Tumor Immune Dysfunction and Exclusion (TIDE) analysis indicated that immunotherapy is likely to be more beneficial for patients in the low-risk group. Drug analysis results revealed that several drugs, including ML210, ML162, dasatinib, idronoxil, and temsirolimus, may contribute to the treatment of GC patients in the high-risk group. Moreover, the risk model genes presented a distinct expression in single-cell levels in the GSE150290 dataset.

CONCLUSION

The two molecular subtypes, with their own individual OS rate, expression patterns, and immune infiltration, lay the foundation for further exploration into the GC molecular mechanism. The eight gene signatures could effectively predict the GC prognosis and can serve as reliable markers for GC patients.

Single-cell histone chaperones patterns guide intercellular communication of tumor microenvironment that contribute to breast cancer metastases

BACKGROUND

Histone chaperones (HCs) are crucial for governing genome stability and gene expression in multiple cancers. However, the functioning of HCs in the tumor microenvironment (TME) is still not clearly understood.

METHODS

Self-tested single-cell RNA-seq data derived from 6 breast cancer (BC) patients with brain and liver metastases were reanalyzed by nonnegative matrix factorization (NMF) algorithm for 36 HCs. TME subclusters were observed with BC and immunotherapy public cohorts to assess their prognosis and immune response. The biological effect of HSPA8, one of the HCs, was verified by transwell assay and wound-healing assays.

RESULTS

Cells including fibroblasts, macrophages, B cells, and T cells, were classified into various subclusters based on marker genes. Additionally, it showed that HCs might be strongly associated with biological and clinical features of BC metastases, along with the pseudotime trajectory of each TME cell type. Besides, the results of bulk-seq analysis revealed that TME cell subclusters mediated by HCs distinguished significant prognostic value for BC patients and were relevant to patients' immunotherapy responses, especially for B cells and macrophages. In particular, CellChat analysis exhibited that HCs-related TME cell subclusters revealed extensive and diverse interactions with malignant cells. Finally, transwell and wound-healing assays exhibited that HSPA8 deficiency inhibited BC cell migration and invasion.

CONCLUSIONS

Collectively, our study first dissected HCs-guided intercellular communication of TME that contribute to BC metastases.

Single-cell transcriptomics reveals heterogeneity in esophageal squamous epithelial cells and constructs models for predicting patient prognosis and immunotherapy

Background

Esophageal squamous cell carcinoma (ESCC), characterized by its high invasiveness and malignant potential, has long been a formidable challenge in terms of treatment.

Methods

A variety of advanced analytical techniques are employed, including single-cell RNA sequencing (scRNA-seq), cell trajectory inference, transcription factor regulatory network analysis, GSVA enrichment analysis, mutation profile construction, and the inference of potential immunotherapeutic drugs. The purpose is to conduct a more comprehensive exploration of the heterogeneity among malignant squamous epithelial cell subgroups within the ESCC microenvironment and establish a model for predicting the prognosis and immunotherapy outcomes of ESCC patients.

Results

An analysis was conducted through scRNA-seq, and three Cluster of malignant epithelial cells were identified using the infer CNV method. Cluster 0 was found to exhibit high invasiveness, whereas Cluster 1 displayed prominent characteristics associated with epithelial-mesenchymal transition. Confirmation of these findings was provided through cell trajectory analysis, which positioned Cluster 0 at the initiation stage of development and Cluster 1 at the final developmental stage. The abundance of Cluster 0-2 groups in TCGA-LUAD samples was assessed using ssGSEA and subsequently categorized into high and low-expression groups. Notably, it was observed that Cluster 0-1 had a significant impact on survival (p<0.05). Furthermore, GSVA enrichment analysis demonstrated heightened activity in hallmark pathways for Cluster 0, whereas Cluster 1 exhibited notable enrichment in pathways related to cell proliferation. It is noteworthy that a prognostic model was established utilizing feature genes from Cluster 0-1, employing the Lasso and stepwise regression methods. The results revealed that in TCGA and GSE53624 cohorts, the low-risk group demonstrated significantly higher overall survival and increased levels of immune infiltration. An examination of four external immunotherapy cohorts unveiled that the low-risk group exhibited improved immunotherapeutic efficacy. Additionally, more meaningful treatment options were identified for the low-risk group.

Conclusion

The findings revealed distinct interactions between malignant epithelial cells of ESCC and subgroups within the tumor microenvironment. Two cell clusters, strongly linked to survival, were pinpointed, and a signature was formulated. This signature is expected to play a crucial role in identifying and advancing precision medicine approaches for the treatment of ESCC.

Transfer learning for clustering single-cell RNA-seq data crossing-species and batch, case on uterine fibroids

Due to the high dimensionality and sparsity of the gene expression matrix in single-cell RNA-sequencing (scRNA-seq) data, coupled with significant noise generated by shallow sequencing, it poses a great challenge for cell clustering methods. While numerous computational methods have been proposed, the majority of existing approaches center on processing the target dataset itself. This approach disregards the wealth of knowledge present within other species and batches of scRNA-seq data. In light of this, our paper proposes a novel method named graph-based deep embedding clustering (GDEC) that leverages transfer learning across species and batches. GDEC integrates graph convolutional networks, effectively overcoming the challenges posed by sparse gene expression matrices. Additionally, the incorporation of DEC in GDEC enables the partitioning of cell clusters within a lower-dimensional space, thereby mitigating the adverse effects of noise on clustering outcomes. GDEC constructs a model based on existing scRNA-seq datasets and then applying transfer learning techniques to fine-tune the model using a limited amount of prior knowledge gleaned from the target dataset. This empowers GDEC to adeptly cluster scRNA-seq data cross different species and batches. Through cross-species and cross-batch clustering experiments, we conducted a comparative analysis between GDEC and conventional packages. Furthermore, we implemented GDEC on the scRNA-seq data of uterine fibroids. Compared results obtained from the Seurat package, GDEC unveiled a novel cell type (epithelial cells) and identified a notable number of new pathways among various cell types, thus underscoring the enhanced analytical capabilities of GDEC. Availability and implementation: https://github.com/YuzhiSun/GDEC/tree/main.

Integrating bulk and single-cell RNA sequencing data reveals epithelial-mesenchymal transition molecular subtype and signature to predict prognosis, immunotherapy efficacy, and drug candidates in low-grade gliomas

Objective: Epithelial-mesenchymal transition (EMT) is a tightly regulated and dynamic process occurring in both embryonic development and tumor progression. Our study aimed to comprehensively explore the molecular subtypes, immune landscape, and prognostic signature based on EMT-related genes in low-grade gliomas (LGG) in order to facilitate treatment decision-making and drug discovery. Methods: We curated EMT-related genes and performed molecular subtyping with consensus clustering algorithm to determine EMT expression patterns in LGG. The infiltration level of diverse immune cell subsets was evaluated by implementing the single-sample gene set enrichment analysis (ssGSEA) and ESTIMATE algorithms. The distinctions in clinical characteristics, mutation landscape, and immune tumor microenvironment (TME) among the subtypes were subjected to further investigation. Gene Set Variation Analysis (GSVA) was performed to explore the biological pathways that were involved in subtypes. The chemo drug sensitivity and immunotherapy of subtypes were estimated through GDSC database and NTP algorithm. To detect EMT subtype-related prognostic gene modules, the analysis of weighted gene co-expression network (WGCNA) was performed. The LASSO algorithm was utilized to construct a prognostic risk model, and its efficacy was verified through an independent CGGA dataset. Finally, the expression of the hub genes from the prognostic model was evaluated through the single-cell dataset and in-vitro experiment. Results: The TCGA-LGG dataset revealed the creation of two molecular subtypes that presented different prognoses, clinical implications, TME, mutation landscapes, chemotherapy, and immunotherapy. A three-gene signature (SLC39A1, CTSA and CLIC1) based on EMT expression pattern were established through WGCNA analysis. Low-risk patients showed a positive outlook, increased immune cell presence, and higher expression of immune checkpoint proteins. In addition, several promising drugs, including birinapant, fluvastatin, clofarabine, dasatinib, tanespimycin, TAK-733, GDC-0152, AZD8330, trametinib and ingenol-mebutate had great potential to the treatment of high risk patients. Finally, CTSA and CLIC1 were highly expressed in monocyte cell through single-cell RNA sequencing analysis. Conclusion: Our research revealed non-negligible role of EMT in the TME diversity and complexity of LGG. A prognostic signature may contribute to the personalized treatment and prognostic determination.

Multi-omics analyses unravel DNA damage repair-related clusters in breast cancer with experimental validation

Background

As one of the most common malignancies worldwide, breast cancer (BC) exhibits high heterogeneity of molecular phenotypes. The evolving view regarding DNA damage repair (DDR) is that it is context-specific and heterogeneous, but its role in BC remains unclear.

Methods

Multi-dimensional data of transcriptomics, genomics, and single-cell transcriptome profiling were obtained to characterize the DDR-related features of BC. We collected 276 DDR-related genes based on the Molecular Signature Database (MSigDB) database and previous studies. We acquired public datasets included the SCAN-B dataset (GEO: GSE96058), METABRIC database, and TCGA-BRCA database. Corresponding repositories such as transcriptomics, genomics, and clinical information were also downloaded. We selected scRNA-seq data from GEO: GSE176078, GSE114727, GSE161529, and GSE158724. Bulk RNA-seq data from GEO: GSE176078, GSE18728, GSE5462, GSE20181, and GSE130788 were extracted for independent analyses.

Results

The DDR classification was constructed in the SCAN-B dataset (GEO: GSE96058) and METABRIC database, Among BC patients, there were two clusters with distinct clinical and molecular characteristics: the DDR-suppressed cluster and the DDR-active cluster. A superior survival rate is found for tumors in the DDR-suppressed cluster, while those with the DDR-activated cluster tend to have inferior prognoses and clinically aggressive behavior. The DDR classification was validated in the TCGA-BRCA cohort and shown similar results. We also found that two clusters have different pathway activities at the genomic level. Based on the intersection of the different expressed genes among these cohorts, we found that PRAME might play a vital role in DDR. The DDR classification was then enabled by establishing a DDR score, which was verified through multilayer cohort analysis. Furthermore, our results revealed that malignant cells contributed more to the DDR score at the single-cell level than nonmalignant cells. Particularly, immune cells with immunosuppressive properties (such as FOXP3+ CD4+ T cells) displayed higher DDR scores among those with distinguishable characteristics.

Conclusion

Collectively, this study performs general analyses of DDR heterogeneity in BC and provides insight into the understanding of individualized molecular and clinicopathological mechanisms underlying unique DDR profiles.

Construction and validation of a novel prognostic model of neutrophil‑related genes signature of lung adenocarcinoma

Lung adenocarcinoma (LUAD) remains an incurable disease with a poor prognosis. This study aimed to explore neutrophil‑related genes (NRGs) and develop a prognostic signature for predicting the prognosis of LUAD. NRGs were obtained by intersecting modular genes identified by weighted gene co-expression network analysis (WGCNA) using bulk RNA-seq data and the marker genes of neutrophils identified from single-cell RNA-sequencing(scRNA-seq) data. Univariate Cox regression, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analyses were run to construct a prognostic signature, follow by delineation of risk groups, and external validation. Analyses of ESTIMAT, immune function, Tumor Immune Dysfunction and Exclusion (TIDE) scores, Immune cell Proportion Score (IPS), and immune checkpoint genes between high- and low-risk groups were performed, and then analyses of drug sensitivity to screen for sensitive anticancer drugs in high-risk groups. A total of 45 candidate NRGs were identified, of which PLTP, EREG, CD68, CD69, PLAUR, and CYP27A1 were considered to be significantly associated with prognosis in LUAD and were used to construct a prognostic signature. Correlation analysis showed significant differences in the immune landscape between high- and low-risk groups. In addition, our prognostic signature was important for predicting drug sensitivity in the high-risk group. Our study screened for NRGs in LUAD and constructed a novel and effective signature, revealing the immune landscape and providing more appropriate guidance protocols in LUAD treatment.

Single-cell analysis reveals exosome-associated biomarkers for prognostic prediction and immunotherapy in lung adenocarcinoma

BACKGROUND

Exosomes play a crucial role in tumor initiation and progression, yet the precise involvement of exosome-related genes (ERGs) in lung adenocarcinoma (LUAD) remains unclear.

METHODS

We conducted a comprehensive investigation of ERGs within the tumor microenvironment (TME) of LUAD using single-cell RNA sequencing (scRNA-seq) analysis. Multiple scoring methods were employed to assess exosome activity (EA). Differences in cell communication were examined between high and low EA groups, utilizing the "CellChat" R package. Subsequently, we leveraged multiple bulk RNA-seq datasets to develop and validate exosome-associated signatures (EAS), enabling a multifaceted exploration of prognosis and immunotherapy outcomes between high- and low-risk groups.

RESULTS

In the LUAD TME, epithelial cells demonstrated the highest EA, with even more elevated levels observed in advanced LUAD epithelial cells. The high-EA group exhibited enhanced intercellular interactions. EAS were established through the analysis of multiple bulk RNA-seq datasets. Patients in the high-risk group exhibited poorer overall survival (OS), reduced immune infiltration, and decreased expression of immune checkpoint genes. Finally, we experimentally validated the high expression of SEC61G in LUAD cell lines and demonstrated that knockdown of SEC61G reduced the proliferative capacity of LUAD cells using colony formation assays.

CONCLUSION

The integration of single-cell and bulk RNA-seq analyses culminated in the development of the profound and significant EAS, which imparts invaluable insights for the clinical diagnosis and therapeutic management of LUAD patients.

The integrated single-cell analysis developed an immunogenic cell death signature to predict lung adenocarcinoma prognosis and immunotherapy

BACKGROUND

Research on immunogenic cell death (ICD) in lung adenocarcinoma (LUAD) has been relatively limited. This study aims to create ICD-related signatures for accurate survival prognosis prediction in LUAD patients, addressing the challenge of lacking reliable early prognostic indicators for this type of cancer.

METHODS

Using single-cell RNA sequencing (scRNA-seq) analysis, ICD activity in cells was calculated by AUCell algorithm, divided into high- and low-ICD groups according to median values, and key ICD regulatory genes were identified through differential analysis, and these genes were integrated into TCGA data to construct prognostic signatures using LASSO and COX regression analysis, and multi-dimensional analysis of ICD-related signatures in terms of prognosis, immunotherapy, tumor microenvironment (TME), and mutational landscape.

RESULTS

The constructed signature reveals a pronounced disparity in prognosis between the high- and low-risk groups of LUAD patients. The statistical discrepancies in survival times among LUAD patients from both the TCGA and GEO databases further corroborate this observation. Additionally, heightened levels of immune cell infiltration expression are evidenced in the low-risk group, suggesting a potential benefit from immunotherapeutic interventions for these patients. The expression levels of pivotal risk-associated genes in tissue samples were assessed utilizing qRT-PCR, thereby unveiling PITX3 as a plausible therapeutic target in the context of LUAD.

CONCLUSIONS

Our constructed ICD-related signatures provide help in predicting the prognosis and immunotherapy of LUAD patients, and to some extent guide the clinical treatment of LUAD patients.

Deciphering Treg cell roles in esophageal squamous cell carcinoma: a comprehensive prognostic and immunotherapeutic analysis

Background: Esophageal squamous cell carcinoma (ESCC) is a prevalent and aggressive form of cancer that poses significant challenges in terms of prognosis and treatment. Regulatory T cells (Treg cells) have gained attention due to their influential role in immune modulation within the tumor microenvironment (TME). Understanding the intricate interactions between Treg cells and the tumor microenvironment is essential for unraveling the mechanisms underlying ESCC progression and for developing effective prognostic models and immunotherapeutic strategies. Methods: A combination of single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq analysis was utilized to explore the role of Treg cells within the TME of ESCC. The accuracy and applicability of the prognostic model were assessed through multi-dimensional evaluations, encompassing an examination of the model's performance across various dimensions, such as the mutation landscape, clinical relevance, enrichment analysis, and potential implications for immunotherapy strategies. Results: The pivotal role of the macrophage migration inhibitory factor (MIF) signaling pathway within the ESCC TME was investigated, with a focus on its impact on Treg cells and other subpopulations. Through comprehensive integration of bulk sequencing data, a Treg-associated signature (TAS) was constructed, revealing that ESCC patients with elevated TAS (referred to as high-TAS individuals) experienced significantly improved prognoses. Heightened immune infiltration and increased expression of immune checkpoint markers were observed in high-TAS specimens. The model's validity was established through the IMvigor210 dataset, demonstrating its robustness in predicting prognosis and responsiveness to immunotherapy. Heightened therapeutic benefits were observed in immune-based interventions for high-TAS ESCC patients. Noteworthy differences in pathway enrichment patterns emerged between high and low-TAS cohorts, highlighting potential avenues for therapeutic exploration. Furthermore, the clinical relevance of key model genes was substantiated by analyzing clinical samples from ten paired tumor and adjacent tissues, revealing differential expression levels. Conclusion: The study established a TAS that enables accurate prediction of patient prognosis and responsiveness to immunotherapy. This achievement holds significant implications for the clinical management of ESCC, offering valuable insights for informed therapeutic interventions.

Multimodal deep learning approaches for single-cell multi-omics data integration

Integrating single-cell multi-omics data is a challenging task that has led to new insights into complex cellular systems. Various computational methods have been proposed to effectively integrate these rapidly accumulating datasets, including deep learning. However, despite the proven success of deep learning in integrating multi-omics data and its better performance over classical computational methods, there has been no systematic study of its application to single-cell multi-omics data integration. To fill this gap, we conducted a literature review to explore the use of multimodal deep learning techniques in single-cell multi-omics data integration, taking into account recent studies from multiple perspectives. Specifically, we first summarized different modalities found in single-cell multi-omics data. We then reviewed current deep learning techniques for processing multimodal data and categorized deep learning-based integration methods for single-cell multi-omics data according to data modality, deep learning architecture, fusion strategy, key tasks and downstream analysis. Finally, we provided insights into using these deep learning models to integrate multi-omics data and better understand single-cell biological mechanisms.

Revealing the role of regulatory T cells in the tumor microenvironment of lung adenocarcinoma: a novel prognostic and immunotherapeutic signature

Background

Regulatory T cells (Tregs), are a key class of cell types in the immune system. In the tumor microenvironment (TME), the presence of Tregs has important implications for immune response and tumor development. Relatively little is known about the role of Tregs in lung adenocarcinoma (LUAD).

Methods

Tregs were identified using but single-cell RNA sequencing (scRNA-seq) analysis and interactions between Tregs and other cells in the TME were investigated. Next, we used multiple bulk RNA-seq datasets to construct risk models based on marker genes of Tregs and explored differences in prognosis, mutational landscape, immune cell infiltration and immunotherapy between high- and low-risk groups, and finally, qRT-PCR and cell function experiments were performed to validate the model genes.

Results

The cellchat analysis showed that MIF-(CD74+CXCR4) pairs play a key role in the interaction of Tregs with other cell subpopulations, and the Tregs-associated signatures (TRAS) could well classify multiple LUAD cohorts into high- and low-risk groups. Immunotherapy may offer greater potential benefits to the low-risk group, as indicated by their superior survival, increased infiltration of immune cells, and heightened expression of immune checkpoints. Finally, the experiment verified that the model genes LTB and PTTG1 were relatively highly expressed in cancer tissues, while PTPRC was relatively highly expressed in paracancerous tissues. Colony Formation assay confirmed that knockdown of PTTG1 reduced the proliferation ability of LUAD cells.

Conclusion

TRAS were constructed using scRNA-seq and bulk RNA-seq to distinguish patient risk subgroups, which may provide assistance in the clinical management of LUAD patients.

A novel signature predicts prognosis and immunotherapy in lung adenocarcinoma based on cancer-associated fibroblasts

BACKGROUND

Extensive research has established the significant correlations between cancer-associated fibroblasts (CAFs) and various stages of cancer development, including initiation, angiogenesis, progression, and resistance to therapy. In this study, we aimed to investigate the characteristics of CAFs in lung adenocarcinoma (LUAD) and develop a risk signature to predict the prognosis of patients with LUAD.

METHODS

We obtained single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the public database. The Seurat R package was used to process the scRNA-seq data and identify CAF clusters based on several biomarkers. CAF-related prognostic genes were further identified using univariate Cox regression analysis. To reduce the number of genes, Lasso regression was performed, and a risk signature was established. A novel nomogram that incorporated the risk signature and clinicopathological features was developed to predict the clinical applicability of the model. Additionally, we conducted immune landscape and immunotherapy responsiveness analyses. Finally, we performed in vitro experiments to verify the functions of EXO1 in LUAD.

RESULTS

We identified 5 CAF clusters in LUAD using scRNA-seq data, of which 3 clusters were significantly associated with prognosis in LUAD. A total of 492 genes were found to be significantly linked to CAF clusters from 1731 DEGs and were used to construct a risk signature. Moreover, our immune landscape exploration revealed that the risk signature was significantly related to immune scores, and its ability to predict responsiveness to immunotherapy was confirmed. Furthermore, a novel nomogram incorporating the risk signature and clinicopathological features showed excellent clinical applicability. Finally, we verified the functions of EXP1 in LUAD through in vitro experiments.

CONCLUSIONS

The risk signature has proven to be an excellent predictor of LUAD prognosis, stratifying patients more appropriately and precisely predicting immunotherapy responsiveness. The comprehensive characterization of LUAD based on the CAF signature can predict the response of LUAD to immunotherapy, thus offering fresh perspectives into the management of LUAD patients. Our study ultimately confirms the role of EXP1 in facilitating the invasion and growth of tumor cells in LUAD. Nevertheless, further validation can be achieved by conducting in vivo experiments.

Integrating multiple machine learning methods to construct glutamine metabolism-related signatures in lung adenocarcinoma

Background

Glutamine metabolism (GM) is known to play a critical role in cancer development, including in lung adenocarcinoma (LUAD), although the exact contribution of GM to LUAD remains incompletely understood. In this study, we aimed to discover new targets for the treatment of LUAD patients by using machine learning algorithms to establish prognostic models based on GM-related genes (GMRGs).

Methods

We used the AUCell and WGCNA algorithms, along with single-cell and bulk RNA-seq data, to identify the most prominent GMRGs associated with LUAD. Multiple machine learning algorithms were employed to develop risk models with optimal predictive performance. We validated our models using multiple external datasets and investigated disparities in the tumor microenvironment (TME), mutation landscape, enriched pathways, and response to immunotherapy across various risk groups. Additionally, we conducted in vitro and in vivo experiments to confirm the role of LGALS3 in LUAD.

Results

We identified 173 GMRGs strongly associated with GM activity and selected the Random Survival Forest (RSF) and Supervised Principal Components (SuperPC) methods to develop a prognostic model. Our model's performance was validated using multiple external datasets. Our analysis revealed that the low-risk group had higher immune cell infiltration and increased expression of immune checkpoints, indicating that this group may be more receptive to immunotherapy. Moreover, our experimental results confirmed that LGALS3 promoted the proliferation, invasion, and migration of LUAD cells.

Conclusion

Our study established a prognostic model based on GMRGs that can predict the effectiveness of immunotherapy and provide novel approaches for the treatment of LUAD. Our findings also suggest that LGALS3 may be a potential therapeutic target for LUAD.

Mast cell marker gene signature: prognosis and immunotherapy response prediction in lung adenocarcinoma through integrated scRNA-seq and bulk RNA-seq

Background

Mast cells, comprising a crucial component of the tumor immune milieu, modulate neoplastic progression by secreting an array of pro- and antitumorigenic factors. Numerous extant studies have produced conflicting conclusions regarding the impact of mast cells on the prognosis of patients afflicted with lung adenocarcinoma (LUAD).

Methods

Employing single-cell RNA sequencing (scRNA-seq) analysis, mast cell-specific marker genes in LUAD were ascertained. Subsequently, a mast cell-related genes (MRGs) signature was devised to stratify LUAD patients into high- and low-risk cohorts based on the median risk value. Further investigations were conducted to assess the influence of distinct risk categories on the tumor microenvironment. The prognostic import and capacity to prognosticate immunotherapy benefits of the MRGs signature were corroborated using four external cohorts. Ultimately, the functional roles of SYAP1 were validated through in vitro experimentation.

Results

After scRNA-seq and bulk RNA-seq data analysis, we established a prognostic signature consisting of nine MRGs. This profile effectively distinguished favorable survival outcomes in both the training and validation cohorts. In addition, we identified the low-risk group as a population more effective for immunotherapy. In cellular experiments, we found that silencing SYAP1 significantly reduced the proliferation, invasion and migratory capacity of LUAD cells while increasing apoptosis.

Conclusion

Our MRGs signature offers valuable insights into the involvement of mast cells in determining the prognosis of LUAD and may prove instrumental as a navigational aid for immunotherapy selection, as well as a predictor of immunotherapy response in LUAD patients.