Gene Network Reconstruction Research Articles

AI-Assisted Identification of Primary and Secondary Metabolomic Markers for Postoperative Delirium.

Despite considerable investigative efforts, the molecular mechanisms of postoperative delirium (POD) remain unresolved. The present investigation employs innovative methodologies for identifying potential primary and secondary metabolic markers of POD by analyzing serum metabolomic profiles utilizing the genetic algorithm and artificial neural networks. The primary metabolomic markers constitute a combination of metabolites that optimally distinguish between POD and non-POD groups of patients. Our analysis revealed L-lactic acid, inositol, and methylcysteine as the most salient primary markers upon which the prediction accuracy of POD manifestation achieved AUC = 99%. The secondary metabolomic markers represent metabolites that exhibit perturbed correlational patterns within the POD group. We identified 54 metabolites as the secondary markers of POD, incorporating neurotransmitters such as gamma-aminobutyric acid (GABA) and serotonin. These findings imply a systemic disruption in metabolic processes in patients with POD. The deployment of gene network reconstruction techniques facilitated the postulation of hypotheses describing the role of established genomic POD markers in the molecular-genetic mechanisms of metabolic pathways dysregulation, and involving the identified primary and secondary metabolomic markers. This study not only expands the understanding of POD pathogenesis but also introduces a novel technology for the bioinformatic analysis of metabolomic data that could aid in uncovering potential primary and secondary markers in diverse research domains.

In Vitro vs. In Vivo Transcriptomic Approach Revealed Core Pathways of Nitrogen Deficiency Response in Tea Plant (Camellia sinensis (L.) Kuntze).

For the first time, we used an in vitro vs. in vivo experimental design to reveal core pathways under nitrogen deficiency (ND) in an evergreen tree crop. These pathways were related to lignin biosynthesis, cell redox homeostasis, the defense response to fungus, the response to Karrikin, amino acid transmembrane transport, the extracellular region, the cellular protein catabolic process, and aspartic-type endopeptidase activity. In addition, the mitogen-activated protein kinase pathway and ATP synthase (ATP)-binding cassette transporters were significantly upregulated under nitrogen deficiency in vitro and in vivo. Most of the MAPK downstream genes were related to calcium signaling (818 genes) rather than hormone signaling (157 genes). Moreover, the hormone signaling pathway predominantly contained auxin- and abscisic acid-related genes, indicating the crucial role of these hormones in ND response. Overall, 45 transcription factors were upregulated in both experiments, 5 WRKYs, 3 NACs, 2 MYBs, 2 ERFs, HD-Zip, RLP12, bHLH25, RADIALIS-like, and others, suggesting their ND regulation is independent from the presence of a root system. Gene network reconstruction displayed that these transcription factors participate in response to fungus/chitin, suggesting that nitrogen response and pathogen response have common regulation. The upregulation of lignin biosynthesis genes, cytochrome genes, and strigalactone response genes was much more pronounced under in vitro ND as compared to in vivo ND. Several cell wall-related genes were closely associated with cytochromes, indicating their important role in flavanols biosynthesis in tea plant. These results clarify the signaling mechanisms and regulation of the response to nitrogen deficiency in evergreen tree crops.

Concerted transcriptional regulation of the morphogenesis of hypothalamic neurons by ONECUT3

Acquisition of specialized cellular features is controlled by the ordered expression of transcription factors (TFs) along differentiation trajectories. Here, we find a member of the Onecut TF family, ONECUT3, expressed in postmitotic neurons that leave their Ascl1+/Onecut1/2+ proliferative domain in the vertebrate hypothalamus to instruct neuronal differentiation. We combined single-cell RNA-seq and gain-of-function experiments for gene network reconstruction to show that ONECUT3 affects the polarization and morphogenesis of both hypothalamic GABA-derived dopamine and thyrotropin-releasing hormone (TRH)+ glutamate neurons through neuron navigator-2 (NAV2). In vivo, siRNA-mediated knockdown of ONECUT3 in neonatal mice reduced NAV2 mRNA, as well as neurite complexity in Onecut3-containing neurons, while genetic deletion of Onecut3/ceh-48 in C. elegans impaired neurocircuit wiring, and sensory discrimination-based behaviors. Thus, ONECUT3, conserved across neuronal subtypes and many species, underpins the polarization and morphological plasticity of phenotypically distinct neurons that descend from a common pool of Ascl1+ progenitors in the hypothalamus.

RNA-Binding Proteins as Novel Effectors in Osteoblasts and Osteoclasts: A Systems Biology Approach to Dissect the Transcriptional Landscape.

Bone health is ensured by the coordinated action of two types of cells-the osteoblasts that build up bone structure and the osteoclasts that resorb the bone. The loss of balance in their action results in pathological conditions such as osteoporosis. Central to this study is a class of RNA-binding proteins (RBPs) that regulates the biogenesis of miRNAs. In turn, miRNAs represent a critical level of regulation of gene expression and thus control multiple cellular and biological processes. The impact of miRNAs on the pathobiology of various multifactorial diseases, including osteoporosis, has been demonstrated. However, the role of RBPs in bone remodeling is yet to be elucidated. The aim of this study is to dissect the transcriptional landscape of genes encoding the compendium of 180 RBPs in bone cells. We developed and applied a multi-modular integrative analysis algorithm. The core methodology is gene expression analysis using the GENEVESTIGATOR platform, which is a database and analysis tool for manually curated and publicly available transcriptomic data sets, and gene network reconstruction using the Ingenuity Pathway Analysis platform. In this work, comparative insights into gene expression patterns of RBPs in osteoblasts and osteoclasts were obtained, resulting in the identification of 24 differentially expressed genes. Furthermore, the regulation patterns upon different treatment conditions revealed 20 genes as being significantly up- or down-regulated. Next, novel gene-gene associations were dissected and gene networks were reconstructed. Additively, a set of osteoblast- and osteoclast-specific gene signatures were identified. The consolidation of data and information gained from each individual analytical module allowed nominating novel promising candidate genes encoding RBPs in osteoblasts and osteoclasts and will significantly enhance the understanding of potential regulatory mechanisms directing intracellular processes in the course of (patho)physiological bone turnover.

Experimental evolution for the recovery of growth loss due to genome reduction.

As the genome encodes the information crucial for cell growth, a sizeable genomic deficiency often causes a significant decrease in growth fitness. Whether and how the decreased growth fitness caused by genome reduction could be compensated by evolution was investigated here. Experimental evolution with an Escherichia coli strain carrying a reduced genome was conducted in multiple lineages for approximately 1000 generations. The growth rate, which largely declined due to genome reduction, was considerably recovered, associated with the improved carrying capacity. Genome mutations accumulated during evolution were significantly varied across the evolutionary lineages and were randomly localized on the reduced genome. Transcriptome reorganization showed a common evolutionary direction and conserved the chromosomal periodicity, regardless of highly diversified gene categories, regulons, and pathways enriched in the differentially expressed genes. Genome mutations and transcriptome reorganization caused by evolution, which were found to be dissimilar to those caused by genome reduction, must have followed divergent mechanisms in individual evolutionary lineages. Gene network reconstruction successfully identified three gene modules functionally differentiated, which were responsible for the evolutionary changes of the reduced genome in growth fitness, genome mutation, and gene expression, respectively. The diversity in evolutionary approaches improved the growth fitness associated with the homeostatic transcriptome architecture as if the evolutionary compensation for genome reduction was like all roads leading to Rome.

Experimental evolution for the recovery of growth loss due to genome reduction

Experimental evolution for the recovery of growth loss due to genome reduction

Prioritization of potential pharmacological targets for the development of anti-hepatocarcinoma drugs modulating the extrinsic apoptosis pathway: the reconstruction and analysis of associative gene networks help

Prioritization of potential pharmacological targets for the development of anti-hepatocarcinoma drugs modulating the extrinsic apoptosis pathway: the reconstruction and analysis of associative gene networks help

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 in adequately 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) me thods. 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.

Writers and readers of H3K9me2 form distinct protein networks during the cell cycle that include candidates for H3K9 mimicry.

Histone H3 lysine 9 methylation (H3K9me), which is written by the Euchromatic Histone Lysine Methyltransferases EHMT1 and EHMT2 and read by the heterochromatin protein 1 (HP1) chromobox (CBX) protein family, is dysregulated in many types of cancers. Approaches to inhibit regulators of this pathway are currently being evaluated for therapeutic purposes. Thus, knowledge of the complexes supporting the function of these writers and readers during the process of cell proliferation is critical for our understanding of their role in carcinogenesis. Here, we immunopurified each of these proteins and used mass spectrometry to define their associated non-histone proteins, individually and at two different phases of the cell cycle, namely G1/S and G2/M. Our findings identify novel binding proteins for these writers and readers, as well as corroborate known interactors, to show the formation of distinct protein complex networks in a cell cycle phase-specific manner. Furthermore, there is an organizational switch between cell cycle phases for interactions among specific writer-reader pairs. Through a multi-tiered bioinformatics-based approach, we reveal that many interacting proteins exhibit histone mimicry, based on an H3K9-like linear motif. Gene ontology analyses, pathway enrichment, and network reconstruction inferred that these comprehensive EHMT and CBX-associated interacting protein networks participate in various functions, including transcription, DNA repair, splicing, and membrane disassembly. Combined, our data reveals novel complexes that provide insight into key functions of cell cycle-associated epigenomic processes that are highly relevant for better understanding these chromatin-modifying proteins during cell cycle and carcinogenesis.

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

BackgroundAttention 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.ResultsUsing 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.ConclusionsThe 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.

MCPNet: a parallel maximum capacity-based genome-scale gene network construction framework.

Gene network reconstruction from gene expression profiles is a compute- and data-intensive problem. Numerous methods based on diverse approaches including mutual information, random forests, Bayesian networks, correlation measures, as well as their transforms and filters such as data processing inequality, have been proposed. However, an effective gene network reconstruction method that performs well in all three aspects of computational efficiency, data size scalability, and output quality remains elusive. Simple techniques such as Pearson correlation are fast to compute but ignore indirect interactions, while more robust methods such as Bayesian networks are prohibitively time consuming to apply to tens of thousands of genes. We developed maximum capacity path (MCP) score, a novel maximum-capacity-path-based metric to quantify the relative strengths of direct and indirect gene-gene interactions. We further present MCPNet, an efficient, parallelized gene network reconstruction software based on MCP score, to reverse engineer networks in unsupervised and ensemble manners. Using synthetic and real Saccharomyces cervisiae datasets as well as real Arabidopsis thaliana datasets, we demonstrate that MCPNet produces better quality networks as measured by AUPRC, is significantly faster than all other gene network reconstruction software, and also scales well to tens of thousands of genes and hundreds of CPU cores. Thus, MCPNet represents a new gene network reconstruction tool that simultaneously achieves quality, performance, and scalability requirements. Source code freely available for download at https://doi.org/10.5281/zenodo.6499747 and https://github.com/AluruLab/MCPNet, implemented in C++ and supported on Linux.

Apoptosis Genes as a Key to Identification of Inverse Comorbidity of Huntington's Disease and Cancer.

Cancer and neurodegenerative disorders present overwhelming challenges for healthcare worldwide. Epidemiological studies showed a decrease in cancer rates in patients with neurodegenerative disorders, including the Huntington disease (HD). Apoptosis is one of the most important processes for both cancer and neurodegeneration. We suggest that genes closely connected with apoptosis and associated with HD may affect carcinogenesis. We applied reconstruction and analysis of gene networks associated with HD and apoptosis and identified potentially important genes for inverse comorbidity of cancer and HD. The top 10 high-priority candidate genes included APOE, PSEN1, INS, IL6, SQSTM1, SP1, HTT, LEP, HSPA4, and BDNF. Functional analysis of these genes was carried out using gene ontology and KEGG pathways. By exploring genome-wide association study results, we identified genes associated with neurodegenerative and oncological disorders, as well as their endophenotypes and risk factors. We used publicly available datasets of HD and breast and prostate cancers to analyze the expression of the identified genes. Functional modules of these genes were characterized according to disease-specific tissues. This integrative approach revealed that these genes predominantly exert similar functions in different tissues. Apoptosis along with lipid metabolism dysregulation and cell homeostasis maintenance in the response to environmental stimulus and drugs are likely key processes in inverse comorbidity of cancer in patients with HD. Overall, the identified genes represent the promising targets for studying molecular relations of cancer and HD.

RNA-seq analysis and reconstruction of gene networks involved in response to salinity stress in quinoa (cv. Titicaca)

To better understand the mechanisms involved in salinity stress, the adaptability of quinoa cv. Titicaca—a halophytic plant—was investigated at the transcriptome level under saline and non-saline conditions. RNA-sequencing analysis of leaf tissue at the four-leaf stage by Illumina paired—end method was used to compare salt stress treatment (four days after stress at 13.8 dsm−1) and control. Among the obtained 30,846,354 transcripts sequenced, 30,303 differentially expressed genes from the control and stress treatment samples were identified, with 3363 genes expressed ≥ 2 and false discovery rate (FDR) of < 0.001. Six differential expression genes were then selected and qRT-PCR was used to confirm the RNA-seq results. Some of the genes (Include; CML39, CBSX5, TRX1, GRXC9, SnRKγ1 and BAG6) and signaling pathways discussed in this paper not been previously studied in quinoa. Genes with ≥ 2 were used to design the gene interaction network using Cytoscape software, and AgriGO software and STRING database were used for gene ontology. The results led to the identification of 14 key genes involved in salt stress. The most effective hub genes involved in salt tolerance were the heat shock protein gene family. The transcription factors that showed a significant increase in expression under stress conditions mainly belonged to the WRKY, bZIP and MYB families. Ontology analysis of salt stress-responsive genes and hub genes revealed that metabolic pathways, binding, cellular processes and cellular anatomical entity are among the most effective processes involved in salt stress.

Noncoding RNA of Zika Virus Affects Interplay between Wnt-Signaling and Pro-Apoptotic Pathways in the Developing Brain Tissue.

Zika virus (ZIKV) has a unique ability among flaviviruses to cross the placental barrier and infect the fetal brain causing severe abnormalities of neurodevelopment known collectively as congenital Zika syndrome. In our recent study, we demonstrated that the viral noncoding RNA (subgenomic flaviviral RNA, sfRNA) of the Zika virus induces apoptosis of neural progenitors and is required for ZIKV pathogenesis in the developing brain. Herein, we expanded on our initial findings and identified biological processes and signaling pathways affected by the production of ZIKV sfRNA in the developing brain tissue. We employed 3D brain organoids generated from induced human pluripotent stem cells (ihPSC) as an ex vivo model of viral infection in the developing brain and utilized wild type (WT) ZIKV (producing sfRNA) and mutant ZIKV (deficient in the production of sfRNA). Global transcriptome profiling by RNA-Seq revealed that the production of sfRNA affects the expression of >1000 genes. We uncovered that in addition to the activation of pro-apoptotic pathways, organoids infected with sfRNA-producing WT, but not sfRNA-deficient mutant ZIKV, which exhibited a strong down-regulation of genes involved in signaling pathways that control neuron differentiation and brain development, indicating the requirement of sfRNA for the suppression of neurodevelopment associated with the ZIKV infection. Using gene set enrichment analysis and gene network reconstruction, we demonstrated that the effect of sfRNA on pathways that control brain development occurs via crosstalk between Wnt-signaling and proapoptotic pathways.

BioGraph: Data Model for Linking and Querying Diverse Biological Metadata.

Studying the association of gene function, diseases, and regulatory gene network reconstruction demands data compatibility. Data from different databases follow distinct schemas and are accessible in heterogenic ways. Although the experiments differ, data may still be related to the same biological entities. Some entities may not be strictly biological, such as geolocations of habitats or paper references, but they provide a broader context for other entities. The same entities from different datasets can share similar properties, which may or may not be found within other datasets. Joint, simultaneous data fetching from multiple data sources is complicated for the end-user or, in many cases, unsupported and inefficient due to differences in data structures and ways of accessing the data. We propose BioGraph-a new model that enables connecting and retrieving information from the linked biological data that originated from diverse datasets. We have tested the model on metadata collected from five diverse public datasets and successfully constructed a knowledge graph containing more than 17 million model objects, of which 2.5 million are individual biological entity objects. The model enables the selection of complex patterns and retrieval of matched results that can be discovered only by joining the data from multiple sources.

Computer analysis of regulation of hepatocarcinoma marker genes hypermethylated by HCV proteins.

Hepatitis C virus (HCV) is a risk factor that leads to hepatocellular carcinoma (HCC) development. Epigenetic changes are known to play an important role in the molecular genetic mechanisms of virus-induced oncogenesis. Aberrant DNA methylation is a mediator of epigenetic changes that are closely associated with the HCC pathogenesis and considered a biomarker for its early diagnosis. The ANDSystem software package was used to reconstruct and evaluate the statistical significance of the pathways HCV could potentially use to regulate 32 hypermethylated genes in HCC, including both oncosuppressor and protumorigenic ones identified by genome-wide analysis of DNA methylation. The reconstructed pathways included those affecting protein-protein interactions (PPI), gene expression, protein activity, stability, and transport regulations, the expression regulation pathways being statistically significant. It has been shown that 8 out of 10 HCV proteins were involved in these pathways, the HCV NS3 protein being implicated in the largest number of regulatory pathways. NS3 was associated with the regulation of 5 tumor-suppressor genes, which may be the evidence of its central role in HCC pathogenesis. Analysis of the reconstructed pathways has demonstrated that following the transcription factor inhibition caused by binding to viral proteins, the expression of a number of oncosuppressors (WT1, MGMT, SOCS1, P53) was suppressed, while the expression of others (RASF1, RUNX3, WIF1, DAPK1) was activated. Thus, the performed gene-network reconstruction has shown that HCV proteins can influence not only the methylation status of oncosuppressor genes, but also their transcriptional regulation. The results obtained can be used in the search for pharmacological targets to develop new drugs against HCV-induced HCC.

Imputation of single-cell transcriptome data enables the reconstruction of networks predictive of breast cancer metastasis

Single-cell transcriptome data provide a unique opportunity to explore the gene networks of a particular cell type. However, insufficient capture rate and high dimensionality of single-cell RNA sequencing (scRNA-seq) data challenge cell-type-specific gene network (CGN) reconstruction. Here, we demonstrated that the imputation of scRNA-seq data enables reconstruction of CGNs by effective retrieval of gene functional associations. We reconstructed CGNs for seven primary and nine metastatic breast cancer cell lines using scRNA-seq data with imputation. Key genes for primary or metastatic cell lines were prioritized based on network centrality measures and CGN hub genes that were presumed to be the major determinant of cell type characteristics. To identify novel genes in breast cancer metastasis, we used the average rank difference of centrality between the primary and metastatic cell lines. Genes predicted using CGN centrality analysis were more enriched for known breast cancer metastatic genes than those predicted using differential expression. The molecular chaperone CCT2 was identified as a novel gene for breast metastasis during knockdown assays of several candidate genes. Overall, our study demonstrated an effective CGN reconstruction technique with imputation of scRNA-seq data and the feasibility of identifying key genes for particular cell subsets using single-cell network analysis.

Plasma metabolomics and gene regulatory networks analysis reveal the role of nonstructural SARS-CoV-2 viral proteins in metabolic dysregulation in COVID-19 patients

Metabolomic analysis of blood plasma samples from COVID-19 patients is a promising approach allowing for the evaluation of disease progression. We performed the metabolomic analysis of plasma samples of 30 COVID-19 patients and the 19 controls using the high-performance liquid chromatography (HPLC) coupled with tandem mass spectrometric detection (LC–MS/MS). In our analysis, we identified 103 metabolites enriched in KEGG metabolic pathways such as amino acid metabolism and the biosynthesis of aminoacyl-tRNAs, which differed significantly between the COVID-19 patients and the controls. Using ANDSystem software, we performed the reconstruction of gene networks describing the potential genetic regulation of metabolic pathways perturbed in COVID-19 patients by SARS-CoV-2 proteins. The nonstructural proteins of SARS-CoV-2 (orf8 and nsp5) and structural protein E were involved in the greater number of regulatory pathways. The reconstructed gene networks suggest the hypotheses on the molecular mechanisms of virus-host interactions in COVID-19 pathology and provide a basis for the further experimental and computer studies of the regulation of metabolic pathways by SARS-CoV-2 proteins. Our metabolomic analysis suggests the need for nonstructural protein-based vaccines and the control strategy to reduce the disease progression of COVID-19.

Gene Regulatory Network Characterization of Gastric Cancer's Histological Subtypes: Distinctive Biological and Clinically Relevant Master Regulators.

Simple SummaryThe prognosis of advanced gastric cancer patients remains unfavorable. Molecular heterogeneity has proven to be a major determinant of clinical outcomes. We characterized the transcriptome of the two major subgroups by highlighting the different biological and molecular pathways. We explored their association with clinicopathological features and survival. This comparative study aimed to define a reproducible in silico analysis so that the molecular mechanisms underlying carcinogenesis, disease natural history and the identification of new therapeutic targets can be traced.Gastric cancer (GC) molecular heterogeneity represents a major determinant for clinical outcomes, and although new molecular classifications have been introduced, they are not easy to translate from bench to bedside. We explored the data from GC public databases by performing differential gene expression analysis (DEGs) and gene network reconstruction to identify master regulators (MRs), as well as a gene set analysis (GSA) to reveal their biological features. Moreover, we evaluated the association of MRs with clinicopathological parameters. According to the GSA, the Diffuse group was characterized by an epithelial-mesenchymal transition (EMT) and inflammatory response, while the Intestinal group was associated with a cell cycle and drug resistance pathways. In particular, the regulons of Diffuse MRs, such as Vgll3 and Ciita, overlapped with the EMT and interferon-gamma response, while the regulons Top2a and Foxm1 were shared with the cell cycle pathways in the Intestinal group. We also found a strict association between MR activity and several clinicopathological features, such as survival. Our approach led to the identification of genes and pathways differentially regulated in the Intestinal and Diffuse GC histotypes, highlighting biologically interesting MRs and subnetworks associated with clinical features and prognosis, suggesting putative actionable candidates.

RECONSTRUCTION OF SCHIZOPHRENIA GENE NETWORK IN SEARCH FOR TARGET GENES

Nowadays, schizophrenia is a poorly understood disease with a variety of symptoms attributed to other malconditions, and controversial diagnosis without well-defined treatment. Target therapy implies disease gene network reconstruction, gene clustering, identification of gene ontology categories and genes with the largest number of network contacts.&#x0D; The aim of the study is to analyze schizophrenia-associated genes, determine their position in the gene network, establish their correlation, identify key genes related to the disease, and evaluate them as target genes for drug therapy.&#x0D; Materials and Methods. The authors analyzed currently relevant data on schizophrenia using such online databases as OMIM, PANTHER, DAVID, GeneMANIA, STRING-DB, and GeneCards. They calculated categories of gene ontologies for 200 genes, such as biological processes, molecular functions and cellular compartments that reflect schizophrenia impact on the transmission of neuronal impulses. The authors also visualized and built gene networks containing the identified key objects and their interaction, identified the most relevant schizophrenia genes (COMT, DISC1, HTR2A, NRXN1) and a strongly connected cluster, including such genes as BDNF, SLC6A4, HTR2A, HTR2C, CHRM1, SRC, AKT, YWHAE, DISC1, DRD2, COMT, NDEL1, NOS1, CAMK28, etc.&#x0D; Results. The biological interpretation of the results obtained is still a great challenge, since schizophrenia is a genetically complex disease with numerous causes and triggering events. Analysis of schizophrenia-associated genes, and identification of their position in the gene network (connectivity) makes it possible to find out their interaction, determine the key genes of the disease, and evaluate their prospects as target genes for drug therapy.