Using graph theory to analyze biological networks

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Quantitative Proteomic Analysis of Serum Reveals MST1 as a Potential Candidate Biomarker in Spontaneously Diabetic Cynomolgus Monkeys

The prevalence of type 2 diabetes (T2DM) is increasing globally, creating essential demands for T2DM animal models for the study of disease pathogenesis, prevention, and therapy. A non-human primate model such as cynomolgus monkeys can develop T2DM spontaneously in an age-dependent way similar to humans. In this study, a data-independent acquisition-based quantitative proteomics strategy was employed to investigate the serum proteomic profiles of spontaneously diabetic cynomolgus monkeys compared with healthy controls. The results revealed significant differences in protein abundances. A total of 95 differentially expressed proteins (DEPs) were quantitatively identified in the current study, among which 31 and 64 proteins were significantly upregulated and downregulated, respectively. Bioinformatic analysis revealed that carbohydrate digestion and absorption was the top enriched pathway by the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Protein-protein interaction network analysis demonstrated that MST1 was identified as the most connected protein in the network and could be considered as the hub protein. MST1 was significantly and inversely associated with FSG and HbA1c. Furthermore, recent lines of evidence also indicate that MST1 acts as a crucial regulator in regulating hepatic gluconeogenesis to maintain metabolic homeostasis while simultaneously suppressing the inflammatory processes. In conclusion, our study provides novel insights into serum proteome changes in spontaneously diabetic cynomolgus monkeys and points out that the dysregulation of several DEPs may play an important role in the pathogenesis of T2DM.

Unpacking COVID-19 Systems Biology in Lung and Whole Blood with Transcriptomics and miRNA Regulators

COVID-19 is a systemic disease affecting tissues and organs, including and beyond the lung. Apart from the current pandemic context, we also have vastly inadequate knowledge of consequences of repeated exposures to SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the virus causing COVID-19, in multiple organ systems and the whole organism scales when the disease evolves from a pandemic to an endemic state. This calls for a systems biology and systems medicine approach and unpacking the effects of COVID-19 in lung as well as other tissues. We report here original findings from transcriptomics analyses and differentially expressed genes (DEGs) in lung samples from 60 patients and 27 healthy controls, and in whole blood samples from 255 patients and 103 healthy individuals. A total of 11 datasets with RNA-seq transcriptomic data were obtained from the Gene Expression Omnibus and the European Nucleotide Archive. The identified DEGs were used to construct protein interaction and functional networks and to identify related pathways and miRNAs. We found 35 DEGs common between lung and the whole blood, and importantly, 2 novel genes, namely CYP1B1 and TNFAIP6, which have not been previously implicated with COVID-19. We also identified four novel miRNA potential regulators, hsa-mir-192-5p, hsa-mir-221-3p, hsa-mir-4756-3p, and hsa-mir-10a-5p, implicated in lung or other diseases induced by coronaviruses. In summary, these findings offer new molecular leads and insights to unpack COVID-19 systems biology in a whole organism context and might inform future antiviral drug, diagnostics, and vaccine discovery efforts.

Nucleotide-based genetic networks: Methods and applications

Genomic variations have been acclaimed as among the key players in understanding the biological mechanisms behind migration, evolution, and adaptation to extreme conditions. Due to stochastic evolutionary forces, the frequency of polymorphisms is affected by changes in the frequency of nearby polymorphisms in the same DNA sample, making them connected in terms of evolution. This article presents all the ingredients to understand the cumulative effects and complex behaviors of genetic variations in the human mitochondrial genome by analyzing co-occurrence networks of nucleotides, and shows key results obtained from such analyses. The article emphasizes recent investigations of these co-occurrence networks, describing the role of interactions between nucleotides in fundamental processes of human migration and viral evolution. The corresponding co-mutation-based genetic networks revealed genetic signatures of human adaptation in extreme environments. This article provides the methods of constructing such networks in detail, along with their graph-theoretical properties, and applications of the genomic networks in understanding the role of nucleotide co-evolution in evolution of the whole genome.

A comprehensive framework for preference aggregation Parrondo's paradox

Individuals can make choices for themselves that are beneficial or detrimental to the entire group. Consider two losing choices that some individuals have to make on behalf of the group. Is it possible that the losing choices combine to give a winning outcome? We show that it is possible through a variant of Parrondo's paradox-the preference aggregation Parrondo's paradox (PAPP). This new variant of Parrondo's paradox makes use of an aggregate rule that combines with a decision-making heuristic that can be applied to individuals or parts of the social group. The aim of this work is to discuss this PAPP framework and exemplify it on a social network. This work enhances existing research by constructing a feedback loop that allows individuals in the social network to adapt its behavior according to the outcome of the Parrondo's games played.

Smart ECG Biosensor Design with an Improved ANN Performance Based on the Taguchi Optimizer

This paper aims to design a smart biosensor to predict electrocardiogram (ECG) signals in a specific auscultation site from other ECG signals measured from other measurement sites. The proposed design is based on a hybrid architecture using the Artificial Neural Networks (ANNs) model and Taguchi optimizer to avoid the ANN issues related to hyperparameters and to improve its accuracy. The proposed approach aims to optimize the number and type of inputs to be considered for the ANN model. Indeed, different combinations are considered in order to find the optimal input combination for the best prediction quality. By identifying the factors that influence a model’s prediction and their degree of importance via the modified Taguchi optimizer, the developed biosensor improves the prediction accuracy of ECG signals collected from different auscultation sites compared to the ANN-based biosensor. Based on an actual database, the simulation results show that this improvement is significant; it can reach more than 94% accuracy.

Microbiome structure in biofilms from a volcanic island in Maritime Antarctica investigated by genome-centric metagenomics and metatranscriptomics

Antarctica is the coldest and driest continent on Earth, characterized by polyextreme environmental conditions, where species adapted form complex networks of interactions. Microbial communities growing in these harsh environments can form biofilms that help the associated species to survive and thrive. A rich body of knowledge describes environmental biofilm communities; however, most studies have focused on dominant community members rather than functional complexity and metabolic potential. To overcome these limitations, the present study used genome-centric metagenomics to describe two biofilm samples subjected to different temperature collected in Deception Island, Maritime Antarctica. The results unraveled a complex biofilm microbiome represented by 180 metagenome-assembled genomes. The potential metabolic interactions were investigated using metabolic flux balance analysis and revealed that purple bacteria are the community members with the highest correlations with other bacteria. Due to their predicted mixotrophic behavior, they may play a crucial role in the microbiome, likely supporting the heterotrophic species in biofilms. Metatranscriptomics results revealed that the chaperone system and proteins counteracting ROS and toxic compounds have a major role in maintaining bacterial cell homeostasis in sediments of volcanic origin.

Inferring the underlying multivariate structure from bivariate networks with highly correlated nodes

Complex systems are often described mathematically as networks. Inferring the actual interactions from observed dynamics of the nodes of the networks is a challenging inverse task. It is crucial to distinguish direct and indirect interactions to allow for a robust identification of the underlying network. If strong and weak links are simultaneously present in the observed network, typical multivariate approaches to address this challenge fail. By means of correlation and partial correlation, we illustrate the challenges that arise and demonstrate how to overcome these. The challenge of strong and weak links translates into ill-conditioned matrices that need to be inverted to obtain the partial correlations, and therefore the correct network topology. Our novel procedure enables robust identification of multivariate network topologies in the presence of highly correlated processes. In applications, this is crucial to avoid erroneous conclusions about network structures and characteristics. Our novel approach applies to other types of interaction measures between processes in a network.

An integrated approach to understand fluid shear stress-driven and reactive oxygen species-mediated metastasis of colon adenocarcinoma through mRNA-miRNA-lncRNA-circRNA networks

Development of colon adenocarcinoma (COAD) metastasis involves several mediators including fluid shear stress (FSS), intracellular ROS levels, and non-coding RNAs. In our present study, we identified and investigated the role of regulatory non-coding RNA molecules specifically involved in COAD metastasis and their association with FSS and ROS. Interactions between the mRNAs associated with FSS and ROS, the corresponding microRNAs (miRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) in COAD metastasis were used to generate the mRNA-miRNA-lncRNA-circRNA network. Experimental validation of the identified RNA hubs using quantitative real-time PCR demonstrated a direct effect of the FSS on their expression levels in cancer cells. FSS resulted in the downregulation of HMGA1 and RAN, as well as the upregulation of HSP90AA1, PMAIP1 and BIRC5. Application of shear stress also led to downregulation of hsa-miR-26b-5p and hsa-miR-34a-5p levels in HCT116 cells. Further, functional enrichment and survival analysis of the significant miRNAs, as well as the OncoPrint and the survival analyses of the selected mRNAs were performed. Subsequently, their functional role was also corroborated with existing literature. Ten significant miRNA hubs were identified, out of which hsa-miR-17-5p and hsa-miR-20a-5p were found to interact with lncRNA (CCAT2) while hsa-miR-335 was found to interact with four circRNAs. Fifteen significant miRNAs were identified in 10 different modules suggesting their importance in FSS and ROS-mediated COAD metastasis. Finally, 10 miRNAs and 3 mRNAs associated with FSS and/or ROS were identified as significant overall survival markers; 33 mRNAs were also identified as metastasis-free survival markers whereas 15 mRNAs showed > 10% gene alterations in TCGA-COAD data and may serve as promising therapeutic biomarkers in the COAD metastasis.

Limbic covariance network alterations in patients with transient global amnesia

BACKGROUND

We compared limbic structure volumes and graph theory parameters of the limbic covariance network between patients with transient global amnesia (TGA) and healthy controls, and between patients with single and recurrent TGA events.

METHODS

We retrospectively enrolled 122 patients with TGA (single event, n = 107; recurrent events, n = 15) and 50 healthy controls who underwent three-dimensional T1-weighted MRI imaging of the brain. Volumetric analysis of the subcortical limbic structures, including the hippocampus, amygdala, thalamus, mammillary body, hypothalamus, basal forebrain, septal nuclei, fornix, and nucleus accumbens, was performed. We examined the limbic covariance network using a graph theory.

RESULTS

Limbic structure volumes did not differ between patients with TGA and healthy controls, and between patients with a single event and those with recurrent events. However, the radius of the limbic covariance network was significantly greater in patients with TGA than in healthy controls (6.595 vs. 4.564, p = 0.040). Furthermore, the radius, diameter, eccentricity, and characteristics path length were greater (4.066 vs. 2.000, p = 0.009; 7.062 vs. 3.645, p = 0.029; 5.633 vs. 2.774, p = 0.013; 3.373 vs. 1.688, p = 0.004; respectively), whereas the average strength, global efficiency, local efficiency, mean clustering coefficient, transitivity, and small-worldness index were lower (5.595 vs. 10.831, p = 0.004; 0.350 vs. 0.642, p = 0.002; 0.531 vs. 1.724, p = 0.004; 0.304 vs. 0.624, p = 0.006; 0.456 vs. 0.935, p = 0.003; 0.913 vs. 0.993, p = 0.017; respectively), in patients with recurrent events than in those with a single event.

CONCLUSION

The limbic covariance network shows significant alterations in patients with TGA, as well as differences between patients with recurrent events and those with a single event. These findings suggest that changes in the limbic covariance network could be related to the pathogenesis of TGA.

Approaches in Gene Coexpression Analysis in Eukaryotes

Simple Summary

Genes whose expression levels rise and fall similarly in a large set of samples, may be considered coexpressed. Gene coexpression analysis refers to the en masse discovery of coexpressed genes from a large variety of transcriptomic experiments. The type of biological networks that studies gene coexpression, known as Gene Coexpression Networks, consist of an undirected graph depicting genes and their coexpression relationships. Coexpressed genes are clustered in smaller subnetworks, the predominant biological roles of which can be determined through enrichment analysis. By studying well-annotated gene partners, the attribution of new roles to genes of unknown function or assumption for participation in common metabolic pathways can be achieved, through a guilt-by-association approach. In this review, we present key issues in gene coexpression analysis, as well as the most popular tools that perform it.

Abstract

Gene coexpression analysis constitutes a widely used practice for gene partner identification and gene function prediction, consisting of many intricate procedures. The analysis begins with the collection of primary transcriptomic data and their preprocessing, continues with the calculation of the similarity between genes based on their expression values in the selected sample dataset and results in the construction and visualisation of a gene coexpression network (GCN) and its evaluation using biological term enrichment analysis. As gene coexpression analysis has been studied extensively, we present most parts of the methodology in a clear manner and the reasoning behind the selection of some of the techniques. In this review, we offer a comprehensive and comprehensible account of the steps required for performing a complete gene coexpression analysis in eukaryotic organisms. We comment on the use of RNA-Seq vs. microarrays, as well as the best practices for GCN construction. Furthermore, we recount the most popular webtools and standalone applications performing gene coexpression analysis, with details on their methods, features and outputs.

Mapping DEHP to the adverse outcome pathway network for human female reproductive toxicity

Adverse outcome pathways (AOPs) and AOP networks are tools for mechanistic presentation of toxicological effects across different levels of biological organization. These tools are used to better understand how chemicals impact human health. In this study, a four-step workflow was used to derive the AOP network of human female reproductive toxicity (HFRT-AOP) from five AOPs available in the AOP-Wiki and ten AOPs obtained from the literature. Standard network analysis identified key events (KEs) that are point of convergence and divergence, upstream and downstream KEs, and bottlenecks across the network. To map di-(2-ethylhexyl) phthalate (DEHP) to the HFRT-AOP network, we extracted DEHP target genes and proteins from the Comparative Toxicogenomic and the CompTox Chemicals Dashboard databases. Enriched GO terms analysis was used to identify relevant biological processes in the ovary that are DEHP targets, whereas screening of scientific literature was performed manually and automatically using AOP-helpFinder. We combined this information to map DEHP to HFRT-AOP network to provide insight on the KEs and system-level perturbations caused by this endocrine disruptor and the emergent paths. This approach can enable better understanding of the toxic mechanism of DEHP-induced human female reproductive toxicity and reveal potential novel DEHP female reproductive targets for experimental studies.

A network biology approach to identify crucial host targets for COVID-19

The ongoing global pandemic of COVID-19, caused by SARS-CoV-2 has killed more than 5.9 million individuals out of ∼43 million confirmed infections. At present, several parts of the world are encountering the 3rd wave. Mass vaccination has been started in several countries but they are less likely to be broadly available for the current pandemic, repurposing of the existing drugs has drawn highest attention for an immediate solution. A recent publication has mapped the physical interactions of SARS-CoV-2 and human proteins by affinity-purification mass spectrometry (AP-MS) and identified 332 high-confidence SARS-CoV-2-human protein-protein interactions (PPIs). Here, we taken a network biology approach and constructed a human protein-protein interaction network (PPIN) with the above SARS-CoV-2 targeted proteins. We utilized a combination of essential network centrality measures and functional properties of the human proteins to identify the critical human targets of SARS-CoV-2. Four human proteins, namely PRKACA, RHOA, CDK5RAP2, and CEP250 have emerged as the best therapeutic targets, of which PRKACA and CEP250 were also found by another group as potential candidates for drug targets in COVID-19. We further found candidate drugs/compounds, such as guanosine triphosphate, remdesivir, adenosine monophosphate, MgATP, and H-89 dihydrochloride that bind the target human proteins. The urgency to prevent the spread of infection and the death of diseased individuals has prompted the search for agents from the pool of approved drugs to repurpose them for COVID-19. Our results indicate that host targeting therapy with the repurposed drugs may be a useful strategy for the treatment of SARS-CoV-2 infection.

Enhancing Data Integration, Interoperability, and Reuse to Address Complex and Emerging Environmental Health Problems

Environmental health sciences (EHS) span many diverse disciplines. Within the EHS community, the National Institute of Environmental Health Sciences Superfund Research Program (SRP) funds multidisciplinary research aimed to address pressing and complex issues on how people are exposed to hazardous substances and their related health consequences with the goal of identifying strategies to reduce exposures and protect human health. While disentangling the interrelationships that contribute to environmental exposures and their effects on human health over the course of life remains difficult, advances in data science and data sharing offer a path forward to explore data across disciplines to reveal new insights. Multidisciplinary SRP-funded teams are well-positioned to examine how to best integrate EHS data across diverse research domains to address multifaceted environmental health problems. As such, SRP supported collaborative research projects designed to foster and enhance the interoperability and reuse of diverse and complex data streams. This perspective synthesizes those experiences as a landscape view of the challenges identified while working to increase the FAIR-ness (Findable, Accessible, Interoperable, and Reusable) of EHS data and opportunities to address them.

Drug repositioning in non-small cell lung cancer (NSCLC) using gene co-expression and drug–gene interaction networks analysis

Lung cancer is the most common cancer in men and women. This cancer is divided into two main types, namely non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC). Around 85 to 90 percent of lung cancers are NSCLC. Repositioning potent candidate drugs in NSCLC treatment is one of the important topics in cancer studies. Drug repositioning (DR) or drug repurposing is a method for identifying new therapeutic uses of existing drugs. The current study applies a computational drug repositioning method to identify candidate drugs to treat NSCLC patients. To this end, at first, the transcriptomics profile of NSCLC and healthy (control) samples was obtained from the GEO database with the accession number GSE21933. Then, the gene co-expression network was reconstructed for NSCLC samples using the WGCNA, and two significant purple and magenta gene modules were extracted. Next, a list of transcription factor genes that regulate purple and magenta modules' genes was extracted from the TRRUST V2.0 online database, and the TF–TG (transcription factors–target genes) network was drawn. Afterward, a list of drugs targeting TF–TG genes was obtained from the DGIdb V4.0 database, and two drug–gene interaction networks, including drug-TG and drug-TF, were drawn. After analyzing gene co-expression TF–TG, and drug–gene interaction networks, 16 drugs were selected as potent candidates for NSCLC treatment. Out of 16 selected drugs, nine drugs, namely Methotrexate, Olanzapine, Haloperidol, Fluorouracil, Nifedipine, Paclitaxel, Verapamil, Dexamethasone, and Docetaxel, were chosen from the drug-TG sub-network. In addition, nine drugs, including Cisplatin, Daunorubicin, Dexamethasone, Methotrexate, Hydrocortisone, Doxorubicin, Azacitidine, Vorinostat, and Doxorubicin Hydrochloride, were selected from the drug-TF sub-network. Methotrexate and Dexamethasone are common in drug-TG and drug-TF sub-networks. In conclusion, this study proposed 16 drugs as potent candidates for NSCLC treatment through analyzing gene co-expression, TF–TG, and drug–gene interaction networks.

Constraint-Based, Score-Based and Hybrid Algorithms to Construct Bayesian Gene Networks in the Bovine Transcriptome

Simple Summary

In this study, we investigated and compared six different Bayesian network algorithms from three different categories to identify hub genes critical to gene expression networks activated in response to progesterone in the bovine uterus. We observed many common hub genes identified between constraint-based algorithms (CBAs) and hybrid algorithms (HAs), while it appeared that score-based algorithm (SBA) methods led to more accurate and relevant predictions of core genes. The results revealed that the identification of hub genes was affected by the type of network reconstruction and by the subsequently used topological parameters. Two identified genes known to have roles during pregnancy are ISG15 and DGAT2. The identified hub genes are associated with biological processes such as amino acid metabolism, hormonal signaling pathways and the immune system. Our analysis revealed a role for miRNAs in the regulation of this system. The biological and physiological roles (enzymatic and hormonal effects) of unannotated identified hub genes should be functionally validated by further studies.

Abstract

Bayesian gene networks are powerful for modelling causal relationships and incorporating prior knowledge for making inferences about relationships. We used three algorithms to construct Bayesian gene networks around genes expressed in the bovine uterus and compared the efficacies of the algorithms. Dataset GSE33030 from the Gene Expression Omnibus (GEO) repository was analyzed using different algorithms for hub gene expression due to the effect of progesterone on bovine endometrial tissue following conception. Six different algorithms (grow-shrink, max-min parent children, tabu search, hill-climbing, max-min hill-climbing and restricted maximum) were compared in three higher categories, including constraint-based, score-based and hybrid algorithms. Gene network parameters were estimated using the bnlearn bundle, which is a Bayesian network structure learning toolbox implemented in R. The results obtained indicated the tabu search algorithm identified the highest degree between genes (390), Markov blankets (25.64), neighborhood sizes (8.76) and branching factors (4.38). The results showed that the highest number of shared hub genes (e.g., proline dehydrogenase 1 (PRODH), Sam-pointed domain containing Ets transcription factor (SPDEF), monocyte-to-macrophage differentiation associated 2 (MMD2), semaphorin 3E (SEMA3E), solute carrier family 27 member 6 (SLC27A6) and actin gamma 2 (ACTG2)) was seen between the hybrid and the constraint-based algorithms, and these genes could be recommended as central to the GSE33030 data series. Functional annotation of the hub genes in uterine tissue during progesterone treatment in the pregnancy period showed that the predicted hub genes were involved in extracellular pathways, lipid and protein metabolism, protein structure and post-translational processes. The identified hub genes obtained by the score-based algorithms had a role in 2-arachidonoylglycerol and enzyme modulation. In conclusion, different algorithms and subsequent topological parameters were used to identify hub genes to better illuminate pathways acting in response to progesterone treatment in the bovine uterus, which should help with our understanding of gene regulatory networks in complex trait expression.

Elastic network modeling of cellular networks unveils sensor and effector genes that control information flow

The high-level organization of the cell is embedded in indirect relationships that connect distinct cellular processes. Existing computational approaches for detecting indirect relationships between genes typically consist of propagating abstract information through network representations of the cell. However, the selection of genes to serve as the source of propagation is inherently biased by prior knowledge. Here, we sought to derive an unbiased view of the high-level organization of the cell by identifying the genes that propagate and receive information most effectively in the cell, and the indirect relationships between these genes. To this aim, we adapted a perturbation-response scanning strategy initially developed for identifying allosteric interactions within proteins. We deployed this strategy onto an elastic network model of the yeast genetic interaction profile similarity network. This network revealed a superior propensity for information propagation relative to simulated networks with similar topology. Perturbation-response scanning identified the major distributors and receivers of information in the network, named effector and sensor genes, respectively. Effectors formed dense clusters centrally integrated into the network, whereas sensors formed loosely connected antenna-shaped clusters and contained genes with previously characterized involvement in signal transduction. We propose that indirect relationships between effector and sensor clusters represent major paths of information flow between distinct cellular processes. Genetic similarity networks for fission yeast and human displayed similarly strong propensities for information propagation and clusters of effector and sensor genes, suggesting that the global architecture enabling indirect relationships is evolutionarily conserved across species. Our results demonstrate that elastic network modeling of cellular networks constitutes a promising strategy to probe the high-level organization and cooperativity in the cell.

Emerging landscape of molecular interaction networks: Opportunities, challenges and prospects

Network biology finds application in interpreting molecular interaction networks and providing insightful inferences using graph theoretical analysis of biological systems. The integration of computational bio-modelling approaches with different hybrid network-based techniques provides additional information about the behaviour of complex systems. With increasing advances in high-throughput technologies in biological research, attempts have been made to incorporate this information into network structures, which has led to a continuous update of network biology approaches over time. The newly minted centrality measures accommodate the details of omics data and regulatory network structure information. The unification of graph network properties with classical mathematical and computational modelling approaches and technologically advanced approaches like machine-learning- and artificial intelligence-based algorithms leverages the potential application of these techniques. These computational advances prove beneficial and serve various applications such as essential gene prediction, identification of drug–disease interaction and gene prioritization. Hence, in this review, we have provided a comprehensive overview of the emerging landscape of molecular interaction networks using graph theoretical approaches. With the aim to provide information on the wide range of applications of network biology approaches in understanding the interaction and regulation of genes, proteins, enzymes and metabolites at different molecular levels, we have reviewed the methods that utilize network topological properties, emerging hybrid network-based approaches and applications that integrate machine learning techniques to analyse molecular interaction networks. Further, we have discussed the applications of these approaches in biomedical research with a note on future prospects.

Comparing protein-protein interaction networks of SARS-CoV-2 and (H1N1) influenza using topological features

SARS-CoV-2 pandemic first emerged in late 2019 in China. It has since infected more than 298 million individuals and caused over 5 million deaths globally. The identification of essential proteins in a protein–protein interaction network (PPIN) is not only crucial in understanding the process of cellular life but also useful in drug discovery. There are many centrality measures to detect influential nodes in complex networks. Since SARS-CoV-2 and (H1N1) influenza PPINs pose 553 common human proteins. Analyzing influential proteins and comparing these networks together can be an effective step in helping biologists for drug-target prediction. We used 21 centrality measures on SARS-CoV-2 and (H1N1) influenza PPINs to identify essential proteins. We applied principal component analysis and unsupervised machine learning methods to reveal the most informative measures. Appealingly, some measures had a high level of contribution in comparison to others in both PPINs, namely Decay, Residual closeness, Markov, Degree, closeness (Latora), Barycenter, Closeness (Freeman), and Lin centralities. We also investigated some graph theory-based properties like the power law, exponential distribution, and robustness. Both PPINs tended to properties of scale-free networks that expose their nature of heterogeneity. Dimensionality reduction and unsupervised learning methods were so effective to uncover appropriate centrality measures.

Darling: A Web Application for Detecting Disease-Related Biomedical Entity Associations with Literature Mining

Finding, exploring and filtering frequent sentence-based associations between a disease and a biomedical entity, co-mentioned in disease-related PubMed literature, is a challenge, as the volume of publications increases. Darling is a web application, which utilizes Name Entity Recognition to identify human-related biomedical terms in PubMed articles, mentioned in OMIM, DisGeNET and Human Phenotype Ontology (HPO) disease records, and generates an interactive biomedical entity association network. Nodes in this network represent genes, proteins, chemicals, functions, tissues, diseases, environments and phenotypes. Users can search by identifiers, terms/entities or free text and explore the relevant abstracts in an annotated format.

Zooplankton network conditioned by turbidity gradient in small anthropogenic reservoirs

Water turbidity can significantly influence interspecific interactions in aquatic ecosystems. We tested the hypothesis that the turbidity gradient significantly differentiates the dynamics, significance and type of relationships in the structure of zooplankton communities colonizing mine pit reservoirs. The interactions between zooplankton species were evaluated by network graph analysis for three water turbidity classes: high turbidity (HT), moderate turbidity (MT) and low turbidity (LT). The HT network was most cohesive, and it was controlled by taxa grazing on various food sources within one ecological niche (Polyarthra longiremis, Brachionus angularis, Cyclops vicinus, Codonella cratera) and the positive and negative relationships between them were balanced. The MT biocenotic network was composed of three sub-networks connected by nodes with high communication attributes (Polyarthra vulgaris, Bosmina longirostris, C. vicinus), and antagonistic interactions (predation and competition) were less important. The LT network was most heterogeneous, and Daphnia cuculllata exerted the strongest influence on the network's structure by forming numerous positive (coexistence with predators) and negative (interference competition with microphagous rotifers) interspecific relationships. The study provides new information about the ecology of aquatic ecosystems, that are disturbed by changes in water turbidity.

A quantitative AOP of mitochondrial toxicity based on data from three cell lines

Adverse Outcome Pathways (AOPs) are increasingly used to support the integration of in vitro data in hazard assessment for chemicals. Quantitative AOPs (qAOPs) use mathematical models to describe the relationship between key events (KEs). In this paper, data obtained in three cell lines, LHUMES, HepG2 and RPTEC/TERT1, using similar experimental protocols, was used to calibrate a qAOP of mitochondrial toxicity for two chemicals, rotenone and deguelin. The objectives were to determine whether the same qAOP could be used for the three cell types, and to test chemical-independence by cross-validation with a dataset obtained on eight other chemicals in LHUMES cells. Repeating the calibration approach for both chemicals in three cell lines highlighted various practical difficulties. Even when the same readouts of KEs are measured, the mathematical functions used to describe the key event relationships may not be the same. Cross-validation in LHUMES cells was attempted by estimating chemical-specific potency at the molecular initiating events and using the rest of the calibrated qAOP to predict downstream KEs: toxicity of azoxystrobin, carboxine, mepronil and thifluzamide was underestimated. Selection of most relevant readouts and accurate characterization of the molecular initiating event for cross-validation are critical when designing in vitro experiments targeted at calibrating qAOPs.

Network Analysis Reveals Different Cellulose Degradation Strategies Across Trichoderma harzianum Strains Associated With XYR1 and CRE1

Trichoderma harzianum, whose gene expression is tightly controlled by the transcription factors (TFs) XYR1 and CRE1, is a potential candidate for hydrolytic enzyme production. Here, we performed a network analysis of T. harzianum IOC-3844 and T. harzianum CBMAI-0179 to explore how the regulation of these TFs varies between these strains. In addition, we explored the evolutionary relationships of XYR1 and CRE1 protein sequences among Trichoderma spp. The results of the T. harzianum strains were compared with those of Trichoderma atroviride CBMAI-0020, a mycoparasitic species. Although transcripts encoding carbohydrate-active enzymes (CAZymes), TFs, transporters, and proteins with unknown functions were coexpressed with cre1 or xyr1, other proteins indirectly related to cellulose degradation were identified. The enriched GO terms describing the transcripts of these groups differed across all strains, and several metabolic pathways with high similarity between both regulators but strain-specific differences were identified. In addition, the CRE1 and XYR1 subnetworks presented different topology profiles in each strain, likely indicating differences in the influences of these regulators according to the fungi. The hubs of the cre1 and xyr1 groups included transcripts not yet characterized or described as being related to cellulose degradation. The first-neighbor analyses confirmed the results of the profile of the coexpressed transcripts in cre1 and xyr1. The analyses of the shortest paths revealed that CAZymes upregulated under cellulose degradation conditions are most closely related to both regulators, and new targets between such signaling pathways were discovered. Although the evaluated T. harzianum strains are phylogenetically close and their amino acid sequences related to XYR1 and CRE1 are very similar, the set of transcripts related to xyr1 and cre1 differed, suggesting that each T. harzianum strain used a specific regulation strategy for cellulose degradation. More interestingly, our findings may suggest that XYR1 and CRE1 indirectly regulate genes encoding proteins related to cellulose degradation in the evaluated T. harzianum strains. An improved understanding of the basic biology of fungi during the cellulose degradation process can contribute to the use of their enzymes in several biotechnological applications and pave the way for further studies on the differences across strains of the same species.

Multidimensional analysis of the host response reveals prognostic and pathogen-driven immune subtypes among adults with sepsis in Uganda

Background

The global burden of sepsis is concentrated in sub-Saharan Africa, where severe infections disproportionately affect young, HIV-infected adults and high-burden pathogens are unique. In this context, poor understanding of sepsis immunopathology represents a crucial barrier to development of locally-effective treatment strategies. We sought to determine inter-individual immunologic heterogeneity among adults hospitalized with sepsis in a sub-Saharan African setting, and characterize associations between immune subtypes, infecting pathogens, and clinical outcomes.

Methods

Among a prospective observational cohort of 288 adults hospitalized with suspected sepsis in Uganda, we applied machine learning methods to 14 soluble host immune mediators, reflective of key domains of sepsis immunopathology (innate and adaptive immune activation, endothelial dysfunction, fibrinolysis), to identify immune subtypes in randomly-split discovery (N = 201) and internal validation (N = 87) sub-cohorts. In parallel, we applied similar methods to whole-blood RNA-sequencing data from a consecutive subset of patients (N = 128) to identify transcriptional subtypes, which we characterized using biological pathway and immune cell-type deconvolution analyses.

Results

Unsupervised clustering consistently identified two immune subtypes defined by differential activation of pro-inflammatory innate and adaptive immune pathways, with transcriptional evidence of concomitant CD56(-)/CD16( +) NK-cell expansion, T-cell exhaustion, and oxidative-stress and hypoxia-induced metabolic and cell-cycle reprogramming in the hyperinflammatory subtype. Immune subtypes defined by greater pro-inflammatory immune activation, T-cell exhaustion, and metabolic reprogramming were consistently associated with a high-prevalence of severe and often disseminated HIV-associated tuberculosis, as well as more extensive organ dysfunction, worse functional outcomes, and higher 30-day mortality.

Conclusions

Our results highlight unique host- and pathogen-driven features of sepsis immunopathology in sub-Saharan Africa, including the importance of severe HIV-associated tuberculosis, and reinforce the need to develop more biologically-informed treatment strategies in the region, particularly those incorporating immunomodulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13054-022-03907-3.

Role of mitochondrial genetic interactions in determining adaptation to high altitude human population

Physiological and haplogroup studies performed to understand high-altitude adaptation in humans are limited to individual genes and polymorphic sites. Due to stochastic evolutionary forces, the frequency of a polymorphism is affected by changes in the frequency of a near-by polymorphism on the same DNA sample making them connected in terms of evolution. Here, first, we provide a method to model these mitochondrial polymorphisms as "co-mutation networks" for three high-altitude populations, Tibetan, Ethiopian and Andean. Then, by transforming these co-mutation networks into weighted and undirected gene-gene interaction (GGI) networks, we were able to identify functionally enriched genetic interactions of CYB and CO3 genes in Tibetan and Andean populations, while NADH dehydrogenase genes in the Ethiopian population playing a significant role in high altitude adaptation. These co-mutation based genetic networks provide insights into the role of different set of genes in high-altitude adaptation in human sub-populations.

A practical guide to interpreting and generating bottom-up proteomics data visualizations

Mass-spectrometry based bottom-up proteomics is the main method to analyze proteomes comprehensively and the rapid evolution of instrumentation and data analysis has made the technology widely available. Data visualization is an integral part of the analysis process and it is crucial for the communication of results. This is a major challenge due to the immense complexity of MS data. In this review, we provide an overview of commonly used visualizations, starting with raw data of traditional and novel MS technologies, then basic peptide and protein level analyses, and finally visualization of highly complex datasets and networks. We specifically provide guidance on how to critically interpret and discuss the multitude of different proteomics data visualizations. Furthermore, we highlight Python-based libraries and other open science tools that can be applied for independent and transparent generation of customized visualizations. To further encourage programmatic data visualization, we provide the Python code used to generate all data Figures in this review on GitHub (https://github.com/MannLabs/ProteomicsVisualization). This article is protected by copyright. All rights reserved.

MathFeature: feature extraction package for DNA, RNA and protein sequences based on mathematical descriptors

One of the main challenges in applying machine learning algorithms to biological sequence data is how to numerically represent a sequence in a numeric input vector. Feature extraction techniques capable of extracting numerical information from biological sequences have been reported in the literature. However, many of these techniques are not available in existing packages, such as mathematical descriptors. This paper presents a new package, MathFeature, which implements mathematical descriptors able to extract relevant numerical information from biological sequences, i.e. DNA, RNA and proteins (prediction of structural features along the primary sequence of amino acids). MathFeature makes available 20 numerical feature extraction descriptors based on approaches found in the literature, e.g. multiple numeric mappings, genomic signal processing, chaos game theory, entropy and complex networks. MathFeature also allows the extraction of alternative features, complementing the existing packages. To ensure that our descriptors are robust and to assess their relevance, experimental results are presented in nine case studies. According to these results, the features extracted by MathFeature showed high performance (0.6350-0.9897, accuracy), both applying only mathematical descriptors, but also hybridization with well-known descriptors in the literature. Finally, through MathFeature, we overcame several studies in eight benchmark datasets, exemplifying the robustness and viability of the proposed package. MathFeature has advanced in the area by bringing descriptors not available in other packages, as well as allowing non-experts to use feature extraction techniques.

The Genomic Physics of COVID-19 Pathogenesis and Spread

Coronavirus disease (COVID-19) spreads mainly through close contact of infected persons, but the molecular mechanisms underlying its pathogenesis and transmission remain unknown. Here, we propose a statistical physics model to coalesce all molecular entities into a cohesive network in which the roadmap of how each entity mediates the disease can be characterized. We argue that the process of how a transmitter transforms the virus into a recipient constitutes a triad unit that propagates COVID-19 along reticulate paths. Intrinsically, person-to-person transmissibility may be mediated by how genes interact transversely across transmitter, recipient, and viral genomes. We integrate quantitative genetic theory into hypergraph theory to code the main effects of the three genomes as nodes, pairwise cross-genome epistasis as edges, and high-order cross-genome epistasis as hyperedges in a series of mobile hypergraphs. Charting a genome-wide atlas of horizontally epistatic hypergraphs can facilitate the systematic characterization of the community genetic mechanisms underlying COVID-19 spread. This atlas can typically help design effective containment and mitigation strategies and screen and triage those more susceptible persons and those asymptomatic carriers who are incubation virus transmitters.

Combinatorial CRISPR/Cas9 Screening Reveals Epistatic Networks of Interacting Tumor Suppressor Genes and Therapeutic Targets in Human Breast Cancer

The majority of cancers are driven by multiple genetic alterations, but how these changes collaborate during tumorigenesis remains largely unknown. To gain mechanistic insights into tumor-promoting genetic interactions among tumor suppressor genes (TSG), we conducted combinatorial CRISPR screening coupled with single-cell transcriptomic profiling in human mammary epithelial cells. As expected, different driver gene alterations in mammary epithelial cells influenced the repertoire of tumor suppressor alterations capable of inducing tumor formation. More surprisingly, TSG interaction networks were comprised of numerous cliques-sets of three or four genes such that each TSG within the clique showed oncogenic cooperation with all other genes in the clique. Genetic interaction profiling indicated that the predominant cooperating TSGs shared overlapping functions rather than distinct or complementary functions. Single-cell transcriptomic profiling of CRISPR double knockouts revealed that cooperating TSGs that synergized in promoting tumorigenesis and growth factor independence showed transcriptional epistasis, whereas noncooperating TSGs did not. These epistatic transcriptional changes, both buffering and synergistic, affected expression of oncogenic mediators and therapeutic targets, including CDK4, SRPK1, and DNMT1. Importantly, the epistatic expression alterations caused by dual inactivation of TSGs in this system, such as PTEN and TP53, were also observed in patient tumors, establishing the relevance of these findings to human breast cancer. An estimated 50% of differentially expressed genes in breast cancer are controlled by epistatic interactions. Overall, our study indicates that transcriptional epistasis is a central aspect of multigenic breast cancer progression and outlines methodologies to uncover driver gene epistatic networks in other human cancers. SIGNIFICANCE: This study provides a roadmap for moving beyond discovery and development of therapeutic strategies based on single driver gene analysis to discovery based on interactions between multiple driver genes.See related commentary by Fong et al., p. 6078.

GPCRs Are Optimal Regulators of Complex Biological Systems and Orchestrate the Interface between Health and Disease

GPCRs arguably represent the most effective current therapeutic targets for a plethora of diseases. GPCRs also possess a pivotal role in the regulation of the physiological balance between healthy and pathological conditions; thus, their importance in systems biology cannot be underestimated. The molecular diversity of GPCR signaling systems is likely to be closely associated with disease-associated changes in organismal tissue complexity and compartmentalization, thus enabling a nuanced GPCR-based capacity to interdict multiple disease pathomechanisms at a systemic level. GPCRs have been long considered as controllers of communication between tissues and cells. This communication involves the ligand-mediated control of cell surface receptors that then direct their stimuli to impact cell physiology. Given the tremendous success of GPCRs as therapeutic targets, considerable focus has been placed on the ability of these therapeutics to modulate diseases by acting at cell surface receptors. In the past decade, however, attention has focused upon how stable multiprotein GPCR superstructures, termed receptorsomes, both at the cell surface membrane and in the intracellular domain dictate and condition long-term GPCR activities associated with the regulation of protein expression patterns, cellular stress responses and DNA integrity management. The ability of these receptorsomes (often in the absence of typical cell surface ligands) to control complex cellular activities implicates them as key controllers of the functional balance between health and disease. A greater understanding of this function of GPCRs is likely to significantly augment our ability to further employ these proteins in a multitude of diseases.

Multi-omics approach to COVID-19: a domain-based literature review

Background

Omics data, driven by rapid advances in laboratory techniques, have been generated very quickly during the COVID-19 pandemic. Our aim is to use omics data to highlight the involvement of specific pathways, as well as that of cell types and organs, in the pathophysiology of COVID-19, and to highlight their links with clinical phenotypes of SARS-CoV-2 infection.

Methods

The analysis was based on the domain model, where for domain it is intended a conceptual repository, useful to summarize multiple biological pathways involved at different levels. The relevant domains considered in the analysis were: virus, pathways and phenotypes. An interdisciplinary expert working group was defined for each domain, to carry out an independent literature scoping review.

Results

The analysis revealed that dysregulated pathways of innate immune responses, (i.e., complement activation, inflammatory responses, neutrophil activation and degranulation, platelet degranulation) can affect COVID-19 progression and outcomes. These results are consistent with several clinical studies.

Conclusions

Multi-omics approach may help to further investigate unknown aspects of the disease. However, the disease mechanisms are too complex to be explained by a single molecular signature and it is necessary to consider an integrated approach to identify hallmarks of severity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03168-8.

Identification of Key Genes in 'Luang Pratahn', Thai Salt-Tolerant Rice, Based on Time-Course Data and Weighted Co-expression Networks

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

Long non-coding RNAs and their potential impact on diagnosis, prognosis, and therapy in prostate cancer: racial, ethnic, and geographical considerations

INTRODUCTION

Advances in high-throughput sequencing have greatly advanced our understanding of long non-coding RNAs (lncRNAs) in a relatively short period of time. This has expanded our knowledge of cancer, particularly how lncRNAs drive many important cancer phenotypes via their regulation of gene expression.

AREAS COVERED

Men of African descent are disproportionately affected by PC in terms of incidence, morbidity, and mortality. LncRNAs could serve as biomarkers to differentiate low-risk from high-risk diseases. Additionally, they may represent therapeutic targets for advanced and castrate-resistant cancer. We review current research surrounding lncRNAs and their association with PC. We discuss how lncRNAs can provide new insights and diagnostic biomarkers for African American men. Finally, we review advances in computational approaches that predict the regulatory effects of lncRNAs in cancer.

EXPERT OPINION

PC diagnostic biomarkers that offer high specificity and sensitivity are urgently needed. PC specific lncRNAs are compelling as diagnostic biomarkers owing to their high tissue and tumor specificity and presence in bodily fluids. Recent studies indicate that PCA3 clinical utility might be restricted to men of European descent. Further work is required to develop lncRNA biomarkers tailored for men of African descent.

Challenges and opportunities in network-based solutions for biological questions

Network biology is useful for modeling complex biological phenomena; it has attracted attention with the advent of novel graph-based machine learning methods. However, biological applications of network methods often suffer from inadequate follow-up. In this perspective, we discuss obstacles for contemporary network approaches—particularly focusing on challenges representing biological concepts, applying machine learning methods, and interpreting and validating computational findings about biology—in an effort to catalyze actionable biological discovery.

Mechanical couplings of protein backbone and side chains exhibit scale-free network properties and specific hotspots for function

•

Statistical learning from protein dynamics unravels rigidities in interaction network.

•

Backbone and side-chain mechanical couplings exhibit scale-free network properties.

•

Graphical depiction of network rigidities captures sequence co-evolution patterns.

•

Functional sites at secondary structure peripheries are mechanical hotspots.

•

Our rigidity scores are compelling metrics for residue biological significance.

A backbone-side-chain elastic network model (bsENM) is devised in this contribution to decipher the network of molecular interactions during protein dynamics. The chemical details in 5 μs all-atom molecular dynamics (MD) simulation are mapped onto the bsENM spring constants by self-consistent iterations. The elastic parameters obtained by this structure-mechanics statistical learning are then used to construct inter-residue rigidity graphs for the chemical components in protein amino acids. A key discovery is that the mechanical coupling strengths of both backbone and side chains exhibit heavy-tailed distributions and scale-free network properties. In both rat trypsin and PDZ3 proteins, the statistically prominent modes of rigidity graphs uncover the sequence-specific coupling patterns and mechanical hotspots. Based on the contributions to graphical modes, our residue rigidity scores in backbone and side chains are found to be very useful metrics for the biological significance. Most functional sites have high residue rigidity scores in side chains while the biologically important glycines are generally next to mechanical hotspots. Furthermore, prominent modes in the rigidity graphs involving side chains oftentimes coincide with the co-evolution patterns due to evolutionary restraints. The bsENM specifically devised to resolve the protein chemical character thus provides useful means for extracting functional information from all-atom MD.

Delineating infection strategies of Leishmania donovani secretory proteins in Human through host-pathogen protein Interactome prediction

Interactions of Leishmania donovani secretory virulence factors with the host proteins and their interplay during the infection process in humans is poorly studied in Visceral Leishmaniasis. Lack of a holistic study of pathway level de-regulations caused due to these virulence factors leads to a poor understanding of the parasite strategies to subvert the host immune responses, secure its survival inside the host and further the spread of infection to the visceral organs. In this study, we propose a computational workflow to predict host-pathogen protein interactome of L.donovani secretory virulence factors with human proteins combining sequence-based Interolog mapping and structure-based Domain Interaction mapping techniques. We further employ graph theoretical approaches and shortest path methods to analyze the interactome. Our study deciphers the infection paths involving some unique and understudied disease-associated signaling pathways influencing the cellular phenotypic responses in the host. Our statistical analysis based in silico knockout study unveils for the first time UBC, 1433Z and HS90A mediator proteins as potential immunomodulatory candidates through which the virulence factors employ the infection paths. These identified pathways and novel mediator proteins can be effectively used as possible targets to control and modulate the infection process further aiding in the treatment of Visceral Leishmaniasis.

Network Biology Approaches to Achieve Precision Medicine in Inflammatory Bowel Disease

Inflammatory bowel disease (IBD) is a chronic immune-mediated condition arising due to complex interactions between multiple genetic and environmental factors. Despite recent advances, the pathogenesis of the condition is not fully understood and patients still experience suboptimal clinical outcomes. Over the past few years, investigators are increasingly capturing multi-omics data from patient cohorts to better characterise the disease. However, reaching clinically translatable endpoints from these complex multi-omics datasets is an arduous task. Network biology, a branch of systems biology that utilises mathematical graph theory to represent, integrate and analyse biological data through networks, will be key to addressing this challenge. In this narrative review, we provide an overview of various types of network biology approaches that have been utilised in IBD including protein-protein interaction networks, metabolic networks, gene regulatory networks and gene co-expression networks. We also include examples of multi-layered networks that have combined various network types to gain deeper insights into IBD pathogenesis. Finally, we discuss the need to incorporate other data sources including metabolomic, histopathological, and high-quality clinical meta-data. Together with more robust network data integration and analysis frameworks, such efforts have the potential to realise the key goal of precision medicine in IBD.

Using BioPAX-Parser (BiP) to enrich lists of genes or proteins with pathway data

BACKGROUND

Pathway enrichment analysis (PEA) is a well-established methodology for interpreting a list of genes and proteins of interest related to a condition under investigation. This paper aims to extend our previous work in which we introduced a preliminary comparative analysis of pathway enrichment analysis tools. We extended the earlier work by providing more case studies, comparing BiP enrichment performance with other well-known PEA software tools.

METHODS

PEA uses pathway information to discover connections between a list of genes and proteins as well as biological mechanisms, helping researchers to overcome the problem of explaining biological entity lists of interest disconnected from the biological context.

RESULTS

We compared the results of BiP with some existing pathway enrichment analysis tools comprising Centrality-based Pathway Enrichment, pathDIP, and Signaling Pathway Impact Analysis, considering three cancer types (colorectal, endometrial, and thyroid), for a total of six datasets (that is, two datasets per cancer type) obtained from the The Cancer Genome Atlas and Gene Expression Omnibus databases. We measured the similarities between the overlap of the enrichment results obtained using each couple of cancer datasets related to the same cancer.

CONCLUSION

As a result, BiP identified some well-known pathways related to the investigated cancer type, validated by the available literature. We also used the Jaccard and meet-min indices to evaluate the stability and the similarity between the enrichment results obtained from each couple of cancer datasets. The obtained results show that BiP provides more stable enrichment results than other tools.

Cell-specific gene association network construction from single-cell RNA sequence

The recent development of a high throughput single-cell RNA sequence devises the opportunity to study entire transcriptomes in the smallest detail. It also leads to the characterization of molecules and subtypes of a cell. Cancer epigenetics induced not only from individual molecules but also from the dysfunction of the system and the coupling effect of genes. While rapid advances are being made in the development of tools for single-cell RNA-seq data analysis, few slants are noticed in the potential advantages of single-cell network construction.Here, we used network perturbation theory with significant analysis to develop a cell-specific network that provides an insight into gene-gene association based on molecular expressions in a single-cell resolution. Besides, using this method, we can characterize each cell by inspecting how genes are connected and can identify the hub genes using network degree theory. Pathway & Gene enrichment analysis of the identified cell-specific high network degree genes supported the effectiveness of this method. This method could be beneficial for personalized drug design and even therapeutics.

NetCom: A Network-Based Tool for Predicting Metabolic Activities of Microbial Communities Based on Interpretation of Metagenomics Data

The study of microbial activity can be viewed as a triangle with three sides: environment (dominant resources in a specific habitat), community (species dictating a repertoire of metabolic conversions) and function (production and/or utilization of resources and compounds). Advances in metagenomics enable a high-resolution description of complex microbial communities in their natural environments and support a systematic study of environment-community-function associations. NetCom is a web-tool for predicting metabolic activities of microbial communities based on network-based interpretation of assembled and annotated metagenomics data. The algorithm takes as an input, lists of differentially abundant enzymatic reactions and generates the following outputs: (i) pathway associations of differently abundant enzymes; (ii) prediction of environmental resources that are unique to each treatment, and their pathway associations; (iii) prediction of compounds that are produced by the microbial community, and pathway association of compounds that are treatment-specific; (iv) network visualization of enzymes, environmental resources and produced compounds, that are treatment specific (2 and 3D). The tool is demonstrated on metagenomic data from rhizosphere and bulk soil samples. By predicting root-specific activities, we illustrate the relevance of our framework for forecasting the impact of soil amendments on the corresponding microbial communities. NetCom is available online.

Cost-Effective Mapping of Genetic Interactions in Mammalian Cells

Comprehensive maps of genetic interactions in mammalian cells are daunting to construct because of the large number of potential interactions, ~ 2 × 10⁸ for protein coding genes. We previously used co-inheritance of distant genes from published radiation hybrid (RH) datasets to identify genetic interactions. However, it was necessary to combine six legacy datasets from four species to obtain adequate statistical power. Mapping resolution was also limited by the low density PCR genotyping. Here, we employ shallow sequencing of nascent human RH clones as an economical approach to constructing interaction maps. In this initial study, 15 clones were analyzed, enabling construction of a network with 225 genes and 2,359 interactions (FDR < 0.05). Despite its small size, the network showed significant overlap with the previous RH network and with a protein-protein interaction network. Consumables were ≲$50 per clone, showing that affordable, high quality genetic interaction maps are feasible in mammalian cells.

Biomolecule and Bioentity Interaction Databases in Systems Biology: A Comprehensive Review

Technological advances in high-throughput techniques have resulted in tremendous growth of complex biological datasets providing evidence regarding various biomolecular interactions. To cope with this data flood, computational approaches, web services, and databases have been implemented to deal with issues such as data integration, visualization, exploration, organization, scalability, and complexity. Nevertheless, as the number of such sets increases, it is becoming more and more difficult for an end user to know what the scope and focus of each repository is and how redundant the information between them is. Several repositories have a more general scope, while others focus on specialized aspects, such as specific organisms or biological systems. Unfortunately, many of these databases are self-contained or poorly documented and maintained. For a clearer view, in this article we provide a comprehensive categorization, comparison and evaluation of such repositories for different bioentity interaction types. We discuss most of the publicly available services based on their content, sources of information, data representation methods, user-friendliness, scope and interconnectivity, and we comment on their strengths and weaknesses. We aim for this review to reach a broad readership varying from biomedical beginners to experts and serve as a reference article in the field of Network Biology.

The molecular profile of mucosal melanoma

Herein, we wanted to explore the molecular landscape of mucosal melanoma from different sites and identify potential molecular targets for future therapy. Mucosal melanomas (N = 40) from different sites (conjunctiva, sinonasal cavity, rectum, and vagina) were investigated. Targeted next-generation sequencing along with Nanostring gene expression profiling was performed. Genetically, conjunctival melanoma was characterized by BRAF-V600E (30%) and NF1 mutations (17%). Mucosal melanomas at nonsun-exposed sites harbored alterations in NRAS, KIT, NF1, along with atypical BRAF mutations. When comparing the gene expression profile of conjunctival melanoma and nonsun-exposed mucosal melanoma, 41 genes were found to be significantly deregulated. Programmed death-ligand 1 (PD-L1) presented a significant sixfold upregulation in conjunctival melanoma compared to the other mucosal melanomas. While melanomas of the sinonasal cavity, vagina, and rectum are molecularly similar, conjunctival melanoma is characterized by a higher frequency of BRAF-V600E mutations and differential expression of several genes involved in the immune response.