RNA interactome capture in Escherichia coli globally identifies RNA-binding proteins.

RNA-binding proteins (RPBs) are deeply involved in fundamental cellular processes in bacteria and are vital for their survival. Despite this, few studies have so far been dedicated to direct and global identification of bacterial RBPs. We have adapted the RNA interactome capture (RIC) technique, originally developed for eukaryotic systems, to globally identify RBPs in bacteria. RIC takes advantage of the base pairing potential of poly(A) tails to pull-down RNA-protein complexes. Overexpressing poly(A) polymerase I in Escherichia coli drastically increased transcriptome-wide RNA polyadenylation, enabling pull-down of crosslinked RNA-protein complexes using immobilized oligo(dT) as bait. With this approach, we identified 169 putative RBPs, roughly half of which are already annotated as RNA-binding. We experimentally verified the RNA-binding ability of a number of uncharacterized RBPs, including YhgF, which is exceptionally well conserved not only in bacteria, but also in archaea and eukaryotes. We identified YhgF RNA targets in vivo using CLIP-seq, verified specific binding in vitro, and reveal a putative role for YhgF in regulation of gene expression. Our findings present a simple and robust strategy for RBP identification in bacteria, provide a resource of new bacterial RBPs, and lay the foundation for further studies of the highly conserved RBP YhgF.

RHS-elements function as type II toxin-antitoxin modules that regulate intra-macrophage replication of Salmonella Typhimurium.

RHS elements are components of conserved toxin-delivery systems, wide-spread within the bacterial kingdom and some of the most positively selected genes known. However, very little is known about how Rhs toxins affect bacterial biology. Salmonella Typhimurium contains a full-length rhs gene and an adjacent orphan rhs gene, which lacks the conserved delivery part of the Rhs protein. Here we show that, in addition to the conventional delivery, Rhs toxin-antitoxin pairs encode for functional type-II toxin-antitoxin (TA) loci that regulate S. Typhimurium proliferation within macrophages. Mutant S. Typhimurium cells lacking both Rhs toxins proliferate 2-times better within macrophages, mainly because of an increased growth rate. Thus, in addition to providing strong positive selection for the rhs loci under conditions when there is little or no toxin delivery, internal expression of the toxin-antitoxin system regulates growth in the stressful environment found inside macrophages.

Contact-dependent growth inhibition induces high levels of antibiotic-tolerant persister cells in clonal bacterial populations.

Bacterial populations can use bet-hedging strategies to cope with rapidly changing environments. One example is non-growing cells in clonal bacterial populations that are able to persist antibiotic treatment. Previous studies suggest that persisters arise in bacterial populations either stochastically through variation in levels of global signalling molecules between individual cells, or in response to various stresses. Here, we show that toxins used in contact-dependent growth inhibition (CDI) create persisters upon direct contact with cells lacking sufficient levels of CdiI immunity protein, which would otherwise bind to and neutralize toxin activity. CDI-mediated persisters form through a feedforward cycle where the toxic activity of the CdiA toxin increases cellular (p)ppGpp levels, which results in Lon-mediated degradation of the immunity protein and more free toxin. Thus, CDI systems mediate a population density-dependent bet-hedging strategy, where the fraction of non-growing cells is increased only when there are many cells of the same genotype. This may be one of the mechanisms of how CDI systems increase the fitness of their hosts.

Class II contact-dependent growth inhibition (CDI) systems allow for broad-range cross-species toxin delivery within the Enterobacteriaceae family.

Contact-dependent growth inhibition (CDI) allows bacteria to recognize kin cells in mixed bacterial populations. In Escherichia coli, CDI mediated effector delivery has been shown to be species-specific, with a preference for the own strain over others. This specificity is achieved through an interaction between a receptor-binding domain in the CdiA protein and its cognate receptor protein on the target cell. But how conserved this specificity is has not previously been investigated in detail. Here, we show that class II CdiA receptor-binding domains and their Enterobacter cloacae analog are highly promiscuous, and can allow for efficient effector delivery into several different Enterobacteriaceae species, including Escherichia, Enterobacter, Klebsiella and Salmonella spp. In addition, although we observe a preference for the own receptors over others for two of the receptor-binding domains, this did not limit cross-species effector delivery in all experimental conditions. These results suggest that class II CdiA proteins could allow for broad-range and cross-species growth inhibition in mixed bacterial populations.

Functional plasticity of antibacterial EndoU toxins.

Bacteria use several different secretion systems to deliver toxic EndoU ribonucleases into neighboring cells. Here, we present the first structure of a prokaryotic EndoU toxin in complex with its cognate immunity protein. The contact-dependent growth inhibition toxin CdiA-CTSTECO31 from Escherichia coli STEC_O31 adopts the eukaryotic EndoU fold and shares greatest structural homology with the nuclease domain of coronavirus Nsp15. The toxin contains a canonical His-His-Lys catalytic triad in the same arrangement as eukaryotic EndoU domains, but lacks the uridylate-specific ribonuclease activity that characterizes the superfamily. Comparative sequence analysis indicates that bacterial EndoU domains segregate into at least three major clades based on structural variations in the N-terminal subdomain. Representative EndoU nucleases from clades I and II degrade tRNA molecules with little specificity. In contrast, CdiA-CTSTECO31 and other clade III toxins are specific anticodon nucleases that cleave tRNAGlu between nucleotides C37 and m2 A38. These findings suggest that the EndoU fold is a versatile scaffold for the evolution of novel substrate specificities. Such functional plasticity may account for the widespread use of EndoU effectors by diverse inter-bacterial toxin delivery systems.

Selection of orphan Rhs toxin expression in evolved Salmonella enterica serovar Typhimurium.

Clonally derived bacterial populations exhibit significant genotypic and phenotypic diversity that contribute to fitness in rapidly changing environments. Here, we show that serial passage of Salmonella enterica serovar Typhimurium LT2 (StLT2) in broth, or within a mouse host, results in selection of an evolved population that inhibits the growth of ancestral cells by direct contact. Cells within each evolved population gain the ability to express and deploy a cryptic "orphan" toxin encoded within the rearrangement hotspot (rhs) locus. The Rhs orphan toxin is encoded by a gene fragment located downstream of the "main" rhs gene in the ancestral strain StLT2. The Rhs orphan coding sequence is linked to an immunity gene, which encodes an immunity protein that specifically blocks Rhs orphan toxin activity. Expression of the Rhs orphan immunity protein protects ancestral cells from the evolved lineages, indicating that orphan toxin activity is responsible for the observed growth inhibition. Because the Rhs orphan toxin is encoded by a fragmented reading frame, it lacks translation initiation and protein export signals. We provide evidence that evolved cells undergo recombination between the main rhs gene and the rhs orphan toxin gene fragment, yielding a fusion that enables expression and delivery of the orphan toxin. In this manner, rhs locus rearrangement provides a selective advantage to a subpopulation of cells. These observations suggest that rhs genes play important roles in intra-species competition and bacterial evolution.

Rhs proteins from diverse bacteria mediate intercellular competition.

Rearrangement hotspot (Rhs) and related YD-peptide repeat proteins are widely distributed in bacteria and eukaryotes, but their functions are poorly understood. Here, we show that Gram-negative Rhs proteins and the distantly related wall-associated protein A (WapA) from Gram-positive bacteria mediate intercellular competition. Rhs and WapA carry polymorphic C-terminal toxin domains (Rhs-CT/WapA-CT), which are deployed to inhibit the growth of neighboring cells. These systems also encode sequence-diverse immunity proteins (RhsI/WapI) that specifically neutralize cognate toxins to protect rhs(+)/wapA(+) cells from autoinhibition. RhsA and RhsB from Dickeya dadantii 3937 carry nuclease domains that degrade target cell DNA. D. dadantii 3937 rhs genes do not encode secretion signal sequences but are linked to hemolysin-coregulated protein and valine-glycine repeat protein G genes from type VI secretion systems. Valine-glycine repeat protein G is required for inhibitor cell function, suggesting that Rhs may be exported from D. dadantii 3937 through a type VI secretion mechanism. In contrast, WapA proteins from Bacillus subtilis strains appear to be exported through the general secretory pathway and deliver a variety of tRNase toxins into neighboring target cells. These findings demonstrate that YD-repeat proteins from phylogenetically diverse bacteria share a common function in contact-dependent growth inhibition.

Selection-driven gene loss in bacteria.

Gene loss by deletion is a common evolutionary process in bacteria, as exemplified by bacteria with small genomes that have evolved from bacteria with larger genomes by reductive processes. The driving force(s) for genome reduction remains unclear, and here we examined the hypothesis that gene loss is selected because carriage of superfluous genes confers a fitness cost to the bacterium. In the bacterium Salmonella enterica, we measured deletion rates at 11 chromosomal positions and the fitness effects of several spontaneous deletions. Deletion rates varied over 200-fold between different regions with the replication terminus region showing the highest rates. Approximately 25% of the examined deletions caused an increase in fitness under one or several growth conditions, and after serial passage of wild-type bacteria in rich medium for 1,000 generations we observed fixation of deletions that substantially increased bacterial fitness when reconstructed in a non-evolved bacterium. These results suggest that selection could be a significant driver of gene loss and reductive genome evolution.

Chromosome-level genome assembly and annotation of the social amoeba Dictyostelium firmibasis.

Dicytostelium firmibasis is a member of Dictyostelia, a group of social amoebae that upon starvation display aggregative multicellularity where the amoebae transition from uni- to multicellular life. The D. firmibasis genome assembly that is currently available is of limited use due to its low contiguity, large number of undetermined bases, and lack of annotations. Here we used Nanopore long read sequencing, complemented with Illumina sequencing, and developmental transcriptomics as well as small RNA-sequencing, to present a new, fully annotated, chromosome-level D. firmibasis genome assembly. The new assembly contains no undetermined bases, and consists mainly of six large contigs representing the chromosomes, as well as a complete mitochondrial genome. This new genome assembly will be a valuable tool, allowing comprehensive comparison to Dictyostelium discoideum, the dictyostelid genetically tractable model. Further, the new genome will be important for studies of evolutionary processes governing the transition from unicellular to multicellular organisms and will aid in the sequencing and annotation of other dictyostelids genomes, many of which are currently of poor quality.

Colicins and T6SS-based competition systems enhance enterotoxigenic E. coli (ETEC) competitiveness.

Diarrheal diseases are still a significant problem for humankind, causing approximately half a million deaths annually. To cause diarrhea, enteric bacterial pathogens must first colonize the gut, which is a niche occupied by the normal bacterial microbiota. Therefore, the ability of pathogenic bacteria to inhibit the growth of other bacteria can facilitate the colonization process. Although enterotoxigenic Escherichia coli (ETEC) is one of the major causative agents of diarrheal diseases, little is known about the competition systems found in and used by ETEC and how they contribute to the ability of ETEC to colonize a host. Here, we collected a set of 94 fully assembled ETEC genomes by performing whole-genome sequencing and mining the NCBI RefSeq database. Using this set, we performed a comprehensive search for delivered bacterial toxins and investigated how these toxins contribute to ETEC competitiveness in vitro. We found that type VI secretion systems (T6SS) were widespread among ETEC (n = 47). In addition, several closely related ETEC strains were found to encode Colicin Ia and T6SS (n = 8). These toxins provide ETEC competitive advantages during in vitro competition against other E. coli, suggesting that the role of T6SS as well as colicins in ETEC biology has until now been underappreciated.

OSA Is Associated With the Human Gut Microbiota Composition and Functional Potential in the Population-Based Swedish CardioPulmonary bioImage Study.

BACKGROUND: OSA is a common sleep-breathing disorder linked to increased risk of cardiovascular disease. Intermittent upper airway obstruction and hypoxia, hallmarks of OSA, have been shown in animal models to induce substantial changes to the gut microbiota composition, and subsequent transplantation of fecal matter to other animals induced changes in BP and glucose metabolism. RESEARCH QUESTION: Does OSA in adults associate with the composition and functional potential of the human gut microbiota? STUDY DESIGN AND METHODS: We used respiratory polygraphy data from up to 3,570 individuals 50 to 64 years of age from the population-based Swedish Cardiopulmonary bioimage Study combined with deep shotgun metagenomics of fecal samples to identify cross-sectional associations between three OSA parameters covering apneas and hypopneas, cumulative sleep time in hypoxia, and number of oxygen desaturation events with gut microbiota composition. Data collection about potential confounders was based on questionnaires, onsite anthropometric measurements, plasma metabolomics, and linkage with the Swedish Prescribed Drug Register. RESULTS: We found that all three OSA parameters were associated with lower diversity of species in the gut. Furthermore, in multivariable-adjusted analysis, the OSA-related hypoxia parameters were associated with the relative abundance of 128 gut bacterial species, including higher abundance of Blautia obeum and Collinsella aerofaciens. The latter species was also independently associated with increased systolic BP. Furthermore, the cumulative time in hypoxia during sleep was associated with the abundance of genes involved in nine gut microbiota metabolic pathways, including propionate production from lactate. Finally, we observed two heterogeneous sets of plasma metabolites with opposite association with species positively and negatively associated with hypoxia parameters, respectively. INTERPRETATION: OSA-related hypoxia, but not the number of apneas/hypopneas, is associated with specific gut microbiota species and functions. Our findings lay the foundation for future research on the gut microbiota-mediated health effects of OSA.

Activation of cryptic aminoglycoside resistance in Salmonella enterica.

Aminoglycoside resistance in bacteria can be acquired by several mechanisms, including drug modification, target alteration, reduced uptake and increased efflux. Here we demonstrate that increased resistance to the aminoglycosides streptomycin and spectinomycin in Salmonella enterica can be conferred by increased expression of an aminoglycoside adenyl transferase encoded by the cryptic, chromosomally located aadA gene. During growth in rich medium the wild-type strain was susceptible but mutations that impaired electron transport and conferred a small colony variant (SCV) phenotype or growth in glucose/glycerol minimal media resulted in activation of the aadA gene and aminoglycoside resistance. Expression of the aadA gene was positively regulated by the stringent response regulator guanosine penta/tetraphosphate ((p)ppGpp). SCV mutants carrying stop codon mutations in the hemA and ubiA genes showed a streptomycin pseudo-dependent phenotype, where growth was stimulated by streptomycin. Our data suggest that this phenotype is due to streptomycin-induced readthrough of the stop codons, a resulting increase in HemA/UbiA levels and improved electron transport and growth. Our results demonstrate that environmental and mutational activation of a cryptic resistance gene can confer clinically significant resistance and that a streptomycin-pseudo-dependent phenotype can be generated via a novel mechanism that does not involve the classical rpsL mutations.

Translesion DNA polymerases are required for spontaneous deletion formation in Salmonella typhimurium.

How spontaneous deletions form in bacteria is still a partly unresolved problem. Here, we show that deletion formation in Salmonella typhimurium requires the presence of functional translesion polymerases. First, in wild-type bacteria, removal of the known translesion DNA polymerases, PolII (polB), PolIV (dinB), PolV (umuDC), and SamAB (samAB), resulted in a 10-fold decrease in the deletion rate, indicating that 90% of all spontaneous deletions require these polymerases for their formation. Second, overexpression of these polymerases by derepression of the DNA damage-inducible LexA regulon caused a 25-fold increase in deletion rate that depended on the presence of functional translesion polymerases. Third, overexpression of the polymerases PolII and PolIV from a plasmid increased the deletion rate 12- to 30-fold, respectively. Last, in a recBC(-) mutant where dsDNA ends are stabilized due to the lack of the end-processing nuclease RecBC, the deletion rate was increased 20-fold. This increase depended on the translesion polymerases. In lexA(def) mutant cells with constitutive SOS expression, a 10-fold increase in DNA breaks was observed. Inactivation of all 4 translesion polymerases in the lexA(def) mutant reduced the deletion rate 250-fold without any concomitant reduction in the amount of DNA breaks. Mutational inactivation of 3 endonucleases under LexA control reduced the number of DNA breaks to the wild-type level in a lexA(def) mutant with a concomitant 50-fold reduction in deletion rate. These findings suggest that the translesion polymerases are not involved in forming the DNA breaks, but that they require them to stimulate deletion formation.

Biological roles of translesion synthesis DNA polymerases in eubacteria.

Biological systems are strongly selected to maintain the integrity of their genomes by prevention and repair of external and internal DNA damages. However, some types of DNA lesions persist and might block the replication apparatus. The universal existence of specialized translesion synthesis DNA polymerases (TLS polymerases) that can bypass such lesions in DNA implies that replication blockage is a general biological problem. We suggest that the primary function for which translesion synthesis polymerases are selected is to rescue cells from replication arrest at lesions in DNA, a situation that, if not amended, is likely to cause an immediate and severe reduction in cell fitness and survival. We will argue that the mutagenesis observed during translesion synthesis is an unavoidable secondary consequence of this primary function and not, as has been suggested, an evolved mechanism to increase mutation rates in response to various stresses. Finally, we will discuss recent data on additional roles for translesion synthesis polymerases in the formation of spontaneous deletions and in transcription-coupled TLS, where the coupling of transcription to TLS is proposed to allow the rescue of the transcription machinery arrested at DNA lesions.

Genetic Evidence for SecY Translocon-Mediated Import of Two Contact-Dependent Growth Inhibition (CDI) Toxins.

The C-terminal (CT) toxin domains of contact-dependent growth inhibition (CDI) CdiA proteins target Gram-negative bacteria and must breach both the outer and inner membranes of target cells to exert growth inhibitory activity. Here, we examine two CdiA-CT toxins that exploit the bacterial general protein secretion machinery after delivery into the periplasm. A Ser281Phe amino acid substitution in transmembrane segment 7 of SecY, the universally conserved channel-forming subunit of the Sec translocon, decreases the cytotoxicity of the membrane depolarizing orphan10 toxin from enterohemorrhagic Escherichia coli EC869. Target cells expressing secYS281F and lacking either PpiD or YfgM, two SecY auxiliary factors, are fully protected from CDI-mediated inhibition either by CdiA-CTo10EC869 or by CdiA-CTGN05224, the latter being an EndoU RNase CdiA toxin from Klebsiella aerogenes GN05224 that has a related cytoplasm entry domain. RNase activity of CdiA-CTGN05224 was reduced in secYS281F target cells and absent in secYS281F ΔppiD or secYS281F ΔyfgM target cells during competition co-cultures. Importantly, an allele-specific mutation in secY (secYG313W ) renders ΔppiD or ΔyfgM target cells specifically resistant to CdiA-CTGN05224 but not to CdiA-CTo10EC869, further suggesting a direct interaction between SecY and the CDI toxins. Our results provide genetic evidence of a unique confluence between the primary cellular export route for unfolded polypeptides and the import pathways of two CDI toxins.IMPORTANCE Many bacterial species interact via direct cell-to-cell contact using CDI systems, which provide a mechanism to inject toxins that inhibit bacterial growth into one another. Here, we find that two CDI toxins, one that depolarizes membranes and another that degrades RNA, exploit the universally conserved SecY translocon machinery used to export proteins for target cell entry. Mutations in genes coding for members of the Sec translocon render cells resistant to these CDI toxins by blocking their movement into and through target cell membranes. This work lays the foundation for understanding how CDI toxins interact with the protein export machinery and has direct relevance to development of new antibiotics that can penetrate bacterial cell envelopes.

Escherichia coli EC93 deploys two plasmid-encoded class I contact-dependent growth inhibition systems for antagonistic bacterial interactions.

The phenomenon of contact-dependent growth inhibition (CDI) and the genes required for CDI (cdiBAI) were identified and isolated in 2005 from an Escherichia coli isolate (EC93) from rats. Although the cdiBAIEC93 locus has been the focus of extensive research during the past 15 years, little is known about the EC93 isolate from which it originates. Here we sequenced the EC93 genome and find two complete and functional cdiBAI loci (including the previously identified cdi locus), both carried on a large 127 kb plasmid. These cdiBAI systems are differentially expressed in laboratory media, enabling EC93 to outcompete E. coli cells lacking cognate cdiI immunity genes. The two CDI systems deliver distinct effector peptides that each dissipate the membrane potential of target cells, although the two toxins display different toxic potencies. Despite the differential expression and toxic potencies of these CDI systems, both yielded similar competitive advantages against E. coli cells lacking immunity. This can be explained by the fact that the less expressed cdiBAI system (cdiBAIEC93-2) delivers a more potent toxin than the highly expressed cdiBAIEC93-1 system. Moreover, our results indicate that unlike most sequenced CDI+ bacterial isolates, the two cdi loci of E. coli EC93 are located on a plasmid and are expressed in laboratory media.

Evolution of high-level resistance during low-level antibiotic exposure.

It has become increasingly clear that low levels of antibiotics present in many environments can select for resistant bacteria, yet the evolutionary pathways for resistance development during exposure to low amounts of antibiotics remain poorly defined. Here we show that Salmonella enterica exposed to sub-MIC levels of streptomycin evolved high-level resistance via novel mechanisms that are different from those observed during lethal selections. During lethal selection only rpsL mutations are found, whereas at sub-MIC selection resistance is generated by several small-effect resistance mutations that combined confer high-level resistance via three different mechanisms: (i) alteration of the ribosomal RNA target (gidB mutations), (ii) reduction in aminoglycoside uptake (cyoB, nuoG, and trkH mutations), and (iii) induction of the aminoglycoside-modifying enzyme AadA (znuA mutations). These results demonstrate how the strength of the selective pressure influences evolutionary trajectories and that even weak selective pressures can cause evolution of high-level resistance.

Fluorescent CRISPR Adaptation Reporter for rapid quantification of spacer acquisition.

CRISPR-Cas systems are adaptive prokaryotic immune systems protecting against horizontally transferred DNA or RNA such as viruses and other mobile genetic elements. Memory of past invaders is stored as spacers in CRISPR loci in a process called adaptation. Here we developed a novel assay where spacer integration results in fluorescence, enabling detection of memory formation in single cells and quantification of as few as 0.05% cells with expanded CRISPR arrays in a bacterial population. Using this fluorescent CRISPR Adaptation Reporter (f-CAR), we quantified adaptation of the two CRISPR arrays of the type I-E CRISPR-Cas system in Escherichia coli, and confirmed that more integration events are targeted to CRISPR-II than to CRISPR-I. The f-CAR conveniently analyzes and compares many samples, allowing new insights into adaptation. For instance, we show that in an E. coli culture the majority of acquisition events occur in late exponential phase.

CdiA Effectors Use Modular Receptor-Binding Domains To Recognize Target Bacteria.

Contact-dependent growth inhibition (CDI) systems encode CdiA effectors, which bind to specific receptors on neighboring bacteria and deliver C-terminal toxin domains to suppress target cell growth. Two classes of CdiA effectors that bind distinct cell surface receptors have been identified, but the molecular basis of receptor specificity is not understood. Alignment of BamA-specific CdiAEC93 from Escherichia coli EC93 and OmpC-specific CdiAEC536 from E. coli 536 suggests that the receptor-binding domain resides within a central region that varies between the two effectors. In support of this hypothesis, we find that CdiAEC93 fragments containing residues Arg1358 to Phe1646 bind specifically to purified BamA. Moreover, chimeric CdiAEC93 that carries the corresponding sequence from CdiAEC536 is endowed with OmpC-binding activity, demonstrating that this region dictates receptor specificity. A survey of E. coli CdiA proteins reveals two additional effector classes, which presumably recognize distinct receptors. Using a genetic approach, we identify the outer membrane nucleoside transporter Tsx as the receptor for a third class of CdiA effectors. Thus, CDI systems exploit multiple outer membrane proteins to identify and engage target cells. These results underscore the modularity of CdiA proteins and suggest that novel effectors can be constructed through genetic recombination to interchange different receptor-binding domains and toxic payloads.IMPORTANCE CdiB/CdiA two-partner secretion proteins mediate interbacterial competition through the delivery of polymorphic toxin domains. This process, known as contact-dependent growth inhibition (CDI), requires stable interactions between the CdiA effector protein and specific receptors on the surface of target bacteria. Here, we localize the receptor-binding domain to the central region of E. coli CdiA. Receptor-binding domains vary between CdiA proteins, and E. coli strains collectively encode at least four distinct effector classes. Further, we show that receptor specificity can be altered by exchanging receptor-binding regions, demonstrating the modularity of this domain. We propose that novel CdiA effectors are naturally generated through genetic recombination to interchange different receptor-binding domains and toxin payloads.