Associating Genes and Protein Complexes with Disease via Network Propagation

A fundamental challenge in human health is the identification of disease-causing genes. Recently, several studies have tackled this challenge via a network-based approach, motivated by the observation that genes causing the same or similar diseases tend to lie close to one another in a network of protein-protein or functional interactions. However, most of these approaches use only local network information in the inference process and are restricted to inferring single gene associations. Here, we provide a global, network-based method for prioritizing disease genes and inferring protein complex associations, which we call PRINCE. The method is based on formulating constraints on the prioritization function that relate to its smoothness over the network and usage of prior information. We exploit this function to predict not only genes but also protein complex associations with a disease of interest. We test our method on gene-disease association data, evaluating both the prioritization achieved and the protein complexes inferred. We show that our method outperforms extant approaches in both tasks. Using data on 1,369 diseases from the OMIM knowledgebase, our method is able (in a cross validation setting) to rank the true causal gene first for 34% of the diseases, and infer 139 disease-related complexes that are highly coherent in terms of the function, expression and conservation of their member proteins. Importantly, we apply our method to study three multi-factorial diseases for which some causal genes have been found already: prostate cancer, alzheimer and type 2 diabetes mellitus. PRINCE''s predictions for these diseases highly match the known literature, suggesting several novel causal genes and protein complexes for further investigation.     

Ontogeny of the GnRH systems in zebrafish brain: in situ hybridization and promoter-reporter expression analyses in intact animals

The ontogeny of two gonadotropin-releasing-hormone (GnRH) systems, salmon GnRH (sGnRH) and chicken GnRH-II (cGnRH-II), was investigated in zebrafish (Danio rerio). In situ hybridization (ISH) first detected sGnRH mRNA-expressing cells at 1 day post-fertilization (pf) anterior to the developing olfactory organs. Subsequently, cells were seen along the ventral olfactory organs and the olfactory bulbs, reaching the terminal nerve (TN) ganglion at 5–6 days pf. Some cells were detected passing posteriorly through the ventral telencephalon (10–25 days pf), and by 25–30 days pf, sGnRH cells were found in the hypothalamic/preoptic area. Continuous documentation in live zebrafish was achieved by a promoter-reporter expression system. The expression of enhanced green fluorescent protein (EGFP) driven by the sGnRH promoter allowed the earlier detection of cells and projections and the migration of sGnRH neurons. This expression system revealed that long leading processes, presumably axons, preceded the migration of the sGnRH neuron somata. cGnRH-II mRNA expressing cells were initially detected (1 day pf) by ISH analysis at lateral aspects of the midbrain and later on (starting at 5 days pf) at the midline of the midbrain tegmentum. Detection of red fluorescent protein (DsRed) driven by the cGnRH-II promoter confirmed the midbrain expression domain and identified specific hindbrain and forebrain cGnRH-II-cells that were not identified by ISH. The forebrain DsRed-expressing cells seemed to emerge from the same site as the sGnRH-EGFP-expressing cells, as revealed by co-injection of both constructs. These studies indicate that zebrafish TN and hypothalamic sGnRH cell populations share a common embryonic origin and migratory path, and that midbrain cGnRH-II cells originate within the midbrain. This study was supported by the US-Israel Bi-national Agricultural Research and Development (BARD) Foundation (grant 3428-03).     

Gene loss rate: a probabilistic measure for the conservation of eukaryotic genes

The rate of conservation of a gene in evolution is believed to be correlated with its biological importance. Recent studies have devised various conservation measures for genes and have shown that they are correlated with several biological characteristics of functional importance. Specifically, the state-of-the-art propensity for gene loss (PGL) measure was shown to be strongly correlated with gene essentiality and its number of protein–protein interactions (PPIs). The observed correlation between conservation and functional importance varies however between conservation measures, underscoring the need for accurate and general measures for the rate of gene conservation. Here we develop a novel maximum-likelihood approach to computing the rate in which a gene is lost in evolution, motivated by the same principles as those underlying PGL. However, in difference to PGL which considers only the most parsimonious ancestral states of the internal nodes of the phylogenetic tree relating the species, our approach weighs in a probabilistic manner all possible ancestral states, and includes the branch length information as part of the probabilistic model. In application to data of 16 eukaryotic genomes, our approach shows higher correlations with experimental data than PGL, including data on gene lethality, level of connectivity in a PPI network and coherence within functionally related genes.     

Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification

MOTIVATION: Genes are often characterized dichotomously as either housekeeping or single-tissue specific. We conjectured that crucial functional information resides in genes with midrange profiles of expression. RESULTS: To obtain such novel information genome-wide, we have determined the mRNA expression levels for one of the largest hitherto analyzed set of 62 839 probesets in 12 representative normal human tissues. Indeed, when using a newly defined graded tissue specificity index tau, valued between 0 for housekeeping genes and 1 for tissue-specific genes, genes with midrange profiles having 0.15< tau<0.85 were found to constitute >50% of all expression patterns. We developed a binary classification, indicating for every gene the I(B) tissues in which it is overly expressed, and the 12-I(B) tissues in which it shows low expression. The 85 dominant midrange patterns with I(B)=2-11 were found to be bimodally distributed, and to contribute most significantly to the definition of tissue specification dendrograms. Our analyses provide a novel route to infer expression profiles for presumed ancestral nodes in the tissue dendrogram. Such definition has uncovered an unsuspected correlation, whereby de novo enhancement and diminution of gene expression go hand in hand. These findings highlight the importance of gene suppression events, with implications to the course of tissue specification in ontogeny and phylogeny. AVAILABILITY: All data and analyses are publically available at the GeneNote website, http://genecards.weizmann.ac.il/genenote/ and, GEO accession GSE803. CONTACT: doron.lancet@weizmann.ac.il SUPPLEMENTARY INFORMATION: Four tables available at the above site.     

Purification of Na+,K+-ATPase expressed in Pichia pastoris reveals an essential role of phospholipid-protein interactions

Na+,K+-ATPase (porcine alpha/his10-beta) has been expressed in Pichia Pastoris, solubilized in n-dodecyl-beta-maltoside and purified to 70-80% purity by nickel-nitrilotriacetic acid chromatography combined with size exclusion chromatography. The recombinant protein is inactive if the purification is done without added phospholipids. The neutral phospholipid, dioleoylphosphatidylcholine, preserves Na+,K+-ATPase activity of protein prepared in a Na+-containing medium, but activity is lost in a K+-containing medium. By contrast, the acid phospholipid, dioleoylphosphatidylserine, preserves activity in either Na+- or K+-containing media. In optimal conditions activity is preserved for about 2 weeks at 0 degrees C. Both recombinant Na+,K+-ATPase and native pig kidney Na+,K+-ATPase, dissolved in n-dodecyl-beta-maltoside, appear to be mainly stable monomers (alpha/beta) as judged by size exclusion chromatography and sedimentation velocity. Na+,K+-ATPase activities at 37 degrees C of the size exclusion chromatography-purified recombinant and renal Na+,K+-ATPase are comparable but are lower than that of membrane-bound renal Na+,K+-ATPase. The beta subunit is expressed in Pichia Pastoris as two lightly glycosylated polypeptides and is quantitatively deglycosylated by endoglycosidase-H at 0 degrees C, to a single polypeptide. Deglycosylation inactivates Na+,K+-ATPase prepared with dioleoylphosphatidylcholine, whereas dioleoylphosphatidylserine protects after deglycosylation, and Na+,K+-ATPase activity is preserved. This work demonstrates an essential role of phospholipid interactions with Na+,K+-ATPase, including a direct interaction of dioleoylphosphatidylserine, and possibly another interaction of either the neutral or acid phospholipid. Additional lipid effects are likely. A role for the beta subunit in stabilizing conformations of Na+,K+-ATPase (or H+,K+-ATPase) with occluded K+ ions can also be inferred. Purified recombinant Na+,K+-ATPase could become an important experimental tool for various purposes, including, hopefully, structural work.     

The effect of phenotypic plasticity on evolution in multipeaked fitness landscapes

When facing the challenge of developing an individual that best fits its environment, nature demonstrates an interesting combination of two fundamentally different adaptive mechanisms: genetic evolution and phenotypic plasticity. Following numerous computational models, it has become the accepted wisdom that lifetime acclimation (e.g. via learning) smooths the fitness landscape and consequently accelerates evolution. However, analytical studies, focusing on the effect of phenotypic plasticity on evolution in simple unimodal landscapes, have often found that learning hinders the evolutionary process rather than accelerating it. Here, we provide a general framework for studying the effect of plasticity on evolution in multipeaked landscapes and introduce a rigorous mathematical analysis of these dynamics. We show that the convergence rate of the evolutionary process in a given arbitrary one-dimensional fitness landscape is dominated by the largest descent (drawdown) in the landscape and provide numerical evidence to support an analogous dominance also in multidimensional landscapes. We consider several schemes of phenotypic plasticity and examine their effect on the landscape drawdown, identifying the conditions under which phenotypic plasticity is advantageous. The lack of such a drawdown in unimodal landscapes vs. its dominance in multipeaked landscapes accounts for the seemingly contradictory findings of previous studies.     

Gene expression of Caenorhabditis elegans neurons carries information on their synaptic connectivity

The claim that genetic properties of neurons significantly influence their synaptic network structure is a common notion in neuroscience. The nematode Caenorhabditis elegans provides an exciting opportunity to approach this question in a large-scale quantitative manner. Its synaptic connectivity network has been identified, and, combined with cellular studies, we currently have characteristic connectivity and gene expression signatures for most of its neurons. By using two complementary analysis assays we show that the expression signature of a neuron carries significant information about its synaptic connectivity signature, and identify a list of putative genes predicting neural connectivity. The current study rigorously quantifies the relation between gene expression and synaptic connectivity signatures in the C. elegans nervous system and identifies subsets of neurons where this relation is highly marked. The results presented and the genes identified provide a promising starting point for further, more detailed computational and experimental investigations.