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. An interpretable machine‐learning framework predicts disease genes for tissue‐selective hereditary diseases. The framework highlights known and novel tissue‐selectivity features and enhances genetic diagnosis of rare‐disease patients. An interpretable machine‐learning framework predicts disease genes for tissue‐selective hereditary diseases. The framework highlights known and novel tissue‐selectivity features and enhances genetic diagnosis of rare‐disease patients.

ABSTRACT Genetic studies of Mendelian and rare diseases face the critical challenges of identifying pathogenic gene variants and their modes-of-action. Previous efforts rarely utilized the tissue-selective manifestation of these diseases for their elucidation. Here we introduce an interpretable machine learning (ML) platform that utilizes heterogeneous and large-scale tissue-aware datasets of human genes, and rigorously, concurrently and quantitatively assesses hundreds of candidate mechanisms per disease. The resulting tissue-aware ML platform is applicable in gene-specific, tissue-specific, or patient-specific modes. Application of the platform to selected Mendelian disease genes pinpointed mechanisms that lead to tissue-specific disease manifestation. When applied jointly to diseases that manifest in the same tissue, the models revealed common known and previously underappreciated factors that underlie tissue-selective disease manifestation. Lastly, we harnessed our ML platform toward genetic diagnosis of tissue-selective rare diseases. Patient-specific models of candidate disease-causing genes from 50 patients successfully prioritized the pathogenic gene in 86% of the cases, implying that the tissue-selectivity of rare diseases aids in filtering out unlikely candidate genes. Thus, interpretable tissue-aware ML models can boost mechanistic understanding and genetic diagnosis of tissue-selective heritable diseases. A webserver supporting gene prioritization is available at https://netbio.bgu.ac.il/trace/ .

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.

Exploring the Phenotypic Consequences of Tissue Specific Gene Expression Variation Inferred GWAS Summary Statistics. Scalable, integrative methods to understand mechanisms that link genetic variants with phenotypes are needed. Here we derive a mathematical expression to compute PrediXcan (a gene mapping approach) results using summary data (S-PrediXcan) and show its accuracy and general robustness to misspeci fi ed reference sets. We apply this framework to 44 GTEx tissues and 100 + phenotypes from GWAS and meta-analysis studies, creating a growing public catalog of associations that seeks to capture the effects of gene expression variation on human phenotypes. Replication in an independent cohort is shown. Most of the associations are tissue speci fi c, suggesting context speci fi city of the trait etiology. Colocalized signi fi cant associations in unexpected tissues underscore the need for an agnostic scanning of multiple contexts to improve our ability to detect causal regulatory mechanisms. Monogenic disease genes are enriched among signi fi cant associations for related traits, suggesting that smaller alterations of these genes may cause a spectrum of milder phenotypes. in concordance with COLOC. A few genes fall in the “ colocalized ” region, in disagreement with COLOC classi fi cation. Unlike COLOC results, HEIDI does not partition the genes into distinct clusters and an arbitrary cutoff p -value has to be chosen. e Shows genes with large HEIDI p -value (no evidence of heterogeneity) which fall in large part in the “ colocalized ” region. However a substantial number fall in “ independent signal ” region, disagreeing with COLOC ’ s classi fi cation signals) from S-TWAS signi fi cant vs. S-PrediXcan signi fi cant results. For most phenotypes, TWAS has lower proportion of colocalized signals compared to S-PrediXcan, as estimated via COLOC. Phenotype abbreviations are as follows: FNBD Femoral Neck Bone Density, LSBD Lumbar Spine Bone Density, BMI Body Mass Index, HEIGHT Height, LDL Low-Density Lipoprotein Cholesterol, HDL High-Density Lipoprotein Cholesterol, TRYG Tryglicerides, CROHN Crohn ’ s Disease, INFBOWEL In fl ammatory Bowel ’ s Disease, ULCERC Ulcerative Colitis, HBA1C Hemogoblin Levels, HOMA-IR HOMA Insulin Response, SCZ Schizophrenia, RA Rheumatoid Arthritis, COLLEGE College Completion, EDUCYEARS Education Years and genes with low some Bone Low-Density Cholesterol, HDL High-Density Lipoprotein Tryglicerides, In Bowel HBA1C Hemogoblin Levels, HOMA-IR HOMA

Global insights into cellular organization and function require comprehensive understanding of interactome networks. Similar to how a reference genome sequence revolutionized human genetics, a reference map of the human interactome network is critical to fully understand genotype-phenotype relationships. Here we present the first human “all-by-all” binary reference interactome map, or “HuRI”. With ~53,000 high-quality protein-protein interactions (PPIs), HuRI is approximately four times larger than the information curated from small-scale studies available in the literature. Integrating HuRI with genome, transcriptome and proteome data enables the study of cellular function within essentially any physiological or pathological cellular context. We demonstrate the use of HuRI in identifying specific subcellular roles of PPIs and protein function modulation via splicing during brain development. 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 thus represents an unprecedented, systematic reference linking genomic variation to phenotypic outcomes.

ABSTRACTMotivationDifferential network analysis, designed to highlight interaction changes between conditions, is an important paradigm in network biology. However, network analysis methods have been typically designed to compare between few conditions, were rarely applied to protein interaction networks (interactomes). Moreover, large-scale benchmarks for their evaluation have been lacking.ResultsHere, we assess five network analysis methods by applying them to 34 human tissues interactomes. For this, we created a manually-curated benchmark of 6,499 tissue-specific, gene ontology biological processes, and analyzed the ability of each method to expose these tissue-process associations. The four differential network analysis methods outperformed the non-differential, expression-based method (AUCs of 0.82-0.9 versus 0.69, respectively). We then created another benchmark, of 1,527 tissue-specific disease cases, and analyzed the ability of differential network analysis methods to highlight additional disease-related genes. Compared to a non-differential subnetworks surrounding a known disease-causing gene, the extremely-differential subnetwork (top 1%) was significantly enriched for additional disease-causing genes in 18.6% of the cases (p≤10e-3). In 5/10 tissues tested, including Muscle, nerve and heart tissues (p = 2.54E-05, 2.71E-04, 3.63E-19), such enrichments were highly significant.SummaryAltogether, our study demonstrates that differential network analysis of human tissue interactomes is a powerful tool for highlighting processes and genes with tissue-selective functionality and clinical impact. Moreover, it offers expansive manually-curated datasets of tissue-selective processes and diseases that could serve for benchmark and for analyses in many other studies.Contactestiyl@bgu.ac.il

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.

RSRC1, whose polymorphism is associated with altered brain function in schizophrenia, is a member of the serine and arginine rich-related protein family. Through homozygosity mapping and whole exome sequencing we show that RSRC1 mutation causes an autosomal recessive syndrome of intellectual disability, aberrant behaviour, hypotonia and mild facial dysmorphism with normal brain MRI. Further, we show that RSRC1 is ubiquitously expressed, and that the RSRC1 mutation triggers nonsense-mediated mRNA decay of the RSRC1 transcript in patients' fibroblasts. Short hairpin RNA (shRNA)-mediated lentiviral silencing and overexpression of RSRC1 in SH-SY5Y cells demonstrated that RSRC1 has a role in alternative splicing and transcription regulation. Transcriptome profiling of RSRC1-silenced cells unravelled specific differentially expressed genes previously associated with intellectual disability, hypotonia and schizophrenia, relevant to the disease phenotype. Protein-protein interaction network modelling suggested possible intermediate interactions by which RSRC1 affects gene-specific differential expression. Patient-derived induced pluripotent stem cells, differentiated into neural progenitor cells, showed expression dynamics similar to the RSRC1-silenced SH-SY5Y model. Notably, patient neural progenitor cells had 9.6-fold downregulated expression of IGFBP3, whose brain expression is affected by MECP2, aberrant in Rett syndrome. Interestingly, Igfbp3-null mice have behavioural impairment, abnormal synaptic function and monoaminergic neurotransmission, likely correlating with the disease phenotype.

The Saccharomyces cerevisiae kinase/adenosine triphosphatase Rio1 regulates rDNA transcription and segregation, pre-rRNA processing and small ribosomal subunit maturation. Other roles are unknown. When overexpressed, human ortholog RIOK1 drives tumor growth and metastasis. Likewise, RIOK1 promotes 40S ribosomal subunit biogenesis and has not been characterized globally. We show that Rio1 manages directly and via a series of regulators, an essential signaling network at the protein, chromatin and RNA levels. Rio1 orchestrates growth and division depending on resource availability, in parallel to the nutrient-activated Tor1 kinase. To define the Rio1 network, we identified its physical interactors, profiled its target genes/transcripts, mapped its chromatin-binding sites and integrated our data with yeast's protein-protein and protein-DNA interaction catalogs using network computation. We experimentally confirmed network components and localized Rio1 also to mitochondria and vacuoles. Via its network, Rio1 commands protein synthesis (ribosomal gene expression, assembly and activity) and turnover (26S proteasome expression), and impinges on metabolic, energy-production and cell-cycle programs. We find that Rio1 activity is conserved to humans and propose that pathological RIOK1 may fuel promiscuous transcription, ribosome production, chromosomal instability, unrestrained metabolism and proliferation; established contributors to cancer. Our study will advance the understanding of numerous processes, here revealed to depend on Rio1 activity.

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.

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.

As challenges to food security increase, the demand for lead genes for improving crop production is growing. However, genetic screens of plant mutants typically yield very low frequencies of desired phenotypes. Here, we present a powerful computational approach for selecting candidate genes for screening insertion mutants. We combined ranking of Arabidopsis thaliana regulatory genes according to their expression in response to multiple abiotic stresses (Multiple Stress [MST] score), with stress-responsive RNA co-expression network analysis to select candidate multiple stress regulatory (MSTR) genes. Screening of 62 T-DNA insertion mutants defective in candidate MSTR genes, for abiotic stress germination phenotypes yielded a remarkable hit rate of up to 62%; this gene discovery rate is 48-fold greater than that of other large-scale insertional mutant screens. Moreover, the MST score of these genes could be used to prioritize them for screening. To evaluate the contribution of the co-expression analysis, we screened 64 additional mutant lines of MST-scored genes that did not appear in the RNA co-expression network. The screening of these MST-scored genes yielded a gene discovery rate of 36%, which is much higher than that of classic mutant screens but not as high as when picking candidate genes from the co-expression network. The MSTR co-expression network that we created, AraSTressRegNet is publicly available at http://netbio.bgu.ac.il/arnet. This systems biology-based screening approach combining gene ranking and network analysis could be generally applicable to enhancing identification of genes regulating additional processes in plants and other organisms provided that suitable transcriptome data are available.

The identification of the molecular pathways active in specific contexts, such as disease states or drug responses, often requires an extensive view of the potential interactions between a subset of proteins. This view is not easily obtained: it requires the integration of context-specific protein list or expression data with up-to-date data of protein interactions that are typically spread across multiple databases. The MyProteinNet web server allows users to easily create such context-sensitive protein interaction networks. Users can automatically gather and consolidate data from up to 11 different databases to create a generic protein interaction network (interactome). They can score the interactions based on reliability and filter them by user-defined contexts including molecular expression and protein annotation. The output of MyProteinNet includes the generic and filtered interactome files, together with a summary of their network attributes. MyProteinNet is particularly geared toward building human tissue interactomes, by maintaining tissue expression profiles from multiple resources. The ability of MyProteinNet to facilitate the construction of up-to-date, context-specific interactomes and its applicability to 11 different organisms and to tens of human tissues, make it a powerful tool in meaningful analysis of protein networks. MyProteinNet is available at http://netbio.bgu.ac.il/myproteinnet.

Genome sequencing and transcriptomic profiling are two widely used approaches for the identification of human disease pathways. However, each approach typically provides a limited view of disease pathways: Genome sequencing can identify disease-related mutations but rarely reveals their mode-of-action, while transcriptomic assays do not reveal the series of events that lead to the transcriptomic change. ResponseNet is an integrative network-optimization approach that we developed to fill these gaps by highlighting major signaling and regulatory molecular interaction paths that connect disease-related mutations and genes. The ResponseNet web-server provides a user-friendly interface to ResponseNet. Specifically, users can upload weighted lists of proteins and genes and obtain a sparse, weighted, molecular interaction subnetwork connecting them, that is biased toward regulatory and signaling pathways. ResponseNet2.0 enhances the functionality of the ResponseNet web-server in two important ways. First, it supports analysis of human data by offering a human interactome composed of proteins, genes and micro-RNAs. Second, it offers a new informative view of the output, including a randomization analysis, to help users assess the biological relevance of the output subnetwork. ResponseNet2.0 is available at http://netbio.bgu.ac.il/respnet .

Knowledge of protein-protein interactions (PPIs) is important for identifying the functions of proteins and the processes they are involved in. Although data of human PPIs are easily accessible through several public databases, these databases do not specify the human tissues in which these PPIs take place. The TissueNet database of human tissue PPIs (http://netbio.bgu.ac.il/tissuenet/) associates each interaction with human tissues that express both pair mates. This was achieved by integrating current data of experimentally detected PPIs with extensive data of gene and protein expression across 16 main human tissues. Users can query TissueNet using a protein and retrieve its PPI partners per tissue, or using a PPI and retrieve the tissues expressing both pair mates. The graphical representation of the output highlights tissue-specific and tissue-wide PPIs. Thus, TissueNet provides a unique platform for assessing the roles of human proteins and their interactions across tissues.