Gene networks for use in metabolomic data analysis of blood plasma from patients with postoperative delirium

Postoperative delirium (POD) is considered one of the most severe complications, resulting in impaired cognitive function, extended hospitalization, and higher treatment costs. The challenge of early POD diagnosis becomes particularly significant in cardiac surgery cases, as the incidence of this complication exceeds 50 % in certain patient categories. While it is known that neuroinflammation, neurotransmitter imbalances, disruptions in neuroendocrine regulation, and interneuronal connections contribute significantly to the development of POD, the molecular, genetic mechanisms of POD in cardiac surgery patients, along with potential metabolomic diagnostic markers, remain inadequately understood. In this study, blood plasma was collected from a group of patients over 65 years old after cardiac surgery involving artificial circulation. The collected samples were analyzed for sphingomyelin content and quantity using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS/MS) methods. The analysis revealed four significantly different sphingomyelin contents in patients with POD compared to those who did not develop POD (control group). Employing gene network reconstruction, we perceived a set of 82 regulatory enzymes affiliated with the genetic coordination of the sphingolipid metabolism pathway. Within this set, 47 are assumed to be regulators of gene expression, governing the transcription of enzymes pivotal to the metabolic cascade. Complementing this, an additional assembly of 35 regulators are considered to be regulators of activity, degradation, and translocation dynamics of enzymes integral to the aforementioned pathway. Analysis of the overrepresentation of diseases with which these regulatory proteins are associated showed that the regulators can be categorized into two groups, associated with cardiovascular pathologies (CVP) and neuropsychiatric diseases (NPD), respectively. The regulators associated with CVP are expectedly related to the effects on myocardial tissue during surgery. It is hypothesized that dysfunction of NPD-associated regulators may specifically account for the development of POD after cardiac surgery. Thus, the identified regulatory genes may provide a basis for planning further experiments, in order to study disorders at the level of expression of these genes, as well as impaired function of proteins encoded by them in patients with POD. The identified significant sphingolipids can be considered as potential markers of POD.

Associative gene networks reveal novel candidates important for ADHD and dyslexia comorbidity

BACKGROUND

Attention deficit hyperactivity disorder (ADHD) is commonly associated with developmental dyslexia (DD), which are both prevalent and complicated pediatric neurodevelopmental disorders that have a significant influence on children's learning and development. Clinically, the comorbidity incidence of DD and ADHD is between 25 and 48%. Children with DD and ADHD may have more severe cognitive deficiencies, a poorer level of schooling, and a higher risk of social and emotional management disorders. Furthermore, patients with this comorbidity are frequently treated for a single condition in clinical settings, and the therapeutic outcome is poor. The development of effective treatment approaches against these diseases is complicated by their comorbidity features. This is often a major problem in diagnosis and treatment. In this study, we developed bioinformatical methodology for the analysis of the comorbidity of these two diseases. As such, the search for candidate genes related to the comorbid conditions of ADHD and DD can help in elucidating the molecular mechanisms underlying the comorbid condition, and can also be useful for genotyping and identifying new drug targets.

RESULTS

Using the ANDSystem tool, the reconstruction and analysis of gene networks associated with ADHD and dyslexia was carried out. The gene network of ADHD included 599 genes/proteins and 148,978 interactions, while that of dyslexia included 167 genes/proteins and 27,083 interactions. When the ANDSystem and GeneCards data were combined, a total of 213 genes/proteins for ADHD and dyslexia were found. An approach for ranking genes implicated in the comorbid condition of the two diseases was proposed. The approach is based on ten criteria for ranking genes by their importance, including relevance scores of association between disease and genes, standard methods of gene prioritization, as well as original criteria that take into account the characteristics of an associative gene network and the presence of known polymorphisms in the analyzed genes. Among the top 20 genes with the highest priority DRD2, DRD4, CNTNAP2 and GRIN2B are mentioned in the literature as directly linked with the comorbidity of ADHD and dyslexia. According to the proposed approach, the genes OPRM1, CHRNA4 and SNCA had the highest priority in the development of comorbidity of these two diseases. Additionally, it was revealed that the most relevant genes are involved in biological processes related to signal transduction, positive regulation of transcription from RNA polymerase II promoters, chemical synaptic transmission, response to drugs, ion transmembrane transport, nervous system development, cell adhesion, and neuron migration.

CONCLUSIONS

The application of methods of reconstruction and analysis of gene networks is a powerful tool for studying the molecular mechanisms of comorbid conditions. The method put forth to rank genes by their importance for the comorbid condition of ADHD and dyslexia was employed to predict genes that play key roles in the development of the comorbid condition. The results can be utilized to plan experiments for the identification of novel candidate genes and search for novel pharmacological targets.

[Primary and Secondary micro-RNA Modulation the Extrinsic Pathway of Apoptosis in Hepatocellular Carcinoma]

One of the most common malignant liver diseases is hepatocellular carcinoma, which has a high recurrence rate and a low five-year survival rate. It is very heterogeneous both in structure and between patients, which complicates the diagnosis, prognosis and response to treatment. In this regard, an individualized, patient-centered approach becomes important, in which the use of mimetics and hsa-miRNA inhibitors involved in the pathogenesis of the disease may be determinative. From this point of view hsa-miRNAs are of interest, their aberrant expression is associated with poor prognosis for patients and is associated with tumor progression due to dysregulation of programmed cell death (apoptosis). However, the effect of hsa-miRNA on tumor development depends not only on its direct effect on expression of genes, the primary targets, but also on secondary targets mediated by regulatory pathways. While the former are actively studied, the role of secondary targets of these hsa-miRNAs in modulating apoptosis is still unclear. The present work summarizes data on hsa-miRNAs whose primary targets are key genes of the extrinsic pathway of apoptosis. Their aberrant expression is associated with early disease relapse and poor patient outcome. For these hsa-miRNAs, using the software package ANDSystem, we reconstructed the regulation of the expression of secondary targets and analyzed their impact on the activity of the extrinsic pathway of apoptosis. The potential effect of hsa-miRNAs mediated by action on secondary targets is shown to negatively correlate with the number of primary targets. It is also shown that hsa-miR-373, hsa-miR-106b and hsa-miR-96 have the highest priority as markers of hepatocellular carcinoma, whose action on secondary targets enhances their anti-apoptotic effect.

Primary and Secondary micro-RNA Modulation the Extrinsic Pathway of Apoptosis in Hepatocellular Carcinoma

Abstract-One of the most common malignant liver diseases is hepatocellular carcinoma, which has a high recurrence rate and a low five-year survival rate. It is very heterogeneous both in structure and between patients, which complicates the diagnosis, prognosis and response to treatment. In this regard, an individualized, patient-centered approach becomes important, in which the use of mimetics and hsa-miRNA inhibitors involved in the pathogenesis of the disease may be determinative. From this point of view hsa-miRNAs are of interest, their aberrant expression is associated with poor prognosis for patients and is associated with tumor progression due to dysregulation of programmed cell death (apoptosis). However, the effect of hsa-miRNA on tumor development depends not only on its direct effect on expression of genes, the primary targets, but also on secondary targets mediated by regulatory pathways. While the former are actively studied, the role of secondary targets of these hsa-miRNAs in modulating apoptosis is still unclear. The present work summarizes data on hsa-miRNAs whose primary targets are key genes of the extrinsic pathway of apoptosis. Their aberrant expression is associated with early disease relapse and poor patient outcome. For these hsa-miRNAs, using the software package ANDSystem, we reconstructed the regulation of the expression of secondary targets and analyzed their impact on the activity of the extrinsic pathway of apoptosis. The potential effect of hsa-miRNAs mediated by action on secondary targets is shown to negatively correlate with the number of primary targets. It is also shown that hsa-miR-373, hsa-miR-106b and hsa-miR-96 have the highest priority as markers of hepatocellular carcinoma, whose action on secondary targets enhances their anti-apoptotic effect.

The New Version of the ANDDigest Tool with Improved AI-Based Short Names Recognition

The body of scientific literature continues to grow annually. Over 1.5 million abstracts of biomedical publications were added to the PubMed database in 2021. Therefore, developing cognitive systems that provide a specialized search for information in scientific publications based on subject area ontology and modern artificial intelligence methods is urgently needed. We previously developed a web-based information retrieval system, ANDDigest, designed to search and analyze information in the PubMed database using a customized domain ontology. This paper presents an improved ANDDigest version that uses fine-tuned PubMedBERT classifiers to enhance the quality of short name recognition for molecular-genetics entities in PubMed abstracts on eight biological object types: cell components, diseases, side effects, genes, proteins, pathways, drugs, and metabolites. This approach increased average short name recognition accuracy by 13%.

Bioinformatic Reconstruction and Analysis of Gene Networks Related to Glucose Variability in Diabetes and Its Complications

Glucose variability (GV) has been recognized recently as a promoter of complications and therapeutic targets in diabetes. The aim of this study was to reconstruct and analyze gene networks related to GV in diabetes and its complications. For network analysis, we used the ANDSystem that provides automatic network reconstruction and analysis based on text mining. The network of GV consisted of 37 genes/proteins associated with both hyperglycemia and hypoglycemia. Cardiovascular system, pancreas, adipose and muscle tissues, gastrointestinal tract, and kidney were recognized as the loci with the highest expression of GV-related genes. According to Gene Ontology enrichment analysis, these genes are associated with insulin secretion, glucose metabolism, glycogen biosynthesis, gluconeogenesis, MAPK and JAK-STAT cascades, protein kinase B signaling, cell proliferation, nitric oxide biosynthesis, etc. GV-related genes were found to occupy central positions in the networks of diabetes complications (cardiovascular disease, diabetic nephropathy, retinopathy, and neuropathy) and were associated with response to hypoxia. Gene prioritization analysis identified new gene candidates (THBS1, FN1, HSP90AA1, EGFR, MAPK1, STAT3, TP53, EGF, GSK3B, and PTEN) potentially involved in GV. The results expand the understanding of the molecular mechanisms of the GV phenomenon in diabetes and provide molecular markers and therapeutic targets for future research.

Prioritization of genes involved in endothelial cell apoptosis by their implication in lymphedema using an analysis of associative gene networks with ANDSystem

BACKGROUND

Currently, more than 150 million people worldwide suffer from lymphedema. It is a chronic progressive disease characterized by high-protein edema of various parts of the body due to defects in lymphatic drainage. Molecular-genetic mechanisms of the disease are still poorly understood. Beginning of a clinical manifestation of primary lymphedema in middle age and the development of secondary lymphedema after treatment of breast cancer can be genetically determined. Disruption of endothelial cell apoptosis can be considered as one of the factors contributing to the development of lymphedema. However, a study of the relationship between genes associated with lymphedema and genes involved in endothelial apoptosis, in the associative gene network was not previously conducted.

METHODS

In the current work, we used well-known methods (ToppGene and Endeavour), as well as methods previously developed by us, to prioritize genes involved in endothelial apoptosis and to find potential participants of molecular-genetic mechanisms of lymphedema among them. Original methods of prioritization took into account the overrepresented Gene Ontology biological processes, the centrality of vertices in the associative gene network, describing the interactions of endothelial apoptosis genes with genes associated with lymphedema, and the association of the analyzed genes with diseases that are comorbid to lymphedema.

RESULTS

An assessment of the quality of prioritization was performed using criteria, which involved an analysis of the enrichment of the top-most priority genes by genes, which are known to have simultaneous interactions with lymphedema and endothelial cell apoptosis, as well as by genes differentially expressed in murine model of lymphedema. In particular, among genes involved in endothelial apoptosis, KDR, TNF, TEK, BMPR2, SERPINE1, IL10, CD40LG, CCL2, FASLG and ABL1 had the highest priority. The identified priority genes can be considered as candidates for genotyping in the studies involving the search for associations with lymphedema.

CONCLUSIONS

Analysis of interactions of these genes in the associative gene network of lymphedema can improve understanding of mechanisms of interaction between endothelial apoptosis and lymphangiogenesis, and shed light on the role of disturbance of these processes in the development of edema, chronic inflammation and connective tissue transformation during the progression of the disease.

A new version of the ANDSystem tool for automatic extraction of knowledge from scientific publications with expanded functionality for reconstruction of associative gene networks by considering tissue-specific gene expression

BACKGROUND

Consideration of tissue-specific gene expression in reconstruction and analysis of molecular genetic networks is necessary for a proper description of the processes occurring in a specified tissue. Currently, there are a number of computer systems that allow the user to reconstruct molecular-genetic networks using the data automatically extracted from the texts of scientific publications. Examples of such systems are STRING, Pathway Commons, MetaCore and Ingenuity. The MetaCore and Ingenuity systems permit taking into account tissue-specific gene expression during the reconstruction of gene networks. Previously, we developed the ANDSystem tool, which also provides an automated extraction of knowledge from scientific texts and allows the reconstruction of gene networks. The main difference between our system and other tools is in the different types of interactions between objects, which makes the ANDSystem complementary to existing well-known systems. However, previous versions of the ANDSystem did not contain any information on tissue-specific expression.

RESULTS

A new version of the ANDSystem has been developed. It offers the reconstruction of associative gene networks while taking into account the tissue-specific gene expression. The ANDSystem knowledge base features information on tissue-specific expression for 272 tissues. The system allows the reconstruction of combined gene networks, as well as performing the filtering of genes from such networks using the information on their tissue-specific expression. As an example of the application of such filtering, the gene network of the extrinsic apoptotic signaling pathway was analyzed. It was shown that considering different tissues can lead to changes in gene network structure, including changes in such indicators as betweenness centrality of vertices, clustering coefficient, network centralization, network density, etc. CONCLUSIONS: The consideration of tissue specificity can play an important role in the analysis of gene networks, in particular solving the problem of finding the most significant central genes. Thus, the new version of ANDSystem can be employed for a wide range of tasks related to biomedical studies of individual tissues. It is available at http://www-bionet.sscc.ru/and/cell /.

Search for New Candidate Genes Involved in the Comorbidity of Asthma and Hypertension Based on Automatic Analysis of Scientific Literature

Comorbid states of diseases significantly complicate diagnosis and treatment. Molecular mechanisms of comorbid states of asthma and hypertension are still poorly understood. Prioritization is a way for identifying genes involved in complex phenotypic traits. Existing methods of prioritization consider genetic, expression and evolutionary data, molecular-genetic networks and other. In the case of molecular-genetic networks, as a rule, protein-protein interactions and KEGG networks are used. ANDSystem allows reconstructing associative gene networks, which include more than 20 types of interactions, including protein-protein interactions, expression regulation, transport, catalysis, etc. In this work, a set of genes has been prioritized to find genes potentially involved in asthma and hypertension comorbidity. The prioritization was carried out using well-known methods (ToppGene and Endeavor) and a cross-talk centrality criterion, calculated by analysis of associative gene networks from ANDSystem. The identified genes, including IL1A, CD40LG, STAT3, IL15, FAS, APP, TLR2, C3, IL13 and CXCL10, may be involved in the molecular mechanisms of comorbid asthma/hypertension. An analysis of the dynamics of the frequency of mentioning the most priority genes in scientific publications revealed that the top 100 priority genes are significantly enriched with genes with increased positive dynamics, which may be a positive sign for further studies of these genes.

Novel candidate genes important for asthma and hypertension comorbidity revealed from associative gene networks

BACKGROUND

Hypertension and bronchial asthma are a major issue for people's health. As of 2014, approximately one billion adults, or ~ 22% of the world population, have had hypertension. As of 2011, 235-330 million people globally have been affected by asthma and approximately 250,000-345,000 people have died each year from the disease. The development of the effective treatment therapies against these diseases is complicated by their comorbidity features. This is often a major problem in diagnosis and their treatment. Hence, in this study the bioinformatical methodology for the analysis of the comorbidity of these two diseases have been developed. As such, the search for candidate genes related to the comorbid conditions of asthma and hypertension can help in elucidating the molecular mechanisms underlying the comorbid condition of these two diseases, and can also be useful for genotyping and identifying new drug targets.

RESULTS

Using ANDSystem, the reconstruction and analysis of gene networks associated with asthma and hypertension was carried out. The gene network of asthma included 755 genes/proteins and 62,603 interactions, while the gene network of hypertension - 713 genes/proteins and 45,479 interactions. Two hundred and five genes/proteins and 9638 interactions were shared between asthma and hypertension. An approach for ranking genes implicated in the comorbid condition of two diseases was proposed. The approach is based on nine criteria for ranking genes by their importance, including standard methods of gene prioritization (Endeavor, ToppGene) as well as original criteria that take into account the characteristics of an associative gene network and the presence of known polymorphisms in the analysed genes. According to the proposed approach, the genes IL10, TLR4, and CAT had the highest priority in the development of comorbidity of these two diseases. Additionally, it was revealed that the list of top genes is enriched with apoptotic genes and genes involved in biological processes related to the functioning of central nervous system.

CONCLUSIONS

The application of methods of reconstruction and analysis of gene networks is a productive tool for studying the molecular mechanisms of comorbid conditions. The method put forth to rank genes by their importance to the comorbid condition of asthma and hypertension was employed that resulted in prediction of 10 genes, playing the key role in the development of the comorbid condition. The results can be utilised to plan experiments for identification of novel candidate genes along with searching for novel pharmacological targets.

NACE: A web-based tool for prediction of intercompartmental efficiency of human molecular genetic networks

Molecular genetic processes generally involve proteins from distinct intracellular localisations. Reactions that follow the same process are distributed among various compartments within the cell. In this regard, the reaction rate and the efficiency of biological processes can depend on the subcellular localisation of proteins. Previously, the authors proposed a method of evaluating the efficiency of biological processes based on the analysis of the distribution of protein subcellular localisation (Popik et al., 2014). Here, NACE is presented, which is an open access web-oriented program that implements this method and allows the user to evaluate the intercompartmental efficiency of human molecular genetic networks. The method has been extended by a new feature that provides the evaluation of the tissue-specific efficiency of networks for more than 2800 anatomical structures. Such assessments are important in cases when molecular genetic pathways in different tissues proceed with the participation of various proteins with a number of intracellular localisations. For example, an analysis of KEGG pathways, conducted using the developed program, showed that the efficiencies of many KEGG pathways are tissue-specific. Analysis of efficiencies of regulatory pathways in the liver, linking proteins of the hepatitis C virus with human proteins involved in the KEGG apoptosis pathway, showed that intercompartmental efficiency might play an important role in host-pathogen interactions. Thus, the developed tool can be useful in the study of the effectiveness of functioning of various molecular genetic networks, including metabolic, regulatory, host-pathogen interactions and others taking into account tissue-specific gene expression. The tool is available via the following link: http://www-bionet.sscc.ru/nace/.

Interactome of the hepatitis C virus: Literature mining with ANDSystem

A study of the molecular genetics mechanisms of host-pathogen interactions is of paramount importance in developing drugs against viral diseases. Currently, the literature contains a huge amount of information that describes interactions between HCV and human proteins. In addition, there are many factual databases that contain experimentally verified data on HCV-host interactions. The sources of such data are the original data along with the data manually extracted from the literature. However, the manual analysis of scientific publications is time consuming and, because of this, databases created with such an approach often do not have complete information. One of the most promising methods to provide actualisation and completeness of information is text mining. Here, with the use of a previously developed method by the authors using ANDSystem, an automated extraction of information on the interactions between HCV and human proteins was conducted. As a data source for the text mining approach, PubMed abstracts and full text articles were used. Additionally, external factual databases were analyzed. On the basis of this analysis, a special version of ANDSystem, extended with the HCV interactome, was created. The HCV interactome contains information about the interactions between 969 human and 11 HCV proteins. Among the 969 proteins, 153 'new' proteins were found not previously referred to in any external databases of protein-protein interactions for HCV-host interactions. Thus, the extended ANDSystem possesses a more comprehensive detailing of HCV-host interactions versus other existing databases. It was interesting that HCV proteins more preferably interact with human proteins that were already involved in a large number of protein-protein interactions as well as those associated with many diseases. Among human proteins of the HCV interactome, there were a large number of proteins regulated by microRNAs. It turned out that the results obtained for protein-protein interactions and microRNA-regulation did not depend on how well the proteins were studied, while protein-disease interactions appeared to be dependent on the level of study. In particular, the mean number of diseases linked to well-studied proteins (proteins were considered well-studied if they were mentioned in 50 or more PubMed publications) from the HCV interactome was 20.8, significantly exceeding the mean number of associations with diseases (10.1) for the total set of well-studied human proteins present in ANDSystem. For proteins not highly poorly-studied investigated, proteins from the HCV interactome (each protein was referred to in less than 50 publications) distribution of the number of diseases associated with them had no statistically significant differences from the distribution of the number of diseases associated with poorly-studied proteins based on the total set of human proteins stored in ANDSystem. With this, the average number of associations with diseases for the HCV interactome and the total set of human proteins were 0.3 and 0.2, respectively. Thus, ANDSystem, extended with the HCV interactome, can be helpful in a wide range of issues related to analyzing HCV-host interactions in the search for anti-HCV drug targets. The demo version of the extended ANDSystem covered here containing only interactions between human proteins, genes, metabolites, diseases, miRNAs and molecular-genetic pathways, as well as interactions between human proteins/genes and HCV proteins, is freely available at the following web address: http://www-bionet.sscc.ru/psd/andhcv/.