Supervised machine learning models and protein-protein interaction network analysis of gene expression profiles induced by omega-3 polyunsaturated fatty acids

Background: Omega-3 polyunsaturated fatty acids (PUFAs), such as eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids, have beneficial effects on human health, but their effect on gene expression in elderly individuals (age ≥ 65) is largely unknown. In order to examine this, the gene expression profiles were analyzed in the healthy subjects (n = 96) at baseline and after 26 weeks of supplementation with EPA+DHA to determine up-regulated and down-regulated dif-ferentially expressed genes (DEGs) triggered by PUFAs. The protein-protein interaction (PPI) networks were constructed by mapping these DEGs to a human interactome and linking them to the specific pathways. Objective: This study aimed to implement supervised machine learning models and protein-protein interaction network analysis of gene expression profiles induced by PUFAs. Methods: The transcriptional profile of GSE12375 was obtained from the Gene Expression Om-nibus database, which is based on the Affymetrix NuGO array. The probe cell intensity data were converted into the gene expression values, and the background correction was performed by the multi-array average algorithm. The LIMMA (Linear Models for Microarray Data) algo-rithm was implemented to identify relevant DEGs at baseline and after 26 weeks of supplemen-tation with a p-value < 0.05. The DAVID web server was used to identify and construct the en-riched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Finally, the construction of machine learning (ML) models, including logistic regression, naïve Bayes, and deep neural networks, were implemented for the analyzed DEGs associated with the specific pathways. Results: The results revealed that up-regulated DEGs were associated with neurotrophin/MAPK signaling, whereas the down-regulated DEGs were linked to cancer, acute myeloid leukemia, and long-term depression pathways. Additionally, ML approaches were able to cluster the EPA/DHA-treated and control groups by the logistic regression performing the best. Conclusion: Overall, this study highlights the pivotal changes in DEGs induced by PUFAs and provides the rationale for the implementation of ML algorithms as predictive models for this type of biomedical data.

Co-expression analysis of placental genes in the search for key signaling pathways and biomarkers of the great obstetrical syndromes

Objective. To study the molecular mechanisms responsible for the development of diseases grouped within the great obstetrical syndromes (GOS) at the level of the transcriptome of human maternal placenta.Material and Methods. We gathered the results of genome-wide transcriptome studies of the human placental tissue using Gene Expression Omnibus (GEO) data repository for the following phenotypes: physiological pregnancy, preeclampsia (PE), premature birth, and intrauterine growth restriction (IUGR). Eleven data sets were selected and supplemented with our experimental data; a total of 481 samples of human placental tissue were included in the integrative analysis. Bioinformatic data processing and statistical analyses were performed in the R v3.6.1 software environment using the Bioconductor packages. The pooled dataset was used to search for common molecular targets for GOS via weighted gene co-expression network analysis (WGCNA). The functional annotation of genes and the resulting clusters was carried out with the DAVID database; protein-protein interaction network was built using the STRING software; and the hub genes for the network were identified using the MCC analysis with plugin cytoHubba in Cytoscape software 3.7.2.Results. We obtained a table of expression levels for 15,167 genes in 246 samples. Hierarchical clustering of this network allowed to find 55 modules of co-expressed genes in the group with PE, 109 modules in the group with PB, 75 modules in patients with IUGR, and 56 modules in the control group. The preservation analysis of co-expressed modules for the studied phenotypes suggested the presence of a common cluster comprising eight genes specific only for patients with PE and IUGR, as well as the module of 23 co-expressed genes typical only for patients with PB and IUGR. Protein-protein interaction network was built for these gene sets, and the SOD1, TXNRD1, and UBB genes were the central nodes in the network. Based on network topology evaluation with cytoHubba, six hub genes (rank ˂ 5) were identified as follows: SOD1, TKT, TXNRD1, GCLM, GOT1, and ACO1.Conclusion. The obtained results allowed to identify promising genetic markers for preeclampsia, intrauterine growth restriction, and miscarriage. Moreover, the study also made it possible to identify the most important overlapping molecular mechanisms of these diseases occurring in the placental tissue.

Integrated analysis and the identification of a circRNA-miRNA-mRNA network in the progression of abdominal aortic aneurysm

Background Abdominal aortic aneurysm (AAA) is a disease commonly seen in the elderly. The aneurysm diameter increases yearly, and the larger the AAA the higher the risk of rupture, increasing the risk of death. However, there are no current effective interventions in the early stages of AAA. Methods Four gene expression profiling datasets, including 23 normal artery (NOR) tissue samples and 97 AAA tissue samples, were integrated in order to explore potential molecular biological targets for early intervention. After preprocessing, differentially expressed genes (DEGs) between AAA and NOR were identified using LIMMA package. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis were conducted using the DAVID database. The protein-protein interaction network was constructed and hub genes were identified using the STRING database and plugins in Cytoscape. A circular RNA (circRNA) profile of four NOR tissues versus four AAA tissues was then reanalyzed. A circRNA-miRNA-mRNA interaction network was constructed after predictions were made using the Targetscan and Circinteractome databases. Results A total of 440 DEGs (263 up-regulated and 177 down-regulated) were identified in the AAA group, compared with the NOR group. The majority were associated with the extracellular matrix, tumor necrosis factor-α, and transforming growth factor-β. Ten hub gene-encoded proteins (namely IL6, RPS27A, JUN, UBC, UBA52, FOS, IL1B, MMP9, SPP1 and CCL2) coupled with a higher degree of connectivity hub were identified after protein‐protein interaction network analysis. Our results, in combination with the results of previous studies revealed that miR-635, miR-527, miR-520h, miR-938 and miR-518a-5p may be affected by circ_0005073 and impact the expression of hub genes such as CCL2, SPP1 and UBA52. The miR-1206 may also be affected by circ_0090069 and impact RPS27A expression. Conclusions This circRNA-miRNA-mRNA network may perform critical roles in AAA and may be a novel target for early intervention.

Endocytic proteins with prion-like domains form viscoelastic condensates that enable membrane remodeling

Membrane invagination and vesicle formation are key steps in endocytosis and cellular trafficking. Here, we show that endocytic coat proteins with prion-like domains (PLDs) form hemispherical puncta in the budding yeast, Saccharomyces cerevisiae. These puncta have the hallmarks of biomolecular condensates and organize proteins at the membrane for actin-dependent endocytosis. They also enable membrane remodeling to drive actin-independent endocytosis. The puncta, which we refer to as endocytic condensates, form and dissolve reversibly in response to changes in temperature and solution conditions. We find that endocytic condensates are organized around dynamic protein–protein interaction networks, which involve interactions among PLDs with high glutamine contents. The endocytic coat protein Sla1 is at the hub of the protein–protein interaction network. Using active rheology, we inferred the material properties of endocytic condensates. These experiments show that endocytic condensates are akin to viscoelastic materials. We use these characterizations to estimate the interfacial tension between endocytic condensates and their surroundings. We then adapt the physics of contact mechanics, specifically modifications of Hertz theory, to develop a quantitative framework for describing how interfacial tensions among condensates, the membrane, and the cytosol can deform the plasma membrane to enable actin-independent endocytosis.

Identifying potential key genes and existing drugs for Multiple sclerosis, Schizophrenia, and Autism- an in silico approach

Nowadays, neurological conditions are a major concern as it not only preys on a patients health but also is a huge economic burden that is placed on the patients family. The diagnosis and treatment of disease sometimes cause methodological limitations. This is mainly common for individuals who have the signs of MS and schizophrenia (SZ). Patients suffering from multiple sclerosis are more likely to develop schizophrenia. Besides, a significant portion of patients who have been diagnosed with Autism Spectrum Disorder (ASD) later acquire the symptoms of Schizophrenia. In this study, we used bioinformatics tools to determine differentially expressed genes (DEGs) in all these diseases, and then we created a protein-protein interaction network using the online software STRING and identified 15 significant genes with the help of Cytohubba a plug-in tool in Cytoscape, the offline software (version3.8.2). We then used a drug-gene interaction database to conduct a drug-gene interaction study of the 15 hub genes and from there we showed that there are 37 existing FDA-approved drugs were obtained. These findings may provide a new and common therapeutic approach for MS, SZ, and ASD therapy.

Interpretable network propagation with application to expanding the repertoire of human proteins that interact with SARS-CoV-2

Background Network propagation has been widely used for nearly 20 years to predict gene functions and phenotypes. Despite the popularity of this approach, little attention has been paid to the question of provenance tracing in this context, e.g., determining how much any experimental observation in the input contributes to the score of every prediction. Results We design a network propagation framework with 2 novel components and apply it to predict human proteins that directly or indirectly interact with SARS-CoV-2 proteins. First, we trace the provenance of each prediction to its experimentally validated sources, which in our case are human proteins experimentally determined to interact with viral proteins. Second, we design a technique that helps to reduce the manual adjustment of parameters by users. We find that for every top-ranking prediction, the highest contribution to its score arises from a direct neighbor in a human protein-protein interaction network. We further analyze these results to develop functional insights on SARS-CoV-2 that expand on known biology such as the connection between endoplasmic reticulum stress, HSPA5, and anti-clotting agents. Conclusions We examine how our provenance-tracing method can be generalized to a broad class of network-based algorithms. We provide a useful resource for the SARS-CoV-2 community that implicates many previously undocumented proteins with putative functional relationships to viral infection. This resource includes potential drugs that can be opportunistically repositioned to target these proteins. We also discuss how our overall framework can be extended to other, newly emerging viruses.

Comparative analysis of aneurysm subtypes associated genes based on protein–protein interaction network

The arterial aneurysm refers to localized dilation of blood vessel wall and is common in general population. The majority of aneurysm cases remains asymptomatic until a sudden rupture which is usually fatal and of extremely high mortality (~ 50–60%). Therefore, early diagnosis, prevention and management of aneurysm are in urgent need. Unfortunately, current understanding of disease driver genes of various aneurysm subtypes is still limited, and without appropriate biomarkers and drug targets no specialized drug has been developed for aneurysm treatment. In this research, aneurysm subtypes were analyzed based on protein–protein interaction network to better understand aneurysm pathogenesis. By measuring network-based proximity of aneurysm subtypes, we identified a relevant closest relationship between aortic aneurysm and aortic dissection. An improved random walk method was performed to prioritize candidate driver genes of each aneurysm subtype. Thereafter, transcriptomes of 6 human aneurysm subtypes were collected and differential expression genes were identified to further filter potential driver genes. Functional enrichment of above driver genes indicated a general role of ubiquitination and programmed cell death in aneurysm pathogenesis. Especially, we further observed participation of BCL-2-mediated apoptosis pathway and caspase-1 related pyroptosis in the development of cerebral aneurysm and aneurysmal subarachnoid hemorrhage in corresponding transcriptomes.

Transcriptome sequencing revealed the influence of blue light on the expression levels of light-stress response genes in Centella asiatica

Centella asiatica is rich in medical and cosmetic properties. While physiological responses of C. asiatica to light have been widely reported, the knowledge of the effects of light on its gene expression is sparse. In this study, we used RNA sequencing (RNA-seq) to investigate the expression of the C. asiatica genes in response to monochromatic red and blue light. Most of the differentially expressed genes (DEGs) under blue light were up-regulated but those under red light were down-regulated. The DEGs encoded for CRY-DASH and UVR3 were among up-regulated genes that play significant roles in responses under blue light. The DEGs involved in the response to photosystem II photodamages and in the biosynthesis of photoprotective xanthophylls were also up-regulated. The expression of flavonoid biosynthetic DEGs under blue light was up-regulated but that under red light was down-regulated. Correspondingly, total flavonoid content under blue light was higher than that under red light. The ABI5, MYB4, and HYH transcription factors appeared as hub nodes in the protein-protein interaction network of the DEGs under blue light while ERF38 was a hub node among the DEGs under red light. In summary, stress-responsive genes were predominantly up-regulated under blue light to respond to stresses that could be induced under high energy light. The information obtained from this study can be useful to better understand the responses of C. asiatica to different light qualities.