Comparative analysis of human tissue interactomes reveals factors leading to tissue-specific manifestation of hereditary diseases

Affected cell types for hundreds of Mendelian diseases revealed by analysis of human and mouse single-cell data

Mendelian diseases tend to manifest clinically in certain tissues, yet their affected cell types typically remain elusive. Single-cell expression studies showed that overexpression of disease-associated genes may point to the affected cell types. Here, we developed a method that infers disease-affected cell types from the preferential expression of disease-associated genes in cell types (PrEDiCT). We applied PrEDiCT to single-cell expression data of six human tissues, to infer the cell types affected in Mendelian diseases. Overall, we inferred the likely affected cell types for 328 diseases. We corroborated our findings by literature text-mining, expert validation, and recapitulation in mouse corresponding tissues. Based on these findings, we explored characteristics of disease-affected cell types, showed that diseases manifesting in multiple tissues tend to affect similar cell types, and highlighted cases where gene functions could be used to refine inference. Together, these findings expand the molecular understanding of disease mechanisms and cellular vulnerability.

Gene Interactions in Human Disease Studies-Evidence Is Mounting

Despite monumental advances in molecular technology to generate genome sequence data at scale, there is still a considerable proportion of heritability in most complex diseases that remains unexplained. Because many of the discoveries have been single-nucleotide variants with small to moderate effects on disease, the functional implication of many of the variants is still unknown and, thus, we have limited new drug targets and therapeutics. We, and many others, posit that one primary factor that has limited our ability to identify novel drug targets from genome-wide association studies may be due to gene interactions (epistasis), gene-environment interactions, network/pathway effects, or multiomic relationships. We propose that many of these complex models explain much of the underlying genetic architecture of complex disease. In this review, we discuss the evidence from multiple research avenues, ranging from pairs of alleles to multiomic integration studies and pharmacogenomics, that supports the need for further investigation of gene interactions (or epistasis) in genetic and genomic studies of human disease. Our goal is to catalog the mounting evidence for epistasis in genetic studies and the connections between genetic interactions and human health and disease that could enable precision medicine of the future.

Predicting molecular mechanisms of hereditary diseases by using their tissue-selective manifestation

How do aberrations in widely expressed genes lead to tissue-selective hereditary diseases? Previous attempts to answer this question were limited to testing a few candidate mechanisms. To answer this question at a larger scale, we developed "Tissue Risk Assessment of Causality by Expression" (TRACE), a machine learning approach to predict genes that underlie tissue-selective diseases and selectivity-related features. TRACE utilized 4,744 biologically interpretable tissue-specific gene features that were inferred from heterogeneous omics datasets. Application of TRACE to 1,031 disease genes uncovered known and novel selectivity-related features, the most common of which was previously overlooked. Next, we created a catalog of tissue-associated risks for 18,927 protein-coding genes (https://netbio.bgu.ac.il/trace/). As proof-of-concept, we prioritized candidate disease genes identified in 48 rare-disease patients. TRACE ranked the verified disease gene among the patient's candidate genes significantly better than gene prioritization methods that rank by gene constraint or tissue expression. Thus, tissue selectivity combined with machine learning enhances genetic and clinical understanding of hereditary diseases.

Splicing Modulators Are Involved in Human Polyglutamine Diversification via Protein Complexes Shuttling between Nucleus and Cytoplasm

Length polymorphisms of polyglutamine (polyQs) in triplet-repeat-disease-causing genes have diversified during primate evolution despite them conferring a risk of human-specific diseases. To explain the evolutionary process of this diversification, there is a need to focus on mechanisms by which rapid evolutionary changes can occur, such as alternative splicing. Proteins that can bind polyQs are known to act as splicing factors and may provide clues about the rapid evolutionary process. PolyQs are also characterized by the formation of intrinsically disordered (ID) regions, so I hypothesized that polyQs are involved in the transportation of various molecules between the nucleus and cytoplasm to regulate mechanisms characteristic of humans such as neural development. To determine target molecules for empirical research to understand the evolutionary change, I explored protein-protein interactions (PPIs) involving the relevant proteins. This study identified pathways related to polyQ binding as hub proteins scattered across various regulatory systems, including regulation via PQBP1, VCP, or CREBBP. Nine ID hub proteins with both nuclear and cytoplasmic localization were found. Functional annotations suggested that ID proteins containing polyQs are involved in regulating transcription and ubiquitination by flexibly changing PPI formation. These findings explain the relationships among splicing complex, polyQ length variations, and modifications in neural development.

Homo cerevisiae-Leveraging Yeast for Investigating Protein-Protein Interactions and Their Role in Human Disease

Understanding how genetic variation affects phenotypes represents a major challenge, particularly in the context of human disease. Although numerous disease-associated genes have been identified, the clinical significance of most human variants remains unknown. Despite unparalleled advances in genomics, functional assays often lack sufficient throughput, hindering efficient variant functionalization. There is a critical need for the development of more potent, high-throughput methods for characterizing human genetic variants. Here, we review how yeast helps tackle this challenge, both as a valuable model organism and as an experimental tool for investigating the molecular basis of phenotypic perturbation upon genetic variation. In systems biology, yeast has played a pivotal role as a highly scalable platform which has allowed us to gain extensive genetic and molecular knowledge, including the construction of comprehensive interactome maps at the proteome scale for various organisms. By leveraging interactome networks, one can view biology from a systems perspective, unravel the molecular mechanisms underlying genetic diseases, and identify therapeutic targets. The use of yeast to assess the molecular impacts of genetic variants, including those associated with viral interactions, cancer, and rare and complex diseases, has the potential to bridge the gap between genotype and phenotype, opening the door for precision medicine approaches and therapeutic development.

Biological interacting units identified in human protein networks reveal tissue-functional diversification and its impact on disease

•

Biological processes are exerted by groups of physically interacting proteins.

•

Proteins display variable biological roles depending on tissue-interactomic context.

•

Tissue-specific protein-protein interaction networks reveal functional diversification.

•

Most disease associated genes/proteins display tissue-specific phenotypes.

•

Protein interaction network analysis is a valuable resource to identify disease genes.

Protein-protein interactions (PPI) play an essential role in the biological processes that occur in the cell. Therefore, the dissection of PPI networks becomes decisive to model functional coordination and predict pathological de-regulation. Cellular networks are dynamic and proteins display varying roles depending on the tissue-interactomic context. Thus, the use of centrality measures in individual proteins fall short to dissect the functional properties of the cell. For this reason, there is a need for more comprehensive, relational, and context-specific ways to analyze the multiple actions of proteins in different cells and identify specific functional assemblies within global biomolecular networks. Under this framework, we define Biological Interacting units (BioInt-U) as groups of proteins that interact physically and are enriched in a common Gene Ontology. A search strategy was applied on 33 tissue-specific (TS) PPI networks to generate BioInt libraries associated with each particular human tissue. The cross-tissue comparison showed that housekeeping assemblies incorporate different proteins and exhibit distinct network properties depending on the tissue. Furthermore, disease genes (DGs) of tissue-associated pathologies preferentially accumulate in units in the expected tissues, which in turn were more central in the TS networks. Overall, the study reveals a tissue-specific functional diversification based on the identification of specific protein units and suggests vulnerabilities specific of each tissue network, which can be applied to refine protein-disease association methods.

Balancing the transcriptome: leveraging sample similarity to improve measures of gene specificity

The spatial and temporal domain of a gene's expression can range from ubiquitous to highly specific. Quantifying the degree to which this expression is unique to a specific tissue or developmental timepoint can provide insight into the etiology of genetic diseases. However, quantifying specificity remains challenging as measures of specificity are sensitive to similarity between samples in the sample set. For example, in the Gene-Tissue Expression project (GTEx), brain subregions are overrepresented at 13 of 54 (24%) unique tissues sampled. In this dataset, existing specificity measures have a decreased ability to identify genes specific to the brain relative to other organs. To solve this problem, we leverage sample similarity information to weight samples such that overrepresented tissues do not have an outsized effect on specificity estimates. We test this reweighting procedure on 4 measures of specificity, Z-score, Tau, Tsi and Gini, in the GTEx data and in single cell datasets for zebrafish and mouse. For all of these measures, incorporating sample similarity information to weight samples results in greater stability of sets of genes called as specific and decreases the overall variance in the change of specificity estimates as sample sets become more unbalanced. Furthermore, the genes with the largest improvement in their specificity estimate's stability are those with functions related to the overrepresented sample types. Our results demonstrate that incorporating similarity information improves specificity estimates' stability to the choice of the sample set used to define the transcriptome, providing more robust and reproducible measures of specificity for downstream analyses.

The Organ-Disease Annotations (ODiseA) database of hereditary diseases and inflicted tissues

Hereditary diseases tend to manifest clinically in few selected tissues. Knowledge of those tissues is important for better understanding of disease mechanisms, which often remain elusive. However, information on the tissues inflicted by each disease is not easily obtainable. Well-established resources, such as the Online Mendelian Inheritance in Man (OMIM) database and Human Phenotype Ontology (HPO), report on a spectrum of disease manifestations, yet do not highlight the main inflicted tissues. The Organ-Disease Annotations (ODiseA) database contains 4,357 thoroughly-curated annotations for 2,181 hereditary diseases and 45 inflicted tissues. Additionally, ODiseA reports 692 annotations of 635 diseases and the pathogenic tissues where they emerge. ODiseA can be queried by disease, disease gene, or inflicted tissue. Owing to its expansive, high-quality annotations, ODiseA serves as a valuable and unique tool for biomedical and computational researchers studying genotype-phenotype relationships of hereditary diseases. ODiseA is available at https://netbio.bgu.ac.il/odisea.

The differential activity of biological processes in tissues and cell subsets can illuminate disease-related processes and cell-type identities

MOTIVATION

The distinct functionalities of human tissues and cell types underlie complex phenotype-genotype relationships, yet often remain elusive. Harnessing the multitude of bulk and single-cell human transcriptomes while focusing on processes can help reveal these distinct functionalities.

RESULTS

The Tissue-Process Activity (TiPA) method aims to identify processes that are preferentially active or under-expressed in specific contexts, by comparing the expression levels of process genes between contexts. We tested TiPA on 1579 tissue-specific processes and bulk tissue transcriptomes, finding that it performed better than another method. Next, we used TiPA to ask whether the activity of certain processes could underlie the tissue-specific manifestation of 1233 hereditary diseases. We found that 21% of the disease-causing genes indeed participated in such processes, thereby illuminating their genotype-phenotype relationships. Lastly, we applied TiPA to single-cell transcriptomes of 108 human cell types, revealing that process activities often match cell-type identities and can thus aid annotation efforts. Hence, differential activity of processes can highlight the distinct functionality of tissues and cells in a robust and meaningful manner.

AVAILABILITY AND IMPLEMENTATION

TiPA code is available in GitHub (https://github.com/moranshar/TiPA). In addition, all data are available as part of the Supplementary Material.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcription factor protein interactomes reveal genetic determinants in heart disease

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

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.

An atlas of protein-protein interactions across mouse tissues

Cellular processes arise from the dynamic organization of proteins in networks of physical interactions. Mapping the interactome has therefore been a central objective of high-throughput biology. However, the dynamics of protein interactions across physiological contexts remain poorly understood. Here, we develop a quantitative proteomic approach combining protein correlation profiling with stable isotope labeling of mammals (PCP-SILAM) to map the interactomes of seven mouse tissues. The resulting maps provide a proteome-scale survey of interactome rewiring across mammalian tissues, revealing more than 125,000 unique interactions at a quality comparable to the highest-quality human screens. We identify systematic suppression of cross-talk between the evolutionarily ancient housekeeping interactome and younger, tissue-specific modules. Rewired proteins are tightly regulated by multiple cellular mechanisms and are implicated in disease. Our study opens up new avenues to uncover regulatory mechanisms that shape in vivo interactome responses to physiological and pathophysiological stimuli in mammalian systems.

Disease-specific interactome alterations via epichaperomics: the case for Alzheimer's disease

The increasingly appreciated prevalence of complicated stressor-to-phenotype associations in human disease requires a greater understanding of how specific stressors affect systems or interactome properties. Many currently untreatable diseases arise due to variations in, and through a combination of, multiple stressors of genetic, epigenetic, and environmental nature. Unfortunately, how such stressors lead to a specific disease phenotype or inflict a vulnerability to some cells and tissues but not others remains largely unknown and unsatisfactorily addressed. Analysis of cell- and tissue-specific interactome networks may shed light on organization of biological systems and subsequently to disease vulnerabilities. However, deriving human interactomes across different cell and disease contexts remains a challenge. To this end, this opinion article links stressor-induced protein interactome network perturbations to the formation of pathologic scaffolds termed epichaperomes, revealing a viable and reproducible experimental solution to obtaining rigorous context-dependent interactomes. This article presents our views on how a specialized 'omics platform called epichaperomics may complement and enhance the currently available conventional approaches and aid the scientific community in defining, understanding, and ultimately controlling interactome networks of complex diseases such as Alzheimer's disease. Ultimately, this approach may aid the transition from a limited single-alteration perspective in disease to a comprehensive network-based mindset, which we posit will result in precision medicine paradigms for disease diagnosis and treatment.

The landscape of molecular chaperones across human tissues reveals a layered architecture of core and variable chaperones

The sensitivity of the protein-folding environment to chaperone disruption can be highly tissue-specific. Yet, the organization of the chaperone system across physiological human tissues has received little attention. Through computational analyses of large-scale tissue transcriptomes, we unveil that the chaperone system is composed of core elements that are uniformly expressed across tissues, and variable elements that are differentially expressed to fit with tissue-specific requirements. We demonstrate via a proteomic analysis that the muscle-specific signature is functional and conserved. Core chaperones are significantly more abundant across tissues and more important for cell survival than variable chaperones. Together with variable chaperones, they form tissue-specific functional networks. Analysis of human organ development and aging brain transcriptomes reveals that these functional networks are established in development and decline with age. In this work, we expand the known functional organization of de novo versus stress-inducible eukaryotic chaperones into a layered core-variable architecture in multi-cellular organisms.

FunCoup 5: Functional Association Networks in All Domains of Life, Supporting Directed Links and Tissue-Specificity

FunCoup (https://funcoup.sbc.su.se) is one of the most comprehensive functional association networks of genes/proteins available. Functional associations are inferred by integrating different types of evidence using a redundancy-weighted naïve Bayesian approach, combined with orthology transfer. FunCoup's high coverage comes from using eleven different types of evidence, and extensive transfer of information between species. Since the latest update of the database, the availability of source data has improved drastically, and user expectations on a tool for functional associations have grown. To meet these requirements, we have made a new release of FunCoup with updated source data and improved functionality. FunCoup 5 now includes 22 species from all domains of life, and the source data for evidences, gold standards, and genomes have been updated to the latest available versions. In this new release, directed regulatory links inferred from transcription factor binding can be visualized in the network viewer for the human interactome. Another new feature is the possibility to filter by genes expressed in a certain tissue in the network viewer. FunCoup 5 further includes the SARS-CoV-2 proteome, allowing users to visualize and analyze interactions between SARS-CoV-2 and human proteins in order to better understand COVID-19. This new release of FunCoup constitutes a major advance for the users, with updated sources, new species and improved functionality for analysis of the networks.

Integrative Epigenomic Analysis of Transcriptional Regulation of Human CircRNAs

Circular RNAs (circRNAs) are evolutionarily conserved and abundant non-coding RNAs whose functions and regulatory mechanisms remain largely unknown. Here, we identify and characterize an epigenomically distinct group of circRNAs (TAH-circRNAs), which are transcribed to a higher level than their host genes. By integrative analysis of cistromic and transcriptomic data, we find that compared with other circRNAs, TAH-circRNAs are expressed more abundantly and have more transcription factors (TFs) binding sites and lower DNA methylation levels. Concordantly, TAH-circRNAs are enriched in open and active chromatin regions. Importantly, ChIA-PET results showed that 23–52% of transcription start sites (TSSs) of TAH-circRNAs have direct interactions with cis-regulatory regions, strongly suggesting their independent transcriptional regulation from host genes. In addition, we characterize molecular features of super-enhancer-driven circRNAs in cancer biology. Together, this study comprehensively analyzes epigenomic characteristics of circRNAs and identifies a distinct group of TAH-circRNAs that are independently transcribed via enhancers and super-enhancers by TFs. These findings substantially advance our understanding of the regulatory mechanism of circRNAs and may have important implications for future investigations of this class of non-coding RNAs.

An analytical method for the identification of cell type-specific disease gene modules

Background

Genome-wide association studies have identified genetic variants associated with the risk of brain-related diseases, such as neurological and psychiatric disorders, while the causal variants and the specific vulnerable cell types are often needed to be studied. Many disease-associated genes are expressed in multiple cell types of human brains, while the pathologic variants affect primarily specific cell types. We hypothesize a model in which what determines the manifestation of a disease in a cell type is the presence of disease module comprised of disease-associated genes, instead of individual genes. Therefore, it is essential to identify the presence/absence of disease gene modules in cells.

Methods

To characterize the cell type-specificity of brain-related diseases, we construct human brain cell type-specific gene interaction networks integrating human brain nucleus gene expression data with a referenced tissue-specific gene interaction network. Then from the cell type-specific gene interaction networks, we identify significant cell type-specific disease gene modules by performing statistical tests.

Results

Between neurons and glia cells, the constructed cell type-specific gene networks and their gene functions are distinct. Then we identify cell type-specific disease gene modules associated with autism spectrum disorder and find that different gene modules are formed and distinct gene functions may be dysregulated in different cells. We also study the similarity and dissimilarity in cell type-specific disease gene modules among autism spectrum disorder, schizophrenia and bipolar disorder. The functions of neurons-specific disease gene modules are associated with synapse for all three diseases, while those in glia cells are different. To facilitate the use of our method, we develop an R package, CtsDGM, for the identification of cell type-specific disease gene modules.

Conclusions

The results support our hypothesis that a disease manifests itself in a cell type through forming a statistically significant disease gene module. The identification of cell type-specific disease gene modules can promote the development of more targeted biomarkers and treatments for the disease. Our method can be applied for depicting the cell type heterogeneity of a given disease, and also for studying the similarity and dissimilarity between different disorders, providing new insights into the molecular mechanisms underlying the pathogenesis and progression of diseases.

Dosage-sensitive molecular mechanisms are associated with the tissue-specificity of traits and diseases

Hereditary diseases and complex traits often manifest in specific tissues, whereas their causal genes are expressed in many tissues that remain unaffected. Among the mechanisms that have been suggested for this enigmatic phenomenon is dosage-sensitive compensation by paralogs of causal genes. Accordingly, tissue-selectivity stems from dosage imbalance between causal genes and paralogs that occurs particularly in disease-susceptible tissues. Here, we used a large-scale dataset of thousands of tissue transcriptomes and applied a linear mixed model (LMM) framework to assess this and other dosage-sensitive mechanisms. LMM analysis of 382 hereditary diseases consistently showed evidence for dosage-sensitive compensation by paralogs across diseases subsets and susceptible tissues. LMM analysis of 135 candidate genes that are strongly associated with 16 tissue-selective complex traits revealed a similar tendency among half of the trait-associated genes. This suggests that dosage-sensitive compensation by paralogs affects the tissue-selectivity of complex traits, and can be used to illuminate candidate genes' modes of action. Next, we applied LMM to analyze dosage imbalance between causal genes and three classes of genetic modifiers, including regulatory micro-RNAs, pseudogenes, and genetic interactors. Our results propose modifiers as a fundamental axis in tissue-selectivity of diseases and traits, and demonstrates the power of LMM as a statistical framework for discovering treatment avenues.

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.

Precision medicine - networks to the rescue

Genetic variants are often not predictive of the phenotypic outcome. Individuals carrying the same pathogenic variant, associated with Mendelian or complex disease, can manifest to different extents, from severe-to-mild to no disease. Improving the accuracy of predicted clinical manifestations of genetic variants has emerged as one of the biggest challenges in precision medicine, which can only be addressed by understanding the mechanisms underlying genotype-phenotype relationships. Efforts to understand the molecular basis of these relationships have identified complex systems of interacting biomolecules that underlie cellular function. Here, we review recent advances in how modeling cellular systems as networks of interacting proteins has fueled identification of disease-associated processes, delineation of underlying molecular mechanisms, and prediction of the pathogenicity of variants. This review is intended to be inspiring for clinicians, geneticists, and network biologists alike who aim to jointly advance our understanding of human disease and accelerate progress toward precision medicine.

ResponseNet v.3: revealing signaling and regulatory pathways connecting your proteins and genes across human tissues

ResponseNet v.3 is an enhanced version of ResponseNet, a web server that is designed to highlight signaling and regulatory pathways connecting user-defined proteins and genes by using the ResponseNet network optimization approach (http://netbio.bgu.ac.il/respnet). Users run ResponseNet by defining source and target sets of proteins, genes and/or microRNAs, and by specifying a molecular interaction network (interactome). The output of ResponseNet is a sparse, high-probability interactome subnetwork that connects the two sets, thereby revealing additional molecules and interactions that are involved in the studied condition. In recent years, massive efforts were invested in profiling the transcriptomes of human tissues, enabling the inference of human tissue interactomes. ResponseNet v.3 expands ResponseNet2.0 by harnessing ∼11,600 RNA-sequenced human tissue profiles made available by the Genotype-Tissue Expression consortium, to support context-specific analysis of 44 human tissues. Thus, ResponseNet v.3 allows users to illuminate the signaling and regulatory pathways potentially active in the context of a specific tissue, and to compare them with active pathways in other tissues. In the era of precision medicine, such analyses open the door for tissue- and patient-specific analyses of pathways and diseases.

Rare Genetic Diseases: Nature's Experiments on Human Development

Rare genetic diseases are the result of a continuous forward genetic screen that nature is conducting on humans. Here, we present epistemological and systems biology arguments highlighting the importance of studying these rare genetic diseases. We contend that the expanding catalog of mutations in ∼4,000 genes, which cause ∼6,500 diseases and their annotated phenotypes, offer a wide landscape for discovering fundamental mechanisms required for human development and involved in common diseases. Rare afflictions disproportionately affect the nervous system in children, but paradoxically, the majority of these disease-causing genes are evolutionarily ancient and ubiquitously expressed in human tissues. We propose that the biased prevalence of childhood rare diseases affecting nervous tissue results from the topological complexity of the protein interaction networks formed by ubiquitous and ancient proteins encoded by childhood disease genes. Finally, we illustrate these principles discussing Menkes disease, an example of the discovery power afforded by rare diseases.

A reference map of the human binary protein interactome

Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype-phenotype relationships^1,2. Here, we present a human “all-by-all” reference interactome map of human binary protein interactions, or “HuRI”. With ~53,000 high-quality protein-protein interactions (PPIs), HuRI has approximately four times more such interactions than high-quality curated interactions from small-scale studies. Integrating HuRI with genome³, transcriptome⁴, and proteome⁵ data enables the study of cellular function within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying specific subcellular roles of PPIs. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms underlying tissue-specific phenotypes of Mendelian diseases. HuRI represents a systematic proteome-wide reference linking genomic variation to phenotypic outcomes.

Differential network analysis of multiple human tissue interactomes highlights tissue-selective processes and genetic disorder genes

Motivation

Differential network analysis, designed to highlight network changes between conditions, is an important paradigm in network biology. However, differential network analysis methods have been typically designed to compare between two conditions and were rarely applied to multiple protein interaction networks (interactomes). Importantly, large-scale benchmarks for their evaluation have been lacking.

Results

Here, we present a framework for assessing the ability of differential network analysis of multiple human tissue interactomes to highlight tissue-selective processes and disorders. For this, we created a benchmark of 6499 curated tissue-specific Gene Ontology biological processes. We applied five methods, including four differential network analysis methods, to construct weighted interactomes for 34 tissues. Rigorous assessment of this benchmark revealed that differential analysis methods perform well in revealing tissue-selective processes (AUCs of 0.82–0.9). Next, we applied differential network analysis to illuminate the genes underlying tissue-selective hereditary disorders. For this, we curated a dataset of 1305 tissue-specific hereditary disorders and their manifesting tissues. Focusing on subnetworks containing the top 1% differential interactions in disease-relevant tissue interactomes revealed significant enrichment for disorder-causing genes in 18.6% of the cases, with a significantly high success rate for blood, nerve, muscle and heart diseases.

Summary

Altogether, we offer a framework that includes expansive manually curated datasets of tissue-selective processes and disorders to be used as benchmarks or to illuminate tissue-selective processes and genes. Our results demonstrate that differential analysis of multiple human tissue interactomes is a powerful tool for highlighting processes and genes with tissue-selective functionality and clinical impact.

Availability and implementation

Datasets are available as part of the Supplementary data.

Supplementary information

Supplementary data are available at Bioinformatics online.

DESE: estimating driver tissues by selective expression of genes associated with complex diseases or traits

The driver tissues or cell types in which susceptibility genes initiate diseases remain elusive. We develop a unified framework to detect the causal tissues of complex diseases or traits according to selective expression of disease-associated genes in genome-wide association studies (GWASs). This framework consists of three components which run iteratively to produce a converged prioritization list of driver tissues. Additionally, this framework also outputs a list of prioritized genes as a byproduct. We apply the framework to six representative complex diseases or traits with GWAS summary statistics, which leads to the estimation of the lung as an associated tissue of rheumatoid arthritis.

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.

The first Autumn School on Proteostasis: from molecular mechanisms to organismal consequences

The first Autumn School on Proteostasis was held at the Mediterranean Institute for Life Sciences (MedILS) in Split, Croatia, from November 12th-16th, 2018, bringing together 44 graduate students and postdoctoral fellows and 22 principal investigators from around the world. This meeting was geared towards providing students with an in-depth understanding of the field of proteostasis, with the aim of broadening their perspectives of the field. Session topics covered multiple aspects of cellular and organismal proteostasis, including fundamental principles, responses to heat shock, aging and disease, and protein folding, misfolding, and degradation. The structure of the meeting and the restricted number of participants afforded the students and postdocs the opportunity to interact with principal investigators to discuss not only their latest research, but also their career prospects and progression in a close, supportive environment.

Role of duplicate genes in determining the tissue-selectivity of hereditary diseases

A longstanding puzzle in human genetics is what limits the clinical manifestation of hundreds of hereditary diseases to certain tissues, while their causal genes are expressed throughout the human body. A general conception is that tissue-selective disease phenotypes emerge when masking factors operate in unaffected tissues, but are specifically absent or insufficient in disease-manifesting tissues. Although this conception has critical impact on the understanding of disease manifestation, it was never challenged in a systematic manner across a variety of hereditary diseases and affected tissues. Here, we address this gap in our understanding via rigorous analysis of the susceptibility of over 30 tissues to 112 tissue-selective hereditary diseases. We focused on the roles of paralogs of causal genes, which are presumably capable of compensating for their aberration. We show for the first time at large-scale via quantitative analysis of omics datasets that, preferentially in the disease-manifesting tissues, paralogs are under-expressed relative to causal genes in more than half of the diseases. This was observed for several susceptible tissues and for causal genes with varying number of paralogs, suggesting that imbalanced expression of paralogs increases tissue susceptibility. While for many diseases this imbalance stemmed from up-regulation of the causal gene in the disease-manifesting tissue relative to other tissues, it was often combined with down-regulation of its paralog. Notably in roughly 20% of the cases, this imbalance stemmed only from significant down-regulation of the paralog. Thus, dosage relationships between paralogs appear as important, yet currently under-appreciated, modifiers of disease manifestation.

The DifferentialNet database of differential protein-protein interactions in human tissues

DifferentialNet is a novel database that provides users with differential interactome analysis of human tissues (http://netbio.bgu.ac.il/diffnet/). Users query DifferentialNet by protein, and retrieve its differential protein-protein interactions (PPIs) per tissue via an interactive graphical interface. To compute differential PPIs, we integrated available data of experimentally detected PPIs with RNA-sequencing profiles of tens of human tissues gathered by the Genotype-Tissue Expression consortium (GTEx) and by the Human Protein Atlas (HPA). We associated each PPI with a score that reflects whether its corresponding genes were expressed similarly across tissues, or were up- or down-regulated in the selected tissue. By this, users can identify tissue-specific interactions, filter out PPIs that are relatively stable across tissues, and highlight PPIs that show relative changes across tissues. The differential PPIs can be used to identify tissue-specific processes and to decipher tissue-specific phenotypes. Moreover, they unravel processes that are tissue-wide yet tailored to the specific demands of each tissue.

Induced Pluripotent Stem Cells in Disease Modeling and Gene Identification

Experimental modeling of human inherited disorders provides insight into the cellular and molecular mechanisms involved, and the underlying genetic component influencing, the disease phenotype. The breakthrough development of induced pluripotent stem cell (iPSC) technology represents a quantum leap in experimental modeling of human diseases, providing investigators with a self-renewing and, thus, unlimited source of pluripotent cells for targeted differentiation. In principle, the entire range of cell types found in the human body can be interrogated using an iPSC approach. Therefore, iPSC technology, and the increasingly refined abilities to differentiate iPSCs into disease-relevant target cells, has far-reaching implications for understanding disease pathophysiology, identifying disease-causing genes, and developing more precise therapeutics, including advances in regenerative medicine. In this chapter, we discuss the technological perspectives and recent developments in the application of patient-derived iPSC lines for human disease modeling and disease gene identification.

Distinct distributions of genomic features of the 5’ and 3’ partners of coding somatic cancer gene fusions: arising mechanisms and functional implications

The genomic features and arising mechanisms of coding cancer somatic gene fusions (CSGFs) largely remain elusive. In this study, we show the gene origin stratification pattern of CSGF partners that fusion partners in human cancers are significantly enriched for genes with the gene age ofEuteleostomes and with the gene family age of Bilateria. GC skew (a measurement of G, C nucleotide content bias, (G-C)/(G+C)) is a useful measurement to indicate the DNA leading strand, lagging strand, replication origin, and replication terminal and DNA-RNA R-loop formation. We find that GC skew bias at the 5 prime (5′) but not the 3 prime (3’) partners of CSGFs, coincident with the polarity feature of gene expression breadth that the 5’ partners are more ubiquitous while the 3’ fusion partners are more tissue specific in general. We reveal distinct length and composition distributions of 5’ and 3’ of CSGFs, including sequence features corresponded to the 5’ untranslated regions (UTRs), 3’ UTRs, and the N-terminal sequences of the encoded proteins. Oncogenic somatic gene fusions are most enriched for the 5’ and 3’ genes’ somatic amplification alongside a substantial proportion of other types of combinations. At the function level, 5’ partners of CSGFs appear more likely to be tumour suppressor genes while many 3’ partners appear to be proto-oncogene. Such distinct polarities of CSGFs at the evolutionary, structural, genomic and functional levels indicate the heterogeneous arsing mechanisms of CSGFs including R-loops and suggest potential novel targeted therapeutics specific to CSGF functional categories.

Comprehensive Analysis of Tissue-wide Gene Expression and Phenotype Data Reveals Tissues Affected in Rare Genetic Disorders

Linking putatively pathogenic variants to the tissues they affect is necessary for determining the correct diagnostic workup and therapeutic regime in undiagnosed patients. Here, we explored how gene expression across healthy tissues can be used to infer this link. We integrated 6,665 tissue-wide transcriptomes with genetic disorder knowledge bases covering 3,397 diseases. Receiver-operating characteristics (ROC) analysis using expression levels in each tissue and across tissues indicated significant but modest associations between elevated expression and phenotype for most tissues (maximum area under ROC curve = 0.69). At extreme elevation, associations were marked. Upregulation of disease genes in affected tissues was pronounced for genes associated with autosomal dominant over recessive disorders. Pathways enriched for genes expressed and associated with phenotypes highlighted tissue functionality, including lipid metabolism in spleen and DNA repair in adipose tissue. These results suggest features useful for evaluating the likelihood of particular tissue manifestations in genetic disorders. The web address of an interactive platform integrating these data is provided.

Rewiring of the inferred protein interactome during blood development studied with the tool PPICompare

BACKGROUND

Differential analysis of cellular conditions is a key approach towards understanding the consequences and driving causes behind biological processes such as developmental transitions or diseases. The progress of whole-genome expression profiling enabled to conveniently capture the state of a cell's transcriptome and to detect the characteristic features that distinguish cells in specific conditions. In contrast, mapping the physical protein interactome for many samples is experimentally infeasible at the moment. For the understanding of the whole system, however, it is equally important how the interactions of proteins are rewired between cellular states. To overcome this deficiency, we recently showed how condition-specific protein interaction networks that even consider alternative splicing can be inferred from transcript expression data. Here, we present the differential network analysis tool PPICompare that was specifically designed for isoform-sensitive protein interaction networks.

RESULTS

Besides detecting significant rewiring events between the interactomes of grouped samples, PPICompare infers which alterations to the transcriptome caused each rewiring event and what is the minimal set of alterations necessary to explain all between-group changes. When applied to the development of blood cells, we verified that a reasonable amount of rewiring events were reported by the tool and found that differential gene expression was the major determinant of cellular adjustments to the interactome. Alternative splicing events were consistently necessary in each developmental step to explain all significant alterations and were especially important for rewiring in the context of transcriptional control.

CONCLUSIONS

Applying PPICompare enabled us to investigate the dynamics of the human protein interactome during developmental transitions. A platform-independent implementation of the tool PPICompare is available at https://sourceforge.net/projects/ppicompare/ .

A Differentiation Transcription Factor Establishes Muscle-Specific Proteostasis in Caenorhabditis elegans

Safeguarding the proteome is central to the health of the cell. In multi-cellular organisms, the composition of the proteome, and by extension, protein-folding requirements, varies between cells. In agreement, chaperone network composition differs between tissues. Here, we ask how chaperone expression is regulated in a cell type-specific manner and whether cellular differentiation affects chaperone expression. Our bioinformatics analyses show that the myogenic transcription factor HLH-1 (MyoD) can bind to the promoters of chaperone genes expressed or required for the folding of muscle proteins. To test this experimentally, we employed HLH-1 myogenic potential to genetically modulate cellular differentiation of Caenorhabditis elegans embryonic cells by ectopically expressing HLH-1 in all cells of the embryo and monitoring chaperone expression. We found that HLH-1-dependent myogenic conversion specifically induced the expression of putative HLH-1-regulated chaperones in differentiating muscle cells. Moreover, disrupting the putative HLH-1-binding sites on ubiquitously expressed daf-21(Hsp90) and muscle-enriched hsp-12.2(sHsp) promoters abolished their myogenic-dependent expression. Disrupting HLH-1 function in muscle cells reduced the expression of putative HLH-1-regulated chaperones and compromised muscle proteostasis during and after embryogenesis. In turn, we found that modulating the expression of muscle chaperones disrupted the folding and assembly of muscle proteins and thus, myogenesis. Moreover, muscle-specific over-expression of the DNAJB6 homolog DNJ-24, a limb-girdle muscular dystrophy-associated chaperone, disrupted the muscle chaperone network and exposed synthetic motility defects. We propose that cellular differentiation could establish a proteostasis network dedicated to the folding and maintenance of the muscle proteome. Such cell-specific proteostasis networks can explain the selective vulnerability that many diseases of protein misfolding exhibit even when the misfolded protein is ubiquitously expressed.

Molecular chaperones protect proteins from misfolding and aggregation. In multi-cellular organisms, the composition and expression levels of chaperones vary between tissues. However, little is known of how such differential expression is regulated. We hypothesized that the cellular differentiation that regulates the cell-type specific expression program may be involved in establishing a cell-type specific chaperone network. To test this possibility, we addressed the myogenic commitment transcription factor HLH-1 (CeMyoD) that converts embryonic cells to muscle cells in Caenorhabditis elegans. We demonstrated that HLH-1 regulates the expression of muscle chaperones during muscle differentiation. Moreover, we showed that HLH-1-dependent expression of chaperones is required for embryonic development and muscle function. We propose that cellular differentiation results in cell-specific differences in the chaperone network that may be detrimental in terms of the susceptibility of neurons and muscle cells to protein misfolding diseases.

The TissueNet v.2 database: A quantitative view of protein-protein interactions across human tissues

Knowledge of the molecular interactions of human proteins within tissues is important for identifying their tissue-specific roles and for shedding light on tissue phenotypes. However, many protein–protein interactions (PPIs) have no tissue-contexts. The TissueNet database bridges this gap by associating experimentally-identified PPIs with human tissues that were shown to express both pair-mates. Users can select a protein and a tissue, and obtain a network view of the query protein and its tissue-associated PPIs. TissueNet v.2 is an updated version of the TissueNet database previously featured in NAR. It includes over 40 human tissues profiled via RNA-sequencing or protein-based assays. Users can select their preferred expression data source and interactively set the expression threshold for determining tissue-association. The output of TissueNet v.2 emphasizes qualitative and quantitative features of query proteins and their PPIs. The tissue-specificity view highlights tissue-specific and globally-expressed proteins, and the quantitative view highlights proteins that were differentially expressed in the selected tissue relative to all other tissues. Together, these views allow users to quickly assess the unique versus global functionality of query proteins. Thus, TissueNet v.2 offers an extensive, quantitative and user-friendly interface to study the roles of human proteins across tissues. TissueNet v.2 is available at http://netbio.bgu.ac.il/tissuenet.

HIPPIE v2.0: enhancing meaningfulness and reliability of protein-protein interaction networks

The increasing number of experimentally detected interactions between proteins makes it difficult for researchers to extract the interactions relevant for specific biological processes or diseases. This makes it necessary to accompany the large-scale detection of protein–protein interactions (PPIs) with strategies and tools to generate meaningful PPI subnetworks. To this end, we generated the Human Integrated Protein–Protein Interaction rEference or HIPPIE (http://cbdm.uni-mainz.de/hippie/). HIPPIE is a one-stop resource for the generation and interpretation of PPI networks relevant to a specific research question. We provide means to generate highly reliable, context-specific PPI networks and to make sense out of them. We just released the second major update of HIPPIE, implementing various new features. HIPPIE grew substantially over the last years and now contains more than 270 000 confidence scored and annotated PPIs. We integrated different types of experimental information for the confidence scoring and the construction of context-specific networks. We implemented basic graph algorithms that highlight important proteins and interactions. HIPPIE's graphical interface implements several ways for wet lab and computational scientists alike to access the PPI data.

Interactome-transcriptome analysis discovers signatures complementary to GWAS Loci of Type 2 Diabetes

Protein interactions play significant roles in complex diseases. We analyzed peripheral blood mononuclear cells (PBMC) transcriptome using a multi-method strategy. We constructed a tissue-specific interactome (T2Di) and identified 420 molecular signatures associated with T2D-related comorbidity and symptoms, mainly implicated in inflammation, adipogenesis, protein phosphorylation and hormonal secretion. Apart from explaining the residual associations within the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) study, the T2Di signatures were enriched in pathogenic cell type-specific regulatory elements related to fetal development, immunity and expression quantitative trait loci (eQTL). The T2Di revealed a novel locus near a well-established GWAS loci AChE, in which SRRT interacts with JAZF1, a T2D-GWAS gene implicated in pancreatic function. The T2Di also included known anti-diabetic drug targets (e.g. PPARD, MAOB) and identified possible druggable targets (e.g. NCOR2, PDGFR). These T2Di signatures were validated by an independent computational method, and by expression data of pancreatic islet, muscle and liver with some of the signatures (CEBPB, SREBF1, MLST8, SRF, SRRT and SLC12A9) confirmed in PBMC from an independent cohort of 66 T2D and 66 control subjects. By combining prior knowledge and transcriptome analysis, we have constructed an interactome to explain the multi-layered regulatory pathways in T2D.

Tissue Specificity of Human Disease Module

Genes carrying mutations associated with genetic diseases are present in all human cells; yet, clinical manifestations of genetic diseases are usually highly tissue-specific. Although some disease genes are expressed only in selected tissues, the expression patterns of disease genes alone cannot explain the observed tissue specificity of human diseases. Here we hypothesize that for a disease to manifest itself in a particular tissue, a whole functional subnetwork of genes (disease module) needs to be expressed in that tissue. Driven by this hypothesis, we conducted a systematic study of the expression patterns of disease genes within the human interactome. We find that genes expressed in a specific tissue tend to be localized in the same neighborhood of the interactome. By contrast, genes expressed in different tissues are segregated in distinct network neighborhoods. Most important, we show that it is the integrity and the completeness of the expression of the disease module that determines disease manifestation in selected tissues. This approach allows us to construct a disease-tissue network that confirms known and predicts unexpected disease-tissue associations.

Comparative analysis of housekeeping and tissue-specific driver nodes in human protein interaction networks

Background

Several recent studies have used the Minimum Dominating Set (MDS) model to identify driver nodes, which provide the control of the underlying networks, in protein interaction networks. There may exist multiple MDS configurations in a given network, thus it is difficult to determine which one represents the real set of driver nodes. Because these previous studies only focus on static networks and ignore the contextual information on particular tissues, their findings could be insufficient or even be misleading.

Results

In this study, we develop a Collective-Influence-corrected Minimum Dominating Set (CI-MDS) model which takes into account the collective influence of proteins. By integrating molecular expression profiles and static protein interactions, 16 tissue-specific networks are established as well. We then apply the CI-MDS model to each tissue-specific network to detect MDS proteins. It generates almost the same MDSs when it is solved using different optimization algorithms. In addition, we classify MDS proteins into Tissue-Specific MDS (TS-MDS) proteins and HouseKeeping MDS (HK-MDS) proteins based on the number of tissues in which they are expressed and identified as MDS proteins. Notably, we find that TS-MDS proteins and HK-MDS proteins have significantly different topological and functional properties. HK-MDS proteins are more central in protein interaction networks, associated with more functions, evolving more slowly and subjected to a greater number of post-translational modifications than TS-MDS proteins. Unlike TS-MDS proteins, HK-MDS proteins significantly correspond to essential genes, ageing genes, virus-targeted proteins, transcription factors and protein kinases. Moreover, we find that besides HK-MDS proteins, many TS-MDS proteins are also linked to disease related genes, suggesting the tissue specificity of human diseases. Furthermore, functional enrichment analysis reveals that HK-MDS proteins carry out universally necessary biological processes and TS-MDS proteins usually involve in tissue-dependent functions.

Conclusions

Our study uncovers key features of TS-MDS proteins and HK-MDS proteins, and is a step forward towards a better understanding of the controllability of human interactomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1233-0) contains supplementary material, which is available to authorized users.

Integrative analysis of human omics data using biomolecular networks

High-throughput '-omics' technologies have given rise to an increasing abundance of genome-scale data detailing human biology at the molecular level. Although these datasets have already made substantial contributions to a more comprehensive understanding of human physiology and diseases, their interpretation becomes increasingly cryptic and nontrivial as they continue to expand in size and complexity. Systems biology networks offer a scaffold upon which omics data can be integrated, facilitating the extraction of new and physiologically relevant information from the data. Two of the most prevalent networks that have been used for such integrative analyses of omics data are genome-scale metabolic models (GEMs) and protein-protein interaction (PPI) networks, both of which have demonstrated success among many different omics and sample types. This integrative approach seeks to unite 'top-down' omics data with 'bottom-up' biological networks in a synergistic fashion that draws on the strengths of both strategies. As the volume and resolution of high-throughput omics data continue to grow, integrative network-based analyses are expected to play an increasingly important role in their interpretation.

Identification of Tissue-Specific Protein-Coding and Noncoding Transcripts across 14 Human Tissues Using RNA-seq

Many diseases and adverse drug reactions exhibit tissue specificity. To better understand the tissue-specific expression characteristics of transcripts in different human tissues, we deeply sequenced RNA samples from 14 different human tissues. After filtering many lowly expressed transcripts, 24,729 protein-coding transcripts and 1,653 noncoding transcripts were identified. By analyzing highly expressed tissue-specific protein-coding transcripts (TSCTs) and noncoding transcripts (TSNTs), we found that testis expressed the highest numbers of TSCTs and TSNTs. Brain, monocytes, ovary, and heart expressed more TSCTs than the rest tissues, whereas brain, placenta, heart, and monocytes expressed more TSNTs than other tissues. Co-expression network constructed based on the TSCTs and TSNTs showed that each hub TSNT was co-expressed with several TSCTs, allowing functional annotation of TSNTs. Important biological processes and KEGG pathways highly related to the specific functions or diseases of each tissue were enriched with the corresponding TSCTs. These TSCTs and TSNTs may participate in the tissue-specific physiological or pathological processes. Our study provided a unique data set and systematic analysis of expression characteristics and functions of both TSCTs and TSNTs based on 14 distinct human tissues, and could facilitate future investigation of the mechanisms behind tissue-specific diseases and adverse drug reactions.

The analyses of human inherited disease and tissue-specific proteins in the interaction network

OBJECTIVES

With the announcement of human proteome and interaction data, it becomes possible to comprehensively investigate the tissue-expression and network properties of inherited disease proteins. In this study, our goal was to develop methods to map the disease and expression data and analyze the disorder-tissue associations.

METHODS

In this paper, we manually classified the human disease proteins into 22 disorder classes and systematically analyzed the properties of disease proteins in different disorder classes. Then, we investigated the similarity of different disorder classes by computing the overlap of different disorder proteins and networks. We proposed two novel measures, Enrichment Ratio and P-value for comparative analysis of disease proteins across tissues and revealed the associations between disorder classes and tissues/cells.

RESULTS

Compared with non-disease proteins, disease proteins tend to express in more tissues, have higher expression levels and interact with more other proteins in the network. The overlap percentages of networks are much higher than those of proteins, implying that different disorder classes usually influence each other by means of their interacting neighbors. The metabolic, muscular and hematologic proteins are related with most tissues/cells, and cancer proteins are closely associated with the disorders in immune cells.

CONCLUSION

This paper provided novel methods to investigate proteome-wide disease proteins and their interacting networks in order to understand different disease's associations.

Cell type-specific properties and environment shape tissue specificity of cancer genes

One of the biggest mysteries in cancer research remains why mutations in certain genes cause cancer only at specific sites in the human body. The poor correlation between the expression level of a cancer gene and the tissues in which it causes malignant transformations raises the question of which factors determine the tissue-specific effects of a mutation. Here, we explore why some cancer genes are associated only with few different cancer types (i.e., are specific), while others are found mutated in a large number of different types of cancer (i.e., are general). We do so by contrasting cellular functions of specific-cancer genes with those of general ones to identify properties that determine where in the body a gene mutation is causing malignant transformations. We identified different groups of cancer genes that did not behave as expected (i.e., DNA repair genes being tissue specific, immune response genes showing a bimodal specificity function or strong association of generally expressed genes to particular cancers). Analysis of these three groups demonstrates the importance of environmental impact for understanding why certain cancer genes are only involved in the development of some cancer types but are rarely found mutated in other types of cancer.

PPIXpress: construction of condition-specific protein interaction networks based on transcript expression

Protein-protein interaction networks are an important component of modern systems biology. Yet, comparatively few efforts have been made to tailor their topology to the actual cellular condition being studied. Here, we present a network construction method that exploits expression data at the transcript-level and thus reveals alterations in protein connectivity not only caused by differential gene expression but also by alternative splicing. We achieved this by establishing a direct correspondence between individual protein interactions and underlying domain interactions in a complete but condition-unspecific protein interaction network. This knowledge was then used to infer the condition-specific presence of interactions from the dominant protein isoforms. When we compared contextualized interaction networks of matched normal and tumor samples in breast cancer, our transcript-based construction identified more significant alterations that affected proteins associated with cancerogenesis than a method that only uses gene expression data. The approach is provided as the user-friendly tool PPIXpress.

AVAILABILITY AND IMPLEMENTATION

PPIXpress is available at https://sourceforge.net/projects/ppixpress/.

Exploring the cellular basis of human disease through a large-scale mapping of deleterious genes to cell types

BACKGROUND

Each cell type found within the human body performs a diverse and unique set of functions, the disruption of which can lead to disease. However, there currently exists no systematic mapping between cell types and the diseases they can cause.

METHODS

In this study, we integrate protein-protein interaction data with high-quality cell-type-specific gene expression data from the FANTOM5 project to build the largest collection of cell-type-specific interactomes created to date. We develop a novel method, called gene set compactness (GSC), that contrasts the relative positions of disease-associated genes across 73 cell-type-specific interactomes to map genes associated with 196 diseases to the cell types they affect. We conduct text-mining of the PubMed database to produce an independent resource of disease-associated cell types, which we use to validate our method.

RESULTS

The GSC method successfully identifies known disease-cell-type associations, as well as highlighting associations that warrant further study. This includes mast cells and multiple sclerosis, a cell population currently being targeted in a multiple sclerosis phase 2 clinical trial. Furthermore, we build a cell-type-based diseasome using the cell types identified as manifesting each disease, offering insight into diseases linked through etiology.

CONCLUSIONS

The data set produced in this study represents the first large-scale mapping of diseases to the cell types in which they are manifested and will therefore be useful in the study of disease systems. Overall, we demonstrate that our approach links disease-associated genes to the phenotypes they produce, a key goal within systems medicine.

The mRNA related ceRNA-ceRNA landscape and significance across 20 major cancer types

Cross-talk between competitive endogenous RNAs (ceRNAs) through shared miRNAs represents a novel layer of gene regulation that plays important roles in the physiology and development of cancers. However, a global view of their system-level properties across various types of cancers is still unknown. Here, we constructed the mRNA related ceRNA–ceRNA interaction landscape across 20 cancer types by systematically analyzing molecular profiles of 5203 tumors and miRNA regulations. Our study highlights the conserved features shared by pan-cancer and higher similarity within similar origin cell type. Moreover, a core ceRNA network was identified. Function analysis identified a common theme of cancer hallmarks, however they exhibit phenotype-specific connectivity patterns. Besides, we found a marked rewiring in the ceRNA program between various cancers, and further revealed conserved and rewired network ceRNA hubs in each cancer, which were tensely competitive interactions to constitute conserved and cancer-specific modules. By providing mechanistic linkage between known cancer miRNAs, their mediated ceRNA–ceRNA interactions, and the associations with known cancer hallmarks, the inferred cancer ceRNA–ceRNA interaction landscape will serve as a powerful public resource for further biological discoveries of tumorigenesis.

Human protein interaction networks across tissues and diseases

Protein interaction networks are an important framework for studying protein function, cellular processes, and genotype-to-phenotype relationships. While our view of the human interaction network is constantly expanding, less is known about networks that form in biologically important contexts such as within distinct tissues or in disease conditions. Here we review efforts to characterize these networks and to harness them to gain insights into the molecular mechanisms underlying human disease.