This study presents a novel method to assess the pathogenicity of pyrin protein mutations by using mutual information (MI) as a measure to quantify the correlation between residue motions or fluctuations and associated changes affecting the phenotype. The concept of MI profile shift is presented to quantify changes in MI upon mutation, revealing insights into residue-residue interactions at critical positions. We apply this method to the pyrin protein variants, which are associated with an autosomal recessively inherited disease called familial Mediterranean fever (FMF) since the available tools do not help predict the pathogenicity of the most penetrant variants. We demonstrate the utility of MI profile shifts in assessing the effects of mutations on protein stability, function, and disease phenotype. The importance of MI shifts, particularly the negative shifts observed in the pyrin example, as indicators of severe functional effects is emphasized. Additionally, the exploration of potential compensatory mechanisms suggested by positive MI shifts, which are otherwise random and inconsequential, is highlighted. The study also discusses challenges in relating MI profile changes to disease severity and advocates for comprehensive analysis considering genetic, environmental, and stochastic factors. Overall, this study provides insights into the molecular mechanisms underlying the pathogenesis of FMF and offers a framework for identifying potential therapeutic targets based on MI profile changes induced by mutations.

Serotonin-receptor binding plays a key role in several neurological and biological processes, including mood, sleep, hunger, cognition, learning, and memory. In this article, we performed molecular dynamics simulation to examine the key residues that play an essential role in the binding of serotonin to the G-protein-coupled 5-HT1B receptor (5HT1BR) via electrostatic interactions. Key residues for electrostatic interactions were identified via bond distance analysis and frustration analysis methods. An end-point free energy calculation method determines the stability of the 5-HT1BR due to serotonin binding. The single-point mutation of the polar/charged amino acid residues (Asp129, Thr134) on the binding sites and the calculation of binding free energy validate the quantitative contribution of these residues to the stability of the serotonin-receptor complex. The principal component analysis reflects that the serotonin-bound 5-HT1BR is more stabilized than the apo-receptor regarding dynamical changes. The difference dynamic cross-correlations map shows the correlation between the transmembranes and mini-Go, which indicates that the signal transduction happens between mini-Go and the receptor. Allosteric pathway analysis reveals the key nodes and key pathways for signal transduction in 5-HT1BR. These results provide useful insights into the study of signal transduction pathways and mutagenesis to regulate the binding and functionality of the complex. The developed protocols can be applied to study local non-covalent interactions and long-range allosteric communications in any protein-ligand system for computer-aided drug design.

To apply a network medicine-based approach to analyze the phenome of the prenatal fetal MRI and biometric findings in the Chiari II malformation (CM II) to detect specific patterns and co-occurrences.

A single-center retrospective review of fetal MRI scans obtained in fetuses with CM II was performed. Co-occurrence analysis was utilized to generate a phenotypic comorbidity matrix and visualized by Gephi software. Traditional univariate regression and geometric thin-plate spline methodology were used to elucidate the mechanisms underlying the relationships between morphometric measurements and geometric landmarks of the spine, skull, and brain deformations.

The CM II phenome consists of 35 nodes interconnected by 979 edges with a density of 0.828. Key "hubs" identified within this network include spinal bony defects, reduced posterior fossa dimensions, and vermis ectopia. The brain edema phenotype appearing only in the fetal stage but disappearing after postnatal surgery, links to increased postnatal morbidity and demonstrates distinct shape patterns by geometric analysis. Traditional univariate regression reveals correlations among spinal defects, posterior fossa dimensions, and caudal extent of vermis ectopia. The degree of brain rearrangement versus spinal bony rearrangement shows a correlation (r = 0.721, p = 0.0023) by partial least-squares analysis.

The CM II prenatal phenome is a multifaceted network centered around three key elements-spinal bony defects, small posterior fossa, and vermis ectopia-with strong interconnections. Fetal brain edema emerged as an exclusively prenatally detectable and transient phenotype of prognostic relevance.

Lung cancer is one of the most common cancer and the leading cause of cancer-related death worldwide. Early detection of lung cancer can help reduce the death rate; therefore, the identification of potential biomarkers is crucial. Thus, this study aimed to identify potential biomarkers for lung cancer by integrating bioinformatics analysis and machine learning (ML)-based approaches. Data were normalized using the robust multiarray average method and batch effect were corrected using the ComBat method. Differentially expressed genes were identified by the LIMMA approach and carcinoma-associated genes were selected using Enrichr, based on the DisGeNET database. Protein-protein interaction (PPI) network analysis was performed using STRING, and the PPI network was visualized using Cytoscape. The core hub genes were identified by overlapping genes obtained from degree, betweenness, closeness, and MNC. Moreover, the MCODE plugin for Cytoscape was used to perform module analysis, and optimal modules were selected based on MCODE scores along with their associated genes. Subsequently, Boruta-based ML approach was utilized to identify the important genes. Consequently, the core genes were identified by the overlapping genes obtained from PPI networks, module analysis, and ML-based approach. The prognostic and discriminative power analysis of the core genes was assessed through survival and ROC analysis. We extracted five datasets from USA cohort and three datasets from Taiwan cohort and performed same experimental protocols to determine potential biomarkers. Four genes (LPL, CLDN18, EDNRB, MME) were identified from USA cohort, while three genes (DNRB, MME, ROBO4) were from Taiwan cohort. Finally, two biomarkers (EDNRB and MME) were identified by intersecting genes, obtained from USA and Taiwan cohorts. The proposed biomarkers can significantly improve patient outcomes by enabling earlier detection, precise diagnosis, and tailored treatment, ultimately contributing to better survival rates and quality of life for patients.

Multi-omic studies provide comprehensive insight into biological systems by evaluating cellular changes between normal and pathological conditions at multiple levels of measurement. Biological networks, which represent interactions or associations between biomolecules, have been highly effective in facilitating omic analysis. However, current network-based methods lack generalizability to accommodate multiple data types across a range of diverse experiments.

We present AMEND 2.0, an updated active module identification method which can analyze multiplex and/or heterogeneous networks integrated with multi-omic data in a highly generalizable framework, in contrast to existing methods, which are mostly appropriate for at most two specific omic types. It is powered by Random Walk with Restart for multiplex-heterogeneous networks, with additional capabilities including degree bias adjustment and biased random walk for multi-objective module identification. AMEND was applied to two real-world multi-omic datasets: renal cell carcinoma data from The cancer genome atlas and an O-GlcNAc Transferase knockout study. Additional analyses investigate the performance of various subroutines of AMEND on tasks of node ranking and degree bias adjustment.

While the analysis of multi-omic datasets in a network context is poised to provide deeper understanding of health and disease, new methods are required to fully take advantage of this increasingly complex data. The current study combines several network analysis techniques into a single versatile method for analyzing biological networks with multi-omic data that can be applied in many diverse scenarios. Software is freely available in the R programming language at https://github.com/samboyd0/AMEND .

Colorectal cancer (CRC) holds the position of being the third most prevalent cancer and the second primary cause of cancer-related fatalities on a global scale. Approximately 65% of CRC patients survive for 5 years following diagnosis. Metastasis and recurrence frequently occur in half of CRC patients diagnosed at the late stage. This study used bioinformatics analysis to identify key signaling pathways, hub genes, transcription factors, and protein kinases involved in transforming primary CRC with liver metastasis potential. Prognostic markers in CRC were also identified.

The GSE81582 dataset was re-analyzed to identify differentially expressed genes (DEGs) in early CRC compared to non-tumoral tissues. A protein interaction network (PIN) was constructed, revealing significant modules and hub genes. Prognostic markers, transcription factors, and protein kinases were determined. Boxplot and gene set enrichment analyses were performed.

This study identified 1113 DEGs in primary CRC compared to healthy controls. PIN analysis revealed 75 hub genes and 8 significant clusters associated with early CRC. The down-regulation of SUCLG2 and KPNA2 correlated with poor prognosis. SIN3A and CDK6 played crucial roles in early CRC transformation, affecting rRNA processing pathways.

This study demonstrated several pathways, biological processes, and genes mediating the malignant transformation of healthy colorectal tissues to primary CRC and may help the prognosis and treatment of patients with early CRC.

Coral reefs are degrading at an accelerating rate owing to climate change. Understanding the heat stress tolerance of corals is vital for their sustainability. However, this tolerance varies substantially geographically, and information regarding coral responses across latitudes is lacking. In this study, we conducted a high temperature (34 °C) stress experiment on Porites lutea from tropical Xisha Islands (XS) and subtropical Daya Bay (DY) in the South China Sea (SCS). We compared physiological levels, antioxidant activities, and transcriptome sequencing to explore heat tolerance mechanisms and adaptive potential. At 34 °C, both XS and DY corals experienced significant bleaching and the physiological/biochemical index decreased, with XS corals exhibiting greater changes than DY corals. Transcriptome analysis revealed that coral hosts respond to heat stress mainly by boosting metabolic activity. The subtle transcriptional responses of zooxanthellae C15 underscored the host's pivotal role in thermal stress responses. DY coral hosts showed lower bleaching, stronger physiological plasticity, and higher temperature tolerance thresholds than XS, indicating superior heat tolerance. This superiority is linked to negative feedback transcriptional regulation strategies, including active environmental stress response and genetic information damage repair. The differences in thermal adaptability between tropical and subtropical P. lutea in the SCS may be attributed to their genetic differences and native habitat environments, suggesting that subtropical P. lutea may have the potential to adapt to future climate change. This study provides novel insights for predicting the fate of corals at different latitudes in terms of global warming and provides instructive guidance for coral reef ecological restoration.

Genome merging is a common phenomenon causing a wide range of consequences on phenotype, adaptation, and gene expression, yet its broader implications are not well-understood. Two consequences of genome merging on gene expression remain particularly poorly understood: dosage effects and evolution of expression. We employed Chlamydomonas reinhardtii as a model to investigate the effects of asymmetric genome merging by crossing a diploid with a haploid strain to create a novel triploid line. Five independent clonal lineages derived from this triploid line were evolved for 425 asexual generations in a laboratory natural selection experiment. Utilizing fitness assays, flow cytometry, and RNA-Seq, we assessed the immediate consequences of genome merging and subsequent evolution. Our findings reveal substantial alterations in genome size, gene expression, protein homeostasis, and cytonuclear stoichiometry. Gene expression exhibited expression-level dominance and transgressivity (i.e. expression level higher or lower than either parent). Ongoing expression-level dominance and a pattern of 'functional dominance' from the haploid parent was observed. Despite major genomic and nucleo-cytoplasmic disruptions, enhanced fitness was detected in the triploid strain. By comparing gene expression across generations, our results indicate that proteostasis restoration is a critical component of rapid adaptation following genome merging in Chlamydomonas reinhardtii and possibly other systems.

Amyloidosis diseases are characterized by protein misfolding, which forms insoluble beta-sheet fibrils progressively deposited in tissues. Deposition in the form of amyloid aggregates can occur in various organs, damaging their structure and function. The hallmark of amyloidosis is aberrant interactions leading to protein aggregation and proteotoxicity. Accordingly, amyloidosis-related samples represent a valuable source of information to generate new knowledge useful for diagnostic, prognostic, and therapeutic purposes. In this scenario, we outline the path to apply computational methods and strategies based on the combination of proteomics and systems biology approaches. In addition to algorithms useful for subtyping amyloid deposits or assessing proteome recovery after drug treatment, our chapter provides workflows based on protein-protein interaction and protein co-expression network models. In particular, the main steps to reconstruct and analyze them at both functional and topological levels are described. Our chapter aims to provide tools and instructions to identify and monitor prognostic, diagnostic, and therapeutic markers and to shed light on the processes, pathways, and functions affected by amyloidogenic proteins.

Understanding what drives protein abundance is essential to biology, medicine, and biotechnology. Driven by evolutionary selection, an amino acid sequence is tailored to meet the required abundance of a proteome, underscoring the intricate relationship between sequence and functional demand. Yet, the specific role of amino acid sequences in determining proteome abundance remains elusive. Here we show that the amino acid sequence alone encodes over half of protein abundance variation across all domains of life, ranging from bacteria to mouse and human. With an attempt to go beyond predictions, we trained a manageable-size Transformer model to interpret latent factors predictive of protein abundances. Intuitively, the model's attention focused on the protein's structural features linked to stability and metabolic costs related to protein synthesis. To probe these relationships, we introduce MGEM (Mutation Guided by an Embedded Manifold), a methodology for guiding protein abundance through sequence modifications. We find that mutations which increase predicted abundance have significantly altered protein polarity and hydrophobicity, underscoring a connection between protein structural features and abundance. Through molecular dynamics simulations we revealed that abundance-enhancing mutations possibly contribute to protein thermostability by increasing rigidity, which occurs at a lower synthesis cost.

The breakthrough discovery of novel biomarkers with prognostic and diagnostic value enables timely medical intervention for the survival of patients diagnosed with gastric cancer (GC). Typically, in studies focused on biomarker analysis, highly connected nodes (hubs) within the protein-protein interaction network (PPIN) are proposed as potential biomarkers. However, this study revealed an unexpected finding following the clustering of network nodes. Consequently, it is essential not to overlook weakly connected nodes (nonhubs) when determining suitable biomarkers from PPIN.

In this study, several potential biomarkers for GC were proposed based on the findings from RNA-sequencing (RNA-Seq) datasets, along with differential gene expression (DGE) analysis, PPINs, and weighted gene co-expression network analysis (WGCNA). Considering the overall survival (OS) analysis and the evaluation of expression levels alongside statistical parameters of the PPIN cluster nodes, it is plausible to suggest that THY1, CDH17, TGIF1, and AEBP1, categorized as nonhub nodes, along with ITGA5, COL1A1, FN1, and MMP2, identified as hub nodes, possess characteristics that render them applicable as biomarkers for the GC. Additionally, insulin-like growth factor (IGF)-binding protein-2 (IGFBP2), classified as a nonhub node, demonstrates a significant negative correlation with both groups within the same cluster. This observation underscores the conflicting findings regarding IGFBP2 in various cancer studies and enhances the potential of this gene to serve as a biomarker.

The findings of the current study not only identified the hubs and nonhubs that may serve as potential biomarkers for GC but also revealed a PPIN cluster that includes both hubs and nonhubs in conjunction with IGFBP2, thereby enhancing the understanding of the complex behavior associated with IGFBP2.

We employed LC-MS/MS to investigate the metabolic profiles of polysaccharide compounds in white, orange, and purple sweet potato flesh. Comparisons between Orange vs White, Purple vs Orange, and Purple vs White identified 69 polysaccharide metabolites, including 23, 36, and 44 differential metabolites, respectively, with distinct differentiation. Among the three sample groups, 14 polysaccharide compounds and 2 anthocyanins exhibited significant differences. Our further analysis indicated that anthocyanins occupy a central position in the related network diagram and are interconnected with polysaccharides. In metabolic pathways, sucrose and the anthocyanin precursor UDP-glucose were upregulated in purple sweet potatoes, along with elevated levels of pelargonidin 3-O-β-D-sambubioside and delphinidin 3,5-diglucoside. Conversely, sucrose was downregulated in purple sweet potatoes while increasing in white and orange varieties. Therefore, we hypothesize that the competition between sugars and anthocyanins for shared biosynthesis precursors is attributed to differential polysaccharide metabolites among sweet potato tubers of three colors.

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, significantly contributes to the phenotypic variance of complex traits. Traditionally, epistasis has been modeled using the Cartesian epistatic model, a multiplicative approach based on standard statistical regression. However, a recent study investigating epistasis in obesity-related traits has identified potential limitations of the Cartesian epistatic model, revealing that it likely only detects a fraction of the genetic interactions occurring in natural systems. In contrast, the exclusive-or (XOR) epistatic model has shown promise in detecting a broader range of epistatic interactions and revealing more biologically relevant functions associated with interacting variants. To investigate whether the XOR epistatic model also forms distinct network structures compared to the Cartesian model, we applied network science to examine genetic interactions underlying body mass index (BMI) in rats (Rattus norvegicus).

Our comparative analysis of XOR and Cartesian epistatic models in rats reveals distinct topological characteristics. The XOR model exhibits enhanced sensitivity to epistatic interactions between the network communities found in the Cartesian epistatic network, facilitating the identification of novel trait-related biological functions via community-based enrichment analysis. Additionally, the XOR network features triangle network motifs, indicative of higher-order epistatic interactions. This research also evaluates the impact of linkage disequilibrium (LD)-based edge pruning on network-based epistasis analysis, finding that LD-based edge pruning may lead to increased network fragmentation, which may hinder the effectiveness of network analysis for the investigation of epistasis. We confirmed through network permutation analysis that most XOR and Cartesian epistatic networks derived from the data display distinct structural properties compared to randomly shuffled networks.

Collectively, these findings highlight the XOR model's ability to uncover meaningful biological associations and higher-order epistasis derived from lower-order network topologies. The introduction of community-based enrichment analysis and motif-based epistatic discovery emphasize network science as a critical approach for advancing epistasis research and understanding complex genetic architectures.

The explosion of next-generation sequencing technologies has allowed researchers to move from studying single genes, to thousands of genes, and thereby to also consider the relationships within gene networks. Like others, we are interested in understanding how developmental and evolutionary forces shape the expression of individual genes, as well as the interactions among genes. To this end, we characterized the effects of genetic background and developmental environment on brain gene coexpression in two parallel, independent evolutionary lineages of Trinidadian guppies (Poecilia reticulata). We asked whether connectivity patterns among genes differed based on genetic background and rearing environment, and whether a gene's connectivity predicted its propensity for expression divergence. In pursuing these questions, we confronted the central challenge that standard approaches fail to control the Type I error and/or have low power in the presence of high dimensionality (i.e., large number of genes) and small sample size, as in many gene expression studies. Using our data as a case study, we detail central challenges, discuss sample size guidelines, and provide rigorous statistical approaches for exploring coexpression differences with small sample sizes. Using these approaches, we find evidence that coexpression relationships differ based on both genetic background and rearing environment. We report greater expression divergence in less connected genes and suggest this pattern may arise and be reinforced by selection.

Degree distributions in protein-protein interaction (PPI) networks are believed to follow a power law (PL). However, technical and study biases affect the experimental procedures for detecting PPIs. For instance, cancer-associated proteins have received disproportional attention. Moreover, bait proteins in large-scale experiments tend to have many false-positive interaction partners. Studying the degree distributions of thousands of PPI networks of controlled provenance, we address the question if PL distributions in observed PPI networks could be explained by these biases alone. Our findings are supported by mathematical models and extensive simulations, and indicate that study bias and technical bias suffice to produce the observed PL distribution. It is, hence, problematic to derive hypotheses about the topology of the true biological interactome from the PL distributions in observed PPI networks. Our study casts doubt on the use of the PL property of biological networks as a modeling assumption or quality criterion in network biology.

As an early indicator of dementia, mild cognitive impairment (MCI) requires specialized treatment according to its subtypes for the effective prevention and management of dementia progression. Based on the neuropathological characteristics, MCI can be classified into Alzheimer's disease (AD)-related cognitive impairment (ADCI) and subcortical vascular cognitive impairment (SVCI), being more likely to progress to AD and subcortical vascular dementia (SVD), respectively. For identifying MCI subtypes, plasma protein biomarkers are recently seen as promising tools due to their minimal invasiveness and cost-effectiveness in diagnostic procedures. Furthermore, the application of machine learning (ML) has led the preciseness in the biomarker discovery and the resulting diagnostics. Nevertheless, previous ML-based studies often fail to consider interactions between proteins, which are essential in complex neurodegenerative disorders such as MCI and dementia. Although protein-protein interactions (PPIs) have been employed in network models, these models frequently do not fully capture the diverse properties of PPIs due to their local awareness. This limitation increases the likelihood of overlooking critical components and amplifying the impact of noisy interactions. In this study, we introduce a new graph-based ML model for classifying MCI subtypes, called eXplainable Graph Propagational Network (XGPN). The proposed method extracts the globally interactive effects between proteins by propagating the independent effect of plasma proteins on the PPI network, and thereby, MCI subtypes are predicted by estimation of the risk effect of each protein. Moreover, the process of model training and the outcome of subtype classification are fully explainable due to the simplicity and transparency of XGPN's architecture. The experimental results indicated that the interactive effect between proteins significantly contributed to the distinct differences between MCI subtype groups, resulting in an enhanced classification performance with an average improvement of 10.0 % compared to existing methods, also identifying key biomarkers and their impact on ADCI and SVCI.

Identifying essential proteins is of utmost importance in the field of biomedical research due to their essential functions in cellular activities and their involvement in mechanisms related to diseases. In this research, a novel approach called AttentionEP for predicting essential proteins (EP) is introduced by attention mechanisms. This method leverages both cross-attention and self-attention frameworks, focusing on enhancing prediction accuracy through the integration of features across diverse scales. Spatial characteristics of proteins are obtained from the protein-protein interaction (PPI) network by employing Graph Convolutional Networks (GCN) and Graph Attention Networks (GAT). Following this, Bidirectional Long Short-Term Memory networks (BiLSTM) are employed to derive temporal features from gene expression datasets. Furthermore, spatial characteristics are derived by integrating data on subcellular localization with the application of Deep Neural Networks (DNN). In order to effectively integrate features across multiple scales, initial steps involve the application of self-attention techniques to derive essential insights from each unique data set. Following this, mechanisms involving self-attention and cross-attention are employed to enhance the interaction between diverse information sources. To identify essential proteins, a classifier based on the ResNet architecture is developed. The findings from the experiments indicate that the method introduced here shows superior performance in identifying essential proteins, recording an Area Under the Curve (AUC) value of 0.9433. This approach shows a considerable advantage over established techniques. The findings of this study provide a significant advancement in the comprehension of critical proteins, revealing promising potential for applications in the development of therapeutics and addressing various diseases.

In response to distinct cellular stresses, the p53 exhibits distinct dynamics. These p53 dynamics subsequently control cell fate. However, different stresses can generate the same p53 dynamics with different cell fate outcomes, suggesting that the integration of dynamic information from other pathways is important for cell fate regulation. The interactions between miRNA-125b, p53, and reactive oxygen species (ROS) are significant in the context of cellular stress responses and apoptosis. However, the regulating mechanism of miR-125b with p53 is not fully studied. The dynamics of p53 and its response to the miR-125b regulation are still open questions. In the present study, we try to answer some of these fundamental questions based on basic model built from available experimental reports. The miR-125b-p53 regulatory network is modeled using a set of 11 molecular species variables. The biochemical network of miR-125b-p53, described by 22 reaction channels, is represented by coupled ordinary differential equations (ODEs) using the mass action law of chemical kinetics. These ODEs are solved numerically using the standard fourth-order Runge-Kutta method to analyze the dynamical behavior of the system. The biochemical network model we designed is based on both experimental and theoretical reported data. The p53 dynamics driven by miR-125b exhibit five distinct dynamical states: first and second stable states, first and second dynamical states, and a sustained oscillation state. These different p53 dynamical states may correspond to various cellular conditions. If the stress induced by miR-125b is weak, the system will be weakly activated, favoring a return to normal functioning. However, if the stress is significantly strong, the system will move to an active state. To sustain this active state, which is far from equilibrium with little scope for returning to normal conditions, the system may transition to an apoptotic state by crossing through other intermediate states, as it is unlikely to regain normal functioning. The p53 dynamical states show a multifractal nature, contributed by both short- and long-range correlations. The networks illustrated from these dynamical states follow hierarchical scale-free features, exhibiting an assortative nature with an absence of the centrality-lethality rule. Furthermore, the active dynamical state is generally closer to hierarchical characteristics and is self-organized. Our research study reveals that significant activity of miR-125b on the p53 regulatory network and its dynamics can only be observed when the system is slightly activated by ROS. However, this process does not necessarily require the direct study of ROS activity. These findings elucidate the mechanisms by which cells integrate signaling pathways with distinct temporal activity patterns to encode stress specificity and direct diverse cell fate decisions.

Pecan scab, caused by the fungal pathogen Venturia effusa, is the most devastating disease of pecan (Carya illinoinensis) in the southeastern United States. Resistance to this pathogen is determined by a complex interaction between host genetics and disease pathotype with even field-susceptible cultivars being resistant to most scab isolates. To understand the underlying molecular mechanisms of scab resistance in pecan, we performed a transcriptome analysis of the pecan cultivar, 'Desirable', in response to inoculation with a pathogenic and a non-pathogenic scab isolate at three different time points (24, 48, and 96 hrs. post-inoculation). Differential gene expression and gene ontology enrichment analyses showed contrasting gene expression patterns and pathway enrichment in response to the contrasting isolates with varying pathogenicity. The weighted gene co-expression network analysis of differentially expressed genes detected 11 gene modules. Among them, two modules had significant enrichment of genes involved with defense responses. These genes were particularly upregulated in the resistant reaction at the early stage of fungal infection (24 h) compared to the susceptible reaction. Hub genes in these modules were predominantly related to receptor-like protein kinase activity, signal reception, signal transduction, biosynthesis and transport of plant secondary metabolites, and oxidoreductase activity. Results of this study suggest that the early response of pathogen-related signal transduction and development of cellular barriers against the invading fungus are likely defense mechanisms employed by pecan cultivars against non-virulent scab isolates. The transcriptomic data generated here provide the foundation for identifying candidate resistance genes in pecan against V. effusa and for exploring the molecular mechanisms of disease resistance.

Spatial Analysis of Functional Enrichment (SAFE) is a popular tool for biologists to investigate the functional organization of biological networks via highly intuitive 2D functional maps. To create these maps, SAFE uses Spring embedding to project a given network into a 2D space in which nodes connected in the network are near each other in space. However, many biological networks are scale-free, containing highly connected hub nodes. Because Spring embedding fails to separate hub nodes, it provides uninformative embeddings that resemble a 'hairball'. In addition, Spring embedding only captures direct node connectivity in the network and does not consider higher-order node wiring patterns, which are best captured by graphlets, small, connected, nonisomorphic, induced subgraphs. The scale-free structure of biological networks is hypothesized to stem from an underlying low-dimensional hyperbolic geometry, which novel hyperbolic embedding methods try to uncover. These include coalescent embedding, which projects a network onto a 2D disk.

To better capture the functional organization of scale-free biological networks, whilst also going beyond simple direct connectivity patterns, we introduce Graphlet Coalescent (GraCoal) embedding, which embeds nodes nearby on a disk if they frequently co-occur on a given graphlet together. We use GraCoal to extend SAFE-based network analysis. Through SAFE-enabled enrichment analysis, we show that GraCoal outperforms graphlet-based Spring embedding in capturing the functional organization of the genetic interaction networks of fruit fly, budding yeast, fission yeast and Escherichia coli. We show that depending on the underlying graphlet, GraCoal embeddings capture different topology-function relationships. We show that triangle-based GraCoal embedding captures functional redundancies between paralogs.

https://gitlab.bsc.es/swindels/gracoal_embedding.

Construction of minimal metabolic networks (MMNs) contributes both to our understanding of the origins of metabolism and to the efficiency of biotechnological processes by preventing the diversion of flux away from product formation. We have designed MMNs using a novel in silico synthetic biology pipeline that removes genes encoding enzymes and transporters from genome-scale metabolic models. The resulting minimal gene-set still ensures both viability and high growth rates. The composition of these MMNs has defined a new functional class of genes termed Network Efficiency Determinants (NEDs). These genes, whilst not essential, are very rarely eliminated in constructing an MMN, suggesting that it is difficult for metabolism to be re-routed to obviate the need for such genes. Moreover, the removal of NED genes from an MMN significantly reduces its global efficiency. Bioinformatic analyses of the NED genes have revealed that not only do these genes have more genetic interactions than the bulk of metabolic genes but their protein products also show more protein-protein interactions. In yeast, NED genes are predominantly single-copy and are highly conserved across evolutionarily distant organisms. These features confirm the importance of the NED genes to the metabolic network, including why they are so rarely excluded during minimisation.

The use of graph centrality measures applied to biological networks, such as protein interaction networks, underpins much research into identifying key players within biological processes. This approach however is restricted to dyadic interactions and it is well-known that in many instances interactions are polyadic. In this study we illustrate the merit of using hypergraph centrality applied to a hypernetwork as an alternative. Specifically, we review and propose an extension to a recently introduced node and edge nonlinear hypergraph centrality model which provides mutually dependent node and edge centralities. A Saccharomyces Cerevisiae protein complex hypernetwork is used as an example application with nodes representing proteins and hyperedges representing protein complexes. The resulting rankings of the nodes and edges are considered to see if they provide insight into the essentiality of the proteins and complexes. We find that certain variations of the model predict essentiality more accurately and that the degree-based variation illustrates that the centrality-lethality rule extends to a hypergraph setting. In particular, through exploitation of the models flexibility, we identify small sets of proteins densely populated with essential proteins. One of the key advantages of applying this model to a protein complex hypernetwork is that it also provides a classification method for protein complexes, unlike previous approaches which are only concerned with classifying proteins.

The Katz centrality of a node in a complex network is a measure of the node's importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its "global" character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability [Formula: see text] that a random node in the ensemble of undirected random graphs has centrality [Formula: see text] satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank-[Formula: see text] projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

Accurately identifying essential proteins is vital for drug research and disease diagnosis. Traditional centrality methods and machine learning approaches often face challenges in accurately discerning essential proteins, primarily relying on information derived from protein-protein interaction (PPI) networks. Despite attempts by some researchers to integrate biological data and PPI networks for predicting essential proteins, designing effective integration methods remains a challenge. In response to these challenges, this paper presents the ACDMBI model, specifically designed to overcome the aforementioned issues. ACDMBI is comprised of two key modules: feature extraction and classification. In terms of capturing relevant information, we draw insights from three distinct data sources. Initially, structural features of proteins are extracted from the PPI network through community division. Subsequently, these features are further optimized using Graph Convolutional Networks (GCN) and Graph Attention Networks (GAT). Moving forward, protein features are extracted from gene expression data utilizing Bidirectional Long Short-Term Memory networks (BiLSTM) and a multi-head self-attention mechanism. Finally, protein features are derived by mapping subcellular localization data to a one-dimensional vector and processing it through fully connected layers. In the classification phase, we integrate features extracted from three different data sources, crafting a multi-layer deep neural network (DNN) for protein classification prediction. Experimental results on brewing yeast data showcase the ACDMBI model's superior performance, with AUC reaching 0.9533 and AUPR reaching 0.9153. Ablation experiments further reveal that the effective integration of features from diverse biological information significantly boosts the model's performance.

Despite its medical relevance, there is no commercial vaccine that protects the population at risk from multidrug-resistant (MDR) Klebsiella pneumoniae infections. The availability of massive omic data and novel algorithms may improve antigen selection to develop effective prophylactic strategies. Up to 133 exposed proteins in the core proteomes, between 516 and 8666 genome samples, of the six most relevant MDR clonal groups (CGs) carried conserved B-cell epitopes, suggesting minimized future evasion if utilized for vaccination. Antigens showed a range of epitopicity, functional constraints, and potential side effects. Eleven antigens, including three sugar porins, were represented in all MDR-CGs, constitutively expressed, and showed limited reactivity with gut microbiota. Some of these antigens had important interactomic interactions and may elicit adhesion-neutralizing antibodies. Synergistic bivalent to pentavalent combinations that address expression conditions, interactome location, virulence activities, and clone-specific proteins may overcome the limiting protection of univalent vaccines. The combination of five central antigens accounted for 41% of all non-redundant interacting partners of the antigen dataset. Specific antigen mixtures represented in a few or just one MDR-CG further reduced the chance of microbiota interference. Rational antigen selection schemes facilitate the design of high-coverage and "magic bullet" multivalent vaccines against recalcitrant K. pneumoniae lineages.

Acupuncture is frequently used to manage pediatric cerebral palsy (CP), yet updated evidence is needed to guide future research and clinical practice.

Seven databases were searched from 1994 to 26 June 2023. Randomized controlled trials (RCTs) involving body, scalp, or ear acupuncture for managing CP, excluding acupoint injection, catgut embedding, electro-acupuncture, or laser acupuncture, were included.

Twenty RCTs with 1797 participants were analyzed. Acupuncture groups had better improvements in gross motor function measure (GMFM) scores by 5% (mean difference: 5.93, 95% CI: 3.67-8.19, p < 0.001, I2 = 57%); a 16% higher probability to yield prominent improvement in effectiveness rate (ER) (risk ratio: 1.16, 95% CI: 1.08-1.25, p < 0.001, I2 = 0%); and better outcomes in the Modified Ashworth Scale (MAS) (standardized mean difference [SMD]: 0.3, 95%, CI: 0.11-0.49, p < 0.001, I2 = 0%), the Berg Balance Scale (BBS) (SMD: 2.48; 95% CI: 2.00-2.97, p < 0.001, I2 = 72%) and ADL (SMD: 1.66; 95% CI: 1.23-2.08, p < 0.001, I2 = 91%). Studies with eight core acupoints identified from all ninety-five acupoints had better ER.

Acupuncture, especially using core acupoints, may be effective for managing symptoms in children with CP.

Epistasis is caused by genetic interactions among mutations that affect fitness. To characterize properties and potential mechanisms of epistasis, we engineered eight double mutants that combined mutations from the rho and rpoB genes of Escherichia coli. The two genes encode essential functions for transcription, and the mutations in each gene were chosen because they were beneficial for adaptation to thermal stress (42.2 °C). The double mutants exhibited patterns of fitness epistasis that included diminishing returns epistasis at 42.2 °C, stronger diminishing returns between mutations with larger beneficial effects and both negative and positive (sign) epistasis across environments (20.0 °C and 37.0 °C). By assessing gene expression between single and double mutants, we detected hundreds of genes with gene expression epistasis. Previous work postulated that highly connected hub genes in coexpression networks have low epistasis, but we found the opposite: hub genes had high epistasis values in both coexpression and protein-protein interaction networks. We hypothesized that elevated epistasis in hub genes reflected that they were enriched for targets of Rho termination but that was not the case. Altogether, gene expression and coexpression analyses revealed that thermal adaptation occurred in modules, through modulation of ribonucleotide biosynthetic processes and ribosome assembly, the attenuation of expression in genes related to heat shock and stress responses, and with an overall trend toward restoring gene expression toward the unstressed state.

Root metaxylems are phenotypically diverse structures whose function is particularly important under drought stress. Significant research has dissected the genetic machinery underlying metaxylem phenotypes in dicots, but that of monocots are relatively underexplored. In maize (Zea mays), a robust pipeline integrated a genome-wide association study (GWAS) of root metaxylem phenes under well-watered and water-stress conditions with a gene co-expression network to prioritize the strongest gene candidates. We identified 244 candidate genes by GWAS, of which 103 reside in gene co-expression modules most relevant to xylem development. Several candidate genes may be involved in biosynthetic processes related to the cell wall, hormone signaling, oxidative stress responses, and drought responses. Of those, six gene candidates were detected in multiple root metaxylem phenes in both well-watered and water-stress conditions. We posit that candidate genes that are more essential to network function based on gene co-expression (i.e., hubs or bottlenecks) should be prioritized and classify 33 essential genes for further investigation. Our study demonstrates a new strategy for identifying promising gene candidates and presents several gene candidates that may enhance our understanding of vascular development and responses to drought in cereals.

Forest management by thinning can mitigate the detrimental impact of increasing drought caused by global warming. Growing evidence shows that the soil microbiota can coordinate the dynamic relationship between forest functions and drought intensity, but how they function as a cohesive whole remains elusive. We outline a statistical topology model to chart the roadmap of how each microbe acts and interacts with every other microbe to shape the dynamic changes of microbial communities under forest management. To demonstrate its utility, we analyze a soil microbiota data collected from a two-way longitudinal factorial experiment involving three stand densities and three levels of rainfall over a growing season in artificial plantations of a forest tree - larix (Larix kaempferi). We reconstruct the most sophisticated soil microbiota networks that code maximally informative microbial interactions and trace their dynamic trajectories across time, space, and environmental signals. By integrating GLMY homology theory, we dissect the topological architecture of these so-called omnidirectional networks and identify key microbial interaction pathways that play a pivotal role in mediating the structure and function of soil microbial communities. The statistical topological model described provides a systems tool for studying how microbial community assembly alters its structure, function and evolution under climate change.

Type 2 diabetes (T2D) and Clear-cell renal cell carcinoma (ccRCC) are both complicated diseases which incidence rates gradually increasing. Population based studies show that severity of ccRCC might be associated with T2D. However, so far, no researcher yet investigated about the molecular mechanisms of their association. This study explored T2D and ccRCC causing shared key genes (sKGs) from multiple transcriptomics profiles to investigate their common pathogenetic processes and associated drug molecules. We identified 259 shared differentially expressed genes (sDEGs) that can separate both T2D and ccRCC patients from control samples. Local correlation analysis based on the expressions of sDEGs indicated significant association between T2D and ccRCC. Then ten sDEGs (CDC42, SCARB1, GOT2, CXCL8, FN1, IL1B, JUN, TLR2, TLR4, and VIM) were selected as the sKGs through the protein-protein interaction (PPI) network analysis. These sKGs were found significantly associated with different CpG sites of DNA methylation that might be the cause of ccRCC. The sKGs-set enrichment analysis with Gene Ontology (GO) terms and KEGG pathways revealed some crucial shared molecular functions, biological process, cellular components and KEGG pathways that might be associated with development of both T2D and ccRCC. The regulatory network analysis of sKGs identified six post-transcriptional regulators (hsa-mir-93-5p, hsa-mir-203a-3p, hsa-mir-204-5p, hsa-mir-335-5p, hsa-mir-26b-5p, and hsa-mir-1-3p) and five transcriptional regulators (YY1, FOXL1, FOXC1, NR2F1 and GATA2) of sKGs. Finally, sKGs-guided top-ranked three repurposable drug molecules (Digoxin, Imatinib, and Dovitinib) were recommended as the common treatment for both T2D and ccRCC by molecular docking and ADME/T analysis. Therefore, the results of this study may be useful for diagnosis and therapies of ccRCC patients who are also suffering from T2D.

Alzheimer's disease (AD) is a complex multifaceted disease. Recently approved anti-amyloid monoclonal antibodies slow disease progression by approximately 30%, and combination therapy appears necessary to prevent the onset of AD or produce greater slowing of cognitive and functional decline. Combination therapies may address core features, non-specific co-pathology commonly occurring in patients with AD (e.g., inflammation), or non-AD pathologies that may co-occur with AD (e.g., α-synuclein). Combination therapies may be advanced through co-development of more than one new molecular entity or through add-on strategies including an approved agent plus a new molecular entity. Addressing add-on combination therapy is currently urgent since patients on anti-amyloid monoclonal antibodies may be included in clinical trials for experimental agents. Phase 1 information must be generated for each agent in combination drug development. Phase 2 and Phase 3 of add-on therapies may contrast the new molecular entity, the approved agent as standard of care, and the combination. More complex development programs including standard or modified combinatorial designs are required for co-development of two or more new molecular entities. Biomarkers are markedly affected by anti-amyloid monoclonal antibodies, and these effects must be anticipated in add-on trials. Examining target engagement biomarkers and comparing the magnitude and sequence of biomarker changes in those receiving more than one therapy, compared with those on monotherapy, may be informative. Using network-based medicine approaches, computational strategies may identify rational combinations using disease and drug effect network mapping.

Proteins execute numerous cell functions in concert with one another in protein-protein interactions (PPI). While essential in each cell, such interactions are not identical from cell to cell. Instead, PPI heterogeneity contributes to cellular phenotypic heterogeneity in health and diseases such as cancer. Understanding cellular phenotypic heterogeneity thus requires measurements of properties of PPIs such as abundance, stoichiometry, and kinetics at the single-cell level. Here, we review recent, exciting progress in single-cell PPI measurements. Novel technology in this area is enabled by microscale and microfluidic approaches that control analyte concentration in timescales needed to outpace PPI disassembly kinetics. We describe microscale innovations, needed technical capabilities, and methods poised to be adapted for single-cell analysis in the near future.

The current study focuses on the importance of Protein-Protein Interactions (PPIs) in biological processes and the potential of targeting PPIs as a new treatment strategy for diseases. Specifically, the study explores the cross-links of PPIs network associated with obesity, type 1 diabetes mellitus (T1DM), and cardiac disease (CD), which is an unexplored area of research. The research aimed to understand the role of highly connected proteins in the network and their potential as drug targets. The methodology for this research involves retrieving genes from the NCBI online gene database, intersecting genes among three diseases (type 1 diabetes, obesity, and cardiovascular) using Interactivenn, determining suitable drug molecules using NetworkAnalyst, and performing various bioinformatics analyses such as Generic Protein-Protein Interactions, topological properties analysis, function enrichment analysis in terms of GO, and Kyoto Encyclopedia of Genes and Genomes (KEGG), gene co-expression network, and protein drug as well as protein chemical interaction network. The study focuses on human subjects. The results of this study identified 12 genes [VEGFA (Vascular Endothelial Growth Factor A), IL6 (Interleukin 6), MTHFR (Methylenetetrahydrofolate reductase), NPPB (Natriuretic Peptide B), RAC1 (Rac Family Small GTPase 1), LMNA (Lamin A/C), UGT1A1 (UDP-glucuronosyltransferase family 1 membrane A1), RETN (Resistin), GCG (Glucagon), NPPA (Natriuretic Peptide A), RYR2 (Ryanodine receptor 2), and PRKAG2 (Protein Kinase AMP-Activated Non-Catalytic Subunit Gamma 2)] that were shared across the three diseases and could be used as key proteins for protein-drug/chemical interaction. Additionally, the study provides an in-depth understanding of the complex molecular and biological relationships between the three diseases and the cellular mechanisms that lead to their development. Potentially significant implications for the therapy and management of various disorders are highlighted by the findings of this study by improving treatment efficacy, simplifying treatment regimens, cost-effectiveness, better understanding of the underlying mechanism of these diseases, early diagnosis, and introducing personalized medicine. In conclusion, the current study provides new insights into the cross-links of PPIs network associated with obesity, T1DM, and CD, and highlights the potential of targeting PPIs as a new treatment strategy for these prevalent diseases.

Chronic atrophic gastritis (CAG), characterized by inflammation and erosion of the gastric lining, is a prevalent digestive disorder and considered a precursor to gastric cancer (GC). Coptis chinensis France (CCF) is renowned for its potent heat-clearing, detoxification, and anti-inflammatory properties. Zuojin Pill (ZJP), a classic Chinese medicine primarily composed of CCF, has demonstrated effectiveness in CAG treatment. This study aims to elucidate the potential mechanism of CCF treatment for CAG through a multifaceted approach encompassing network pharmacology, molecular docking, molecular dynamics simulation and experimental verification. The study identified three major active compounds of CCF and elucidated key pathways, such as TNF signaling, PI3K-Akt signaling and p53 signaling. Molecular docking revealed interactions between these active compounds and pivotal targets like PTGS2, TNF, MTOR, and TP53. Additionally, molecular dynamics simulation validated berberine as the primary active compound of CCF, which was further confirmed through experimental verification. This study not only identified berberine as the primary active compound of CCF but also provided valuable insights into the molecular mechanisms underlying CCF's efficacy in treating CAG. Furthermore, it offers a reference for refining therapeutic strategies for CAG management.

Cardiovascular disease is a chronic inflammatory disease with high mortality rates. TNF-alpha is pro-inflammatory and associated with the disease, but current medications have adverse effects. Therefore, efficient inhibitors are urgently needed as alternatives. This study represents a structural-activity relationship investigation of TNF-alpha, curated from the ChEMBL database. Exploratory data analysis was performed to visualize the physicochemical properties of different bioactivity groups. The extracted molecules were subjected to PubChem and SubStructure fingerprints, and a QSAR-based Random Forest (QSAR-RF) model was generated using the WEKA tool. The QSAR random Forest model was built based on the SubStructure fingerprint with a correlation coefficient of 0.992 and 0.716 as the respective tenfold cross-validation scores. The variance important plot (VIP) method was used to extract the important features for TNF-alpha inhibition. The Substructure-based QSAR-RF (SS-QSAR-RF) model was validated using molecules from PubChem and ZINC databases. The generated model also predicts the pIC50 value of the molecules selected from the docking study followed by molecular dynamic simulation with the time step of 100 ns. Through virtual reverse pharmacology, we determined the main drug targets from the top four hit compounds obtained via molecular docking study. Our analysis included an integrated bioinformatics approach to pinpoint crucial targets like EGRF, HSP900A1, STAT3, PSEN1, AKT1, and MDM2. Further, GO and KEGG pathways analysis identified relevant cardiovascular disease-related pathways for the hub gene involved. However, this study provides valuable insights, it is important to note that it lacks experimental application. Future research may benefit from conducting in-vitro and in-vivo studies.

Plasma protein biomarkers have been considered promising tools for diagnosing dementia subtypes due to their low variability, cost-effectiveness, and minimal invasiveness in diagnostic procedures. Machine learning (ML) methods have been applied to enhance accuracy of the biomarker discovery. However, previous ML-based studies often overlook interactions between proteins, which are crucial in complex disorders like dementia. While protein-protein interactions (PPIs) have been used in network models, these models often fail to fully capture the diverse properties of PPIs due to their local awareness. This drawback increases the chance of neglecting critical components and magnifying the impact of noisy interactions. In this study, we propose a novel graph-based ML model for dementia subtype diagnosis, the graph propagational network (GPN). By propagating the independent effect of plasma proteins on PPI network, the GPN extracts the globally interactive effects between proteins. Experimental results showed that the interactive effect between proteins yielded to further clarify the differences between dementia subtype groups and contributed to the performance improvement where the GPN outperformed existing methods by 10.4% on average.

Genome-wide association studies (GWASs) have uncovered over 75 genomic loci associated with risk for late-onset Alzheimer's disease (LOAD), but identification of the underlying causal genes remains challenging. Studies of induced pluripotent stem cell (iPSC)-derived neurons from LOAD patients have demonstrated the existence of neuronal cell-intrinsic functional defects. Here, we searched for genetic contributions to neuronal dysfunction in LOAD using an integrative systems approach that incorporated multi-evidence-based gene mapping and network-analysis-based prioritization. A systematic perturbation screening of candidate risk genes in Caenorhabditis elegans (C. elegans) revealed that neuronal knockdown of the LOAD risk gene orthologs vha-10 (ATP6V1G2), cmd-1 (CALM3), amph-1 (BIN1), ephx-1 (NGEF), and pho-5 (ACP2) alters short-/intermediate-term memory function, the cognitive domain affected earliest during LOAD progression. These results highlight the impact of LOAD risk genes on evolutionarily conserved memory function, as mediated through neuronal endosomal dysfunction, and identify new targets for further mechanistic interrogation.

The selective pressure of pathogen-host symbiosis drives adaptations. How these interactions shape the metabolism of pathogens is largely unknown. Here, we use comparative genomics to systematically analyze the metabolic networks of oomycetes, a diverse group of eukaryotes that includes saprotrophs as well as animal and plant pathogens, with the latter causing devastating diseases with significant economic and/or ecological impacts. In our analyses of 44 oomycete species, we uncover considerable variation in metabolism that can be linked to lifestyle differences. Comparisons of metabolic gene content reveal that plant pathogenic oomycetes have a bipartite metabolism consisting of a conserved core and an accessory set. The accessory set can be associated with the degradation of defense compounds produced by plants when challenged by pathogens. Obligate biotrophic oomycetes have smaller metabolic networks, and taxonomically distantly related biotrophic lineages display convergent evolution by repeated gene losses in both the conserved as well as the accessory set of metabolisms. When investigating to what extent the metabolic networks in obligate biotrophs differ from those in hemibiotrophic plant pathogens, we observe that the losses of metabolic enzymes in obligate biotrophs are not random and that gene losses predominantly influence the terminal branches of the metabolic networks. Our analyses represent the first metabolism-focused comparison of oomycetes at this scale and will contribute to a better understanding of the evolution of oomycete metabolism in relation to lifestyle adaptation. Numerous oomycete species are devastating plant pathogens that cause major damage in crops and natural ecosystems. Their interactions with hosts are shaped by strong selection, but how selection affects adaptation of the primary metabolism to a pathogenic lifestyle is not yet well established. By pan-genome and metabolic network analyses of distantly related oomycete pathogens and their nonpathogenic relatives, we reveal considerable lifestyle- and lineage-specific adaptations. This study contributes to a better understanding of metabolic adaptations in pathogenic oomycetes in relation to lifestyle, host, and environment, and the findings will help in pinpointing potential targets for disease control. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Spinal cord injury (SCI) is a severe neurological complication following spinal fracture, which has long posed a challenge for clinicians. Microglia play a dual role in the pathophysiological process after SCI, both beneficial and detrimental. The underlying mechanisms of microglial actions following SCI require further exploration. The present study combined three different machine learning algorithms, namely weighted gene co-expression network analysis, random forest analysis and least absolute shrinkage and selection operator analysis, to screen for differentially expressed genes in the GSE96055 microglia dataset after SCI. It then used protein-protein interaction networks and gene set enrichment analysis with single genes to investigate the key genes and signaling pathways involved in microglial function following SCI. The results indicated that microglia not only participate in neuroinflammation but also serve a significant role in the clearance mechanism of apoptotic cells following SCI. Notably, bioinformatics analysis and lipopolysaccharide + UNC569 (a MerTK-specific inhibitor) stimulation of BV2 cell experiments showed that the expression levels of Anxa2, Myo1e and Spp1 in microglia were significantly upregulated following SCI, thus potentially involved in regulating the clearance mechanism of apoptotic cells. The present study suggested that Anxa2, Myo1e and Spp1 may serve as potential targets for the future treatment of SCI and provided a theoretical basis for the development of new methods and drugs for treating SCI.

The metastasis suppressor, N-myc downstream regulated gene-1 (NDRG1), inhibits pro-oncogenic signaling in pancreatic cancer (PC). This investigation dissected a novel mechanism induced by NDRG1 on WNT/β-catenin signaling in multiple PC cell types. NDRG1 overexpression decreased β-catenin and downregulated glycogen synthase kinase-3β (GSK-3β) protein levels and its activation. However, β-catenin phosphorylation at Ser33, Ser37, and Thr41 are classically induced by GSK-3β was significantly increased after NDRG1 overexpression, suggesting a GSK-3β-independent mechanism. Intriguingly, NDRG1 overexpression upregulated protein kinase Cα (PKCα), with PKCα silencing preventing β-catenin phosphorylation at Ser33, Ser37, and Thr41, and decreasing β-catenin expression. Further, NDRG1 and PKCα were demonstrated to associate, with PKCα stabilization occurring after NDRG1 overexpression. PKCα half-life increased from 1.5 ± 0.8 h (3) in control cells to 11.0 ± 2.5 h (3) after NDRG1 overexpression. Thus, NDRG1 overexpression leads to the association of NDRG1 with PKCα and PKCα stabilization, resulting in β-catenin phosphorylation at Ser33, Ser37, and Thr41. The association between PKCα, NDRG1, and β-catenin was identified, with the formation of a potential metabolon that promotes the latter β-catenin phosphorylation. This anti-oncogenic activity of NDRG1 was multi-modal, with the above mechanism accompanied by the downregulation of the nucleo-cytoplasmic shuttling protein, p21-activated kinase 4 (PAK4), which is involved in β-catenin nuclear translocation, inhibition of AKT phosphorylation (Ser473), and decreased β-catenin phosphorylation at Ser552 that suppresses its transcriptional activity. These mechanisms of NDRG1 activity are important to dissect to understand the marked anti-cancer efficacy of NDRG1-inducing thiosemicarbazones that upregulate PKCα and inhibit WNT signaling.

Bow-tie architecture is a layered network structure that has a narrow middle layer with multiple inputs and outputs. Such structures are widely seen in the molecular networks in cells, suggesting that a universal evolutionary mechanism underlies the emergence of bow-tie architecture. The previous theoretical studies have implemented evolutionary simulations of the feedforward network to satisfy a given input-output goal and proposed that the bow-tie architecture emerges when the ideal input-output relation is given as a rank-deficient matrix with mutations in network link intensities in a multiplicative manner. Here, we report that the bow-tie network inevitably appears when the link intensities representing molecular interactions are small at the initial condition of the evolutionary simulation, regardless of the rank of the goal matrix. Our dynamical system analysis clarifies the mechanisms underlying the emergence of the bow-tie structure. Further, we demonstrate that the increase in the input-output matrix reduces the width of the middle layer, resulting in the emergence of bow-tie architecture, even when evolution starts from large link intensities. Our data suggest that bow-tie architecture emerges as a side effect of evolution rather than as a result of evolutionary adaptation.