ZapG (YhcB/DUF1043), a novel cell division protein in gamma-proteobacteria linking the Z-ring to septal peptidoglycan synthesis.

YhcB, a poorly understood protein conserved across gamma-proteobacteria, contains a domain of unknown function (DUF1043) and an N-terminal transmembrane domain. Here, we used an integrated approach including X-ray crystallography, genetics, and molecular biology to investigate the function and structure of YhcB. The Escherichia coli yhcB KO strain does not grow at 45 °C and is hypersensitive to cell wall-acting antibiotics, even in the stationary phase. The deletion of yhcB leads to filamentation, abnormal FtsZ ring formation, and aberrant septum development. The Z-ring is essential for the positioning of the septa and the initiation of cell division. We found that YhcB interacts with proteins of the divisome (e.g., FtsI, FtsQ) and elongasome (e.g., RodZ, RodA). Seven of these interactions are also conserved in Yersinia pestis and/or Vibrio cholerae. Furthermore, we mapped the amino acid residues likely involved in the interactions of YhcB with FtsI and RodZ. The 2.8 Å crystal structure of the cytosolic domain of Haemophilus ducreyi YhcB shows a unique tetrameric α-helical coiled-coil structure likely to be involved in linking the Z-ring to the septal peptidoglycan-synthesizing complexes. In summary, YhcB is a conserved and conditionally essential protein that plays a role in cell division and consequently affects envelope biogenesis. Based on these findings, we propose to rename YhcB to ZapG (Z-ring-associated protein G). This study will serve as a starting point for future studies on this protein family and on how cells transit from exponential to stationary survival.

A Gaussian process-based definition reveals new and bona fide genetic interactions compared to a multiplicative model in the Gram-negative Escherichia coli.

MOTIVATION: A digenic genetic interaction (GI) is observed when mutations in two genes within the same organism yield a phenotype that is different from the expected, given each mutation's individual effects. While multiplicative scoring is widely applied to define GIs, revealing underlying gene functions, it remains unclear if it is the most suitable choice for scoring GIs in Escherichia coli. Here, we assess many different definitions, including the multiplicative model, for mapping functional links between genes and pathways in E.coli. RESULTS: Using our published E.coli GI datasets, we show computationally that a machine learning Gaussian process (GP)-based definition better identifies functional associations among genes than a multiplicative model, which we have experimentally confirmed on a set of gene pairs. Overall, the GP definition improves the detection of GIs, biological reasoning of epistatic connectivity, as well as the quality of GI maps in E.coli, and, potentially, other microbes. AVAILABILITY AND IMPLEMENTATION: The source code and parameters used to generate the machine learning models in WEKA software were provided in the Supplementary information. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Legionella pneumophila effector Lem4 is a membrane-associated protein tyrosine phosphatase.

Legionella pneumophila is a Gram-negative pathogenic bacterium that causes severe pneumonia in humans. It establishes a replicative niche called Legionella-containing vacuole (LCV) that allows bacteria to survive and replicate inside pulmonary macrophages. To hijack host cell defense systems, L. pneumophila injects over 300 effector proteins into the host cell cytosol. The Lem4 effector (lpg1101) consists of two domains: an N-terminal haloacid dehalogenase (HAD) domain with unknown function and a C-terminal phosphatidylinositol 4-phosphate-binding domain that anchors Lem4 to the membrane of early LCVs. Herein, we demonstrate that the HAD domain (Lem4-N) is structurally similar to mouse MDP-1 phosphatase and displays phosphotyrosine phosphatase activity. Substrate specificity of Lem4 was probed using a tyrosine phosphatase substrate set, which contained a selection of 360 phosphopeptides derived from human phosphorylation sites. This assay allowed us to identify a consensus pTyr-containing motif. Based on the localization of Lem4 to lysosomes and to some extent to plasma membrane when expressed in human cells, we hypothesize that this protein is involved in protein-protein interactions with an LCV or plasma membrane-associated tyrosine-phosphorylated host target.

Insights from protein-protein interaction studies on bacterial pathogenesis.

INTRODUCTION: The threat bacterial pathogens pose to human health is increasing with the number and distribution of antibiotic-resistant bacteria, while the rate of discovery of new antimicrobials dwindles. Proteomics is playing key roles in understanding the molecular mechanisms of bacterial pathogenesis, and in identifying disease outcome determinants. The physical associations identified by proteomics can provide the means to develop pathogen-specific treatment methods that reduce the spread of antibiotic resistance and alleviate the negative effects of broad-spectrum antibiotics on beneficial bacteria. Areas covered: This review discusses recent trends in proteomics and introduces new and developing approaches that can be applied to the study of protein-protein interactions (PPIs) underlying bacterial pathogenesis. The approaches examined encompass options for mapping proteomes as well as stable and transient interactions in vivo and in vitro. We also explored the coverage of bacterial and human-bacterial PPIs, knowledge gaps in this area, and how they can be filled. Expert commentary: Identifying potential antimicrobial candidates is confounded by the complex molecular biology of bacterial pathogenesis and the lack of knowledge about PPIs underlying this process. Proteomics approaches can offer new perspectives for mechanistic insights and identify essential targets for guiding the discovery of next generation antimicrobials.

Quantitative genome-wide genetic interaction screens reveal global epistatic relationships of protein complexes in Escherichia coli.

Large-scale proteomic analyses in Escherichia coli have documented the composition and physical relationships of multiprotein complexes, but not their functional organization into biological pathways and processes. Conversely, genetic interaction (GI) screens can provide insights into the biological role(s) of individual gene and higher order associations. Combining the information from both approaches should elucidate how complexes and pathways intersect functionally at a systems level. However, such integrative analysis has been hindered due to the lack of relevant GI data. Here we present a systematic, unbiased, and quantitative synthetic genetic array screen in E. coli describing the genetic dependencies and functional cross-talk among over 600,000 digenic mutant combinations. Combining this epistasis information with putative functional modules derived from previous proteomic data and genomic context-based methods revealed unexpected associations, including new components required for the biogenesis of iron-sulphur and ribosome integrity, and the interplay between molecular chaperones and proteases. We find that functionally-linked genes co-conserved among γ-proteobacteria are far more likely to have correlated GI profiles than genes with divergent patterns of evolution. Overall, examining bacterial GIs in the context of protein complexes provides avenues for a deeper mechanistic understanding of core microbial systems.

Novel function discovery with GeneMANIA: a new integrated resource for gene function prediction in Escherichia coli.

MOTIVATION: The model bacterium Escherichia coli is among the best studied prokaryotes, yet nearly half of its proteins are still of unknown biological function. This is despite a wealth of available large-scale physical and genetic interaction data. To address this, we extended the GeneMANIA function prediction web application developed for model eukaryotes to support E.coli. RESULTS: We integrated 48 distinct E.coli functional interaction datasets and used the GeneMANIA algorithm to produce thousands of novel functional predictions and prioritize genes for further functional assays. Our analysis achieved cross-validation performance comparable to that reported for eukaryotic model organisms, and revealed new functions for previously uncharacterized genes in specific bioprocesses, including components required for cell adhesion, iron-sulphur complex assembly and ribosome biogenesis. The GeneMANIA approach for network-based function prediction provides an innovative new tool for probing mechanisms underlying bacterial bioprocesses. CONTACT: gary.bader@utoronto.ca; mohan.babu@uregina.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Investigating Bacterial Protein Synthesis Using Systems Biology Approaches.

Protein synthesis is essential for bacterial growth and survival. Its study in Escherichia coli helped uncover features conserved among bacteria as well as universally. The pattern of discovery and the identification of some of the longest-known components of the protein synthesis machinery, including the ribosome itself, tRNAs, and translation factors proceeded through many stages of successively more refined biochemical purifications, finally culminating in the isolation to homogeneity, identification, and mapping of the smallest unit required for performing the given function. These early studies produced a wealth of information. However, many unknowns remained. Systems biology approaches provide an opportunity to investigate protein synthesis from a global perspective, overcoming the limitations of earlier ad hoc methods to gain unprecedented insights. This chapter reviews innovative systems biology approaches, with an emphasis on those designed specifically for investigating the protein synthesis machinery in E. coli.

Genetic interaction maps in Escherichia coli reveal functional crosstalk among cell envelope biogenesis pathways.

As the interface between a microbe and its environment, the bacterial cell envelope has broad biological and clinical significance. While numerous biosynthesis genes and pathways have been identified and studied in isolation, how these intersect functionally to ensure envelope integrity during adaptive responses to environmental challenge remains unclear. To this end, we performed high-density synthetic genetic screens to generate quantitative functional association maps encompassing virtually the entire cell envelope biosynthetic machinery of Escherichia coli under both auxotrophic (rich medium) and prototrophic (minimal medium) culture conditions. The differential patterns of genetic interactions detected among > 235,000 digenic mutant combinations tested reveal unexpected condition-specific functional crosstalk and genetic backup mechanisms that ensure stress-resistant envelope assembly and maintenance. These networks also provide insights into the global systems connectivity and dynamic functional reorganization of a universal bacterial structure that is both broadly conserved among eubacteria (including pathogens) and an important target.

Investigating the Ustilago maydis/Zea mays pathosystem: transcriptional responses and novel functional aspects of a fungal calcineurin regulatory B subunit.

The sustainable control of basidiomycete biotrophic plant pathogenesis requires an understanding of host responses to infection, as well as the identification and functional analysis of fungal genes involved in disease development. The creation and analysis of a suppressive subtractive hybridization (SSH) cDNA library from Ustilago maydis-infected Zea mays seedlings enabled the identification of fungal and plant genes expressed during disease development, and uncovered new insights into the interactions of this model system. Candidate U. maydis pathogenesis genes were identified by using the current SSH cDNA library analysis, and by knowledge generated from previous cDNA microarray and comparative genomic analyses. These identifications were supported by the independent determination of transcript level changes in different cell-types and during pathogenic development. The basidiomycete specific um01632, the highly in planta expressed um03046 (zig1), and the calcineurin regulatory B subunit (um10226, cnb1), were chosen for deletion experiments. um01632 and zig1 mutants showed no difference in morphology and did not have a statistically significant impact on pathogenesis. cnb1 mutants had a distinct cell division phenotype and reduced virulence in seedling assays. Infections with reciprocal wild-type×Δcnb1 haploid strain crosses revealed that the wild-type allele was unable to fully compensate for the lack of a second cnb1 allele. This haploinsufficiency was undetected in other fungal cnb1 mutational analyses. The reported data improves U. maydis genome annotation and expands on the current understanding of pathogenesis genes in this model basidiomycete.

Co-fractionation-mass spectrometry to characterize native mitochondrial protein assemblies in mammalian neurons and brain.

Human mitochondrial (mt) protein assemblies are vital for neuronal and brain function, and their alteration contributes to many human disorders, e.g., neurodegenerative diseases resulting from abnormal protein-protein interactions (PPIs). Knowledge of the composition of mt protein complexes is, however, still limited. Affinity purification mass spectrometry (MS) and proximity-dependent biotinylation MS have defined protein partners of some mt proteins, but are too technically challenging and laborious to be practical for analyzing large numbers of samples at the proteome level, e.g., for the study of neuronal or brain-specific mt assemblies, as well as altered mtPPIs on a proteome-wide scale for a disease of interest in brain regions, disease tissues or neurons derived from patients. To address this challenge, we adapted a co-fractionation-MS platform to survey native mt assemblies in adult mouse brain and in human NTERA-2 embryonal carcinoma stem cells or differentiated neuronal-like cells. The workflow consists of orthogonal separations of mt extracts isolated from chemically cross-linked samples to stabilize PPIs, data-dependent acquisition MS to identify co-eluted mt protein profiles from collected fractions and a computational scoring pipeline to predict mtPPIs, followed by network partitioning to define complexes linked to mt functions as well as those essential for neuronal and brain physiological homeostasis. We developed an R/CRAN software package, Macromolecular Assemblies from Co-elution Profiles for automated scoring of co-fractionation-MS data to define complexes from mtPPI networks. Presently, the co-fractionation-MS procedure takes 1.5-3.5 d of proteomic sample preparation, 31 d of MS data acquisition and 8.5 d of data analyses to produce meaningful biological insights.

A dual function of the CRISPR-Cas system in bacterial antivirus immunity and DNA repair.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) and the associated proteins (Cas) comprise a system of adaptive immunity against viruses and plasmids in prokaryotes. Cas1 is a CRISPR-associated protein that is common to all CRISPR-containing prokaryotes but its function remains obscure. Here we show that the purified Cas1 protein of Escherichia coli (YgbT) exhibits nuclease activity against single-stranded and branched DNAs including Holliday junctions, replication forks and 5'-flaps. The crystal structure of YgbT and site-directed mutagenesis have revealed the potential active site. Genome-wide screens show that YgbT physically and genetically interacts with key components of DNA repair systems, including recB, recC and ruvB. Consistent with these findings, the ygbT deletion strain showed increased sensitivity to DNA damage and impaired chromosomal segregation. Similar phenotypes were observed in strains with deletion of CRISPR clusters, suggesting that the function of YgbT in repair involves interaction with the CRISPRs. These results show that YgbT belongs to a novel, structurally distinct family of nucleases acting on branched DNAs and suggest that, in addition to antiviral immunity, at least some components of the CRISPR-Cas system have a function in DNA repair.

Structural biology of the invasion arsenal of Gram-negative bacterial pathogens.

In the last several years, there has been a tremendous progress in the understanding of host-pathogen interactions and the mechanisms by which bacterial pathogens modulate behavior of the host cell. Pathogens use secretion systems to inject a set of proteins, called effectors, into the cytosol of the host cell. These effectors are secreted in a highly regulated, temporal manner and interact with host proteins to modify a multitude of cellular processes. The number of effectors varies between pathogens from ~ 30 to as many as ~ 350. The functional redundancy of effectors encoded by each pathogen makes it difficult to determine the cellular effects or function of individual effectors, since their individual knockouts frequently produce no easily detectable phenotypes. Structural biology of effector proteins and their interactions with host proteins, in conjunction with cell biology approaches, has provided invaluable information about the cellular function of effectors and underlying molecular mechanisms of their modes of action. Many bacterial effectors are functionally equivalent to host proteins while being structurally divergent from them. Other effector proteins display new, previously unobserved functionalities. Here, we summarize the contribution of the structural characterization of effectors and effector-host protein complexes to our understanding of host subversion mechanisms used by the most commonly investigated Gram-negative bacterial pathogens. We describe in some detail the enzymatic activities discovered among effector proteins and how they affect various cellular processes.

HPiP: an R/Bioconductor package for predicting host-pathogen protein-protein interactions from protein sequences using ensemble machine learning approach.

Motivation: Despite arduous and time-consuming experimental efforts, protein-protein interactions (PPIs) for many pathogenic microbes with their human host are still unknown, limiting our understanding of the intricate interactions during infection and the identification of therapeutic targets. Since computational tools offer a promising alternative, we developed an R/Bioconductor package, HPiP (Host-Pathogen Interaction Prediction) software with a series of amino acid sequence property descriptors and an ensemble machine learning classifiers to predict the yet unmapped interactions between pathogen and host proteins. Results: Using severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1) or the novel SARS-CoV-2 coronavirus-human PPI training sets as a case study, we show that HPiP achieves a good performance with PPI predictions between SARS-CoV-2 and human proteins, which we confirmed experimentally in human monocyte THP-1 cells, and with several quality control metrics. HPiP also exhibited strong performance in accurately predicting the previously reported PPIs when tested against the sequences of pathogenic bacteria, Mycobacterium tuberculosis and human proteins. Collectively, our fully documented HPiP software will hasten the exploration of PPIs for a systems-level understanding of many understudied pathogens and uncover molecular targets for repurposing existing drugs. Availability and implementation: HPiP is released as an open-source code under the MIT license that is freely available on GitHub (https://github.com/BabuLab-UofR/HPiP) as well as on Bioconductor (http://bioconductor.org/packages/devel/bioc/html/HPiP.html). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Auxotrophic and prototrophic conditional genetic networks reveal the rewiring of transcription factors in Escherichia coli.

Bacterial transcription factors (TFs) are widely studied in Escherichia coli. Yet it remains unclear how individual genes in the underlying pathways of TF machinery operate together during environmental challenge. Here, we address this by applying an unbiased, quantitative synthetic genetic interaction (GI) approach to measure pairwise GIs among all TF genes in E. coli under auxotrophic (rich medium) and prototrophic (minimal medium) static growth conditions. The resulting static and differential GI networks reveal condition-dependent GIs, widespread changes among TF genes in metabolism, and new roles for uncharacterized TFs (yjdC, yneJ, ydiP) as regulators of cell division, putrescine utilization pathway, and cold shock adaptation. Pan-bacterial conservation suggests TF genes with GIs are co-conserved in evolution. Together, our results illuminate the global organization of E. coli TFs, and remodeling of genetic backup systems for TFs under environmental change, which is essential for controlling the bacterial transcriptional regulatory circuits.

Bioinformatic Analysis of Temporal and Spatial Proteome Alternations During Infections.

Microbial pathogens have evolved numerous mechanisms to hijack host's systems, thus causing disease. This is mediated by alterations in the combined host-pathogen proteome in time and space. Mass spectrometry-based proteomics approaches have been developed and tailored to map disease progression. The result is complex multidimensional data that pose numerous analytic challenges for downstream interpretation. However, a systematic review of approaches for the downstream analysis of such data has been lacking in the field. In this review, we detail the steps of a typical temporal and spatial analysis, including data pre-processing steps (i.e., quality control, data normalization, the imputation of missing values, and dimensionality reduction), different statistical and machine learning approaches, validation, interpretation, and the extraction of biological information from mass spectrometry data. We also discuss current best practices for these steps based on a collection of independent studies to guide users in selecting the most suitable strategies for their dataset and analysis objectives. Moreover, we also compiled the list of commonly used R software packages for each step of the analysis. These could be easily integrated into one's analysis pipeline. Furthermore, we guide readers through various analysis steps by applying these workflows to mock and host-pathogen interaction data from public datasets. The workflows presented in this review will serve as an introduction for data analysis novices, while also helping established users update their data analysis pipelines. We conclude the review by discussing future directions and developments in temporal and spatial proteomics and data analysis approaches. Data analysis codes, prepared for this review are available from https://github.com/BabuLab-UofR/TempSpac, where guidelines and sample datasets are also offered for testing purposes.

Quantitative Genetic Screens for Mapping Bacterial Pathways and Functional Networks.

Escherichia coli synthetic genetic array (eSGA) screening procedure enables high-throughput systematic mapping of pairwise genetic interactions in E. coli. The eSGA method exploits E. coli's rapid growth, its ease of genetic manipulation, and efficient genetic exchange via conjugation. Replica pinning is used to grow and mate arrayed sets of single gene mutant strains as well as to select double mutants en masse. Strain fitness, which is the eSGA readout, is determined by the digital imaging of the plates and subsequent colony size measurements. Comparing single and double mutant colony sizes then allows for identifying interacting genes. Using eSGA on a global or a smaller process-centric scale can help reveal gene functions and reconstruct genetic interaction networks with known and novel connections between genes and pathways.

Structural and functional study of Legionella pneumophila effector RavA.

Legionella pneumophila is an intracellular pathogen that causes Legionnaire's disease in humans. This bacterium can be found in freshwater environments as a free-living organism, but it is also an intracellular parasite of protozoa. Human infection occurs when inhaled aerosolized pathogen comes into contact with the alveolar mucosa and replicates in alveolar macrophages. Legionella enters the host cell by phagocytosis and redirects the Legionella-containing phagosomes from the phagocytic maturation pathway. These nascent phagosomes fuse with ER-derived secretory vesicles and membranes forming the Legionella-containing vacuole. Legionella subverts many host cellular processes by secreting over 300 effector proteins into the host cell via the Dot/Icm type IV secretion system. The cellular function for many Dot/Icm effectors is still unknown. Here, we present a structural and functional study of L. pneumophila effector RavA (Lpg0008). Structural analysis revealed that the RavA consists of four ~85 residue long α-helical domains with similar folds, which show only a low level of structural similarity to other protein domains. The ~90 residues long C-terminal segment is predicted to be natively unfolded. We show that during L. pneumophila infection of human cells, RavA localizes to the Golgi apparatus and to the plasma membrane. The same localization is observed when RavA is expressed in human cells. The localization signal resides within the C-terminal sequence C409 WTSFCGLF417 . Yeast-two-hybrid screen using RavA as bait identified RAB11A as a potential binding partner. RavA is present in L. pneumophila strains but only distant homologs are found in other Legionella species, where the number of repeats varies.

Chemical-genetic profiling reveals limited cross-resistance between antimicrobial peptides with different modes of action.

Antimicrobial peptides (AMPs) are key effectors of the innate immune system and promising therapeutic agents. Yet, knowledge on how to design AMPs with minimal cross-resistance to human host-defense peptides remains limited. Here, we systematically assess the resistance determinants of Escherichia coli against 15 different AMPs using chemical-genetics and compare to the cross-resistance spectra of laboratory-evolved AMP-resistant strains. Although generalizations about AMP resistance are common in the literature, we find that AMPs with different physicochemical properties and cellular targets vary considerably in their resistance determinants. As a consequence, cross-resistance is prevalent only between AMPs with similar modes of action. Finally, our screen reveals several genes that shape susceptibility to membrane- and intracellular-targeting AMPs in an antagonistic manner. We anticipate that chemical-genetic approaches could inform future efforts to minimize cross-resistance between therapeutic and human host AMPs.

Rewiring of the Human Mitochondrial Interactome during Neuronal Reprogramming Reveals Regulators of the Respirasome and Neurogenesis.

Mitochondrial protein (MP) assemblies undergo alterations during neurogenesis, a complex process vital in brain homeostasis and disease. Yet which MP assemblies remodel during differentiation remains unclear. Here, using mass spectrometry-based co-fractionation profiles and phosphoproteomics, we generated mitochondrial interaction maps of human pluripotent embryonal carcinoma stem cells and differentiated neuronal-like cells, which presented as two discrete cell populations by single-cell RNA sequencing. The resulting networks, encompassing 6,442 high-quality associations among 600 MPs, revealed widespread changes in mitochondrial interactions and site-specific phosphorylation during neuronal differentiation. By leveraging the networks, we show the orphan C20orf24 as a respirasome assembly factor whose disruption markedly reduces respiratory chain activity in patients deficient in complex IV. We also find that a heme-containing neurotrophic factor, neuron-derived neurotrophic factor [NENF], couples with Parkinson disease-related proteins to promote neurotrophic activity. Our results provide insights into the dynamic reorganization of mitochondrial networks during neuronal differentiation and highlights mechanisms for MPs in respirasome, neuronal function, and mitochondrial diseases.

Identification and molecular characterization of two naturally occurring Soybean mosaic virus isolates that are closely related but differ in their ability to overcome Rsv4 resistance.

A naturally occurring Rsv4 resistance-breaking isolate (L-RB) and a closely related non-resistance-breaking isolate (L) of Soybean mosaic virus (SMV) were identified in soybean fields in London, Ontario, Canada. The viral genomes of L and L-RB were completely sequenced. Each isolate has a 9585-nucleotide genome with a single open reading frame encoding a polyprotein of approximately 350 kDa. L-RB and L have a very high sequence similarity (99.6%) at both the nucleotide and amino acid levels. Phylogenetic analysis showed that the two isolates belong to the G2 pathotype. Pathogenicity predictions of all virus/soybean combinations, based on the phylogenetic profile, were confirmed by pathogenicity tests using L and L-RB isolates and soybeans carrying different resistance genes, with an exception that L-RB infected a soybean cultivar carrying Rsv4 resistance. The temporal and spatial proximity of L and L-RB and their high sequence similarity suggest L-RB was likely derived from the SMV-L quasispecies. Recombination analysis did not reveal the evidence of genetic recombination for the emergence of L-RB. Mutations introduced by virus-encoded RNA-dependent RNA polymerase during viral genome replication and selection pressure probably contributed to the occurrence of L-RB.