Search

Abstract Background Flux Balance Analysis (FBA) is a key metabolic modeling method used to simulate cellular metabolism under steady-state conditions. Its simplicity and versatility have led to various strategies incorporating transcriptomic and proteomic data into FBA, successfully predicting flux distribution and phenotypic results. However, despite these advances, the untapped potential lies in leveraging gene-related connections like co-expression patterns for valuable insights. Results To fill this gap, we introduce ICON-GEMs, an innovative constraint-based model to incorporate gene co-expression network into the FBA model, facilitating more precise determination of flux distributions and functional pathways. In this study, transcriptomic data from both Escherichia coli and Saccharomyces cerevisiae were integrated into their respective genome-scale metabolic models. A comprehensive gene co-expression network was constructed as a global view of metabolic mechanism of the cell. By leveraging quadratic programming, we maximized the alignment between pairs of reaction fluxes and the correlation of their corresponding genes in the co-expression network. The outcomes notably demonstrated that ICON-GEMs outperformed existing methodologies in predictive accuracy. Flux variabilities over subsystems and functional modules also demonstrate promising results. Furthermore, a comparison involving different types of biological networks, including protein–protein interactions and random networks, reveals insights into the utilization of the co-expression network in genome-scale metabolic engineering. Conclusion ICON-GEMs introduce an innovative constrained model capable of simultaneous integration of gene co-expression networks, ready for board application across diverse transcriptomic data sets and multiple organisms. It is freely available as open-source at https://github.com/ThummaratPaklao/ICOM-GEMs.git .

Breast cancer has the highest diagnosis rate among all cancers. Tumor budding (TB) is recognized as a recent prognostic marker. Identifying genes specific to high-TB samples is crucial for hindering tumor progression and metastasis. In this study, we utilized an RNA sequencing technique, called TempO-Seq, to profile transcriptomic data from breast cancer samples, aiming to identify biomarkers for high-TB cases. Through differential expression analysis and mutual information, we identified seven genes (NOL4, STAR, C8G, NEIL1, SLC46A3, FRMD6, and SCARF2) that are potential biomarkers in breast cancer. To gain more relevant proteins, further investigation based on a protein–protein interaction network and the network diffusion technique revealed enrichment in the Hippo signaling and Wnt signaling pathways, promoting tumor initiation, invasion, and metastasis in several cancer types. In conclusion, these novel genes, recognized as overexpressed in high-TB samples, along with their associated pathways, offer promising therapeutic targets, thus advancing treatment and diagnosis for breast cancer.

Background Identifying the genes responsible for diseases requires precise prioritization of significant genes. Gene expression analysis enables differentiation between gene expressions in disease and normal samples. Increasing the number of high-quality samples enhances the strength of evidence regarding gene involvement in diseases. This process has led to the discovery of disease biomarkers through the collection of diverse gene expression data. Methods This study presents GeneCompete, a web-based tool that integrates gene expression data from multiple platforms and experiments to identify the most promising biomarkers. GeneCompete incorporates a novel union strategy and eight well-established ranking methods, including Win-Loss, Massey, Colley, Keener, Elo, Markov, PageRank, and Bi-directional PageRank algorithms, to prioritize genes across multiple gene expression datasets. Each gene in the competition is assigned a score based on log-fold change values, and significant genes are determined as winners. Results We tested the tool on the expression datasets of Hypertrophic cardiomyopathy (HCM) and the datasets from Microarray Quality Control (MAQC) project, which include both microarray and RNA-Sequencing techniques. The results demonstrate that all ranking scores have more power to predict new occurrence datasets than the classical method. Moreover, the PageRank method with a union strategy delivers the best performance for both up-regulated and down-regulated genes. Furthermore, the top-ranking genes exhibit a strong association with the disease. For MAQC, the two-sides ranking score shows a high relationship with TaqMan validation set in all log-fold change thresholds. Conclusion GeneCompete is a powerful web-based tool that revolutionizes the identification of disease-causing genes through the integration of gene expression data from multiple platforms and experiments.

Background: Despite advances in autism spectrum disorder (ASD) research and the vast genomic, transcriptomic, and proteomic data available, there are still controversies regarding the pathways and molecular signatures underlying the neurodevelopmental disorders leading to ASD. Purpose: To delineate these underpinning signatures, we examined the two largest gene expression meta-analysis datasets obtained from the brain and peripheral blood mononuclear cells (PBMCs) of 1355 ASD patients and 1110 controls. Methods: We performed network, enrichment, and annotation analyses using the differentially expressed genes, transcripts, and proteins identified in ASD patients. Results: Transcription factor network analyses in up- and down-regulated genes in brain tissue and PBMCs in ASD showed eight main transcription factors, namely: BCL3, CEBPB, IRF1, IRF8, KAT2A, NELFE, RELA, and TRIM28. The upregulated gene networks in PBMCs of ASD patients are strongly associated with activated immune-inflammatory pathways, including interferon-α signaling, and cellular responses to DNA repair. Enrichment analyses of the upregulated CNS gene networks indicate involvement of immune-inflammatory pathways, cytokine production, Toll-Like Receptor signalling, with a major involvement of the PI3K-Akt pathway. Analyses of the downregulated CNS genes suggest electron transport chain dysfunctions at multiple levels. Network topological analyses revealed that the consequent aberrations in axonogenesis, neurogenesis, synaptic transmission, and regulation of transsynaptic signalling affect neurodevelopment with subsequent impairments in social behaviours and neurocognition. The results suggest a defense response against viral infection. Conclusions: Peripheral activation of immune-inflammatory pathways, most likely induced by viral infections, may result in CNS neuroinflammation and mitochondrial dysfunction, leading to abnormalities in transsynaptic transmission, and brain neurodevelopment.

Topological data analysis (TDA) methods have recently emerged as powerful tools for uncovering intricate patterns and relationships in complex biological data, demonstrating their effectiveness in identifying key genes in breast, lung, and blood cancer. In this study, we applied a TDA technique, specifically persistent homology (PH), to identify key pathways for early detection of hepatocellular carcinoma (HCC). Recognizing the limitations of current strategies for this purpose, we meticulously used PH to analyze RNA sequencing (RNA-seq) data from peripheral blood of both HCC patients and normal controls. This approach enabled us to gain nuanced insights by detecting significant differences between control and disease sample classes. By leveraging topological descriptors crucial for capturing subtle changes between these classes, our study identified 23 noteworthy pathways, including the apelin signaling pathway, the IL-17 signaling pathway, and the p53 signaling pathway. Subsequently, we performed a comparative analysis with a classical enrichment-based pathway analysis method which revealed both shared and unique findings. Notably, while the IL-17 signaling pathway was identified by both methods, the HCC-related apelin signaling and p53 signaling pathways emerged exclusively through our topological approach. In summary, our study underscores the potential of PH to complement traditional pathway analysis approaches, potentially providing additional knowledge for the development of innovative early detection strategies of HCC from blood samples.

Hyperspectral imaging has emerged as a powerful tool for the non-destructive assessment of plant properties, including the quantification of phytochemical contents. Traditional methods for antioxidant analysis in holy basil (Ocimum tenuiflorum L.) are time-consuming, while hyperspectral imaging has the potential to rapidly observe holy basil. In this study, we employed hyperspectral imaging combined with machine learning techniques to determine the levels of total phenolic contents in Thai holy basil. Spectral data were acquired from 26 holy basil cultivars at different growth stages, and the total phenolic contents of the samples were measured. To extract the characteristics of the spectral data, we used 22 statistical features in both time and frequency domains. Relevant features were selected and combined with the corresponding total phenolic content values to develop a neural network model for classifying the phenolic content levels into 'low' and 'normal-to-high' categories. The neural network model demonstrated high performance, achieving an area under the receiver operating characteristic curve of 0.8113, highlighting its effectiveness in predicting phenolic content levels based on the spectral data. Comparative analysis with other machine learning techniques confirmed the superior performance of the neural network approach. Further investigation revealed that the model exhibited increased confidence in predicting the phenolic content levels of older holy basil samples. This study exhibits the potential of integrating hyperspectral imaging, feature extraction, and machine learning techniques for the rapid and non-destructive assessment of phenolic content levels in holy basil. The demonstrated effectiveness of this approach opens new possibilities for screening antioxidant properties in plants, facilitating efficient decision-making processes for researchers based on comprehensive spectral data. [ABSTRACT FROM AUTHOR]

Holy basil (Ocimum tenuiflorum L.) is a widely renowned herb for its abundance of bioactive compounds and medicinal applications. Nevertheless, there exists a dearth of knowledge regarding the variability among holy basil cultivars capable of yielding substantial bioactive compounds. This study aims to address this gap by shedding light on the diversity of antioxidant capacities within different accessions of Thai holy basil by employing a density-based clustering algorithm to categorize the holy basil cultivars that demonstrate notable antioxidant potential. The study involves the analysis of the anthocyanin, flavonoid, phenolic, and terpenoid content, as well as DPPH antioxidant activity, in 26 Thai holy basil accessions collected from diverse locations in Thailand. Among the 26 tested Thai holy basil cultivars, terpenoids were found to be the dominant class of compounds, with average values of 707 mg/gDW, while the levels of flavonoids and phenolic compounds remained below 65 mg rutin/gDW and 46 mg GAE/gDW, respectively. The DPPH assay in holy basil cultivars demonstrated that the antioxidant activity ranged between 50% and 93%. After standardizing the data, the clustering results revealed four distinct groups of cultivars: the first group, with low antioxidant levels; the second group, with high terpenoid content; the third group, with high flavonoid, DPPH antioxidant activity, and phenolic content; and the fourth group, with elevated levels of anthocyanin, DPPH antioxidant activity, and phenolic content. A strong positive correlation was observed among DPPH antioxidant activity, flavonoids, and phenolics. Specific cultivars: The Red, OC108, and OC106 holy basil cultivars in cluster 4 exhibited high anthocyanin and phenolic production. In cluster 3, the accessions OC113, OC057, OC063, and OC059 showed high DPPH antioxidant activity, flavonoids, and phenolics, while, in cluster 2, only accessions from Udon Thani, Thailand—namely OC194 and OC195—displayed high terpenoid content. Ultimately, this study significantly contributes to the inherent diversity in the antioxidant capacities among various Thai holy basil cultivars. It lays the foundation for targeted breeding strategies and informed choices regarding consumption. The comprehensive insights from this analysis hold the potential to accurately identify holy basil cultivars with promising applications in medicine, functional foods, and the nutraceutical industry.

Identification of drug–target interaction (DTI) is a crucial step to reduce time and cost in the drug discovery and development process. Since various biological data are publicly available, DTIs have been identified computationally. To predict DTIs, most existing methods focus on a single similarity measure of drugs and target proteins, whereas some recent methods integrate a particular set of drug and target similarity measures by a single integration function. Therefore, many DTIs are still missing. In this study, we propose heterogeneous network propagation with the forward similarity integration (FSI) algorithm, which systematically selects the optimal integration of multiple similarity measures of drugs and target proteins. Seven drug–drug and nine target–target similarity measures are applied with four distinct integration methods to finally create an optimal heterogeneous network model. Consequently, the optimal model uses the target similarity based on protein sequences and the fused drug similarity, which combines the similarity measures based on chemical structures, the Jaccard scores of drug–disease associations, and the cosine scores of drug–drug interactions. With an accuracy of 99.8%, this model significantly outperforms others that utilize different similarity measures of drugs and target proteins. In addition, the validation of the DTI predictions of this model demonstrates the ability of our method to discover missing potential DTIs.

Identifying genes associated with autism spectrum disorder (ASD) is crucial for understanding the underlying mechanisms of the disorder. However, ASD is a complex condition involving multiple mechanisms, and this has resulted in an unclear understanding of the disease and a lack of precise knowledge concerning the genes associated with ASD. To address these challenges, we conducted a systematic analysis that integrated multiple data sources, including associations among ASD-associated genes and gene expression data from ASD studies. With these data, we generated both a gene embedding profile that captured the complex relationships between genes and a differential gene expression profile (built from the gene expression data). We utilized the XGBoost classifier and leveraged these profiles to identify novel ASD associations. This approach revealed 10,848 potential gene–gene associations and inferred 125 candidate genes, with DNA Topoisomerase I, ATP Synthase F1 Subunit Gamma, and Neuronal Calcium Sensor 1 being the top three candidates. We conducted a statistical analysis to assess the relevance of candidate genes to specific functions and pathways. Additionally, we identified sub-networks within the candidate network to uncover sub-groups of associations that could facilitate the identification of potential ASD-related genes. Overall, our systematic analysis, which integrated multiple data sources, represents a significant step towards unraveling the complexities of ASD. By combining network-based gene associations, gene expression data, and machine learning, we contribute to ASD research and facilitate the discovery of new targets for molecularly targeted therapies.

Disease-related gene prioritization is one of the most well-established pharmaceutical techniques used to identify genes that are important to a biological process relevant to a disease. In identifying these essential genes, the network diffusion (ND) approach is a widely used technique applied in gene prioritization. However, there is still a large number of candidate genes that need to be evaluated experimentally. Therefore, it would be of great value to develop a new strategy to improve the precision of the prioritization. Given the efficiency and simplicity of centrality measures in capturing a gene that might be important to the network structure, herein, we propose a technique that extends the scope of ND through a centrality measure to identify new disease-related genes. Five common centrality measures with different aspects were examined for integration in the traditional ND model. A total of 40 diseases were used to test our developed approach and to find new genes that might be related to a disease. Results indicated that the best measure to combine with the diffusion is closeness centrality. The novel candidate genes identified by the model for all 40 diseases were provided along with supporting evidence. In conclusion, the integration of network centrality in ND is a simple but effective technique to discover more precise disease-related genes, which is extremely useful for biomedical science.

COVID-19 is a respiratory disease that can spread rapidly. Controlling the spread through vaccination is one of the measures for activating immunization that helps to reduce the number of infected people. Different types of vaccines are effective in preventing and alleviating the symptoms of the disease in different ways. In this study, a mathematical model, SVIHR, was developed to assess the behavior of disease transmission in Thailand by considering the vaccine efficacy of different vaccine types and the vaccination rate. The equilibrium points were investigated and the basic reproduction number R0 was calculated using a next-generation matrix to determine the stability of the equilibrium. We found that the disease-free equilibrium point was asymptotically stable if, and only if, R0<1, and the endemic equilibrium was asymptotically stable if, and only if, R0>1. The simulation results and the estimation of the parameters applied to the actual data in Thailand are reported. The sensitivity of parameters related to the basic reproduction number was compared with estimates of the effectiveness of pandemic controls. The simulations of different vaccine efficacies for different vaccine types were compared and the average mixing of vaccine types was reported to assess the vaccination policies. Finally, the trade-off between the vaccine efficacy and the vaccination rate was investigated, resulting in the essentiality of vaccine efficacy to restrict the spread of COVID-19.

Breast cancer has the highest diagnosis rate among all cancers. Tumor budding (TB) is recognized as a recent prognostic marker. Identifying genes specific to high-TB samples is crucial for hindering tumor progression and metastasis. In this study, we utilized an RNA sequencing technique, called TempO-Seq, to profile transcriptomic data from breast cancer samples, aiming to identify biomarkers for high-TB cases. Through differential expression analysis and mutual information, we identified seven genes (NOL4, STAR, C8G, NEIL1, SLC46A3, FRMD6, and SCARF2) that are potential biomarkers in breast cancer. To gain more relevant proteins, further investigation based on a protein–protein interaction network and the network diffusion technique revealed enrichment in the Hippo signaling and Wnt signaling pathways, promoting tumor initiation, invasion, and metastasis in several cancer types. In conclusion, these novel genes, recognized as overexpressed in high-TB samples, along with their associated pathways, offer promising therapeutic targets, thus advancing treatment and diagnosis for breast cancer. [ABSTRACT FROM AUTHOR]

Drug repositioning, discovering new indications for existing drugs, is known to solve the bottleneck of drug discovery and development. To support a task of drug repositioning, many in silico methods have been proposed for predicting drug-disease associations. A meta-path based approach, which extracts network-based information through paths from a drug to a disease, can produce comparable performance with less required information when compared to other approaches. However, existing meta-path based methods typically use counts of extracted paths and discard information of intermediate nodes in those paths although they are very important indicators, such as drug- and disease-associated proteins. Herein, we propose an ensemble learning method with Meta-path based Gene ontology Profiles for predicting Drug-Disease Associations (MGP-DDA). We exploit gene ontology (GO) terms to link drugs and diseases to their associated functions and act as intermediate nodes in a drug-GO-disease tripartite network. For each drug-disease pair, MGP-DDA utilizes meta-paths to generate novel profiles of GO functions, termed as meta-path based GO profiles. We train bagging and boosting classifiers with those novel features to recognize known (positive) from unknown (unlabeled) drug-disease associations. Consequently, MGP-DDA outperforms the state-of-the-art methods and yields the precision of 88.6%. By MGP-DDA, the eminent number of new drug-disease associations with supporting evidence in ClinicalTrials.gov (37.7%) ensures the practicality of our method in drug repositioning.

At present, a large number of people worldwide have been infected by coronavirus 2019 (COVID-19). When the outbreak of the COVID-19 pandemic begins in a country, its impact is disastrous to both the country and its neighbors. In early 2020, the spread of COVID-19 was associated with global aviation. More recently, COVID-19 infections due to illegal or undocumented immigration have played a significant role in spreading the disease in Southeast Asia countries. Therefore, the spread of COVID-19 of all countries' border should be curbed. Many countries closed their borders to all nations, causing an unprecedented decline in global travel, especially cross-border travel. This restriction affects social and economic trade-offs. Therefore, immigration policies are essential to control the COVID-19 pandemic. To understand and simulate the spread of the disease under different immigration conditions, we developed a novel mathematical model called the Legal immigration and Undocumented immigration from natural borders for Susceptible-Infected-Hospitalized and Recovered people (LUSIHR). The purpose of the model was to simulate the number of infected people under various policies, including uncontrolled, fully controlled, and partially controlled countries. The infection rate was parameterized using the collected data from the Department of Disease Control, Ministry of Public Health, Thailand. We demonstrated that the model possesses nonnegative solutions for favorable initial conditions. The analysis of numerical experiments showed that we could control the virus spread and maintain the number of infected people by increasing the control rate of undocumented immigration across the unprotected natural borders. Next, the obtained parameters were used to visualize the effect of the control rate on immigration at the natural border. Overall, the model was well-suited to explaining and building the simulation. The parameters were used to simulate the trends in the number of people infected from COVID-19.

Salinity is an important environmental factor causing a negative effect on rice production. To prevent salinity effects on rice yields, genetic diversity concerning salt tolerance must be evaluated. In this study, we investigated the salinity responses of rice (Oryza sativa) to determine the critical genes. The transcriptomes of ‘Luang Pratahn’ rice, a local Thai rice variety with high salt tolerance, were used as a model for analyzing and identifying the key genes responsible for salt-stress tolerance. Based on 3' Tag-Seq data from the time course of salt-stress treatment, weighted gene co-expression network analysis was used to identify key genes in gene modules. We obtained 1,386 significantly differentially expressed genes in eight modules. Among them, six modules indicated a significant correlation within 6, 12, or 48h after salt stress. Functional and pathway enrichment analysis was performed on the co-expressed genes of interesting modules to reveal which genes were mainly enriched within important functions for salt-stress responses. To identify the key genes in salt-stress responses, we considered the two-state co-expression networks, normal growth conditions, and salt stress to investigate which genes were less important in a normal situation but gained more impact under stress. We identified key genes for the response to biotic and abiotic stimuli and tolerance to salt stress. Thus, these novel genes may play important roles in salinity tolerance and serve as potential biomarkers to improve salt tolerance cultivars.

Malaria caused by Plasmodium vivax can lead to severe morbidity and death. In addition, resistance has been reported to existing drugs in treating this malaria. Therefore, the identification of new human proteins associated with malaria is urgently needed for the development of additional drugs. In this study, we established an analysis framework to predict human- P. vivax protein associations using network topological profiles from a heterogeneous network structure of human and P. vivax , machine-learning techniques and statistical analysis. Novel associations were predicted and ranked to determine the importance of human proteins associated with malaria. With the best-ranking score, 411 human proteins were identified as promising proteins. Their regulations and functions were statistically analyzed, which led to the identification of proteins involved in the regulation of membrane and vesicle formation, and proteasome complexes as potential targets for the treatment of P. vivax malaria. In conclusion, by integrating related data, our analysis was efficient in identifying potential targets providing an insight into human-parasite protein associations. Furthermore, generalizing this model could allow researchers to gain further insights into other diseases and enhance the field of biomedical science.

Coronavirus disease 2019 (COVID-19) is still an active global public health issue. Although vaccines and therapeutic options are available, some patients experience severe conditions and need critical care support. Hence, identifying key genes or proteins involved in immune-related severe COVID-19 is necessary to find or develop the targeted therapies. This study proposed a novel construction of an immune-related protein interaction network (IPIN) in severe cases with the use of a network diffusion technique on a human interactome network and transcriptomic data. Enrichment analysis revealed that the IPIN was mainly associated with antiviral, innate immune, apoptosis, cell division, and cell cycle regulation signaling pathways. Twenty-three proteins were identified as key proteins to find associated drugs. Finally, poly (I:C), mitomycin C, decitabine, gemcitabine, hydroxyurea, tamoxifen, and curcumin were the potential drugs interacting with the key proteins to heal severe COVID-19. In conclusion, IPIN can be a good representative network for the immune system that integrates the protein interaction network and transcriptomic data. Thus, the key proteins and target drugs in IPIN help to find a new treatment with the use of existing drugs to treat the disease apart from vaccination and conventional antiviral therapy.

Major depressive disorder and major depressive episodes (MDD/MDE) are characterized by the activation of the immune–inflammatory response system (IRS) and the compensatory immune–regulatory system (CIRS). Cannabidiol (CBD) is a phytocannabinoid isolated from the cannabis plant, which is reported to have antidepressant-like and anti-inflammatory effects. The aim of the present study is to examine the effects of CBD on IRS, CIRS, M1, T helper (Th)-1, Th-2, Th-17, T regulatory (Treg) profiles, and growth factors in depression and healthy controls. Culture supernatant of stimulated (5 μg/mL of PHA and 25 μg/mL of LPS) whole blood of 30 depressed patients and 20 controls was assayed for cytokines using the LUMINEX assay. The effects of three CBD concentrations (0.1 µg/mL, 1 µg/mL, and 10 µg/mL) were examined. Depression was characterized by significantly increased PHA + LPS-stimulated Th-1, Th-2, Th-17, Treg, IRS, CIRS, and neurotoxicity profiles. CBD 0.1 µg/mL did not have any immune effects. CBD 1.0 µg/mL decreased CIRS activities but increased growth factor production, while CBD 10.0 µg/mL suppressed Th-1, Th-17, IRS, CIRS, and a neurotoxicity profile and enhanced T cell growth and growth factor production. CBD 1.0 to 10.0 µg/mL dose-dependently decreased sIL-1RA, IL-8, IL-9, IL-10, IL-13, CCL11, G-CSF, IFN-γ, CCL2, CCL4, and CCL5, and increased IL-1β, IL-4, IL-15, IL-17, GM-CSF, TNF-α, FGF, and VEGF. In summary, in this experiment, there was no beneficial effect of CBD on the activated immune profile of depression and higher CBD concentrations can worsen inflammatory processes.

This study used established biomarkers of death from ischemic stroke (IS) versus stroke survival to perform network, enrichment, and annotation analyses. Protein-protein interaction (PPI) network analysis revealed that the backbone of the highly connective network of IS death consisted of IL6, ALB, TNF, SERPINE1, VWF, VCAM1, TGFB1, and SELE. Cluster analysis revealed immune and hemostasis subnetworks, which were strongly interconnected through the major switches ALB and VWF. Enrichment analysis revealed that the PPI immune subnetwork of death due to IS was highly associated with TLR2/4, TNF, JAK-STAT, NOD, IL10, IL13, IL4, and TGF-β1/SMAD pathways. The top biological and molecular functions and pathways enriched in the hemostasis network of death due to IS were platelet degranulation and activation, the intrinsic pathway of fibrin clot formation, the urokinase-type plasminogen activator pathway, post-translational protein phosphorylation, integrin cell-surface interactions, and the proteoglycan-integrin extracellular matrix complex (ECM). Regulation Explorer analysis of transcriptional factors shows: (a) that NFKB1, RELA and SP1 were the major regulating actors of the PPI network; and (b) hsa-mir-26-5p and hsa-16-5p were the major regulating microRNA actors. In conclusion, prevention of death due to IS should consider that current IS treatments may be improved by targeting VWF, the proteoglycan-integrin-ECM complex, TGF-β1/SMAD, NF-κB/RELA and SP1.

There is evidence that schizophrenia is characterized by activation of the immune-inflammatory response (IRS) and compensatory immune-regulatory systems (CIRS) and lowered neuroprotection. Studies performed on antipsychotic-naïve first episode psychosis (AN-FEP) and schizophrenia (FES) patients are important as they may disclose the pathogenesis of FES. However, the protein–protein interaction (PPI) network of FEP/FES is not established. The aim of the current study was to delineate a) the characteristics of the PPI network of AN-FEP and its transition to FES; and b) the biological functions, pathways, and molecular patterns, which are over-represented in FEP/FES. Toward this end, we used PPI network, enrichment, and annotation analyses. FEP and FEP/FES are strongly associated with a response to a bacterium, alterations in Toll-Like Receptor-4 and nuclear factor-κB signaling, and the Janus kinases/signal transducer and activator of the transcription proteins pathway. Specific molecular complexes of the peripheral immune response are associated with microglial activation, neuroinflammation, and gliogenesis. FEP/FES is accompanied by lowered protection against inflammation, in part attributable to dysfunctional miRNA maturation, deficits in neurotrophin and Wnt/catenin signaling, and adherens junction organization. Multiple interactions between reduced brain derived neurotrophic factor, E-cadherin, and β-catenin and disrupted schizophrenia-1 (DISC1) expression increase the vulnerability to the neurotoxic effects of immune molecules, including cytokines and complement factors. In summary: FEP and FES are systemic neuro-immune disorders that are probably triggered by a bacterial stimulus which induces neuro-immune toxicity cascades that are overexpressed in people with reduced anti-inflammatory and miRNA protections, cell–cell junction organization, and neurotrophin and Wnt/catenin signaling.

Functional annotation of unknown function genes reveals unidentified functions that can enhance our understanding of complex genome communications. A common approach for inferring gene function involves the ortholog-based method. However, genetic data alone are often not enough to provide information for function annotation. Thus, integrating other sources of data can potentially increase the possibility of retrieving annotations. Network-based methods are efficient techniques for exploring interactions among genes and can be used for functional inference. In this study, we present an analysis framework for inferring the functions of Plasmodium falciparum genes based on connection profiles in a heterogeneous network between human and Plasmodium falciparum proteins. These profiles were fed into a hybrid deep learning algorithm to predict the orthologs of unknown function genes. The results show high performance of the model’s predictions, with an AUC of 0.89. One hundred and twenty-one predicted pairs with high prediction scores were selected for inferring the functions using statistical enrichment analysis. Using this method, PF3D7_1248700 and PF3D7_0401800 were found to be involved with muscle contraction and striated muscle tissue development, while PF3D7_1303800 and PF3D7_1201000 were found to be related to protein dephosphorylation. In conclusion, combining a heterogeneous network and a hybrid deep learning technique can allow us to identify unknown gene functions of malaria parasites. This approach is generalized and can be applied to other diseases that enhance the field of biomedical science.

The first publication demonstrating that major depressive disorder (MDD) is associated with alterations in the gut microbiota appeared in 2008 (Maes et al., 2008). The purpose of the present study is to delineate a) the microbiome signature of the phenome of depression, including suicidal behaviours (SB) and cognitive deficits; the effects of adverse childhood experiences (ACEs) and recurrence of illness index (ROI) on the microbiome; and the microbiome signature of lowered high-density lipoprotein cholesterol (HDLc). We determined isometric log-ratio abundances or prevalences of gut microbiome phyla, genera, and species by analysing stool samples from 37 healthy Thai controls and 32 MDD patients using 16S rDNA sequencing. Six microbiome taxa accounted for 36% of the variance in the depression phenome, namely Hungatella and Fusicatenibacter (positive associations) and Butyricicoccus, Clostridium, Parabacteroides merdae , and Desulfovibrio piger (inverse association). This profile (labelled enterotype 1) indicates compositional dysbiosis, is strongly predicted by ACE and ROI, and is linked to SB. A second enterotype was developed that predicted a decrease in HDLc and an increase in the atherogenic index of plasma (Bifidobacterium, P. merdae , and Romboutsi a were positively associated, while Proteobacteria and Clostridium sensu stricto were negatively associated). Together, enterotypes 1 and 2 explained 40.4% of the variance in the depression phenome, and enterotype 1 in conjunction with HDLc explained 39.9% of the variance in current SB. In conclusion, the microimmuneoxysome is a potential new drug target for the treatment of severe depression and SB and possibly for the prevention of future episodes. [ABSTRACT FROM AUTHOR]

Drug repositioning has been proposed to develop drugs for diseases. However, the similarity in a single aspect may not be sufficient to reveal hidden information. Therefore, we established protein–protein similarity vectors (PPSVs) based on potential similarities in various types of biological information associated with proteins, including their network topology, proteomic data, functional analysis, and druggable property. Based on the proposed PPSVs, a separate drug–disease matrix was constructed for individual to prevent characteristics from being obscured between diseases. The classification technique was employed for prediction. The results showed that more than half of the tested disease models exhibited high performance, with overall F1 scores of more than 80%. Furthermore, comparing all diseases using traditional methods in one run, we obtained an (area under the curve) AUC of 98.9%. All candidate drugs were then tested in clinical trials (p-value < 2.2 × 10−16) and were known drugs based on their functions (p-value < 0.05). An analysis revealed that, in the functional aspect, the confidence value of an interaction in the protein–protein interaction network and the functional pathway score were the best descriptors for prediction. Based on the learning processes of PPSVs with an isolated disease, the classifier exhibited high performance in predicting and identifying new potential drugs for that disease.

Drug target prediction is an important method for drug discovery and design, can disclose the potential inhibitory effect of active compounds, and is particularly relevant to many diseases that have the potential to kill, such as dengue, but lack any healing agent. An antiviral drug is urgently required for dengue treatment. Some potential antiviral agents are still in the process of drug discovery, but the development of more effective active molecules is in critical demand. Herein, we aimed to provide an efficient technique for target prediction using homopharma and network-based methods, which is reliable and expeditious to hunt for the possible human targets of three phenolic lipids (anarcardic acid, cardol, and cardanol) related to dengue viral (DENV) infection as a case study. Using several databases, the similarity search and network-based analyses were applied on the three phenolic lipids resulting in the identification of seven possible targets as follows. Based on protein annotation, three phenolic lipids may interrupt or disturb the human proteins, namely KAT5, GAPDH, ACTB, and HSP90AA1, whose biological functions have been previously reported to be involved with viruses in the family Flaviviridae. In addition, these phenolic lipids might inhibit the mechanism of the viral proteins: NS3, NS5, and E proteins. The DENV and human proteins obtained from this study could be potential targets for further molecular optimization on compounds with a phenolic lipid core structure in anti-dengue drug discovery. As such, this pipeline could be a valuable tool to identify possible targets of active compounds.

Integration of multiple sources and data levels provides a great insight into the complex associations between human and malaria systems. In this study, a meta-analysis framework was developed based on a heterogeneous network model for integrating human-malaria protein similarities, a human protein interaction network, and a Plasmodium vivax protein interaction network. An iterative network propagation was performed on the heterogeneous network until we obtained stabilized weights. The association scores were calculated for qualifying a novel potential human-malaria protein association. This method provided a better performance compared to random experiments. After that, the stabilized network was clustered into association modules. The potential association candidates were then thoroughly analyzed by statistical enrichment analysis with protein complexes and known drug targets. The most promising target proteins were the succinate dehydrogenase protein complex in the human citrate (TCA) cycle pathway and the nicotinic acetylcholine receptor in the human central nervous system. Promising associations and potential drug targets were also provided for further studies and designs in therapeutic approaches for malaria at a systematic level. In conclusion, this method is efficient to identify new human-malaria protein associations and can be generalized to infer other types of association studies to further advance biomedical science.

Protein-protein interactions integrated with disease-gene associations represent important information for revealing protein functions under disease conditions to improve the prevention, diagnosis, and treatment of complex diseases. Although several studies have attempted to identify disease-gene associations, the number of possible disease-gene associations is very small. High-throughput technologies have been established experimentally to identify the association between genes and diseases. However, these techniques are still quite expensive, time consuming, and even difficult to perform. Thus, based on currently available data and knowledge, computational methods have served as alternatives to provide more possible associations to increase our understanding of disease mechanisms. Here, a new network-based algorithm, namely, Disease-Gene Association (DGA), was developed to calculate the association score of a query gene to a new possible set of diseases. First, a large-scale protein interaction network was constructed, and the relationship between two interacting proteins was calculated with regard to the disease relationship. Novel plausible disease-gene pairs were identified and statistically scored by our algorithm using neighboring protein information. The results yielded high performance for disease-gene prediction, with an F-measure of 0.78 and an AUC of 0.86. To identify promising candidates of disease-gene associations, the association coverage of genes and diseases were calculated and used with the association score to perform gene and disease selection. Based on gene selection, we identified promising pairs that exhibited evidence related to several important diseases, e.g., inflammation, lipid metabolism, inborn errors, xanthomatosis, cerebellar ataxia, cognitive deterioration, malignant neoplasms of the skin and malignant tumors of the cervix. Focusing on disease selection, we identified target genes that were important to blistering skin diseases and muscular dystrophy. In summary, our developed algorithm is simple, efficiently identifies disease-gene associations in the protein-protein interaction network and provides additional knowledge regarding disease-gene associations. This method can be generalized to other association studies to further advance biomedical science.