Development of an Autophagy-Related Gene Prognostic Model and Nomogram for Estimating Renal Clear Cell Carcinoma Survival

m7G-related lncRNAs are potential biomarkers for predicting prognosis and immune responses in patients with oral squamous cell carcinoma

Among head and neck cancers, oral squamous cell carcinoma (OSCC) is the most common malignant tumor. N-7-methylguanosine (m7G) and lncRNAs are both related to the development and progression of tumors. Therefore, this study aims to explore and establish the prognostic signal of OSCC based on m7G-related lncRNAs. In this study, RNA sequencing transcriptome data of OSCC patients were downloaded from The Cancer Genome Atlas (TCGA) database. Therefore, m7G-related lncRNAs were identified as differentially expressed in OSCC. Then, univariate Cox regression analysis and LASSO regression analysis were conducted to evaluate the prognostic significance of differentially expressed lncRNAs. Consequently, the abovementioned lncRNAs were assigned five OSCC patient risk scores, with high-risk and low-risk patients assigned to each group. Different signaling pathways were significantly enriched between the two groups as determined by set enrichment analysis (GSEA). Multivariate Cox regression analysis confirmed the factors used to construct the nomogram model. Then, the prognosis of the nomogram model was evaluated. Consequently, high-risk individuals had higher immune infiltration levels. According to the results of a study that evaluated the sensitivity of different risk subgroups to antitumour drugs, the high-risk group had a high sensitivity to doxorubicin. By performing real-time polymerase chain reaction (RT‒PCR), we verified the expression of these five m7G lncRNAs. Therefore, the model based on five m7G-related lncRNAs was able to predict the overall survival rates of OSCC patients and guide their treatment. It can also spur new ideas about how to prevent and treat OSCC.

Molecular Subtyping Based on Cuproptosis-Related Genes and Characterization of Tumor Microenvironment Infiltration in Kidney Renal Clear Cell Carcinoma

The incidence of kidney renal clear cell carcinoma (KIRC) is rising worldwide, and the prognosis is poor. Cuproptosis is a new form of cell death that is dependent on and regulated by copper ions. The relationship between cuproptosis and KIRC remains unclear. In the current study, changes in cuproptosis-related genes (CRGs) in TCGA-KIRC transcriptional datasets were characterized, and the expression patterns of these genes were analyzed. We identified three main molecular subtypes and discovered that multilayer CRG changes were associated with patient clinicopathological traits, prognosis, elesclomol sensitivity, and tumor microenvironment (TME) cell infiltration characteristics. Then, a CRG score was created to predict overall survival (OS). The CRG score was found to be strongly linked to the TME. These findings may help elucidate the roles of CRGs in KIRC, potentially enhancing understanding of cuproptosis and supporting the development of more effective immunotherapy strategies.

Identifying the role of transient receptor potential channels (TRPs) in kidney renal clear cell carcinoma and their potential therapeutic significances using genomic and transcriptome analyses

Kidney renal clear cell carcinoma (KIRC) is among the major causes of cancer-caused mortality around the world. Transient receptor potential channels (TRPs), due to their role in various human diseases, might become potential drug targets in cancer. The mRNA expression, copy number variation, single-nucleotide variation, prognostic values, drug sensitivity, and pathway regulation of TRPs were studied across cancer types. The ArrayExpress and The Cancer Genome Atlas (TCGA) databases were used to retrieve KIRC samples. Simultaneously, training, internal, and external cohorts were grouped. In KIRC, a prognostic signature with superior survival prediction in contrast with other well-established signatures was created after a stepwise screening of optimized genes linked to TRPs using univariate Cox, weighted gene co-expression network analysis, multivariate Cox, and least absolute shrinkage and selection operator regression analyses. Subsequent to the determination of risk levels, the variations in the expression of immune checkpoint genes, tumor mutation burden, and immune subtypes and response between low-risk and high-risk subgroups were studied using a variety of bioinformatics algorithms, including ESTIMATE, XCELL, EPIC, CIBERSORT-ABS, CIBERSORT, MCPCOUNTER, TIMER, and QUANTISEQ. Gene set enrichment analysis helped in the identification of abnormal pathways across the low- and high-risk subgroups. Besides, high-risk KIRC patients might benefit from ABT888, AZD6244, AZD7762, Bosutinib, Camptothecin, CI1040, JNK inhibitor VIII, KU55933, Lenalidomide, Nilotinib, PLX4720, RO3306, Vinblastine, and ZM.447439; however, low-risk populations might benefit from Bicalutamide, FH535, and OSI906. Finally, calibration curves were used to validate the nomogram with a satisfactory predictive survival probability. In conclusion, this research provides useful insight that can aid and guide clinical practice and scientific research.

Identification of Diagnostic Biomarkers, Immune Infiltration Characteristics, and Potential Compounds in Rheumatoid Arthritis

Aims

This study is aimed at investigating the pathogenesis of rheumatoid arthritis (RA) by identifying key biomarkers, associated immune infiltration, and small-molecule compounds using bioinformatic analysis.

Methods

Six datasets were obtained from the Gene Expression Omnibus database, and the batch effect was adjusted. Functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to analyse differentially expressed genes (DEGs). Furthermore, candidate small-molecule drugs associated with RA were selected from the Connectivity Map (CMap) database. The least absolute shrinkage and selection operator regression, support vector machine recursive feature elimination, and multivariate logistic regression analyses were performed on DEGs to screen for RA diagnostic markers. The receiver operating characteristic curve, concordance index, and GiViTi calibration band were the metrics used to assess the diagnostic markers of RA identified in this analysis. The single-sample gene set enrichment analysis was performed to calculate the scores of infiltrating immune cells and evaluate the activities of immune-related pathways. Finally, the correlation between screening markers and RA diagnosis was determined.

Results

A total of 227 DEGs were identified. Functional enrichment analysis and KEGG revealed that DEGs were enriched by the immune response. CMap analysis identified 11 small-molecule compounds with therapeutic potential for RA. In gene expression, the activities of 13 immune cells and 12 immune-related pathways significantly differed between patients with RA and healthy controls. DPYSL3 and SPP1 had the potential to diagnose RA. SPP1 expression was positively correlated with DPYSL3 in 11 immune cells and 10 immune-related pathways.

Conclusion

This study comprehensively analysed DEGs and immune infiltration and screened for potential diagnostic markers and small-molecule compounds of RA.

Identification of Hypoxia-Immune-Related Gene Signatures and Construction of a Prognostic Model in Kidney Renal Clear Cell Carcinoma

Introduction: Kidney renal clear cell carcinoma (KIRC), a kind of malignant disease, is a severe threat to public health. Tracking the information of tumor progression and conducting a related dynamic prognosis model are necessary for KIRC. It is crucial to identify hypoxia-immune-related genes and construct a prognostic model due to immune interaction and the influence of hypoxia in the prognosis of patients with KIRC. Methods: The hypoxia and immune status of KIRC patients were identified by utilizing t-SNE and ImmuCellAI for gene expression data. COX and Lasso regression were used to identify some hypoxia-immune-related signature genes and further construct a prognostic risk model based on these genes. Internal and external validations were also conducted to construct a prognostic model. Finally, some potentially effective drugs were screened by the CMap dataset. Results: We found that high-hypoxia and low-immune status tend to induce poor overall survival (OS). Six genes, including PLAUR, UCN, PABPC1L, SLC16A12, NFE2L3, and KCNAB1, were identified and involved in our hypoxia-immune-related prognostic risk model. Internal verification showed that the area under the curve (AUC) for the constructed models for 1-, 3-, 4-, and 5-year OS were 0.768, 0.754, 0.775, and 0.792, respectively. For the external verification, the AUC for 1-, 3-, 4-, and 5-year OS were 0.768, 0.739, 0.763, and 0.643 respectively. Furthermore, the decision curve analysis findings demonstrated excellent clinical effectiveness. Finally, we found that four drugs (including vorinostat, fludroxycortide, oxolinic acid, and flutamide) might be effective and efficient in alleviating or reversing the status of severe hypoxia and poor infiltration of immune cells. Conclusion: Our constructed prognostic model, based on hypoxia-immune-related genes, has excellent effectiveness and clinical application value. Moreover, some small-molecule drugs are screened to alleviate severe hypoxia and poor infiltration of immune cells.

Identification of a Four-Gene Signature for Diagnosing Paediatric Sepsis

Aim

Early diagnosis of paediatric sepsis is crucial for the proper treatment of children and reduction of hospitalization and mortality. Biomarkers are a convenient and effective method for diagnosing any disease. However, huge differences among the studies reporting biomarkers for diagnosing sepsis have limited their clinical application. Therefore, in this study, we aimed to evaluate the diagnostic value of key genes involved in paediatric sepsis based on the data of the Gene Expression Omnibus database.

Methods

We used the GSE119217 dataset to identify differentially expressed genes (DEGs) between patients with and without paediatric sepsis. The most relevant gene modules of paediatric sepsis were screened through the weighted gene coexpression network analysis (WGCNA). Common genes (CGs) were found between DEGs and WGCNA. Genes with a potential diagnostic value in paediatric sepsis were selected from the CGs using least absolute shrinkage and selection operator regression and support vector machine recursive feature elimination. The principal component analysis, receiver operating characteristic curves, and C-index were used to verify the diagnostic value of the identified genes in six other independent sepsis datasets. Subsequently, a meta-analysis of the selected genes was performed to evaluate the value of these genes as biomarkers in paediatric sepsis.

Results

A total of 41 CGs were selected from the GSE119217 dataset. A four-gene signature composed of ANXA3, CD177, GRAMD1C, and TIGD3 effectively distinguished patients with paediatric sepsis from those in the control group. The signature was verified using six other independent datasets. In addition, the meta-analysis results showed that the pooled sensitivity, specificity, and area under the curve values were 1.00, 0.98, and 1.00, respectively.

Conclusion

The four-gene signature can be used as new biomarkers to distinguish patients with paediatric sepsis from healthy individuals.

The Construction of a Prognostic Model Based on a Peptidyl Prolyl Cis-Trans Isomerase Gene Signature in Hepatocellular Carcinoma

Objective: The aim of the present study was to construct a prognostic model based on the peptidyl prolyl cis–trans isomerase gene signature and explore the prognostic value of this model in patients with hepatocellular carcinoma.

Methods: The transcriptome and clinical data of hepatocellular carcinoma patients were downloaded from The Cancer Genome Atlas and the International Cancer Genome Consortium database as the training set and validation set, respectively. Peptidyl prolyl cis–trans isomerase gene sets were obtained from the Molecular Signatures Database. The differential expression of peptidyl prolyl cis–trans isomerase genes was analyzed by R software. A prognostic model based on the peptidyl prolyl cis–trans isomerase signature was established by Cox, Lasso, and stepwise regression methods. Kaplan–Meier survival analysis was used to evaluate the prognostic value of the model and validate it with an independent external data. Finally, nomogram and calibration curves were developed in combination with clinical staging and risk score.

Results: Differential gene expression analysis of hepatocellular carcinoma and adjacent tissues showed that there were 16 upregulated genes. A prognostic model of hepatocellular carcinoma was constructed based on three gene signatures by Cox, Lasso, and stepwise regression analysis. The Kaplan–Meier curve showed that hepatocellular carcinoma patients in high-risk score group had a worse prognosis (p < 0.05). The receiver operating characteristic curve revealed that the area under curve values of predicting the survival rate at 1, 2, 3, 4, and 5 years were 0.725, 0.680, 0.644, 0.630, and 0.639, respectively. In addition, the evaluation results of the model by the validation set were basically consistent with those of the training set. A nomogram incorporating clinical stage and risk score was established, and the calibration curve matched well with the diagonal.

Conclusion: A prognostic model based on 3 peptidyl prolyl cis–trans isomerase gene signatures is expected to provide reference for prognostic risk stratification in patients with hepatocellular carcinoma.