Integrated COVID-19 Predictor: Differential expression analysis to reveal potential biomarkers and prediction of coronavirus using RNA-Seq profile data

A chitosan-based hydrogel with ultrasound-driven immuno-sonodynamic therapeutic effect for accelerated bacterial infected wound healing

Sonodynamic therapy has attracted much attention as a noninvasive treatment for deep infections. However, it is challenging to achieve high antibacterial activity for hydrogels under ultrasonic irradiation due to the relatively weak sono-catalysis capability of sonosensitizers. Herein, an ultrasound-responsive antibacterial hydrogel (Fe3O4/HA/Ber-LA) composed of Fe3O4-grafted-Berberine, chitosan molecules modified with L-arginine and poly (vinyl alcohol) is prepared for enhanced sonodynamic therapy and immunoregulation. The formation of heterojunction between berberine and Fe3O4 with different work function promotes the charge separation and electron flow and disrupts the conjugated structure of berberine, causing a significant decrease in the band gap, eventually enhancing the sonocatalytic activity. The combination of berberine with Fe3O4 also significantly improves the oxygen adsorption energy, enabling more O2 molecules to react with the electron-rich regions on the surface of Fe3O4 to generate more reactive oxygen species (ROS). L-arginine grafted in the hydrogel is catalyzed by the ROS to release nitric oxide, which not only possesses antibacterial activity, but also positively affects macrophage M1 polarization to display potent phagocytosis to Staphylococcus aureus, thus achieving immuno-sonodynamic therapy. Hence, Fe3O4/HA/Ber-LA hydrogel under ultrasound irradiation shows excellent antibacterial activity. Furthermore, the antioxidative activity and anti-inflammatory effect of berberine released from the hybrid hydrogel induces macrophages to polarize towards the anti-inflammatory M2 status as infection comes under control, thus accelerating the wound healing. The hybrid hydrogel based on the immuno-sonodynamic therapy may be an extraordinary candidate for the treatment of deep infections.

Comparative transcriptome analysis of canola carrying a single vs stacked resistance genes against clubroot

Pyramiding resistance genes may expand the efficacy and scope of a canola variety against clubroot (Plasmodiophora brassicae), a serious threat to canola production in western Canada. However, the mechanism(s) of multigenic resistance, especially the potential interaction among clubroot resistance (CR) genes, are not well understood. In this study, transcriptome was compared over three canola (Brassica napus L.) inbred/hybrid lines carrying a single CR gene in chromosome A03 (CRaM, Line 16) or A08 (Crr1rutb, Line 20), and both genes (CRaM+Crr1rutb, Line 15) inoculated with a field population (L-G2) of P. brassicae pathotype X, a new variant found in western Canada recently. The line16 was susceptible, while lines 15 and 20 were partially resistant. Functional annotation identified differential expression of genes (DEGs) involved in biosynthetic processes responsive to stress and regulation of cellular process; The Venn diagram showed that the partially resistant lines 15 and 20 shared 1,896 differentially expressed genes relative to the susceptible line 16, and many of these DEGs are involved in defense responses, activation of innate immunity, hormone biosynthesis and programmed cell death. The transcription of genes involved in Pathogen-Associated Molecular Pattern (PAMP)-Triggered and Effector-Triggered Immunity (PTI and ETI) was particularly up-regulated, and the transcription level was higher in line 15 (CRaM + Crr1rutb) than in line 20 (Crr1rutb only) for most of the DEGs. These results indicated that the partial resistance to the pathotype X was likely conferred by the CR gene Crr1rutb for both lines 15 and 20 that functioned via the activation of both PTI and ETI signaling pathways. Additionally, these two CR genes might have synergistic effects against the pathotype X, based on the higher transcription levels of defense-related DEGs expressed by inoculated line 15, highlighting the benefit of gene stacking for improved canola resistance as opposed to a single CR gene alone.

Integrated analysis of RNA-seq datasets reveals novel targets and regulators of COVID-19 severity

During the COVID-19 pandemic, RNA-seq datasets were produced to investigate the virus-host relationship. However, much of these data remains underexplored. To improve the search for molecular targets and biomarkers, we performed an integrated analysis of multiple RNA-seq datasets, expanding the cohort and including patients from different countries, encompassing severe and mild COVID-19 patients. Our analysis revealed that severe COVID-19 patients exhibit overexpression of genes coding for proteins of extracellular exosomes, endomembrane system, and neutrophil granules (e.g., S100A9, LY96, and RAB1B), which may play an essential role in the cellular response to infection. Concurrently, these patients exhibit down-regulation of genes encoding components of the T cell receptor complex and nucleolus, including TP53, IL2RB, and NCL Finally, SPI1 may emerge as a central transcriptional factor associated with the up-regulated genes, whereas TP53, MYC, and MAX were associated with the down-regulated genes during COVID-19. This study identified targets and transcriptional factors, lighting on the molecular pathophysiology of syndrome coronavirus 2 infection.

SinCWIm: An imputation method for single-cell RNA sequence dropouts using weighted alternating least squares

BACKGROUND AND OBJECTIVES

Single-cell RNA sequencing (scRNA-seq) provides a powerful tool for exploring cellular heterogeneity, discovering novel or rare cell types, distinguishing between tissue-specific cellular composition, and understanding cell differentiation during development. However, due to technological limitations, dropout events in scRNA-seq can mistakenly convert some entries in the real data to zero. This is equivalent to introducing noise into the data of cell gene expression entries. The data is contaminated, which affects the performance of downstream analyses, including clustering, cell annotation, differential gene expression analysis, and so on. Therefore, it is a crucial work to accurately determine which zeros are due to dropout events and perform imputation operations on them.

METHODS

Considering the different confidence levels of different zeros in the gene expression matrix, this paper proposes a SinCWIm method for dropout events in scRNA-seq based on weighted alternating least squares (WALS). The method utilizes Pearson correlation coefficient and hierarchical clustering to quantify the confidence of zero entries. It is then combined with WALS for matrix decomposition. And the imputation result is made close to the actual number by outlier removal and data correction operations.

RESULTS

A total of eight single-cell sequencing datasets were used for comparative experiments to demonstrate the overall superiority of SinCWIm over state-of-the-art models. SinCWIm was applied to cluster the data to obtain an adjusted RAND index evaluation, and the Usoskin, Pollen and Bladder datasets scored 94.46%, 96.48% and 76.74%, respectively. In addition, significant improvements were made in the retention of differential expression genes and visualization.

CONCLUSIONS

SinCWIm provides a valuable imputation method for handling dropout events in single-cell sequencing data. In comparison to advanced methods, SinCWIm demonstrates excellent performance in clustering, visualization and other aspects. It is applicable to various single-cell sequencing datasets.

Identifying OGN as a Biomarker Covering Multiple Pathogenic Pathways for Diagnosing Heart Failure: From Machine Learning to Mechanism Interpretation

BACKGROUND

The pathophysiologic heterogeneity of heart failure (HF) necessitates a more detailed identification of diagnostic biomarkers that can reflect its diverse pathogenic pathways.

METHODS

We conducted weighted gene and multiscale embedded gene co-expression network analysis on differentially expressed genes obtained from HF and non-HF specimens. We employed a machine learning integration framework and protein-protein interaction network to identify diagnostic biomarkers. Additionally, we integrated gene set variation analysis, gene set enrichment analysis (GSEA), and transcription factor (TF)-target analysis to unravel the biomarker-dominant pathways. Leveraging single-sample GSEA and molecular docking, we predicted immune cells and therapeutic drugs related to biomarkers. Quantitative polymerase chain reaction validated the expressions of biomarkers in the plasma of HF patients. A two-sample Mendelian randomization analysis was implemented to investigate the causal impact of biomarkers on HF.

RESULTS

We first identified COL14A1, OGN, MFAP4, and SFRP4 as candidate biomarkers with robust diagnostic performance. We revealed that regulating biomarkers in HF pathogenesis involves TFs (BNC2, MEOX2) and pathways (cell adhesion molecules, chemokine signaling pathway, cytokine-cytokine receptor interaction, oxidative phosphorylation). Moreover, we observed the elevated infiltration of effector memory CD4+ T cells in HF, which was highly related to biomarkers and could impact immune pathways. Captopril, aldosterone antagonist, cyclopenthiazide, estradiol, tolazoline, and genistein were predicted as therapeutic drugs alleviating HF via interactions with biomarkers. In vitro study confirmed the up-regulation of OGN as a plasma biomarker of HF. Mendelian randomization analysis suggested that genetic predisposition toward higher plasma OGN promoted the risk of HF.

CONCLUSIONS

We propose OGN as a diagnostic biomarker for HF, which may advance our understanding of the diagnosis and pathogenesis of HF.

Undifferentially Expressed CXXC5 as a Transcriptionally Regulatory Biomarker of Breast Cancer

This work hypothesizes that some genes undergo radically changed transcription regulations (TRs) in breast cancer (BC), but don't show differential expressions for unknown reasons. The TR of a gene is quantitatively formulated by a regression model between the expression of this gene and multiple transcription factors (TFs). The difference between the predicted and real expression levels of a gene in a query sample is defined as the mqTrans value of this gene, which quantitatively reflects its regulatory changes. This work systematically screens the undifferentially expressed genes with differentially expressed mqTrans values in 1036 samples across five datasets and three ethnic groups. This study calls the 25 genes satisfying the above hypothesis in at least four datasets as dark biomarkers, and the strong dark biomarker gene CXXC5 (CXXC Finger Protein 5) is even supported by all the five independent BC datasets. Although CXXC5 does not show differential expressions in BC, its transcription regulations show quantitative associations with BCs in diversified cohorts. The overlapping long noncoding RNAs (lncRNAs) may have contributed their transcripts to the expression miscalculations of dark biomarkers. The mqTrans analysis serves as a complementary view of the transcriptome-based detections of biomarkers that are ignored by many existing studies.

Monkeypox diagnosis using ensemble classification

The world has recently been exposed to a fierce attack from many viral diseases, such as Covid-19, that exhausted medical systems around the world. Such attack had a negative impact not only on the health status of people or the high death rate, but also had a bad impact on the economic situation, which affected all countries of the world especially the poor and the developing ones. Monkeypox is one of the latest viral diseases that may cause a pandemic in the near future if not dealt and diagnosed with appropriately. This paper provides a new strategy for diagnosing monkeypox, which is called; Accurate Monkeypox Diagnosing Strategy (AMDS). The proposed AMDS consists of two phases, which are; (i) pre-processing and (ii) classification. During the pre-processing phase, the most effective feature are selected using Binary Tiki-Taka Algorithm (BTTA). On the other hand, in the classification phase, ensemble classification is used for diagnosing new cases, which combines evidence from three different new classifiers, namely; (a) Layered K-Nearest Neighbors (LKNN), (b) Statistical Naïve Bayes (SNB), and (c) Deep Learning Classifier (DLC). Moreover, the decisions of the proposed classifiers are merged in a new voting scheme called Fuzzified Voting Scheme (FVS). AMDS has been compared against recent diagnostic strategies. Experimental results have proven that AMDS outperforms other monkeypox diagnostic strategies as it introduces the most accurate diagnosis according to two different datasets.

Computational identification of differentially-expressed genes as suggested novel COVID-19 biomarkers: A bioinformatics analysis of expression profiles

COVID-19 was declared a pandemic in March 2020, and since then, it has not stopped spreading like wildfire in almost every corner of the world, despite the many efforts made to stem its spread. SARS-CoV-2 has one of the biggest genomes among RNA viruses and presents unique characteristics that differentiate it from other coronaviruses, making it even more challenging to find a cure or vaccine that is efficient enough. This work aims, using RNA sequencing (RNA-Seq) data, to evaluate whether the expression of specific human genes in the host can vary in different grades of disease severity and to determine the molecular origins of the differences in response to SARS-CoV-2 infection in different patients. In addition to quantifying gene expression, data coming from RNA-Seq allow for the discovery of new transcripts, the identification of alternative splicing events, the detection of allele-specific expression, and the detection of post-transcriptional alterations. For this reason, we performed differential expression analysis on different expression profiles of COVID-19 patients, using RNA-Seq data coming from NCBI public repository, and we obtained the lists of all differentially expressed genes (DEGs) emerging from 7 experimental conditions. We performed a Gene Set Enrichment Analysis (GSEA) on these genes to find possible correlations between DEGs and known disease phenotypes. We mainly focused on DEGs coming out from the analysis of the contrasts involving severe conditions to infer any possible relation between a worsening of the clinical picture and an over-representation of specific genes. Based on the obtained results, this study indicates a small group of genes that result up-regulated in the severe form of the disease. EXOSC5, MESD, REXO2, and TRMT2A genes are not differentially expressed or not present in the other conditions, being for that reason, good biomarkers candidates for the severe form of COVID-19 disease. The use of specific over-expressed genes, whether up-regulated or down-regulated, which have an individual role in each different condition of COVID-19 as a biomarker, can assist in early diagnosis.

A systematic review of artificial intelligence-based COVID-19 modeling on multimodal genetic information

This study systematically reviews the Artificial Intelligence (AI) methods developed to resolve the critical process of COVID-19 gene data analysis, including diagnosis, prognosis, biomarker discovery, drug responsiveness, and vaccine efficacy. This systematic review follows the guidelines of Preferred Reporting for Systematic Reviews and Meta-Analyses (PRISMA). We searched PubMed, Embase, Web of Science, and Scopus databases to identify the relevant articles from January 2020 to June 2022. It includes the published studies of AI-based COVID-19 gene modeling extracted through relevant keyword searches in academic databases. This study included 48 articles discussing AI-based genetic studies for several objectives. Ten articles confer about the COVID-19 gene modeling with computational tools, and five articles evaluated ML-based diagnosis with observed accuracy of 97% on SARS-CoV-2 classification. Gene-based prognosis study reviewed three articles and found host biomarkers detecting COVID-19 progression with 90% accuracy. Twelve manuscripts reviewed the prediction models with various genome analysis studies, nine articles examined the gene-based in silico drug discovery, and another nine investigated the AI-based vaccine development models. This study compiled the novel coronavirus gene biomarkers and targeted drugs identified through ML approaches from published clinical studies. This review provided sufficient evidence to delineate the potential of AI in analyzing complex gene information for COVID-19 modeling on multiple aspects like diagnosis, drug discovery, and disease dynamics. AI models entrenched a substantial positive impact by enhancing the efficiency of the healthcare system during the COVID-19 pandemic.

Smart systems and data-driven services in healthcare

The modern development of Medicine and Healthcare is primarily based on the automation of various processes to support the correct and timely medical decisions. A doctor or other medical staff's fast, accurate, reliable diagnosis, prevention, and treatment processes are the keys to quality patient care. This aim can be achieved by developing reliable, intelligent Healthcare and medicine service systems for various purposes. The Special Issue covers science-intensive real-time solutions that underlie smart systems and data-driven services in Healthcare and Medicine.

A Novel 3-Gene Signature for Identifying COVID-19 Patients Based on Bioinformatics and Machine Learning

Although many biomarkers associated with coronavirus disease 2019 (COVID-19) were found, a novel signature relevant to immune cells has not been developed. In this work, the “CIBERSORT” algorithm was used to assess the fraction of immune infiltrating cells in GSE152641 and GSE171110. Key modules associated with important immune cells were selected by the “WGCNA” package. The “GO” enrichment analysis was used to reveal the biological function associated with COVID-19. The “Boruta” algorithm was used to screen candidate genes, and the “LASSO” algorithm was used for collinearity reduction. A novel gene signature was developed based on multivariate logistic regression analysis. Subsequently, M0 macrophages (PRAUC = 0.948 in GSE152641 and PRAUC = 0.981 in GSE171110) and neutrophils (PRAUC = 0.892 in GSE152641 and PRAUC = 0.960 in GSE171110) were considered as important immune cells. Forty-three intersected genes from two modules were selected, which mainly participated in some immune-related activities. Finally, a three-gene signature comprising CLEC4D, DUSP13, and UNC5A that can accurately distinguish COVID-19 patients and healthy controls in three datasets was constructed. The ROCAUC was 0.974 in the training set, 0.946 in the internal test set, and 0.709 in the external test set. In conclusion, we constructed a three-gene signature to identify COVID-19, and CLEC4D, DUSP13, and UNC5A may be potential biomarkers for COVID-19 patients.