Gene Essentiality Analyzed by In Vivo Transposon Mutagenesis and Machine Learning in a Stable Haploid Isolate of Candida albicans

This work was supported by European Research Council Advanced Award 340087 (RAPLODAPT) to J.B., the Dahlem Centre of Plant Sciences (DCPS) of the Freie Universität Berlin (R.K.), Israel Science Foundation grant no. 715/18 (R.S.), the Wellcome Trust (grants 086827, 075470, 101873, and 200208) and the MRC Centre for Medical Mycology (N006364/1) (N.A.R.G.). Data availability.All of the code and required dependencies for analysis of the TnSeq data are available at https://github.com/berman-lab/transposon-pipeline. Library insertion sequences are available at NCBI under project PRJNA490565 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490565). Datasets S1 through S9 are available at https://doi.org/10.6084/m9.figshare.c.4251182.Peer reviewedPublisher PD

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

Gene Expression in the Rodent Brain is Associated with Its Regional Connectivity

The putative link between gene expression of brain regions and their neural connectivity patterns is a fundamental question in neuroscience. Here this question is addressed in the first large scale study of a prototypical mammalian rodent brain, using a combination of rat brain regional connectivity data with gene expression of the mouse brain. Remarkably, even though this study uses data from two different rodent species (due to the data limitations), we still find that the connectivity of the majority of brain regions is highly predictable from their gene expression levels–the outgoing (incoming) connectivity is successfully predicted for 73% (56%) of brain regions, with an overall fairly marked accuracy level of 0.79 (0.83). Many genes are found to play a part in predicting both the incoming and outgoing connectivity (241 out of the 500 top selected genes, p-value<1e-5). Reassuringly, the genes previously known from the literature to be involved in axon guidance do carry significant information about regional brain connectivity. Surveying the genes known to be associated with the pathogenesis of several brain disorders, we find that those associated with schizophrenia, autism and attention deficit disorder are the most highly enriched in the connectivity-related genes identified here. Finally, we find that the profile of functional annotation groups that are associated with regional connectivity in the rodent is significantly correlated with the annotation profile of genes previously found to determine neural connectivity in C. elegans (Pearson correlation of 0.24, p<1e-6 for the outgoing connections and 0.27, p<1e-5 for the incoming). Overall, the association between connectivity and gene expression in a specific extant rodent species' brain is likely to be even stronger than found here, given the limitations of current data

Genome-Scale Metabolic Modeling Elucidates the Role of Proliferative Adaptation in Causing the Warburg Effect

The Warburg effect - a classical hallmark of cancer metabolism - is a counter-intuitive phenomenon in which rapidly proliferating cancer cells resort to inefficient ATP production via glycolysis leading to lactate secretion, instead of relying primarily on more efficient energy production through mitochondrial oxidative phosphorylation, as most normal cells do. The causes for the Warburg effect have remained a subject of considerable controversy since its discovery over 80 years ago, with several competing hypotheses. Here, utilizing a genome-scale human metabolic network model accounting for stoichiometric and enzyme solvent capacity considerations, we show that the Warburg effect is a direct consequence of the metabolic adaptation of cancer cells to increase biomass production rate. The analysis is shown to accurately capture a three phase metabolic behavior that is observed experimentally during oncogenic progression, as well as a prominent characteristic of cancer cells involving their preference for glutamine uptake over other amino acids

Decoupling Environment-Dependent and Independent Genetic Robustness across Bacterial Species

The evolutionary origins of genetic robustness are still under debate: it may arise as a consequence of requirements imposed by varying environmental conditions, due to intrinsic factors such as metabolic requirements, or directly due to an adaptive selection in favor of genes that allow a species to endure genetic perturbations. Stratifying the individual effects of each origin requires one to study the pertaining evolutionary forces across many species under diverse conditions. Here we conduct the first large-scale computational study charting the level of robustness of metabolic networks of hundreds of bacterial species across many simulated growth environments. We provide evidence that variations among species in their level of robustness reflect ecological adaptations. We decouple metabolic robustness into two components and quantify the extents of each: the first, environmental-dependent, is responsible for at least 20% of the non-essential reactions and its extent is associated with the species' lifestyle (specialized/generalist); the second, environmental-independent, is associated (correlation = ∼0.6) with the intrinsic metabolic capacities of a species—higher robustness is observed in fast growers or in organisms with an extensive production of secondary metabolites. Finally, we identify reactions that are uniquely susceptible to perturbations in human pathogens, potentially serving as novel drug-targets

Transcriptional Regulation by CHIP/LDB Complexes

It is increasingly clear that transcription factors play versatile roles in turning genes “on” or “off” depending on cellular context via the various transcription complexes they form. This poses a major challenge in unraveling combinatorial transcription complex codes. Here we use the powerful genetics of Drosophila combined with microarray and bioinformatics analyses to tackle this challenge. The nuclear adaptor CHIP/LDB is a major developmental regulator capable of forming tissue-specific transcription complexes with various types of transcription factors and cofactors, making it a valuable model to study the intricacies of gene regulation. To date only few CHIP/LDB complexes target genes have been identified, and possible tissue-dependent crosstalk between these complexes has not been rigorously explored. SSDP proteins protect CHIP/LDB complexes from proteasome dependent degradation and are rate-limiting cofactors for these complexes. By using mutations in SSDP, we identified 189 down-stream targets of CHIP/LDB and show that these genes are enriched for the binding sites of APTEROUS (AP) and PANNIER (PNR), two well studied transcription factors associated with CHIP/LDB complexes. We performed extensive genetic screens and identified target genes that genetically interact with components of CHIP/LDB complexes in directing the development of the wings (28 genes) and thoracic bristles (23 genes). Moreover, by in vivo RNAi silencing we uncovered novel roles for two of the target genes, xbp1 and Gs-alpha, in early development of these structures. Taken together, our results suggest that loss of SSDP disrupts the normal balance between the CHIP-AP and the CHIP-PNR transcription complexes, resulting in down-regulation of CHIP-AP target genes and the concomitant up-regulation of CHIP-PNR target genes. Understanding the combinatorial nature of transcription complexes as presented here is crucial to the study of transcription regulation of gene batteries required for development

PathBLAST: a tool for alignment of protein interaction networks

PathBLAST is a network alignment and search tool for comparing protein interaction networks across species to identify protein pathways and complexes that have been conserved by evolution. The basic method searches for high-scoring alignments between pairs of protein interaction paths, for which proteins of the first path are paired with putative orthologs occurring in the same order in the second path. This technique discriminates between true- and false-positive interactions and allows for functional annotation of protein interaction pathways based on similarity to the network of another, well-characterized species. PathBLAST is now available at http://www.pathblast.org/ as a web-based query. In this implementation, the user specifies a short protein interaction path for query against a target protein–protein interaction network selected from a network database. PathBLAST returns a ranked list of matching paths from the target network along with a graphical view of these paths and the overlap among them. Target protein–protein interaction networks are currently available for Helicobacter pylori, Saccharomyces cerevisiae, Caenorhabditis elegans and Drosophila melanogaster. Just as BLAST enables rapid comparison of protein sequences between genomes, tools such as PathBLAST are enabling comparative genomics at the network level