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.

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.

The bacterium E. coli has around 300 transcriptional factors, but the functions of many of them, and the interactions between their respective regulatory networks, are unclear. Here, the authors study genetic interactions among all transcription factor genes in E. coli, revealing condition-dependent interactions and roles for uncharacterized transcription factors.

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.

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.

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.

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 C⁴⁰⁹ WTSFCGLF⁴¹⁷ . 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.

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.

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.

•

Rewiring of mitochondrial (mt) protein interaction network in distinct cell states

•

Dramatic changes in site-specific phosphorylation during neuronal differentiation

•

C20orf24 is a respirasome assembly factor depleted in patients deficient in CIV

•

NENF binding with DJ-1/PINK1 promotes neurotrophic activity and neuronal 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.

Systematic Genetic Screens Reveal the Dynamic Global Functional Organization of the Bacterial Translation Machinery

Bacterial protein synthesis is an essential, conserved, and environmentally responsive process. Yet, many of its components and dependencies remain unidentified. To address this gap, we used quantitative synthetic genetic arrays to map functional relationships among >48,000 gene pairs in Escherichia coli under four culture conditions differing in temperature and nutrient availability. The resulting data provide global functional insights into the roles and associations of genes, pathways, and processes important for efficient translation, growth, and environmental adaptation. We predict and independently verify the requirement of unannotated genes for normal translation, including a previously unappreciated role of YhbY in 30S biogenesis. Dynamic changes in the patterns of genetic dependencies across the four growth conditions and data projections onto other species reveal overarching functional and evolutionary pressures impacting the translation system and bacterial fitness, underscoring the utility of systematic screens for investigating protein synthesis, adaptation, and evolution.

Conditional Epistatic Interaction Maps Reveal Global Functional Rewiring of Genome Integrity Pathways in Escherichia coli

As antibiotic resistance is increasingly becoming a public health concern, an improved understanding of the bacterial DNA damage response (DDR), which is commonly targeted by antibiotics, could be of tremendous therapeutic value. Although the genetic components of the bacterial DDR have been studied extensively in isolation, how the underlying biological pathways interact functionally remains unclear. Here, we address this by performing systematic, unbiased, quantitative synthetic genetic interaction (GI) screens and uncover widespread changes in the GI network of the entire genomic integrity apparatus of Escherichia coli under standard and DNA-damaging growth conditions. The GI patterns of untreated cultures implicated two previously uncharacterized proteins (YhbQ and YqgF) as nucleases, whereas reorganization of the GI network after DNA damage revealed DDR roles for both annotated and uncharacterized genes. Analyses of pan-bacterial conservation patterns suggest that DDR mechanisms and functional relationships are near universal, highlighting a modular and highly adaptive genomic stress response.

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.

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.

Genome-wide genetic interaction (GI) screens have been performed in yeast, but no analogous large-scale studies have yet been reported for bacteria. Here, we have used E. coli synthetic genetic array (eSGA) technology developed by our group to quantitatively map GIs to reveal epistatic dependencies and functional cross-talk among ∼600,000 digenic mutant combinations. By combining this epistasis information with functional modules derived by our group's earlier efforts from proteomic and genomic context (GC)-based methods, we identify several unexpected pathway-level dependencies, functional links between protein complexes, and biological roles of uncharacterized bacterial gene products. As part of the study, two of our pathway predictions from GI screens were validated experimentally, where we confirmed the role of these new components in iron-sulphur biogenesis and ribosome integrity. We also extrapolated the epistatic connectivity diagram of E. coli to 233 distantly related γ-proteobacterial species lacking GI information, and identified co-conserved genes and functional modules important for bacterial pathogenesis. Overall, this study describes the first genome-scale map of GIs in gram-negative bacterium, and through integrative analysis with previously derived protein-protein and GC-based interaction networks presents a number of novel insights into the architecture of bacterial pathways that could not have been discerned through either network alone.

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.

Mapping bacterial functional networks and pathways in Escherichia Coli using synthetic genetic arrays

Phenotypes are determined by a complex series of physical (e.g. protein-protein) and functional (e.g. gene-gene or genetic) interactions (GI)(1). While physical interactions can indicate which bacterial proteins are associated as complexes, they do not necessarily reveal pathway-level functional relationships1. GI screens, in which the growth of double mutants bearing two deleted or inactivated genes is measured and compared to the corresponding single mutants, can illuminate epistatic dependencies between loci and hence provide a means to query and discover novel functional relationships(2). Large-scale GI maps have been reported for eukaryotic organisms like yeast(3-7), but GI information remains sparse for prokaryotes(8), which hinders the functional annotation of bacterial genomes. To this end, we and others have developed high-throughput quantitative bacterial GI screening methods(9, 10). Here, we present the key steps required to perform quantitative E. coli Synthetic Genetic Array (eSGA) screening procedure on a genome-scale(9), using natural bacterial conjugation and homologous recombination to systemically generate and measure the fitness of large numbers of double mutants in a colony array format. Briefly, a robot is used to transfer, through conjugation, chloramphenicol (Cm) - marked mutant alleles from engineered Hfr (High frequency of recombination) 'donor strains' into an ordered array of kanamycin (Kan) - marked F- recipient strains. Typically, we use loss-of-function single mutants bearing non-essential gene deletions (e.g. the 'Keio' collection(11)) and essential gene hypomorphic mutations (i.e. alleles conferring reduced protein expression, stability, or activity(9, 12, 13)) to query the functional associations of non-essential and essential genes, respectively. After conjugation and ensuing genetic exchange mediated by homologous recombination, the resulting double mutants are selected on solid medium containing both antibiotics. After outgrowth, the plates are digitally imaged and colony sizes are quantitatively scored using an in-house automated image processing system(14). GIs are revealed when the growth rate of a double mutant is either significantly better or worse than expected(9). Aggravating (or negative) GIs often result between loss-of-function mutations in pairs of genes from compensatory pathways that impinge on the same essential process(2). Here, the loss of a single gene is buffered, such that either single mutant is viable. However, the loss of both pathways is deleterious and results in synthetic lethality or sickness (i.e. slow growth). Conversely, alleviating (or positive) interactions can occur between genes in the same pathway or protein complex(2) as the deletion of either gene alone is often sufficient to perturb the normal function of the pathway or complex such that additional perturbations do not reduce activity, and hence growth, further. Overall, systematically identifying and analyzing GI networks can provide unbiased, global maps of the functional relationships between large numbers of genes, from which pathway-level information missed by other approaches can be inferred(9).

Genome-scale genetic manipulation methods for exploring bacterial molecular biology

Bacteria are diverse and abundant, playing key roles in human health and disease, the environment, and biotechnology. Despite progress in genome sequencing and bioengineering, much remains unknown about the functional organization of prokaryotes. For instance, roughly a third of the protein-coding genes of the best-studied model bacterium, Escherichia coli, currently lack experimental annotations. Systems-level experimental approaches for investigating the functional associations of bacterial genes and genetic structures are essential for defining the fundamental molecular biology of microbes, preventing the spread of antibacterial resistance in the clinic, and driving the development of future biotechnological applications. This review highlights recently introduced large-scale genetic manipulation and screening procedures for the systematic exploration of bacterial gene functions, molecular relationships, and the global organization of bacteria at the gene, pathway, and genome levels.

Ribosome-dependent ATPase interacts with conserved membrane protein in Escherichia coli to modulate protein synthesis and oxidative phosphorylation

Elongation factor RbbA is required for ATP-dependent deacyl-tRNA release presumably after each peptide bond formation; however, there is no information about the cellular role. Proteomic analysis in Escherichia coli revealed that RbbA reciprocally co-purified with a conserved inner membrane protein of unknown function, YhjD. Both proteins are also physically associated with the 30S ribosome and with members of the lipopolysaccharide transport machinery. Genome-wide genetic screens of rbbA and yhjD deletion mutants revealed aggravating genetic interactions with mutants deficient in the electron transport chain. Cells lacking both rbbA and yhjD exhibited reduced cell division, respiration and global protein synthesis as well as increased sensitivity to antibiotics targeting the ETC and the accuracy of protein synthesis. Our results suggest that RbbA appears to function together with YhjD as part of a regulatory network that impacts bacterial oxidative phosphorylation and translation efficiency.

Array-based synthetic genetic screens to map bacterial pathways and functional networks in Escherichia coli

Cellular processes are carried out through a series of molecular interactions. Various experimental approaches can be used to investigate these functional relationships on a large-scale. Recently, the power of investigating biological systems from the perspective of genetic (gene-gene or epistatic) interactions has been evidenced by the ability to elucidate novel functional relationships. Examples of functionally related genes include genes that buffer each other's function or impinge on the same biological process. Genetic interactions have traditionally been investigated in bacteria by combining pairs of mutations (e.g., gene deletions) and assessing deviation of the phenotype of each double mutant from an expected neutral (or no interaction) phenotype. Fitness is a particularly convenient phenotype to measure: when the double mutant grows faster or slower than expected, the two mutated genes are said to show alleviating or aggravating interactions, respectively. The most commonly used neutral model assumes that the fitness of the double mutant is equal to the product of individual single mutant fitness. A striking genetic interaction is exemplified by the loss of two nonessential genes that buffer each other in performing an essential biological function: deleting only one of these genes produces no detectable fitness defect; however, loss of both genes simultaneously results in systems failure, leading to synthetic sickness or lethality. Systematic large-scale genetic interaction screens have been used to generate functional maps for model eukaryotic organisms, such as yeast, to describe the functional organization of gene products into pathways and protein complexes within a cell. They also reveal the modular arrangement and cross talk of pathways and complexes within broader functional neighborhoods (Dixon et al., Annu Rev Genet 43:601-625, 2009). Here, we present a high-throughput quantitative Escherichia coli Synthetic Genetic Array (eSGA) screening procedure, which we developed to systematically infer genetic interactions by scoring growth defects among large numbers of double mutants in a classic Gram-negative bacterium. The eSGA method exploits the rapid colony growth, ease of genetic manipulation, and natural efficient genetic exchange via conjugation of laboratory E. coli strains. Replica pinning is used to grow and mate arrayed sets of single gene mutant strains and to select double mutants en masse. Strain fitness, which is used as the eSGA readout, is quantified by the digital imaging of the plates and subsequent measuring and comparing single and double mutant colony sizes. While eSGA can be used to screen select mutants to probe the functions of individual genes, using eSGA more broadly to collect genetic interaction data for many combinations of genes can help reconstruct a functional interaction network to reveal novel links and components of biological pathways as well as unexpected connections between pathways. A variety of bacterial systems can be investigated, wherein the genes impinge on a essential biological process (e.g., cell wall assembly, ribosome biogenesis, chromosome replication) that are of interest from the perspective of drug development (Babu et al., Mol Biosyst 12:1439-1455, 2009). We also show how genetic interactions generated by high-throughput eSGA screens can be validated by manual small-scale genetic crosses and by genetic complementation and gene rescue experiments.

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.