scPanel: a tool for automatic identification of sparse gene panels for generalizable patient classification using scRNA-seq datasets.

Single-cell RNA sequencing (scRNA-seq) technologies can generate transcriptomic profiles at a single-cell resolution in large patient cohorts, facilitating discovery of gene and cellular biomarkers for disease. Yet, when the number of biomarker genes is large, the translation to clinical applications is challenging due to prohibitive sequencing costs. Here, we introduce scPanel, a computational framework designed to bridge the gap between biomarker discovery and clinical application by identifying a sparse gene panel for patient classification from the cell population(s) most responsive to perturbations (e.g. diseases/drugs). scPanel incorporates a data-driven way to automatically determine a minimal number of informative biomarker genes. Patient-level classification is achieved by aggregating the prediction probabilities of cells associated with a patient using the area under the curve score. Application of scPanel to scleroderma, colorectal cancer, and COVID-19 datasets resulted in high patient classification accuracy using only a small number of genes (<20), automatically selected from the entire transcriptome. In the COVID-19 case study, we demonstrated cross-dataset generalizability in predicting disease state in an external patient cohort. scPanel outperforms other state-of-the-art gene selection methods for patient classification and can be used to identify parsimonious sets of reliable biomarker candidates for clinical translation.

SAE-Impute: imputation for single-cell data via subspace regression and auto-encoders.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) technology has emerged as a crucial tool for studying cellular heterogeneity. However, dropouts are inherent to the sequencing process, known as dropout events, posing challenges in downstream analysis and interpretation. Imputing dropout data becomes a critical concern in scRNA-seq data analysis. Present imputation methods predominantly rely on statistical or machine learning approaches, often overlooking inter-sample correlations. RESULTS: To address this limitation, We introduced SAE-Impute, a new computational method for imputing single-cell data by combining subspace regression and auto-encoders for enhancing the accuracy and reliability of the imputation process. Specifically, SAE-Impute assesses sample correlations via subspace regression, predicts potential dropout values, and then leverages these predictions within an autoencoder framework for interpolation. To validate the performance of SAE-Impute, we systematically conducted experiments on both simulated and real scRNA-seq datasets. These results highlight that SAE-Impute effectively reduces false negative signals in single-cell data and enhances the retrieval of dropout values, gene-gene and cell-cell correlations. Finally, We also conducted several downstream analyses on the imputed single-cell RNA sequencing (scRNA-seq) data, including the identification of differential gene expression, cell clustering and visualization, and cell trajectory construction. CONCLUSIONS: These results once again demonstrate that SAE-Impute is able to effectively reduce the droupouts in single-cell dataset, thereby improving the functional interpretability of the data.

Identification of functional heterogeneity of immune cells and tubular-immune cellular interplay action in diabetic kidney disease.

Background: Renal inflammation plays key roles in the pathogenesis of diabetic kidney disease (DKD). Immune cell infiltration is the main pathological feature in the progression of DKD. Sodium glucose cotransporter 2 inhibitor (SGLT2i) were reported to have antiinflammatory effects on DKD. While the heterogeneity and molecular basis of the pathogenesis and treatment with SGLT2i in DKD remains poorly understood. Methods: To address this question, we performed a single-cell transcriptomics data analysis and cell cross-talk analysis based on the database (GSE181382). The single-cell transcriptome analysis findings were validated using multiplex immunostaining. Results: A total of 58760 cells are categorized into 25 distinct cell types. A subset of macrophages with anti-inflammatory potential was identified. We found that Ccl3+ (S100a8/a9 high) macrophages with anti-inflammatory and antimicrobial in the pathogenesis of DKD decreased and reversed the dapagliflozin treatment. Besides, dapagliflozin treatment enhanced the accumulation of Pck1+ macrophage, characterized by gluconeogenesis signaling pathway. Cell-cross talk analysis showed the GRN/SORT1 pair and CD74 related signaling pathways were enriched in the interactions between tubular epithelial cells and immune cells. Conclusions: Our study depicts the heterogeneity of macrophages and clarifies a new possible explanation of dapagliflozin treatment, showing the metabolism shifts toward gluconeogenesis in macrophages, fueling the anti-inflammatory function of M2 macrophages, highlighting the new molecular features and signaling pathways and potential therapeutic targets, which has provided an important reference for the study of immune-related mechanisms in the progression of the disease.

Dissecting the Single-Cell Diversity and Heterogeneity Underlying Cervical Precancerous Lesions and Cancer Tissues.

The underlying cellular diversity and heterogeneity from cervix precancerous lesions to cervical squamous cell carcinoma (CSCC) is investigated. Four single-cell datasets including normal tissues, normal adjacent tissues, precancerous lesions, and cervical tumors were integrated to perform disease stage analysis. Single-cell compositional data analysis (scCODA) was utilized to reveal the compositional changes of each cell type. Differentially expressed genes (DEGs) among cell types were annotated using BioCarta. An assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis was performed to correlate epigenetic alterations with gene expression profiles. Lastly, a logistic regression model was used to assess the similarity between the original and new cohort data (HRA001742). After global annotation, seven distinct cell types were categorized. Eight consensus-upregulated DEGs were identified in B cells among different disease statuses, which could be utilized to predict the overall survival of CSCC patients. Inferred copy number variation (CNV) analysis of epithelial cells guided disease progression classification. Trajectory and ATAC-seq integration analysis identified 95 key transcription factors (TF) and one immunohistochemistry (IHC) testified key-node TF (YY1) involved in epithelial cells from CSCC initiation to progression. The consistency of epithelial cell subpopulation markers was revealed with single-cell sequencing, bulk sequencing, and RT-qPCR detection. KRT8 and KRT15, markers of Epi6, showed progressively higher expression with disease progression as revealed by IHC detection. The logistic regression model testified the robustness of the resemblance of clusters among the various datasets utilized in this study. Valuable insights into CSCC cellular diversity and heterogeneity provide a foundation for future targeted therapy.

Integrated Single-cell RNA-seq and Bulk RNA-seq Identify Diagnostic Biomarkers for Postmenopausal Osteoporosis.

AIM: We aimed to explore diagnostic biomarkers of postmenopausal osteoporosis (PMOP). BACKGROUND: PMOP brings enormous physical and economic burden to elderly women. OBJECTIVES: This study aims to screen new biomarkers for osteoporosis, providing insights for early diagnosis and therapeutic targets of osteoporosis. METHODS: Weighted gene co-expression network analysis (WGCNA) was applied to identify osteoporosis-related hub genes. Single-cell transcriptomic atlas of osteoporosis was depicted and the heterogeneity of monocytes was analyzed, based on which the biomarkers for osteoporosis were screened. Gene set enrichment analysis (GSEA) was conducted on the biomarkers. The diagnostic model (nomogram) was established and evaluated based on the expression levels of biomarkers. Additionally, the transcription factor (TF) regulatory network was constructed to predict the potential TF and targeted miRNA of biomarkers. The drugs with significant correlation with biomarkers were identified by Spearman correlation analysis. RESULTS: We obtained 30 osteoporosis-associated hub genes. 9 cell types were identified, and the monocytes were subdivided to 4 subtypes. Three biomarkers, DHX29, LSM5, and UBE2V2, were screened. DHX29 and UBE2V2 were highly expressed in non-classical monocytes, while LSM5 exhibited the highest expression in other monocytes, followed by non-classical monocytes. GSEA indicated that osteoporosis may be correlated with vascular calcification and the biomarkers may be involved in the formation of immune cells. Then, nomogram was constructed and exhibited good robustness. In addition, MYC and SETDB1 were the shared IF in three biomarkers, which may play critical regulatory roles in the progression of osteoporosis. Moreover, 41, 49, and 68 drugs appeared significant correlations with DHX29, LSM5, and UBE2V2, respectively. CONCLUSION: This study provided a basis for early diagnosis and targeted treatment of osteoporosis.

Identification and validation of SHC1 and FGFR1 as novel immune-related oxidative stress biomarkers of non-obstructive azoospermia.

Background: Non-obstructive azoospermia (NOA) is a major contributor of male infertility. Herein, we used existing datasets to identify novel biomarkers for the diagnosis and prognosis of NOA, which could have great significance in the field of male infertility. Methods: NOA datasets were obtained from the Gene Expression Omnibus (GEO) database. CIBERSORT was utilized to analyze the distributions of 22 immune cell populations. Hub genes were identified by applying weighted gene co-expression network analysis (WGCNA), machine learning methods, and protein-protein interaction (PPI) network analysis. The expression of hub genes was verified in external datasets and was assessed by receiver operating characteristic (ROC) curve analysis. Gene set enrichment analysis (GSEA) was applied to explore the important functions and pathways of hub genes. The mRNA-microRNA (miRNA)-transcription factors (TFs) regulatory network and potential drugs were predicted based on hub genes. Single-cell RNA sequencing data from the testes of patients with NOA were applied for analyzing the distribution of hub genes in single-cell clusters. Furthermore, testis tissue samples were obtained from patients with NOA and obstructive azoospermia (OA) who underwent testicular biopsy. RT-PCR and Western blot were used to validate hub gene expression. Results: Two immune-related oxidative stress hub genes (SHC1 and FGFR1) were identified. Both hub genes were highly expressed in NOA samples compared to control samples. ROC curve analysis showed a remarkable prediction ability (AUCs > 0.8). GSEA revealed that hub genes were predominantly enriched in toll-like receptor and Wnt signaling pathways. A total of 24 TFs, 82 miRNAs, and 111 potential drugs were predicted based on two hub genes. Single-cell RNA sequencing data in NOA patients indicated that SHC1 and FGFR1 were highly expressed in endothelial cells and Leydig cells, respectively. RT-PCR and Western blot results showed that mRNA and protein levels of both hub genes were significantly upregulated in NOA testis tissue samples, which agree with the findings from analysis of the microarray data. Conclusion: It appears that SHC1 and FGFR1 could be significant immune-related oxidative stress biomarkers for detecting and managing patients with NOA. Our findings provide a novel viewpoint for illustrating potential pathogenesis in men suffering from infertility.

Dynamics of alternative polyadenylation in single root cells of Arabidopsis thaliana.

Introduction: Single-cell RNA-seq (scRNA-seq) technologies have been widely used to reveal the diversity and complexity of cells, and pioneering studies on scRNA-seq in plants began to emerge since 2019. However, existing studies on plants utilized scRNA-seq focused only on the gene expression regulation. As an essential post-transcriptional mechanism for regulating gene expression, alternative polyadenylation (APA) generates diverse mRNA isoforms with distinct 3' ends through the selective use of different polyadenylation sites in a gene. APA plays important roles in regulating multiple developmental processes in plants, such as flowering time and stress response. Methods: In this study, we developed a pipeline to identify and integrate APA sites from different scRNA-seq data and analyze APA dynamics in single cells. First, high-confidence poly(A) sites in single root cells were identified and quantified. Second, three kinds of APA markers were identified for exploring APA dynamics in single cells, including differentially expressed poly(A) sites based on APA site expression, APA markers based on APA usages, and APA switching genes based on 3' UTR (untranslated region) length change. Moreover, cell type annotations of single root cells were refined by integrating both the APA information and the gene expression profile. Results: We comprehensively compiled a single-cell APA atlas from five scRNA-seq studies, covering over 150,000 cells spanning four major tissue branches, twelve cell types, and three developmental stages. Moreover, we quantified the dynamic APA usages in single cells and identified APA markers across tissues and cell types. Further, we integrated complementary information of gene expression and APA profiles to annotate cell types and reveal subtle differences between cell types. Discussion: This study reveals that APA provides an additional layer of information for determining cell identity and provides a landscape of APA dynamics during Arabidopsis root development.

Detecting expressed genes in cell populations at the single-cell level with scGeneXpress.

Determining whether genes are expressed or not remains a challenge in single-cell RNAseq experiments due to their different expression spectra, which are influenced by genetics, the microenvironment and gene length. Current approaches for addressing this issue fail to provide a comprehensive landscape of expressed genes, since they neglect the inherent differences in the expression ranges and distributions of genes. Here, we present scGeneXpress, a method for detecting expressed genes in cell populations of single-cell RNAseq samples based on gene-specific reference distributions. We demonstrate that scGeneXpress accurately detects expressed cell markers and identity genes in 34 human and mouse tissues and can be employed to improve differential expression analysis of single-cell RNAseq data.

Single-cell transcriptome analysis of the mouse lungs during the injury and recovery periods after lipopolysaccharide administration.

OBJECTIVE: This study sought to investigate the cellular and molecular alterations during the injury and recovery periods of ALI and develop effective treatments for ALI. METHODS: Pulmonary histology at 1, 3, 6, and 9 days after lipopolysaccharide administration mice were assessed. An unbiased single-cell RNA sequencing was performed in alveoli tissues from injury (day 3) and recovery (day 6) mice after lipopolysaccharide administration. The roles of Fpr2 and Dpp4 in ALI were assessed. RESULTS: The most severe lung injury occurred on day 3, followed by recovery entirely on day 9 after lipopolysaccharide administration. The numbers of Il1a+ neutrophils, monocytes/macrophages, and Cd4+ and Cd8+ T cells significantly increased at day 3 after LPS administration; subsequently, the number of Il1a+ neutrophils greatly decreased, the numbers of monocytes/macrophages and Cd4+ and Cd8+ T cells continuously increased, and the number of resident alveolar macrophages significantly increased at day 6. The interactions between monocytes/macrophages and pneumocytes during the injury period were enhanced by the Cxcl10/Dpp4 pair, and inhibiting Dpp4 improved ALI significantly, while inhibiting Fpr2 did not. CONCLUSIONS: Our results offer valuable insights into the cellular and molecular mechanisms underlying its progression and identify Dpp4 as an effective therapeutic target for ALI.

Germinal center B-cell subgroups in the tumor microenvironment cannot be overlooked: Their involvement in prognosis, immunotherapy response, and treatment resistance in head and neck squamous carcinoma.

Background: More than 60 % of patients with head and neck squamous carcinoma (HNSCC) are diagnosed at advanced stages and miss radical treatment. This has prompted the need to find new biomarkers to achieve early diagnosis and predict early recurrence and metastasis of tumors. Methods: Single-cell RNA sequencing (scRNA-seq) data from HNSCC tissues and peripheral blood samples were obtained through the Gene Expression Omnibus (GEO) database (GSE164690) to characterize the B-cell subgroups, differentiation trajectories, and intercellular communication networks in HNSCC and to construct a prognostic model of the associated risks. In addition, this study analyzed the differences in clinical features, immune cell infiltration, functional enrichment, tumor mutational burden (TMB), and drug sensitivity between the high- and low-risk groups. Results: Using scRNA-seq of HNSCC, we classified B and plasma cells into a total of four subgroups: naive B cells (NBs), germinal center B cells (GCBs), memory B cells (MBs), and plasma cells (PCs). Pseudotemporal trajectory analysis revealed that NBs and GCBs were at the early stage of B cell differentiation, while MBs and PCs were at the end. Cellular communication revealed that GCBs acted on tumor cells through the CD99 and SEMA4 signaling pathways. The independent prognostic value, immune cell infiltration, TMB and drug sensitivity assays were validated for the MEF2B+ GCB score groups. Conclusions: We identified GCBs as B cell-specific prognostic biomarkers for the first time. The MEF2B+ GCB score fills the research gap in the genetic prognostic prediction model of HNSCC and is expected to provide a theoretical basis for finding new therapeutic targets for HNSCC.

Leveraging explainable deep learning methodologies to elucidate the biological underpinnings of Huntington's disease using single-cell RNA sequencing data.

BACKGROUND: Huntington's disease (HD) is a hereditary neurological disorder caused by mutations in HTT, leading to neuronal degeneration. Traditionally, HD is associated with the misfolding and aggregation of mutant huntingtin due to an extended polyglutamine domain encoded by an expanded CAG tract. However, recent research has also highlighted the role of global transcriptional dysregulation in HD pathology. However, understanding the intricate relationship between mRNA expression and HD at the cellular level remains challenging. Our study aimed to elucidate the underlying mechanisms of HD pathology using single-cell sequencing data. RESULTS: We used single-cell RNA sequencing analysis to determine differential gene expression patterns between healthy and HD cells. HD cells were effectively modeled using a residual neural network (ResNet), which outperformed traditional and convolutional neural networks. Despite the efficacy of our approach, the F1 score for the test set was 96.53%. Using the SHapley Additive exPlanations (SHAP) algorithm, we identified genes influencing HD prediction and revealed their roles in HD pathobiology, such as in the regulation of cellular iron metabolism and mitochondrial function. SHAP analysis also revealed low-abundance genes that were overlooked by traditional differential expression analysis, emphasizing its effectiveness in identifying biologically relevant genes for distinguishing between healthy and HD cells. Overall, the integration of single-cell RNA sequencing data and deep learning models provides valuable insights into HD pathology. CONCLUSION: We developed the model capable of analyzing HD at single-cell transcriptomic level.

Impact of Mitophagy-Related Genes on the Diagnosis and Development of Esophageal Squamous Cell Carcinoma via Single-Cell RNA-seq Analysis and Machine Learning Algorithms.

As a treatment for esophageal squamous cell carcinoma (ESCC), which is common and fatal, mitophagy is a conserved cellular mechanism that selectively removes damaged mitochondria and is crucial for cellular homeostasis. While tumor development and resistance to anticancer therapies are related to ESCC, their role in ESCC remains unclear. Here, we investigated the relationship between mitophagy-related genes (MRGs) and ESCC to provide novel insights into the role of mitophagy in ESCC prognosis and diagnosis prediction. First, we identified MRGs from the GeneCards database and examined them at both the single-cell and transcriptome levels. Key genes were selected and a prognostic model was constructed using least absolute shrinkage and selection operator analysis. External validation was performed using the GSE53624 dataset and Kaplan-Meier survival analysis was performed to identify PYCARD as a gene significantly associated with survival in ESCC. We then examined the effect of PYCARD on ESCC cell proliferation and migration and identified 169 MRGs at the single-cell and transcriptome levels, as well as the high-risk groups associated with cancer-related pathways. Thirteen key genes were selected for model construction via multiple machine learning algorithms. PYCARD, which is upregulated in patients with ESCC, was negatively correlated with prognosis and its knockdown inhibited ESCC cell proliferation and migration. Our ESCC prediction model based on mitophagy-related genes demonstrated promising results and provides more options for the management and clinical treatment of ESCC patients. Moreover, targeting or regulating PYCARD levels might offer new therapeutic strategies for ESCC patients in clinical settings.

MSIsensor-RNA: Microsatellite Instability Detection for Bulk and Single-cell Gene Expression Data.

Microsatellite instability (MSI) is an indispensable biomarker in cancer immunotherapy. Currently, MSI scoring methods by high-throughput omics methods have gained popularity and demonstrated better performance than the gold standard method for MSI detection. However, the MSI detection method on expression data, especially single-cell expression data, is still lacking, limiting the scope of clinical application and prohibiting the investigation of MSI at a single-cell level. Herein, we developed MSIsensor-RNA, an accurate, robust, adaptable, and standalone software to detect MSI status based on expression values of MSI-associated genes. We demonstrated the favorable performance and promise of MSIsensor-RNA in both bulk and single-cell gene expression data in multiplatform technologies including RNA sequencing (RNA-seq), microarray, and single-cell RNA-seq. MSIsensor-RNA is a versatile, efficient, and robust method for MSI status detection from both bulk and single-cell gene expression data in clinical studies and applications. MSIsensor-RNA is available at https://github.com/xjtu-omics/msisensor-rna.

Pan-cancer integrative analyses dissect the remodeling of endothelial cells in human cancers.

Therapeutics targeting tumor endothelial cells (TECs) have been explored for decades, with only suboptimal efficacy achieved, partly due to an insufficient understanding of the TEC heterogeneity across cancer patients. We integrated single-cell RNA-seq data of 575 cancer patients from 19 solid tumor types, comprehensively charting the TEC phenotypic diversities. Our analyses uncovered underappreciated compositional and functional heterogeneity in TECs from a pan-cancer perspective. Two subsets, CXCR4 + tip cells and SELE + veins, represented the prominent angiogenic and proinflammatory phenotypes of TECs, respectively. They exhibited distinct spatial organization patterns, and compared to adjacent non-tumor tissues, tumor tissue showed an increased prevalence of CXCR4 + tip cells, yet with SELE + veins depleted. Such functional and spatial characteristics underlie their differential associations with the response of anti-angiogenic therapies and immunotherapies. Our integrative resources and findings open new avenues to understand and clinically intervene in the tumor vasculature.

Characteristics of PD-1+CD4+ T cells in peripheral blood and synovium of rheumatoid arthritis patients.

Objectives: PD-1 plays a crucial role in the immune dysregulation of rheumatoid arthritis (RA), but the specific characteristics of PD-1+CD4+ T cells remain unclear and require further investigation. Methods: Circulating PD-1+CD4+ T cells from RA patients were analysed using flow cytometry. Plasma levels of soluble PD-1 (sPD-1) were measured using enzyme-linked immunosorbent assay (ELISA). Single-cell RNA sequence data from peripheral blood mononuclear cells (PBMCs) and synovial tissue of patients were obtained from the GEO and the ImmPort databases. Bioinformatics analyses were performed in the R studio to characterise PD-1+CD4+ T cells. Expression of CCR7, KLF2 and IL32 in PD-1+CD4+ T cells was validated by flow cytometry. Results: RA patients showed an elevated proportion of PD-1+CD4+ T cells in peripheral blood, along with increased plasma sPD-1 levels, which positively correlated with TNF-α and erythrocyte sedimentation rate. Bioinformatic analysis revealed PD-1 expression on CCR7+CD4+ T cells in PBMCs, and on both CCR7+CD4+ T cells and CXCL13+CD4+ T cells in RA synovium. PD-1 was co-expressed with CCR7, KLF2, and IL32 in peripheral CD4+ T cells. In synovium, PD-1+CCR7+CD4+ T cells had higher expression of TNF and LCP2, while PD-1+CXCL13+CD4+ T cells showed elevated levels of ARID5A and DUSP2. PD-1+CD4+ T cells in synovium also appeared to interact with B cells and fibroblasts through BTLA and TNFSF signalling pathways. Conclusion: This study highlights the increased proportion of PD-1+CD4+ T cells and elevated sPD-1 levels in RA. The transcriptomic profiles and signalling networks of PD-1+CD4+ T cells offer new insights into their role in RA pathogenesis.

HPV-driven heterogeneity in cervical cancer: study on the role of epithelial cells and myofibroblasts in the tumor progression based on single-cell RNA sequencing analysis.

Background: Cervical cancer (CC) is a neoplasia with a high heterogeneity. We aimed to explore the characteristics of tumor microenvironment (TME) for CC treatment. Methods: HPV positive (+) and negative (-) samples from cervical cancer (CC) patients were sourced from the Gene Expression Omnibus (GEO) database. The single-cell RNA sequencing (scRNA-seq) data were processed and annotated for cell types utilizing the Seurat package. Following this, the expression levels and biological roles of the marker genes were analyzed applying real-time PCR (RT-PCR) and transwell assays. Furthermore, the enrichment of genes with significantly differential expressions and copy number variations was assessed by the ClusterProlifer and inferCNV software packages. Results: Seven main cell clusters were classified based on a total of 12,431 cells. The HPV- CC samples exhibited a higher immune cell infiltration level, while epithelial cells and myofibroblasts had higher proportion in the HPV+ CC samples with extensive heterogeneity. Immune pathways including antigen treatment and presentation, immunoglobulin production and T cell mediated immunity were significantly activated in the HPV- CC group with lower cell cycle and proliferation activity. However, the anti-tumor immunity of these cells was inhibited in HPV+ CC group with higher cell proliferation activity. Moreover, the amplification and loss of CNVs also supported that these cells in HPV- CC samples were prone to anti-tumor activation. Further cell validation results showed that except GZMA, the levels of APOC1, CEACAM6, FOXP3, SFRP4 and TFF3 were all higher in CC cells Hela, and that silencing TFF3 could inhibit the migration and invasion of CC cells in-vitro. Conclusion: This study highlighted the critical role of HPV infection in CC progression, providing a novel molecular basis for optimizing the current preventive screening and personalized treatment for the cancer.

Causal effects of systemic inflammatory proteins on Guillain-Barre Syndrome: insights from genome-wide Mendelian randomization, single-cell RNA sequencing analysis, and network pharmacology.

Background: Evidence from observational studies indicates that inflammatory proteins play a vital role in Guillain-Barre Syndrome (GBS). Nevertheless, it is unclear how circulating inflammatory proteins are causally associated with GBS. Herein, we conducted a two-sample Mendelian randomization (MR) analysis to systematically explore the causal links of genetically determined systemic inflammatory proteins on GBS. Methods: A total of 8,293 participants of European ancestry were included in a genome-wide association study of 41 inflammatory proteins as instrumental variables. Five MR approaches, encompassing inverse-variance weighted, weighted median, MR-Egger, simple model, and weighted model were employed to explore the causal links between inflammatory proteins and GBS. MR-Egger regression was utilized to explore the pleiotropy. Cochran's Q statistic was implemented to quantify the heterogeneity. Furthermore, we performed single-cell RNA sequencing analysis and predicted potential drug targets through molecular docking technology. Results: By applying MR analysis, four inflammatory proteins causally associated with GBS were identified, encompassing IFN-γ (OR:1.96, 95%CI: 1.02-3.78, PIVW=0.045), IL-7 (OR:1.86, 95%CI: 1.07-3.23, PIVW=0.029), SCGF-ß (OR:1.56, 95%CI: 1.11-2.19, PIVW=0.011), and Eotaxin (OR:1.99, 95%CI: 1.01-3.90, PIVW=0.046). The sensitivity analysis revealed no evidence of pleiotropy or heterogeneity. Additionally, significant genes were found through single-cell RNA sequencing analysis and several anti-inflammatory or neuroprotective small molecular compounds were identified by utilizing molecular docking technology. Conclusions: Our MR analysis suggested that IFN-γ, IL-7, SCGF-ß, and Eotaxin were causally linked to the occurrence and development of GBS. These findings elucidated potential causal associations and highlighted the significance of these inflammatory proteins in the pathogenesis and prospective therapeutic targets for GBS.

Upregulation of CKS2 in immunosuppressive cells is associated with metastasis and poor prognosis in prostate cancer: a single-cell RNA-sequencing analysis.

Background: Metastasis worsens prostate cancer (PCa) prognosis, with the immunosuppressive microenvironment playing a key role in bone metastasis. This study aimed to investigate how an immunosuppressive environment promotes PCa metastasis and worsens prognosis of patients with PCa. Methods: Candidate oncogenes were identified through analysis of the Gene Expression Omnibus (GEO) database. A prognostic model was developed for the purpose of identifying target genes. A single-cell RNA sequencing data from GEO database was used to analyze the localization of target genes in the tumor microenvironment. A pan-cancer analysis was conducted to study the cancer-causing potential of target genes across different types of tumors. Results: Fifty-one genes were found to be differentially expressed in bone metastasis compared to non-metastatic PCa, with CKS2 identified as the most significant gene associated with poor prognosis. CKS2 was shown to be linked to an immunosuppressive microenvironment and osteoclastic bone metastases, as shown by its negative correlation with immune cell infiltration and osteoblast-related gene expression. Moreover, CKS2 was found in immunosuppressive cells and was linked to bone metastasis in PCa. It was also overexpressed in different types of tumors, making it as an oncogenic gene. Conclusions: This research offers a new perspective on the potential utility of CKS2 as a therapeutic target for the prevention of metastatic PCa.

scHyper: reconstructing cell-cell communication through hypergraph neural networks.

Cell-cell communications is crucial for the regulation of cellular life and the establishment of cellular relationships. Most approaches of inferring intercellular communications from single-cell RNA sequencing (scRNA-seq) data lack a comprehensive global network view of multilayered communications. In this context, we propose scHyper, a new method that can infer intercellular communications from a global network perspective and identify the potential impact of all cells, ligand, and receptor expression on the communication score. scHyper designed a new way to represent tripartite relationships, by extracting a heterogeneous hypergraph that includes the source (ligand expression), the target (receptor expression), and the relevant ligand-receptor (L-R) pairs. scHyper is based on hypergraph representation learning, which measures the degree of match between the intrinsic attributes (static embeddings) of nodes and their observed behaviors (dynamic embeddings) in the context (hyperedges), quantifies the probability of forming hyperedges, and thus reconstructs the cell-cell communication score. Additionally, to effectively mine the key mechanisms of signal transmission, we collect a rich dataset of multisubunit complex L-R pairs and propose a nonparametric test to determine significant intercellular communications. Comparing with other tools indicates that scHyper exhibits superior performance and functionality. Experimental results on the human tumor microenvironment and immune cells demonstrate that scHyper offers reliable and unique capabilities for analyzing intercellular communication networks. Therefore, we introduced an effective strategy that can build high-order interaction patterns, surpassing the limitations of most methods that can only handle low-order interactions, thus more accurately interpreting the complexity of intercellular communications.

Deciphering single-cell gene expression variability and its role in drug response.

The effectiveness of drug treatments is profoundly influenced by individual responses, which are shaped by gene expression variability, particularly within pharmacogenes. Leveraging single-cell RNA sequencing (scRNA-seq) data, our study explores the extent of expression variability among pharmacogenes in a wide array of cell types across eight different human tissues, shedding light on their impact on drug responses. Our findings broaden the established link between variability in pharmacogene expression and drug efficacy to encompass variability at the cellular level. Moreover, we unveil a promising approach to enhance drug efficacy prediction. This is achieved by leveraging a combination of cross-cell and cross-individual pharmacogene expression variation measurements. Our study opens avenues for more precise forecasting of drug performance, facilitating tailored and more effective treatments in the future.