Search

The distinct functions and phenotypes of human tissues and cells derive from the activity of biological processes that varies in a context-dependent manner. Here, we present the Process Activity (ProAct) webserver that estimates the preferential activity of biological processes in tissues, cells, and other contexts. Users can upload a differential gene expression matrix measured across contexts or cells, or use a built-in matrix of differential gene expression in 34 human tissues. Per context, ProAct associates gene ontology (GO) biological processes with estimated preferential activity scores, which are inferred from the input matrix. ProAct visualizes these scores across processes, contexts, and process-associated genes. ProAct also offers potential cell-type annotations for cell subsets, by inferring them from the preferential activity of 2001 cell-type-specific processes. Thus, ProAct output can highlight the distinct functions of tissues and cell types in various contexts, and can enhance cell-type annotation efforts. The ProAct webserver is available at https://netbio.bgu.ac.il/ProAct/.

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.

Hereditary diseases manifest clinically in certain tissues, however 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 1,113 hereditary diseases. Overall, we identified 110 cell types affected by 714 diseases. We corroborated our findings by literature text-mining and recapitulation in mouse corresponding tissues. Based on these findings, we explored features of disease-affected cell types and cell classes, highlighted cell types affected by mitochondrial diseases and heritable cancers, and identified diseases that perturb intercellular communication. This study expands our understanding of disease mechanisms and cellular vulnerability.

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.

Hundreds of heritable traits and diseases that are caused by germline aberrations in ubiquitously expressed genes manifest in a remarkably limited number of cell types and tissues across the body. Unravelling mechanisms that govern their tissue-specific manifestations is critical for our understanding of disease aetiologies and may direct efforts to develop treatments. Owing to recent advances in high-throughput technologies and open resources, data and tools are now available to approach this enigmatic phenomenon at large scales, both computationally and experimentally. Here, we discuss the large prevalence of tissue-selective traits and diseases, describe common molecular mechanisms underlying their tissue-selective manifestation and present computational strategies and publicly available resources for elucidating the molecular basis of their genotype–phenotype relationships. The pathology of heritable human traits and diseases often affects a narrow range of tissues, even when causal genes are expressed widely across the body. In this Review, Hekselman and Yeger-Lotem discuss the latest understanding of tissue specificity in human traits and disease, including the diverse underlying molecular mechanisms, experimental and bioinformatics resources to leverage omics data, and implications for understanding disease aetiology.

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.

Tissue contexts are extremely valuable when studying protein functions and their associated phenotypes. Recently, the study of proteins in tissue contexts was greatly facilitated by the availability of thousands of tissue transcriptomes. To provide access to these data we developed the TissueNet integrative database that displays protein-protein interactions (PPIs) in tissue contexts. Through TissueNet, users can create tissue-sensitive network views of the PPI landscape of query proteins. Unlike other tools, TissueNet output networks highlight tissue-specific and broadly expressed proteins, as well as over- and under-expressed proteins per tissue. The TissueNet v.3 upgrade has a much larger dataset of proteins and PPIs, and represents 125 adult tissues and seven embryonic tissues. Thus, TissueNet provides an extensive, quantitative, and user-friendly interface to study the roles of human proteins in adulthood and embryonic stages. TissueNet v3 is freely available at https://netbio.bgu.ac.il/tissuenet3.

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.

Significance Bioinformatic analysis revealed that approximately 40% of human messenger RNAs contain upstream open reading frames (uORFs) in their 5′ untranslated regions. Some of these sequences are translated, but the function of the encoded peptides remains unknown. Our study revealed a uORF encoding for a peptide exhibiting kinase inhibitory activity. This uORF, upstream of a PKC family member, possess the typical pseudosubstrate motif, which autoinhibits the catalytic activity of all PKCs. Using mouse models and human cells, we show that this peptide inhibits cancer cell survival, tumor progression, invasion, and metastasis and synergizes with chemotherapy by interfering with DNA damage response. Together, we point to a previously unrecognized function of a uORF-encoded peptide as a kinase inhibitor, pertinent to cancer therapy., Approximately 40% of human messenger RNAs (mRNAs) contain upstream open reading frames (uORFs) in their 5′ untranslated regions. Some of these uORF sequences, thought to attenuate scanning ribosomes or lead to mRNA degradation, were recently shown to be translated, although the function of the encoded peptides remains unknown. Here, we show a uORF-encoded peptide that exhibits kinase inhibitory functions. This uORF, upstream of the protein kinase C-eta (PKC-η) main ORF, encodes a peptide (uPEP2) containing the typical PKC pseudosubstrate motif present in all PKCs that autoinhibits their kinase activity. We show that uPEP2 directly binds to and selectively inhibits the catalytic activity of novel PKCs but not of classical or atypical PKCs. The endogenous deletion of uORF2 or its overexpression in MCF-7 cells revealed that the endogenously translated uPEP2 reduces the protein levels of PKC-η and other novel PKCs and restricts cell proliferation. Functionally, treatment of breast cancer cells with uPEP2 diminished cell survival and their migration and synergized with chemotherapy by interfering with the response to DNA damage. Furthermore, in a xenograft of MDA-MB-231 breast cancer tumor in mice models, uPEP2 suppressed tumor progression, invasion, and metastasis. Tumor histology showed reduced proliferation, enhanced cell death, and lower protein expression levels of novel PKCs along with diminished phosphorylation of PKC substrates. Hence, our study demonstrates that uORFs may encode biologically active peptides beyond their role as translation regulators of their downstream ORFs. Together, we point to a unique function of a uORF-encoded peptide as a kinase inhibitor, pertinent to cancer therapy.

Motivation The effectiveness of drugs tends to vary between patients. One of the well-known reasons for this phenomenon is genetic polymorphisms in drug target genes among patients. Here, we propose that differences in expression levels of drug target genes across individuals can also contribute to this phenomenon. Results To explore this hypothesis, we analyzed the expression variability of protein-coding genes, and particularly drug target genes, across individuals. For this, we developed a novel variability measure, termed local coefficient of variation (LCV), which ranks the expression variability of each gene relative to genes with similar expression levels. Unlike commonly used methods, LCV neutralizes expression levels biases without imposing any distribution over the variation and is robust to data incompleteness. Application of LCV to RNA-sequencing profiles of 19 human tissues and to target genes of 1076 approved drugs revealed that drug target genes were significantly more variable than protein-coding genes. Analysis of 113 drugs with available effectiveness scores showed that drugs targeting highly variable genes tended to be less effective in the population. Furthermore, comparison of approved drugs to drugs that were withdrawn from the market showed that withdrawn drugs targeted significantly more variable genes than approved drugs. Last, upon analyzing gender differences we found that the variability of drug target genes was similar between men and women. Altogether, our results suggest that expression variability of drug target genes could contribute to the variable responsiveness and effectiveness of drugs, and is worth considering during drug treatment and development. Availability and implementation LCV is available as a python script in GitHub (https://github.com/eyalsim/LCV). Supplementary information Supplementary data are available at Bioinformatics online.

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., Tissue-specific differences in protein folding capacities are poorly understood. Here, the authors show that the human chaperone system consists of ubiquitous core chaperones and tissue-specific variable chaperones, perturbation of which leads to tissue-specific phenotypes.

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 role of different cell types and their interactions in Alzheimer’s disease (AD) is an open question. Here we pursued it by assembling a high-resolution cellular map of the aging frontal cortex by single nucleus RNA-seq of 24 individuals with different clinicopathologic characteristics. We used the map to infer the neocortical cellular architecture of 638 individuals profiled by bulk RNA-seq, providing the sample size necessary for identifying statistically robust associations. We uncovered diverse cell populations associated with AD, including inhibitory neuronal subtypes and oligodendroglial states. We further recovered a network of multicellular communities, each composed of coordinated subpopulations of neuronal, glial and endothelial cells, and found that two of these communities are altered in AD. Finally, we used mediation analyses to prioritize cellular changes that might contribute to cognitive decline. Thus, our deconstruction of the aging neocortex provides a roadmap for evaluating the cellular microenvironments underlying AD and dementia.

The Genotype-Tissue Expression (GTEx) project dissects how genetic variation affects gene expression and splicing. Some human genetic variants affect the amount of RNA produced and the splicing of gene transcripts, crucial steps in development and maintaining a healthy individual. However, some of these changes only occur in a small number of tissues within the body. The Genotype-Tissue Expression (GTEx) project has been expanded over time, and, looking at the final data in version 8, Aguet et al. present a deep characterization of genetic associations and gene expression and splicing in 838 individuals over 49 tissues (see the Perspective by Wilson). This large study was able to characterize the details underlying many aspects of gene expression and provides a resource with which to better understand the fundamental molecular mechanisms of how genetic variants affect gene regulation and complex traits in humans. Science, this issue p. 1318; see also p. 1298 The Genotype-Tissue Expression (GTEx) project was established to characterize genetic effects on the transcriptome across human tissues and to link these regulatory mechanisms to trait and disease associations. Here, we present analyses of the version 8 data, examining 15,201 RNA-sequencing samples from 49 tissues of 838 postmortem donors. We comprehensively characterize genetic associations for gene expression and splicing in cis and trans, showing that regulatory associations are found for almost all genes, and describe the underlying molecular mechanisms and their contribution to allelic heterogeneity and pleiotropy of complex traits. Leveraging the large diversity of tissues, we provide insights into the tissue specificity of genetic effects and show that cell type composition is a key factor in understanding gene regulatory mechanisms in human tissues. We thank the donors and their families for their generous gifts of organ donation for transplantation and tissue donations for the GTEx research project; the Genomics Platform at the Broad Institute for data generation; J. Struewing for support and leadership of the GTEx project; M. Khan and C. Stolte for the illustrations in Fig. 1; and R. Do, D. Jordan, and M. Verbanck for providing GWAS pleiotropy scores. Funding: This work was supported by the Common Fund of the Office of the Director, U.S. National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, NIA, NIAID, and NINDS through NIH contracts HHSN261200800001E (Leidos Prime contract with NCI: A.M.S., D.E.T., N.V.R., J.A.M., L.S., M.E.B., L.Q., T.K., D.B., K.R., and A.U.), 10XS170 (NDRI: W.F.L., J.A.T., G.K., A.M., S.S., R.H., G.Wa., M.J., M.Wa., L.E.B., C.J., J.W., B.R., M.Hu., K.M., L.A.S., H.M.G., M.Mo., and L.K.B.), 10XS171 (Roswell Park Cancer Institute: B.A.F., M.T.M., E.K., B.M.G., K.D.R., and J.B.), 10X172 (Science Care Inc.), 12ST1039 (IDOX), 10ST1035 (Van Andel Institute: S.D.J., D.C.R., and D.R.V.), HHSN268201000029C (Broad Institute: F.A., G.G., K.G.A., A.V.S., X.Li., E.T., S.G., A.G., S.A., K.H.H., D.T.N., K.H., S.R.M., and J.L.N.), 5U41HG009494 (F.A., G.G., and K.G.A.), and through NIH grants R01 DA006227-17 (University of Miami Brain Bank: D.C.M. and D.A.D.), Supplement to University of Miami grant DA006227 (D.C.M. and D.A.D.), R01 MH090941 (University of Geneva), R01 MH090951 and R01 MH090937 (University of Chicago), R01 MH090936 (University of North Carolina–Chapel Hill), R01MH101814 (M.M.-A., V.W., S.B.M., R.G., E.T.D., D.G.-M., and A.V.), U01HG007593 (S.B.M.), R01MH101822 (C.D.B.), U01HG007598 (M.O. and B.E.S.), U01MH104393 (A.P.F.), extension H002371 to 5U41HG002371 (W.J.K.), as well as other funding sources: R01MH106842 (T.L., P.M., E.F., and P.J.H.), R01HL142028 (T.L., Si.Ka., and P.J.H.), R01GM122924 (T.L. and S.E.C.), R01MH107666 (H.K.I.), P30DK020595 (H.K.I.), UM1HG008901 (T.L.), R01GM124486 (T.L.), R01HG010067 (Y.Pa.), R01HG002585 (G.Wa. and M.St.), Gordon and Betty Moore Foundation GBMF 4559 (G.Wa. and M.St.), 1K99HG009916-01 (S.E.C.), R01HG006855 (Se.Ka. and R.E.H.), BIO2015-70777-P, Ministerio de Economia y Competitividad and FEDER funds (M.M.-A., V.W., R.G., and D.G.-M.), la Caixa Foundation ID 100010434 under agreement LCF/BQ/SO15/52260001 (D.G.-M.), NIH CTSA grant UL1TR002550-01 (P.M.), Marie-Skłodowska Curie fellowship H2020 Grant 706636 (S.K.-H.), R35HG010718 (E.R.G.), FPU15/03635, Ministerio de Educación, Cultura y Deporte (M.M.-A.),R01MH109905, 1R01HG010480 (A.Ba.), Searle Scholar Program (A.Ba.), R01HG008150 (S.B.M.), 5T32HG000044-22, NHGRI Institutional Training Grant in Genome Science (N.R.G.), EU IMI program (UE7-DIRECT-115317-1) (E.T.D. and A.V.), FNS funded project RNA1 (31003A_149984) (E.T.D. and A.V.), DK110919 (F.H.), F32HG009987 (F.H.), Massachusetts Lions Eye Research Fund Grant (A.R.H.), Wellcome grant WT108749/Z/15/Z (P.F.), and European Molecular Biology Laboratory (P.F. and D.Z.). Peer Reviewed "Article signat per 1 autors/es del BSC membres del THE GTEX CONSORTIUM: Marta Mele Messeguer"

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. Here, we used human tissue RNA-sequencing profiles to analyze the expression and organization of chaperones across 29 main tissues. We found that relative to protein-coding genes, chaperones were significantly more ubiquitously and highly expressed across all tissues. Nevertheless, differential expression analysis revealed that most chaperones were up- or down-regulated in certain tissues, suggesting that they have tissue-specific roles. In agreement, chaperones that were upregulated in skeletal muscle were highly enriched in mouse myoblasts and in nematode’s muscle tissue, and overlapped significantly with chaperones that are causal for muscle diseases. We also identified a distinct subset of chaperones that formed a uniformly-expressed, cross-family core group conducting basic cellular functions that was significantly more essential for cell survival. Altogether, this suggests a layered architecture of chaperones across tissues that is composed of shared core elements that are complemented by variable elements which give rise to tissue-specific functions and sensitivities, thereby contributing to the tissue-specificity of protein misfolding diseases.Significance StatementProtein misfolding diseases, such as neurodegenerative disorders and myopathies, are often manifested in a specific tissue or even a specific cell type. Enigmatically, however, they are typically caused by mutations in widely expressed proteins. Here we focused on chaperones, the main and basic components of the protein-folding machinery of cells. Computational analyses of large scale tissue transcriptomes unveils that the chaperone system is composed of core essential elements that are uniformly expressed across tissues, and of variable elements that are differentially expressed in a tissue-specific manner. This organization allows each tissue to fit the quality control system to its specific requirements and illuminates the mechanisms that underlie a tissue’s susceptibility to protein-misfolding diseases.

Allele expression (AE) analysis robustly measures cis-regulatory effects. Here, we present and demonstrate the utility of a vast AE resource generated from the GTEx v8 release, containing 15,253 samples spanning 54 human tissues for a total of 431 million measurements of AE at the SNP level and 153 million measurements at the haplotype level. In addition, we develop an extension of our tool phASER that allows effect sizes of cis-regulatory variants to be estimated using haplotype-level AE data. This AE resource is the largest to date, and we are able to make haplotype-level data publicly available. We anticipate that the availability of this resource will enable future studies of regulatory variation across human tissues.

et al., Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype–phenotype relationships1,2. Here we present a human ‘all-by-all’ reference interactome map of human binary protein interactions, or ‘HuRI’. With approximately 53,000 protein–protein interactions, HuRI has approximately four times as many such interactions as there are high-quality curated interactions from small-scale studies. The integration of HuRI with genome3, transcriptome4 and proteome5 data enables cellular function to be studied within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying the specific subcellular roles of protein–protein interactions. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms that might underlie tissue-specific phenotypes of Mendelian diseases. HuRI is a systematic proteome-wide reference that links genomic variation to phenotypic outcomes., This work was primarily supported by the National Institutes of Health (NIH) National Human Genome Research Institute (NHGRI) grant U41HG001715 (M.V., F.P.R., D.E.H., M.A.C., G.D.B. and J.T.) with additional support from NIH grants P50HG004233 (M.V. and F.P.R.), U01HL098166 (M.V.), U01HG007690 (M.V.), R01GM109199 (M.A.C.), Canadian Institute for Health Research (CIHR) Foundation Grants (F.P.R. and J. Rak), the Canada Excellence Research Chairs Program (F.P.R.) and an American Heart Association grant 15CVGPS23430000 (M.V.). D.-K.K. was supported by a Banting Postdoctoral Fellowship through the Natural Sciences and Engineering Research Council (NSERC) of Canada and by the Basic Science Research Program through the National Research Foundation (NRF) of Korea funded by the Ministry of Education (2017R1A6A3A03004385). C. Pons was supported by a Ramon Cajal fellowship (RYC-2017-22959). G.M.S. was supported by NIH Training Grant T32CA009361. M.V. is a Chercheurv Qualifié Honoraire from the Fonds de la Recherche Scientifique (FRS-FNRS, Wallonia-Brussels Federation, Belgium).

The impact of inherited genetic variation on gene expression in humans is well-established. The majority of known expression quantitative trait loci (eQTLs) impact expression of local genes (cis-eQTLs). More research is needed to identify effects of genetic variation on distant genes (trans-eQTLs) and understand their biological mechanisms. One common trans-eQTLs mechanism is “mediation” by a local (cis) transcript. Thus, mediation analysis can be applied to genome-wide SNP and expression data in order to identify transcripts that are “cis-mediators” of trans-eQTLs, including those “cis-hubs” involved in regulation of many trans-genes. Identifying such mediators helps us understand regulatory networks and suggests biological mechanisms underlying trans-eQTLs, both of which are relevant for understanding susceptibility to complex diseases. The multitissue expression data from the Genotype-Tissue Expression (GTEx) program provides a unique opportunity to study cis-mediation across human tissue types. However, the presence of complex hidden confounding effects in biological systems can make mediation analyses challenging and prone to confounding bias, particularly when conducted among diverse samples. To address this problem, we propose a new method: Genomic Mediation analysis with Adaptive Confounding adjustment (GMAC). It enables the search of a very large pool of variables, and adaptively selects potential confounding variables for each mediation test. Analyses of simulated data and GTEx data demonstrate that the adaptive selection of confounders by GMAC improves the power and precision of mediation analysis. Application of GMAC to GTEx data provides new insights into the observed patterns of cis-hubs and trans-eQTL regulation across tissue types.

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.

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

Single-cell analysis of CD4 T cells reveals a population structure associated with inflammation and declined immunity in aging., Age-associated changes in CD4 T-cell functionality have been linked to chronic inflammation and decreased immunity. However, a detailed characterization of CD4 T cell phenotypes that could explain these dysregulated functional properties is lacking. We used single-cell RNA sequencing and multidimensional protein analyses to profile thousands of CD4 T cells obtained from young and old mice. We found that the landscape of CD4 T cell subsets differs markedly between young and old mice, such that three cell subsets—exhausted, cytotoxic, and activated regulatory T cells (aTregs)—appear rarely in young mice but gradually accumulate with age. Most unexpected were the extreme pro- and anti-inflammatory phenotypes of cytotoxic CD4 T cells and aTregs, respectively. These findings provide a comprehensive view of the dynamic reorganization of the CD4 T cell milieu with age and illuminate dominant subsets associated with chronic inflammation and immunity decline, suggesting new therapeutic avenues for age-related diseases.

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.

MotivationDifferential 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

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.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

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.

Determining protein levels in each tissue and how they compare with RNA levels is important for understanding human biology and disease as well as regulatory processes that control protein levels. We quantified the relative protein levels from 12,627 genes across 32 normal human tissue types prepared by the GTEx project. Known and new tissue specific or enriched proteins (5,499) were identified and compared to transcriptome data. Many ubiquitous transcripts are found to encode highly tissue specific proteins. Discordance in the sites of RNA expression and protein detection also revealed potential sites of synthesis and action of protein signaling molecules. Overall, these results provide an extraordinary resource, and demonstrate that understanding protein levels can provide insights into metabolism, regulation, secretome, and human diseases.SummaryQuantitative proteome study of 32 human tissues and integrated analysis with transcriptome data revealed that understanding protein levels could provide in-depth knowledge to post transcriptional or translational regulations, human metabolism, secretome, and diseases.

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.

Gene co-expression networks capture biologically important patterns in gene expression data, enabling functional analyses of genes, discovery of biomarkers, and interpretation of regulatory genetic variants. Most network analyses to date have been limited to assessing correlation between total gene expression levels in a single or small sets of tissues. Here, we have reconstructed networks that capture a much more complete set of regulatory relationships, specifically including regulation of relative isoform abundance and splicing, and tissue-specific connections unique to each of a diverse set of tissues. Using the Genotype-Tissue Expression (GTEx) project v6 RNA-sequencing data across 44 tissues in 449 individuals, we evaluated shared and tissue-specific network relationships. First, we developed a framework called Transcriptome Wide Networks (TWNs) for combining total expression and relative isoform levels into a single sparse network, capturing the complex interplay between the regulation of splicing and transcription. We built TWNs for sixteen tissues, and found that hubs with isoform node neighbors in these networks were strongly enriched for splicing and RNA binding genes, demonstrating their utility in unraveling regulation of splicing in the human transcriptome, and providing a set of candidate shared and tissue-specific regulatory hub genes. Next, we used a Bayesian biclustering model that identifies network edges between genes with co-expression in a single tissue to reconstruct tissue-specific networks (TSNs) for 27 distinct GTEx tissues and for four subsets of related tissues. Using both TWNs and TSNs, we characterized gene co-expression patterns shared across tissues. Finally, we found genetic variants associated with multiple neighboring nodes in our networks, supporting the estimated network structures and identifying 33 genetic variants with distant regulatory impact on transcription and splicing. Our networks provide an improved understanding of the complex relationships between genes in the human transcriptome, including tissue-specificity of gene co-expression, regulation of splicing, and the coordinated impact of genetic variation on transcription.

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.

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., Author summary A longstanding enigma in human genetics is what limits the clinical manifestation of hundreds of hereditary diseases to certain tissues or cell types, while their causal genes are present and expressed throughout the human body. A general conception was that the tissue-wide robustness to the causal aberration is achieved owing to the presence of a compensatory factor, and that disease phenotypes emerge wherever this factor is limited. Here, we tested this general conception at large-scale for the first time. We focused on paralogs of disease-causing genes, which share their functionality and may compensate for their aberration. Based on quantitative analyses of several types of omics data, we show that paralogs of causal genes are down-regulated relative to the disease-causing gene preferentially in the respective disease-manifesting tissue. This tendency is common across various subsets of causal genes, diseases, and tissues. Thus, paralogs of causal genes appear to contribute to the tissue-wide robustness against causal aberrations, and serve as important, yet currently under-appreciated, modifiers of disease manifestation.

Rare genetic variants are abundant in humans and are expected to contribute to individual disease risk1,2,3,4. While genetic association studies have successfully identified common genetic variants associated with susceptibility, these studies are not practical for identifying rare variants1,5. Efforts to distinguish pathogenic variants from benign rare variants have leveraged the genetic code to identify deleterious protein-coding alleles1,6,7, but no analogous code exists for non-coding variants. Therefore, ascertaining which rare variants have phenotypic effects remains a major challenge. Rare non-coding variants have been associated with extreme gene expression in studies using single tissues8,9,10,11, but their effects across tissues are unknown. Here we identify gene expression outliers, or individuals showing extreme expression levels for a particular gene, across 44 human tissues by using combined analyses of whole genomes and multi-tissue RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project v6p release12. We find that 58% of underexpression and 28% of overexpression outliers have nearby conserved rare variants compared to 8% of non-outliers. Additionally, we developed RIVER (RNA-informed variant effect on regulation), a Bayesian statistical model that incorporates expression data to predict a regulatory effect for rare variants with higher accuracy than models using genomic annotations alone. Overall, we demonstrate that rare variants contribute to large gene expression changes across tissues and provide an integrative method for interpretation of rare variants in individual genomes.

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., Author Summary 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.