An experimental study of the intrinsic stability of random forest variable importance measures

Integrating machine learning and single-cell transcriptomic analysis to identify potential biomarkers and analyze immune features of ischemic stroke

This study employs machine learning and single-cell transcriptome sequencing (scRNA-seq) analysis to unearth novel biomarkers and delineate the immune characteristics of ischemic stroke (IS), thereby contributing fresh insights into IS treatment strategies.Our research leverages gene expression data sourced from the GEO database. We undertake weighted gene co-expression network analysis (WGCNA) to filter pertinent genes and subsequently employ machine learning algorithms for the identification of feature genes. Concurrently, we rigorously execute quality control measures, dimensionality reduction techniques, and cell annotation on the scRNA-seq data to pinpoint differentially expressed genes (DEGs). The identification of core genes, denoted as Hub genes, among the feature genes and DEGs, is achieved through meticulous overlapping analysis. We illuminate the immune characteristics of these Hub genes using a suite of analytical tools, encompassing CIBERSORT, MCPcounter, and pseudotemporal analysis, all based on immune cell annotations and single-cell transcriptome data.Subsequently, we harness the CMap database to prognosticate potential therapeutic drugs and scrutinize their associations with the identified Hub genes. Our findings unveil robust linkages between three pivotal Hub genes-namely, RNF13, VASP, and CD163-and specific immune cell types such as T cells and neutrophils. These Hub genes predominantly manifest in macrophages and microglial cells within the scRNA-seq immune cell population, exhibiting variances across different stages of cellular differentiation. In conclusion, this study unearths highly pertinent biomarkers for IS diagnosis and elucidates IS-induced immune infiltration characteristics, thus providing a firm foundation for a comprehensive exploration of potential immune mechanisms and the identification of novel therapeutic targets for IS.

Understanding of Wetting Mechanism Toward the Sticky Powder and Machine Learning in Predicting Granule Size Distribution Under High Shear Wet Granulation

The granulation of traditional Chinese medicine (TCM) has attracted widespread attention, there is limited research on the high shear wet granulation (HSWG) and wetting mechanisms of sticky TCM powders, which profoundly impact the granule size distribution (GSD). Here we investigate the wetting mechanism of binders and the influence of various parameters on the GSD of HSWG and establish a GSD prediction model. Permeability and contact angle experiments combined with molecular dynamics (MD) simulations were used to explore the wetting mechanism of hydroalcoholic solutions with TCM powder. Machine learning (ML) was employed to build a GSD prediction model, feature importance explained the influence of features on the predictive performance of the model, and correlation analysis was used to assess the influence of various parameters on GSD. The results show that water increases powder viscosity, forming high-viscosity aggregates, while ethanol primarily acted as a wetting agent. The contact angle of water on the powder bed was the largest and decreased with an increase in ethanol concentration. Extreme Gradient Boosting (XGBoost) outperformed other models in overall prediction accuracy in GSD prediction, the binder had the greatest impact on the predictions and GSD, adjusting the amount and concentration of adhesive can control the adhesion and growth of granules while the impeller speed had the least influence on granulation. The study elucidates the wetting mechanism and provides a GSD prediction model, along with the impact of material properties, formulation, and process parameters obtained, aiding the intelligent manufacturing and formulation development of TMC.

Integrated machine learning and bioinformatic analysis of mitochondrial-related signature in chronic rhinosinusitis with nasal polyps

Background

Chronic rhinosinusitis with nasal polyps (CRSwNP) is a prevalent inflammatory disorder affecting the upper respiratory tract. Recent studies have indicated an association between CRSwNP and mitochondrial metabolic disorder characterized by impaired metabolic pathways; however, the precise mechanisms remain unclear. This study aims to investigate the mitochondrial-related signature in individuals diagnosed with CRSwNP.

Methods

Through the integration of differentially expressed genes (DEGs) with the mitochondrial gene set, differentially expressed mitochondrial-related genes (DEMRGs) were identified. Subsequently, the hub DEMRGs were selected using 4 integrated machine learning algorithms. Immune and mitochondrial characteristics were estimated based on CIBERSORT and ssGSEA algorithms. Bioinformatic findings were confirmed through RT-qPCR, immunohistochemistry, and ELISA for nasal tissues, as well as Western blotting analysis for human nasal epithelial cells (hNECs). The relationship between hub DEMRGs and disease severity was assessed using Spearman correlation analysis.

Results

A total of 24 DEMRGs were screened, most of which exhibited lower expression levels in CRSwNP samples. Five hub DEMRGs (ALDH1L1, BCKDHB, CBR3, HMGCS2, and OXR1) were consistently downregulated in both the discovery and validation cohorts. The hub genes showed a high diagnostic performance and were positively correlated with the infiltration of M2 macrophages and resting mast cells. Experimental results confirmed that the 5 genes were downregulated at both the mRNA and protein levels within nasal polyp tissues. Finally, a significant and inverse relationship was identified between the expression levels of these genes and both the Lund-Mackay and Lund-Kennedy scores.

Conclusion

Our findings systematically unraveled 5 hub markers correlated with mitochondrial metabolism and immune cell infiltration in CRSwNP, suggesting their potential to be based to design diagnostic and therapeutic strategies for the disease.

Prognostic model incorporating immune checkpoint genes to predict the immunotherapy efficacy for lung adenocarcinoma: a cohort study integrating machine learning algorithms

This study aimed to develop and validate a nomogram based on immune checkpoint genes (ICGs) for predicting prognosis and immune checkpoint blockade (ICB) efficacy in lung adenocarcinoma (LUAD) patients. A total of 385 LUAD patients from the TCGA database and 269 LUAD patients in the combined dataset (GSE41272 + GSE50081) were divided into training and validation cohorts, respectively. Three different machine learning algorithms including random forest (RF), least absolute shrinkage and selection operator (LASSO) logistic regression analysis, and support vector machine (SVM) were employed to select the predictive markers from 82 ICGs to construct the prognostic nomogram. The X-tile software was used to stratify patients into high- and low-risk subgroups based on the nomogram-derived risk scores. Differences in functional enrichment and immune infiltration between the two subgroups were assessed using gene set variation analysis (GSVA) and various algorithms. Additionally, three lung cancer cohorts receiving ICB therapy were utilized to evaluate the ability of the model to predict ICB efficacy in the real world. Five ICGs were identified as predictive markers across all three machine learning algorithms, leading to the construction of a nomogram with strong potential for prognosis prediction in both the training and validation cohorts (all AUC values close to 0.800). The patients were divided into high- (risk score ≥ 185.0) and low-risk subgroups (risk score < 185.0). Compared to the high-risk subgroup, the low-risk subgroup exhibited enrichment in immune activation pathways and increased infiltration of activated immune cells, such as CD8 + T cells and M1 macrophages (P < 0.05). Furthermore, the low-risk subgroup had a greater likelihood of benefiting from ICB therapy and longer progression-free survival (PFS) than did the high-risk subgroup (P < 0.05) in the two cohorts receiving ICB therapy. A nomogram based on ICGs was constructed and validated to aid in predicting prognosis and ICB treatment efficacy in LUAD patients.

Machine learning and bioinformatics analysis of diagnostic biomarkers associated with the occurrence and development of lung adenocarcinoma

Objective

Lung adenocarcinoma poses a major global health challenge and is a leading cause of cancer-related deaths worldwide. This study is a review of three molecular biomarkers screened by machine learning that are not only important in the occurrence and progression of lung adenocarcinoma but also have the potential to serve as biomarkers for clinical diagnosis, prognosis evaluation and treatment guidance.

Methods

Differentially expressed genes (DEGs) were identified using comprehensive GSE1987 and GSE18842 gene expression databases. A comprehensive bioinformatics analysis of these DEGs was conducted to explore enriched functions and pathways, relative expression levels, and interaction networks. Random Forest and LASSO regression analysis techniques were used to identify the three most significant target genes. The TCGA database and quantitative polymerase chain reaction (qPCR) experiments were used to verify the expression levels and receiver operating characteristic (ROC) curves of these three target genes. Furthermore, immune invasiveness, pan-cancer, and mRNA-miRNA interaction network analyses were performed.

Results

Eighty-nine genes showed increased expression and 190 genes showed decreased expression. Notably, the upregulated DEGs were predominantly associated with organelle fission and nuclear division, whereas the downregulated DEGs were mainly associated with genitourinary system development and cell-substrate adhesion. The construction of the DEG protein-protein interaction network revealed 32 and 19 hub genes with the highest moderate values among the upregulated and downregulated genes, respectively. Using random forest and LASSO regression analyses, the hub genes were employed to identify three most significant target genes.TCGA database and qPCR experiments were used to verify the expression levels and ROC curves of these three target genes, and immunoinvasive analysis, pan-cancer analysis and mRNA-miRNA interaction network analysis were performed.

Conclusion

Three target genes identified by machine learning: BUB1B, CENPF, and PLK1 play key roles in LUAD development of lung adenocarcinoma.

MUC16 as a serum-based prognostic indicator of prometastatic gastric cancer

Metastatic gastric cancer (GC) presents significant clinical challenges due to its poor prognosis and limited treatment options. To address this, we conducted a targeted protein biomarker discovery study to identify markers predictive of metastasis in advanced GC (AGC). Serum samples from 176 AGC patients (T stage 3 or higher) were analyzed using the Olink Proteomics Target panels. Patients were retrospectively categorized into nonmetastatic, metastatic, and recurrence groups, and differential protein expression was assessed. Machine learning and gene set enrichment analysis (GSEA) methods were applied to discover biomarkers and predict prognosis. Four proteins (MUC16, CAIX, 5'-NT, and CD8A) were significantly elevated in metastatic GC patients compared to the control group. Additionally, GSEA indicated that the response to interleukin-4 and hypoxia-related pathways were enriched in metastatic patients. Random forest classification and decision-tree modeling showed that MUC16 could be a predictive marker for metastasis in GC patients. Additionally, ELISA validation confirmed elevated MUC16 levels in metastatic patients. Notably, high MUC16 levels were independently associated with metastatic progression in T3 or higher GC. These findings suggest the potential of MUC16 as a clinically relevant biomarker for identifying GC patients at high risk of metastasis.

Exploring the Roles of Key Mediators IKBKE and HSPA1A in Alzheimer's Disease and Hepatocellular Carcinoma through Bioinformatics Analysis

Recent studies have hinted at a potential link between Alzheimer's Disease (AD) and cancer. Thus, our study focused on finding genes common to AD and Liver Hepatocellular Carcinoma (LIHC), assessing their promise as diagnostic indicators and guiding future treatment approaches for both conditions. Our research utilized a broad methodology, including differential gene expression analysis, Weighted Gene Co-expression Network Analysis (WGCNA), gene enrichment analysis, Receiver Operating Characteristic (ROC) curves, and Kaplan-Meier plots, supplemented with immunohistochemistry data from the Human Protein Atlas (HPA) and machine learning techniques, to identify critical genes and significant pathways shared between AD and LIHC. Through differential gene expression analysis, WGCNA, and machine learning methods, we identified nine key genes associated with AD, which served as entry points for LIHC analysis. Subsequent analyses revealed IKBKE and HSPA1A as shared pivotal genes in patients with AD and LIHC, suggesting these genes as potential targets for intervention in both conditions. Our study indicates that IKBKE and HSPA1A could influence the onset and progression of AD and LIHC by modulating the infiltration levels of immune cells. This lays a foundation for future research into targeted therapies based on their shared mechanisms.

Hsa-miR-342-3p and hsa-miR-360 may be the key molecules that promote periodontitis in type 2 diabetes mellitus

Background

Periodontitis (PD) has been acknowledged as a complication associated with type 2 diabetes mellitus (T2DM). However, the precise mechanism through which T2DM fosters the development of PD remains elusive. Our objective is to elucidate the connection between these two conditions by conducting bioinformatics analysis.

Methods

In this study, we analyzed miRNA datasets pertaining to T2DM and PD sourced from GEO. Through differential expression analysis, we identified common differentially expressed miRNAs (DE-miRNAs) and subsequently analyzed the functional enrichment of these common DE-miRNAs. We further leveraged the PD transcriptome database to select DE-miRNA-targeted mRNAs and examined their association with immune infiltration. Finally, machine learning was used to further screen hub DE-miRNA-targeted mRNAs and validate our data in external datasets.

Results

Two common DE-miRNAs, namely hsa-miR-342-3p and hsa-miR-360, were identified from the miRNA datasets of PD and T2DM. Functional enrichment analysis indicated that these two common DE-miRNAs predominantly participate in Ras, PI3K-Akt, p53, and MAPK signaling pathways. Integration of the PD transcriptome dataset revealed a total of 21 DE-miRNA-targeted mRNAs in PD, with strong correlations observed with plasma cells and dendritic cells. Finally, three hub DE-miRNA-targeted mRNAs (hsa-miR-342-3p-/hsa-miR-360-RASAL2, hsa-miR-360-ENTPD1/PLXDC2) were identified. ENTPD1 exhibited a robust positive correlation with plasma cells and a negative correlation with resting dendritic cells.

Conclusions

Therefore, hsa-miR-342-3p-/hsa-miR-360-RASAL2, as well as hsa-miR-360-ENTPD1/PLXDC2, may serve as diagnostic and therapeutic targets for T2DM-associated PD.

LC-MS metabolomics analysis of serum metabolites during neoadjuvant chemoradiotherapy in locally advanced rectal cancer

BACKGROUND

This study aimed to investigate the serum metabolite profiles during neoadjuvant chemoradiotherapy (NCRT) in locally advanced rectal cancer (LARC) using liquid chromatography-mass spectrometry (LC-MS) metabolomics analysis.

METHODS

60 serum samples were collected from 20 patients with LARC before, during, and after radiotherapy. LC-MS metabolomics analysis was performed to identify the metabolite variations. Functional annotation was applied to discover altered metabolic pathways. The key metabolites were screened and their ability to predict sensitivity to radiotherapy was calculated using random forests and ROC curves.

RESULTS

The results showed that NCRT led to significant changes in the serum metabolite profiles. The serum metabolic profiles showed an apparent separation between different time points and different sensitivity groups. Moreover, the functional annotation showed that the differential metabolites were associated with a series of important metabolic pathways. Pre-radiotherapy (3Z,6Z)-3,6-Nonadiena and pro-radiotherapy 1-Hydroxyibuprofen showed good predictive performance in discriminating the sensitive and non-sensitive group to NCRT, with an AUC of 0.812 and 0.75, respectively. Importantly, the combination of different metabolites significantly increased the predictive ability.

CONCLUSION

This study demonstrated the potential of LC-MS metabolomics for revealing the serum metabolite profiles during NCRT in LARC. The identified metabolites may serve as potential biomarkers and therapeutic targets for the management of this disease. Furthermore, the understanding of the affected metabolic pathways may help design more personalized therapeutic strategies for LARC patients.

Predictive modeling of plastic pyrolysis process for the evaluation of activation energy: Explainable artificial intelligence based comprehensive insights

Pyrolysis, a thermochemical conversion approach of transforming plastic waste to energy has tremendous potential to manage the exponentially increasing plastic waste. However, understanding the process kinetics is fundamental to engineering a sustainable process. Conventional analysis techniques do not provide insights into the influence of characteristics of feedstock on the process kinetics. Present study exemplifies the efficacy of using machine learning for predictive modeling of pyrolysis of waste plastics to understand the complexities of the interrelations of predictor variables and their influence on activation energy. The activation energy for pyrolysis of waste plastics was evaluated using machine learning models namely Random Forest, XGBoost, CatBoost, and AdaBoost regression models. Feature selection based on the multicollinearity of data and hyperparameter tuning of the models utilizing RandomizedSearchCV was conducted. Random forest model outperformed the other models with coefficient of determination (R2) value of 0.941, root mean square error (RMSE) value of 14.69 and mean absolute error (MAE) value of 8.66 for the testing dataset. The explainable artificial intelligence-based feature importance plot and the summary plot of the shapely additive explanations projected fixed carbon content, ash content, conversion value, and carbon content as significant parameters of the model in the order; fixed carbon > carbon > ash content > degree of conversion. Present study highlighted the potential of machine learning as a powerful tool to understand the influence of the characteristics of plastic waste and the degree of conversion on the activation energy of a process that is essential for designing the large-scale operations and future scale-up of the process.

Mining bone metastasis related key genes of prostate cancer from the STING pathway based on machine learning

Background

Prostate cancer (PCa) is the second most prevalent malignant tumor in male, and bone metastasis occurs in about 70% of patients with advanced disease. The STING pathway, an innate immune signaling mechanism, has been shown to play a key role in tumorigenesis, metastasis, and cancerous bone pain. Hence, exploring regulatory mechanism of STING in PCa bone metastasis will bring novel opportunities for treating PCa bone metastasis.

Methods

First, key genes were screened from STING-related genes (SRGs) based on random forest algorithm and their predictive performance was evaluated. Subsequently, a comprehensive analysis of key genes was performed to explore their roles in prostate carcinogenesis, metastasis and tumor immunity. Next, cellular experiments were performed to verify the role of RELA in proliferation and migration in PCa cells, meanwhile, based on immunohistochemistry, we verified the difference of RELA expression between PCa primary foci and bone metastasis. Finally, based on the key genes to construct an accurate and reliable nomogram, and mined targeting drugs of key genes.

Results

In this study, three key genes for bone metastasis were mined from SRGs based on the random forest algorithm. Evaluation analysis showed that the key genes had excellent prediction performance, and it also showed that the key genes played a key role in carcinogenesis, metastasis and tumor immunity in PCa by comprehensive analysis. In addition, cellular experiments and immunohistochemistry confirmed that overexpression of RELA significantly inhibited the proliferation and migration of PCa cells, and RELA was significantly low-expression in bone metastasis. Finally, the constructed nomogram showed excellent predictive performance in Receiver Operating Characteristic (ROC, AUC = 0.99) curve, calibration curve, and Decision Curve Analysis (DCA) curve; and the targeted drugs showed good molecular docking effects.

Conclusion

In sum, this study not only provides a new theoretical basis for the mechanism of PCa bone metastasis, but also provides novel therapeutic targets and novel diagnostic tools for advanced PCa treatment.

A comparative study of feature selection and feature extraction methods for financial distress identification

Financial distress identification remains an essential topic in the scientific literature due to its importance for society and the economy. The advancements in information technology and the escalating volume of stored data have led to the emergence of financial distress that transcends the realm of financial statements and its' indicators (ratios). The feature space could be expanded by incorporating new perspectives on feature data categories such as macroeconomics, sectors, social, board, management, judicial incident, etc. However, the increased dimensionality results in sparse data and overfitted models. This study proposes a new approach for efficient financial distress classification assessment by combining dimensionality reduction and machine learning techniques. The proposed framework aims to identify a subset of features leading to the minimization of the loss function describing the financial distress in an enterprise. During the study, 15 dimensionality reduction techniques with different numbers of features and 17 machine-learning models were compared. Overall, 1,432 experiments were performed using Lithuanian enterprise data covering the period from 2015 to 2022. Results revealed that the artificial neural network (ANN) model with 30 ranked features identified using the Random Forest mean decreasing Gini (RF_MDG) feature selection technique provided the highest AUC score. Moreover, this study has introduced a novel approach for feature extraction, which could improve financial distress classification models.

Predictive Model to Identify the Long Time Survivor in Patients with Glioblastoma: A Cohort Study Integrating Machine Learning Algorithms

We aimed to develop and validate a predictive model for identifying long-term survivors (LTS) among glioblastoma (GB) patients, defined as those with an overall survival (OS) of more than 3 years. A total of 293 GB patients from CGGA and 169 from TCGA database were assigned to training and validation cohort, respectively. The differences in expression of immune checkpoint genes (ICGs) and immune infiltration landscape were compared between LTS and short time survivor (STS) (OS<1.5 years). The differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) were used to identify the genes differentially expressed between LTS and STS. Three different machine learning algorithms were employed to select the predictive genes from the overlapping region of DEGs and WGCNA to construct the nomogram. The comparison between LTS and STS revealed that STS exhibited an immune-resistant status, with higher expression of ICGs (P<0.05) and greater infiltration of immune suppression cells compared to LTS (P<0.05). Four genes, namely, OSMR, FMOD, CXCL14, and TIMP1, were identified and incorporated into the nomogram, which possessed good potential in predicting LTS probability among GB patients both in the training (C-index, 0.791; 0.772-0.817) and validation cohort (C-index, 0.770; 0.751-0.806). STS was found to be more likely to exhibit an immune-cold phenotype. The identified predictive genes were used to construct the nomogram with potential to identify LTS among GB patients.

Potential therapeutic targets for COVID-19 complicated with pulmonary hypertension: a bioinformatics and early validation study

Coronavirus disease (COVID-19) and pulmonary hypertension (PH) are closely correlated. However, the mechanism is still poorly understood. In this article, we analyzed the molecular action network driving the emergence of this event. Two datasets (GSE113439 and GSE147507) from the GEO database were used for the identification of differentially expressed genes (DEGs).Common DEGs were selected by VennDiagram and their enrichment in biological pathways was analyzed. Candidate gene biomarkers were selected using three different machine-learning algorithms (SVM-RFE, LASSO, RF).The diagnostic efficacy of these foundational genes was validated using independent datasets. Eventually, we validated molecular docking and medication prediction. We found 62 common DEGs, including several ones that could be enriched for Immune Response and Inflammation. Two DEGs (SELE and CCL20) could be identified by machine-learning algorithms. They performed well in diagnostic tests on independent datasets. In particular, we observed an upregulation of functions associated with the adaptive immune response, the leukocyte-lymphocyte-driven immunological response, and the proinflammatory response. Moreover, by ssGSEA, natural killer T cells, activated dendritic cells, activated CD4 T cells, neutrophils, and plasmacytoid dendritic cells were correlated with COVID-19 and PH, with SELE and CCL20 showing the strongest correlation with dendritic cells. Potential therapeutic compounds like FENRETI-NIDE, AFLATOXIN B1 and 1-nitropyrene were predicted. Further molecular docking and molecular dynamics simulations showed that 1-nitropyrene had the most stable binding with SELE and CCL20.The findings indicated that SELE and CCL20 were identified as novel diagnostic biomarkers for COVID-19 complicated with PH, and the target of these two key genes, FENRETI-NIDE and 1-nitropyrene, was predicted to be a potential therapeutic target, thus providing new insights into the prediction and treatment of COVID-19 complicated with PH in clinical practice.

Contextualizing injury severity from occupational accident reports using an optimized deep learning prediction model

Background

This study introduced a novel approach for predicting occupational injury severity by leveraging deep learning-based text classification techniques to analyze unstructured narratives. Unlike conventional methods that rely on structured data, our approach recognizes the richness of information within injury narrative descriptions with the aim of extracting valuable insights for improved occupational injury severity assessment.

Methods

Natural language processing (NLP) techniques were harnessed to preprocess the occupational injury narratives obtained from the US Occupational Safety and Health Administration (OSHA) from January 2015 to June 2023. The methodology involved meticulous preprocessing of textual narratives to standardize text and eliminate noise, followed by the innovative integration of Term Frequency-Inverse Document Frequency (TF-IDF) and Global Vector (GloVe) word embeddings for effective text representation. The proposed predictive model adopts a novel Bidirectional Long Short-Term Memory (Bi-LSTM) architecture and is further refined through model optimization, including random search hyperparameters and in-depth feature importance analysis. The optimized Bi-LSTM model has been compared and validated against other machine learning classifiers which are naïve Bayes, support vector machine, random forest, decision trees, and K-nearest neighbor.

Results

The proposed optimized Bi-LSTM models' superior predictability, boasted an accuracy of 0.95 for hospitalization and 0.98 for amputation cases with faster model processing times. Interestingly, the feature importance analysis revealed predictive keywords related to the causal factors of occupational injuries thereby providing valuable insights to enhance model interpretability.

Conclusion

Our proposed optimized Bi-LSTM model offers safety and health practitioners an effective tool to empower workplace safety proactive measures, thereby contributing to business productivity and sustainability. This study lays the foundation for further exploration of predictive analytics in the occupational safety and health domain.

CFLAR: A novel diagnostic and prognostic biomarker in soft tissue sarcoma, which positively modulates the immune response in the tumor microenvironment

Anoikis is highly associated with tumor cell apoptosis and tumor prognosis; however, the specific role of anoikis-related genes (ARGs) in soft tissue sarcoma (STS) remains to be fully elucidated. The present study aimed to use a variety of bioinformatics methods to determine differentially expressed anoikis-related genes in STS and healthy tissues. Subsequently, three machine learning algorithms, Least Absolute Shrinkage and Selection Operator, Support Vector Machine and Random Forest, were used to screen genes with the highest importance score. The results of the bioinformatics analyses demonstrated that CASP8 and FADD-like apoptosis regulator (CFLAR) exhibited the highest importance score. Subsequently, the diagnostic and prognostic value of CFLAR in STS development was determined using multiple public and in-house cohorts. The results of the present study demonstrated that CFLAR may be considered a diagnostic and prognostic marker of STS, which acts as an independent prognostic factor of STS development. The present study also aimed to explore the potential role of CFLAR in the STS tumor microenvironment, and the results demonstrated that CFLAR significantly enhanced the immune response of STS, and exerted a positive effect on the infiltration of CD8+ T cells and M1 macrophages in the STS immune microenvironment. Notably, the aforementioned results were verified using multiplex immunofluorescence analysis. Collectively, the results of the present study demonstrated that CFLAR may act as a novel diagnostic and prognostic marker for STS, and may positively regulate the immune response of STS. Thus, the present study provided a novel theoretical basis for the use of CFLAR in STS diagnosis, in predicting clinical outcomes and in tailoring individualized treatment options.

Phthalate metabolites and sex steroid hormones in relation to obesity in US adults: NHANES 2013-2016

Background

Obesity and metabolic syndrome pose significant health challenges in the United States (US), with connections to disruptions in sex hormone regulation. The increasing prevalence of obesity and metabolic syndrome might be associated with exposure to phthalates (PAEs). Further exploration of the impact of PAEs on obesity is crucial, particularly from a sex hormone perspective.

Methods

A total of 7780 adult participants in the National Health and Nutrition Examination Survey (NHANES) from 2013 to 2016 were included in the study. Principal component analysis (PCA) coupled with multinomial logistic regression was employed to elucidate the association between urinary PAEs metabolite concentrations and the likelihood of obesity. Weighted quartiles sum (WQS) regression was utilized to consolidate the impact of mixed PAEs exposure on sex hormone levels (total testosterone (TT), estradiol and sex hormone-binding globulin (SHBG)). We also delved into machine learning models to accurately discern obesity status and identify the key variables contributing most to these models.

Results

Principal Component 1 (PC1), characterized by mono(2-ethyl-5-carboxypentyl) phthalate (MECPP), mono(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), and mono(2-ethyl-5-oxohexyl) phthalate (MEOHP) as major contributors, exhibited a negative association with obesity. Conversely, PC2, with monocarboxyononyl phthalate (MCNP), monocarboxyoctyl phthalate (MCOP), and mono(3-carboxypropyl) phthalate (MCPP) as major contributors, showed a positive association with obesity. Mixed exposure to PAEs was associated with decreased TT levels and increased estradiol and SHBG. During the exploration of the interrelations among obesity, sex hormones, and PAEs, models based on Random Forest (RF) and eXtreme Gradient Boosting (XGBoost) algorithms demonstrated the best classification efficacy. In both models, sex hormones exhibited the highest variable importance, and certain phthalate metabolites made significant contributions to the model's performance.

Conclusions

Individuals with obesity exhibit lower levels of TT and SHBG, accompanied by elevated estradiol levels. Exposure to PAEs disrupts sex hormone levels, contributing to an increased risk of obesity in US adults. In the exploration of the interrelationships among these three factors, the RF and XGBoost algorithm models demonstrated superior performance, with sex hormones displaying higher variable importance.

Accumulation mechanism of metabolite markers identified by machine learning between Qingyuan and Xiushui counties in Polygonatum cyrtonema Hua

Polygonatum cyrtonema Hua is a traditional Chinese medicinal plant acclaimed for its therapeutic potential in diabetes and various chronic diseases. Its rhizomes are the main functional parts rich in secondary metabolites, such as flavonoids and saponins. But their quality varies by region, posing challenges for industrial and medicinal application of P. cyrtonema. In this study, 482 metabolites were identified in P. cyrtonema rhizome from Qingyuan and Xiushui counties. Cluster analysis showed that samples between these two regions had distinct secondary metabolite profiles. Machine learning methods, specifically support vector machine-recursive feature elimination and random forest, were utilized to further identify metabolite markers including flavonoids, phenolic acids, and lignans. Comparative transcriptomics and weighted gene co-expression analysis were performed to uncover potential candidate genes including CHI, UGT1, and PcOMT10/11/12/13 associated with these compounds. Functional assays using tobacco transient expression system revealed that PcOMT10/11/12/13 indeed impacted metabolic fluxes of the phenylpropanoid pathway and phenylpropanoid-related metabolites such as chrysoeriol-6,8-di-C-glucoside, syringaresinol-4'-O-glucopyranosid, and 1-O-Sinapoyl-D-glucose. These findings identified metabolite markers between these two regions and provided valuable genetic insights for engineering the biosynthesis of these compounds.

Identification of CXC Chemokine Receptor 2 (CXCR2) as a Novel Eosinophils-Independent Diagnostic Biomarker of Pediatric Eosinophilic Esophagitis by Integrated Bioinformatic and Machine-Learning Analysis

Background

Eosinophilic esophagitis (EoE) is a complex allergic condition frequently accompanied by various atopic comorbidities in children, which significantly affects their life qualities. Therefore, this study aimed to evaluate pivotal molecular markers that may facilitate the diagnosis of EoE in pediatric patients.

Methods

Three available EoE-associated gene expression datasets in children: GSE184182, GSE 197702, GSE55794, along with GSE173895 were downloaded from the GEO database. Differentially expressed genes (DEGs) identified by "limma" were intersected with key module genes identified by weighted gene co-expression network analysis (WGCNA), and the shared genes went through functional enrichment analysis. The protein-protein interaction (PPI) network and the machine learning algorithms: least absolute shrinkage and selection operator (LASSO), random forest (RF), and XGBoost were used to reveal candidate diagnostic markers for EoE. The receiver operating characteristic (ROC) curve showed the efficacy of differential diagnosis of this marker, along with online databases predicting its molecular regulatory network. Finally, we performed gene set enrichment analysis (GSEA) and assessed immune cell infiltration of EoE/control samples by using the CIBERSORT algorithm. The correlations between the key diagnostic biomarker and immune cells were also investigated.

Results

The intersection of 936 DEGs and 1446 key module genes in EoE generated 567 genes, which were primarily enriched in immune regulation. Following the construction of the PPI network and filtration by machine learning, CXCR2 served as a potential diagnostic biomarker of pediatric EoE with a perfect diagnostic efficacy (AUC = ~1.00) in regional tissue/peripheral whole blood samples. Multiple infiltrated immune cells were observed to participate in disrupting the homeostasis of esophageal epithelium to varying degrees.

Conclusion

The immune-correlated CXCR2 gene was proved to be a promising diagnostic indicator for EoE, and dysregulated regulatory T cells (Tregs)/neutrophils might play a crucial role in the pathogenesis of EoE in children.

PINK1 dominated mitochondria associated genes signature predicts abdominal aortic aneurysm with metabolic syndrome

Abdominal aortic aneurysm (AAA) is typically asymptomatic but a devastating cardiovascular disorder, with overall mortality exceeding 80 % once it ruptures. Some patients with AAA may also have comorbid metabolic syndrome (MS), suggesting a potential common underlying pathogenesis. Mitochondrial dysfunction has been reported as a key factor contributing to the deterioration of both AAA and MS. However, the intricate interplay between metabolism and mitochondrial function, both contributing to the development of AAA, has not been thoroughly explored. In this study, we identified candidate genes related to mitochondrial function in AAA and MS. Subsequently, we developed a nomoscore model comprising hub genes (PINK1, ACSL1, CYP27A1, and SLC25A11), identified through the application of two machine learning algorithms, to predict AAA. We observed a marked disparity in immune infiltration profiles between high- and low-nomoscore groups. Furthermore, we confirmed a significant upregulation of the expression of the four hub genes in AAA tissues. Among these, ACSL1 showed relatively higher expression in LPS-treated RAW264.7 cell lines, while CYP27A1 exhibited a notable decrease. Moreover, SLC25A11 displayed a significant upregulation in AngII-treated VSMCs. Conversely, the expression level of PINK1 declined in LPS-stimulated RAW264.7 cell lines but significantly increased in AngII-treated VSMCs. In vivo experiments revealed that the activation of PINK1-mediated mitophagy inhibited the development of AAA in mice. In this current study, we have innovatively identified four mitochondrial function-related genes through integrated bioinformatic analysis. This discovery sheds light on the regulatory mechanisms and unveils promising therapeutic targets for the comorbidity of AAA and MS.

Integrative analysis of potential diagnostic markers and therapeutic targets for glomerulus-associated diabetic nephropathy based on cellular senescence

Introduction

Diabetic nephropathy (DN), distinguished by detrimental changes in the renal glomeruli, is regarded as the leading cause of death from end-stage renal disease among diabetics. Cellular senescence plays a paramount role, profoundly affecting the onset and progression of chronic kidney disease (CKD) and acute kidney injuries. This study was designed to delve deeply into the pathological mechanisms between glomerulus-associated DN and cellular senescence.

Methods

Glomerulus-associated DN datasets and cellular senescence-related genes were acquired from the Gene Expression Omnibus (GEO) and CellAge database respectively. By integrating bioinformatics and machine learning methodologies including the LASSO regression analysis and Random Forest, we screened out four signature genes. The receiver operating characteristic (ROC) curve was performed to evaluate the diagnostic performance of the selected genes. Rigorous experimental validations were subsequently conducted in the mouse model to corroborate the identification of three signature genes, namely LOX, FOXD1 and GJA1. Molecular docking with chlorogenic acids (CGA) was further established not only to validate LOX, FOXD1 and GJA1 as diagnostic markers but also reveal their potential therapeutic effects.

Results and discussion

In conclusion, our findings pinpointed three diagnostic markers of glomerulus-associated DN on the basis of cellular senescence. These markers could not only predict an increased risk of DN progression but also present promising therapeutic targets, potentially ushering in innovative treatments for DN in the elderly population.

Identification and validation of key biomarkers associated with macrophages in nonalcoholic fatty liver disease based on hdWGCNA and machine learning

BACKGROUND

NAFLD has attracted increasing attention because of its high prevalence and risk of progression to cirrhosis or even hepatocellular carcinoma. Therefore, research into the root causes and molecular indicators of NAFLD is crucial.

METHODS

We analyzed scRNA-seq data and RNA-seq data from normal and NAFLD liver samples. We utilized hdWGCNA to find module-related genes associated with the phenotype. Multiple machine learning algorithms were used to validate the model diagnostics and further screen for genes that are characteristic of NAFLD. The NAFLD mouse model was constructed using the MCD diet to validate the diagnostic effect of the genes.

RESULTS

We identified a specific macrophage population called NASH-macrophages by single-cell sequencing analysis. Cell communication analysis and Pseudo-time trajectory analysis revealed the specific role and temporal distribution of NASH-macrophages in NAFLD. The hdWGCNA screening yielded 30 genes associated with NASH-macrophages, and machine learning algorithms screened and obtained two genes characterizing NAFLD. The immune infiltration indicated that these genes were highly associated with macrophages. Notably, we verified by RT-qPCR, IHC, and WB that MAFB and CX3CR1 are highly expressed in the MCD mouse model and may play important roles.

CONCLUSIONS

Our study revealed a macrophage population that is closely associated with NAFLD. Using hdWGCNA analysis and multiple machine learning algorithms, we identified two NAFLD signature genes that are highly correlated with macrophages. Our findings may provide potential feature markers and therapeutic targets for NAFLD.

Analysis and validation of diagnostic biomarkers and immune cell infiltration characteristics in pediatric sepsis by integrating bioinformatics and machine learning

BACKGROUND

Pediatric sepsis is a complicated condition characterized by life-threatening organ failure resulting from a dysregulated host response to infection in children. It is associated with high rates of morbidity and mortality, and rapid detection and administration of antimicrobials have been emphasized. The objective of this study was to evaluate the diagnostic biomarkers of pediatric sepsis and the function of immune cell infiltration in the development of this illness.

METHODS

Three gene expression datasets were available from the Gene Expression Omnibus collection. First, the differentially expressed genes (DEGs) were found with the use of the R program, and then gene set enrichment analysis was carried out. Subsequently, the DEGs were combined with the major module genes chosen using the weighted gene co-expression network. The hub genes were identified by the use of three machine-learning algorithms: random forest, support vector machine-recursive feature elimination, and least absolute shrinkage and selection operator. The receiver operating characteristic curve and nomogram model were used to verify the discrimination and efficacy of the hub genes. In addition, the inflammatory and immune status of pediatric sepsis was assessed using cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT). The relationship between the diagnostic markers and infiltrating immune cells was further studied.

RESULTS

Overall, after overlapping key module genes and DEGs, we detected 402 overlapping genes. As pediatric sepsis diagnostic indicators, CYSTM1 (AUC = 0.988), MMP8 (AUC = 0.973), and CD177 (AUC = 0.986) were investigated and demonstrated statistically significant differences (P < 0.05) and diagnostic efficacy in the validation set. As indicated by the immune cell infiltration analysis, multiple immune cells may be involved in the development of pediatric sepsis. Additionally, all diagnostic characteristics may correlate with immune cells to varying degrees.

CONCLUSIONS

The candidate hub genes (CD177, CYSTM1, and MMP8) were identified, and the nomogram was constructed for pediatric sepsis diagnosis. Our study could provide potential peripheral blood diagnostic candidate genes for pediatric sepsis patients.

Identifying potential biomarkers of idiopathic pulmonary fibrosis through machine learning analysis

Idiopathic pulmonary fibrosis (IPF) is the most common and serious type of idiopathic interstitial pneumonia, characterized by chronic, progressive, and low survival rates, while unknown disease etiology. Until recently, patients with idiopathic pulmonary fibrosis have a poor prognosis, high mortality, and limited treatment options, due to the lack of effective early diagnostic and prognostic tools. Therefore, we aimed to identify biomarkers for idiopathic pulmonary fibrosis based on multiple machine-learning approaches and to evaluate the role of immune infiltration in the disease. The gene expression profile and its corresponding clinical data of idiopathic pulmonary fibrosis patients were downloaded from Gene Expression Omnibus (GEO) database. Next, the differentially expressed genes (DEGs) with the threshold of FDR < 0.05 and |log2 foldchange (FC)| > 0.585 were analyzed via R package "DESeq2" and GO enrichment and KEGG pathways were run in R software. Then, least absolute shrinkage and selection operator (LASSO) logistic regression, support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF) algorithms were combined to screen the key potential biomarkers of idiopathic pulmonary fibrosis. The diagnostic performance of these biomarkers was evaluated through receiver operating characteristic (ROC) curves. Moreover, the CIBERSORT algorithm was employed to assess the infiltration of immune cells and the relationship between the infiltrating immune cells and the biomarkers. Finally, we sought to understand the potential pathogenic role of the biomarker (SLAIN1) in idiopathic pulmonary fibrosis using a mouse model and cellular model. A total of 3658 differentially expressed genes of idiopathic pulmonary fibrosis were identified, including 2359 upregulated genes and 1299 downregulated genes. FHL2, HPCAL1, RNF182, and SLAIN1 were identified as biomarkers of idiopathic pulmonary fibrosis using LASSO logistic regression, RF, and SVM-RFE algorithms. The ROC curves confirmed the predictive accuracy of these biomarkers both in the training set and test set. Immune cell infiltration analysis suggested that patients with idiopathic pulmonary fibrosis had a higher level of B cells memory, Plasma cells, T cells CD8, T cells follicular helper, T cells regulatory (Tregs), Macrophages M0, and Mast cells resting compared with the control group. Correlation analysis demonstrated that FHL2 was significantly associated with the infiltrating immune cells. qPCR and western blotting analysis suggested that SLAIN1 might be a signature for the diagnosis of idiopathic pulmonary fibrosis. In this study, we identified four potential biomarkers (FHL2, HPCAL1, RNF182, and SLAIN1) and evaluated the potential pathogenic role of SLAIN1 in idiopathic pulmonary fibrosis. These findings may have great significance in guiding the understanding of disease mechanisms and potential therapeutic targets in idiopathic pulmonary fibrosis.

Public health factors help explain cross country heterogeneity in excess death during the COVID19 pandemic

The COVID-19 pandemic has taken a devastating toll around the world. Since January 2020, the World Health Organization estimates 14.9 million excess deaths have occurred globally. Despite this grim number quantifying the deadly impact, the underlying factors contributing to COVID-19 deaths at the population level remain unclear. Prior studies indicate that demographic factors like proportion of population older than 65 and population health explain the cross-country difference in COVID-19 deaths. However, there has not been a comprehensive analysis including variables describing government policies and COVID-19 vaccination rate. Furthermore, prior studies focus on COVID-19 death rather than excess death to assess the impact of the pandemic. Through a robust statistical modeling framework, we analyze 80 countries and show that actionable public health efforts beyond just the factors intrinsic to each country are important for explaining the cross-country heterogeneity in excess death.

Hub gene identification and molecular subtype construction for Helicobacter pylori in gastric cancer via machine learning methods and NMF algorithm

Helicobacter pylori (HP) is a gram-negative and spiral-shaped bacterium colonizing the human stomach and has been recognized as the risk factor of gastritis, peptic ulcer disease, and gastric cancer (GC). Moreover, it was recently identified as a class I carcinogen, which affects the occurrence and progression of GC via inducing various oncogenic pathways. Therefore, identifying the HP-related key genes is crucial for understanding the oncogenic mechanisms and improving the outcomes of GC patients. We retrieved the list of HP-related gene sets from the Molecular Signatures Database. Based on the HP-related genes, unsupervised non-negative matrix factorization (NMF) clustering method was conducted to stratify TCGA-STAD, GSE15459, GSE84433 samples into two clusters with distinct clinical outcomes and immune infiltration characterization. Subsequently, two machine learning (ML) strategies, including support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF), were employed to determine twelve hub HP-related genes. Beyond that, receiver operating characteristic and Kaplan-Meier curves further confirmed the diagnostic value and prognostic significance of hub genes. Finally, expression of HP-related hub genes was tested by qRT-PCR array and immunohistochemical images. Additionally, functional pathway enrichment analysis indicated that these hub genes were implicated in the genesis and progression of GC by activating or inhibiting the classical cancer-associated pathways, such as epithelial-mesenchymal transition, cell cycle, apoptosis, RAS/MAPK, etc. In the present study, we constructed a novel HP-related tumor classification in different datasets, and screened out twelve hub genes via performing the ML algorithms, which may contribute to the molecular diagnosis and personalized therapy of GC.

Filter and Wrapper Stacking Ensemble (FWSE): a robust approach for reliable biomarker discovery in high-dimensional omics data

Selecting informative features, such as accurate biomarkers for disease diagnosis, prognosis and response to treatment, is an essential task in the field of bioinformatics. Medical data often contain thousands of features and identifying potential biomarkers is challenging due to small number of samples in the data, method dependence and non-reproducibility. This paper proposes a novel ensemble feature selection method, named Filter and Wrapper Stacking Ensemble (FWSE), to identify reproducible biomarkers from high-dimensional omics data. In FWSE, filter feature selection methods are run on numerous subsets of the data to eliminate irrelevant features, and then wrapper feature selection methods are applied to rank the top features. The method was validated on four high-dimensional medical datasets related to mental illnesses and cancer. The results indicate that the features selected by FWSE are stable and statistically more significant than the ones obtained by existing methods while also demonstrating biological relevance. Furthermore, FWSE is a generic method, applicable to various high-dimensional datasets in the fields of machine intelligence and bioinformatics.

Thrombus radiomics in patients with anterior circulation acute ischemic stroke undergoing endovascular treatment

BACKGROUND

Thrombus radiomics (TR) describe complex shape and textural thrombus imaging features. We aimed to study the relationship of TR extracted from non-contrast CT with procedural and functional outcome in endovascular-treated patients with acute ischemic stroke.

METHODS

Thrombi were segmented on thin-slice non-contrast CT (≤1 mm) from 699 patients included in the MR CLEAN Registry. In a pilot study, we selected 51 TR with consistent values across two raters' segmentations (ICC >0.75). Random forest models using TR in addition or as a substitute to baseline clinical variables (CV) and manual thrombus measurements (MTM) were trained with 499 patients and evaluated on 200 patients for predicting successful reperfusion (extended Thrombolysis in Cerebral Ischemia (eTICI) ≥2B), first attempt reperfusion, reperfusion within three attempts, and functional independence (modified Rankin Scale (mRS) ≤2). Three texture and shape features were selected based on feature importance and related to eTICI ≥2B, number of attempts to eTICI ≥2B, and 90-day mRS with ordinal logistic regression.

RESULTS

Random forest models using TR, CV or MTM had comparable predictive performance. Thrombus texture (inverse difference moment normalized) was independently associated with reperfusion (adjusted common OR (acOR) 0.85, 95% CI 0.72 to 0.99). Thrombus volume and texture were also independently associated with the number of attempts to successful reperfusion (acOR 1.36, 95% CI 1.03 to 1.88 and acOR 1.24, 95% CI 1.04 to 1.49).

CONCLUSIONS

TR describing thrombus volume and texture were associated with more attempts to successful reperfusion. Compared with models using CV and MTM, TR had no added value for predicting procedural and functional outcome.

Recurrent Implantation Failure: Bioinformatic Discovery of Biomarkers and Identification of Metabolic Subtypes

Recurrent implantation failure (RIF) is a challenging scenario from different standpoints. This study aimed to investigate its correlation with the endometrial metabolic characteristics. Transcriptomics data of 70 RIF and 99 normal endometrium tissues were retrieved from the Gene Expression Omnibus database. Common differentially expressed metabolism-related genes were extracted and various enrichment analyses were applied. Then, RIF was classified using a consensus clustering approach. Three machine learning methods were employed for screening key genes, and they were validated through the RT-qPCR experiment in the endometrium of 10 RIF and 10 healthy individuals. Receiver operator characteristic (ROC) curves were generated and validated by 20 RIF and 20 healthy individuals from Peking University People's Hospital. We uncovered 109 RIF-related metabolic genes and proposed a novel two-subtype RIF classification according to their metabolic features. Eight characteristic genes (SRD5A1, POLR3E, PPA2, PAPSS1, PRUNE, CA12, PDE6D, and RBKS) were identified, and the area under curve (AUC) was 0.902 and the external validated AUC was 0.867. Higher immune cell infiltration levels were found in RIF patients and a metabolism-related regulatory network was constructed. Our work has explored the metabolic and immune characteristics of RIF, which paves a new road to future investigation of the related pathogenic mechanisms.

Machine learning uncovers accumulation mechanism of flavonoid compounds in Polygonatum cyrtonema Hua

The compositions and yield of flavonoid compounds of Polygonatum cyrtonema Hua (PC) are important indices of the quality of medicinal materials. However, the flavonoids compositions and accumulation mechanism are still unclear in PC. Here, we identified 22 flavonoids using widely-targeted metabolome analysis in 15 genotypes of PC. Then weighted gene co-expression network analysis based on 45 transcriptome samples was performed to construct 12 co-expressed modules, in which blue module highly correlated with flavonoids was identified. Furthermore, 4 feature genes including PcCHS1, PcCHI, PcCHS2 and PcCHR5 were identified from 94 hub genes in blue module via machine learning methods support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF), and their functions on metabolic flux of flavonoids pathway were confirmed by tobacco transient expression system. Our findings identified representative flavonoids and key enzymes in PC that provided new insight for elite breeding rich in flavonoids, and thus will be beneficial for rapid development of great potential economic and medicinal value of PC.

The activity of cuproptosis pathway calculated by AUCell algorithm was employed to construct cuproptosis landscape in lung adenocarcinoma

Cuproptosis is a recently described copper-dependent cell death pathway. Consequently, there are still few studies on lung adenocarcinoma (LUAD)-related cuproptosis, and we aimed to deepen in this matter. In this study, data from 503 patients with lung cancer from the TCGA-LUAD cohort data collection and 11 LUAD single-cells from GSE131907 as well as from 10 genes associated with cuproptosis were analyzed. The AUCell R package was used to determine the copper-dependent cell death pathway activity for each cell subpopulation, calculate the CellChat score, and display cell communication for each cell subpopulation. The PROGENy score was calculated to show the scores of tumor-related pathways in different cell populations. GO and KEGG analyses were used to calculate pathway activity. Univariate COX and random forest analyses were used to screen prognosis-associated genes and construct models. The ssGSEA and xCell algorithms were used to calculate the immunocyte infiltration score. Based on data from the GDSC database, the drug sensitivity score was calculated using oncoPredict. Finally, in vitro experiments were performed to determine the role of TLE1, the most important gene in the prognostic model. The 11 LUAD single-cell samples were classified into 8 different cell populations, from which epithelial cells showed the highest copper-dependent cell death pathway activity. Epithelial cell subsets were significantly positively correlated with MAKP, hypoxia, and other pathways. In addition, cell subgroup communication showed highly active collagen and APP pathways. Using the Findmark algorithm, differentially expressed genes (DEGs) between epithelial and other cell types were identified. Combined with the bulk data in the TCGA-LUAD database, DEGs were enriched in pathways such as EGFR tyrosine kinase inhibitor resistance, Hippo signaling pathway, and tight junction. Subsequently, we selected 4 genes (out of 112) with prognostic significance, ANKRD29, RHOV, TLE1, and NPAS2, and used them to construct a prognostic model. The high- and low-risk groups, distinguished by the median risk score, showed significantly different prognoses. Finally, we chose TLE1 as a biomarker based on the relative importance score in the prognostic model. In vitro experiments showed that TLE1 promotes tumor proliferation and migration and inhibits apoptosis.

Identifying hub genes and common biological pathways between COVID-19 and benign prostatic hyperplasia by machine learning algorithms

Background

COVID-19, a serious respiratory disease that has the potential to affect numerous organs, is a serious threat to the health of people around the world. The objective of this article is to investigate the potential biological targets and mechanisms by which SARS-CoV-2 affects benign prostatic hyperplasia (BPH) and related symptoms.

Methods

We downloaded the COVID-19 datasets (GSE157103 and GSE166253) and the BPH datasets (GSE7307 and GSE132714) from the Gene Expression Omnibus (GEO) database. In GSE157103 and GSE7307, differentially expressed genes (DEGs) were found using the "Limma" package, and the intersection was utilized to obtain common DEGs. Further analyses followed, including those using Protein-Protein Interaction (PPI), Gene Ontology (GO) function enrichment analysis, and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Potential hub genes were screened using three machine learning methods, and they were later verified using GSE132714 and GSE166253. The CIBERSORT analysis and the identification of transcription factors, miRNAs, and drugs as candidates were among the subsequent analyses.

Results

We identified 97 common DEGs from GSE157103 and GSE7307. According to the GO and KEGG analyses, the primary gene enrichment pathways were immune-related pathways. Machine learning methods were used to identify five hub genes (BIRC5, DNAJC4, DTL, LILRB2, and NDC80). They had good diagnostic properties in the training sets and were validated in the validation sets. According to CIBERSORT analysis, hub genes were closely related to CD4 memory activated of T cells, T cells regulatory and NK cells activated. The top 10 drug candidates (lucanthone, phytoestrogens, etoposide, dasatinib, piroxicam, pyrvinium, rapamycin, niclosamide, genistein, and testosterone) will also be evaluated by the P value, which is expected to be helpful for the treatment of COVID-19-infected patients with BPH.

Conclusion

Our findings reveal common signaling pathways, possible biological targets, and promising small molecule drugs for BPH and COVID-19. This is crucial to understand the potential common pathogenic and susceptibility pathways between them.

Identification of necroptosis-related gene TRAF5 as potential target of diagnosing atherosclerosis and assessing its stability

BACKGROUND

Atherosclerosis (AS) is a leading cause of morbidity and mortality in older patients and features progressive formation of plaques in vascular tissues. With the progression of atherosclerosis, plaque rupture may occur and cause stroke, myocardial infarction, etc. Different forms of cell death promote the formation of a necrotic core of the plaque, leading to rupture. Necroptosis is a type of programmed cell death that contributes to the development of cardiovascular disease. However, the role of necroptosis in AS has not yet been investigated.

METHODS

The Gene Expression Omnibus (GEO) database was used to obtain gene expression profiles. Differentially expressed genes (DEGs) and necroptosis gene sets were used to identify necroptosis-related differentially expressed genes (NRDEGs). The NRDEGs were used to construct a diagnostic model and were further screened using least absolute shrinkage selection operator (LASSO) regression and random forest (RF) analysis. The discriminatory capacity of the NRDEGs was evaluated using receiver operating characteristic (ROC) curves. Immune infiltration levels were estimated based on CIBERSORTx analysis. The GSE21545 dataset, containing survival information, was used to determine prognosis-associated genes. Univariate and multivariate Cox regression analyses combined with survival analysis determined gene prognostic values. RNA and protein levels were detected by RT-qPCR and western blotting in arteriosclerosis obliterans(ASO) and normal vascular tissues. Vascular smooth muscle cells (VSMCs) were treated with oxidized low-density lipoprotein (ox-LDL) to develop cell models of advanced AS. The effects of protein knockdown on necroptosis were assessed by western blotting and flow cytometry. EdU and Cell Counting Kit-8 assays were used to examine cell proliferation.

RESULTS

TNF Receptor Associated Factor 5 (TRAF5) was identified as a diagnostic marker for AS based on the AUC value in both the GSE20129 and GSE43292 datasets. According to differential expression analysis, LASSO regression analysis, RF analysis, univariate analysis, multivariate analysis, and gene-level survival analysis, TRAF5 was markedly associated with necroptosis in AS. Silencing TRAF5 promotes necroptosis and attenuates the proliferation of ox-LDL-induced cell models of advanced AS.

CONCLUSIONS

This study identified a diagnostic marker of necroptosis-related atherosclerosis, TRAF5, which can also be used to diagnose and assess atherosclerotic plaque stability. This novel finding has important implications in the diagnosis and assessment of plaque stability in atherosclerosis.

Nonnegative matrix factorization analysis and multiple machine learning methods identified IL17C and ACOXL as novel diagnostic biomarkers for atherosclerosis

BACKGROUND

Atherosclerosis is the common pathological basis for many cardiovascular and cerebrovascular diseases. The purpose of this study is to identify the diagnostic biomarkers related to atherosclerosis through machine learning algorithm.

METHODS

Clinicopathological parameters and transcriptomics data were obtained from 4 datasets (GSE21545, GSE20129, GSE43292, GSE100927). A nonnegative matrix factorization algorithm was used to classify arteriosclerosis patients in GSE21545 dataset. Then, we identified prognosis-related differentially expressed genes (DEGs) between the subtypes. Multiple machine learning methods to detect pivotal markers. Discrimination, calibration and clinical usefulness of the predicting model were assessed using area under curve, calibration plot and decision curve analysis respectively. The expression level of the feature genes was validated in GSE20129, GSE43292, GSE100927.

RESULTS

2 molecular subtypes of atherosclerosis was identified, and 223 prognosis-related DEGs between the 2 subtypes were identified. These genes are not only related to epithelial cell proliferation, mitochondrial dysfunction, but also to immune related pathways. Least absolute shrinkage and selection operator, random forest, support vector machine- recursive feature elimination show that IL17C and ACOXL were identified as diagnostic markers of atherosclerosis. The prediction model displayed good discrimination and good calibration. Decision curve analysis showed that this model was clinically useful. Moreover, IL17C and ACOXL were verified in other 3 GEO datasets, and also have good predictive performance.

CONCLUSION

IL17C and ACOXL were diagnostic genes of atherosclerosis and associated with higher incidence of ischemic events.

Integrated bioinformatics and machine learning algorithms reveal the critical cellular senescence-associated genes and immune infiltration in heart failure due to ischemic cardiomyopathy

Heart failure (HF) is the final stage of many cardiovascular illnesses and the leading cause of death worldwide. At the same time, ischemic cardiomyopathy has replaced valvular heart disease and hypertension as the primary causes of heart failure. Cellular senescence in heart failure is currently receiving more attention. In this paper, we investigated the correlation between the immunological properties of myocardial tissue and the pathological mechanisms of cellular senescence during ischemic cardiomyopathy leading to heart failure (ICM-HF) using bioinformatics and machine learning methodologies. Our goals were to clarify the pathogenic causes of heart failure and find new treatment options. First, after obtaining GSE5406 from the Gene Expression Omnibus (GEO) database and doing limma analysis, differential genes (DEGs) among the ICM-HF and control groups were identified. We intersected these differential genes with cellular senescence-associated genes (CSAG) via the CellAge database to obtain 39 cellular senescence-associated DEGs (CSA-DEGs). Then, a functional enrichment analysis was performed to elucidate the precise biological processes by which the hub genes control cellular senescence and immunological pathways. Then, the respective key genes were identified by Random Forest (RF) method, LASSO (Least Absolute Shrinkage and Selection Operator) algorithms, and Cytoscape's MCODE plug-in. Three sets of key genes were taken to intersect to obtain three CSA-signature genes (including MYC, MAP2K1, and STAT3), and these three CSA-signature genes were validated in the test gene set (GSE57345), and Nomogram analysis was done. In addition, we assessed the relationship between these three CSA- signature genes and the immunological landscape of heart failure encompassing immunological infiltration expression profiles. This work implies that cellular senescence may have a crucial role in the pathogenesis of ICM-HF, which may be closely tied to its effect on the immune microenvironment. Exploring the molecular underpinnings of cellular senescence during ICM-HF is anticipated to yield significant advances in the disease's diagnosis and therapy.

Identification of NETs-related biomarkers and molecular clusters in systemic lupus erythematosus

Neutrophil extracellular traps (NETs) is an important process involved in the pathogenesis of systemic lupus erythematosus (SLE), but the potential mechanisms of NETs contributing to SLE at the genetic level have not been clearly investigated. This investigation aimed to explore the molecular characteristics of NETs-related genes (NRGs) in SLE based on bioinformatics analysis, and identify associated reliable biomarkers and molecular clusters. Dataset GSE45291 was acquired from the Gene Expression Omnibus repository and used as a training set for subsequent analysis. A total of 1006 differentially expressed genes (DEGs) were obtained, most of which were associated with multiple viral infections. The interaction of DEGs with NRGs revealed 8 differentially expressed NRGs (DE-NRGs). The correlation and protein-protein interaction analyses of these DE-NRGs were performed. Among them, HMGB1, ITGB2, and CREB5 were selected as hub genes by random forest, support vector machine, and least absolute shrinkage and selection operator algorithms. The significant diagnostic value for SLE was confirmed in the training set and three validation sets (GSE81622, GSE61635, and GSE122459). Additionally, three NETs-related sub-clusters were identified based on the hub genes' expression profiles analyzed by unsupervised consensus cluster assessment. Functional enrichment was performed among the three NETs subgroups, and the data revealed that cluster 1 highly expressed DEGs were prevalent in innate immune response pathways while that of cluster 3 were enriched in adaptive immune response pathways. Moreover, immune infiltration analysis also revealed that innate immune cells were markedly infiltrated in cluster 1 while the adaptive immune cells were upregulated in cluster 3. As per our knowledge, this investigation is the first to explore the molecular characteristics of NRGs in SLE, identify three potential biomarkers (HMGB1, ITGB2, and CREB5), and three distinct clusters based on these hub biomarkers.

Identification of shared fatty acid metabolism related signatures in dilated cardiomyopathy and myocardial infarction

Aim: It is to be elucidated the risk-predictive role of differentially expressed fatty acid metabolism related genes (DE-FRGs) in dilated cardiomyopathy (DCM) and myocardial infarction. Materials & methods: Four gene enrichment analyses defined DE-FRGs’ biological functions and pathways. Three strategies were applied to identify risk biomarkers and construct a nomogram. The 4-DE-FRG correlation with immune cell infiltration, drugs, and ceRNA was explored. Results: DE-FRGs were enriched in lipid metabolism. A risk nomogram was established by ACSL1, ALDH2, CYP27A1 and PPARA, demonstrating a good ability for DCM and myocardial infarction prediction. PPARA was positively correlated with adaptive immunocytes. Thirty-five drugs are candidate therapeutic targets. Conclusion: A nomogram and new biological targets for early diagnosis and treatment of DCM and myocardial infarction were provided.

Fecal Metabolites as Biomarkers for Predicting Food Intake by Healthy Adults

BACKGROUND

The fecal metabolome is affected by diet and includes metabolites generated by human and microbial metabolism. Advances in -omics technologies and analytic approaches have allowed researchers to identify metabolites and better utilize large data sets to generate usable information. One promising aspect of these advancements is the ability to determine objective biomarkers of food intake.

OBJECTIVES

We aimed to utilize a multivariate, machine learning approach to identify metabolite biomarkers that accurately predict food intake.

METHODS

Data were aggregated from 5 controlled feeding studies in adults that tested the impact of specific foods (almonds, avocados, broccoli, walnuts, barley, and oats) on the gastrointestinal microbiota. Fecal samples underwent GC-MS metabolomic analysis; 344 metabolites were detected in preintervention samples, whereas 307 metabolites were detected postintervention. After removing metabolites that were only detected in either pre- or postintervention and those undetectable in ≥80% of samples in all study groups, changes in 96 metabolites relative concentrations (treatment postintervention minus preintervention) were utilized in random forest models to 1) examine the relation between food consumption and fecal metabolome changes and 2) rank the fecal metabolites by their predictive power (i.e., feature importance score).

RESULTS

Using the change in relative concentration of 96 fecal metabolites, 6 single-food random forest models for almond, avocado, broccoli, walnuts, whole-grain barley, and whole-grain oats revealed prediction accuracies between 47% and 89%. When comparing foods with one another, almond intake was differentiated from walnut intake with 91% classification accuracy.

CONCLUSIONS

Our findings reveal promise in utilizing fecal metabolites as objective complements to certain self-reported food intake estimates. Future research on other foods at different doses and dietary patterns is needed to identify biomarkers that can be applied in feeding study compliance and clinical settings.

Network analyses of upper and lower airway transcriptomes identify shared mechanisms among children with recurrent wheezing and school-age asthma

Background

Predicting which preschool children with recurrent wheezing (RW) will develop school-age asthma (SA) is difficult, highlighting the critical need to clarify the pathogenesis of RW and the mechanistic relationship between RW and SA. Despite shared environmental exposures and genetic determinants, RW and SA are usually studied in isolation. Based on network analysis of nasal and tracheal transcriptomes, we aimed to identify convergent transcriptomic mechanisms in RW and SA.

Methods

RNA-sequencing data from nasal and tracheal brushing samples were acquired from the Gene Expression Omnibus. Combined with single-cell transcriptome data, cell deconvolution was used to infer the composition of 18 cellular components within the airway. Consensus weighted gene co-expression network analysis was performed to identify consensus modules closely related to both RW and SA. Shared pathways underlying consensus modules between RW and SA were explored by enrichment analysis. Hub genes between RW and SA were identified using machine learning strategies and validated using external datasets and quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Finally, the potential value of hub genes in defining RW subsets was determined using nasal and tracheal transcriptome data.

Results

Co-expression network analysis revealed similarities in the transcriptional networks of RW and SA in the upper and lower airways. Cell deconvolution analysis revealed an increase in mast cell fraction but decrease in club cell fraction in both RW and SA airways compared to controls. Consensus network analysis identified two consensus modules highly associated with both RW and SA. Enrichment analysis of the two consensus modules indicated that fatty acid metabolism-related pathways were shared key signals between RW and SA. Furthermore, machine learning strategies identified five hub genes, i.e., CST1, CST2, CST4, POSTN, and NRTK2, with the up-regulated hub genes in RW and SA validated using three independent external datasets and qRT-PCR. The gene signatures of the five hub genes could potentially be used to determine type 2 (T2)-high and T2-low subsets in preschoolers with RW.

Conclusions

These findings improve our understanding of the molecular pathogenesis of RW and provide a rationale for future exploration of the mechanistic relationship between RW and SA.

Identification of crucial genes for predicting the risk of atherosclerosis with system lupus erythematosus based on comprehensive bioinformatics analysis and machine learning

BACKGROUND

Systemic lupus erythematosus (SLE) has become a major public health problem over the years, and atherosclerosis (AS) is one of the main complications of SLE associated with serious cardiovascular consequences in this patient population. The present study aimed to identify potential biomarkers for SLE patients with AS.

METHODS

Five microarray datasets (GSE50772, GSE81622, GSE100927, GSE28829, GSE37356) were downloaded from the NCBI Gene Expression Omnibus database. The Limma package was used to identify differentially expressed genes (DEGs) in AS. Weighted gene coexpression network analysis (WGCNA) was used to identify significant module genes associated with SLE. Functional enrichment analysis, protein-protein interaction (PPI) network construction, and machine learning algorithms (least absolute shrinkage and selection operator (Lasso, Support Vector Machine-Recursive Feature Elimination (SVM-RFE), and random forest) were applied to identify hub genes. Subsequently, we generated a nomogram and receiver operating characteristic curve (ROC) for predicting the risk of AS in SLE patients. Finally, immune cell infiltrations were analyzed, and Consensus Cluster Analysis was conducted based on Single Sample Gene Set Enrichment Analysis (ssGSEA) scores.

RESULTS

Five hub genes (SPI1, MMP9, C1QA, CX3CR1, and MNDA) were identified and used to establish a nomogram that yielded a high predictive performance (area under the curve 0.900-0.981). Dysregulated immune cell infiltrations were found in AS, with positive correlations with the five hub genes. Consensus clustering showed that the optimal number of subtypes was 3. Compared to subtypes A and B, subtype C presented higher expression of the five hub genes, immune cell infiltration levels and immune checkpoint expression.

CONCLUSION

Our study systematically identified five candidate hub genes (SPI1, MMP9, C1QA, CX3CR1, MNDA) and established a nomogram that could predict the risk of AS with SLE using various bioinformatic analyses and machine learning algorithms. Our findings provide the foothold for future studies on potential crucial genes for AS in SLE patients. Additionally, the dysregulated immune cell proportions and immune checkpoint expressions in AS with SLE were identified.

Examining the Human Activity-Intensity Change at Different Stages of the COVID-19 Pandemic across Chinese Working, Residential and Entertainment Areas

The COVID-19 pandemic has already resulted in more than 6 million deaths worldwide as of December 2022. The COVID-19 has also been greatly affecting the activity of the human population in China and the world. It remains unclear how the human activity-intensity changes have been affected by the COVID-19 spread in China at its different stages along with the lockdown and relaxation policies. We used four days of Location-based services data from Tencent across China to capture the real-time changes in human activity intensity in three stages of COVID-19-namely, during the lockdown, at the first stage of work resuming and at the stage of total work resuming-and observed the changes in different land use categories. We applied the mean decrease Gini (MDG) approach in random forest to examine how these changes are influenced by land attributes, relying on the CART algorithm in Python. This approach was also compared with Geographically Weighted Regression (GWR). Our analysis revealed that the human activity intensity decreased by 22-35%, 9-16% and 6-15%, respectively, in relation to the normal conditions before the spread of COVID-19 during the three periods. The human activity intensity associated with commercial sites, sports facilities/gyms and tourism experienced the relatively largest contraction during the lockdown. During the relaxations of restrictions, government institutions showed a 13.89% rise in intensity at the first stage of work resuming, which was the highest rate among all the working sectors. Furthermore, the GDP and road junction density were more influenced by the change in human activity intensity for all land use categories. The bus stop density was importantly associated with mixed-use land recovery during the relaxing stages, while the coefficient of density of population in entertainment land were relatively higher at these two stages. This study aims to provide additional support to investigate the human activity changes due to the spread of COVID-19 at different stages across different sectors.

Exploring the Mechanism and Suggesting Combination Therapies for HDAC Inhibitors in Androgen Receptor-Null Prostate Cancer Using Multivariate Statistical Analysis and Data Mining Techniques

Previously, we showed that novel histone deacetylase (HDAC) inhibitors, N1-hydroxy-N ⁸-(4-(pyridine-2-carbothioamido)phenyl)octanediamide (Jazz90) and [chlorido(η⁵-pentamethylcyclopentadienyl)(N1-hydroxy-N8-(4-(pyridine-2-carbothioamido-κ² N, S)phenyl)octanediamide)rhodium(III)] chloride (Jazz167), have cytostatic and anti-angiogenic effects in androgen receptor-negative prostate cancer cells and are also non-toxic in BALB/c mice. However, only univariate statistical analysis was carried out to determine the role of individual proteins. In this study, multivariate statistical analyses (MVAs) and data mining procedures were carried out with the objective of determining the molecular networks that explain the growth inhibitory potential of Jazz90 and Jazz167 in PC3 cells and to determine potential inhibitors that can be used in combination with these HDAC inhibitors. Lasso regression revealed that angiogenic factors, vascular endothelial growth factor-A (VEGF-A), and vascular endothelial growth factor receptor-2 (VEGFR-2), alongside HDAC inhibition, predicted the reduction in cell number with an adjusted R ² value of 0.99 following Jazz90 treatment, whereas VEGFR-2, acetylation of histone-3, and HDAC inhibition predicted cell number with an adjusted R ² value of 0.84 following Jazz167 treatment. These results were further followed up with ridge regression, hierarchical cluster analysis, random forest classification (RFC), and support vector machines. RFC and support vector machines also predicted the treatment groups with a 100% accuracy. MVAs also revealed that Jazz90 should be examined in combination with epithelial to mesenchymal transitioning inhibitors, such as simvastatin and olaparib, whereas Jazz167 should be examined with venetoclax or navitoclax. Future studies should also address the roles of VEGF-A and VEGFR-2 in cellular proliferation, whereas p27 function should be examined for its role in PC3 cell migration.

Development and analysis of a comprehensive diagnostic model for aortic valve calcification using machine learning methods and artificial neural networks

BACKGROUND

Although advanced surgical and interventional treatments are available for advanced aortic valve calcification (AVC) with severe clinical symptoms, early diagnosis, and intervention is critical in order to reduce calcification progression and improve patient prognosis. The aim of this study was to develop therapeutic targets for improving outcomes for patients with AVC.

MATERIALS AND METHODS

We used the public expression profiles of individuals with AVC (GSE12644 and GSE51472) to identify potential diagnostic markers. First, the R software was used to identify differentially expressed genes (DEGs) and perform functional enrichment analysis. Next, we combined bioinformatics techniques with machine learning methodologies such as random forest algorithms and support vector machines to screen for and identify diagnostic markers of AVC. Subsequently, artificial neural networks were employed to filter and model the diagnostic characteristics for AVC incidence. The diagnostic values were determined using the receiver operating characteristic (ROC) curves. Furthermore, CIBERSORT immune infiltration analysis was used to determine the expression of different immune cells in the AVC. Finally, the CMap database was used to predict candidate small compounds as prospective AVC therapeutics.

RESULTS

A total of 78 strong DEGs were identified. The leukocyte migration and pid integrin 1 pathways were highly enriched for AVC-specific DEGs. CXCL16, GPM6A, BEX2, S100A9, and SCARA5 genes were all regarded diagnostic markers for AVC. The model was effectively constructed using a molecular diagnostic score system with significant diagnostic value (AUC = 0.987) and verified using the independent dataset GSE83453 (AUC = 0.986). Immune cell infiltration research revealed that B cell naive, B cell memory, plasma cells, NK cell activated, monocytes, and macrophage M0 may be involved in the development of AVC. Additionally, all diagnostic characteristics may have varying degrees of correlation with immune cells. The most promising small molecule medicines for reversing AVC gene expression are Doxazosin and Terfenadine.

CONCLUSION

It was identified that CXCL16, GPM6A, BEX2, S100A9, and SCARA5 are potentially beneficial for diagnosing and treating AVC. A diagnostic model was constructed based on a molecular prognostic score system using machine learning. The aforementioned immune cell infiltration may have a significant influence on the development and incidence of AVC.

FAM171B as a Novel Biomarker Mediates Tissue Immune Microenvironment in Pulmonary Arterial Hypertension

The purpose of this study was to uncover potential diagnostic indicators of pulmonary arterial hypertension (PAH), evaluate the function of immune cells in the pathogenesis of the disease, and find innovative treatment targets and medicines with the potential to enhance prognosis. Gene Expression Omnibus was utilized to acquire the PAH datasets. We recognized differentially expressed genes (DEGs) and investigated their functions utilizing R software. Weighted gene coexpression network analysis, least absolute shrinkage and selection operators, and support vector machines were used to identify biomarkers. The extent of immune cell infiltration in the normal and PAH tissues was determined using CIBERSORT. Additionally, the association between diagnostic markers and immune cells was analyzed. In this study, 258DEGs were used to analyze the disease ontology. Most DEGs were linked with atherosclerosis, arteriosclerotic cardiovascular disease, and lung disease, including obstructive lung disease. Gene set enrichment analysis revealed that compared to normal samples, results from PAH patients were mostly associated with ECM-receptor interaction, arrhythmogenic right ventricular cardiomyopathy, the Wnt signaling pathway, and focal adhesion. FAM171B was identified as a biomarker for PAH (area under the curve = 0.873). The mechanism underlying PAH may be mediated by nave CD4 T cells, resting memory CD4 T cells, resting NK cells, monocytes, activated dendritic cells, resting mast cells, and neutrophils, according to an investigation of immune cell infiltration. FAM171B expression was also associated with resting mast cells, monocytes, and CD8 T cells. The results suggest that PAH may be closely related to FAM171B with high diagnostic performance and associated with immune cell infiltration, suggesting that FAM171B may promote the progression of PAH by stimulating immune infiltration and immune response. This study provides valuable insights into the pathogenesis and treatment of PAH.

Advanced machine learning algorithms for flood susceptibility modeling - performance comparison: Red Sea, Egypt

Floods are among the most devastating environmental hazards that directly and indirectly affect people's lives and activities. In many countries, sustainable environmental management requires the assessment of floods and the likely flood-prone areas to avoid potential hazards. In this study, the performance and capabilities of seven machine learning algorithms (MLAs) for flood susceptibility mapping were tested, evaluated, and compared. These MLAs, including support vector machine (SVM), random forest (RF), multivariate adaptive regression spline (MARS), boosted regression tree (BRT), functional data analysis (FDA), general linear model (GLM), and multivariate discriminant analysis (MDA), were tested for the area between Safaga and Ras Gharib cities, Red Sea, Egypt. A geospatial database was developed with eleven flood-related factors, namely altitude, slope aspect, lithology, land use/land cover (LULC), slope length (LS), topographic wetness index (TWI), slope angle, profile curvature, plan curvature, stream power index (SPI), and hydrolithology units. In addition, 420 actual flooded areas were recorded from the study area to create a flood inventory map. The inventory data were randomly divided into training group with 70% and validation group with 30%. The flood-related factors were tested with a multicollinearity test, the variance inflation factor (VIF) was less than 2.135, the tolerance (TOL) was more than 0.468, and their importance was evaluated with a partial least squares (PLS) method. The results show that RF performed the best with the highest AUC (area under curve) value of 0.813, followed by GLM with 0.802, MARS with 0.801, BRT with 0.777, MDA with 0.768%, FDA with 0.763, and SVM with 0.733. The results of this study and the flood susceptibility maps could be useful for environmental mitigation, future development activities in the area, and flood control areas.

MX2: Identification and systematic mechanistic analysis of a novel immune-related biomarker for systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE) is an autoimmune disease that involves multiple organs. However, the current SLE-related biomarkers still lack sufficient sensitivity, specificity and predictive power for clinical application. Thus, it is significant to explore new immune-related biomarkers for SLE diagnosis and development.

Methods

We obtained seven SLE gene expression profile microarrays (GSE121239/11907/81622/65391/100163/45291/49454) from the GEO database. First, differentially expressed genes (DEGs) were screened using GEO2R, and SLE biomarkers were screened by performing WGCNA, Random Forest, SVM-REF, correlation with SLEDAI and differential gene analysis. Receiver operating characteristic curves (ROCs) and AUC values were used to determine the clinical value. The expression level of the biomarker was verified by RT‒qPCR. Subsequently, functional enrichment analysis was utilized to identify biomarker-associated pathways. ssGSEA, CIBERSORT, xCell and ImmuCellAI algorithms were applied to calculate the sample immune cell infiltration abundance. Single-cell data were analyzed for gene expression specificity in immune cells. Finally, the transcriptional regulatory network of the biomarker was constructed, and the corresponding therapeutic drugs were predicted.

Results

Multiple algorithms were screened together for a unique marker gene, MX2, and expression analysis of multiple datasets revealed that MX2 was highly expressed in SLE compared to the normal group (all P < 0.05), with the same trend validated by RT‒qPCR (P = 0.026). Functional enrichment analysis identified the main pathway of MX2 promotion in SLE as the NOD-like receptor signaling pathway (NES=2.492, P < 0.001, etc.). Immuno-infiltration analysis showed that MX2 was closely associated with neutrophils, and single-cell and transcriptomic data revealed that MX2 was specifically expressed in neutrophils. The NOD-like receptor signaling pathway was also remarkably correlated with neutrophils (r >0.3, P < 0.001, etc.). Most of the MX2-related interacting proteins were associated with SLE, and potential transcription factors of MX2 and its related genes were also significantly associated with the immune response.

Conclusion

Our study found that MX2 can serve as an immune-related biomarker for predicting the diagnosis and disease activity of SLE. It activates the NOD-like receptor signaling pathway and promotes neutrophil infiltration to aggravate SLE.

Radiomic assessment of oesophageal adenocarcinoma: a critical review of 18F-FDG PET/CT, PET/MRI and CT

OBJECTIVES

Radiomic models present an avenue to improve oesophageal adenocarcinoma assessment through quantitative medical image analysis. However, model selection is complicated by the abundance of available predictors and the uncertainty of their relevance and reproducibility. This analysis reviews recent research to facilitate precedent-based model selection for prospective validation studies.

METHODS

This analysis reviews research on 18F-FDG PET/CT, PET/MRI and CT radiomics in oesophageal adenocarcinoma between 2016 and 2021. Model design, testing and reporting are evaluated according to the Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis (TRIPOD) score and Radiomics Quality Score (RQS). Key results and limitations are analysed to identify opportunities for future research in the area.

RESULTS

Radiomic models of stage and therapeutic response demonstrated discriminative capacity, though clinical applications require greater sensitivity. Although radiomic models predict survival within institutions, generalisability is limited. Few radiomic features have been recommended independently by multiple studies.

CONCLUSIONS

Future research must prioritise prospective validation of previously proposed models to further clinical translation.

Single-Cell Sequencing Analysis and Multiple Machine Learning Methods Identified G0S2 and HPSE as Novel Biomarkers for Abdominal Aortic Aneurysm

Identifying biomarkers for abdominal aortic aneurysms (AAA) is key to understanding their pathogenesis, developing novel targeted therapeutics, and possibly improving patients outcomes and risk of rupture. Here, we identified AAA biomarkers from public databases using single-cell RNA-sequencing, weighted co-expression network (WGCNA), and differential expression analyses. Additionally, we used the multiple machine learning methods to identify biomarkers that differentiated large AAA from small AAA. Biomarkers were validated using GEO datasets. CIBERSORT was used to assess immune cell infiltration into AAA tissues and investigate the relationship between biomarkers and infiltrating immune cells. Therefore, 288 differentially expressed genes (DEGs) were screened for AAA and normal samples. The identified DEGs were mostly related to inflammatory responses, lipids, and atherosclerosis. For the large and small AAA samples, 17 DEGs, mostly related to necroptosis, were screened. As biomarkers for AAA, G0/G1 switch 2 (G0S2) (Area under the curve [AUC] = 0.861, 0.875, and 0.911, in GSE57691, GSE47472, and GSE7284, respectively) and for large AAA, heparinase (HPSE) (AUC = 0.669 and 0.754, in GSE57691 and GSE98278, respectively) were identified and further verified by qRT-PCR. Immune cell infiltration analysis revealed that the AAA process may be mediated by T follicular helper (Tfh) cells and the large AAA process may also be mediated by Tfh cells, M1, and M2 macrophages. Additionally, G0S2 expression was associated with neutrophils, activated and resting mast cells, M0 and M1 macrophages, regulatory T cells (Tregs), resting dendritic cells, and resting CD4 memory T cells. Moreover, HPSE expression was associated with M0 and M1 macrophages, activated and resting mast cells, Tregs, and resting CD4 memory T cells. Additional, G0S2 may be an effective diagnostic biomarker for AAA, whereas HPSE may be used to confer risk of rupture in large AAAs. Immune cells play a role in the onset and progression of AAA, which may improve its diagnosis and treatment.

Analysis of immune related gene expression profiles and immune cell components in patients with Barrett esophagus

Barrett's esophagus (BE) is a well-known precancerous condition of esophageal adenocarcinoma. However, the immune cells and immune related genes involved in BE development and progression are not fully understood. Therefore, our study attempted to investigate the roles of immune cells and immune related genes in BE patients. The raw gene expression data were downloaded from the GEO database. The limma package in R was used to screen differentially expressed genes (DEGs). Then we performed the least absolute shrinkage and selection operator (LASSO) and random forest (RF) analyses to screen key genes. The proportion of infiltrated immune cells was evaluated using the CIBERSORT algorithm between BE and normal esophagus (NE) samples. The spearman index was used to show the correlations of immune genes and immune cells. Receiver operating characteristic (ROC) curves were used to assess the diagnostic value of key genes in BE. A total of 103 differentially expressed immune-related genes were identified between BE samples and normal samples. Then, 7 genes (CD1A, LTF, FABP4, PGC, TCF7L2, INSR,SEMA3C) were obtained after Lasso analysis and RF modeling. CIBERSORT analysis revealed that resting CD4 T memory cells and gamma delta T cells were present at significantly lower levels in BE samples. Moreover, plasma cell and regulatory T cells were present at significantly higher levels in BE samples than in NE samples. INSR had the highest AUC values in ROC analysis. We identified 7 immune related genes and 4 different immune cells in our study, that may play vital roles in the occurrence and development of BE. Our findings improve the understanding of the molecular mechanisms of BE.

Revealing Potential Spinal Cord Injury Biomarkers and Immune Cell Infiltration Characteristics in Mice

Spinal cord injury (SCI) is a disabling condition with significant morbidity and mortality. Currently, no effective SCI treatment exists. This study aimed to identify potential biomarkers and characterize the properties of immune cell infiltration during this pathological event. To eliminate batch effects, we concurrently analyzed two mouse SCI datasets (GSE5296, GSE47681) from the GEO database. First, we identified differentially expressed genes (DEGs) using linear models for microarray data (LIMMA) and performed functional enrichment studies on those DEGs. Next, we employed bioinformatics and machine-learning methods to identify and define the characteristic genes of SCI. Finally, we validated them using immunofluorescence and qRT-PCR. Additionally, this study assessed the inflammatory status of SCI by identifying cell types using CIBERSORT. Furthermore, we investigated the link between key markers and infiltrating immune cells. In total, we identified 561 robust DEGs. We identified Rab20 and Klf6 as SCI-specific biomarkers and demonstrated their significance using qRT-PCR in the mouse model. According to the examination of immune cell infiltration, M0, M1, and M2 macrophages, along with naive CD8, dendritic cell-activated, and CD4 Follicular T cells may have a role in the progression of SCI. Therefore, Rab20 and Klf6 could be accessible targets for diagnosing and treating SCI. Moreover, as previously stated, immune cell infiltration may significantly impact the development and progression of SCI.