FlowAtlas: an interactive tool for high-dimensional immunophenotyping analysis bridging FlowJo with computational tools in Julia.

As the dimensionality, throughput and complexity of cytometry data increases, so does the demand for user-friendly, interactive analysis tools that leverage high-performance machine learning frameworks. Here we introduce FlowAtlas: an interactive web application that enables dimensionality reduction of cytometry data without down-sampling and that is compatible with datasets stained with non-identical panels. FlowAtlas bridges the user-friendly environment of FlowJo and computational tools in Julia developed by the scientific machine learning community, eliminating the need for coding and bioinformatics expertise. New population discovery and detection of rare populations in FlowAtlas is intuitive and rapid. We demonstrate the capabilities of FlowAtlas using a human multi-tissue, multi-donor immune cell dataset, highlighting key immunological findings. FlowAtlas is available at https://github.com/gszep/FlowAtlas.jl.git.

Systems biology-based analysis exploring shared biomarkers and pathogenesis of myocardial infarction combined with osteoarthritis.

Background: More and more evidence supports the association between myocardial infarction (MI) and osteoarthritis (OA). The purpose of this study is to explore the shared biomarkers and pathogenesis of MI complicated with OA by systems biology. Methods: Gene expression profiles of MI and OA were downloaded from the Gene Expression Omnibus (GEO) database. The Weighted Gene Co-Expression Network Analysis (WGCNA) and differentially expressed genes (DEGs) analysis were used to identify the common DEGs. The shared genes related to diseases were screened by three public databases, and the protein-protein interaction (PPI) network was built. GO and KEGG enrichment analyses were performed on the two parts of the genes respectively. The hub genes were intersected and verified by Least absolute shrinkage and selection operator (LASSO) analysis, receiver operating characteristic (ROC) curves, and single-cell RNA sequencing analysis. Finally, the hub genes differentially expressed in primary cardiomyocytes and chondrocytes were verified by RT-qPCR. The immune cell infiltration analysis, subtypes analysis, and transcription factors (TFs) prediction were carried out. Results: In this study, 23 common DEGs were obtained by WGCNA and DEGs analysis. In addition, 199 common genes were acquired from three public databases by PPI. Inflammation and immunity may be the common pathogenic mechanisms, and the MAPK signaling pathway may play a key role in both disorders. DUSP1, FOS, and THBS1 were identified as shared biomarkers, which is entirely consistent with the results of single-cell RNA sequencing analysis, and furher confirmed by RT-qPCR. Immune infiltration analysis illustrated that many types of immune cells were closely associated with MI and OA. Two potential subtypes were identified in both datasets. Furthermore, FOXC1 may be the crucial TF, and the relationship of TFs-hub genes-immune cells was visualized by the Sankey diagram, which could help discover the pathogenesis between MI and OA. Conclusion: In summary, this study first revealed 3 (DUSP1, FOS, and THBS1) novel shared biomarkers and signaling pathways underlying both MI and OA. Additionally, immune cells and key TFs related to 3 hub genes were examined to further clarify the regulation mechanism. Our study provides new insights into shared molecular mechanisms between MI and OA.

Towards an enhanced understanding of osteoanabolic effects of PTH-induced microRNAs on osteoblasts using a bioinformatic approach.

In this study, we used a bioinformatic approach to construct a miRNA-target gene interaction network potentially involved in the anabolic effect of parathyroid hormone analogue teriparatide [PTH (1-34)] on osteoblasts. We extracted a dataset of 26 microRNAs (miRNAs) from previously published studies and predicted miRNA target interactions (MTIs) using four software tools: DIANA, miRWalk, miRDB, and TargetScan. By constructing an interactome of PTH-regulated miRNAs and their predicted target genes, we elucidated signaling pathways regulating pluripotency of stem cells, the Hippo signaling pathway, and the TGF-beta signaling pathway as the most significant pathways in the effects of PTH on osteoblasts. Furthermore, we constructed intersection of MTI networks for these three pathways and added validated interactions. There are 8 genes present in all three selected pathways and a set of 18 miRNAs are predicted to target these genes, according to literature data. The most important genes in all three pathways were BMPR1A, BMPR2 and SMAD2 having the most interactions with miRNAs. Among these miRNAs, only miR-146a-5p and miR-346 have validated interactions in these pathways and were shown to be important regulators of these pathways. In addition, we also propose miR-551b-5p and miR-338-5p for further experimental validation, as they have been predicted to target important genes in these pathways but none of their target interactions have yet been verified. Our wet-lab experiment on miRNAs differentially expressed between PTH (1-34) treated and untreated mesenchymal stem cells supports miR-186-5p from the literature obtained data as another prominent miRNA. The meticulous selection of miRNAs outlined will significantly support and guide future research aimed at discovering and understanding the crucial pathways of osteoanabolic PTH-epigenetic effects on osteoblasts. Additionally, they hold potential for the discovery of new PTH target genes, innovative biomarkers for the effectiveness and safety of osteoporosis-affected treatment, as well as novel therapeutic targets.

Microbiome-host interactions in the pathogenesis of acute exacerbation of chronic obstructive pulmonary disease.

Objective: To explore the underlying mechanisms the airway microbiome contributes to Acute Exacerbation of Chronic Obstructive Pulmonary Disease(AECOPD). Methods: We enrolled 31 AECOPD patients and 26 stable COPD patients, their sputum samples were collected for metagenomic and RNA sequencing, and then subjected to bioinformatic analyses. The expression of host genes was validated by Quantitative Real-time PCR(qPCR) using the same batch of specimens. Results: Our results indicated a higher expression of Rothia mucilaginosa(p=0.015) in the AECOPD group and Haemophilus influenzae(p=0.005) in the COPD group. The Different expressed genes(DEGs) detected were significantly enriched in "type I interferon signaling pathway"(p<0.001, q=0.001) in gene function annotation, and "Cytosolic DNA-sensing pathway"(p=0.002, q=0.024), "Toll-like receptor signaling pathway"(p=0.006, q=0.045), and "TNF signaling pathway"(p=0.006, q=0.045) in KEGG enrichment analysis. qPCR amplification experiment verified that the expression of OASL and IL6 increased significantly in the AECOPD group. Conclusion: Pulmonary bacteria dysbiosis may regulate the pathogenesis of AECOPD through innate immune system pathways like type I interferon signaling pathway and Toll-like receptor signaling pathway.

Holimap: an accurate and efficient method for solving stochastic gene network dynamics.

Gene-gene interactions are crucial to the control of sub-cellular processes but our understanding of their stochastic dynamics is hindered by the lack of simulation methods that can accurately and efficiently predict how the distributions of gene product numbers vary across parameter space. To overcome these difficulties, here we present Holimap (high-order linear-mapping approximation), an approach that approximates the protein or mRNA number distributions of a complex gene regulatory network by the distributions of a much simpler reaction system. We demonstrate Holimap's computational advantages over conventional methods by applying it to predict the stochastic time-dependent dynamics of various gene networks, including transcriptional networks ranging from simple autoregulatory loops to complex randomly connected networks, post-transcriptional networks, and post-translational networks. Holimap is ideally suited to study how the intricate network of gene-gene interactions results in precise coordination and control of gene expression.

OmicScope unravels systems-level insights from quantitative proteomics data.

Shotgun proteomics analysis presents multifaceted challenges, demanding diverse tool integration for insights. Addressing this complexity, OmicScope emerges as an innovative solution for quantitative proteomics data analysis. Engineered to handle various data formats, it performs data pre-processing - including joining replicates, normalization, data imputation - and conducts differential proteomics analysis for both static and longitudinal experimental designs. Empowered by Enrichr with over 224 databases, OmicScope performs Over Representation Analysis (ORA) and Gene Set Enrichment Analysis (GSEA). Additionally, its Nebula module facilitates meta-analysis from independent datasets, providing a systems biology approach for enriched insights. Complete with a data visualization toolkit and accessible as Python package and a web application, OmicScope democratizes proteomics analysis, offering an efficient and high-quality pipeline for researchers.

Tissue characterization at an enhanced resolution across spatial omics platforms with deep generative model.

Recent advances in spatial omics have expanded the spectrum of profiled molecular categories beyond transcriptomics. However, many of these technologies are constrained by limited spatial resolution, hindering our ability to deeply characterize intricate tissue architectures. Existing computational methods primarily focus on the resolution enhancement of transcriptomics data, lacking the adaptability to address the emerging spatial omics technologies that profile various omics types. Here, we introduce soScope, a unified generative framework designed to enhance data quality and spatial resolution for molecular profiles obtained from diverse spatial technologies. soScope aggregates multimodal tissue information from omics, spatial relations and images, and jointly infers omics profiles at enhanced resolutions with omics-specific modeling through distribution priors. With comprehensive evaluations on diverse spatial omics platforms, including Visium, Xenium, spatial-CUT&Tag, and slide-DNA/RNA-seq, soScope improves performances in identifying biologically meaningful intestine and kidney architectures, revealing embryonic heart structure that cannot be resolved at the original resolution and correcting sample and technical biases arising from sequencing and sample processing. Furthermore, soScope extends to spatial multiomics technology spatial-CITE-seq and spatial ATAC-RNA-seq, leveraging cross-omics reference for simultaneous multiomics enhancement. soScope provides a versatile tool to improve the utilization of continually expanding spatial omics technologies and resources.

Phenotype prediction using biologically interpretable neural networks on multi-cohort multi-omics data.

Integrating multi-omics data into predictive models has the potential to enhance accuracy, which is essential for precision medicine. In this study, we developed interpretable predictive models for multi-omics data by employing neural networks informed by prior biological knowledge, referred to as visible networks. These neural networks offer insights into the decision-making process and can unveil novel perspectives on the underlying biological mechanisms associated with traits and complex diseases. We tested the performance, interpretability and generalizability for inferring smoking status, subject age and LDL levels using genome-wide RNA expression and CpG methylation data from the blood of the BIOS consortium (four population cohorts, Ntotal = 2940). In a cohort-wise cross-validation setting, the consistency of the diagnostic performance and interpretation was assessed. Performance was consistently high for predicting smoking status with an overall mean AUC of 0.95 (95% CI: 0.90-1.00) and interpretation revealed the involvement of well-replicated genes such as AHRR, GPR15 and LRRN3. LDL-level predictions were only generalized in a single cohort with an R2 of 0.07 (95% CI: 0.05-0.08). Age was inferred with a mean error of 5.16 (95% CI: 3.97-6.35) years with the genes COL11A2, AFAP1, OTUD7A, PTPRN2, ADARB2 and CD34 consistently predictive. For both regression tasks, we found that using multi-omics networks improved performance, stability and generalizability compared to interpretable single omic networks. We believe that visible neural networks have great potential for multi-omics analysis; they combine multi-omic data elegantly, are interpretable, and generalize well to data from different cohorts.

Immunoinformatics design and synthesis of a multi-epitope vaccine against Helicobacter pylori based on lipid nanoparticles.

Helicobacter pylori (H. pylori) is responsible for various chronic or acute diseases, such as stomach ulcers, dyspepsia, peptic ulcers, gastroesophageal reflux, gastritis, lymphoma, and stomach cancers. Although specific drugs are available to treat the bacterium's harmful effects, there is an urgent need to develop a preventive or therapeutic vaccine. Therefore, the current study aims to create a multi-epitope vaccine against H. pylori using lipid nanoparticles. Five epitopes from five target proteins of H. pylori, namely, Urease, CagA, HopE, SabA, and BabA, were used. Immunogenicity, MHC (Major Histocompatibility Complex) bonding, allergenicity, toxicity, physicochemical analysis, and global population coverage of the entire epitopes and final construct were carefully examined. The study involved using various bioinformatic web tools to accomplish the following tasks: modeling the three-dimensional structure of a set of epitopes and the final construct and docking them with Toll-Like Receptor 4 (TLR4). In the experimental phase, the final multi-epitope construct was synthesized using the solid phase method, and it was then enclosed in lipid nanoparticles. After synthesizing the construct, its loading, average size distribution, and nanoliposome shape were checked using Nanodrop at 280 nm, dynamic light scattering (DLS), and atomic force microscope (AFM). The designed vaccine has been confirmed to be non-toxic and anti-allergic. It can bind with different MHC alleles at a rate of 99.05%. The construct loading was determined to be about 91%, with an average size of 54 nm. Spherical shapes were also observed in the AFM images. Further laboratory tests are necessary to confirm the safety and immunogenicity of the multi-epitope vaccine.

Investigating the WNT and TGF-beta pathways alterations and tumor mutant burden in young-onset colorectal cancer.

Colorectal cancer (CRC) is the third most common cancer in the United States. Recent epidemiological evidence demonstrates an increasing incidence of young-onset CRC cases, defined as CRC cases in individuals 50 years old or younger. Studies have established that alterations in both the WNT and TGF-Beta signaling pathways have contributed to CRC development. While this is well understood, the comprehensive analysis of WNT and TGF-Beta pathway alterations in young-onset CRC cases has yet to be investigated. Here, we conducted a comprehensive bioinformatics analysis of mutations associated with each of the WNT and TGF-Beta signaling pathways according to age (≤ 50 years old versus > 50 years old) utilizing published genomic data from the cBioPortal. Chi-square results demonstrated no significant difference in WNT alterations between young-onset CRC and those > 50 years old. However, across all age groups, WNT alterations were frequently found in rectal cancers. We also found that WNT alterations were associated with better outcomes. The mutations associated with TGF-beta were observed at a higher rate in older CRC patients when compared to those ≤ 50 years old. Additionally, these mutations were found more frequently in colon primaries.

Identification and verification of immune-related genes for diagnosing the progression of atherosclerosis and metabolic syndrome.

BACKGROUND: Atherosclerosis and metabolic syndrome are the main causes of cardiovascular events, but their underlying mechanisms are not clear. In this study, we focused on identifying genes associated with diagnostic biomarkers and effective therapeutic targets associated with these two diseases. METHODS: Transcriptional data sets of atherosclerosis and metabolic syndrome were obtained from GEO database. The differentially expressed genes were analyzed by RStudio software, and the function-rich and protein-protein interactions of the common differentially expressed genes were analyzed.Furthermore, the hub gene was screened by Cytoscape software, and the immune infiltration of hub gens was analyzed. Finally, relevant clinical blood samples were collected for qRT-PCR verification of the three most important hub genes. RESULTS: A total of 1242 differential genes (778 up-regulated genes and 464 down-regulated genes) were screened from GSE28829 data set. A total of 1021 differential genes (492 up-regulated genes and 529 down-regulated genes) were screened from the data set GSE98895. Then 23 up-regulated genes and 11 down-regulated genes were screened by venn diagram. Functional enrichment analysis showed that cytokines and immune activation were involved in the occurrence and development of these two diseases. Through the construction of the Protein-Protein Interaction(PPI) network and Cytoscape software analysis, we finally screened 10 hub genes. The immune infiltration analysis was further improved. The results showed that the infiltration scores of 7 kinds of immune cells in GSE28829 were significantly different among groups (Wilcoxon Test < 0.05), while in GSE98895, the infiltration scores of 4 kinds of immune cells were significantly different between groups (Wilcoxon Test < 0.05). Spearman method was used to analyze the correlation between the expression of 10 key genes and 22 kinds of immune cell infiltration scores in two data sets. The results showed that there were 42 pairs of significant correlations between 10 genes and 22 kinds of immune cells in GSE28829 (|Cor| > 0.3 & P < 0.05). There were 41 pairs of significant correlations between 10 genes and 22 kinds of immune cells in GSE98895 (|Cor| > 0.3 & P < 0.05). Finally, our results identified 10 small molecules with the highest absolute enrichment value, and the three most significant key genes (CX3CR1, TLR5, IL32) were further verified in the data expression matrix and clinical blood samples. CONCLUSION: We have established a co-expression network between atherosclerotic progression and metabolic syndrome, and identified key genes between the two diseases. Through the method of bioinformatics, we finally obtained 10 hub genes in As and MS, and selected 3 of the most significant genes (CX3CR1, IL32, TLR5) for blood PCR verification. This may be helpful to provide new research ideas for the diagnosis and treatment of AS complicated with MS.

Interactive molecular causal networks of hypertension using a fast machine learning algorithm MRdualPC.

BACKGROUND: Understanding the complex interactions between genes and their causal effects on diseases is crucial for developing targeted treatments and gaining insight into biological mechanisms. However, the analysis of molecular networks, especially in the context of high-dimensional data, presents significant challenges. METHODS: This study introduces MRdualPC, a computationally tractable algorithm based on the MRPC approach, to infer large-scale causal molecular networks. We apply MRdualPC to investigate the upstream causal transcriptomics influencing hypertension using a comprehensive dataset of kidney genome and transcriptome data. RESULTS: Our algorithm proves to be 100 times faster than MRPC on average in identifying transcriptomics drivers of hypertension. Through clustering, we identify 63 modules with causal driver genes, including 17 modules with extensive causal networks. Notably, we find that genes within one of the causal networks are associated with the electron transport chain and oxidative phosphorylation, previously linked to hypertension. Moreover, the identified causal ancestor genes show an over-representation of blood pressure-related genes. CONCLUSIONS: MRdualPC has the potential for broader applications beyond gene expression data, including multi-omics integration. While there are limitations, such as the need for clustering in large gene expression datasets, our study represents a significant advancement in building causal molecular networks, offering researchers a valuable tool for analyzing big data and investigating complex diseases.

PLEKv2: predicting lncRNAs and mRNAs based on intrinsic sequence features and the coding-net model.

BACKGROUND: Long non-coding RNAs (lncRNAs) are RNA transcripts of more than 200 nucleotides that do not encode canonical proteins. Their biological structure is similar to messenger RNAs (mRNAs). To distinguish between lncRNA and mRNA transcripts quickly and accurately, we upgraded the PLEK alignment-free tool to its next version, PLEKv2, and constructed models tailored for both animals and plants. RESULTS: PLEKv2 can achieve 98.7% prediction accuracy for human datasets. Compared with classical tools and deep learning-based models, this is 8.1%, 3.7%, 16.6%, 1.4%, 4.9%, and 48.9% higher than CPC2, CNCI, Wen et al.'s CNN, LncADeep, PLEK, and NcResNet, respectively. The accuracy of PLEKv2 was > 90% for cross-species prediction. PLEKv2 is more effective and robust than CPC2, CNCI, LncADeep, PLEK, and NcResNet for primate datasets (including chimpanzees, macaques, and gorillas). Moreover, PLEKv2 is not only suitable for non-human primates that are closely related to humans, but can also predict the coding ability of RNA sequences in plants such as Arabidopsis. CONCLUSIONS: The experimental results illustrate that the model constructed by PLEKv2 can distinguish lncRNAs and mRNAs better than PLEK. The PLEKv2 software is freely available at https://sourceforge.net/projects/plek2/ .

Integrating bioinformatics and machine learning methods to analyze diagnostic biomarkers for HBV-induced hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a malignant tumor. It is estimated that approximately 50-80% of HCC cases worldwide are caused by hepatitis b virus (HBV) infection, and other pathogenic factors have been shown to promote the development of HCC when coexisting with HBV. Understanding the molecular mechanisms of HBV-induced hepatocellular carcinoma (HBV-HCC) is crucial for the prevention, diagnosis, and treatment of the disease. In this study, we analyzed the molecular mechanisms of HBV-induced HCC by combining bioinformatics and deep learning methods. Firstly, we collected a gene set related to HBV-HCC from the GEO database, performed differential analysis and WGCNA analysis to identify genes with abnormal expression in tumors and high relevance to tumors. We used three deep learning methods, Lasso, random forest, and SVM, to identify key genes RACGAP1, ECT2, and NDC80. By establishing a diagnostic model, we determined the accuracy of key genes in diagnosing HBV-HCC. In the training set, RACGAP1(AUC:0.976), ECT2(AUC:0.969), and NDC80 (AUC: 0.976) showed high accuracy. They also exhibited good accuracy in the validation set: RACGAP1(AUC:0.878), ECT2(AUC:0.731), and NDC80(AUC:0.915). The key genes were found to be highly expressed in liver cancer tissues compared to normal liver tissues, and survival analysis indicated that high expression of key genes was associated with poor prognosis in liver cancer patients. This suggests a close relationship between key genes RACGAP1, ECT2, and NDC80 and the occurrence and progression of HBV-HCC. Molecular docking results showed that the key genes could spontaneously bind to the anti-hepatocellular carcinoma drugs Lenvatinib, Regorafenib, and Sorafenib with strong binding activity. Therefore, ECT2, NDC80, and RACGAP1 may serve as potential biomarkers for the diagnosis of HBV-HCC and as targets for the development of targeted therapeutic drugs.

Single-cell tumor heterogeneity landscape of hepatocellular carcinoma: unraveling the pro-metastatic subtype and its interaction loop with fibroblasts.

BACKGROUND: Tumor heterogeneity presents a formidable challenge in understanding the mechanisms driving tumor progression and metastasis. The heterogeneity of hepatocellular carcinoma (HCC) in cellular level is not clear. METHODS: Integration analysis of single-cell RNA sequencing data and spatial transcriptomics data was performed. Multiple methods were applied to investigate the subtype of HCC tumor cells. The functional characteristics, translation factors, clinical implications and microenvironment associations of different subtypes of tumor cells were analyzed. The interaction of subtype and fibroblasts were analyzed. RESULTS: We established a heterogeneity landscape of HCC malignant cells by integrated 52 single-cell RNA sequencing data and 5 spatial transcriptomics data. We identified three subtypes in tumor cells, including ARG1+ metabolism subtype (Metab-subtype), TOP2A+ proliferation phenotype (Prol-phenotype), and S100A6+ pro-metastatic subtype (EMT-subtype). Enrichment analysis found that the three subtypes harbored different features, that is metabolism, proliferating, and epithelial-mesenchymal transition. Trajectory analysis revealed that both Metab-subtype and EMT-subtype originated from the Prol-phenotype. Translation factor analysis found that EMT-subtype showed exclusive activation of SMAD3 and TGF-ß signaling pathway. HCC dominated by EMT-subtype cells harbored an unfavorable prognosis and a deserted microenvironment. We uncovered a positive loop between tumor cells and fibroblasts mediated by SPP1-CD44 and CCN2/TGF-ß-TGFBR1 interaction pairs. Inhibiting CCN2 disrupted the loop, mitigated the transformation to EMT-subtype, and suppressed metastasis. CONCLUSION: By establishing a heterogeneity landscape of malignant cells, we identified a three-subtype classification in HCC. Among them, S100A6+ tumor cells play a crucial role in metastasis. Targeting the feedback loop between tumor cells and fibroblasts is a promising anti-metastatic strategy.

The mechanism of geniposide in patients with COVID-19 and atherosclerosis: A pharmacological and bioinformatics analysis.

In patients with severe acute respiratory syndrome coronavirus 2 (which causes coronavirus disease 2019 [COVID-19]), oxidative stress (OS) is associated with disease severity and death. OS is also involved in the pathogenesis of atherosclerosis (AS). Previous studies have shown that geniposide has anti-inflammatory and anti-viral properties, and can protect cells against OS. However, the potential target(s) of geniposide in patients with COVID-19 and AS, as well as the mechanism it uses, are unclear. We combined pharmacology and bioinformatics analysis to obtain geniposide against COVID-19/AS targets, and build protein-protein interaction network to filter hub genes. The hub genes were performed an enrichment analysis by ClueGO, including Gene Ontology and KEGG. The Enrichr database and the target microRNAs (miRNAs) of hub genes were predicted through the MiRTarBase via Enrichr. The common miRNAs were used to construct the miRNAs-mRNAs regulated network, and the miRNAs' function was evaluated by mirPath v3.0 software. Two hundred forty-seven targets of geniposide were identified in patients with COVID-19/AS comorbidity by observing the overlap between the genes modulated by geniposide, COVID-19, and AS. A protein-protein interaction network of geniposide in patients with COVID-19/AS was constructed, and 27 hub genes were identified. The results of enrichment analysis suggested that geniposide may be involved in regulating the OS via the FoxO signaling pathway. MiRNA-mRNA network revealed that hsa-miR-34a-5p may play an important role in the therapeutic mechanism of geniposide in COVID-19/AS patients. Our study found that geniposide represents a promising therapy for patients with COVID-19 and AS comorbidity. Furthermore, the target genes and miRNAs that we identified may aid the development of new treatment strategies against COVID-19/AS.

Investigation of the relationship between COVID-19 and pancreatic cancer using bioinformatics and systems biology approaches.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, poses a huge threat to human health. Pancreatic cancer (PC) is a malignant tumor with high mortality. Research suggests that infection with SARS-CoV-2 may increase disease severity and risk of death in patients with pancreatic cancer, while pancreatic cancer may also increase the likelihood of contracting SARS-CoV-2, but the link is unclear. METHODS: This study investigated the transcriptional profiles of COVID-19 and PC patients, along with their respective healthy controls, using bioinformatics and systems biology approaches to uncover the molecular mechanisms linking the 2 diseases. Specifically, gene expression data for COVID-19 and PC patients were obtained from the Gene Expression Omnibus datasets, and common differentially expressed genes (DEGs) were identified. Gene ontology and pathway enrichment analyses were performed on the common DEGs to elucidate the regulatory relationships between the diseases. Additionally, hub genes were identified by constructing a protein-protein interaction network from the shared DEGs. Using these hub genes, we conducted regulatory network analyses of microRNA/transcription factors-genes relationships, and predicted potential drugs for treating COVID-19 and PC. RESULTS: A total of 1722 and 2979 DEGs were identified from the transcriptome data of PC (GSE119794) and COVID-19 (GSE196822), respectively. Among these, 236 common DEGs were found between COVID-19 and PC based on protein-protein interaction analysis. Functional enrichment analysis indicated that these shared DEGs were involved in pathways related to viral genome replication and tumorigenesis. Additionally, 10 hub genes, including extra spindle pole bodies like 1, holliday junction recognition protein, marker of proliferation Ki-67, kinesin family member 4A, cyclin-dependent kinase 1, topoisomerase II alpha, cyclin B2, ubiquitin-conjugating enzyme E2 C, aurora kinase B, and targeting protein for Xklp2, were identified. Regulatory network analysis revealed 42 transcription factors and 23 microRNAs as transcriptional regulatory signals. Importantly, lucanthone, etoposide, troglitazone, resveratrol, calcitriol, ciclopirox, dasatinib, enterolactone, methotrexate, and irinotecan emerged as potential therapeutic agents against both COVID-19 and PC. CONCLUSION: This study unveils potential shared pathogenic mechanisms between PC and COVID-19, offering novel insights for future research and therapeutic strategies for the treatment of PC and SARS-CoV-2 infection.

Exploring therapeutic targets and molecular mechanisms for treating diabetes mellitus-associated heart failure with Qishen Yiqi dropping pills: A network pharmacology and bioinformatics approach.

Diabetes mellitus (DM) and heart failure frequently coexist, presenting significant public health challenges. QiShenYiQi Dropping Pills (QSDP) are widely employed in the treatment of diabetes mellitus concomitant with heart failure (DM-HF). Nevertheless, the precise mechanisms underlying their efficacy have yet to be elucidated. Active ingredients and likely targets of QSDP were retrieved from the TCMSP and UniProt databases. Genes associated with DM-HF were pinpointed through searches in the GeneCards, OMIM, DisGeNET, and TTD databases. Differential genes connected to DM-HF were sourced from the GEO database. Enrichment analyses via gene ontology and Kyoto Encyclopedia of Genes and Genomes pathways, as well as immune infiltration assessments, were conducted using R software. Further analysis involved employing molecular docking strategies to explore the interactions between the identified targets and active substances in QSDP that are pertinent to DM-HF treatment. This investigation effectively discerned 108 active compounds and 257 targets relevant to QSDP. A protein-protein interaction network was constructed, highlighting 6 central targets for DM-HF treatment via QSDP. Gene ontology enrichment analysis predominantly linked these targets with responses to hypoxia, metabolism of reactive oxygen species, and cytokine receptor interactions. Analysis of Kyoto Encyclopedia of Genes and Genomes pathways demonstrated that these targets mainly participate in pathways linked to diabetic complications, such as AGE-RAGE signaling, dyslipidemia, arteriosclerosis, the HIF-1 signaling pathway, and the tumor necrosis factor signaling pathway. Further, immune infiltration analysis implied that QSDP's mechanism in treating DM-HF might involve immune-mediated inflammation and crucial signaling pathways. Additionally, molecular docking studies showed that the active substances in QSDP have strong binding affinities with these identified targets. This research presents a new model for addressing DM-HF through the use of QSDP, providing novel insights into incorporating traditional Chinese medicine (TCM) principles in the clinical treatment of DM-HF. The implications of these findings are substantial for both clinical application and further scientific inquiry.

Bioinformatics and machine learning approaches reveal key genes and underlying molecular mechanisms of atherosclerosis: A review.

Atherosclerosis (AS) causes thickening and hardening of the arterial wall due to accumulation of extracellular matrix, cholesterol, and cells. In this study, we used comprehensive bioinformatics tools and machine learning approaches to explore key genes and molecular network mechanisms underlying AS in multiple data sets. Next, we analyzed the correlation between AS and immune fine cell infiltration, and finally performed drug prediction for the disease. We downloaded GSE20129 and GSE90074 datasets from the Gene expression Omnibus database, then employed the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts algorithm to analyze 22 immune cells. To enrich for functional characteristics, the black module correlated most strongly with T cells was screened with weighted gene co-expression networks analysis. Functional enrichment analysis revealed that the genes were mainly enriched in cell adhesion and T-cell-related pathways, as well as NF-κ B signaling. We employed the Lasso regression and random forest algorithms to screen out 5 intersection genes (CCDC106, RASL11A, RIC3, SPON1, and TMEM144). Pathway analysis in gene set variation analysis and gene set enrichment analysis revealed that the key genes were mainly enriched in inflammation, and immunity, among others. The selected key genes were analyzed by single-cell RNA sequencing technology. We also analyzed differential expression between these 5 key genes and those involved in iron death. We found that ferroptosis genes ACSL4, CBS, FTH1 and TFRC were differentially expressed between AS and the control groups, RIC3 and FTH1 were significantly negatively correlated, whereas SPON1 and VDAC3 were significantly positively correlated. Finally, we used the Connectivity Map database for drug prediction. These results provide new insights into AS genetic regulation.

Accurate prediction of protein function using statistics-informed graph networks.

Understanding protein function is pivotal in comprehending the intricate mechanisms that underlie many crucial biological activities, with far-reaching implications in the fields of medicine, biotechnology, and drug development. However, more than 200 million proteins remain uncharacterized, and computational efforts heavily rely on protein structural information to predict annotations of varying quality. Here, we present a method that utilizes statistics-informed graph networks to predict protein functions solely from its sequence. Our method inherently characterizes evolutionary signatures, allowing for a quantitative assessment of the significance of residues that carry out specific functions. PhiGnet not only demonstrates superior performance compared to alternative approaches but also narrows the sequence-function gap, even in the absence of structural information. Our findings indicate that applying deep learning to evolutionary data can highlight functional sites at the residue level, providing valuable support for interpreting both existing properties and new functionalities of proteins in research and biomedicine.