Spatial Analysis of Functional Enrichment (SAFE) in Large Biological Networks

High-throughput fitness experiments reveal specific vulnerabilities of human-adapted Salmonella during stress and infection

Salmonella enterica is comprised of genetically distinct 'serovars' that together provide an intriguing model for exploring the genetic basis of pathogen evolution. Although the genomes of numerous Salmonella isolates with broad variations in host range and human disease manifestations have been sequenced, the functional links between genetic and phenotypic differences among these serovars remain poorly understood. Here, we conduct high-throughput functional genomics on both generalist (Typhimurium) and human-restricted (Typhi and Paratyphi A) Salmonella at unprecedented scale in the study of this enteric pathogen. Using a comprehensive systems biology approach, we identify gene networks with serovar-specific fitness effects across 25 host-associated stresses encountered at key stages of human infection. By experimentally perturbing these networks, we characterize previously undescribed pseudogenes in human-adapted Salmonella. Overall, this work highlights specific vulnerabilities encoded within human-restricted Salmonella that are linked to the degradation of their genomes, shedding light into the evolution of this enteric pathogen.

Genetic interactions of schizophrenia using gene-based statistical epistasis exclusively identify nervous system-related pathways and key hub genes

Background: The relationship between genotype and phenotype is governed by numerous genetic interactions (GIs), and the mapping of GI networks is of interest for two main reasons: 1) By modelling biological robustness, GIs provide a powerful opportunity to infer compensatory biological mechanisms via the identification of functional relationships between genes, which is of interest for biological discovery and translational research. Biological systems have evolved to compensate for genetic (i.e., variations and mutations) and environmental (i.e., drug efficacy) perturbations by exploiting compensatory relationships between genes, pathways and biological processes; 2) GI facilitates the identification of the direction (alleviating or aggravating interactions) and magnitude of epistatic interactions that influence the phenotypic outcome. The generation of GIs for human diseases is impossible using experimental biology approaches such as systematic deletion analysis. Moreover, the generation of disease-specific GIs has never been undertaken in humans. Methods: We used our Indian schizophrenia case-control (case-816, controls-900) genetic dataset to implement the workflow. Standard GWAS sample quality control procedure was followed. We used the imputed genetic data to increase the SNP coverage to analyse epistatic interactions across the genome comprehensively. Using the odds ratio (OR), we identified the GIs that increase or decrease the risk of a disease phenotype. The SNP-based epistatic results were transformed into gene-based epistatic results. Results: We have developed a novel approach by conducting gene-based statistical epistatic analysis using an Indian schizophrenia case-control genetic dataset and transforming these results to infer GIs that increase the risk of schizophrenia. There were ∼9.5 million GIs with a p-value ≤ 1 × 10-5. Approximately 4.8 million GIs showed an increased risk (OR > 1.0), while ∼4.75 million GIs had a decreased risk (OR <1.0) for schizophrenia. Conclusion: Unlike model organisms, this approach is specifically viable in humans due to the availability of abundant disease-specific genome-wide genotype datasets. The study exclusively identified brain/nervous system-related processes, affirming the findings. This computational approach fills a critical gap by generating practically non-existent heritable disease-specific human GIs from human genetic data. These novel datasets can train innovative deep-learning models, potentially surpassing the limitations of conventional GWAS.

Identification of potential inhibitors of casein kinase 2 alpha of Plasmodium falciparum with potent in vitro activity

Drug resistance to commercially available antimalarials is a major obstacle in malaria control and elimination, creating the need to find new antiparasitic compounds with novel mechanisms of action. The success of kinase inhibitors for oncological treatments has paved the way for the exploitation of protein kinases as drug targets in various diseases, including malaria. Casein kinases are ubiquitous serine/threonine kinases involved in a wide range of cellular processes such as mitotic checkpoint signaling, DNA damage response, and circadian rhythm. In Plasmodium, it is suggested that these protein kinases are essential for both asexual and sexual blood-stage parasites, reinforcing their potential as targets for multi-stage antimalarials. To identify new putative PfCK2α inhibitors, we utilized an in silico chemogenomic strategy involving virtual screening with docking simulations and quantitative structure-activity relationship predictions. Our investigation resulted in the discovery of a new quinazoline molecule (542), which exhibited potent activity against asexual blood stages and a high selectivity index (>100). Subsequently, we conducted chemical-genetic interaction analysis on yeasts with mutations in casein kinases. Our chemical-genetic interaction results are consistent with the hypothesis that 542 inhibits yeast Cka1, which has a hinge region with high similarity to PfCK2α. This finding is in agreement with our in silico results suggesting that 542 inhibits PfCK2α via hinge region interaction.

A proximity-dependent biotinylation map of a human cell

Compartmentalization is a defining characteristic of eukaryotic cells, and partitions distinct biochemical processes into discrete subcellular locations. Microscopy¹ and biochemical fractionation coupled with mass spectrometry^2-4 have defined the proteomes of a variety of different organelles, but many intracellular compartments have remained refractory to such approaches. Proximity-dependent biotinylation techniques such as BioID provide an alternative approach to define the composition of cellular compartments in living cells^5-7. Here we present a BioID-based map of a human cell on the basis of 192 subcellular markers, and define the intracellular locations of 4,145 unique proteins in HEK293 cells. Our localization predictions exceed the specificity of previous approaches, and enabled the discovery of proteins at the interface between the mitochondrial outer membrane and the endoplasmic reticulum that are crucial for mitochondrial homeostasis. On the basis of this dataset, we created humancellmap.org as a community resource that provides online tools for localization analysis of user BioID data, and demonstrate how this resource can be used to understand BioID results better.

τ-SGA: synthetic genetic array analysis for systematically screening and quantifying trigenic interactions in yeast

Systematic complex genetic interaction studies have provided insight into high-order functional redundancies and genetic network wiring of the cell. Here, we describe a method for screening and quantifying trigenic interactions from ordered arrays of yeast strains grown on agar plates as individual colonies. The protocol instructs users on the trigenic synthetic genetic array analysis technique, τ-SGA, for high-throughput screens. The steps describe construction of the double-mutant query strains and the corresponding single-mutant control query strains, which are screened in parallel in two replicates. The screening experimental set-up consists of sequential replica-pinning steps that enable automated mating, meiotic recombination and successive haploid selection steps for the generation of triple mutants, which are scored for colony size as a proxy for fitness, which enables the calculation of trigenic interactions. The procedure described here was used to conduct 422 trigenic interaction screens, which generated ~460,000 yeast triple mutants for trigenic interaction analysis. Users should be familiar with robotic equipment required for high-throughput genetic interaction screens and be proficient at the command line to execute the scoring pipeline. Large-scale screen computational analysis is achieved by using MATLAB pipelines that score raw colony size data to produce τ-SGA interaction scores. Additional recommendations are included for optimizing experimental design and analysis of smaller-scale trigenic interaction screens by using a web-based analysis system, SGAtools. This protocol provides a resource for those who would like to gain a deeper, more practical understanding of trigenic interaction screening and quantification methodology.

Trigenic Synthetic Genetic Array (τ-SGA) Technique for Complex Interaction Analysis

Complex genetic interactions occur when mutant alleles of multiple genes combine to elicit an unexpected phenotype, which could not be predicted given the expectation based on the combination of phenotypes associated with individual mutant alleles. Trigenic Synthetic Genetic Array (τ-SGA) methodology was developed for the systematic analysis of complex interactions involving combinations of three gene perturbations. With a series of replica pinning steps of the τ-SGA procedure, haploid triple mutants are constructed through automated mating and meiotic recombination. For example, a double-mutant query strain carrying two mutant alleles of interest, such as a deletion allele of a nonessential gene and a conditional temperature-sensitive allele of an essential gene, is crossed to an input array of yeast mutants, such as the diagnostic array set of ~1200 mutants, to generate an output array of triple mutants. The colony-size measurements of the resulting triple mutants are used to estimate cellular fitness and quantify trigenic interactions by incorporating corresponding single- and double-mutant fitness estimates. Trigenic interaction networks can be further analyzed for functional modules using various clustering and enrichment analysis tools. Complex genetic interactions are rich in functional information and provide insight into the genotype-to-phenotype relationship, genome size, and speciation.

FMN reduces Amyloid-β toxicity in yeast by regulating redox status and cellular metabolism

Alzheimer's disease (AD) is defined by progressive neurodegeneration, with oligomerization and aggregation of amyloid-β peptides (Aβ) playing a pivotal role in its pathogenesis. In recent years, the yeast Saccharomyces cerevisiae has been successfully used to clarify the roles of different human proteins involved in neurodegeneration. Here, we report a genome-wide synthetic genetic interaction array to identify toxicity modifiers of Aβ42, using yeast as the model organism. We find that FMN1, the gene encoding riboflavin kinase, and its metabolic product flavin mononucleotide (FMN) reduce Aβ42 toxicity. Classic experimental analyses combined with RNAseq show the effects of FMN supplementation to include reducing misfolded protein load, altering cellular metabolism, increasing NADH/(NADH + NAD+) and NADPH/(NADPH + NADP+) ratios and increasing resistance to oxidative stress. Additionally, FMN supplementation modifies Htt103QP toxicity and α-synuclein toxicity in the humanized yeast. Our findings offer insights for reducing cytotoxicity of Aβ42, and potentially other misfolded proteins, via FMN-dependent cellular pathways.

The Proteomic Landscape of Pancreatic Ductal Adenocarcinoma Liver Metastases Identifies Molecular Subtypes and Associations with Clinical Response

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) is a highly metastatic disease that can be separated into distinct subtypes based on molecular signatures. Identifying PDAC subtype-specific therapeutic vulnerabilities is necessary to develop precision medicine approaches to treat PDAC.

EXPERIMENTAL DESIGN

A total of 56 PDAC liver metastases were obtained from the UNMC Rapid Autopsy Program and analyzed with quantitative proteomics. PDAC subtypes were identified by principal component analysis based on protein expression profiling. Proteomic subtypes were further characterized by the associated clinical information, including but not limited to survival analysis, drug treatment response, and smoking and drinking status.

RESULTS

Over 3,960 proteins were identified and used to delineate four distinct PDAC microenvironment subtypes: (i) metabolic; (ii) progenitor-like; (iii) proliferative; and (iv) inflammatory. PDAC risk factors of alcohol and tobacco consumption correlate with subtype classifications. Enhanced survival is observed in FOLFIRINOX treated metabolic and progenitor-like subtypes compared with the proliferative and inflammatory subtypes. In addition, TYMP, PDCD6IP, ERAP1, and STMN showed significant association with patient survival in a subtype-specific manner. Gemcitabine-induced alterations in the proteome identify proteins, such as serine hydroxymethyltransferase 1, associated with drug resistance.

CONCLUSIONS

These data demonstrate that proteomic analysis of clinical PDAC liver metastases can identify molecular signatures unique to disease subtypes and point to opportunities for therapeutic development to improve the treatment of PDAC.

To Embed or Not: Network Embedding as a Paradigm in Computational Biology

Current technology is producing high throughput biomedical data at an ever-growing rate. A common approach to interpreting such data is through network-based analyses. Since biological networks are notoriously complex and hard to decipher, a growing body of work applies graph embedding techniques to simplify, visualize, and facilitate the analysis of the resulting networks. In this review, we survey traditional and new approaches for graph embedding and compare their application to fundamental problems in network biology with using the networks directly. We consider a broad variety of applications including protein network alignment, community detection, and protein function prediction. We find that in all of these domains both types of approaches are of value and their performance depends on the evaluation measures being used and the goal of the project. In particular, network embedding methods outshine direct methods according to some of those measures and are, thus, an essential tool in bioinformatics research.