Tissue-specific Protein-protein Interaction Research Articles

The PENGUIN approach to reconstruct protein interactions at enhancer-promoter regions and its application to prostate cancer

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/.

Immune cell-specific and common molecular signatures in rheumatoid arthritis through molecular network approaches

Rheumatoid arthritis (RA) is an autoimmune disorder and common symptom of RA is chronic synovial inflammation. The pathogenesis of RA is not fully understood. Therefore, we aimed to identify underlying common and distinct molecular signatures and pathways among ten types of tissue and cells obtained from patients with RA. In this study, transcriptomic data including synovial tissues, macrophages, blood, T cells, CD4+T cells, CD8+T cells, natural killer T (NKT), cells natural killer (NK) cells, neutrophils, and monocyte cells were analyzed with an integrative and comparative network biology perspective. Each dataset yielded a list of differentially expressed genes as well as a reconstruction of the tissue-specific protein-protein interaction (PPI) network. Molecular signatures were identified by a statistical test using the hypergeometric probability density function by employing the interactions of transcriptional regulators and PPI. Reporter metabolites of each dataset were determined by using genome-scale metabolic networks. It was defined as the common hub proteins, novel molecular signatures, and metabolites in two or more tissue types while immune cell-specific molecular signatures were identified, too. Importantly, miR-155-5p is found as a common miRNA in all tissues. Moreover, NCOA3, PRKDC and miR-3160 might be novel molecular signatures for RA. Our results establish a novel approach for identifying immune cell-specific molecular signatures of RA and provide insights into the role of common tissue-specific genes, miRNAs, TFs, receptors, and reporter metabolites. Experimental research should be used to validate the corresponding genes, miRNAs, and metabolites.

Affinity purification of in vivo assembled whirlin-associated protein complexes from the zebrafish retina

Mutations in WHRN lead to Usher syndrome type 2d or to non-syndromic hearing impairment. The WHRN-encoded gene product whirlin directly interacts with the intracellular regions of the other two Usher syndrome type 2-associated proteins, usherin and ADGRV1. In photoreceptor cells, this protein complex constitutes fibrous links between the periciliary membrane and the connecting cilium. However, the molecular mechanism(s) of retinal degeneration due to compromised formation and function of the USH2-associated protein complex remains elusive. To unravel this pathogenic mechanism, we isolated and characterized whirlin-associated protein complexes from zebrafish photoreceptor cells. We generated transgenic zebrafish that express Strep/FLAG-tagged Whrna, a zebrafish ortholog of human whirlin, under the control of a photoreceptor-specific promoter. Affinity purification of Strep/FLAG-tagged Whrna and associated proteins from adult transgenic zebrafish retinas followed by mass spectrometry identified 19 novel candidate associated proteins. Pull down experiments and dedicated yeast two-hybrid assays confirmed the association of Whrna with 7 of the co-purified proteins. Several of the co-purified proteins are part of the synaptic proteome, which indicates a role for whirlin in the photoreceptor synapse. Future studies will elucidate which of the newly identified protein-protein interactions contribute to the development of the retinal phenotype observed in USH2d patients. SignificanceSince protein-protein interactions identified using targeted in vitro studies do not always recapitulate interactions that are functionally relevant in vivo, we established a transgenic zebrafish line that stably expresses a Strep/FLAG-tagged ortholog of human whirlin (SF-Whrna) in photoreceptor cells. Affinity purification of in vivo-assembled SF-Whrna-associated protein complexes from retinal lysates followed by mass spectrometry, identified 19 novel candidate interaction partners, many of which are enriched in the synaptic proteome. Two human orthologs of the identified candidate interaction partners, FRMPD4 and Kir2.3, were validated as direct interaction partners of human whirlin using a yeast two-hybrid assay. The strong connection of whirlin with postsynaptic density proteins was not identified in previous in vitro protein-protein interaction assays, presumably due to the absence of a biologically relevant context. Isolation and identification of in vivo-assembled whirlin-associated protein complexes from the tissue of interest is therefore a powerful methodology to obtain novel insight into tissue specific protein-protein interactions and has the potential to improve significantly our understanding of the function of whirlin and the molecular pathogenesis underlying Usher syndrome type 2.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data.

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.

Compound Target Identification in Tissue-Specific Interaction Networks

Characterizing all possible side effects of compounds is not a trivial task due to unknown off-target proteins that might eventually lead lethal reactions. There is still a tremendous need of computational methods to identify protein targets of a new compound. We have performed a comprehensive analysis for identification of protein targets by integrating tissue-specific protein-protein interaction (PPI) networks and compound induced transcriptome data. Several network centrality metrics are computed to suggest the most probable off-targets of a given compound. The effects of interaction types between proteins and tissue-specific PPI networks are evaluated from multiple perspectives. Usage of network centrality metrics on a tissue-specific PPI network enhances the correct prediction rate of known targets of a given compound. The detailed analysis of successfully identified known targets indicated that degree and local radiality metrics are more practical for determining different types of target protein families, such as GPCR, chemokines, proteasome, and protein kinase families. Therefore, the proposed computational pipeline is applicable while investigating proteins from these families that can be especially targeted for treatment of complex diseases. To the best of our knowledge, this study for the first time presents how tissue-specific transcriptome changes alter the topological structure of PPI networks and the relative effects of tissue-specific networks in target identification. Overall, the proposed computational methods are practical tools for choosing more accurate protein targets for later computational analysis and wet-lab experiments.

Exploring Drug Treatment Patterns Based on the Action of Drug and Multilayer Network Model.

Some drugs can be used to treat multiple diseases, suggesting potential patterns in drug treatment. Determination of drug treatment patterns can improve our understanding of the mechanisms of drug action, enabling drug repurposing. A drug can be associated with a multilayer tissue-specific protein–protein interaction (TSPPI) network for the diseases it is used to treat. Proteins usually interact with other proteins to achieve functions that cause diseases. Hence, studying drug treatment patterns is similar to studying common module structures in multilayer TSPPI networks. Therefore, we propose a network-based model to study the treatment patterns of drugs. The method was designated SDTP (studying drug treatment pattern) and was based on drug effects and a multilayer network model. To demonstrate the application of the SDTP method, we focused on analysis of trichostatin A (TSA) in leukemia, breast cancer, and prostate cancer. We constructed a TSPPI multilayer network and obtained candidate drug-target modules from the network. Gene ontology analysis provided insights into the significance of the drug-target modules and co-expression networks. Finally, two modules were obtained as potential treatment patterns for TSA. Through analysis of the significance, composition, and functions of the selected drug-target modules, we validated the feasibility and rationality of our proposed SDTP method for identifying drug treatment patterns. In summary, our novel approach used a multilayer network model to overcome the shortcomings of single-layer networks and combined the network with information on drug activity. Based on the discovered drug treatment patterns, we can predict the potential diseases that the drug can treat. That is, if a disease-related protein module has a similar structure, then the drug is likely to be a potential drug for the treatment of the disease.

Tissue-Specific Subcellular Localization Prediction Using Multi-Label Markov Random Fields.

The understanding of subcellular localization (SCL) of proteins and proteome variation in the different tissues and organs of the human body are two crucial aspects for increasing our knowledge of the dynamic rules of proteins, the cell biology, and the mechanism of diseases. Although there have been tremendous contributions to these two fields independently, the lack of knowledge of the variation of spatial distribution of proteins in the different tissues still exists. Here, we proposed an approach that allows predicting protein SCL on tissue specificity through the use of tissue-specific functional associations and physical protein-protein interactions (PPIs). We applied our previously developed Bayesian collective Markov random fields (BCMRFs) on tissue-specific protein-protein interaction network (PPI network) for nine types of tissues focusing on eight high-level SCL. The evaluated results demonstrate the strength of our approach in predicting tissue-specific SCL. We identified 1,314 proteins that their SCL were previously proven cell line dependent. We predicted 549 novel tissue-specific localized candidate proteins while some of them were validated via text-mining.

INBIA: a boosting methodology for proteomic network inference

BackgroundThe analysis of tissue-specific protein interaction networks and their functional enrichment in pathological and normal tissues provides insights on the etiology of diseases. The Pan-cancer proteomic project, in The Cancer Genome Atlas, collects protein expressions in human cancers and it is a reference resource for the functional study of cancers. However, established protocols to infer interaction networks from protein expressions are still missing.ResultsWe have developed a methodology called Inference Network Based on iRefIndex Analysis (INBIA) to accurately correlate proteomic inferred relations to protein-protein interaction (PPI) networks. INBIA makes use of 14 network inference methods on protein expressions related to 16 cancer types. It uses as reference model the iRefIndex human PPI network. Predictions are validated through non-interacting and tissue specific PPI networks resources. The first, Negatome, takes into account likely non-interacting proteins by combining both structure properties and literature mining. The latter, TissueNet and GIANT, report experimentally verified PPIs in more than 50 human tissues. The reliability of the proposed methodology is assessed by comparing INBIA with PERA, a tool which infers protein interaction networks from Pathway Commons, by both functional and topological analysis.ConclusionResults show that INBIA is a valuable approach to predict proteomic interactions in pathological conditions starting from the current knowledge of human protein interactions.

Integrative network analysis highlights biological processes underlying GLP-1 stimulated insulin secretion: A DIRECT study.

Glucagon-like peptide 1 (GLP-1) stimulated insulin secretion has a considerable heritable component as estimated from twin studies, yet few genetic variants influencing this phenotype have been identified. We performed the first genome-wide association study (GWAS) of GLP-1 stimulated insulin secretion in non-diabetic individuals from the Netherlands Twin register (n = 126). This GWAS was enhanced using a tissue-specific protein-protein interaction network approach. We identified a beta-cell protein-protein interaction module that was significantly enriched for low gene scores based on the GWAS P-values and found support at the network level in an independent cohort from Tübingen, Germany (n = 100). Additionally, a polygenic risk score based on SNPs prioritized from the network was associated (P < 0.05) with glucose-stimulated insulin secretion phenotypes in up to 5,318 individuals in MAGIC cohorts. The network contains both known and novel genes in the context of insulin secretion and is enriched for members of the focal adhesion, extracellular-matrix receptor interaction, actin cytoskeleton regulation, Rap1 and PI3K-Akt signaling pathways. Adipose tissue is, like the beta-cell, one of the target tissues of GLP-1 and we thus hypothesized that similar networks might be functional in both tissues. In order to verify peripheral effects of GLP-1 stimulation, we compared the transcriptome profiling of ob/ob mice treated with liraglutide, a clinically used GLP-1 receptor agonist, versus baseline controls. Some of the upstream regulators of differentially expressed genes in the white adipose tissue of ob/ob mice were also detected in the human beta-cell network of genes associated with GLP-1 stimulated insulin secretion. The findings provide biological insight into the mechanisms through which the effects of GLP-1 may be modulated and highlight a potential role of the beta-cell expressed genes RYR2, GDI2, KIAA0232, COL4A1 and COL4A2 in GLP-1 stimulated insulin secretion.

Integrating phenotypic features and tissue-specific information to prioritize disease genes

Prioritization of candidate disease genes is crucial for improving medical care, and is one of the fundamental challenges in the post-genomic era. In recent years, different network-based methods for gene prioritization are proposed. Previous studies on gene prioritization show that tissue-specific protein-protein interaction (PPI) networks built by integrating PPIs with tissue-specific gene expression profiles can perform better than tissue-na¨ive global PPI network. Based on the observations that diseases with similar phenotypes are likely to have common related genes, and genes associated with the same phenotype tend to interact with each other, we propose a method to prioritize disease genes based on a heterogeneous network built by integrating phenotypic features and tissue-specific information. In this heterogeneous network, the PPI network is built by integrating phenotypic features with a tissue-specific PPI network, and the disease network consists of the diseases that are associated with the same phenotype and tissue as the query disease. To determine the impacts of these two factors on gene prioritization, we test three typical network-based prioritization methods on heterogeneous networks consisting of combinations of different PPIs and disease networks built with or without phenotypic features and tissue-specific information. We also compare the proposed method with other tissuespecific networks. The results of case studies reveals that integrating phenotypic features with a tissue-specific PPI network improves the prioritization results. Moreover, the disease networks generated using our method not only show comparable performance with the widely used disease similarity dataset of 5080 human diseases, but are also effective for diseases that are not in the dataset.

Integrated interactions database: tissue-specific view of the human and model organism interactomes.

IID (Integrated Interactions Database) is the first database providing tissue-specific protein–protein interactions (PPIs) for model organisms and human. IID covers six species (S. cerevisiae (yeast), C. elegans (worm), D. melonogaster (fly), R. norvegicus (rat), M. musculus (mouse) and H. sapiens (human)) and up to 30 tissues per species. Users query IID by providing a set of proteins or PPIs from any of these organisms, and specifying species and tissues where IID should search for interactions. If query proteins are not from the selected species, IID enables searches across species and tissues automatically by using their orthologs; for example, retrieving interactions in a given tissue, conserved in human and mouse. Interaction data in IID comprises three types of PPI networks: experimentally detected PPIs from major databases, orthologous PPIs and high-confidence computationally predicted PPIs. Interactions are assigned to tissues where their proteins pairs or encoding genes are expressed. IID is a major replacement of the I2D interaction database, with larger PPI networks (a total of 1,566,043 PPIs among 68,831 proteins), tissue annotations for interactions, and new query, analysis and data visualization capabilities. IID is available at http://ophid.utoronto.ca/iid.

SPECTRA: An Integrated Knowledge Base for Comparing Tissue and Tumor-Specific PPI Networks in Human.

Protein–protein interaction (PPI) networks available in public repositories usually represent relationships between proteins within the cell. They ignore the specific set of tissues or tumors where the interactions take place. Indeed, proteins can form tissue-selective complexes, while they remain inactive in other tissues. For these reasons, a great attention has been recently paid to tissue-specific PPI networks, in which nodes are proteins of the global PPI network whose corresponding genes are preferentially expressed in specific tissues. In this paper, we present SPECTRA, a knowledge base to build and compare tissue or tumor-specific PPI networks. SPECTRA integrates gene expression and protein interaction data from the most authoritative online repositories. We also provide tools for visualizing and comparing such networks, in order to identify the expression and interaction changes of proteins across tissues, or between the normal and pathological states of the same tissue. SPECTRA is available as a web server at http://alpha.dmi.unict.it/spectra.

Nodes occupying central positions in human tissue specific PPI networks are enriched with many splice variants.

The functional repertoire of genes in the eukaryotic organisms is enhanced by the phenomenon of alternative splicing. Hence, a node in a tissue specific protein-protein interaction (TS PPIN) network can be thought of as an ensemble of various spliced protein products of the corresponding gene expressed in that tissue. Here we demonstrate that the nodes that occupy topologically central positions characterized by high degree, betweenness, closeness, and eigenvector centrality values in TS PPINs of Homo sapiens are associated with high number of splice variants. We also show that the high "centrality" of these genes/nodes could in part be explained by the presence of a large number of promiscuous domains.

Comparative Analysis of Human Tissue Interactomes Reveals Factors Leading to Tissue-Specific Manifestation of Hereditary Diseases

An open question in human genetics is what underlies the tissue-specific manifestation of hereditary diseases, which are caused by genomic aberrations that are present in cells across the human body. Here we analyzed this phenomenon for over 300 hereditary diseases by using comparative network analysis. We created an extensive resource of protein expression and interactions in 16 main human tissues, by integrating recent data of gene and protein expression across tissues with data of protein-protein interactions (PPIs). The resulting tissue interaction networks (interactomes) shared a large fraction of their proteins and PPIs, and only a small fraction of them were tissue-specific. Applying this resource to hereditary diseases, we first show that most of the disease-causing genes are widely expressed across tissues, yet, enigmatically, cause disease phenotypes in few tissues only. Upon testing for factors that could lead to tissue-specific vulnerability, we find that disease-causing genes tend to have elevated transcript levels and increased number of tissue-specific PPIs in their disease tissues compared to unaffected tissues. We demonstrate through several examples that these tissue-specific PPIs can highlight disease mechanisms, and thus, owing to their small number, provide a powerful filter for interrogating disease etiologies. As two thirds of the hereditary diseases are associated with these factors, comparative tissue analysis offers a meaningful and efficient framework for enhancing the understanding of the molecular basis of hereditary diseases.

Global versus Local Hubs in Human Protein–Protein Interaction Network

In this study, we have constructed tissue-specific protein-protein interaction networks for 70 human tissues and have identified three types of hubs based on their expression breadths: (a) tissue-specific hubs (TSHs) (proteins that are expressed in ≤ 10 tissues and also form hubs in ≤ 10 tissues), (b) tissue-preferred hubs (TPHs) (proteins expressed in ≥ 60 tissues but are highly connected in ≤ 10 tissues), and (c) housekeeping hubs (HKHs) (proteins that are expressed in ≥ 60 tissues and also form hubs in ≥ 60 tissues). Comparative analyses revealed significant differences between TSHs and HKHs and also revealed that TPHs behave more like HKHs. TSHs are lengthier, more disordered, and also quickly evolving proteins as compared with HKHs. Despite having a similar number of binding surfaces and interacting domains, TSHs are associated with a lower degree of centrality as compared with HKHs, suggesting that TSHs are "unsaturated" with regard to their binding capability and are perhaps evolving with regard to their interactions. TSHs are less abundantly expressed as compared with HKHs and are enriched with PEST motifs, indicating their tight regulation. All of these properties of TSHs and HKHs correlate with their distinct functional roles; TSHs are involved in tissue-specific functional roles, viz., secretors, receptors, and signaling proteins, whereas HKHs are involved in core-cellular functions such as transcription, translation, and so on. Our study, therefore, brings forth a clear and distinct classification of hubs simply based on their expression breadth and further assumes significance in the light of the highly debated dichotomy of date and party hubs, which is based on the coexpression pattern of hubs with their partners.

Construction and Analysis of Tissue-Specific Protein-Protein Interaction Networks in Humans

We have studied the changes in protein-protein interaction network of 38 different tissues of the human body. 123 gene expression samples from these tissues were used to construct human protein-protein interaction network. This network is then pruned using the gene expression samples of each tissue to construct different protein-protein interaction networks corresponding to different studied tissues of the body. This study is helpful for understanding how human protein interactions change in different tissues. In this way, similar tissues of the body and special functions of each tissue, corresponding to their individualized subnetworks, can be identified. We have calculated graph parameters for studying these protein-protein interaction networks and hubs and non-hubs of the studied protein-protein interaction networks are identified. We found a common subnetwork among protein-protein interaction networks of the studied tissues and a tree of tissue similarities has been constructed. We have also found that average correlation coefficient of hubs in human protein-protein interaction networks obeys a normal-like distribution though it is not possible to separate party and date hubs.

Walking on a Tissue-Specific Disease-Protein-Complex Heterogeneous Network for the Discovery of Disease-Related Protein Complexes

Besides the pinpointing of individual disease-related genes, associating protein complexes to human inherited diseases is also of great importance, because a biological function usually arises from the cooperative behaviour of multiple proteins in a protein complex. Moreover, knowledge about disease-related protein complexes could also enhance the inference of disease genes and pathogenic genetic variants. Here, we have designed a computational systems biology approach to systematically analyse potential relationships between diseases and protein complexes. First, we construct a heterogeneous network which is composed of a disease-disease similarity layer, a tissue-specific protein-protein interaction layer, and a protein complex membership layer. Then, we propose a random walk model on this disease-protein-complex network for identifying protein complexes that are related to a query disease. With a series of leave-one-out cross-validation experiments, we show that our method not only possesses high performance but also demonstrates robustness regarding the parameters and the network structure. We further predict a landscape of associations between human diseases and protein complexes. This landscape can be used to facilitate the inference of disease genes, thereby benefiting studies on pathology of diseases.

Enhancing the Prioritization of Disease-Causing Genes through Tissue Specific Protein Interaction Networks

The prioritization of candidate disease-causing genes is a fundamental challenge in the post-genomic era. Current state of the art methods exploit a protein-protein interaction (PPI) network for this task. They are based on the observation that genes causing phenotypically-similar diseases tend to lie close to one another in a PPI network. However, to date, these methods have used a static picture of human PPIs, while diseases impact specific tissues in which the PPI networks may be dramatically different. Here, for the first time, we perform a large-scale assessment of the contribution of tissue-specific information to gene prioritization. By integrating tissue-specific gene expression data with PPI information, we construct tissue-specific PPI networks for 60 tissues and investigate their prioritization power. We find that tissue-specific PPI networks considerably improve the prioritization results compared to those obtained using a generic PPI network. Furthermore, they allow predicting novel disease-tissue associations, pointing to sub-clinical tissue effects that may escape early detection.

Dynamics of Homologous and Heterologous Calcium Channel α1 and β Subunits in the Skeletal Muscle Triad

Skeletal muscle EC-coupling requires direct interactions of the dihydropyridine receptor (DHPR) with the ryanodine receptor (RyR1) and its organization in tetrads opposite RyR1. Accordingly, tetrad formation and skeletal muscle EC-coupling fail in myotubes expressing the cardiac α1C subunit or lacking the β1a subunit. This indicates that tissue-specific protein-protein interactions between α1S, β1a and RyR1 are important for the structure and function of the skeletal muscle EC-coupling apparatus. Here we apply fluorescence recovery after photobleaching (FRAP) in dysgenic (α1S-/-) myotubes reconstituted with GFP-tagged DHPR subunits to probe the importance of these interactions in stabilizing them in the triad. Consistent with its anchoring in the triad, the α1S subunit is essentially immobile. Surprisingly, also the cardiac α1C subunit, which does not interact with RyR1, is immobile. Thus, the II-III loop - RyR1 interaction is not the determing factor for stabilizing α1S in the triad. When coexpressed with α1S, the β1a subunit showed a higher fluorescence recovery rate than the α1 subunits, suggesting that it can be dynamically exchanged in the complex. Again, the recovery did not increase when β1a was coexpressed with α1C, indicating that its putative interaction with RyR1 does not contribute to the stability of β1a in the complex. Previously we showed that all βisoforms can associate with α1S and α1C in the triads. When the non-muscle β2a or β4b subunits were expressed with α1S in dysgenic myotubes their recovery after photobleaching was significantly higher than that of β1a, demonstrating that the homologous β subunit binds the channel complex with higher affinity than heterologous β subunits. Funding: FWF P20059, W1101, T443-B18.