Analysis of the human protein interactome and comparison with yeast, worm and fly interaction datasets

Structure-guided Evolutionary Analysis of Interactome Network Rewiring at Single Residue Resolution in Yeasts

Protein-protein interactions (PPIs) are known to rewire extensively during evolution leading to lineage-specific and species-specific changes in molecular processes. However, the detailed molecular evolutionary mechanisms underlying interactome network rewiring are not well-understood. Here, we combine high-confidence PPI data, high-resolution three-dimensional structures of protein complexes, and homology-based structural annotation transfer to construct structurally-resolved interactome networks for the two yeasts S. cerevisiae and S. pombe. We then classify PPIs according to whether they are preserved or different between the two yeast species and compare site-specific evolutionary rates of interfacial versus non-interfacial residues for these different categories of PPIs. We find that residues in PPI interfaces evolve significantly more slowly than non-interfacial residues when using lineage-specific measures of evolutionary rate, but not when using non-lineage-specific measures. Furthermore, both lineage-specific and non-lineage-specific evolutionary rate measures can distinguish interfacial residues from non-interfacial residues for preserved PPIs between the two yeasts, but only the lineage-specific measure is appropriate for rewired PPIs. Finally, both lineage-specific and non-lineage-specific evolutionary rate measures are appropriate for elucidating structural determinants of protein evolution for residues outside of PPI interfaces. Overall, our results demonstrate that unlike tertiary structures of single proteins, PPIs and PPI interfaces can be highly volatile in their evolution, thus requiring the use of lineage-specific measures when studying their evolution. These results yield insight into the evolutionary design principles of PPIs and the mechanisms by which interactions are preserved or rewired between species, improving our understanding of the molecular evolution of PPIs and PPI interfaces at the residue level.

Graph Node Classification to Predict Autism Risk in Genes

This study explores the genetic risk associations with autism spectrum disorder (ASD) using graph neural networks (GNNs), leveraging the Sfari dataset and protein interaction network (PIN) data. We built a gene network with genes as nodes, chromosome band location as node features, and gene interactions as edges. Graph models were employed to classify the autism risk associated with newly introduced genes (test set). Three classification tasks were undertaken to test the ability of our models: binary risk association, multi-class risk association, and syndromic gene association. We tested graph convolutional networks, Graph Sage, graph transformer, and Multi-Layer Perceptron (Baseline) architectures on this problem. The Graph Sage model consistently outperformed the other models, showcasing its utility in classifying ASD-related genes. Our ablation studies show that the chromosome band location and protein interactions contain useful information for this problem. The models achieved 85.80% accuracy on the binary risk classification, 81.68% accuracy on the multi-class risk classification, and 90.22% on the syndromic classification.

Predicting cell-type specific disease genes of diabetes with the biological network

Type 2 diabetes (T2D) is a chronic condition that can lead to significant harm, such as heart disease, kidney disease, nerve damage, and blindness. Although T2D-related genes have been identified through Genome-wide association studies (GWAS) and various computational methods, the biological mechanism of T2D at the cell type level remains unclear. Exploring cell type-specific genes related to T2D is essential to understand the cellular mechanisms underlying the disease. To address this issue, we introduce DiGCellNet (predicting Disease Genes with Cell type specificity based on biological Networks), a model that integrates graph convolutional network (GCN) and multi-task learning (MTL) to predict T2D-associated cell type-specific genes based on the biological network. Our work represents the first attempt to predict cell type-specific disease genes using GCN and MTL. We evaluate our approach by predicting genes specific to four cell types and demonstrate that the proposed DiGCellNet outperforms other models that combine node embeddings with traditional machine learning algorithms. Moreover, DiGCellNet successfully identifies CALM1 as a gene specific to beta cell type in T2D cases, and this association is confirmed using an independent dataset. The code is available at https://github.com/23AIBox/23AIBox-DiGCellNet.

Bioinformatics Prediction for Network-Based Integrative Multi-Omics Expression Data Analysis in Hirschsprung Disease

Hirschsprung's disease (HSCR) is a rare developmental disorder in which enteric ganglia are missing along a portion of the intestine. HSCR has a complex inheritance, with RET as the major disease-causing gene. However, the pathogenesis of HSCR is still not completely understood. Therefore, we applied a computational approach based on multi-omics network characterization and clustering analysis for HSCR-related gene/miRNA identification and biomarker discovery. Protein-protein interaction (PPI) and miRNA-target interaction (MTI) networks were analyzed by DPClusO and BiClusO, respectively, and finally, the biomarker potential of miRNAs was computationally screened by miRNA-BD. In this study, a total of 55 significant gene-disease modules were identified, allowing us to propose 178 new HSCR candidate genes and two biological pathways. Moreover, we identified 12 key miRNAs with biomarker potential among 137 predicted HSCR-associated miRNAs. Functional analysis of new candidates showed that enrichment terms related to gene ontology (GO) and pathways were associated with HSCR. In conclusion, this approach has allowed us to decipher new clues of the etiopathogenesis of HSCR, although molecular experiments are further needed for clinical validations.

Network medicine: an approach to complex kidney disease phenotypes

Scientific reductionism has been the basis of disease classification and understanding for more than a century. However, the reductionist approach of characterizing diseases from a limited set of clinical observations and laboratory evaluations has proven insufficient in the face of an exponential growth in data generated from transcriptomics, proteomics, metabolomics and deep phenotyping. A new systematic method is necessary to organize these datasets and build new definitions of what constitutes a disease that incorporates both biological and environmental factors to more precisely describe the ever-growing complexity of phenotypes and their underlying molecular determinants. Network medicine provides such a conceptual framework to bridge these vast quantities of data while providing an individualized understanding of disease. The modern application of network medicine principles is yielding new insights into the pathobiology of chronic kidney diseases and renovascular disorders by expanding the understanding of pathogenic mediators, novel biomarkers and new options for renal therapeutics. These efforts affirm network medicine as a robust paradigm for elucidating new advances in the diagnosis and treatment of kidney disorders.

Tat-hspb1 Suppresses Clear Cell Renal Cell Carcinoma (ccRCC) Growth via Lysosomal Membrane Permeabilization

Clear cell renal cell carcinoma (ccRCC) is the most prevalent kidney cancer, of which the incidence is increasing worldwide with a high mortality rate. Bioactive peptides are considered a significant class of natural medicines. We applied mass spectrometry-based peptidomic analysis to explore the peptide profile of human renal clear cell carcinoma and adjacent normal tissues. A total of 18,031 peptides were identified, of which 105 unique peptides were differentially expressed (44 were up-regulated and 61 were down-regulated in ccRCC tissues). Through bioinformatic analysis, we finally selected one peptide derived from the HSPB1 protein (amino acids 12-35 of the N-terminal region of HSPB1). Next, we fused this peptide to the HIV-Tat, generated a novel peptide named Tat-hspb1, and found that Tat-hspb1 inhibited ccRCC cells' viability while being less cytotoxic to normal epithelial cells. Furthermore, Tat-hspb1 induced apoptosis and inhibited the proliferation and migration of ccRCC cells. Furthermore, we demonstrated that Tat-hspb1 was predominantly localized in lysosomes after entering the ccRCC cell and induced lysosomal membrane permeabilization (LMP) and the release of cathepsin D from lysosomes. Taken together, Tat-hspb1 has the potential to serve as a new anticancer drug candidate.

A novel candidate disease gene prioritization method using deep graph convolutional networks and semi-supervised learning

Background

Selecting and prioritizing candidate disease genes is necessary before conducting laboratory studies as identifying disease genes from a large number of candidate genes using laboratory methods, is a very costly and time-consuming task. There are many machine learning-based gene prioritization methods. These methods differ in various aspects including the feature vectors of genes, the used datasets with different structures, and the learning model. Creating a suitable feature vector for genes and an appropriate learning model on a variety of data with different and non-Euclidean structures, including graphs, as well as the lack of negative data are very important challenges of these methods. The use of graph neural networks has recently emerged in machine learning and other related fields, and they have demonstrated superior performance for a broad range of problems.

Methods

In this study, a new semi-supervised learning method based on graph convolutional networks is presented using the novel constructing feature vector for each gene. In the proposed method, first, we construct three feature vectors for each gene using terms from the Gene Ontology (GO) database. Then, we train a graph convolution network on these vectors using protein–protein interaction (PPI) network data to identify disease candidate genes. Our model discovers hidden layer representations encoding in both local graph structure as well as features of nodes. This method is characterized by the simultaneous consideration of topological information of the biological network (e.g., PPI) and other sources of evidence. Finally, a validation has been done to demonstrate the efficiency of our method.

Results

Several experiments are performed on 16 diseases to evaluate the proposed method's performance. The experiments demonstrate that our proposed method achieves the best results, in terms of precision, the area under the ROC curve (AUCs), and F1-score values, when compared with eight state-of-the-art network and machine learning-based disease gene prioritization methods.

Conclusion

This study shows that the proposed semi-supervised learning method appropriately classifies and ranks candidate disease genes using a graph convolutional network and an innovative method to create three feature vectors for genes based on the molecular function, cellular component, and biological process terms from GO data.

ImitateDB: A database for domain and motif mimicry incorporating host and pathogen protein interactions

Molecular mimicry of host proteins by pathogens constitutes a strategy to hijack the host pathways. At present, there is no dedicated resource for mimicked domains and motifs in the host-pathogen interactome. In this work, the experimental host-pathogen (HP) and host-host (HH) protein-protein interactions (PPIs) were collated. The domains and motifs of these proteins were annotated using CD Search and ScanProsite, respectively. Host and pathogen proteins with a shared host interactor and similar domain/motif constitute a mimicry pair exhibiting global structural similarity (domain mimicry pair; DMP) or local sequence motif similarity (motif mimicry pair; MMP). Mimicry pairs are likely to be co-expressed and co-localized. 1,97,607 DMPs and 32,67,568 MMPs were identified in 49,265 experimental HP-PPIs and organized in a web-based resource, ImitateDB ( http://imitatedb.sblab-nsit.net ) that can be easily queried. The results are externally integrated using hyperlinked domain PSSM ID, motif ID, protein ID and PubMed ID. Kinase, UL36, Smc and DEXDc were frequent DMP domains whereas protein kinase C phosphorylation, casein kinase 2 phosphorylation, glycosylation and myristoylation sites were frequent MMP motifs. Novel DMP domains SANT, Tudor, PhoX and MMP motif microbody C-terminal targeting signal, cornichon signature and lipocalin signature were proposed. ImitateDB is a novel resource for identifying mimicry in interacting host and pathogen proteins.

An insight into major signaling pathways and protein-protein interaction networks involved in the pathogenesis of gestational diabetes mellitus

Gestational diabetes mellitus (GDM) is associated with the increase of glucose in the blood rather than being absorbed by the cells. A better understanding of the signaling pathways is necessary to understand the pathophysiology of GDM. This study provides details about a series of signaling pathways and protein-protein interactions involved in the pathogenesis of GDM and their evaluations in GDM development. Protein-protein interactions were found between proteins of several signaling pathways that suggest interlink between these signaling pathways. Protein-protein interactions were generated with high confidence interaction scores based on textmining, co-occurrence, coexpression, neighborhood, gene fusion, experiments and databases. The dysregulation of signaling pathways may also contribute to the increased risk of complications associated with GDM in the mother and child. Further, studies on signaling pathways involved in the pathogenesis of GDM would help in the development of an effective intervention to prevent GDM along with the identification of key targets for effective therapies in the future. This article is protected by copyright. All rights reserved.

Network-Based Approach to Repurpose Approved Drugs for COVID-19 by Integrating GWAS and Text Mining Data

The coronavirus disease 19 (COVID-19) is a global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has a rapidly increasing prevalence and has caused significant morbidity/mortality. Despite the availability of many vaccines that can offer widespread immunization, it is also important to reach effective treatment for COVID-19 patients. However, the development of novel drug therapeutics is usually a time-consuming and costly process, and therefore, repositioning drugs that were previously approved for other purposes could have a major impact on the fight against COVID-19. Here, we first identified lung-specific gene regulatory/interaction subnetworks (COVID-19-related genes modules) enriched for COVID-19-associated genes obtained from GWAS and text mining. We then screened the targets of 220 approved drugs from DrugBank, obtained their drug-induced gene expression profiles in the LINCS database, and constructed lung-specific drug-related gene modules. By applying an integrated network-based approach to quantify the interactions of the COVID-19-related gene modules and drug-related gene modules, we prioritized 13 approved drugs (e.g., alitretinoin, clocortolone, terazosin, doconexent, and pergolide) that could potentially be repurposed for the treatment of COVID-19. These findings provide important and timely insights into alternative therapeutic options that should be further explored as COVID-19 continues to spread.

Low-density-lipoprotein-receptor-related protein 1 mediates Notch pathway activation

The Notch signaling pathway controls cell growth, differentiation, and fate decisions, and its dysregulation has been linked to various human genetic disorders and cancers. To comprehensively understand the global organization of the Notch pathway and identify potential drug targets for Notch-related diseases, we established a protein interaction landscape for the human Notch pathway. By combining and analyzing genetic and phenotypic data with bioinformatics analysis, we greatly expanded this pathway and identified many key regulators, including low-density-lipoprotein-receptor-related protein 1 (LRP1). We demonstrated that LRP1 mediates the ubiquitination chain linkage switching of Delta ligands, which further affects ligand recycling, membrane localization, and stability. LRP1 inhibition led to Notch signaling inhibition and decreased tumorigenesis in leukemia models. Our study provides a glimpse into the Notch pathway interaction network and uncovers LRP1 as one critical regulator of the Notch pathway, as well as a possible therapeutic target for Notch-related cancers.

Optimized Molecular Interaction Networks for the Study of Skeletal Muscle

BACKGROUND

Molecular interaction networks (MINs) aim to capture the complex relationships between interacting molecules within a biological system. MINs can be constructed from existing knowledge of molecular functional associations, such as protein-protein binding interactions (PPI) or gene co-expression, and these different sources may be combined into a single MIN. A given MIN may be more or less optimal in its representation of the important functional relationships of molecules in a tissue.

OBJECTIVE

The aim of this study was to establish whether a combined MIN derived from different types of functional association could better capture muscle-relevant biology compared to its constituent single-source MINs.

METHODS

MINs were constructed from functional association databases for both protein-binding and gene co-expression. The networks were then compared based on the capture of muscle-relevant genes and gene ontology (GO) terms, tested in two different ways using established biological network clustering algorithms. The top performing MINs were combined to test whether an optimal MIN for skeletal muscle could be constructed.

RESULTS

The STRING PPI network was the best performing single-source MIN among those tested. Combining STRING with interactions from either the MyoMiner or CoXPRESSdb gene co-expression sources resulted in a combined network with improved performance relative to its constituent networks.

CONCLUSION

MINs constructed from multiple types of functional association can better represent the functional relationships of molecules in a given tissue. Such networks may be used to improve the analysis and interpretation of functional genomics data in the study of skeletal muscle and neuromuscular diseases. Networks and clusters described by this study, including the combinations of STRING with MyoMiner or with CoXPRESSdb, are available for download from https://www.sys-myo.com/myominer/download.php.

Protein Interaction Network-based Deep Learning Framework for Identifying Disease-Associated Human Proteins

Infectious diseases in humans appear to be one of the most primary public health issues. Identification of novel disease-associated proteins will furnish an efficient recognition of the novel therapeutic targets. Here, we develop a Graph Convolutional Network (GCN)-based model called PINDeL to identify the disease-associated host proteins by integrating the human Protein Locality Graph and its corresponding topological features. Because of the amalgamation of GCN with the protein interaction network, PINDeL achieves the highest accuracy of 83.45% while AUROC and AUPRC values are 0.90 and 0.88, respectively. With high accuracy, recall, F1-score, specificity, AUROC, and AUPRC, PINDeL outperforms other existing machine-learning and deep-learning techniques for disease gene/protein identification in humans. Application of PINDeL on an independent dataset of 24320 proteins, which are not used for training, validation, or testing purposes, predicts 6448 new disease-protein associations of which we verify 3196 disease-proteins through experimental evidence like disease ontology, Gene Ontology, and KEGG pathway enrichment analyses. Our investigation informs that experimentally-verified 748 proteins are indeed responsible for pathogen-host protein interactions of which 22 disease-proteins share their association with multiple diseases such as cancer, aging, chem-dependency, pharmacogenomics, normal variation, infection, and immune-related diseases. This unique Graph Convolution Network-based prediction model is of utmost use in large-scale disease-protein association prediction and hence, will provide crucial insights on disease pathogenesis and will further aid in developing novel therapeutics.

Integration of peripheral transcriptomics, genomics, and interactomics following trauma identifies causal genes for symptoms of post-traumatic stress and major depression

Posttraumatic stress disorder (PTSD) is a debilitating syndrome with substantial morbidity and mortality that occurs in the aftermath of trauma. Symptoms of major depressive disorder (MDD) are also a frequent consequence of trauma exposure. Identifying novel risk markers in the immediate aftermath of trauma is a critical step for the identification of novel biological targets to understand mechanisms of pathophysiology and prevention, as well as the determination of patients most at risk who may benefit from immediate intervention. Our study utilizes a novel approach to computationally integrate blood-based transcriptomics, genomics, and interactomics to understand the development of risk vs. resilience in the months following trauma exposure. In a two-site longitudinal, observational prospective study, we assessed over 10,000 individuals and enrolled >700 subjects in the immediate aftermath of trauma (average 5.3 h post-trauma (range 0.5-12 h)) in the Grady Memorial Hospital (Atlanta) and Jackson Memorial Hospital (Miami) emergency departments. RNA expression data and 6-month follow-up data were available for 366 individuals, while genotype, transcriptome, and phenotype data were available for 297 patients. To maximize our power and understanding of genes and pathways that predict risk vs. resilience, we utilized a set-cover approach to capture fluctuations of gene expression of PTSD or depression-converting patients and non-converting trauma-exposed controls to find representative sets of disease-relevant dysregulated genes. We annotated such genes with their corresponding expression quantitative trait loci and applied a variant of a current flow algorithm to identify genes that potentially were causal for the observed dysregulation of disease genes involved in the development of depression and PTSD symptoms after trauma exposure. We obtained a final list of 11 driver causal genes related to MDD symptoms, 13 genes for PTSD symptoms, and 22 genes in PTSD and/or MDD. We observed that these individual or combined disorders shared ESR1, RUNX1, PPARA, and WWOX as driver causal genes, while other genes appeared to be causal driver in the PTSD only or MDD only cases. A number of these identified causal pathways have been previously implicated in the biology or genetics of PTSD and MDD, as well as in preclinical models of amygdala function and fear regulation. Our work provides a promising set of initial pathways that may underlie causal mechanisms in the development of PTSD or MDD in the aftermath of trauma.

What Can Machine Learning Approaches in Genomics Tell Us about the Molecular Basis of Amyotrophic Lateral Sclerosis?

Amyotrophic Lateral Sclerosis (ALS) is the most common late-onset motor neuron disorder, but our current knowledge of the molecular mechanisms and pathways underlying this disease remain elusive. This review (1) systematically identifies machine learning studies aimed at the understanding of the genetic architecture of ALS, (2) outlines the main challenges faced and compares the different approaches that have been used to confront them, and (3) compares the experimental designs and results produced by those approaches and describes their reproducibility in terms of biological results and the performances of the machine learning models. The majority of the collected studies incorporated prior knowledge of ALS into their feature selection approaches, and trained their machine learning models using genomic data combined with other types of mined knowledge including functional associations, protein-protein interactions, disease/tissue-specific information, epigenetic data, and known ALS phenotype-genotype associations. The importance of incorporating gene-gene interactions and cis-regulatory elements into the experimental design of future ALS machine learning studies is highlighted. Lastly, it is suggested that future advances in the genomic and machine learning fields will bring about a better understanding of ALS genetic architecture, and enable improved personalized approaches to this and other devastating and complex diseases.

Network medicine in Cardiovascular Research

The ability to generate multi-omics data coupled with deeply characterizing the clinical phenotype of individual patients promises to improve understanding of complex cardiovascular pathobiology. There remains an important disconnection between the magnitude and granularity of these data and our ability to improve phenotype-genotype correlations for complex cardiovascular diseases. This shortcoming may be due to limitations associated with traditional reductionist analytical methods, which tend to emphasize a single molecular event in the pathogenesis of diseases more aptly characterized by crosstalk between overlapping molecular pathways. Network medicine is a rapidly growing discipline that considers diseases as the consequences of perturbed interactions between multiple interconnected biological components. This powerful integrative approach has enabled a number of important discoveries in complex disease mechanisms. In this review, we introduce the basic concepts of network medicine and highlight specific examples by which this approach has accelerated cardiovascular research. We also review how network medicine is well-positioned to promote rational drug design for patients with cardiovascular diseases, with particular emphasis on advancing precision medicine.

Comparison of gene regulatory networks to identify pathogenic genes for lymphoma

Lymphoma is the most complicated cancer that can be divided into several tens of subtypes. It may occur in any part of body that has lymphocytes, and is closely correlated with diverse environmental factors such as the ionizing radiation, chemocarcinogenesis, and virus infection. All the environmental factors affect the lymphoma through genes. Identifying pathogenic genes for lymphoma is consequently an essential task to understand its complexity in a unified framework. In this paper, we propose a new method to expose high-confident edges in gene regulatory networks (GRNs) for a total of 32 organs, called Filtered GRNs (f-GRNs), comparison of which gives us a proper reference for the Lymphoma, i.e. the B-lymphocytes cells, whose f-GRN is closest with that for the Lymphoma. By using the Gene Ontology and Biological Process analysis we display the differences of the two networks' hubs in biological functions. Matching with the Genecards shows that most of the hubs take part in the genetic information transmission and expression, except a specific gene of Retinoic Acid Receptor Alpha (RARA) that encodes the retinoic acid receptor. In the lymphoma, the genes in the RARA ego-network are involved in two cancer pathways, and the RARA is present only in these cancer pathways. For the lymphoid B cells, however, the genes in the RARA ego-network do not participate in cancer-related pathways.

Cryo-EM Reveals Unanchored M1-Ubiquitin Chain Binding at hRpn11 of the 26S Proteasome

The 26S proteasome is specialized for regulated protein degradation and formed by a dynamic regulatory particle (RP) that caps a hollow cylindrical core particle (CP) where substrates are proteolyzed. Its diverse substrates unify as proteasome targets by ubiquitination. We used cryogenic electron microscopy (cryo-EM) to study how human 26S proteasome interacts with M1-linked hexaubiquitin (M1-Ub6) unanchored to a substrate and E3 ubiquitin ligase E6AP/UBE3A. Proteasome structures are available with model substrates extending through the RP ATPase ring and substrate-conjugated K63-linked ubiquitin chains present at inhibited deubiquitinating enzyme hRpn11 and the nearby ATPase hRpt4/hRpt5 coiled coil. In this study, we find M1-Ub6 at the hRpn11 site despite the absence of conjugated substrate, indicating that ubiquitin binding at this location does not require substrate interaction with the RP. Moreover, unanchored M1-Ub6 binds to this hRpn11 site of the proteasome with the CP gating residues in both the closed and opened conformational states.

Analysis of Protein-Protein Interaction Networks through Computational Approaches

The interactions among proteins and genes are extremely important for cellular functions. Molecular interactions at protein or gene levels can be used to construct interaction networks in which the interacting species are categorized based on direct interactions or functional similarities. Compared with the limited experimental techniques, various computational tools make it possible to analyze, filter, and combine the interaction data to get comprehensive information about the biological pathways. By the efficient way of integrating experimental findings in discovering PPIs and computational techniques for prediction, the researchers have been able to gain many valuable data on PPIs, including some advanced databases. Moreover, many useful tools and visualization programs enable the researchers to establish, annotate, and analyze biological networks. We here review and list the computational methods, databases, and tools for protein-protein interaction prediction.

A Machine Learning Approach to Predicting Autism Risk Genes: Validation of Known Genes and Discovery of New Candidates

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic basis. The role of de novo mutations in ASD has been well established, but the set of genes implicated to date is still far from complete. The current study employs a machine learning-based approach to predict ASD risk genes using features from spatiotemporal gene expression patterns in human brain, gene-level constraint metrics, and other gene variation features. The genes identified through our prediction model were enriched for independent sets of ASD risk genes, and tended to be down-expressed in ASD brains, especially in frontal and parietal cortex. The highest-ranked genes not only included those with strong prior evidence for involvement in ASD (for example, NBEA, HERC1, and TCF20), but also indicated potentially novel candidates, such as, MYCBP2 and CAND1, which are involved in protein ubiquitination. We also showed that our method outperformed state-of-the-art scoring systems for ranking curated ASD candidate genes. Gene ontology enrichment analysis of our predicted risk genes revealed biological processes clearly relevant to ASD, including neuronal signaling, neurogenesis, and chromatin remodeling, but also highlighted other potential mechanisms that might underlie ASD, such as regulation of RNA alternative splicing and ubiquitination pathway related to protein degradation. Our study demonstrates that human brain spatiotemporal gene expression patterns and gene-level constraint metrics can help predict ASD risk genes. Our gene ranking system provides a useful resource for prioritizing ASD candidate genes.

An integrated method for identifying essential proteins from multiplex network model of protein-protein interactions

Cell survival requires the presence of essential proteins. Detection of essential proteins is relevant not only because of the critical biological functions they perform but also the role played by them as a drug target against pathogens. Several computational techniques are in place to identify essential proteins based on protein-protein interaction (PPI) network. Essential protein detection using only physical interaction data of proteins is challenging due to its inherent uncertainty. Hence, in this work, we propose a multiplex network-based framework that incorporates multiple protein interaction data from their physical, coexpression and phylogenetic profiles. An extended version termed as multiplex eigenvector centrality (MEC) is used to identify essential proteins from this network. The methodology integrates the score obtained from the multiplex analysis with subcellular localization and Gene Ontology information and is implemented using Saccharomyces cerevisiae datasets. The proposed method outperformed many recent essential protein prediction techniques in the literature.

uKIN Combines New and Prior Information with Guided Network Propagation to Accurately Identify Disease Genes

Protein interaction networks provide a powerful framework for identifying genes causal for complex genetic diseases. Here, we introduce a general framework, uKIN, that uses prior knowledge of disease-associated genes to guide, within known protein-protein interaction networks, random walks that are initiated from newly identified candidate genes. In large-scale testing across 24 cancer types, we demonstrate that our network propagation approach for integrating both prior and new information not only better identifies cancer driver genes than using either source of information alone but also readily outperforms other state-of-the-art network-based approaches. We also apply our approach to genome-wide association data to identify genes functionally relevant for several complex diseases. Overall, our work suggests that guided network propagation approaches that utilize both prior and new data are a powerful means to identify disease genes. uKIN is freely available for download at: https://github.com/Singh-Lab/uKIN.

A Knowledge-Based Machine Learning Approach to Gene Prioritisation in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis is a neurodegenerative disease of the upper and lower motor neurons resulting in death from neuromuscular respiratory failure, typically within two to five years of first symptoms. Several rare disruptive gene variants have been associated with ALS and are responsible for about 15% of all cases. Although our knowledge of the genetic landscape of this disease is improving, it remains limited. Machine learning models trained on the available protein-protein interaction and phenotype-genotype association data can use our current knowledge of the disease genetics for the prediction of novel candidate genes. Here, we describe a knowledge-based machine learning method for this purpose. We trained our model on protein-protein interaction data from IntAct, gene function annotation from Gene Ontology, and known disease-gene associations from DisGeNet. Using several sets of known ALS genes from public databases and a manual review as input, we generated a list of new candidate genes for each input set. We investigated the relevance of the predicted genes in ALS by using the available summary statistics from the largest ALS genome-wide association study and by performing functional and phenotype enrichment analysis. The predicted sets were enriched for genes associated with other neurodegenerative diseases known to overlap with ALS genetically and phenotypically, as well as for biological processes associated with the disease. Moreover, using ALS genes from ClinVar and our manual review as input, the predicted sets were enriched for ALS-associated genes (ClinVar p = 0.038 and manual review p = 0.060) when used for gene prioritisation in a genome-wide association study.

HPOLabeler: improving prediction of human protein–phenotype associations by learning to rank

Motivation

Annotating human proteins by abnormal phenotypes has become an important topic. Human Phenotype Ontology (HPO) is a standardized vocabulary of phenotypic abnormalities encountered in human diseases. As of November 2019, only <4000 proteins have been annotated with HPO. Thus, a computational approach for accurately predicting protein–HPO associations would be important, whereas no methods have outperformed a simple Naive approach in the second Critical Assessment of Functional Annotation, 2013–2014 (CAFA2).

Results

We present HPOLabeler, which is able to use a wide variety of evidence, such as protein–protein interaction (PPI) networks, Gene Ontology, InterPro, trigram frequency and HPO term frequency, in the framework of learning to rank (LTR). LTR has been proved to be powerful for solving large-scale, multi-label ranking problems in bioinformatics. Given an input protein, LTR outputs the ranked list of HPO terms from a series of input scores given to the candidate HPO terms by component learning models (logistic regression, nearest neighbor and a Naive method), which are trained from given multiple evidence. We empirically evaluate HPOLabeler extensively through mainly two experiments of cross validation and temporal validation, for which HPOLabeler significantly outperformed all component models and competing methods including the current state-of-the-art method. We further found that (i) PPI is most informative for prediction among diverse data sources and (ii) low prediction performance of temporal validation might be caused by incomplete annotation of new proteins.

Availability and implementation

http://issubmission.sjtu.edu.cn/hpolabeler/.

Contact

zhusf@fudan.edu.cn

Supplementary information

Supplementary data are available at Bioinformatics online.

Construction, comparison and evolution of networks in life sciences and other disciplines

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

Molecular networks in Network Medicine: Development and applications

Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.

Predicting subcellular localization of proteins using protein-protein interaction data

The knowledge of subcellular localization of proteins can provide useful clues about their functions. The conventional methods to determine the subcellular localization are unable to keep pace with the rate at which the new data is being generated. Thus, though sequence information is available, the localization and function of a number of proteins remains unknown. In this study, we have developed a script that makes use of the physical interactors of a protein and their localization data to predict the subcellular localization. We used the script to predict the localization of yeast proteins for which there is no localization data. Further, we experimentally verified the predicted localization for six arbitrarily chosen proteins and found our predictions to be correct for five of the proteins.

Visual Analysis of Protein-Protein Interaction Docking Models Using COZOID Tool

Networks of protein-protein interactions (PPI) constitute either stable or transient complexes in every cell. Most of the cellular complexes keep their function, and therefore stay similar, during evolution. The evolutionary constraints preserve most cellular functions via preservation of protein structures and interactions. The evolutionary conservation information is utilized in template-based approaches, like protein structure modeling or docking. Here we use the combination of the template-free docking method with conservation-based selection of the best docking model using our newly developed COZOID tool.We describe a step-by-step protocol for visual selection of docking models, based on their similarity to the original protein complex structure. Using the COZOID tool, we first analyze contact zones of the original complex structure and select contact amino acids for docking restraints. Then we model and dock the homologous proteins. Finally, we utilize different analytical modes of our COZOID tool to select the docking models most similar to the original complex structure.

Understanding allergic multimorbidity within the non-eosinophilic interactome

BACKGROUND

The mechanisms explaining multimorbidity between asthma, dermatitis and rhinitis (allergic multimorbidity) are not well known. We investigated these mechanisms and their specificity in distinct cell types by means of an interactome-based analysis of expression data.

METHODS

Genes associated to the diseases were identified using data mining approaches, and their multimorbidity mechanisms in distinct cell types were characterized by means of an in silico analysis of the topology of the human interactome.

RESULTS

We characterized specific pathomechanisms for multimorbidities between asthma, dermatitis and rhinitis for distinct emergent non-eosinophilic cell types. We observed differential roles for cytokine signaling, TLR-mediated signaling and metabolic pathways for multimorbidities across distinct cell types. Furthermore, we also identified individual genes potentially associated to multimorbidity mechanisms.

CONCLUSIONS

Our results support the existence of differentiated multimorbidity mechanisms between asthma, dermatitis and rhinitis at cell type level, as well as mechanisms common to distinct cell types. These results will help understanding the biology underlying allergic multimorbidity, assisting in the design of new clinical studies.

Integration of Heterogeneous Experimental Data Improves Global Map of Human Protein Complexes

Protein complexes play a significant role in the core functionality of cells. These complexes are typically identified by detecting densely connected subgraphs in protein-protein interaction (PPI) networks. Recently, multiple large-scale mass spectrometry-based experiments have significantly increased the availability of PPI data in order to further expand the set of known complexes. However, high-throughput experimental data generally are incomplete, show limited agreement between experiments, and show frequent false positive interactions. There is a need for computational approaches that can address these limitations in order to improve the coverage and accuracy of human protein complexes. Here, we present a new method that integrates data from multiple heterogeneous experiments and sources in order to increase the reliability and coverage of predicted protein complexes. We first fused the heterogeneous data into a feature matrix and trained classifiers to score pairwise protein interactions. We next used graph based methods to combine pairwise interactions into predicted protein complexes. Our approach improves the accuracy and coverage of protein pairwise interactions, accurately identifies known complexes, and suggests both novel additions to known complexes and entirely new complexes. Our results suggest that integration of heterogeneous experimental data helps improve the reliability and coverage of diverse high-throughput mass-spectrometry experiments, leading to an improved global map of human protein complexes.

A dual controllability analysis of influenza virus-host protein-protein interaction networks for antiviral drug target discovery

BACKGROUND

Host factors of influenza virus replication are often found in key topological positions within protein-protein interaction networks. This work explores how protein states can be manipulated through controllability analysis: the determination of the minimum manipulation needed to drive the cell system to any desired state. Here, we complete a two-part controllability analysis of two protein networks: a host network representing the healthy cell state and an influenza A virus-host network representing the infected cell state. In this context, controllability analyses aim to identify key regulating host factors of the infected cell's progression. This knowledge can be utilized in further biological analysis to understand disease dynamics and isolate proteins for study as drug target candidates.

RESULTS

Both topological and controllability analyses provide evidence of wide-reaching network effects stemming from the addition of viral-host protein interactions. Virus interacting and driver host proteins are significant both topologically and in controllability, therefore playing important roles in cell behavior during infection. Functional analysis finds overlap of results with previous siRNA studies of host factors involved in influenza replication, NF-kB pathway and infection relevance, and roles as interferon regulating genes. 24 proteins are identified as holding regulatory roles specific to the infected cell by measures of topology, controllability, and functional role. These proteins are recommended for further study as potential antiviral drug targets.

CONCLUSIONS

Seasonal outbreaks of influenza A virus are a major cause of illness and death around the world each year with a constant threat of pandemic infection. This research aims to increase the efficiency of antiviral drug target discovery using existing protein-protein interaction data and network analysis methods. These results are beneficial to future studies of influenza virus, both experimental and computational, and provide evidence that the combination of topology and controllability analyses may be valuable for future efforts in drug target discovery.

Substrate complex competition is a regulatory motif that allows NFκB RelA to license but not amplify NFκB RelB

Signaling pathways often share molecular components, tying the activity of one pathway to the functioning of another. In the NFκB signaling system, distinct kinases mediate inflammatory and developmental signaling via RelA and RelB, respectively. Although the substrates of the developmental, so-called noncanonical, pathway are induced by inflammatory/canonical signaling, crosstalk is limited. Through dynamical systems modeling, we identified the underlying regulatory mechanism. We found that as the substrate of the noncanonical kinase NIK, the nfkb2 gene product p100, transitions from a monomer to a multimeric complex, it may compete with and inhibit p100 processing to the active p52. Although multimeric complexes of p100 (IκBδ) are known to inhibit preexisting RelA:p50 through sequestration, here we report that p100 complexes can inhibit the enzymatic formation of RelB:p52. We show that the dose-response systems properties of this complex substrate competition motif are poorly accounted for by standard Michaelis-Menten kinetics, but require more detailed mass action formulations. In sum, although tonic inflammatory signaling is required for adequate expression of the noncanonical pathway precursors, the substrate complex competition motif identified here can prevent amplification of the active RelB:p52 dimer in elevated inflammatory conditions to ensure reliable RelB-dependent developmental signaling independent of inflammatory context.

Identification of epistasis loci underlying rice flowering time by controlling population stratification and polygenic effect

Flowering time is an important agronomic trait, attributed by multiple genes, gene-gene interactions and environmental factors. Population stratification and polygenic effects might confound genetic effects of the causal loci underlying this complex trait. We proposed a two-step approach for detecting epistasis interactions underlying rice flowering time by accounting population structure and polygenic effects. Simulation studies showed that the approach used in this study performs better than classical and PC-linear approaches in terms of powers and false discovery rates in the case of population stratification and polygenic effects. Whole genome epistasis analyses identified 589 putative genetic interactions for flowering time. Eighteen of these interactions are located within 10 kilobases of regions of known protein-protein interactions. Thirty-seven SNPs near to twenty-five genes involve in rice or/and Arabidopsis (orthologue) flowering pathway. Bioinformatics analysis showed that 66.55% pairwise genes of the identified interactions (392 out of the 589 interactions) have similarity in various genomic features. Moreover, significant numbers of detected epistatic genes have high expression in different floral tissues. Our findings highlight the importance of epistasis analysis by controlling population stratification and polygenic effect and provided novel insights into the genetic architecture of rice flowering which could assist breeding programmes.

Network-based prediction of protein interactions

Despite exceptional experimental efforts to map out the human interactome, the continued data incompleteness limits our ability to understand the molecular roots of human disease. Computational tools offer a promising alternative, helping identify biologically significant, yet unmapped protein-protein interactions (PPIs). While link prediction methods connect proteins on the basis of biological or network-based similarity, interacting proteins are not necessarily similar and similar proteins do not necessarily interact. Here, we offer structural and evolutionary evidence that proteins interact not if they are similar to each other, but if one of them is similar to the other's partners. This approach, that mathematically relies on network paths of length three (L3), significantly outperforms all existing link prediction methods. Given its high accuracy, we show that L3 can offer mechanistic insights into disease mechanisms and can complement future experimental efforts to complete the human interactome.

Utilizing microphysiological systems and induced pluripotent stem cells for disease modeling: a case study for blood brain barrier research in a pharmaceutical setting

Microphysiological systems (MPS) may be able to provide the pharmaceutical industry models that can reflect human physiological responses to improve drug discovery and translational outcomes. With lack of efficacy being the primary cause for drug attrition, developing MPS disease models would help researchers identify novel targets, study mechanisms in more physiologically-relevant depth, screen for novel biomarkers and test/optimize various therapeutics (small molecules, nanoparticles and biologics). Furthermore, with advances in inducible pluripotent stem cell technology (iPSC), pharmaceutical companies can access cells from patients to help recreate specific disease phenotypes in MPS platforms. Combining iPSC and MPS technologies will contribute to our understanding of the complexities of neurodegenerative diseases and of the blood brain barrier (BBB) leading to development of enhanced therapeutics.

Allosteric Modulation of Intrinsically Disordered Proteins

The allosteric property of globular proteins is applauded as their intrinsic ability to regulate distant sites, and this property further plays a critical role in a wide variety of cellular regulatory mechanisms. Recent advancements and studies have revealed the manifestation of allostery in intrinsically disordered proteins or regions as allosteric sites present within or mediated by IDP/IDRs facilitates the signaling interactions for various biological mechanisms which would otherwise be impossible for globular proteins to regulate. This thematic review has highlighted the biological outcomes that can be achieved by the mechanism of allosteric regulation of intrinsically disordered proteins or regions. The similar mechanism has been implemented on Adenovirus 5 early region 1A and tumor apoptosis protein p53 in correspondence with other partners in binary and ternary complexes, which are the subject of the current review. Both these proteins regulate once they bind to their partners, consequently, forming either a binary or a ternary complex. Allosteric regulation by IDPs is currently a subject undergoing intense study, and the ongoing research work will ensure a better understanding of precision and efficiency of cellular regulation by them. Allosteric regulation mechanism can also be researched by intrinsically disordered protein-specific force field.

Protein-protein interaction network analysis of insecticide resistance molecular mechanism in Drosophila melanogaster

The problem of resistance has not been solved fundamentally at present, because the development speed of new insecticides can not keep pace with the development speed of resistance, and the lack of understanding of molecular mechanism of resistance. Here we collected seed genes and their interacting proteins involved in insecticide resistance molecular mechanism in Drosophila melanogaster by literature mining and the String database. We identified a total of 528 proteins and 13514 protein-protein interactions. The protein interaction network was constructed by String and Pajek, and we analyzed the topological properties, such as degree centrality and eigenvector centrality. Proteasome complexes and drug metabolism-cytochrome P450 were an enrichment by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. This is the first time to explore the insecticide resistance molecular mechanism of D. melanogaster by the methods and tools of network biology, it can provide the bioinformatic foundation for further understanding the mechanisms of insecticide resistance.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

Effect of IAPP on the proteome of cultured Rin-5F cells

Background

Islet amyloid polypeptide (IAPP) or amylin deposits can be found in the islets of type 2 diabetes patients. The peptide is suggested to be involved in the etiology of the disease through formation of amyloid deposits and destruction of β islet cells, though the underlying molecular events leading from IAPP deposition to β cell death are still largely unknown.

Results

We used OFFGEL™ proteomics to study how IAPP exposure affects the proteome of rat pancreatic insulinoma Rin-5F cells. The OFFGEL™ methodology is highly effective at generating quantitative data on hundreds of proteins affected by IAPP, with its accuracy confirmed by In Cell Western and Quantitative Real Time PCR results. Combining data on individual proteins identifies pathways and protein complexes affected by IAPP. IAPP disrupts protein synthesis and degradation, and induces oxidative stress. It causes decreases in protein transport and localization. IAPP disrupts the regulation of ubiquitin-dependent protein degradation and increases catabolic processes. IAPP causes decreases in protein transport and localization, and affects the cytoskeleton, DNA repair and oxidative stress.

Conclusions

Results are consistent with a model where IAPP aggregates overwhelm the ability of a cell to degrade proteins via the ubiquitin system. Ultimately this leads to apoptosis. IAPP aggregates may be also toxic to the cell by causing oxidative stress, leading to DNA damage or by decreasing protein transport. The reversal of any of these effects, perhaps by targeting proteins which alter in response to IAPP, may be beneficial for type II diabetes.

Electronic supplementary material

The online version of this article (10.1186/s12858-018-0099-3) contains supplementary material, which is available to authorized users.

OligoPVP: Phenotype-driven analysis of individual genomic information to prioritize oligogenic disease variants

An increasing number of disorders have been identified for which two or more distinct alleles in two or more genes are required to either cause the disease or to significantly modify its onset, severity or phenotype. It is difficult to discover such interactions using existing approaches. The purpose of our work is to develop and evaluate a system that can identify combinations of alleles underlying digenic and oligogenic diseases in individual whole exome or whole genome sequences. Information that links patient phenotypes to databases of gene-phenotype associations observed in clinical or non-human model organism research can provide useful information and improve variant prioritization for genetic diseases. Additional background knowledge about interactions between genes can be utilized to identify sets of variants in different genes in the same individual which may then contribute to the overall disease phenotype. We have developed OligoPVP, an algorithm that can be used to prioritize causative combinations of variants in digenic and oligogenic diseases, using whole exome or whole genome sequences together with patient phenotypes as input. We demonstrate that OligoPVP has significantly improved performance when compared to state of the art pathogenicity detection methods in the case of digenic diseases. Our results show that OligoPVP can efficiently prioritize sets of variants in digenic diseases using a phenotype-driven approach and identify etiologically important variants in whole genomes. OligoPVP naturally extends to oligogenic disease involving interactions between variants in two or more genes. It can be applied to the identification of multiple interacting candidate variants contributing to phenotype, where the action of modifier genes is suspected from pedigree analysis or failure of traditional causative variant identification.

Covalent Rpn13-Binding Inhibitors for the Treatment of Ovarian Cancer

Substitution of the m,p-chloro groups of bis-benzylidinepiperidone RA190 for p-nitro, generating RA183, enhanced covalent drug binding to Cys88 of RPN13. Treatment of cancer cell lines with RA183 inhibited ubiquitin-mediated protein degradation, resulting in rapid accumulation of high-molecular-weight polyubiquitinated proteins, blockade of NFκB signaling, endoplasmic reticulum stress, an unfolded protein response, production of reactive oxygen species, and apoptotic cell death. High-grade ovarian cancer, triple-negative breast cancer, and multiple myeloma cell lines were particularly vulnerable to RA183. RA183 stabilized a tetraubiquitin-linked firefly luciferase reporter protein in cancer cell lines and mice, demonstrating in vitro and in vivo proteasomal inhibition, respectively. However, RA183 was rapidly cleared from plasma, likely reflecting its rapid degradation to the active compound RA9, as seen in human liver microsomes. Intraperitoneal administration of RA183 inhibited proteasome function and orthotopic tumor growth in mice bearing human ovarian cancer model ES2-luc ascites or syngeneic ID8-luc tumor.

PMLPR: A novel method for predicting subcellular localization based on recommender systems

The importance of protein subcellular localization problem is due to the importance of protein's functions in different cell parts. Moreover, prediction of subcellular locations helps to identify the potential molecular targets for drugs and has an important role in genome annotation. Most of the existing prediction methods assign only one location for each protein. But, since some proteins move between different subcellular locations, they can have multiple locations. In recent years, some multiple location predictors have been introduced. However, their performances are not accurate enough and there is much room for improvement. In this paper, we introduced a method, PMLPR, to predict locations for a protein. PMLPR predicts a list of locations for each protein based on recommender systems and it can properly overcome the multiple location prediction problem. For evaluating the performance of PMLPR, we considered six datasets RAT, FLY, HUMAN, Du et al., DBMLoc and Höglund. The performance of this algorithm is compared with six state-of-the-art algorithms, YLoc, WOLF-PSORT, prediction channel, MDLoc, Du et al. and MultiLoc2-HighRes. The results indicate that our proposed method is significantly superior on RAT and Fly proteins, and decent on HUMAN proteins. Moreover, on the datasets introduced by Du et al., DBMLoc and Höglund, PMLPR has comparable results. For the case study, we applied the algorithms on 8 proteins which are important in cancer research. The results of comparison with other methods indicate the efficiency of PMLPR.

A Comprehensive Analysis of Nuclear-Encoded Mitochondrial Genes in Schizophrenia

BACKGROUND

The genetic risk factors of schizophrenia (SCZ), a severe psychiatric disorder, are not yet fully understood. Multiple lines of evidence suggest that mitochondrial dysfunction may play a role in SCZ, but comprehensive association studies are lacking. We hypothesized that variants in nuclear-encoded mitochondrial genes influence susceptibility to SCZ.

METHODS

We conducted gene-based and gene-set analyses using summary association results from the Psychiatric Genomics Consortium Schizophrenia Phase 2 (PGC-SCZ2) genome-wide association study comprising 35,476 cases and 46,839 control subjects. We applied the MAGMA method to three sets of nuclear-encoded mitochondrial genes: oxidative phosphorylation genes, other nuclear-encoded mitochondrial genes, and genes involved in nucleus-mitochondria crosstalk. Furthermore, we conducted a replication study using the iPSYCH SCZ sample of 2290 cases and 21,621 control subjects.

RESULTS

In the PGC-SCZ2 sample, 1186 mitochondrial genes were analyzed, among which 159 had p values < .05 and 19 remained significant after multiple testing correction. A meta-analysis of 818 genes combining the PGC-SCZ2 and iPSYCH samples resulted in 104 nominally significant and nine significant genes, suggesting a polygenic model for the nuclear-encoded mitochondrial genes. Gene-set analysis, however, did not show significant results. In an in silico protein-protein interaction network analysis, 14 mitochondrial genes interacted directly with 158 SCZ risk genes identified in PGC-SCZ2 (permutation p = .02), and aldosterone signaling in epithelial cells and mitochondrial dysfunction pathways appeared to be overrepresented in this network of mitochondrial and SCZ risk genes.

CONCLUSIONS

This study provides evidence that specific aspects of mitochondrial function may play a role in SCZ, but we did not observe its broad involvement even using a large sample.

Identification of an interaction between calcium-dependent protein kinase 4 (EtCDPK4) and serine protease inhibitor (EtSerpin) in Eimeria tenella

Background

Eimeria tenella is an obligate intracellular apicomplexan protozoan parasite that has a complex life-cycle. Calcium ions, through various calcium-dependent protein kinases (CDPKs), regulate key events in parasite growth and development, including protein secretion, movement, differentiation, and invasion of and escape from host cells. In this study, we identified proteins that interact with EtCDPK4 to lay a foundation for clarifying the role of CDPKs in calcium channels.

Methods

Eimeria tenella merozoites were collected to construct a yeast two-hybrid (Y2H) cDNA library. The Y2H system was used to identify proteins that interact with EtCDPK4. One of interacting proteins was confirmed using bimolecular fluorescence complementation and co-immunoprecipitation in vivo. Co-localization of proteins was performed using immunofluorescence assays.

Results

Eight proteins that interact with EtCDPK4 were identified using the Y2H system. One of the proteins, E. tenella serine protease inhibitor 1 (EtSerpin), was further confirmed.

Conclusion

In this study, we screened for proteins that interact with EtCDPK4. An interaction between EtSerpin and EtCDPK4 was identified that may contribute to the invasion and development of E. tenella in host cells.

Prioritizing disease genes with an improved dual label propagation framework

BACKGROUND

Prioritizing disease genes is trying to identify potential disease causing genes for a given phenotype, which can be applied to reveal the inherited basis of human diseases and facilitate drug development. Our motivation is inspired by label propagation algorithm and the false positive protein-protein interactions that exist in the dataset. To the best of our knowledge, the false positive protein-protein interactions have not been considered before in disease gene prioritization. Label propagation has been successfully applied to prioritize disease causing genes in previous network-based methods. These network-based methods use basic label propagation, i.e. random walk, on networks to prioritize disease genes in different ways. However, all these methods can not deal with the situation in which plenty false positive protein-protein interactions exist in the dataset, because the PPI network is used as a fixed input in previous methods. This important characteristic of data source may cause a large deviation in results.

RESULTS

A novel network-based framework IDLP is proposed to prioritize candidate disease genes. IDLP effectively propagates labels throughout the PPI network and the phenotype similarity network. It avoids the method falling when few disease genes are known. Meanwhile, IDLP models the bias caused by false positive protein interactions and other potential factors by treating the PPI network matrix and the phenotype similarity matrix as the matrices to be learnt. By amending the noises in training matrices, it improves the performance results significantly. We conduct extensive experiments over OMIM datasets, and IDLP has demonstrated its effectiveness compared with eight state-of-the-art approaches. The robustness of IDLP is also validated by doing experiments with disturbed PPI network. Furthermore, We search the literatures to verify the predicted new genes got by IDLP are associated with the given diseases, the high prediction accuracy shows IDLP can be a powerful tool to help biologists discover new disease genes.

CONCLUSIONS

IDLP model is an effective method for disease gene prioritization, particularly for querying phenotypes without known associated genes, which would be greatly helpful for identifying disease genes for less studied phenotypes.

AVAILABILITY

https://github.com/nkiip/IDLP.

Predicted networks of protein-protein interactions in Stegodyphus mimosarum by cross-species comparisons

BACKGROUND

Stegodyphus mimosarum is a candidate model organism belonging to the class Arachnida in the phylum Arthropoda. Studies on the biology of S. mimosarum over the past several decades have consisted of behavioral research and comparison of gene sequences based on the assembled genome sequence. Given the lack of systematic protein analyses and the rich source of information in the genome, we predicted the relationships of proteins in S. mimosarum by bioinformatics comparison with genome-wide proteins from select model organisms using gene mapping.

RESULTS

The protein-protein interactions (PPIs) of 11 organisms were integrated from four databases (BioGrid, InAct, MINT, and DIP). Here, we present comprehensive prediction and analysis of 3810 proteins in S. mimosarum with regard to interactions between proteins using PPI data of organisms. Interestingly, a portion of the protein interactions conserved among Saccharomyces cerevisiae, Homo sapiens, Arabidopsis thaliana, and Drosophila melanogaster were found to be associated with RNA splicing. In addition, overlap of predicted PPIs in reference organisms, Gene Ontology, and topology models in S. mimosarum are also reported.

CONCLUSIONS

Addition of Stegodyphus, a spider representative of interactomic research, provides the possibility of obtaining deeper insights into the evolution of PPI networks among different animal species. This work largely supports the utility of the "stratus clouds" model for predicted PPIs, providing a roadmap for integrative systems biology in S. mimosarum.