Introduction to Bioinformatics

Identification of candidate biomarkers correlated with the pathogenesis of breast cancer patients

Breast cancer (BC) is the second leading cause of cancer-related death in females, followed by lung cancer. Disadvantages exist in conventional diagnostic techniques of BC, such as radiation risk. The present study integrated bioinformatics analysis with machine learning to elucidate potential key candidate genes associated with the tumorigenesis of BC. Eleven datasets were downloaded from the Gene Expression Omnibus (GEO) database and were consolidated into two independent cohorts (training cohort and validation cohort) after batch-effect removal. We employed "limma" package to screen differentially expressed genes (DEGs) between BC and adjacent normal breast samples. Subsequently, the most reliable diagnostic indicators were identified utilizing LASSO-Logistic regression, SVM-RFE and multivariate stepwise Logistic regression analysis. Logistic model and nomogram were created based on these hub genes and applied in external validation cohort to verify the robustness of the model. As a result, a total of six hub genes connected with BC pathogenesis were identified, including CD300LG, IGSF10, FAM83D, MAMDC2, COMP and SEMA3G. Then, a diagnostic model of BC on the basis of these genes was established. ROC analysis of the diagnostic model illustrated that AUC of the training cohort was 0.978 (0.962, 0.995). In the validation cohort, AUC of training set and validation set were 0.936 (0.910, 0.961) and 0.921 (0.870, 0.972), respectively. This indicated that the model was reliable in separating BC patients from healthy individuals. The model may assist in early diagnosis of BC with implications for improving the prognosis of BC patients.

IL7R, GZMA and CD8A serve as potential molecular biomarkers for sepsis based on bioinformatics analysis

Purpose

Sepsis is an unusual systemic reaction to what is sometimes an otherwise ordinary infection, and it probably represents a pattern of response by the immune system to injury. However, the relationship between biomarkers and sepsis remains unclear. This study aimed to find potential molecular biomarkers, which could do some help to patients with sepsis.

Methods

The sepsis dataset GSE28750, GSE57065 was downloaded from the GEO database, and ten patients with or without sepsis from our hospital were admitted for RNA-seq and the differentially expressed genes (DEGs) were screened. The Metascape database was used for functional enrichment analysis and was used to found the differential gene list. Protein-protein interaction network was used and further analyzed by using Cytoscape and STRING. Logistic regression and Correlation analysis were used to find the potential molecular biomarkers.

Results

Taking the intersection of the three datasets yielded 287 differential genes. The enrichment results included Neutrophil degranulation, leukocyte activation, immune effectors process, positive regulation of immune response, regulation of leukocyte activation. The top 10 key genes of PPI connectivity were screened using cytoHubba plugin, which were KLRK1, KLRB1, IL7R, GZMA, CD27, PRF1, CD8A, CD2, IL2RB, and GZMB. All of the hub genes are higher expressed in health group of different databases. Logistic regression showed that IL7R, GZMA and CD8A proteins were analyzed and all of them were statistically significant. Correlation analysis showed that there was a statistically significant correlation between IL7R, GZMA and CD8A.

Conclusion

KLRK1, KLRB1, IL7R, GZMA, CD27, PRF1, CD8A, CD2, IL2RB, GZMB are key genes in sepsis, which associated with the development of sepsis. However, IL7R, GZMA and CD8A may serve as the attractively potential molecular biomarkers for sepsis.

Resistin as a potential diagnostic biomarker for sepsis: insights from DIA and ELISA analyses

PURPOSE

The primary objective of this investigation is to systematically screen and identify differentially expressed proteins (DEPs) within the plasma of individuals afflicted with sepsis. This endeavor employs both Data-Independent Acquisition (DIA) and enzyme-linked immunosorbent assay (ELISA) methodologies. The overarching goal is to furnish accessible and precise serum biomarkers conducive to the diagnostic discernment of sepsis.

METHOD

The study encompasses 53 sepsis patients admitted to the Affiliated Hospital of Southwest Medical University between January 2019 and December 2020, alongside a control cohort consisting of 16 individuals devoid of sepsis pathology. Subsequently, a subset comprising 10 randomly selected subjects from the control group and 22 from the sepsis group undergoes quantitative proteomic analysis via DIA. The acquired data undergoes Gene Ontology (GO) and Kyoto Encyclopedia of Genes (KEGG) analyses, facilitating the construction of a Protein-Protein Interaction (PPI) network to discern potential markers. Validation of core proteins is then accomplished through ELISA. Comparative analysis between the normal and sepsis groups ensues, characterized by Receiver Operating Characteristic (ROC) curve construction to evaluate diagnostic efficacy.

RESULT

A total of 187 DEPs were identified through bioinformatic methodologies. Examination reveals their predominant involvement in biological processes such as wound healing, coagulation, and blood coagulation. Functional pathway analysis further elucidates their engagement in the complement pathway and malaria. Resistin emerges as a candidate plasma biomarker, subsequently validated through ELISA. Notably, the protein exhibits significantly elevated levels in the serum of sepsis patients compared to the normal control group. ROC curve analysis underscores the robust diagnostic capacity of these biomarkers for sepsis.

CONCLUSION

Data-Independent Acquisition (DIA) and Enzyme-Linked Immunosorbent Assay (ELISA) show increased Resistin levels in sepsis patients, suggesting diagnostic potential, warranting further research.

LCN2 and ELANE overexpression induces sepsis

Sepsis is a syndrome characterized by a systemic inflammatory response due to the invasion of pathogenic microorganisms. The relationship between Lipocalin-2 (LCN2), elastase, neutrophil expressed (ELANE) and sepsis remains unclear. The sepsis datasets GSE137340 and GSE154918 profiles were downloaded from gene expression omnibus generated from GPL10558. Batch normalization, differentially expressed Genes (DEGs) screening, weighted gene co-expression network analysis (WGCNA), functional enrichment analysis, Gene Set Enrichment Analysis (GSEA), immune infiltration analysis, construction and analysis of protein-protein interaction (PPI) networks, Comparative Toxicogenomics Database (CTD) analysis were performed. Gene expression heatmaps were generated. TargetScan was used to screen miRNAs of DEGs. 328 DEGs were identified. According to Gene Ontology (GO), in the Biological Process analysis, they were mainly enriched in immune response, apoptosis, inflammatory response, and immune response regulation signaling pathways. In cellular component analysis, they were mainly enriched in vesicles, cytoplasmic vesicles, and secretory granules. In Molecular Function analysis, they were mainly concentrated in hemoglobin binding, Toll-like receptor binding, immunoglobulin binding, and RAGE receptor binding. In Kyoto Encyclopedia of Genes and Genomes (KEGG), they were mainly enriched in NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, TNF signaling pathway, P53 signaling pathway, and legionellosis. Seventeen modules were generated. The PPI network identified 4 core genes (MPO, ELANE, CTSG, LCN2). Gene expression heatmaps revealed that core genes (MPO, ELANE, CTSG, LCN2) were highly expressed in sepsis samples. CTD analysis found that MPO, ELANE, CTSG and LCN2 were associated with sepsis, peritonitis, meningitis, pneumonia, infection, and inflammation. LCN2 and ELANE are highly expressed in sepsis and may serve as molecular targets.

Identification of diagnostic biomarkers in Alzheimer's disease by integrated bioinformatic analysis and machine learning strategies

Background

Alzheimer's disease (AD) is the most prevalent form of dementia, and is becoming one of the most burdening and lethal diseases. More useful biomarkers for diagnosing AD and reflecting the disease progression are in need and of significance.

Methods

The integrated bioinformatic analysis combined with machine-learning strategies was applied for exploring crucial functional pathways and identifying diagnostic biomarkers of AD. Four datasets (GSE5281, GSE131617, GSE48350, and GSE84422) with samples of AD frontal cortex are integrated as experimental datasets, and another two datasets (GSE33000 and GSE44772) with samples of AD frontal cortex were used to perform validation analyses. Functional Correlation enrichment analyses were conducted based on Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Reactome database to reveal AD-associated biological functions and key pathways. Four models were employed to screen the potential diagnostic biomarkers, including one bioinformatic analysis of Weighted gene co-expression network analysis (WGCNA)and three machine-learning algorithms: Least absolute shrinkage and selection operator (LASSO), support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF) analysis. The correlation analysis was performed to explore the correlation between the identified biomarkers with CDR scores and Braak staging.

Results

The pathways of the immune response and oxidative stress were identified as playing a crucial role during AD. Thioredoxin interacting protein (TXNIP), early growth response 1 (EGR1), and insulin-like growth factor binding protein 5 (IGFBP5) were screened as diagnostic markers of AD. The diagnostic efficacy of TXNIP, EGR1, and IGFBP5 was validated with corresponding AUCs of 0.857, 0.888, and 0.856 in dataset GSE33000, 0.867, 0.909, and 0.841 in dataset GSE44770. And the AUCs of the combination of these three biomarkers as a diagnostic tool for AD were 0.954 and 0.938 in the two verification datasets.

Conclusion

The pathways of immune response and oxidative stress can play a crucial role in the pathogenesis of AD. TXNIP, EGR1, and IGFBP5 are useful biomarkers for diagnosing AD and their mRNA level may reflect the development of the disease by correlation with the CDR scores and Breaking staging.

A metabolic reprogramming-related prognostic risk model for clear cell renal cell carcinoma: From construction to preliminary application

Metabolic reprogramming is one of the characteristics of clear cell renal cell carcinoma (ccRCC). Although some treatments associated with the metabolic reprogramming for ccRCC have been identified, remain still lacking. In this study, we identified the differentially expressed genes (DEGs) associated with clinical traits with a total of 965 samples via DEG analysis and weighted correlation network analysis (WGCNA), screened the prognostic metabolism-related genes, and constructed the risk score prognostic models. We took the intersection of DEGs with significant difference coexpression modules and received two groups of intersection genes that were connected with metabolism via functional enrichment analysis. Then we respectively screened prognostic metabolic-related genes from the genes of the two intersection groups and constructed the risk score prognostic models. Compared with the predicted effect of clinical grade and stage for ccRCC patients, finally, we selected the model constructed with genes of ABAT, ALDH6A1, CHDH, EPHX2, ETNK2, and FBP1. The risk scores of the prognostic model were significantly related to overall survival (OS) and could serve as an independent prognostic factor. The Kaplan-Meier analysis and ROC curves revealed that the model efficiently predicts prognosis in the TCGA-KIRC cohort and the validation cohort. Then we investigated the potential underlying mechanism and sensitive drugs between high- and low-risk groups. The six key genes were significantly linked with worse OS and were downregulated in ccRCC, we confirmed the results in clinical samples. These results demonstrated the efficacy and robustness of the risk score prognostic model, based on the characteristics of metabolic reprogramming in ccRCC, and the key genes used in constructing the model also could develop into targets of molecular therapy for ccRCC.

Exploring biomarkers for ischemic stroke through integrated microarray data analysis

Stroke is the third leading cause of disability-adjusted life years worldwide, and drugs available for its treatment are limited. This study aimed to explore high-confidence candidate genes associated with ischemic stroke (IS) through bioinformatics analysis and identify potential diagnostic biomarkers and gene-drug interactions. Weighted gene coexpression network analysis (WGCNA) and differentially expressed genes (DEGs) were integrated to identify overlapping genes. Then, high-confidence candidate genes were screened by least absolute shrinkage and selection operator (LASSO) regression. Receiver operating characteristic (ROC) curves were used to evaluate the diagnostic value of high-confidence candidate genes as biomarkers for IS. The NetworkAnalyst database was used to construct the TF-gene network and miRNA-TF regulatory network of the high-confidence candidate genes. The DGIdb database was used to identified gene-drug interactions. Through the comprehensive analysis of GSE58294 and GSE16561, 10 high-confidence candidate genes were identified by LASSO regression: ARG1, LY96, ABCA1, SLC22A4, CD163, TPM2, SLC25A42, ID3, FAM102A and CD79B. FAM102A had the highest diagnostic value, and the area under curve (AUC), sensitivity and specificity values were 0.974, 0.919 and 0.936, respectively. The HPA database demonstrated that 10 high-confidence candidate genes were expressed in the brain and blood in normal humans. Finally, DGIdb database analysis identified 8 gene-drug interactions. We identified IS-related diagnostic biomarkers and gene-drug interactions that potentially provide new insights into the diagnosis and treatment of IS.

Identification of diagnosis and prognosis gene markers in B-ALL with ETV6-RUNX1 fusion by integrated bioinformatics analysis

B-acute lymphoblastic leukemia (B-ALL) is characterized by clonal expansion of immature B-lymphocytes in the bone marrow, blood, or other tissues. Chromosomal translocations have often been reported in B-ALL, which are important for its prognosis. B-ALL patients with ETV6-RUNX1 fusion have favorable outcomes, but the mechanisms remain to be clarified. In the present study, we crossed the selected WGCNA module genes and differential expression genes to obtain core genes, and random forest algorithm, a type of supervised learning analysis, was conducted to evaluate the importance of those core genes in distinguishing B-ALL samples with ETV6-RUNX2 fusion with extracting 5 genes as gene markers for ETV6-RUNX2 fusion. Moreover, we calculated the immune infiltration profiles and screened out the ETV6-RUNX2 association immune cells using the CIBERSORT algorithm. In conclusion, combined with various solid informatics methods, we depicted the underlying molecular and immune mechanism of ETV6-RUNX2 fusion and providing potential biological targets for diagnosing and treating B-ALL in the future.

Uncovering the Pathogenesis of Orofacial Clefts Using Bioinformatics Analysis

OBJECTIVE

Many genes have been found to be associated with the occurrence of the orofacial clefts (OFC). The links between these pathogenic genes are rarely studied. In this study, bioinformatics analysis were performed in order to find associations between OFC-related genes and provide new ideas for etiology study of OFCs.

METHODS

Orofacial clefts-related genes were searched and identified from the Online Mendelian Inheritance of Man (OMIM.org). These genes were then analyzed by bioinformatics methods, including protein-protein interaction network, functional enrichment analysis, module analysis, and hub genes analysis.

RESULTS

After searching the database of OMIM.org and removing duplicate results, 279 genes were finally obtained. These genes were involved to 369 pathways in biological process, 56 in cell component, 64 in molecular function, and 45 in the Kyoto Encyclopedia of Genes and Genomes. Most identified genes were significantly enriched in embryonic appendage morphogenesis (29.17%), embryonic limb morphogenesis (6.06%), and limb development (4.33%) for biological process (Fig. 5A); ciliary tip (42.86%), MKS complex (28.57%), ciliary basal body (14.29%), and ciliary membrane (14.29%) for cell component. The top 10 hub genes were identified, including SHH, GLI2, PTCH1, SMAD4, FGFR1, BMP4, SOX9, SOX2, RUNX2, and CDH1.

CONCLUSIONS

Bioinformatics methods were used to analyze OFC-related genes in this study, including hub gene identifying and analysis, protein-protein interaction network construction, and functional enrichment analysis. Several potential mechanisms related to occurrence of OFCs were also discussed. These results may be helpful for further studies of the etiology of OFC.

Novel Plasma Biomarker-Based Model for Predicting Acute Kidney Injury After Cardiac Surgery: A Case Control Study

Introduction:

Acute kidney injury (AKI) after cardiac surgery is independently associated with a prolonged hospital stay, increased cost of care, and increased post-operative mortality. Delayed elevation of serum creatinine (SCr) levels requires novel biomarkers to provide a prediction of AKI after cardiac surgery. Our objective was to find a novel blood biomarkers combination to construct a model for predicting AKI after cardiac surgery and risk stratification.

Methods:

This was a case-control study. Weighted Gene Co-expression Network Analysis (WGCNA) was applied to Gene Expression Omnibus (GEO) dataset GSE30718 to seek potential biomarkers associated with AKI. We measured biomarker levels in venous blood samples of 67 patients with AKI after cardiac surgery and 59 control patients in two cohorts. Clinical data were collected. We developed a multi-biomarker model for predicting cardiac-surgery-associated AKI and compared it with a traditional clinical-factor-based model.

Results:

From bioinformatics analysis and previous articles, we found 6 potential plasma biomarkers for the prediction of AKI. Among them, 3 biomarkers, such as growth differentiation factor 15 (GDF15), soluble suppression of tumorigenicity 2 (ST2, IL1RL1), and soluble urokinase plasminogen activator receptor (uPAR) were found to have prediction ability for AKI (area under the curve [AUC] > 0.6) in patients undergoing cardiac surgery. They were then incorporated into a multi-biomarker model for predicting AKI (C-statistic: 0.84, Brier 0.15) which outperformed the traditional clinical-factor-based model (C-statistic: 0.73, Brier 0.16).

Conclusion:

Our research validated a promising plasma multi-biomarker model for predicting AKI after cardiac surgery.

SOX10 Knockdown Inhibits Melanoma Cell Proliferation via Notch Signaling Pathway

Purpose

Melanoma is a serious and malignant disease worldwide. Seeking diagnostic markers and potential therapeutic targets is urgent for melanoma treatment. SOX10, a member of the SoxE family of genes, is a transcription factor which can regulate the transcription of a wide variety of genes in multiple cellular processes.

Methods

The mRNA level and protein expression of SOX10 is confirmed by bioinformatic analysis and IHC staining. MTT, clone formation and EdU analysis showed that SOX10 knockdown (KD) could significantly inhibit melanoma cell proliferation. FACS analysis showed that SOX10 KD could markedly enhance the level of cell apoptosis. The downstream target signaling pathway is predicted by RNA-seq based on the public GEO database. The activation of Notch signaling mediated by SOX10 is tested by qPCR and Western blot.

Results

Ectopic upregulation of SOX10 was found in melanoma patient tissues compared to normal nevus tissues in mRNA and protein levels. Furthermore, both mRNA and protein level of SOX10 were negatively correlated with melanoma patient's prognosis. SOX10 knockdown could obviously suppress the proliferation ability of melanoma cells by inactivating Notch signaling pathway.

Conclusion

Our study confirmed that SOX10 is an oncogene and activate Notch signaling pathway, which suggests the potential treatment for melanoma patients by target SOX10/Notch axis.

Incorporating Machine Learning into Established Bioinformatics Frameworks

The exponential growth of biomedical data in recent years has urged the application of numerous machine learning techniques to address emerging problems in biology and clinical research. By enabling the automatic feature extraction, selection, and generation of predictive models, these methods can be used to efficiently study complex biological systems. Machine learning techniques are frequently integrated with bioinformatic methods, as well as curated databases and biological networks, to enhance training and validation, identify the best interpretable features, and enable feature and model investigation. Here, we review recently developed methods that incorporate machine learning within the same framework with techniques from molecular evolution, protein structure analysis, systems biology, and disease genomics. We outline the challenges posed for machine learning, and, in particular, deep learning in biomedicine, and suggest unique opportunities for machine learning techniques integrated with established bioinformatics approaches to overcome some of these challenges.