Identification of co-expression modules and pathways correlated with osteosarcoma and its metastasis

The prognostic value of ubiquitin/ubiquitin-like-related genes along with immune cell infiltration and clinicopathological features in osteosarcoma

BACKGROUND

Ubiquitin/ubiquitin-like (Ub/UBL)-related genes have been reported to be associated with the survival of osteosarcoma patients but have not yet been systematically explored.

METHODS

The prognostic value of Ub/UBL-related genes, immune cell infiltration and clinicopathological features of patients were explored by Cox and LASSO regression analyses. A prognostic model was established and then validated in the GSE21257 dataset. The differential expression of hub genes in osteosarcoma was confirmed by qRT-PCR, western blotting and immunohistochemistry.

RESULTS

Tripartite Motif Containing 8 (TRIM8) and Ubiquitin Like With PHD And Ring Finger Domains 2 (UHRF2) were screened as genes with prognostic value in osteosarcoma. Kaplan-Meier analysis and scatter plots indicated that patients in the high gene significance score group tended to have a worse prognosis. The concordance index, calibration analysis and receiver operating characteristic analysis suggested that the model had good prediction accuracy and high sensitivity and specificity. Decision curve analysis revealed that patients could obtain greater net benefit from this model. Functional analyses of the differentially expressed genes indicated that they were involved in important functions and pathways. TRIM8 and UHRF2 were confirmed to be highly expressed in osteosarcoma cell lines and tissues.

CONCLUSIONS

TRIM8 and UHRF2 are potential prognostic genes in osteosarcoma, and these results provide insights into the roles of these genes and their implications for patient outcomes.

Using machine learning methods to study the tumour microenvironment and its biomarkers in osteosarcoma metastasis

Background

The long-term prognosis for patients with osteosarcoma (OS) metastasis remains unfavourable, highlighting the urgent need for research that explores potential biomarkers using innovative methodologies.

Methods

This study explored potential biomarkers for OS metastasis by analysing data from the Cancer Genome Atlas Program (TCGA) and Gene Expression Omnibus (GEO) databases. The synthetic minority oversampling technique (SMOTE) was employed to tackle class imbalances, while genes were selected using four feature selection algorithms (Monte Carlo feature selection [MCFS], Borota, minimum-redundancy maximum-relevance [mRMR], and light gradient-boosting machine [LightGBM]) based on the gene expression matrix. Four machine learning (ML) algorithms (support vector machine [SVM], extreme gradient boosting [XGBoost], random forest [RF], and k-nearest neighbours [kNN]) were utilized to determine the optimal number of genes for building the model. Interpretable machine learning (IML) was applied to construct prediction networks, revealing potential relationships among the selected genes. Additionally, enrichment analysis, survival analysis, and immune infiltration were performed on the featured genes.

Results

In DS1, DS2, and DS3, the IML algorithm identified 53, 45, and 46 features, respectively. Using the merged gene set, we obtained a total of 79 interpretable prediction rules for OS metastasis. We subsequently conducted an in-depth investigation on 39 crucial molecules associated with predicting OS metastasis, elucidating their roles within the tumour microenvironment. Importantly, we found that certain genes act as both predictors and differentially expressed genes. Finally, our study unveiled statistically significant differences in survival between the high and low expression groups of TRIP4, S100A9, SELL and SLC11A1, and there was a certain correlation between these genes and 22 various immune cells.

Conclusions

The biomarkers discovered in this study hold significant implications for personalized therapies, potentially enhancing the clinical prognosis of patients with OS.

Analysis of translesion polymerases in colorectal cancer cells following cetuximab treatment: A network perspective

INTRODUCTION

Adaptive mutagenesis observed in colorectal cancer (CRC) cells upon exposure to EGFR inhibitors contributes to the development of resistance and recurrence. Multiple investigations have indicated a parallel between cancer cells and bacteria in terms of exhibiting adaptive mutagenesis. This phenomenon entails a transient and coordinated escalation of error-prone translesion synthesis polymerases (TLS polymerases), resulting in mutagenesis of a magnitude sufficient to drive the selection of resistant phenotypes.

METHODS

In this study, we conducted a comprehensive pan-transcriptome analysis of the regulatory framework within CRC cells, with the objective of identifying potential transcriptome modules encompassing certain translesion polymerases and the associated transcription factors (TFs) that govern them. Our sampling strategy involved the collection of transcriptomic data from tumors treated with cetuximab, an EGFR inhibitor, untreated CRC tumors, and colorectal-derived cell lines, resulting in a diverse dataset. Subsequently, we identified co-regulated modules using weighted correlation network analysis with a minKMEtostay threshold set at 0.5 to minimize false-positive module identifications and mapped the modules to STRING annotations. Furthermore, we explored the putative TFs influencing these modules using KBoost, a kernel PCA regression model.

RESULTS

Our analysis did not reveal a distinct transcriptional profile specific to cetuximab treatment. Moreover, we elucidated co-expression modules housing genes, for example, POLK, POLI, POLQ, REV1, POLN, and POLM. Specifically, POLK, POLI, and POLQ were assigned to the "blue" module, which also encompassed critical DNA damage response enzymes, for example. BRCA1, BRCA2, MSH6, and MSH2. To delineate the transcriptional control of this module, we investigated associated TFs, highlighting the roles of prominent cancer-associated TFs, such as CENPA, HNF1A, and E2F7.

CONCLUSION

We found that translesion polymerases are co-regulated with DNA mismatch repair and cell cycle-associated factors. We did not, however, identified any networks specific to cetuximab treatment indicating that the response to EGFR inhibitors relates to a general stress response mechanism.

Metastasis-Related Signature for Clinically Predicting Prognosis and Tumor Immune Microenvironment of Osteosarcoma Patients

Osteosarcoma is the most prevalent clinical malignant bone tumor in adolescents. The prognosis of metastatic osteosarcoma is still very poor. The aim of our study was to investigate the clinical diagnosis and prognostic significance of metastasis related genes (MRGs) in patients with osteosarcoma. Clinical information and RNA sequencing data with osteosarcoma patients were obtained and set as the training set from UCSC databases. GSE21257 were downloaded and chosen as the verification cohort. An eight gene metastasis related risk signature including MYC, TAC4, ABCA4, GADD45GIP1, TNFRSF21, HERC5, MAGEA11, and PDE1B was built to predict the overall survival of osteosarcoma patients. Based on risk assessments, patients were classified into high- and low-risk groups. The high-risk patients had higher risk score and shorter survival time. ROC curves revealed that this risk signature can accurately predict survival times of osteosarcoma patients at the 1-, 2-, 3-, 4- and 5- year. GSEA revealed that MYC targets, E2F targets, mTORC1 signaling, Wnt /β-catenin signaling and cell cycle were upregulated, and cell adhesion molecules, and primary immunodeficiency were decreased in high-risk group. MRGs were highly linked with the tumor immune microenvironment and ICB response. These results identified that MRGs as a novel prognostic and diagnostic biomarker in osteosarcoma.

Machine learning algorithms assisted identification of post-stroke depression associated biological features

Objectives

Post-stroke depression (PSD) is a common and serious psychiatric complication which hinders functional recovery and social participation of stroke patients. Stroke is characterized by dynamic changes in metabolism and hemodynamics, however, there is still a lack of metabolism-associated effective and reliable diagnostic markers and therapeutic targets for PSD. Our study was dedicated to the discovery of metabolism related diagnostic and therapeutic biomarkers for PSD.

Methods

Expression profiles of GSE140275, GSE122709, and GSE180470 were obtained from GEO database. Differentially expressed genes (DEGs) were detected in GSE140275 and GSE122709. Functional enrichment analysis was performed for DEGs in GSE140275. Weighted gene co-expression network analysis (WGCNA) was constructed in GSE122709 to identify key module genes. Moreover, correlation analysis was performed to obtain metabolism related genes. Interaction analysis of key module genes, metabolism related genes, and DEGs in GSE122709 was performed to obtain candidate hub genes. Two machine learning algorithms, least absolute shrinkage and selection operator (LASSO) and random forest, were used to identify signature genes. Expression of signature genes was validated in GSE140275, GSE122709, and GSE180470. Gene set enrichment analysis (GSEA) was applied on signature genes. Based on signature genes, a nomogram model was constructed in our PSD cohort (27 PSD patients vs. 54 controls). ROC curves were performed for the estimation of its diagnostic value. Finally, correlation analysis between expression of signature genes and several clinical traits was performed.

Results

Functional enrichment analysis indicated that DEGs in GSE140275 enriched in metabolism pathway. A total of 8,188 metabolism associated genes were identified by correlation analysis. WGCNA analysis was constructed to obtain 3,471 key module genes. A total of 557 candidate hub genes were identified by interaction analysis. Furthermore, two signature genes (SDHD and FERMT3) were selected using LASSO and random forest analysis. GSEA analysis found that two signature genes had major roles in depression. Subsequently, PSD cohort was collected for constructing a PSD diagnosis. Nomogram model showed good reliability and validity. AUC values of receiver operating characteristic (ROC) curve of SDHD and FERMT3 were 0.896 and 0.964. ROC curves showed that two signature genes played a significant role in diagnosis of PSD. Correlation analysis found that SDHD (r = 0.653, P < 0.001) and FERM3 (r = 0.728, P < 0.001) were positively related to the Hamilton Depression Rating Scale 17-item (HAMD) score.

Conclusion

A total of 557 metabolism associated candidate hub genes were obtained by interaction with DEGs in GSE122709, key modules genes, and metabolism related genes. Based on machine learning algorithms, two signature genes (SDHD and FERMT3) were identified, they were proved to be valuable therapeutic and diagnostic biomarkers for PSD. Early diagnosis and prevention of PSD were made possible by our findings.

PKIB involved in the metastasis and survival of osteosarcoma

Osteosarcoma is frequently metastasized at the time of diagnosis in patients. However, the underlying mechanism of osteosarcoma metastasis remains poorly understood. In this study, we evaluated DNA methylation profiles combined with gene expression profiles of 21 patients with metastatic osteosarcoma and 64 patients with non-metastatic osteosarcoma from TARGET database and identified PKIB and AIM2 as hub genes related to the metastasis of osteosarcoma. To verify the effects of PKIB on migration and invasion of osteosarcoma, we performed wound-healing assay and transwell assay. The results showed that PKIB significantly inhibited the migration and invasion of osteosarcoma cells, and the Western blot experiments showed that the protein level of E-cad was upregulated and of VIM was downregulated in 143-B cell recombinant expression PKIB. These results indicate that PKIB inhibit the metastasis of osteosarcoma. CCK-8 assay results showed that PKIB promote the proliferation of osteosarcoma. In addition, the Western blot results showed that the phosphorylation level of Akt was upregulated in 143-B cells overexpressing PKIB, indicating that PKIB promotes the proliferation of osteosarcoma probably through signaling pathway that Akt involved in. These results give us clues that PKIB was a potential target for osteosarcoma therapy. Furthermore, combined clinical profiles analysis showed that the expression of AIM2- and PKIB- related risk scores was significantly related to the overall survival of patients with osteosarcoma. Thus, we constructed a nomogram based on AIM2 and PKIB expression–related risk scores for osteosarcoma prognostic assessment to predict the 1-, 2-, 3-, and 5-year overall survival rate of patients with metastatic osteosarcoma, assisting clinicians in the diagnosis and treatment of metastatic osteosarcoma.

Differential Gene Expression and Weighted Correlation Network Dynamics in High-Throughput Datasets of Prostate Cancer

Precision oncology is an absolute need today due to the emergence of treatment resistance and heterogeneity among cancerous profiles. Target-propelled cancer therapy is one of the treasures of precision oncology which has come together with substantial medical accomplishment. Prostate cancer is one of the most common cancers in males, with tremendous biological heterogeneity in molecular and clinical behavior. The spectrum of molecular abnormalities and varying clinical patterns in prostate cancer suggest substantial heterogeneity among different profiles. To identify novel therapeutic targets and precise biomarkers implicated with prostate cancer, we performed a state-of-the-art bioinformatics study, beginning with analyzing high-throughput genomic datasets from The Cancer Genome Atlas (TCGA). Weighted gene co-expression network analysis (WGCNA) suggests a set of five dysregulated hub genes (MAF, STAT6, SOX2, FOXO1, and WNT3A) that played crucial roles in biological pathways associated with prostate cancer progression. We found overexpressed STAT6 and SOX2 and proposed them as candidate biomarkers and potential targets in prostate cancer. Furthermore, the alteration frequencies in STAT6 and SOX2 and their impact on the patients' survival were explored through the cBioPortal platform. The Kaplan-Meier survival analysis suggested that the alterations in the candidate genes were linked to the decreased overall survival of the patients. Altogether, the results signify that STAT6 and SOX2 and their genomic alterations can be explored in therapeutic interventions of prostate cancer for precision oncology, utilizing early diagnosis and target-propelled therapy.

Transcriptomic Analysis of Canine Osteosarcoma from a Precision Medicine Perspective Reveals Limitations of Differential Gene Expression Studies

Despite significant advances in cancer diagnosis and treatment, osteosarcoma (OSA), an aggressive primary bone tumor, has eluded attempts at improving patient survival for many decades. The difficulty in managing OSA lies in its extreme genetic complexity, drug resistance, and heterogeneity, making it improbable that a single-target treatment would be beneficial for the majority of affected individuals. Precision medicine seeks to fill this gap by addressing the intra- and inter-tumoral heterogeneity to improve patient outcome and survival. The characterization of differentially expressed genes (DEGs) unique to the tumor provides insight into the phenotype and can be useful for informing appropriate therapies as well as the development of novel treatments. Traditional DEG analysis combines patient data to derive statistically inferred genes that are dysregulated in the group; however, the results from this approach are not necessarily consistent across individual patients, thus contradicting the basis of precision medicine. Spontaneously occurring OSA in the dog shares remarkably similar clinical, histological, and molecular characteristics to the human disease and therefore serves as an excellent model. In this study, we use transcriptomic sequencing of RNA isolated from primary OSA tumor and patient-matched normal bone from seven dogs prior to chemotherapy to identify DEGs in the group. We then evaluate the universality of these changes in transcript levels across patients to identify DEGs at the individual level. These results can be useful for reframing our perspective of transcriptomic analysis from a precision medicine perspective by identifying variations in DEGs among individuals.

The Identification of Childhood Asthma Progression-Related lncRNAs and mRNAs Suitable as Biomarkers Using Weighted Gene Coexpression Network Analysis

Background

Asthma is a common chronic respiratory disease in children, seriously affecting children's health and growth. This bioinformatics study aimed to identify potential RNA candidates closely associated with childhood asthma development within current gene databases.

Methods

GSE65204 and GSE19187 datasets were screened and downloaded from the NCBI GEO database. Differentially expressed long noncoding RNAs (DE-lncRNAs) and mRNAs (DE-mRNAs) were identified using the Bioconductor limma package in R, and these DE-mRNAs were used to perform biological process (BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Thereafter, weighted gene coexpression network analysis (WGCNA) was utilized to screen the modules directly related to childhood asthma, and a coexpression network of DE-lncRNAs and DE-mRNAs was built. Finally, principal component analysis (PCA) was performed.

Results

In total, 7 DE-lncRNAs and 1060 DE-mRNAs, as well as 7 DE-lncRNAs and 1027 DE-mRNAs, were identified in GSE65204 and GSE19187, respectively. After comparison, 336 overlapping genes had the same trend of expression, including 2 overlapped DE-lncRNAs and 334 overlapped DE-mRNAs. These overlapped DE-mRNAs were enriched in 28 BP and 12 KEGG pathways. Eleven modules were obtained in GSE65204, and it was found that the purple, black, and yellow modules were significantly positively correlated with asthma development. Subsequently, a coexpression network including 63 DE-mRNAs and 2 DE-lncRNAs was built, and five KEGG pathways, containing 8 genes, were found to be directly associated with childhood asthma. The PCA further verified these results.

Conclusion

LncRNAs LINC01559 and SNHG8 and mRNAs VWF, LAMB3, LAMA4, CAV1, ALDH1A3, SMOX, GNG4, and PPARG were identified as biomarkers associated with the progression of childhood asthma.

Coexpression Module Construction by Weighted Gene Coexpression Network Analysis and Identify Potential Prognostic Markers of Breast Cancer

Background: Breast cancer (BC) is a malignant tumor with the highest morbidity among women, disrupting millions of their lives worldwide each year. However, the molecular mechanisms underlying remain unclear. Methods: The RNA-Sequencing and clinical data of BC patients from The Cancer Genome Atlas (TCGA) database were analyzed by weighted gene coexpression network analysis (WGCNA). Additionally, coexpressed modules were used to detect their correlation with the clinical traits of BC. Next, nodes of the most significant coexpression modules were used for Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, mRNA-lncRNA coexpression network and survival analyses. Results: In total, 2056 differentially expressed mRNAs (DEmRNAs) and 297 differentially expressed lncRNAs (DElncRNAs) were identified and subjected to WGCNA analysis, and 12 coexpression modules were generated. The top five significant modules (turquoise, green, red, brown, and blue modules) were related to one or more clinical traits of BC. In particular, the turquoise and green modules were chosen for further analysis. Next, by lncRNA-mRNA coexpression analysis of the turquoise and green modules, 12 DEmRNAs and 2 DElncRNAs were identified as hub nodes. The lncRNA-associated mRNAs of the networks were commonly related to several cancer-related pathways. Moreover, these networks also revealed central roles for RP11-389C8.2 and TGFBR2 in the turquoise module and MYLK, KIT, and RP11-394O4.5 in the green module. Furthermore, 16 DEmRNAs and 3 DElncRNAs in these two modules were significantly correlated with the overall survival of BC patients. Conclusions: The authors' study identified some prognostic biomarkers that might play important roles in the development and treatment of BC. In particular, lncRNAs AC016995.3, RP1-193H18.2, and RP11-166D19.1 were novel biomarkers for BC.

Metastasis-associated gene MAPK15 promotes the migration and invasion of osteosarcoma cells via the c-Jun/MMPs pathway

Osteosarcoma (OS) is the most common and destructive primary bone malignancy to affect children and adolescents. Metastases remain the primary cause of death in patients with OS. In the present study, weight gene co-expressed network analysis (WGCNA) and differentially-expressed gene analysis were used to identify key genes associated with the metastasis of OS. Reverse transcription-quantitative PCR and immunohistochemical staining were then used to detect the expression levels of these key genes in OS tissues, and to determine the hub genes of interest. Wound-healing and transwell assays, in addition to a lung metastasis model, were used to detect the effects of the hub genes on OS cell proliferation and metastasis in vitro and in vivo. Using WGCNA and differential expression analysis, deleted in lung and esophageal cancer protein 1 (DLEC1), Forkhead box J1 (FOXJ1) and mitogen-activated protein kinase 15 (MAPK15) were predicted to be key metastasis-associated genes, and highly expressed in metastatic OS tissues; among them, the protein and mRNA expression levels of MAPK15 were most significantly increased in our OS tissues from patients who exhibited metastases at diagnosis, and thus MAPK15 was determined to be a metastasis-associated hub gene to further study. Furthermore, inhibiting MAPK15 expression significantly decreased OS cell metastasis in vitro and in vivo, as well as suppressing c-Jun/matrix metalloproteinase (MMP)-associated pathways. Overexpression of MAPK15 activated the c-Jun/MMPs pathway and promoted OS cell metastasis, while inhibition of c-Jun blocked this effect. Taken together, MAPK15 was indicated to be an OS metastasis-associated gene, and was confirmed to promote the migration and invasion of OS cells via the c-Jun/MMP pathway. MAPK15 may therefore be an effective target for the treatment of OS.

Key genes with prognostic values in suppression of osteosarcoma metastasis using comprehensive analysis

BACKGROUND

Osteosarcoma is a primary malignant tumor originating from mesenchymal tissue, with a poor distant metastasis prognosis. The molecular mechanisms of osteosarcoma metastasis are extremely complicated.

METHODS

A public data series (GSE21257) was used to identify differentially expressed genes (DEGs) in osteosarcoma patients that did, or did not, develop metastases. Functional enrichment analysis, a protein-protein interaction network, and survival analysis of DEGs were performed. DEGs with a prognostic value were considered as candidate genes and their functional predictions, different expression in normal and malignant tissues, and immune infiltration were analyzed.

RESULTS

The DEGs were mainly enriched in the immune response. Three candidate genes (ALOX5AP, CD74, and FCGR2A) were found, all of which were expressed at higher levels in lungs and lymph nodes than in matched cancer tissues and were probably expressed in the microenvironment.

CONCLUSIONS

Candidate genes can help us understand the molecular mechanisms underlying osteosarcoma metastasis and provide targets for future research.

Characterization of long non-coding RNA and messenger RNA profiles in laryngeal cancer by weighted gene co-expression network analysis

Laryngeal cancer (LC) is a malignant tumor in the head and neck region. It was recently elucidated that long non-coding RNAs (lncRNAs) participate in the pathogenesis of LC. However, the detailed mechanism of lncRNA in LC and whether long non-coding RNAs serve as effective biomarkers remains unclear. Ribonucleic acid (RNA) sequence data of LC and 11 patient clinical traits were extracted from The Cancer Genome Atlas (TCGA) database and analyzed by weighted gene co-expression network analysis (WGCNA). A total of 9 co-expression modules were identified. The co-expression Pink module significantly correlated with four clinical traits, including history of smoking, lymph node count, tumor status, and the success of follow-up treatment. Based on the co-expression Pink module, lncRNA-microRNA (miRNA)-messenger RNA (mRNA) and lncRNA-RNA binding protein-mRNA networks were constructed. We found that 8 lncRNAs significantly impacted overall survival (OS) in LC patients. These identified lncRNA and hub gene biomarkers were also validated in multiple LC cells in vitro via qPCR. Taken together, this study provided the framework of co-expression gene modules of LC and identified some important biomarkers in LC development and disease progression.