A Metabolism-Related Gene Signature Predicts the Prognosis of Breast Cancer Patients: Combined Analysis of High-Throughput Sequencing and Gene Chip Data Sets

Construction of mitochondrial quality regulation genes-related prognostic model based on bulk-RNA-seq analysis in multiple myeloma

Mitochondrial quality regulation plays an important role in affecting the treatment sensitivity of multiple myeloma (MM). We aimed to develop a mitochondrial quality regulation genes (MQRGs)-related prognostic model for MM patients. The Genomic Data Commons-MM of bulk RNA-seq, mutation, and single-cell RNA-seq (scRNA-seq) dataset were downloaded, and the MQRGs gene set was collected previous study. "maftools" and CIBERSORT were used for mutation and immune-infiltration analysis. Subsequently, the "ConsensusClusterPlus" was used to perform the unsupervised clustering analysis, "survminer" and "ssGSEA" R package was used for the Kaplan-Meier survival and enrichment analysis, "limma" R, univariate and Least Absolute Shrinkage and Selection Operator Cox were used for RiskScore model. The "timeROC" R package was used for Receiver Operating Characteristic Curve analysis. Finally, the "Seurat" R package was used for scRNA-seq analysis. These MQRGs are mainly located on chromosome-1,2,3,7, and 22 and had significant expression differences among age, gender, and stage groups, in which PPARGC1A and PPARG are the high mutation genes. Most MQRGs expression are closely associated with the plasma cells infiltration and can divide the patients into 2 different prognostic clusters (C1, C2). Then, 8 risk models were screened from 60 DEGs for RiskScore, which is an independent prognostic factor and effectively divided the patients into high and low risk groups with significant difference of immune checkpoint expression. Nomogram containing RiskScore can accurately predict patient prognosis, and a series of specific transcription factor PRDM1 and IRF1 were identified. We described the based molecular features and developed a high effective MQRGs-related prognostic model in MM.

Identification of novel diagnostic biomarkers associated with liver metastasis in colon adenocarcinoma by machine learning

BACKGROUND

Liver metastasis is one of the primary causes of poor prognosis in colon adenocarcinoma (COAD) patients, but there are few studies on its biomarkers.

METHODS

The Cancer Genome Atlas (TCGA)-COAD, GSE41258, and GSE49355 datasets were acquired from the public database. Differentially expressed genes (DEGs) between liver metastasis and primary tumor samples in COAD were identified by limma, and functional enrichment analysis were performed. MuTect2 and maftools were used to measure somatic mutation rates, while ADTEx was used to measure copy number variations (CNVs). The intersection of three machine learning methods, support vector machine (SVM), Random Forest, and least absolute shrinkage and selection operator (LASSO), is utilized to screen biomarkers, and their diagnostic performance is subsequently validated. The correlation between biomarkers and immune cells infiltration was analyzed by Spearman method.

RESULTS

47 DEGs between liver metastasis and primary tumor samples in COAD were obtained, which were mainly enriched in the complement and coagulation, extracellular matrix (ECM), and peptidase regulator activity, etc. 38 out of 47 DEGs had mutations and exhibited a high frequency of CNV amplification or deletion. Furthermore, 3 biomarkers (MMP3, MAB21L2, and COLEC11) were screened, which showed good diagnostic performance. The proportion of multiple immune cells, such as B cells naive, T cells CD4 naive, Monocytes, and Dendritic cells resting, was higher in liver metastasis samples than that in primary tumor samples. Meanwhile, MMP3, MAB21L2, and COLEC11 exhibited an outstanding correlation with immune cells infiltration.

CONCLUSION

In short, 3 biomarkers with good diagnostic efficacy were identified, providing a new perspective of therapeutic targets for liver metastasis in COAD.

High infiltration of immune cells with lower immune activity mediated the heterogeneity of gastric adenocarcinoma and promoted metastasis

Background

Gastric adenocarcinoma (GA) is a heterogeneous malignancy with high invasion and metastasis. We aimed to explore the metastatic characteristics of GA using single-cell RNA-sequencing (scRNA-seq) analysis.

Methods

The scRNA-seq dataset was downloaded from the GEO database and the "Seurat" package was used to perform the scRNA-seq analysis. The CellMarker2.0 database provided gene markers. Subsequently, differentially expressed genes (DEGs) were identified using the FindMarkers function and subjected to enrichment analysis with the "ClusterProlifer". "GseaVis" package was used for visualizing the gene levels. Finally, the SCENIC analysis was performed for identifying key regulons. The expression level and functionality of the key genes were verified by quantitative real-time PCR (qRT-PCR), wound healing and transwell assays.

Results

A total of 7697 cells were divided into 8 cell subsets, in which the Cytotoxic NK/T cells, Myeloid cells and Myofibroblasts had higher proportion in the metastatic tissues. Further screening of DEGs and enrichment analysis revealed that in the metastatic tissues, NK cells, monocytes and inflammatory fibroblasts with low immune levels contributed to GA metastasis. In addition, this study identified a series of key immune-related regulons that mediated the lower immune activity of immune cells. Further in vitro experiment verified that CXCL8 was a key factor mediating the proliferation and migration of GA cells.

Conclusion

The scRNA-seq analysis showed that high infiltration of immune cells with lower immune activity mediated heterogeneity to contribute to GA metastasis.

Polycystic ovarian syndrome (PCOS) and recurrent spontaneous abortion (RSA) are associated with the PI3K-AKT pathway activation

Aims

We aimed to elucidate the mechanism leading to polycystic ovarian syndrome (PCOS) and recurrent spontaneous abortion (RSA).

Background

PCOS is an endocrine disorder. Patients with RSA also have a high incidence rate of PCOS, implying that PCOS and RSA may share the same pathological mechanism.

Objective

The single-cell RNA-seq datasets of PCOS (GSE168404 and GSE193123) and RSA GSE113790 and GSE178535) were downloaded from the Gene Expression Omnibus (GEO) database.

Methods

Datasets of PSCO and RSA patients were retrieved from the Gene Expression Omnibus (GEO) database. The "WGCNA" package was used to determine the module eigengenes associated with the PCOS and RSA phenotypes and the gene functions were analyzed using the "DAVID" database. The GSEA analysis was performed in "clusterProfiler" package, and key genes in the activated pathways were identified using the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Real-time quantitative PCR (RT-qPCR) was conducted to determine the mRNA level. Cell viability and apoptosis were measured by cell counting kit-8 (CCK-8) and flow cytometry, respectively.

Results

The modules related to PCOS and RSA were sectioned by weighted gene co-expression network analysis (WGCNA) and positive correlation modules of PCOS and RSA were all enriched in angiogenesis and Wnt pathways. The GSEA further revealed that these biological processes of angiogenesis, Wnt and regulation of cell cycle were significantly positively correlated with the PCOS and RSA phenotypes. The intersection of the positive correlation modules of PCOS and RSA contained 80 key genes, which were mainly enriched in kinase-related signal pathways and were significant high-expressed in the disease samples. Subsequently, visualization of these genes including PDGFC, GHR, PRLR and ITGA3 showed that these genes were associated with the PI3K-AKT signal pathway. Moreover, the experimental results showed that PRLR had a higher expression in KGN cells, and that knocking PRLR down suppressed cell viability and promoted apoptosis of KGN cells.

Conclusion

This study revealed the common pathological mechanisms between PCOS and RSA and explored the role of the PI3K-AKT signaling pathway in the two diseases, providing a new direction for the clinical treatment of PCOS and RSA.

Molecular subtype construction and prognosis model for stomach adenocarcinoma characterized by metabolism-related genes

Background

Metabolic reprogramming is implicated in cancer progression. However, the impact of metabolism-associated genes in stomach adenocarcinomas (STAD) has not been thoroughly reviewed. Herein, we characterized metabolic transcription-correlated STAD subtypes and evaluated a metabolic RiskScore for evaluation survival.

Method

Genes related to metabolism were gathered from previous study and metabolic subtypes were screened using ConsensusClusterPlus in TCGA-STAD and GSE66229 dataset. The ssGSEA, MCP-Count, ESTIMATE and CIBERSORT determined the immune infiltration. A RiskScore model was established using the WGCNA and LASSO Cox regression in the TCGA-STAD queue and verified in the GSE66229 datasets. RT-qPCR was employed to measure the mRNA expressions of genes in the model.

Result

Two metabolism-related subtypes (C1 and C2) of STAD were constructed on account of the expression profiles of 113 prognostic metabolism genes with different immune outcomes and apparently distinct metabolic characteristic. The overall survival (OS) of C2 subtype was shorter than that of C1 subtype. Four metabolism-associated genes in turquoise model, which closely associated with C2 subtype, were employed to build the RiskScore (MATN3, OSBPL1A, SERPINE1, CPNE8) in TCGA-train dataset. Patients developed a poorer prognosis if they had a high RiskScore than having a low RiskScore. The promising effect of RiskScore was verified in the TCGA-test, TCGA-STAD and GSE66229 datasets. The prediction reliability of the RiskScore was validated by time-dependent receiver operating characteristic curve (ROC) and nomogram. Moreover, samples with high RiskScore had an enhanced immune status and TIDE score. Moreover, MATN3, OSBPL1A, SERPINE1 and CPNE8 mRNA levels were all elevated in SGC7901 cells. Inhibition of OSBPL1A decreased SGC7901 cells invasion numbers.

Conclusion

This work provided a new perspective into heterogeneity in metabolism and its association with immune escape in STAD. RiskScore was considered to be a strong prognostic label that could help individualize the treatment of STAD patients.

Research and experimental verification on the mechanisms of cellular senescence in triple-negative breast cancer

Background

Triple-negative breast cancer (TNBC) is an aggressive breast cancer subtype with high heterogeneity, poor prognosis, and a low 10-year survival rate of less than 50%. Although cellular senescence displays extensive effects on cancer, the comprehensions of cellular senescence-related characteristics in TNBC patients remains obscure.

Method

Single-cell RNA sequencing (scRNA-seq) data were analyzed by Seurat package. Scores for cellular senescence-related pathways were computed by single-sample gene set enrichment analysis (ssGSEA). Subsequently, unsupervised consensus clustering was performed for molecular cluster identification. Immune scores of patients in The Cancer Genome Atlas (TCGA) dataset and associated immune cell scores were calculated using Estimation of STromal and Immune cells in MAlignantTumours using Expression data (ESTIMATE) and Microenvironment Cell Populations-counter (MCP-counter), Tumor Immune Estimation Resource (TIMER) and Estimating the Proportion of Immune and Cancer cells (EPIC) methods, respectively. Immunotherapy scores were assessed using TIDE. Furthermore, feature genes were identified by univariate Cox and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses; these were used to construct a risk model. Additionally, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and transwell assay were conducted for in vitro validation of hub genes.

Result

TNBC was classified into three subtypes based on cellular senescence-related pathways as clusters 1, 2, and 3. Specifically, cluster 1 showed the best prognosis, followed by cluster 2 and cluster 3. The levels of gene expression in cluster 2 were the lowest, whereas these were the highest in cluster 3. Moreover, clusters 1 and 3 showed a high degree of immune infiltration. TIDE scores were higher for cluster 3, suggesting that immune escape was more likely in patients with the cluster 3 subtype who were less likely to benefit from immunotherapy. Next, the TNBC risk model was constructed and validated. RT-qPCR revealed that prognostic risk genes (MMP28, ACP5 and KRT6A) were up-regulated while protective genes (CT83) were down-regulated in TNBC cell lines, validating the results of the bioinformatics analysis. Meanwhile, cellular experiments revealed that ACP5 could promote the migration and invasion abilities in two TNBC cell lines. Finally, we evaluated the validity of prognostic models for assessing TME characteristics and TNBC chemotherapy response.

Conclusion

In conclusion, these findings help to assess the efficacy of targeted therapies in patients with different molecular subtypes, have practical applications for subtype-specific treatment of TNBC patients, and provide information on prognostic factors, as well as guidance for the revelation of the molecular mechanisms by which senescence-associated genes influence TNBC progression.

Exploring the resistance mechanism of triple-negative breast cancer to paclitaxel through the scRNA-seq analysis

BACKGROUND

The triple negative breast cancer (TNBC) is the most malignant subtype of breast cancer with high aggressiveness. Although paclitaxel-based chemotherapy scenario present the mainstay in TNBC treatment, paclitaxel resistance is still a striking obstacle for cancer cure. So it is imperative to probe new therapeutic targets through illustrating the mechanisms underlying paclitaxel chemoresistance.

METHODS

The Single cell RNA sequencing (scRNA-seq) data of TNBC cells treated with paclitaxel at different points were downloaded from the Gene Expression Omnibus (GEO) database. The Seurat R package was used to filter and integrate the scRNA-seq expression matrix. Cells were further clustered by the FindClusters function, and the gene marker of each subset was defined by FindAllMarkers function. Then, the hallmark score of each cell was calculated by AUCell R package, the biological function of the highly expressed interest genes was analyzed by the DAVID database. Subsequently, we performed pseudotime analysis to explore the change patterns of drug resistance genes and SCENIC analysis to identify the key transcription factors (TFs). Finally, the inhibitors of which were also analyzed by the CTD database.

RESULTS

We finally obtained 6 cell subsets from 2798 cells, which were marked as AKR1C3+, WNT7A+, FAM72B+, RERG+, IDO1+ and HEY1+HCC1143 cell subsets, among which the AKR1C3+, IDO1+ and HEY1+ cell subsets proportions increased with increasing treatment time, and then were regarded as paclitaxel resistance subsets. Hallmark score and pseudotime analysis showed that these paclitaxel resistance subsets were associated with the inflammatory response, virus and interferon response activation. In addition, the gene regulatory networks (GRNs) indicated that 3 key TFs (STAT1, CEBPB and IRF7) played vital role in promoting resistance development, and five common inhibitors targeted these TFs as potential combination therapies of paclitaxel were identified.

CONCLUSION

In this study, we identified 3 paclitaxel resistance relevant IFs and their inhibitors, which offers essential molecular basis for paclitaxel resistance and beneficial guidance for the combination of paclitaxel in clinical TNBC therapy.

Identifying Hub Genes for Glaucoma based on Bulk RNA Sequencing Data and Multi-machine Learning Models

AIMS

The aims of this study were to determine hub genes in glaucoma through multiple machine learning algorithms.

BACKGROUND

Glaucoma has afflicted many patients for many years, with excessive pressure in the eye continuously damaging the nervous system and leading to severe blindness. An effective molecular diagnostic method is currently lacking.

OBJECTIVE

The present study attempted to reveal the molecular mechanism and gene regulatory network of hub genes in glaucoma, followed by an attempt to reveal the drug-gene-disease network regulated by hub genes.

METHODS

A microarray sequencing dataset (GSE9944) was obtained through the Gene Expression Omnibus database. The differentially expressed genes in Glaucoma were identified. Based on these genes, we constructed three machine learning models for feature training, Random Forest model (RF), Least absolute shrinkage and selection operator regression model (LASSO), and Support Vector Machines model (SVM). Meanwhile, Weighted Gene Co-Expression Network Analysis (WGCNA) was performed for GSE9944 expression profiles to identify Glaucoma-related genes. The overlapping genes in the four groups were considered as hub genes of Glaucoma. Based on these genes, we also constructed a molecular diagnostic model of Glaucoma. In this study, we also performed molecular docking analysis to explore the gene-drug network targeting hub genes. In addition, we evaluated the immune cell infiltration landscape in Glaucoma samples and normal samples by applying CIBERSORT method.

RESULTS

8 hub genes were determined: ATP6V0D1, PLEC, SLC25A1, HRSP12, PKN1, RHOD, TMEM158 and GSN. The diagnostic model showed excellent diagnostic performance (area under the curve=1). GSN might positively regulate T cell CD4 naïve as well as negatively regulate T cell regulation (Tregs). In addition, we constructed gene-drug networks in an attempt to explore novel therapeutic agents for Glaucoma.

CONCLUSION

Our results systematically determined 8 hub genes and established a molecular diagnostic model that allowed the diagnosis of Glaucoma. Our study provided a basis for future systematic studies of Glaucoma pathogenesis.