Correlation between gene expression profiles and protein-protein interactions within and across genomes

AVERON notebook to discover actionable cancer vulnerabilities enabled by neomorph protein-protein interactions

Genomic alterations, such as missense mutations, often lead to the activation of oncogenic pathways and cell transformation by rewiring protein-protein interaction (PPI) networks. Understanding how mutant-directed neomorph PPIs (neoPPIs) drive cancer is vital to developing new personalized clinical strategies. However, the experimental interrogation of neoPPI functions in patients with cancer is highly challenging. To address this challenge, we developed a computational platform, termed AVERON for discovering actionable vulnerabilities enabled by rewired oncogenic networks. AVERON enables rapid systematic profiling of the clinical significance of neomorph PPIs across different cancer types, informing molecular mechanisms of neoPPI-driven tumorigenesis, and revealing therapeutically actionable neoPPI-regulated genes. We demonstrated the application of the AVERON platform by evaluating the biological functions and clinical significance of 130 neomorph interactions, experimentally determined for oncogenic BRAFV600E. The AVERON application to broad sets of mutant-directed PPIs may inform new testable biological models and clinical strategies in cancer.

ACDMBI: A deep learning model based on community division and multi-source biological information fusion predicts essential proteins

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.

Predicting evolutionary targets and parameters of gene deletion from expression data

Motivation

Gene deletion is traditionally thought of as a nonadaptive process that removes functional redundancy from genomes, such that it generally receives less attention than duplication in evolutionary turnover studies. Yet, mounting evidence suggests that deletion may promote adaptation via the "less-is-more" evolutionary hypothesis, as it often targets genes harboring unique sequences, expression profiles, and molecular functions. Hence, predicting the relative prevalence of redundant and unique functions among genes targeted by deletion, as well as the parameters underlying their evolution, can shed light on the role of gene deletion in adaptation.

Results

Here, we present CLOUDe, a suite of machine learning methods for predicting evolutionary targets of gene deletion events from expression data. Specifically, CLOUDe models expression evolution as an Ornstein-Uhlenbeck process, and uses multi-layer neural network, extreme gradient boosting, random forest, and support vector machine architectures to predict whether deleted genes are "redundant" or "unique", as well as several parameters underlying their evolution. We show that CLOUDe boasts high power and accuracy in differentiating between classes, and high accuracy and precision in estimating evolutionary parameters, with optimal performance achieved by its neural network architecture. Application of CLOUDe to empirical data from Drosophila suggests that deletion primarily targets genes with unique functions, with further analysis showing these functions to be enriched for protein deubiquitination. Thus, CLOUDe represents a key advance in learning about the role of gene deletion in functional evolution and adaptation.

Availability and implementation

CLOUDe is freely available on GitHub (https://github.com/anddssan/CLOUDe).

Network Analytics Enabled by Generating a Pool of Network Variants from Noisy Data

Mapping network nodes and edges to communities and network functions is crucial to gaining a higher level of understanding of the network structure and functions. Such mappings are particularly challenging to design for covert social networks, which intentionally hide their structure and functions to protect important members from attacks or arrests. Here, we focus on correctly inferring the structures and functions of such networks, but our methodology can be broadly applied. Without the ground truth, knowledge about the allocation of nodes to communities and network functions, no single network based on the noisy data can represent all plausible communities and functions of the true underlying network. To address this limitation, we apply a generative model that randomly distorts the original network based on the noisy data, generating a pool of statistically equivalent networks. Each unique generated network is recorded, while each duplicate of the already recorded network just increases the repetition count of that network. We treat each such network as a variant of the ground truth with the probability of arising in the real world approximated by the ratio of the count of this network's duplicates plus one to the total number of all generated networks. Communities of variants with frequently occurring duplicates contain persistent patterns shared by their structures. Using Shannon entropy, we can find a variant that minimizes the uncertainty for operations planned on the network. Repeatedly generating new pools of networks from the best network of the previous step for several steps lowers the entropy of the best new variant. If the entropy is too high, the network operators can identify nodes, the monitoring of which can achieve the most significant reduction in entropy. Finally, we also present a heuristic for constructing a new variant, which is not randomly generated but has the lowest expected cost of operating on the distorted mappings of network nodes to communities and functions caused by noisy data.

Narrow funnel-like interaction energy distribution is an indicator of specific protein interaction partner

Protein interaction networks underlie countless biological mechanisms. However, most protein interaction predictions are based on biological evidence that are biased to well-known protein interaction or physical evidence that exhibits low accuracy for weak interactions and requires high computational power. In this study, a novel method has been suggested to predict protein interaction partners by investigating narrow funnel-like interaction energy distribution. In this study, it was demonstrated that various protein interactions including kinases and E3 ubiquitin ligases have narrow funnel-like interaction energy distribution. To analyze protein interaction distribution, modified scores of iRMS and TM-score are introduced. Then, using these scores, algorithm and deep learning model for prediction of protein interaction partner and substrate of kinase and E3 ubiquitin ligase were developed. The prediction accuracy was similar to or even better than that of yeast two-hybrid screening. Ultimately, this knowledge-free protein interaction prediction method will broaden our understanding of protein interaction networks.

The protein kinases of Dictyostelia and their incorporation into a signalome

Protein kinases are major regulators of cellular processes, but the roles of most kinases remain unresolved. Dictyostelid social amoebas have been useful in identifying functions for 30% of its kinases in cell migration, cytokinesis, vesicle trafficking, gene regulation and other processes but their upstream regulators and downstream effectors are mostly unknown. Comparative genomics can assist to distinguish between genes involved in deeply conserved core processes and those involved in species-specific innovations, while co-expression of genes as evident from comparative transcriptomics can provide cues to the protein complement of regulatory networks. Genomes and developmental and cell-type specific transcriptomes are available for species that span the 0.5 billion years of evolution of Dictyostelia from their unicellular ancestors. In this work we analysed conservation and change in the abundance, functional domain architecture and developmental regulation of protein kinases across the 4 major taxon groups of Dictyostelia. All data are summarized in annotated phylogenetic trees of the kinase subtypes and accompanied by functional information of all kinases that were experimentally studied. We detected 393 different protein kinase domains across the five studied genomes, of which 212 were fully conserved. Conservation was highest (71%) in the previously defined AGC, CAMK, CK1, CMCG, STE and TKL groups and lowest (26%) in the "other" group of typical protein kinases. This was mostly due to species-specific single gene amplification of "other" kinases. Apart from the AFK and α-kinases, the atypical protein kinases, such as the PIKK and histidine kinases were also almost fully conserved. The phylogeny-wide developmental and cell-type specific expression profiles of the protein kinase genes were combined with profiles from the same transcriptomic experiments for the families of G-protein coupled receptors, small GTPases and their GEFs and GAPs, the transcription factors and for all genes that upon lesion generate a developmental defect. This dataset was subjected to hierarchical clustering to identify clusters of co-expressed genes that potentially act together in a signalling network. The work provides a valuable resource that allows researchers to identify protein kinases and other regulatory proteins that are likely to act as intermediates in a network of interest.

Predicting gene expression divergence between single-copy orthologs in two species

Predicting gene expression divergence is integral to understanding the emergence of new biological functions and associated traits. Whereas several sophisticated methods have been developed for this task, their applications are either limited to duplicate genes or require expression data from more than two species. Thus, here we present PiXi, the first machine learning framework for predicting gene expression divergence between single-copy orthologs in two species. PiXi models gene expression evolution as an Ornstein-Uhlenbeck process, and overlays this model with multi-layer neural network, random forest, and support vector machine architectures for making predictions. It outputs the predicted class "conserved" or "diverged" for each pair of orthologs, as well as their predicted expression optima in the two species. We show that PiXi has high power and accuracy in predicting gene expression divergence between single-copy orthologs, as well as high accuracy and precision in estimating their expression optima in the two species, across a wide range of evolutionary scenarios, with the globally best performance achieved by a multi-layer neural network. Moreover, application of our best performing PiXi predictor to empirical gene expression data from single-copy orthologs residing at different loci in two species of Drosophila reveals that approximately 23% underwent expression divergence after positional relocation. Further analysis shows that several of these "diverged" genes are involved in the electron transport chain of the mitochondrial membrane, suggesting that new chromatin environments may impact energy production in Drosophila. Thus, by providing a toolkit for predicting gene expression divergence between single-copy orthologs in two species, PiXi can shed light on the origins of novel phenotypes across diverse biological processes and study systems.

Protocol to analyze dysregulation of the eIF4F complex in human cancers using R software and large public datasets

Understanding dysregulation of the eukaryotic initiation factor 4F (eIF4F) complex across tumor types is critical to cancer treatment development. We present a protocol and accompanying R package "eIF4F.analysis". We describe analysis of copy number status, gene abundance and stoichiometry, survival probability, expression covariation, correlating genes, mRNA/protein correlation, and protein co-expression. Using publicly available large multi-omics data, eIF4F.analysis permits computationally derived and statistically powerful inferences regarding initiation factor regulation in human cancers and clinical relevance of protein interactions within the eIF4F complex. For complete details on the use and execution of this protocol, please refer to Wu and Wagner (2021).¹.

Gut Microbiota Alterations in Trace Amine-Associated Receptor 9 (TAAR9) Knockout Rats

Trace amine-associated receptors (TAAR1-TAAR9) are a family of G-protein-coupled monoaminergic receptors which might have great pharmacological potential. It has now been well established that TAAR1 plays an important role in the central nervous system. Interestingly, deletion of TAAR9 in rats leads to alterations in the periphery. Previously, we found that knockout of TAAR9 in rats (TAAR9-KO rats) decreased low-density lipoprotein cholesterol levels in the blood. TAAR9 was also identified in intestinal tissues, and it is known that it responds to polyamines. To elucidate the role of TAAR9 in the intestinal epithelium, we analyzed TAAR9-co-expressed gene clusters in public data for cecum samples. As identified by gene ontology enrichment analysis, in the intestine, TAAR9 is co-expressed with genes involved in intestinal mucosa homeostasis and function, including cell organization, differentiation, and death. Additionally, TAAR9 was co-expressed with genes implicated in dopamine signaling, which may suggest a role for this receptor in the regulation of peripheral dopaminergic transmission. To further investigate how TAAR9 might be involved in colonic mucosal homeostasis, we analyzed the fecal microbiome composition in TAAR9-KO rats and their wild-type littermates. We identified a significant difference in the number of observed taxa between the microbiome of TAAR9-KO and wild-type rats. In TAAR9-KO rats, the gut microbial community became more variable compared with the wild-type rats. Furthermore, it was found that the family Saccharimonadaceae, which is one of the top 10 most abundant families in TAAR9-KO rat feces, is almost completely absent in wild-type animal fecal samples. Taken together, these data indicate a role of TAAR9 in intestinal function.

Interacting partners of Brassica juncea Regulator of G-protein Signaling protein suggest its role in cell wall metabolism and cellular signaling

Heterotrimeric G-proteins interact with various upstream and downstream effectors to regulate various aspects of plant growth and development. G-protein effectors have been recently reported in Arabidopsis thaliana; however, less information is available from polyploid crop species having complex networks of G-protein components. Regulator of G-protein signaling (RGS) is a well-characterized GTPase accelerating protein, which plays an important role in the regulation of the G-protein cycle in plants. In the present study, four homologs encoding RGS proteins were isolated from the allotetraploid Brassica juncea, a globally important oilseed, vegetable, and condiment crop. The B. juncea RGS proteins were grouped into distinct BjuRGS1 and BjuRGS2 orthologous clades, and the expression of BjuRGS1 homologs was predominantly higher than BjuRGS2 homologs across the tested tissue types of B. juncea. Utilizing B. juncea Y2H library screening, a total of 30 nonredundant interacting proteins with the RGS-domain of the highly expressed BjuA.RGS1 was identified. Gene ontology analysis indicated that these effectors exerted various molecular, cellular, and physiological functions. Many of them were known to regulate cell wall metabolism (BjuEXP6, Bju-α-MAN, BjuPGU4, BjuRMS3) and phosphorylation-mediated cell signaling (BjuMEK4, BjuDGK3, and BjuKinase). Furthermore, transcript analysis indicated that the identified interacting proteins have a coexpression pattern with the BjuRGS homologs. These findings increase our knowledge about the novel targets of G-protein components from a globally cultivated Brassica crop and provide an important resource for developing a plant G-protein interactome network.

Computationally Reconstructed Interactome of Bradyrhizobium diazoefficiens USDA110 Reveals Novel Functional Modules and Protein Hubs for Symbiotic Nitrogen Fixation

Symbiotic nitrogen fixation is an important part of the nitrogen biogeochemical cycles and the main nitrogen source of the biosphere. As a classical model system for symbiotic nitrogen fixation, rhizobium-legume systems have been studied elaborately for decades. Details about the molecular mechanisms of the communication and coordination between rhizobia and host plants is becoming clearer. For more systematic insights, there is an increasing demand for new studies integrating multiomics information. Here, we present a comprehensive computational framework integrating the reconstructed protein interactome of B. diazoefficiens USDA110 with its transcriptome and proteome data to study the complex protein-protein interaction (PPI) network involved in the symbiosis system. We reconstructed the interactome of B. diazoefficiens USDA110 by computational approaches. Based on the comparison of interactomes between B. diazoefficiens USDA110 and other rhizobia, we inferred that the slow growth of B. diazoefficiens USDA110 may be due to the requirement of more protein modifications, and we further identified 36 conserved functional PPI modules. Integrated with transcriptome and proteome data, interactomes representing free-living cell and symbiotic nitrogen-fixing (SNF) bacteroid were obtained. Based on the SNF interactome, a core-sub-PPI-network for symbiotic nitrogen fixation was determined and nine novel functional modules and eleven key protein hubs playing key roles in symbiosis were identified. The reconstructed interactome of B. diazoefficiens USDA110 may serve as a valuable reference for studying the mechanism underlying the SNF system of rhizobia and legumes.

Imputing single-cell RNA-seq data by combining graph convolution and autoencoder neural networks

Single-cell RNA sequencing technology promotes the profiling of single-cell transcriptomes at an unprecedented throughput and resolution. However, in scRNA-seq studies, only a low amount of sequenced mRNA in each cell leads to missing detection for a portion of mRNA molecules, i.e. the dropout problem which hinders various downstream analyses. Therefore, it is necessary to develop robust and effective imputation methods for the increasing scRNA-seq data. In this study, we have developed an imputation method (GraphSCI) to impute the dropout events in scRNA-seq data based on the graph convolution networks. Extensive experiments demonstrated that GraphSCI outperforms other state-of-the-art methods for imputation on both simulated and real scRNA-seq data. Meanwhile, GraphSCI is able to accurately infer gene-to-gene relationships and the inferred gene-to-gene relationships could also provide powerful assistance for imputation dynamically during the training process, which is a key promotion of GraphSCI compared with other imputation algorithms.

Learning Retention Mechanisms and Evolutionary Parameters of Duplicate Genes from Their Expression Data

Learning about the roles that duplicate genes play in the origins of novel phenotypes requires an understanding of how their functions evolve. A previous method for achieving this goal, CDROM, employs gene expression distances as proxies for functional divergence and then classifies the evolutionary mechanisms retaining duplicate genes from comparisons of these distances in a decision tree framework. However, CDROM does not account for stochastic shifts in gene expression or leverage advances in contemporary statistical learning for performing classification, nor is it capable of predicting the parameters driving duplicate gene evolution. Thus, here we develop CLOUD, a multi-layer neural network built on a model of gene expression evolution that can both classify duplicate gene retention mechanisms and predict their underlying evolutionary parameters. We show that not only is the CLOUD classifier substantially more powerful and accurate than CDROM, but that it also yields accurate parameter predictions, enabling a better understanding of the specific forces driving the evolution and long-term retention of duplicate genes. Further, application of the CLOUD classifier and predictor to empirical data from Drosophila recapitulates many previous findings about gene duplication in this lineage, showing that new functions often emerge rapidly and asymmetrically in younger duplicate gene copies, and that functional divergence is driven by strong natural selection. Hence, CLOUD represents a major advancement in classifying retention mechanisms and predicting evolutionary parameters of duplicate genes, thereby highlighting the utility of incorporating sophisticated statistical learning techniques to address long-standing questions about evolution after gene duplication.

Development and validation of a novel pseudogene pair-based prognostic signature for prediction of overall survival in patients with hepatocellular carcinoma

BACKGROUND

There is growing evidence that pseudogenes may serve as prognostic biomarkers in several cancers. The present study was designed to develop and validate an accurate and robust pseudogene pairs-based signature for the prognosis of hepatocellular carcinoma (HCC).

METHODS

RNA-sequencing data from 374 HCC patients with clinical follow-up information were obtained from the Cancer Genome Atlas (TCGA) database and used in this study. Survival-related pseudogene pairs were identified, and a signature model was constructed by Cox regression analysis (univariate and least absolute shrinkage and selection operator). All individuals were classified into high- and low-risk groups based on the optimal cutoff. Subgroups analysis of the novel signature was conducted and validated in an independent cohort. Pearson correlation analyses were carried out between the included pseudogenes and the protein-coding genes based on their expression levels. Enrichment analysis was performed to predict the possible role of the pseudogenes identified in the signature.

RESULTS

A 19-pseudogene pair signature, which included 21 pseudogenes, was established. Patients in high-risk group demonstrated an increased the risk of adverse prognosis in the TCGA cohort and the external cohort (all P < 0.001). The novel pseudogene signature was independent of other conventional clinical variables used for survival prediction in HCC patients in the two cohorts revealed by the multivariate Cox regression analysis (all P < 0.001). Subgroup analysis further demonstrated the diagnostic value of the signature across different stages, grades, sexes, and age groups. The C-index of the prognostic signature was 0.761, which was not only higher than that of several previous risk models but was also much higher than that of a single age, sex, grade, and stage risk model. Furthermore, functional analysis revealed that the potential biological mechanisms mediated by these pseudogenes are primarily involved in cytokine receptor activity, T cell receptor signaling, chemokine signaling, NF-κB signaling, PD-L1 expression, and the PD-1 checkpoint pathway in cancer.

CONCLUSION

The novel proposed and validated pseudogene pair-based signature may serve as a valuable independent prognostic predictor for predicting survival of patients with HCC.

A pathogenesis-related protein GmPR08-Bet VI promotes a molecular interaction between the GmSHMT08 and GmSNAP18 in resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is the most devastating pest affecting soybean production worldwide. SCN resistance requires both the GmSHMT08 and the GmSNAP18 in 'Peking'-type resistance. Here, we describe the molecular interaction between GmSHMT08 and GmSNAP18, which is potentiated by a pathogenesis-related protein GmPR08-Bet VI. Like GmSNAP18 and GmSHMT08, GmPR08-Bet VI expression was induced in response to SCN and its overexpression decreased SCN cysts by 65% in infected transgenic soybean roots. Overexpression of GmPR08-Bet VI did not have an effect on SCN resistance when the two cytokinin-binding sites in GmPR08-Bet VI were mutated, indicating a new role of GmPR08-Bet VI in SCN resistance. GmPR08-Bet VI was mapped to a QTL for resistance to SCN using different mapping populations. GmSHMT08, GmSNAP18 and GmPR08-Bet VI localize to the cytosol and plasma membrane. GmSNAP18 expression and localization hyper-accumulated at the plasma membrane and was specific to the root cells surrounding the nematode in SCN-resistant soybeans. Genes encoding key components of the salicylic acid signalling pathway were induced under SCN infection. GmSNAP18 and GmPR08-Bet VI were also induced under salicylic acid and cytokinin exogenous treatments, while GmSHMT08 was induced only when the resistant GmSNAP18 was present, pointing to the presence of a molecular crosstalk between SCN-resistant genes and defence genes. Expression analysis of GmSHMT08 and GmSNAP18 identified the need of a minimum expression requirement to trigger the SCN resistance reaction. These results provide insight into a new response mechanism towards plant nematode resistance involving haplotype compatibility, gene dosage and hormone signalling.

Identification of differentially expressed long non-coding RNAs associated with dilated cardiomyopathy using integrated bioinformatics approaches

The aim of this study was to identify novel long non-coding RNA (lncRNA) biomarkers associated with dilated cardiomyopathy (DCM) and reveal the potential molecular mechanisms of DCM development using bioinformatics approaches. The array data of GSE5406, including 108 DCM samples and 16 non-failing control samples, were obtained from the Gene Expression Omnibus database. The differentially expressed lncRNAs were identified using limma package in R. Pearson's correlation analyses were performed between the differentially expressed lncRNAs and protein-coding genes based on their expression levels. Pathway enrichment of these lncRNAs was conducted based on the significantly co-expressed genes. From the receiver operating characteristic (ROC) curve, the area under the ROC curve (AUC) value was obtained and used for evaluating discriminatory ability. IDI2-AS1 and XIST were differentially expressed in DCM patients. A total of 510 co-expressed genes were identified. The enriched functions and pathways of the co-expressed genes mainly included NADH dehydrogenase activity, cardiac muscle contraction, and oxidative phosphorylation. The ROC curve analysis indicated that the two lncRNAs have favorable diagnostic values in DCM. The AUC values of XIST, IDI2-AS1, and the combination of XIST and IDI2-AS1 were 0.733 (95% CI: 0.646-0.809), 0.796 (95% CI: 0.715-0.863), and 0.823 (95% CI: 0.745-0.886), respectively. This study identified IDI2-AS1 and XIST lncRNAs and related pathways involved in the pathogenesis of DCM, thus providing potential diagnostic and therapeutic targets for DCM.

Microarray Analysis of Differential Gene Expression in Alzheimer's Disease Identifies Potential Biomarkers with Diagnostic Value

BACKGROUND Alzheimer disease (AD) is a common and fatal subtype of dementia that remains a challenge to diagnose and treat. This study aimed to identify potential biomarkers that influence the prognosis of AD. MATERIAL AND METHODS A total of 6 gene expression profiles from the Gene Expression Omnibus (GEO) database were assessed for their potential as AD biomarkers. We identified differentially expressed genes (DEGs) using the prediction analysis for microarray (PAM) algorithm and obtained hub genes through the analysis of the protein-protein interaction (PPI) network and module analysis. RESULTS We identified 6 gene expression profiles from the GEO database and assessed their potential as AD biomarkers. Shared gene sets were extracted and integrated into large expression profile matrices. We identified 2514 DEGs including 68 upregulated- and 2446 downregulated genes through analysis of the limma package. We screened 379 significant DEGs including 68 upregulated and 307 downregulated genes for their ability to distinguish AD from control samples using PAM algorithm. Functional enrichment of the 379 target genes was produced from Database for Annotation, Visualization and Integrated Discovery.(DAVID) and included histone function, beta receptor signaling, cell growth, and angiogenesis. The downregulated genes were significantly enriched in MAPK signaling, synaptic signaling, neuronal apoptosis and AD associated pathways. Upon analysis of the PPI network, 32 hub genes including ENO2, CCT2, CALM2, ACACB, ATP5B, MDH1, and PP2CA were screened. Of these hub genes, NFKBIA and ACACB were upregulated and 29 genes were downregulated in AD patients. CONCLUSIONS We screened 379 significant DEGs as potential biomarkers of AD using PAM and obtained 32 hub genes through PPI network and module analysis. These findings reveal new potential AD biomarkers with prognostic and therapeutic value.

Evaluation of colorectal cancer subtypes and cell lines using deep learning

Colorectal cancer (CRC) is a common cancer with a high mortality rate and a rising incidence rate in the developed world. Molecular profiling techniques have been used to better understand the variability between tumors and disease models such as cell lines. To maximize the translatability and clinical relevance of in vitro studies, the selection of optimal cancer models is imperative. We have developed a deep learning-based method to measure the similarity between CRC tumors and disease models such as cancer cell lines. Our method efficiently leverages multiomics data sets containing copy number alterations, gene expression, and point mutations and learns latent factors that describe data in lower dimensions. These latent factors represent the patterns that are clinically relevant and explain the variability of molecular profiles across tumors and cell lines. Using these, we propose refined CRC subtypes and provide best-matching cell lines to different subtypes. These findings are relevant to patient stratification and selection of cell lines for early-stage drug discovery pipelines, biomarker discovery, and target identification.

Yeast two-hybrid screening for proteins that interact with PFT in wheat

Fusarium head blight (FHB) is a devastating disease of wheat worldwide. Fhb1 is the most consistently reported quantitative trait locus (QTL) for FHB resistance breeding. A pore-forming toxin-like (PFT) gene at Fhb1 was first cloned by map-based cloning and found to confer FHB resistance in wheat. Proteins often interact with each other to execute their functions. Characterization of the proteins interacting with PFT might therefore provide information on the molecular mechanisms of PFT functions. In this study, a high-quality yeast two-hybrid (Y2H) library using RNA extracted from Fusarium graminearum (Fg)-infected wheat spikes of Sumai 3 was constructed. The agglutinin domains of PFT exhibited no self-activation and toxicity to yeast cells and were used as bait to screen the Y2H library. Twenty-three proteins that interact with PFT were obtained, which were mainly involved in the ubiquitination process, clathrin coat assembly, the oxidation-reduction process, and protein phosphorylation. The expression pattern of these interacting genes was analyzed by quantitative real-time PCR. This study clarifies the protein interactions of PFT and raises a regulatory network for PFT regarding FHB resistance in wheat.

Network Medicine in the Age of Biomedical Big Data

Network medicine is an emerging area of research dealing with molecular and genetic interactions, network biomarkers of disease, and therapeutic target discovery. Large-scale biomedical data generation offers a unique opportunity to assess the effect and impact of cellular heterogeneity and environmental perturbations on the observed phenotype. Marrying the two, network medicine with biomedical data provides a framework to build meaningful models and extract impactful results at a network level. In this review, we survey existing network types and biomedical data sources. More importantly, we delve into ways in which the network medicine approach, aided by phenotype-specific biomedical data, can be gainfully applied. We provide three paradigms, mainly dealing with three major biological network archetypes: protein-protein interaction, expression-based, and gene regulatory networks. For each of these paradigms, we discuss a broad overview of philosophies under which various network methods work. We also provide a few examples in each paradigm as a test case of its successful application. Finally, we delineate several opportunities and challenges in the field of network medicine. We hope this review provides a lexicon for researchers from biological sciences and network theory to come on the same page to work on research areas that require interdisciplinary expertise. Taken together, the understanding gained from combining biomedical data with networks can be useful for characterizing disease etiologies and identifying therapeutic targets, which, in turn, will lead to better preventive medicine with translational impact on personalized healthcare.

Network-based methods for predicting essential genes or proteins: a survey

Genes that are thought to be critical for the survival of organisms or cells are called essential genes. The prediction of essential genes and their products (essential proteins) is of great value in exploring the mechanism of complex diseases, the study of the minimal required genome for living cells and the development of new drug targets. As laboratory methods are often complicated, costly and time-consuming, a great many of computational methods have been proposed to identify essential genes/proteins from the perspective of the network level with the in-depth understanding of network biology and the rapid development of biotechnologies. Through analyzing the topological characteristics of essential genes/proteins in protein–protein interaction networks (PINs), integrating biological information and considering the dynamic features of PINs, network-based methods have been proved to be effective in the identification of essential genes/proteins. In this paper, we survey the advanced methods for network-based prediction of essential genes/proteins and present the challenges and directions for future research.

Lineage-Specific Expression Divergence in Grasses Is Associated with Male Reproduction, Host-Pathogen Defense, and Domestication

Poaceae (grasses) is an agriculturally important and widely distributed family of plants with extraordinary phenotypic diversity, much of which was generated under recent lineage-specific evolution. Yet, little is known about the genes and functional modules involved in the lineage-specific divergence of grasses. Here, I address this question on a genome-wide scale by applying a novel branch-based statistic of lineage-specific expression divergence, LED, to RNA-seq data from nine tissues of the wild grass Brachypodium distachyon and its domesticated relatives Oryza sativa japonica (rice) and Sorghum bicolor (sorghum). I find that LED is generally smallest in B. distachyon and largest in O. sativa japonica, which underwent domestication earlier than S. bicolor, supporting the hypothesis that domestication may increase the rate of lineage-specific expression divergence in grasses. Moreover, in all three species, LED is positively correlated with protein-coding sequence divergence and tissue specificity, and negatively correlated with network connectivity. Further analysis reveals that genes with large LED are often primarily expressed in anther, implicating lineage-specific expression divergence in the evolution of male reproductive phenotypes. Gene ontology enrichment analysis also identifies an overrepresentation of terms related to male reproduction in the two domesticated grasses, as well as to those involved in host-pathogen defense in all three species. Last, examinations of genes with the largest LED reveal that their lineage-specific expression divergence may have contributed to antimicrobial functions in B. distachyon, to enhanced adaptation and yield during domestication in O. sativa japonica, and to defense against a widespread and devastating fungal pathogen in S. bicolor. Together, these findings suggest that lineage-specific expression divergence in grasses may increase under domestication and preferentially target rapidly evolving genes involved in male reproduction, host-pathogen defense, and the origin of domesticated phenotypes.

PPIExp: A Web-Based Platform for Integration and Visualization of Protein-Protein Interaction Data and Spatiotemporal Proteomics Data

Integrating spatiotemporal proteomics data with protein-protein interaction (PPI) data can help researchers make an in-depth exploration of their proteins of interest in a dynamic manner. However, there is still a lack of proper tools for the biologists who usually have few programming skills to construct a PPI network for a protein list, visualize active PPI subnetworks, and then select key nodes for further study. We propose a web-based platform named PPIExp that can automatically construct a PPI network, perform clustering analysis according to protein abundances, and perform functional enrichment analysis. More importantly, it provides multiple effective visualization interfaces, such as the interface to display the PPI network map, the interface to display a dendrogram and heatmap for the clustering result, and the interface to display the expression pattern of a selected protein. To visualize the active PPI subnetworks in specific space or time, it provides buttons to highlight the differentially expressed proteins under each condition on the network map. Additionally, to help researchers determine which proteins are worth further attention, PPIExp provides extensive one-click interactive operations to map node centrality measures to node size on the network and highlight three types of proteins, that is, the proteins in an enriched functional term, the coexpressed proteins selected from the dendgrogram and heatmap, and the proteins input by users. PPIExp is available at http://www.fgvis.com/expressvis/PPIExp .

Computational Methods for Predicting Protein-Protein Interactions Using Various Protein Features

Understanding protein-protein interactions (PPIs) in a cell is essential for learning protein functions, pathways, and mechanism of diseases. PPIs are also important targets for developing drugs. Experimental methods, both small-scale and large-scale, have identified PPIs in several model organisms. However, results cover only a part of PPIs of organisms; moreover, there are many organisms whose PPIs have not yet been investigated. To complement experimental methods, many computational methods have been developed that predict PPIs from various characteristics of proteins. Here we provide an overview of literature reports to classify computational PPI prediction methods that consider different features of proteins, including protein sequence, genomes, protein structure, function, PPI network topology, and those which integrate multiple methods. © 2018 by John Wiley & Sons, Inc.

Reconstruction of the Protein-Protein Interaction Network for Protein Complexes Identification by Walking on the Protein Pair Fingerprints Similarity Network

Identifying protein complexes from protein-protein interaction networks (PPINs) is important to understand the science of cellular organization and function. However, PPINs produced by high-throughput studies have high false discovery rate and only represent snapshot interaction information. Reconstructing higher quality PPINs is essential for protein complex identification. Here we present a Multi-Level PPINs reconstruction (MLPR) method for protein complexes detection. From existing PPINs, we generated full combinations of every two proteins. These protein pairs are represented as a vector which includes six different sources. Then the protein pairs with same vector are mapped to the same fingerprint ID. A fingerprint similarity network is constructed next, in which a vertex represents a protein pair fingerprint ID and each vertex is connected to its top 10 similar fingerprints by edges. After random walking on the fingerprints similarity network, each vertex got a score at the steady state. According to the score of protein pairs, we considered the top ranked ones as reliable PPI and the score as the weight of edge between two distinct proteins. Finally, we expanded clusters starting from seeded vertexes based on the new weighted reliable PPINs. Applying our method on the yeast PPINs, our algorithm achieved higher F-value in protein complexes detection than the-state-of-the-art methods. The interactions in our reconstructed PPI network have more significant biological relevance than the exiting PPI datasets, assessed by gene ontology. In addition, the performance of existing popular protein complexes detection methods are significantly improved on our reconstructed network.

Estrogen-related receptors alpha, beta and gamma expression and function is associated with transcriptional repressor EZH2 in breast carcinoma

BACKGROUND

Orphan nuclear receptors ERRα, ERRβ and ERRγ that belong to NR3B or type IV nuclear receptor family are well studied for their role in breast cancer pathophysiology. Their homology with the canonical estrogen receptor dictates their possible contributing role in mammary gland development and disease. Although function and regulation of ERRα, ERRγ and less about ERRβ is reported, role of histone methylation in their altered expression in cancer cells is not studied. Transcriptional activity of nuclear receptors depends on co-regulatory proteins. The present study for the first time gives an insight into regulation of estrogen-related receptors by histone methylation specifically through methyltransferase EZH2 in breast cancer.

METHODS

Expression of ERRα, ERRβ, ERRγ and EZH2 was assessed by immunohistochemistry in four identical tissue array slides that were prepared as per the protocol. The array slides were stained with ERRα, ERRβ, ERRγ and EZH2 simultaneously. Array data was correlated with expression in MERAV expression dataset. Pearson correlation coeficient r was calculated from the partial matrix expression values available at MERAV database to study the strength of association between EZH2 and three orphan nuclear receptors under study. By western blot and real time PCR, their correlated expression was studied in breast cancer cell lines MCF-7, MDA-MB-231, T47D and MDA-MB-453 including normal breast epithelial MCF-10A cells at both protein and RNA level. Regulation of ERRα, ERRβ, ERRγ by EZH2 was further investigated upon overexpression and silencing of EZH2. The interaction between ERRs and EZH2 was validated in vivo by CHIP-qPCR.

RESULTS

We found a negative correlation between estrogen-related receptors and Enhancer of Zeste Homolog 2, a global repressor gene. Immunohistochemistry in primary breast tumors of different grades showed a correlated expression of estrogen-related receptors and EZH2. Their correlated expression was further validated using online MERAV expression dataset where a negative correlation of variable strengths was observed in breast cancer. Ectopic expression of EZH2 in low EZH2-expressing normal breast epithelial cells abrogated their expression and at the same time, its silencing enhanced the expression of estrogen-related receptors in cancerous cells. Global occupancy of EZH2 on ERRα and ERRβ was observed in-vivo.

CONCLUSION

Our findings identify EZH2 as a relevant coregulator for estrogen-related receptors in breast carcinoma.

Elevated CO2 induces alteration in lignin accumulation in celery (Apium graveolens L.)

Carbon dioxide (CO₂) is an important regulator of plant growth and development, and its proportion in the atmosphere continues to rise now. Lignin is one of the major secondary products in plants with vital biological functions. However, the relationship between CO₂ level and xylogenesis in celery is still unknown. In order to investigate the effects of increasing CO₂ concentration on lignin accumulation in celery, 'Jinnanshiqin' were exposed to two CO₂ applications, 400 (e₀) and 1000 μmol mol^-1 (e₁), respectively. Plant morphology and lignin distribution in celery plants treated with elevated CO₂ did not change significantly. There was an upward trend on lignin content in celery leaves, and the transcript abundance of 12 genes involved in lignin metabolism has altered in response to elevated CO₂. The effects of high level of CO₂ on different tissues were different. Our works confirmed that CO₂ may play an important role in lignin accumulation in celery leaves. The current study will offer new evidence to understand the regulation mechanism of lignin biosynthesis under elevated CO₂ and provide a reference to improve celery quality by adjusting the growth environment.

netSmooth: Network-smoothing based imputation for single cell RNA-seq

Single cell RNA-seq (scRNA-seq) experiments suffer from a range of characteristic technical biases, such as dropouts (zero or near zero counts) and high variance. Current analysis methods rely on imputing missing values by various means of local averaging or regression, often amplifying biases inherent in the data. We present netSmooth, a network-diffusion based method that uses priors for the covariance structure of gene expression profiles on scRNA-seq experiments in order to smooth expression values. We demonstrate that netSmooth improves clustering results of scRNA-seq experiments from distinct cell populations, time-course experiments, and cancer genomics. We provide an R package for our method, available at: https://github.com/BIMSBbioinfo/netSmooth.

netSmooth: Network-smoothing based imputation for single cell RNA-seq

Single cell RNA-seq (scRNA-seq) experiments suffer from a range of characteristic technical biases, such as dropouts (zero or near zero counts) and high variance. Current analysis methods rely on imputing missing values by various means of local averaging or regression, often amplifying biases inherent in the data. We present netSmooth, a network-diffusion based method that uses priors for the covariance structure of gene expression profiles on scRNA-seq experiments in order to smooth expression values. We demonstrate that netSmooth improves clustering results of scRNA-seq experiments from distinct cell populations, time-course experiments, and cancer genomics. We provide an R package for our method, available at: https://github.com/BIMSBbioinfo/netSmooth.

netSmooth: Network-smoothing based imputation for single cell RNA-seq

Single cell RNA-seq (scRNA-seq) experiments suffer from a range of characteristic technical biases, such as dropouts (zero or near zero counts) and high variance. Current analysis methods rely on imputing missing values by various means of local averaging or regression, often amplifying biases inherent in the data. We present netSmooth, a network-diffusion based method that uses priors for the covariance structure of gene expression profiles on scRNA-seq experiments in order to smooth expression values. We demonstrate that netSmooth improves clustering results of scRNA-seq experiments from distinct cell populations, time-course experiments, and cancer genomics. We provide an R package for our method, available at: https://github.com/BIMSBbioinfo/netSmooth.

Data integration and predictive modeling methods for multi-omics datasets

We provide an overview of opportunities and challenges in multi-omics predictive analytics with particular emphasis on data integration and machine learning methods.

Natural Selection Drives Rapid Functional Evolution of Young Drosophila Duplicate Genes

Gene duplication is thought to play a major role in phenotypic evolution. Yet the forces involved in the functional divergence of young duplicate genes remain unclear. Here, we use population-genetic inference to elucidate the role of natural selection in the functional evolution of young duplicate genes in Drosophila melanogaster. We find that negative selection acts on young duplicates with ancestral functions, and positive selection on those with novel functions, suggesting that natural selection may determine whether and how young duplicate genes are retained. Moreover, evidence of natural selection is strongest in protein-coding regions and 3' UTRs of young duplicates, indicating that selection may primarily target encoded proteins and regulatory sequences specific to 3' UTRs. Further analysis reveals that natural selection acts immediately after duplication and weakens over time, possibly explaining the observed bias toward the acquisition of new functions by young, rather than old, duplicate gene copies. Last, we find an enrichment of testis-related functions in young duplicates that underwent recent positive selection, but not in young duplicates that did not undergo recent positive selection, or in old duplicates that either did or did not undergo recent positive selection. Thus, our findings reveal that natural selection is a key player in the functional evolution of young duplicate genes, acts rapidly and in a region-specific manner, and may underlie the origin of novel testis-specific phenotypes in Drosophila.

GPCRs globally coevolved with receptor activity-modifying proteins, RAMPs

Receptor activity-modifying proteins (RAMPs) are widely expressed in human tissues and, in some cases, have been shown to affect surface expression or ligand specificity of G-protein-coupled receptors (GPCRs). However, whether RAMP-GPCR interactions are widespread, and the nature of their functional consequences, remains largely unknown. In humans, there are three RAMPs and over 800 expressed GPCRs, making direct experimental approaches challenging. We analyzed relevant genomic data from all currently available sequenced organisms. We discovered that RAMPs and GPCRs tend to have orthologs in the same species and have correlated phylogenetic trees to the same extent, or higher than other interacting protein pairs that play key roles in cellular signaling. In addition, the resulting RAMP-GPCR interaction map suggests that RAMP1 and RAMP3 interact with the same set of GPCRs, which implies functional redundancy. We next analyzed human transcriptomes and found expression correlation for GPCRs and RAMPs. Our results suggest global coevolution of GPCRs and RAMPS and support the hypothesis that GPCRs interact globally with RAMPs in cellular signaling pathways.

Protease-Inhibitor Interaction Predictions: Lessons on the Complexity of Protein-Protein Interactions

Protein interactions shape proteome function and thus biology. Identification of protein interactions is a major goal in molecular biology, but biochemical methods, although improving, remain limited in coverage and accuracy. Whereas computational predictions can guide biochemical experiments, low validation rates of predictions remain a major limitation. Here, we investigated computational methods in the prediction of a specific type of interaction, the inhibitory interactions between proteases and their inhibitors. Proteases generate thousands of proteoforms that dynamically shape the functional state of proteomes. Despite the important regulatory role of proteases, knowledge of their inhibitors remains largely incomplete with the vast majority of proteases lacking an annotated inhibitor. To link inhibitors to their target proteases on a large scale, we applied computational methods to predict inhibitory interactions between proteases and their inhibitors based on complementary data, including coexpression, phylogenetic similarity, structural information, co-annotation, and colocalization, and also surveyed general protein interaction networks for potential inhibitory interactions. In testing nine predicted interactions biochemically, we validated the inhibition of kallikrein 5 by serpin B12. Despite the use of a wide array of complementary data, we found a high false positive rate of computational predictions in biochemical follow-up. Based on a protease-specific definition of true negatives derived from the biochemical classification of proteases and inhibitors, we analyzed prediction accuracy of individual features, thereby we identified feature-specific limitations, which also affected general protein interaction prediction methods. Interestingly, proteases were often not coexpressed with most of their functional inhibitors, contrary to what is commonly assumed and extrapolated predominantly from cell culture experiments. Predictions of inhibitory interactions were indeed more challenging than predictions of nonproteolytic and noninhibitory interactions. In summary, we describe a novel and well-defined but difficult protein interaction prediction task and thereby highlight limitations of computational interaction prediction methods.

HAPPI-2: a Comprehensive and High-quality Map of Human Annotated and Predicted Protein Interactions

BACKGROUND

Human protein-protein interaction (PPI) data is essential to network and systems biology studies. PPI data can help biochemists hypothesize how proteins form complexes by binding to each other, how extracellular signals propagate through post-translational modification of de-activated signaling molecules, and how chemical reactions are coupled by enzymes involved in a complex biological process. Our capability to develop good public database resources for human PPI data has a direct impact on the quality of future research on genome biology and medicine.

RESULTS

The database of Human Annotated and Predicted Protein Interactions (HAPPI) version 2.0 is a major update to the original HAPPI 1.0 database. It contains 2,922,202 unique protein-protein interactions (PPI) linked by 23,060 human proteins, making it the most comprehensive database covering human PPI data today. These PPIs contain both physical/direct interactions and high-quality functional/indirect interactions. Compared with the HAPPI 1.0 database release, HAPPI database version 2.0 (HAPPI-2) represents a 485% of human PPI data coverage increase and a 73% protein coverage increase. The revamped HAPPI web portal provides users with a friendly search, curation, and data retrieval interface, allowing them to retrieve human PPIs and available annotation information on the interaction type, interaction quality, interacting partner drug targeting data, and disease information. The updated HAPPI-2 can be freely accessed by Academic users at http://discovery.informatics.uab.edu/HAPPI .

CONCLUSIONS

While the underlying data for HAPPI-2 are integrated from a diverse data sources, the new HAPPI-2 release represents a good balance between data coverage and data quality of human PPIs, making it ideally suited for network biology.

Personalized characterization of diseases using sample-specific networks

A complex disease generally results not from malfunction of individual molecules but from dysfunction of the relevant system or network, which dynamically changes with time and conditions. Thus, estimating a condition-specific network from a single sample is crucial to elucidating the molecular mechanisms of complex diseases at the system level. However, there is currently no effective way to construct such an individual-specific network by expression profiling of a single sample because of the requirement of multiple samples for computing correlations. We developed here with a statistical method, i.e. a sample-specific network (SSN) method, which allows us to construct individual-specific networks based on molecular expressions of a single sample. Using this method, we can characterize various human diseases at a network level. In particular, such SSNs can lead to the identification of individual-specific disease modules as well as driver genes, even without gene sequencing information. Extensive analysis by using the Cancer Genome Atlas data not only demonstrated the effectiveness of the method, but also found new individual-specific driver genes and network patterns for various types of cancer. Biological experiments on drug resistance further validated one important advantage of our method over the traditional methods, i.e. we can even identify such drug resistance genes that actually have no clear differential expression between samples with and without the resistance, due to the additional network information.

A Network Approach of Gene Co-expression in the Zea mays/Aspergillus flavus Pathosystem to Map Host/Pathogen Interaction Pathways

A gene co-expression network (GEN) was generated using a dual RNA-seq study with the fungal pathogen Aspergillus flavus and its plant host Zea mays during the initial 3 days of infection. The analysis deciphered novel pathways and mapped genes of interest in both organisms during the infection. This network revealed a high degree of connectivity in many of the previously recognized pathways in Z. mays such as jasmonic acid, ethylene, and reactive oxygen species (ROS). For the pathogen A. flavus, a link between aflatoxin production and vesicular transport was identified within the network. There was significant interspecies correlation of expression between Z. mays and A. flavus for a subset of 104 Z. mays, and 1942 A. flavus genes. This resulted in an interspecies subnetwork enriched in multiple Z. mays genes involved in the production of ROS. In addition to the ROS from Z. mays, there was enrichment in the vesicular transport pathways and the aflatoxin pathway for A. flavus. Included in these genes, a key aflatoxin cluster regulator, AflS, was found to be co-regulated with multiple Z. mays ROS producing genes within the network, suggesting AflS may be monitoring host ROS levels. The entire GEN for both host and pathogen, and the subset of interspecies correlations, is presented as a tool for hypothesis generation and discovery for events in the early stages of fungal infection of Z. mays by A. flavus.

Network-Guided Key Gene Discovery for a Given Cellular Process

Identification of key genes for a given physiological or pathological process is an essential but still very challenging task for the entire biomedical research community. Statistics-based approaches, such as genome-wide association study (GWAS)- or quantitative trait locus (QTL)-related analysis have already made enormous contributions to identifying key genes associated with a given disease or phenotype, the success of which is however very much dependent on a huge number of samples. Recent advances in network biology, especially network inference directly from genome-scale data and the following-up network analysis, opens up new avenues to predict key genes driving a given biological process or cellular function. Here we review and compare the current approaches in predicting key genes, which have no chances to stand out by classic differential expression analysis, from gene-regulatory, protein-protein interaction, or gene expression correlation networks. We elaborate these network-based approaches mainly in the context of immunology and infection, and urge more usage of correlation network-based predictions. Such network-based key gene discovery approaches driven by information-enriched 'omics' data should be very useful for systematic key gene discoveries for any given biochemical process or cellular function, and also valuable for novel drug target discovery and novel diagnostic, prognostic and therapeutic-efficiency marker prediction for a specific disease or disorder.

A computational interactome and functional annotation for the human proteome

We present a database, PrePPI (Predicting Protein-Protein Interactions), of more than 1.35 million predicted protein-protein interactions (PPIs). Of these at least 127,000 are expected to constitute direct physical interactions although the actual number may be much larger (~500,000). The current PrePPI, which contains predicted interactions for about 85% of the human proteome, is related to an earlier version but is based on additional sources of interaction evidence and is far larger in scope. The use of structural relationships allows PrePPI to infer numerous previously unreported interactions. PrePPI has been subjected to a series of validation tests including reproducing known interactions, recapitulating multi-protein complexes, analysis of disease associated SNPs, and identifying functional relationships between interacting proteins. We show, using Gene Set Enrichment Analysis (GSEA), that predicted interaction partners can be used to annotate a protein's function. We provide annotations for most human proteins, including many annotated as having unknown function.

Predicting physiologically relevant SH3 domain mediated protein-protein interactions in yeast

MOTIVATION

Many intracellular signaling processes are mediated by interactions involving peptide recognition modules such as SH3 domains. These domains bind to small, linear protein sequence motifs which can be identified using high-throughput experimental screens such as phage display. Binding motif patterns can then be used to computationally predict protein interactions mediated by these domains. While many protein-protein interaction prediction methods exist, most do not work with peptide recognition module mediated interactions or do not consider many of the known constraints governing physiologically relevant interactions between two proteins.

RESULTS

A novel method for predicting physiologically relevant SH3 domain-peptide mediated protein-protein interactions in S. cerevisae using phage display data is presented. Like some previous similar methods, this method uses position weight matrix models of protein linear motif preference for individual SH3 domains to scan the proteome for potential hits and then filters these hits using a range of evidence sources related to sequence-based and cellular constraints on protein interactions. The novelty of this approach is the large number of evidence sources used and the method of combination of sequence based and protein pair based evidence sources. By combining different peptide and protein features using multiple Bayesian models we are able to predict high confidence interactions with an overall accuracy of 0.97.

AVAILABILITY AND IMPLEMENTATION

Domain-Motif Mediated Interaction Prediction (DoMo-Pred) command line tool and all relevant datasets are available under GNU LGPL license for download from http://www.baderlab.org/Software/DoMo-Pred The DoMo-Pred command line tool is implemented using Python 2.7 and C ++.

CONTACT

gary.bader@utoronto.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global De Novo Protein-Protein Interactome Elucidates Interactions of Drought-Responsive Proteins in Horse Gram (Macrotyloma uniflorum)

Inspired by the availability of de novo transcriptome of horse gram (Macrotyloma uniflorum) and recent developments in systems biology studies, the first ever global protein-protein interactome (PPI) map was constructed for this highly drought-tolerant legume. Large-scale studies of PPIs and the constructed database would provide rationale behind the interplay at cascading translational levels for drought stress-adaptive mechanisms in horse gram. Using a bidirectional approach (interolog and domain-based), a high-confidence interactome map and database for horse gram was constructed. Available transcriptomic information for shoot and root tissues of a sensitive (M-191; genotype 1) and a drought-tolerant (M-249; genotype 2) genotype of horse gram was utilized to draw comparative PPI subnetworks under drought stress. High-confidence 6804 interactions were predicted among 1812 proteins covering about one-fourth of the horse gram proteome. The highest number of interactions (33.86%) in horse gram interactome matched with Arabidopsis PPI data. The top five hub nodes mostly included ubiquitin and heat-shock-related proteins. Higher numbers of PPIs were found to be responsive in shoot tissue (416) and root tissue (2228) of genotype 2 compared with shoot tissue (136) and root tissue (579) of genotype 1. Characterization of PPIs using gene ontology analysis revealed that kinase and transferase activities involved in signal transduction, cellular processes, nucleocytoplasmic transport, protein ubiquitination, and localization of molecules were most responsive to drought stress. Hence, these could be framed in stress adaptive mechanisms of horse gram. Being the first legume global PPI map, it would provide new insights into gene and protein regulatory networks for drought stress tolerance mechanisms in horse gram. Information compiled in the form of database (MauPIR) will provide the much needed high-confidence systems biology information for horse gram genes, proteins, and involved processes. This information would ease the effort and increase the efficacy for similar studies on other legumes. Public access is available at http://14.139.59.221/MauPIR/ .

Investigating the Combinatory Effects of Biological Networks on Gene Co-expression

Co-expressed genes often share similar functions, and gene co-expression networks have been widely used in studying the functionality of gene modules. Previous analysis indicated that genes are more likely to be co-expressed if they are either regulated by the same transcription factors, forming protein complexes or sharing similar topological properties in protein-protein interaction networks. Here, we reconstructed transcriptional regulatory and protein-protein networks for Saccharomyces cerevisiae using well-established databases, and we evaluated their co-expression activities using publically available gene expression data. Based on our network-dependent analysis, we found that genes that were co-regulated in the transcription regulatory networks and shared similar neighbors in the protein-protein networks were more likely to be co-expressed. Moreover, their biological functions were closely related.

Inferring interaction type in gene regulatory networks using co-expression data

Background

Knowledge of interaction types in biological networks is important for understanding the functional organization of the cell. Currently information-based approaches are widely used for inferring gene regulatory interactions from genomics data, such as gene expression profiles; however, these approaches do not provide evidence about the regulation type (positive or negative sign) of the interaction.

Results

This paper describes a novel algorithm, “Signing of Regulatory Networks” (SIREN), which can infer the regulatory type of interactions in a known gene regulatory network (GRN) given corresponding genome-wide gene expression data. To assess our new approach, we applied it to three different benchmark gene regulatory networks, including Escherichia coli, prostate cancer, and an in silico constructed network. Our new method has approximately 68, 70, and 100 percent accuracy, respectively, for these networks. To showcase the utility of SIREN algorithm, we used it to predict previously unknown regulation types for 454 interactions related to the prostate cancer GRN.

Conclusions

SIREN is an efficient algorithm with low computational complexity; hence, it is applicable to large biological networks. It can serve as a complementary approach for a wide range of network reconstruction methods that do not provide information about the interaction type.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-015-0054-4) contains supplementary material, which is available to authorized users.

A predicted protein interactome identifies conserved global networks and disease resistance subnetworks in maize

Interactomes are genome-wide roadmaps of protein-protein interactions. They have been produced for humans, yeast, the fruit fly, and Arabidopsis thaliana and have become invaluable tools for generating and testing hypotheses. A predicted interactome for Zea mays (PiZeaM) is presented here as an aid to the research community for this valuable crop species. PiZeaM was built using a proven method of interologs (interacting orthologs) that were identified using both one-to-one and many-to-many orthology between genomes of maize and reference species. Where both maize orthologs occurred for an experimentally determined interaction in the reference species, we predicted a likely interaction in maize. A total of 49,026 unique interactions for 6004 maize proteins were predicted. These interactions are enriched for processes that are evolutionarily conserved, but include many otherwise poorly annotated proteins in maize. The predicted maize interactions were further analyzed by comparing annotation of interacting proteins, including different layers of ontology. A map of pairwise gene co-expression was also generated and compared to predicted interactions. Two global subnetworks were constructed for highly conserved interactions. These subnetworks showed clear clustering of proteins by function. Another subnetwork was created for disease response using a bait and prey strategy to capture interacting partners for proteins that respond to other organisms. Closer examination of this subnetwork revealed the connectivity between biotic and abiotic hormone stress pathways. We believe PiZeaM will provide a useful tool for the prediction of protein function and analysis of pathways for Z. mays researchers and is presented in this paper as a reference tool for the exploration of protein interactions in maize.

Constructing minimal models for complex system dynamics

One of the strengths of statistical physics is the ability to reduce macroscopic observations into microscopic models, offering a mechanistic description of a system's dynamics. This paradigm, rooted in Boltzmann's gas theory, has found applications from magnetic phenomena to subcellular processes and epidemic spreading. Yet, each of these advances were the result of decades of meticulous model building and validation, which are impossible to replicate in most complex biological, social or technological systems that lack accurate microscopic models. Here we develop a method to infer the microscopic dynamics of a complex system from observations of its response to external perturbations, allowing us to construct the most general class of nonlinear pairwise dynamics that are guaranteed to recover the observed behaviour. The result, which we test against both numerical and empirical data, is an effective dynamic model that can predict the system's behaviour and provide crucial insights into its inner workings.

Spatial features for Escherichia coli genome organization

Background

In bacterial genomes, the compactly encoded genes and operons are well organized, with genes in the same biological pathway or operons in the same regulon close to each other on the genome sequence. In addition, the linearly close genes have a higher probability of co-expression and their protein products tend to form protein–protein interactions. However, the organization features of bacterial genomes in a three-dimensional space remain elusive. The DNA interaction data of Escherichia coli, measured by the genome conformation capture (GCC) technique, have recently become available, which allowed us to investigate the spatial features of bacterial genome organization.

Results

By renormalizing the GCC data, we compared the interaction frequency of operon pairs in the same regulon with that of random operon pairs. The results showed that arrangements of operons in the E. coli genome tend to minimize the spatial distance between operons in the same regulon. A similar global organization feature exists for genes in biological pathways of E. coli. In addition, the genes close to each other spatially (even if they are far from each other on the genome sequence) tend to be co-expressed and form protein–protein interactions. These results provided new insights into the organization principles of bacterial genomes and support the notion of transcription factory.

Conclusions

This study revealed the organization features of Escherichia coli genomic functional units in the 3D space and furthered our understanding of the link between the three-dimensional structure of chromosomes and biological function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1258-1) contains supplementary material, which is available to authorized users.

A regulatory similarity measure using the location information of transcription factor binding sites in Saccharomyces cerevisiae

Background

Defining a measure for regulatory similarity (RS) of two genes is an important step toward identifying co-regulated genes. To date, transcription factor binding sites (TFBSs) have been widely used to measure the RS of two genes because transcription factors (TFs) binding to TFBSs in promoters is the most crucial and well understood step in gene regulation. However, existing TFBS-based RS measures consider the relation of a TFBS to a gene as a Boolean (either 'presence' or 'absence') without utilizing the information of TFBS locations in promoters.

Results

Functional TFBSs of many TFs in yeast are known to have a strong positional preference to occur in a small region in the promoters. This biological knowledge prompts us to develop a novel RS measure that exploits the TFBS location information. The performances of different RS measures are evaluated by the fraction of gene pairs that are co-regulated (validated by literature evidence) by at least one common TF under different RS scores. The experimental results show that the proposed RS measure is the best co-regulation indicator among the six compared RS measures. In addition, the co-regulated genes identified by the proposed RS measure are also shown to be able to benefit three co-regulation-based applications: detecting gene co-function, gene co-expression and protein-protein interactions.

Conclusions

The proposed RS measure provides a good indicator for gene co-regulation. Besides, its good performance reveals the importance of the location information in TFBS-based RS measures.