Dissection of triple-negative breast cancer microenvironment and identification of potential therapeutic drugs using single-cell RNA sequencing analysis.

Breast cancer remains a leading cause of mortality in women worldwide. Triple-negative breast cancer (TNBC) is a particularly aggressive subtype characterized by rapid progression, poor prognosis, and lack of clear therapeutic targets. In the clinic, delineation of tumor heterogeneity and development of effective drugs continue to pose considerable challenges. Within the scope of our study, high heterogeneity inherent to breast cancer was uncovered based on the landscape constructed from both tumor and healthy breast tissue samples. Notably, TNBC exhibited significant specificity regarding cell proliferation, differentiation, and disease progression. Significant associations between tumor grade, prognosis, and TNBC oncogenes were established via pseudotime trajectory analysis. Consequently, we further performed comprehensive characterization of the TNBC microenvironment. A crucial epithelial subcluster, E8, was identified as highly malignant and strongly associated with tumor cell proliferation in TNBC. Additionally, epithelial-mesenchymal transition (EMT)-associated fibroblast and M2 macrophage subclusters exerted an influence on E8 through cellular interactions, contributing to tumor growth. Characteristic genes in these three cluster cells could therefore serve as potential therapeutic targets for TNBC. The collective findings provided valuable insights that assisted in the screening of a series of therapeutic drugs, such as pelitinib. We further confirmed the anti-cancer effect of pelitinib in an orthotopic 4T1 tumor-bearing mouse model. Overall, our study sheds light on the unique characteristics of TNBC at single-cell resolution and the crucial cell types associated with tumor cell proliferation that may serve as potent tools in the development of effective anti-cancer drugs.

Optimizing drug combination and mechanism analysis based on risk pathway crosstalk in pan cancer.

Combination therapy can greatly improve the efficacy of cancer treatment, so identifying the most effective drug combination and interaction can accelerate the development of combination therapy. Here we developed a computational network biological approach to identify the effective drug which inhibition risk pathway crosstalk of cancer, and then filtrated and optimized the drug combination for cancer treatment. We integrated high-throughput data concerning pan-cancer and drugs to construct miRNA-mediated crosstalk networks among cancer pathways and further construct networks for therapeutic drug. Screening by drug combination method, we obtained 687 optimized drug combinations of 83 first-line anticancer drugs in pan-cancer. Next, we analyzed drug combination mechanism, and confirmed that the targets of cancer-specific crosstalk network in drug combination were closely related to cancer prognosis by survival analysis. Finally, we save all the results to a webpage for query ( http://bio-bigdata.hrbmu.edu.cn/oDrugCP/ ). In conclusion, our study provided an effective method for screening precise drug combinations for various cancer treatments, which may have important scientific significance and clinical application value for tumor treatment.

Multiscale integrative analyses unveil immune-related diagnostic signature for the progression of MASLD.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is a chronic liver disease prevalent worldwide, with an increasing incidence associated with obesity, diabetes, and metabolic syndrome. The progression of MASLD to metabolic dysfunction-associated steatohepatitis (MASH) poses a pressing health concern, highlighting the significance of accurately identifying MASLD and its progression to MASH as a primary challenge in the field. In this study, a systematic integration of 66 immune cell types was conducted. Comprehensive analyses were performed on bulk, single-cell RNA-Seq, and clinical data to investigate the immune cell types implicated in MASLD progression thoroughly. Multiple approaches, including immune infiltration, gene expression trend analysis, weighted gene coexpression network analysis, and 4 machine learning algorithms, were used to examine the dynamic changes in genes and immune cells during MASLD progression. C-X-C motif chemokine receptor 4 and dedicator of cytokinesis 8 have been identified as potential diagnostic biomarkers for MASLD progression. Furthermore, cell communication analysis at the single-cell level revealed that the involvement of C-X-C motif chemokine receptor 4 and dedicator of cytokinesis 8 in MASLD progression is mediated through their influence on T cells. Overall, our study identified vital immune cells and a 2-gene diagnostic signature for the progression of MASLD, providing a new perspective on the diagnosis and immune-related molecular mechanisms of MASLD. These findings have important implications for developing innovative diagnostic tools and therapies for MASLD.

Identifying and characterizing drug sensitivity-related lncRNA-TF-gene regulatory triplets.

Recently, many studies have shown that lncRNA can mediate the regulation of TF-gene in drug sensitivity. However, there is still a lack of systematic identification of lncRNA-TF-gene regulatory triplets for drug sensitivity. In this study, we propose a novel analytic approach to systematically identify the lncRNA-TF-gene regulatory triplets related to the drug sensitivity by integrating transcriptome data and drug sensitivity data. Totally, 1570 drug sensitivity-related lncRNA-TF-gene triplets were identified, and 16 307 relationships were formed between drugs and triplets. Then, a comprehensive characterization was performed. Drug sensitivity-related triplets affect a variety of biological functions including drug response-related pathways. Phenotypic similarity analysis showed that the drugs with many shared triplets had high similarity in their two-dimensional structures and indications. In addition, Network analysis revealed the diverse regulation mechanism of lncRNAs in different drugs. Also, survival analysis indicated that lncRNA-TF-gene triplets related to the drug sensitivity could be candidate prognostic biomarkers for clinical applications. Next, using the random walk algorithm, the results of which we screen therapeutic drugs for patients across three cancer types showed high accuracy in the drug-cell line heterogeneity network based on the identified triplets. Besides, we developed a user-friendly web interface-DrugSETs (http://bio-bigdata.hrbmu.edu.cn/DrugSETs/) available to explore 1570 lncRNA-TF-gene triplets relevant with 282 drugs. It can also submit a patient's expression profile to predict therapeutic drugs conveniently. In summary, our research may promote the study of lncRNAs in the drug resistance mechanism and improve the effectiveness of treatment.

Prognostic Features of the Tumor Immune Microenvironment in Glioma and Their Clinical Applications: Analysis of Multiple Cohorts.

Glioma is the most common malignant tumor of the central nervous system. Tumor purity is a source of important prognostic factor for glioma patients, showing the key roles of the microenvironment in glioma prognosis. In this study, we systematically screened functional characterization related to the tumor immune microenvironment and constructed a risk model named Glioma MicroEnvironment Functional Signature (GMEFS) based on eight cohorts. The prognostic value of the GMEFS model was also verified in another two glioma cohorts, glioblastoma (GBM) and low-grade glioma (LGG) cohorts, from The Cancer Genome Atlas (TCGA). Nomograms were established in the training and testing cohorts to validate the clinical use of this model. Furthermore, the relationships between the risk score, intrinsic molecular subtypes, tumor purity, and tumor-infiltrating immune cell abundance were also evaluated. Meanwhile, the performance of the GMEFS model in glioma formation and glioma recurrence was systematically analyzed based on 16 glioma cohorts from the Gene Expression Omnibus (GEO) database. Based on multiple-cohort integrated analysis, risk subpathway signatures were identified, and a drug-subpathway association network was further constructed to explore candidate therapy target regions. Three subpathways derived from Focal adhesion (path: 04510) were identified and contained known targets including platelet derived growth factor receptor alpha (PDGFRA), epidermal growth factor receptor (EGFR), and erb-b2 receptor tyrosine kinase 2 (ERBB2). In conclusion, the novel functional signatures identified in this study could serve as a robust prognostic biomarker, and this study provided a framework to identify candidate therapeutic target regions, which further guide glioma patients' clinical decision.

Comprehensive analysis of microglia gene and subpathway signatures for glioma prognosis and drug screening: linking microglia to glioma.

Glioma is the most common malignant tumors in the brain. Previous studies have revealed that, as the innate immune cells in nervous system, microglia cells were involved in glioma pathology. And, the resident microglia displayed its specific biological roles which distinguished with peripheral macrophages. In this study, an integrated analysis was performed based on public resource database to explore specific biological of microglia within glioma. Through comprehensive analysis, the biological characterization underlying two conditions, glioma microglia compared to glioma macrophage (MicT/MacT) as well as glioma microglia compared to normal microglia (MicT/MicN), were revealed. Notably, nine core MicT/MicN genes displayed closely associations with glioma recurrence and prognosis, such as P2RY2, which was analyzed in more than 2800 glioma samples from 25 studies. Furthermore, we applied a random walk based strategy to identify microglia specific subpathways and developed SubP28 signature for glioma prognostic analysis. Multiple validation data sets confirmed the predictive performance of SubP28 and involvement in molecular subtypes. The associations between SuP28 score and microglia M1/M2 polarization were also explored for both GBM and LGG types. Finally, a comprehensive drug-subpathway network was established for screening candidate medicable molecules (drugs) and identifying therapeutic subpathway targets. In conclusions, the comprehensive analysis of microglia related gene and functional signatures in glioma pathobiologic events by large-scale data sets displayed a framework to dissect inner connection between microglia and glioma, and identify robust signature for glioma clinical implications.

Identifying individualized risk subpathways reveals pan-cancer molecular classification based on multi-omics data.

Cancer is a highly heterogeneous disease with different functional disorders among individuals. The initiation and progression of cancer is usually related to dysregulation of local regions within pathways. Identification of individualized risk pathways is crucial for revealing the mechanisms of tumorigenesis and heterogeneity. However, approach that focused on mining patient-specific risk subpathway regions is still lacking. Here, we developed an individualized cancer risk subpathway identification method that was referred as InCRiS by integrating multi-omics data. Then, the method was applied to nearly 3000 samples across 9 TCGA cancer types and its robustness and reliability were evaluated. Dissecting dysregulated subpathways in these tumor samples revealed several key regions which may play oncogenic roles in cancer. The construction of risk subpathway dysregulation profile of pan-cancers revealed 11 pan-cancer molecular classification (InCRiS subtypes) with significantly different clinical outcomes. Moreover, subpathway regions that tend to be disordered in individuals of specific subtypes were examined for understanding the pathogenesis of tumor and some key genes such as CTNNB1, EP300 and PRKDC were nominated in different subtypes. In summary, the proposed method and resulting data presented useful resources to study the mechanism of tumor heterogeneity and to discovery novel therapeutic targets for precise treatment of cancer.

Dissecting immune cell stat regulation network reveals biomarkers to predict ICB therapy responders in melanoma.

BACKGROUND: Immunotherapy is a revolutionary strategy in cancer therapy, but the resistance of which is one of the important challenges. Detecting the regulation of immune cells and biomarkers concerning immune checkpoint blockade (ICB) therapy is of great significance. METHODS: Here, we firstly constructed regulation networks for 11 immune cell clusters by integrating biological pathway data and single cell sequencing data in metastatic melanoma with or without ICB therapy. We then dissected these regulation networks and identified differently expressed genes between responders and non-responders. Finally, we trained and validated a logistic regression model based on ligands and receptors in the regulation network to predict ICB therapy response. RESULTS: We discovered the regulation of genes across eleven immune cell stats. Functional analysis indicated that these stat-specific networks consensually enriched in immune response corrected pathways and highlighted antigen processing and presentation as a core pathway in immune cell regulation. Furthermore, some famous ligands like SIRPA, ITGAM, CD247and receptors like CD14, IL2 and HLA-G were differently expressed between cells of responders and non-responders. A predictive model of gene sets containing ligands and receptors performed accuracy prediction with AUCs above 0.7 in a validation dataset suggesting that they may be server as biomarkers for predicting immunotherapy response. CONCLUSIONS: In summary, our study presented the gene-gene regulation landscape across 11 immune cell clusters and analysis of these networks revealed several important aspects and immunotherapy response biomarkers, which may provide novel insights into immune related mechanisms and immunotherapy response prediction.

Identification of a Novel Prognostic Signature for Gastric Cancer Based on Multiple Level Integration and Global Network Optimization.

Gastric Cancer (GC) is a common cancer worldwide with a high morbidity and mortality rate in Asia. Many prognostic signatures from genes and non-coding RNA (ncRNA) levels have been identified by high-throughput expression profiling for GC. To date, there have been no reports on integrated optimization analysis based on the GC global lncRNA-miRNA-mRNA network and the prognostic mechanism has not been studied. In the present work, a Gastric Cancer specific lncRNA-miRNA-mRNA regulatory network (GCsLMM) was constructed based on the ceRNA hypothesis by combining miRNA-target interactions and data on the expression of GC. To mine for novel prognostic signatures associated with GC, we performed topological analysis, a random walk with restart algorithm, in the GCsLMM from three levels, miRNA-, mRNA-, and lncRNA-levels. We further obtained candidate prognostic signatures by calculating the integrated score and analyzed the robustness of these signatures by combination strategy. The biological roles of key candidate signatures were also explored. Finally, we targeted the PHF10 gene and analyzed the expression patterns of PHF10 in independent datasets. The findings of this study will improve our understanding of the competing endogenous RNA (ceRNA) regulatory mechanisms and further facilitate the discovery of novel prognostic biomarkers for GC clinical guidelines.

The identification of six risk genes for ovarian cancer platinum response based on global network algorithm and verification analysis.

Ovarian cancer is the most lethal gynaecological cancer, and resistance of platinum-based chemotherapy is the main reason for treatment failure. The aim of the present study was to identify candidate genes involved in ovarian cancer platinum response by analysing genes from homologous recombination and Fanconi anaemia pathways. Associations between these two functional genes were explored in the study, and we performed a random walk algorithm based on reconstructed gene-gene network, including protein-protein interaction and co-expression relations. Following the random walk, all genes were ranked and GSEA analysis showed that the biological functions focused primarily on autophagy, histone modification and gluconeogenesis. Based on three types of seed nodes, the top two genes were utilized as examples. We selected a total of six candidate genes (FANCA, FANCG, POLD1, KDM1A, BLM and BRCA1) for subsequent verification. The validation results of the six candidate genes have significance in three independent ovarian cancer data sets with platinum-resistant and platinum-sensitive information. To explore the correlation between biomarkers and clinical prognostic factors, we performed differential analysis and multivariate clinical subgroup analysis for six candidate genes at both mRNA and protein levels. And each of the six candidate genes and their neighbouring genes with a mutation rate greater than 10% were also analysed by network construction and functional enrichment analysis. In the meanwhile, the survival analysis for platinum-treated patients was performed in the current study. Finally, the RT-qPCR assay was used to determine the performance of candidate genes in ovarian cancer platinum response. Taken together, this research demonstrated that comprehensive bioinformatics methods could help to understand the molecular mechanism of platinum response and provide new strategies for overcoming platinum resistance in ovarian cancer treatment.

Comprehensive analysis of mRNA-level and miRNA-level subpathway activities for identifying robust ovarian cancer prognostic signatures.

Ovarian cancer (OvCa) causes the highest mortality among all gynaecologic cancers. A large number of mRNA- or miRNA-based signatures were identified for OvCa patient prognosis. However, the comprehensive analysis of function-level prognostic signatures is currently not considered in OvCa. In the present study, we respectively inferred subpathway activities from mRNA and miRNA levels based on high-throughput expression profiles and reconstructed subpathways. Firstly, the activities of two tumour pathways were calculated and the difference between normal and tumour samples were analysed using multiple tumour types. Then, we calculated subpathway activities for OvCa based on the expression profiles from both mRNA and miRNA levels. Furthermore, based on these subpathway activity matrices, we performed bootstrap analysis to obtain sub-training sets and utilized univariate method to identify robust OvCa prognostic subpathways. A comprehensive comparison of subpathway results between these two levels was performed. As a result, we observed subpathway mutual exclusion trend between the levels of mRNA and miRNA, which indicated the necessary of combining mRNA-miRNA levels. Finally, by using ICGC data as testing sets, we utilized two strategies to verify survival predictive power of the mRNA-miRNA combined subpathway signatures and performed comparisons with results from individual levels. It was confirmed that our framework displayed application to identify robust and efficient prognostic signatures for OvCa, and the combined signatures indeed exhibited advantages over individual ones. In the study, we took a step forward in relevant novel integrated functional signatures for OvCa prognosis.

PRSS1 Upregulation Predicts Platinum Resistance in Ovarian Cancer Patients.

Ovarian cancer is the most frequent cause of death among gynecologic malignancies. A total of 80% of patients who have completed platinum-based chemotherapy suffer from relapse and develop resistance within 2 years. In the present study, we obtained patients' complete platinum (cisplatin and carboplatin) medication information from The Cancer Genome Atlas database and then divided them into two categories: resistance and sensitivity. Difference analysis was performed to screen differentially expressed genes (DEgenes) related to platinum response. Subsequently, we annotated DEgenes into the protein-protein interaction network as seed nodes and analyzed them by random walk. Finally, second-ranking protease serine 1 gene (PRSS1) was selected as a candidate gene for verification analysis. PRSS1's expression pattern was continuously studied in Oncomine and cBio Cancer Genomic Portal databases, revealing the key roles of PRSS1 in ovarian cancer formation. Hereafter, we conducted in-depth explorations on PRSS1's platinum response to ovarian cancer through tissue and cytological experiments. Quantitative real-time polymerase chain reaction and Western blot assay results indicated that PRSS1 expression levels in platinum-resistant samples (tissue/cell) were significantly higher than in samples sensitive to platinum. By cell transfection assay, we observed that knockdown of PRSS1 reduced the resistance of ovarian cancer cells to cisplatin. Meanwhile, overexpression of PRSS1 increased the resistance to cisplatin. In conclusion, we identified a novel risk gene PRSS1 related to ovarian cancer platinum response and confirmed its key roles using multiple levels of low-throughput experiments, revealing a new treatment strategy based on a novel target factor for overcoming cisplatin resistance in ovarian cancer.