Comparative interaction networks: bridging genotype to phenotype

A global genetic interaction network by single-cell imaging and machine learning

Cellular and organismal phenotypes are controlled by complex gene regulatory networks. However, reference maps of gene function are still scarce across different organisms. Here, we generated synthetic genetic interaction and cell morphology profiles of more than 6,800 genes in cultured Drosophila cells. The resulting map of genetic interactions was used for machine learning-based gene function discovery, assigning functions to genes in 47 modules. Furthermore, we devised Cytoclass as a method to dissect genetic interactions for discrete cell states at the single-cell resolution. This approach identified an interaction of Cdk2 and the Cop9 signalosome complex, triggering senescence-associated secretory phenotypes and immunogenic conversion in hemocytic cells. Together, our data constitute a genome-scale resource of functional gene profiles to uncover the mechanisms underlying genetic interactions and their plasticity at the single-cell level.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

Testing the Domino Theory of Gene Loss in Buchnera aphidicola: The Relevance of Epistatic Interactions

The domino theory of gene loss states that when some particular gene loses its function and cripples a cellular function, selection will relax in all functionally related genes, which may allow for the non-functionalization and loss of these genes. Here we study the role of epistasis in determining the pattern of gene losses in a set of genes participating in cell envelope biogenesis in the endosymbiotic bacteria Buchnera aphidicola. We provide statistical evidence indicating pairs of genes in B. aphidicola showing correlated gene loss tend to have orthologs in Escherichia coli known to have alleviating epistasis. In contrast, pairs of genes in B. aphidicola not showing correlated gene loss tend to have orthologs in E. coli known to have aggravating epistasis. These results suggest that during the process of genome reduction in B. aphidicola by gene loss, positive or alleviating epistasis facilitates correlated gene losses while negative or aggravating epistasis impairs correlated gene losses. We interpret this as evidence that the reduced proteome of B. aphidicola contains less pathway redundancy and more compensatory interactions, mimicking the situation of E. coli when grown under environmental constrains.

Golgi-Resident Gαo Promotes Protrusive Membrane Dynamics

To form protrusions like neurites, cells must coordinate their induction and growth. The first requires cytoskeletal rearrangements at the plasma membrane (PM), the second requires directed material delivery from cell's insides. We find that the Gαo-subunit of heterotrimeric G proteins localizes dually to PM and Golgi across phyla and cell types. The PM pool of Gαo induces, and the Golgi pool feeds, the growing protrusions by stimulated trafficking. Golgi-residing KDELR binds and activates monomeric Gαo, atypically for G protein-coupled receptors that normally act on heterotrimeric G proteins. Through multidimensional screenings identifying > 250 Gαo interactors, we pinpoint several basic cellular activities, including vesicular trafficking, as being regulated by Gαo. We further find small Golgi-residing GTPases Rab1 and Rab3 as direct effectors of Gαo. This KDELR → Gαo → Rab1/3 signaling axis is conserved from insects to mammals and controls material delivery from Golgi to PM in various cells and tissues.

Strategic Integration of Multiple Bioinformatics Resources for System Level Analysis of Biological Networks

Recent technological advances in genomics allow the production of biological data at unprecedented tera- and petabyte scales. Efficient mining of these vast and complex datasets for the needs of biomedical research critically depends on a seamless integration of the clinical, genomic, and experimental information with prior knowledge about genotype-phenotype relationships. Such experimental data accumulated in publicly available databases should be accessible to a variety of algorithms and analytical pipelines that drive computational analysis and data mining.We present an integrated computational platform Lynx (Sulakhe et al., Nucleic Acids Res 44:D882-D887, 2016) ( http://lynx.cri.uchicago.edu ), a web-based database and knowledge extraction engine. It provides advanced search capabilities and a variety of algorithms for enrichment analysis and network-based gene prioritization. It gives public access to the Lynx integrated knowledge base (LynxKB) and its analytical tools via user-friendly web services and interfaces. The Lynx service-oriented architecture supports annotation and analysis of high-throughput experimental data. Lynx tools assist the user in extracting meaningful knowledge from LynxKB and experimental data, and in the generation of weighted hypotheses regarding the genes and molecular mechanisms contributing to human phenotypes or conditions of interest. The goal of this integrated platform is to support the end-to-end analytical needs of various translational projects.

The BioGRID interaction database: 2015 update

The Biological General Repository for Interaction Datasets (BioGRID: http://thebiogrid.org) is an open access database that houses genetic and protein interactions curated from the primary biomedical literature for all major model organism species and humans. As of September 2014, the BioGRID contains 749 912 interactions as drawn from 43 149 publications that represent 30 model organisms. This interaction count represents a 50% increase compared to our previous 2013 BioGRID update. BioGRID data are freely distributed through partner model organism databases and meta-databases and are directly downloadable in a variety of formats. In addition to general curation of the published literature for the major model species, BioGRID undertakes themed curation projects in areas of particular relevance for biomedical sciences, such as the ubiquitin-proteasome system and various human disease-associated interaction networks. BioGRID curation is coordinated through an Interaction Management System (IMS) that facilitates the compilation interaction records through structured evidence codes, phenotype ontologies, and gene annotation. The BioGRID architecture has been improved in order to support a broader range of interaction and post-translational modification types, to allow the representation of more complex multi-gene/protein interactions, to account for cellular phenotypes through structured ontologies, to expedite curation through semi-automated text-mining approaches, and to enhance curation quality control.

Beyond GWAS in COPD: probing the landscape between gene-set associations, genome-wide associations and protein-protein interaction networks

OBJECTIVES

To use a systems biology approach to integrate genotype and protein-protein interaction (PPI) data to identify disease network modules associated with chronic obstructive pulmonary disease (COPD) and to perform traditional pathway analysis.

METHODS

We utilized a standard gene-set association approach (FORGE) using gene-based association analysis and gene-set definitions from the molecular signatures database (MSigDB). As a discovery step, we analyzed GWAS results from 2 well-characterized COPD cohorts: COPDGene and GenKOLS. We used a third well-characterized COPD case-control cohort for replication: ECLIPSE. Next, we used dmGWAS, a method that integrates GWAS results with PPI, to identify COPD disease modules.

RESULTS

No gene-sets reached experiment-wide significance in either discovery population. We identified a consensus network of 10 genes identified in modules by integrating GWAS results with PPI that replicated in COPDGene, GenKOLS, and ECLIPSE. Members of 4 gene-sets were enriched among these 10 genes: (i) lung adenocarcinoma tumor-sequencing genes, (ii) IL-7 pathway genes, (iii) kidney cell response to arsenic, and (iv) CD4 T-cell responses. Further, several genes have also been associated with pathophysiology relevant to COPD including KCNK3, NEDD4L, and RIN3. In particular, KCNK3 has been associated with pulmonary arterial hypertension, a common complication in advanced COPD.

CONCLUSION

We report a set of new genes that may influence the etiology of COPD that would not have been identified using traditional GWAS and pathway analyses alone.

Predicting the virulence of MRSA from its genome sequence

Microbial virulence is a complex and often multifactorial phenotype, intricately linked to a pathogen’s evolutionary trajectory. Toxicity, the ability to destroy host cell membranes, and adhesion, the ability to adhere to human tissues, are the major virulence factors of many bacterial pathogens, including Staphylococcus aureus. Here, we assayed the toxicity and adhesiveness of 90 MRSA (methicillin resistant S. aureus) isolates and found that while there was remarkably little variation in adhesion, toxicity varied by over an order of magnitude between isolates, suggesting different evolutionary selection pressures acting on these two traits. We performed a genome-wide association study (GWAS) and identified a large number of loci, as well as a putative network of epistatically interacting loci, that significantly associated with toxicity. Despite this apparent complexity in toxicity regulation, a predictive model based on a set of significant single nucleotide polymorphisms (SNPs) and insertion and deletions events (indels) showed a high degree of accuracy in predicting an isolate’s toxicity solely from the genetic signature at these sites. Our results thus highlight the potential of using sequence data to determine clinically relevant parameters and have further implications for understanding the microbial virulence of this opportunistic pathogen.