An antiplasmid system drives antibiotic resistance gene integration in carbapenemase-producing Escherichia coli lineages

Plasmids carrying antibiotic resistance genes (ARG) are the main mechanism of resistance dissemination in Enterobacterales. However, the fitness-resistance trade-off may result in their elimination. Chromosomal integration of ARGs preserves resistance advantage while relieving the selective pressure for keeping costly plasmids. In some bacterial lineages, such as carbapenemase producing sequence type ST38 Escherichia coli, most ARGs are chromosomally integrated. Here we reproduce by experimental evolution the mobilisation of the carbapenemase blaOXA-48 gene from the pOXA-48 plasmid into the chromosome. We demonstrate that this integration depends on a plasmid-induced fitness cost, a mobile genetic structure embedding the ARG and a novel antiplasmid system ApsAB actively involved in pOXA-48 destabilization. We show that ApsAB targets high and low-copy number plasmids. ApsAB combines a nuclease/helicase protein and a novel type of Argonaute-like protein. It belongs to a family of defense systems broadly distributed among bacteria, which might have a strong ecological impact on plasmid diffusion.

Rapid expansion and international spread of M1UK in the post-pandemic UK upsurge of Streptococcus pyogenes

The UK observed a marked increase in scarlet fever and invasive group A streptococcal infection in 2022 with severe outcomes in children and similar trends worldwide. Here we report lineage M1UK to be the dominant source of invasive infections in this upsurge. Compared with ancestral M1global strains, invasive M1UK strains exhibit reduced genomic diversity and fewer mutations in two-component regulator genes covRS. The emergence of M1UK is dated to 2008. Following a bottleneck coinciding with the COVID-19 pandemic, three emergent M1UK clades underwent rapid nationwide expansion, despite lack of detection in previous years. All M1UK isolates thus-far sequenced globally have a phylogenetic origin in the UK, with dispersal of the new clades in Europe. While waning immunity may promote streptococcal epidemics, the genetic features of M1UK point to a fitness advantage in pathogenicity, and a striking ability to persist through population bottlenecks.

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Conjugative type IV secretion systems enable bacterial antagonism that operates independently of plasmid transfer

Bacterial cooperation and antagonism mediated by secretion systems are among the ways in which bacteria interact with one another. Here we report the discovery of an antagonistic property of a type IV secretion system (T4SS) sourced from a conjugative plasmid, RP4, using engineering approaches. We scrutinized the genetic determinants and suggested that this antagonistic activity is independent of molecular cargos, while we also elucidated the resistance genes. We further showed that a range of Gram-negative bacteria and a mixed bacterial population can be eliminated by this T4SS-dependent antagonism. Finally, we showed that such an antagonistic property is not limited to T4SS sourced from RP4, rather it can also be observed in a T4SS originated from another conjugative plasmid, namely R388. Our results are the first demonstration of conjugative T4SS-dependent antagonism between Gram-negative bacteria on the genetic level and provide the foundation for future mechanistic studies.

Genetic barriers more than environmental associations explain Serratia marcescens population structure

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

A role for the stringent response in ciprofloxacin resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a major cause of nosocomial infections and the leading cause of chronic lung infections in cystic fibrosis and chronic obstructive pulmonary disease patients. Antibiotic treatment remains challenging because P. aeruginosa is resistant to high concentrations of antibiotics and has a remarkable ability to acquire mutations conferring resistance to multiple groups of antimicrobial agents. Here we report that when P. aeruginosa is plated on ciprofloxacin (cipro) plates, the majority of cipro-resistant (ciproR) colonies observed at and after 48 h of incubation carry mutations in genes related to the Stringent Response (SR). Mutations in one of the major SR components, spoT, were present in approximately 40% of the ciproR isolates. Compared to the wild-type strain, most of these isolates had decreased growth rate, longer lag phase and altered intracellular ppGpp content. Also, 75% of all sequenced mutations were insertions and deletions, with short deletions being the most frequently occurring mutation type. We present evidence that most of the observed mutations are induced on the selective plates in a subpopulation of cells that are not instantly killed by cipro. Our results suggests that the SR may be an important contributor to antibiotic resistance acquisition in P. aeruginosa.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Chromosomal barcodes for simultaneous tracking of near-isogenic bacterial strains in plant microbiota

DNA-amplicon-based microbiota profiling can estimate species diversity and abundance but cannot resolve genetic differences within individuals of the same species. Here we report the development of modular bacterial tags (MoBacTags) encoding DNA barcodes that enable tracking of near-isogenic bacterial commensals in an array of complex microbiome communities. Chromosomally integrated DNA barcodes are then co-amplified with endogenous marker genes of the community by integrating corresponding primer binding sites into the barcode. We use this approach to assess the contributions of individual bacterial genes to Arabidopsis thaliana root microbiota establishment with synthetic communities that include MoBacTag-labelled strains of Pseudomonas capeferrum. Results show reduced root colonization for certain mutant strains with defects in gluconic-acid-mediated host immunosuppression, which would not be detected with traditional amplicon sequencing. Our work illustrates how MoBacTags can be applied to assess scaling of individual bacterial genetic determinants in the plant microbiota.

Compensatory evolution in NusG improves fitness of drug-resistant M. tuberculosis

Drug-resistant bacteria are emerging as a global threat, despite frequently being less fit than their drug-susceptible ancestors1-8. Here we sought to define the mechanisms that drive or buffer the fitness cost of rifampicin resistance (RifR) in the bacterial pathogen Mycobacterium tuberculosis (Mtb). Rifampicin inhibits RNA polymerase (RNAP) and is a cornerstone of modern short-course tuberculosis therapy9,10. However, RifR Mtb accounts for one-quarter of all deaths due to drug-resistant bacteria11,12. We took a comparative functional genomics approach to define processes that are differentially vulnerable to CRISPR interference (CRISPRi) inhibition in RifR Mtb. Among other hits, we found that the universally conserved transcription factor NusG is crucial for the fitness of RifR Mtb. In contrast to its role in Escherichia coli, Mtb NusG has an essential RNAP pro-pausing function mediated by distinct contacts with RNAP and the DNA13. We find this pro-pausing NusG-RNAP interface to be under positive selection in clinical RifR Mtb isolates. Mutations in the NusG-RNAP interface reduce pro-pausing activity and increase fitness of RifR Mtb. Collectively, these results define excessive RNAP pausing as a molecular mechanism that drives the fitness cost of RifR in Mtb, identify a new mechanism of compensation to overcome this cost, suggest rational approaches to exacerbate the fitness cost, and, more broadly, could inform new therapeutic approaches to develop drug combinations to slow the evolution of RifR in Mtb.

Bacteria can compensate the fitness costs of amplified resistance genes via a bypass mechanism

Antibiotic heteroresistance is a phenotype in which a susceptible bacterial population includes a small subpopulation of cells that are more resistant than the main population. Such resistance can arise by tandem amplification of DNA regions containing resistance genes that in single copy are not sufficient to confer resistance. However, tandem amplifications often carry fitness costs, manifested as reduced growth rates. Here, we investigated if and how these fitness costs can be genetically ameliorated. We evolved four clinical isolates of three bacterial species that show heteroresistance to tobramycin, gentamicin and tetracyclines at increasing antibiotic concentrations above the minimal inhibitory concentration (MIC) of the main susceptible population. This led to a rapid enrichment of resistant cells with up to an 80-fold increase in the resistance gene copy number, an increased MIC, and severely reduced growth rates. When further evolved in the presence of antibiotic, these strains acquired compensatory resistance mutations and showed a reduction in copy number while maintaining high-level resistance. A deterministic model indicated that the loss of amplified units was driven mainly by their fitness costs and that the compensatory mutations did not affect the loss rate of the gene amplifications. Our findings suggest that heteroresistance mediated by copy number changes can facilitate and precede the evolution towards stable resistance.

Inter-species gene flow drives ongoing evolution of Streptococcus pyogenes and Streptococcus dysgalactiae subsp. equisimilis

Streptococcus dysgalactiae subsp. equisimilis (SDSE) is an emerging cause of human infection with invasive disease incidence and clinical manifestations comparable to the closely related species, Streptococcus pyogenes. Through systematic genomic analyses of 501 disseminated SDSE strains, we demonstrate extensive overlap between the genomes of SDSE and S. pyogenes. More than 75% of core genes are shared between the two species with one third demonstrating evidence of cross-species recombination. Twenty-five percent of mobile genetic element (MGE) clusters and 16 of 55 SDSE MGE insertion regions were shared across species. Assessing potential cross-protection from leading S. pyogenes vaccine candidates on SDSE, 12/34 preclinical vaccine antigen genes were shown to be present in >99% of isolates of both species. Relevant to possible vaccine evasion, six vaccine candidate genes demonstrated evidence of inter-species recombination. These findings demonstrate previously unappreciated levels of genomic overlap between these closely related pathogens with implications for streptococcal pathobiology, disease surveillance and prevention.

Dependency on host vitamin B12 has shaped Mycobacterium tuberculosis Complex evolution

Human and animal tuberculosis is caused by the Mycobacterium tuberculosis Complex (MTBC), which has evolved a genomic decay of cobalamin (vitamin B12) biosynthetic genes. Accordingly, and in sharp contrast to environmental, opportunistic and ancestor mycobacteria; we demonstrate that M. tuberculosis (Mtb), M. africanum, and animal-adapted lineages, lack endogenous production of cobalamin, yet they retain the capacity for exogenous uptake. A B12 anemic model in immunocompromised and immunocompetent mice, demonstrates improved survival, and lower bacteria in organs, in B12 anemic animals infected with Mtb relative to non-anemic controls. Conversely, no differences were observed between mice groups infected with M. canettii, an ancestor mycobacterium which retains cobalamin biosynthesis. Interrogation of the B12 transcriptome in three MTBC strains defined L-methionine synthesis by metE and metH genes as a key phenotype. Expression of metE is repressed by a cobalamin riboswitch, while MetH requires the cobalamin cofactor. Thus, deletion of metE predominantly attenuates Mtb in anemic mice; although inactivation of metH exclusively causes attenuation in non-anemic controls. Here, we show how sub-physiological levels of B12 in the host antagonizes Mtb virulence, and describe a yet unknown mechanism of host-pathogen cross-talk with implications for B12 anemic populations.

ProTInSeq: transposon insertion tracking by ultra-deep DNA sequencing to identify translated large and small ORFs

Identifying open reading frames (ORFs) being translated is not a trivial task. ProTInSeq is a technique designed to characterize proteomes by sequencing transposon insertions engineered to express a selection marker when they occur in-frame within a protein-coding gene. In the bacterium Mycoplasma pneumoniae, ProTInSeq identifies 83% of its annotated proteins, along with 5 proteins and 153 small ORF-encoded proteins (SEPs; ≤100 aa) that were not previously annotated. Moreover, ProTInSeq can be utilized for detecting translational noise, as well as for relative quantification and transmembrane topology estimation of fitness and non-essential proteins. By integrating various identification approaches, the number of initially annotated SEPs in this bacterium increases from 27 to 329, with a quarter of them predicted to possess antimicrobial potential. Herein, we describe a methodology complementary to Ribo-Seq and mass spectroscopy that can identify SEPs while providing other insights in a proteome with a flexible and cost-effective DNA ultra-deep sequencing approach.

bacLIFE: a user-friendly computational workflow for genome analysis and prediction of lifestyle-associated genes in bacteria

Bacteria have an extensive adaptive ability to live in close association with eukaryotic hosts, exhibiting detrimental, neutral or beneficial effects on host growth and health. However, the genes involved in niche adaptation are mostly unknown and their functions poorly characterized. Here, we present bacLIFE ( https://github.com/Carrion-lab/bacLIFE ) a streamlined computational workflow for genome annotation, large-scale comparative genomics, and prediction of lifestyle-associated genes (LAGs). As a proof of concept, we analyzed 16,846 genomes from the Burkholderia/Paraburkholderia and Pseudomonas genera, which led to the identification of hundreds of genes potentially associated with a plant pathogenic lifestyle. Site-directed mutagenesis of 14 of these predicted LAGs of unknown function, followed by plant bioassays, showed that 6 predicted LAGs are indeed involved in the phytopathogenic lifestyle of Burkholderia plantarii and Pseudomonas syringae pv. phaseolicola. These 6 LAGs encompassed a glycosyltransferase, extracellular binding proteins, homoserine dehydrogenases and hypothetical proteins. Collectively, our results highlight bacLIFE as an effective computational tool for prediction of LAGs and the generation of hypotheses for a better understanding of bacteria-host interactions.

Transmission and dynamics of mother-infant gut viruses during pregnancy and early life

Early development of the gut ecosystem is crucial for lifelong health. While infant gut bacterial communities have been studied extensively, the infant gut virome remains under-explored. To study the development of the infant gut virome over time and the factors that shape it, we longitudinally assess the composition of gut viruses and their bacterial hosts in 30 women during and after pregnancy and in their 32 infants during their first year of life. Using shotgun metagenomic sequencing applied to dsDNA extracted from Virus-Like Particles (VLPs) and bacteria, we generate 205 VLP metaviromes and 322 total metagenomes. With this data, we show that while the maternal gut virome composition remains stable during late pregnancy and after birth, the infant gut virome is dynamic in the first year of life. Notably, infant gut viromes contain a higher abundance of active temperate phages compared to maternal gut viromes, which decreases over the first year of life. Moreover, we show that the feeding mode and place of delivery influence the gut virome composition of infants. Lastly, we provide evidence of co-transmission of viral and bacterial strains from mothers to infants, demonstrating that infants acquire some of their virome from their mother's gut.

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

FacZ is a GpsB-interacting protein that prevents aberrant division-site placement in Staphylococcus aureus

Staphylococcus aureus is a Gram-positive pathogen responsible for antibiotic-resistant infections. To identify vulnerabilities in cell envelope biogenesis that may overcome resistance, we enriched for S. aureus transposon mutants with defects in cell surface integrity or cell division by sorting for cells that stain with propidium iodide or have increased light-scattering properties, respectively. Transposon sequencing of the sorted populations identified more than 20 previously uncharacterized factors impacting these processes. Cells inactivated for one of these proteins, factor preventing extra Z-rings (FacZ, SAOUHSC_01855), showed aberrant membrane invaginations and multiple FtsZ cytokinetic rings. These phenotypes were suppressed in mutants lacking the conserved cell-division protein GpsB, which forms an interaction hub bridging envelope biogenesis factors with the cytokinetic ring in S. aureus. FacZ was found to interact directly with GpsB in vitro and in vivo. We therefore propose that FacZ is an envelope biogenesis factor that antagonizes GpsB function to prevent aberrant division events in S. aureus.

Synthesis of lipid-linked precursors of the bacterial cell wall is governed by a feedback control mechanism in Pseudomonas aeruginosa

Many bacterial surface glycans such as the peptidoglycan (PG) cell wall are built from monomeric units linked to a polyprenyl lipid carrier. How this limiting carrier is distributed among competing pathways has remained unclear. Here we describe the isolation of hyperactive variants of Pseudomonas aeruginosa MraY, the enzyme that forms the first lipid-linked PG precursor. These variants result in the elevated production of the final PG precursor lipid II in cells and are hyperactive in vitro. The activated MraY variants have substitutions that map to a cavity on the extracellular side of the dimer interface, far from the active site. Our structural and molecular dynamics results suggest that this cavity is a binding site for externalized lipid II. Overall, our results support a model in which excess externalized lipid II allosterically inhibits MraY, providing a feedback mechanism that prevents the sequestration of lipid carrier in the PG biogenesis pathway.

The Deinococcus protease PprI senses DNA damage by directly interacting with single-stranded DNA

Bacteria have evolved various response systems to adapt to environmental stress. A protease-based derepression mechanism in response to DNA damage was characterized in Deinococcus, which is controlled by the specific cleavage of repressor DdrO by metallopeptidase PprI (also called IrrE). Despite the efforts to document the biochemical, physiological, and downstream regulation of PprI-DdrO, the upstream regulatory signal activating this system remains unclear. Here, we show that single-stranded DNA physically interacts with PprI protease, which enhances the PprI-DdrO interactions as well as the DdrO cleavage in a length-dependent manner both in vivo and in vitro. Structures of PprI, in its apo and complexed forms with single-stranded DNA, reveal two DNA-binding interfaces shaping the cleavage site. Moreover, we show that the dynamic monomer-dimer equilibrium of PprI is also important for its cleavage activity. Our data provide evidence that single-stranded DNA could serve as the signal for DNA damage sensing in the metalloprotease/repressor system in bacteria. These results also shed light on the survival and acquired drug resistance of certain bacteria under antimicrobial stress through a SOS-independent pathway.

TcrXY is an acid-sensing two-component transcriptional regulator of Mycobacterium tuberculosis required for persistent infection

The ability of Mycobacterium tuberculosis (Mtb) to persist in the host complicates and prolongs tuberculosis (TB) patient chemotherapy. Here we demonstrate that a neglected two-component system (TCS) of Mtb, TcrXY, is an autoregulated acid-sensing TCS that controls a functionally diverse 70-gene regulon required for bacterial persistence. Characterisation of two representatives of this regulon, Rv3706c and Rv3705A, implicate these genes as key determinants for the survival of Mtb in vivo by serving as important effectors to mitigate redox stress at acidic pH. We show that genetic silencing of the response regulator tcrX using CRISPR interference attenuates the persistence of Mtb during chronic mouse infection and improves treatment with the two front-line anti-TB drugs, rifampicin and isoniazid. We propose that targeting TcrXY signal transduction blocks the ability of Mtb to sense and respond to acid stress, resulting in a disordered program of persistence to render the organism vulnerable to existing TB chemotherapy.

Commensal lifestyle regulated by a negative feedback loop between Arabidopsis ROS and the bacterial T2SS

Despite the plant health-promoting effects of plant microbiota, these assemblages also comprise potentially detrimental microbes. How plant immunity controls its microbiota to promote plant health under these conditions remains largely unknown. We find that commensal bacteria isolated from healthy Arabidopsis plants trigger diverse patterns of reactive oxygen species (ROS) production dependent on the immune receptors and completely on the NADPH oxidase RBOHD that selectively inhibited specific commensals, notably Xanthomonas L148. Through random mutagenesis, we find that L148 gspE, encoding a type II secretion system (T2SS) component, is required for the damaging effects of Xanthomonas L148 on rbohD mutant plants. In planta bacterial transcriptomics reveals that RBOHD suppresses most T2SS gene expression including gspE. L148 colonization protected plants against a bacterial pathogen, when gspE was inhibited by ROS or mutation. Thus, a negative feedback loop between Arabidopsis ROS and the bacterial T2SS tames a potentially detrimental leaf commensal and turns it into a microbe beneficial to the host.

σE of Streptomyces coelicolor can function both as a direct activator or repressor of transcription

σ factors are considered as positive regulators of gene expression. Here we reveal the opposite, inhibitory role of these proteins. We used a combination of molecular biology methods and computational modeling to analyze the regulatory activity of the extracytoplasmic σE factor from Streptomyces coelicolor. The direct activator/repressor function of σE was then explored by experimental analysis of selected promoter regions in vivo. Additionally, the σE interactome was defined. Taken together, the results characterize σE, its regulation, regulon, and suggest its direct inhibitory function (as a repressor) in gene expression, a phenomenon that may be common also to other σ factors and organisms.

Carbapenem-resistant Escherichia coli exhibit diverse spatiotemporal epidemiological characteristics across the globe

Carbapenem-resistant Escherichia coli (CREC) poses a severe global public health risk. This study reveals the worldwide geographic spreading patterns and spatiotemporal distribution characteristics of resistance genes in 7918 CREC isolates belonging to 497 sequence types (ST) and originating from 75 countries. In the last decade, there has been a transition in the prevailing STs from highly virulent ST131 and ST38 to higher antibiotic-resistant ST410 and ST167. The rise of multi-drug resistant strains of CREC carrying plasmids with extended-spectrum beta-lactamase (ESBL) resistance genes could be attributed to three important instances of host-switching events. The spread of CREC was associated with the changing trends in blaNDM-5, blaKPC-2, and blaOXA-48, as well as the plasmids IncFI, IncFII, and IncI. There were intercontinental geographic transfers of major CREC strains. Various crucial transmission hubs and patterns have been identified for ST131 in the United Kingdom, Italy, the United States, and China, ST167 in India, France, Egypt, and the United States, and ST410 in Thailand, Israel, the United Kingdom, France, and the United States. This work is valuable in managing CREC infections and preventing CREC occurrence and transmission inside healthcare settings and among diverse hosts.

CRISPR-Cas-based identification of a sialylated human milk oligosaccharides utilization cluster in the infant gut commensal Bacteroides dorei

The infant gut microbiome is impacted by early-life feeding, as human milk oligosaccharides (HMOs) found in breastmilk cannot be digested by infants and serve as nutrients for their gut bacteria. While the vast majority of HMO-utilization research has focused on Bifidobacterium species, recent studies have suggested additional HMO-utilizers, mostly Bacteroides, yet their utilization mechanism is poorly characterized. Here, we investigate Bacteroides dorei isolates from breastfed-infants and identify that polysaccharide utilization locus (PUL) 33 enables B. dorei to utilize sialylated HMOs. We perform transcriptional profiling and identity upregulated genes when growing on sialylated HMOs. Using CRISPR-Cas12 to knock-out four PUL33 genes, combined with complementation assays, we identify GH33 as the critical gene in PUL33 for sialylated HMO-utilization. This demonstration of an HMO-utilization system by Bacteroides species isolated from infants opens the way to further characterization of additional such systems, to better understand HMO-utilization in the infant gut.

The CRISPR effector Cam1 mediates membrane depolarization for phage defence

Prokaryotic type III CRISPR-Cas systems provide immunity against viruses and plasmids using CRISPR-associated Rossman fold (CARF) protein effectors1-5. Recognition of transcripts of these invaders with sequences that are complementary to CRISPR RNA guides leads to the production of cyclic oligoadenylate second messengers, which bind CARF domains and trigger the activity of an effector domain6,7. Whereas most effectors degrade host and invader nucleic acids, some are predicted to contain transmembrane helices without an enzymatic function. Whether and how these CARF-transmembrane helix fusion proteins facilitate the type III CRISPR-Cas immune response remains unknown. Here we investigate the role of cyclic oligoadenylate-activated membrane protein 1 (Cam1) during type III CRISPR immunity. Structural and biochemical analyses reveal that the CARF domains of a Cam1 dimer bind cyclic tetra-adenylate second messengers. In vivo, Cam1 localizes to the membrane, is predicted to form a tetrameric transmembrane pore, and provides defence against viral infection through the induction of membrane depolarization and growth arrest. These results reveal that CRISPR immunity does not always operate through the degradation of nucleic acids, but is instead mediated via a wider range of cellular responses.

The Helicobacter pylori Genome Project: insights into H. pylori population structure from analysis of a worldwide collection of complete genomes

Helicobacter pylori, a dominant member of the gastric microbiota, shares co-evolutionary history with humans. This has led to the development of genetically distinct H. pylori subpopulations associated with the geographic origin of the host and with differential gastric disease risk. Here, we provide insights into H. pylori population structure as a part of the Helicobacter pylori Genome Project (HpGP), a multi-disciplinary initiative aimed at elucidating H. pylori pathogenesis and identifying new therapeutic targets. We collected 1011 well-characterized clinical strains from 50 countries and generated high-quality genome sequences. We analysed core genome diversity and population structure of the HpGP dataset and 255 worldwide reference genomes to outline the ancestral contribution to Eurasian, African, and American populations. We found evidence of substantial contribution of population hpNorthAsia and subpopulation hspUral in Northern European H. pylori. The genomes of H. pylori isolated from northern and southern Indigenous Americans differed in that bacteria isolated in northern Indigenous communities were more similar to North Asian H. pylori while the southern had higher relatedness to hpEastAsia. Notably, we also found a highly clonal yet geographically dispersed North American subpopulation, which is negative for the cag pathogenicity island, and present in 7% of sequenced US genomes. We expect the HpGP dataset and the corresponding strains to become a major asset for H. pylori genomics.

Post-crotonylation oxidation by a monooxygenase promotes acetyl-CoA synthetase degradation in Streptomyces roseosporus

Protein post-translational modifications (PTMs) with various acyl groups play central roles in Streptomyces. But whether these acyl groups can be further modified, and the influences of these potential modifications on bacterial physiology have not been addressed. Here in Streptomyces roseosporus with rich crotonylation, a luciferase monooxygenase LimB is identified to elaborately regulate the crotonylation level, morphological development and antibiotic production by oxidation on the crotonyl groups of an acetyl-CoA synthetase Acs. This chemical modification on crotonylation leads to Acs degradation via the protease ClpP1/2 pathway and lowered intracellular crotonyl-CoA pool. Thus, we show that acyl groups after PTMs can be further modified, herein named post-PTM modification (PPM), and LimB is a PTM modifier to control the substrate protein turnover for cell development of Streptomyces. These findings expand our understanding of the complexity of chemical modifications on proteins for physiological regulation, and also suggest that PPM would be widespread.

Ongoing diversification of the global fish pathogen Piscirickettsia salmonis through genetic isolation and transposition bursts

The management of bacterial pathogens remains a key challenge of aquaculture. The marine gammaproteobacterium Piscirickettsia salmonis is the etiological agent of piscirickettsiosis and causes multi-systemic infections in different salmon species, resulting in considerable mortality and substantial commercial losses. Here, we elucidate its global diversity, evolution, and selection during human interventions. Our comprehensive analysis of 73 closed, high quality genome sequences covered strains from major outbreaks and was supplemented by an analysis of all P. salmonis 16S rRNA gene sequences and metagenomic reads available in public databases. Genome comparison showed that Piscirickettsia comprises at least three distinct, genetically isolated species of which two showed evidence for continuing speciation. However, at least twice the number of species exist in marine fish or seawater. A hallmark of Piscirickettsia diversification is the unprecedented amount and diversity of transposases which are particularly active in subgroups undergoing rapid speciation and are key to the acquisition of novel genes and to pseudogenization. Several group-specific genes are involved in surface antigen synthesis and may explain the differences in virulence between strains. However, the frequent failure of antibiotic treatment of piscirickettsiosis outbreaks cannot be explained by horizontal acquisition of resistance genes which so far occurred only very rarely. Besides revealing a dynamic diversification of an important pathogen, our study also provides the data for improving its surveillance, predicting the emergence of novel lineages, and adapting aquaculture management, and thereby contributes towards the sustainability of salmon farming.

Soil bacterium manipulates antifungal weapons by sensing intracellular type IVA secretion system effectors of a competitor

Soil beneficial bacteria can effectively inhibit bacterial pathogens by assembling contact-dependent killing weapons, such as the type IVA secretion system (T4ASS). It's not clear whether these antibacterial weapons are involved in biotrophic microbial interactions in soil. Here we showed that an antifungal antibiotic 2,4-DAPG production of the soil bacterium, Pseudomonas protegens can be triggered by another soil bacterium, Lysobacter enzymogenes, via T4ASS by co-culturing on agar plates to mimic cell-to-cell contact. We demonstrated that the induced 2,4-DAPG production of P. protegens is achieved by intracellular detection of the T4ASS effector protein Le1519 translocated from L. enzymogenes. We defined Le1519 as LtaE (Lysobacter T4E triggering antifungal effects), which specifically stimulates the expression of 2,4-DAPG biosynthesis genes in P. protegens, thereby protecting soybean seedlings from infection by the fungus Rhizoctonia solani. We further found that LtaE directly bound to PhlF, a pathway-specific transcriptional repressor of the 2,4-DAPG biosynthesis, then activated the 2,4-DAPG production. Our results highlight a novel pattern of microbial interspecies and interkingdom interactions, providing a unique case for expanding the diversity of soil microbial interactions.

Phylogenetic and ecophysiological novelty of subsurface mercury methylators in mangrove sediments

Mangrove sediment is a crucial component in the global mercury (Hg) cycling and acts as a hotspot for methylmercury (MeHg) production. Early evidence has documented the ubiquity of well-studied Hg methylators in mangrove superficial sediments; however, their diversity and metabolic adaptation in the more anoxic and highly reduced subsurface sediments are lacking. Through MeHg biogeochemical assay and metagenomic sequencing, we found that mangrove subsurface sediments (20-100 cm) showed a less hgcA gene abundance but higher diversity of Hg methylators than superficial sediments (0-20 cm). Regional-scale investigation of mangrove subsurface sediments spanning over 1500 km demonstrated a prevalence and family-level novelty of Hg-methylating microbial lineages (i.e., those affiliated to Anaerolineae, Phycisphaerae, and Desulfobacterales). We proposed the candidate phylum Zixibacteria lineage with sulfate-reducing capacity as a currently understudied Hg methylator across anoxic environments. Unlike other Hg methylators, the Zixibacteria lineage does not use the Wood-Ljungdahl pathway but has unique capabilities of performing methionine synthesis to donate methyl groups. The absence of cobalamin biosynthesis pathway suggests that this Hg-methylating lineage may depend on its syntrophic partners (i.e., Syntrophobacterales members) for energy in subsurface sediments. Our results expand the diversity of subsurface Hg methylators and uncover their unique ecophysiological adaptations in mangrove sediments.

When less is more: shortening the Lpp protein leads to increased vancomycin resistance in Escherichia coli

Vancomycin is a naturally occurring cell-wall-targeting glycopeptide antibiotic. Due to the low potency of this antibiotic against Gram-negative pathogens, such as Escherichia coli, there is a limited knowledge about interactions between vancomycin and this group of bacteria. Here, we show that an in-frame 63 bp deletion of the lpp gene caused a fourfold increase in vancomycin resistance in E. coli. The resulting protein, LppΔ21, is 21 amino acids shorter than the wild-type Lpp, a helical structural lipoprotein that controls the width of the periplasmic space through its length. The mutant remains susceptible to synergistic growth inhibition by combination of furazolidone and vancomycin; with furazolidone decreasing the vancomycin MIC by eightfold. These findings have clinical relevance, given that the vancomycin concentration required to select the lpp mutation is reachable during typical vancomycin oral administration for treating Clostridioides difficile infections. Combination therapy with furazolidone, however, is likely to prevent emergence and outgrowth of the lpp-mutated Gram-negative coliforms, avoiding exacerbation of the patient's condition during the treatment.

Strain-level diversity of symbiont communities between individuals and populations of a bioluminescent fish

The bioluminescent symbiosis involving the urchin cardinalfish, Siphamia tubifer, and Photobacterium mandapamensis, a luminous member of the Vibrionaceae, is highly specific compared to other bioluminescent fish-bacteria associations. Despite this high degree of specificity, patterns of genetic diversity have been observed for the symbionts from hosts sampled over relatively small spatial scales. We characterized and compared sub-species, strain-level symbiont diversity within and between S. tubifer hosts sampled from the Philippines and Japan using PCR fingerprinting. We then carried out whole genome sequencing of the unique symbiont genotypes identified to characterize the genetic diversity of the symbiont community and the symbiont pangenome. We determined that an individual light organ contains six symbiont genotypes on average, but varied between 1-13. Additionally, we found that there were few genotypes shared between hosts from the same location. A phylogenetic analysis of the unique symbiont strains indicated location-specific clades, suggesting some genetic differentiation in the symbionts between host populations. We also identified symbiont genes that were variable between strains, including luxF, a member of the lux operon, which is responsible for light production. We quantified the light emission and growth rate of two strains missing luxF along with the other strains isolated from the same light organs and determined that strains lacking luxF were dimmer but grew faster than most of the other strains, suggesting a potential metabolic trade-off. This study highlights the importance of strain-level diversity in microbial associations and provides new insight into the underlying genetic architecture of intraspecific symbiont communities within a host.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

Diurnal cycles drive rhythmic physiology and promote survival in facultative phototrophic bacteria

Bacteria have evolved many strategies to spare energy when nutrients become scarce. One widespread such strategy is facultative phototrophy, which helps heterotrophs supplement their energy supply using light. Our knowledge of the impact that such behaviors have on bacterial fitness and physiology is, however, still limited. Here, we study how a representative of the genus Porphyrobacter, in which aerobic anoxygenic phototrophy is ancestral, responds to different light regimes under nutrient limitation. We show that bacterial survival in stationary phase relies on functional reaction centers and varies depending on the light regime. Under dark-light alternance, our bacterial model presents a diphasic life history dependent on phototrophy: during dark phases, the cells inhibit DNA replication and part of the population lyses and releases nutrients, while subsequent light phases allow for the recovery and renewed growth of the surviving cells. We correlate these cyclic variations with a pervasive pattern of rhythmic transcription which reflects global changes in diurnal metabolic activity. Finally, we demonstrate that, compared to either a phototrophy mutant or a bacteriochlorophyll a overproducer, the wild type strain is better adapted to natural environments, where regular dark-light cycles are interspersed with additional accidental dark episodes. Overall, our results highlight the importance of light-induced biological rhythms in a new model of aerobic anoxygenic phototroph representative of an ecologically important group of environmental bacteria.

Global pathogenomic analysis identifies known and candidate genetic antimicrobial resistance determinants in twelve species

Surveillance programs for managing antimicrobial resistance (AMR) have yielded thousands of genomes suited for data-driven mechanism discovery. We present a workflow integrating pangenomics, gene annotation, and machine learning to identify AMR genes at scale. When applied to 12 species, 27,155 genomes, and 69 drugs, we 1) find AMR gene transfer mostly confined within related species, with 925 genes in multiple species but just eight in multiple phylogenetic classes, 2) demonstrate that discovery-oriented support vector machines outperform contemporary methods at recovering known AMR genes, recovering 263 genes compared to 145 by Pyseer, and 3) identify 142 AMR gene candidates. Validation of two candidates in E. coli BW25113 reveals cases of conditional resistance: ΔcycA confers ciprofloxacin resistance in minimal media with D-serine, and frdD V111D confers ampicillin resistance in the presence of ampC by modifying the overlapping promoter. We expect this approach to be adaptable to other species and phenotypes.

Distributed genotyping and clustering of Neisseria strains reveal continual emergence of epidemic meningococcus over a century

Core genome multilocus sequence typing (cgMLST) is commonly used to classify bacterial strains into different types, for taxonomical and epidemiological applications. However, cgMLST schemes require central databases for the nomenclature of new alleles and sequence types, which must be synchronized worldwide and involve increasingly intensive calculation and storage demands. Here, we describe a distributed cgMLST (dcgMLST) scheme that does not require a central database of allelic sequences and apply it to study evolutionary patterns of epidemic and endemic strains of the genus Neisseria. We classify 69,994 worldwide Neisseria strains into multi-level clusters that assign species, lineages, and local disease outbreaks. We divide Neisseria meningitidis into 168 endemic lineages and three epidemic lineages responsible for at least 9 epidemics in the past century. According to our analyses, the epidemic and endemic lineages experienced very different population dynamics in the past 100 years. Epidemic lineages repetitively emerged from endemic lineages, disseminated worldwide, and apparently disappeared rapidly afterward. We propose a stepwise model for the evolutionary trajectory of epidemic lineages in Neisseria, and expect that the development of similar dcgMLST schemes will facilitate epidemiological studies of other bacterial pathogens.

Depletion of tRNA CCA-adding enzyme in Mycobacterium tuberculosis leads to polyadenylation of transcripts and precursor tRNAs

In reference to gene annotation, more than half of the tRNA species synthesized by Mycobacterium tuberculosis require the enzymatic addition of the cytosine-cytosine-adenine (CCA) tail, which is indispensable for amino acid charging and tRNA functionality. It makes the mycobacterial CCA-adding enzyme essential for survival of the bacterium and a potential target for novel pipelines in drug discovery avenues. Here, we described the rv3907c gene product, originally annotated as poly(A)polymerase (rv3907c, PcnA) as a functional CCA-adding enzyme (CCAMtb) essential for viability of M. tuberculosis. The depletion of the enzyme affected tRNAs maturation, inhibited bacilli growth, and resulted in abundant accumulation of polyadenylated RNAs. We determined the enzymatic activities displayed by the mycobacterial CCAMtb in vitro and studied the effects of inhibiting of its transcription in bacterial cells. We are the first to properly confirm the existence of RNA polyadenylation in mycobacteria, a previously controversial phenomenon, which we found promoted upon CCA-adding enzyme downexpression.

The genomic epidemiology of shigellosis in South Africa

Shigellosis, a leading cause of diarrhoeal mortality and morbidity globally, predominantly affects children under five years of age living in low- and middle-income countries. While whole genome sequence analysis (WGSA) has been effectively used to further our understanding of shigellosis epidemiology, antimicrobial resistance, and transmission, it has been under-utilised in sub-Saharan Africa. In this study, we applied WGSA to large sub-sample of surveillance isolates from South Africa, collected from 2011 to 2015, focussing on Shigella flexneri 2a and Shigella sonnei. We find each serotype is epidemiologically distinct. The four identified S. flexneri 2a clusters having distinct geographical distributions, and antimicrobial resistance (AMR) and virulence profiles, while the four sub-Clades of S. sonnei varied in virulence plasmid retention. Our results support serotype specific lifestyles as a driver for epidemiological differences, show AMR is not required for epidemiological success in S. flexneri, and that the HIV epidemic may have promoted Shigella population expansion.

Regulation of the physiology and virulence of Ralstonia solanacearum by the second messenger 2',3'-cyclic guanosine monophosphate

Previous studies have demonstrated that bis-(3',5')-cyclic diguanosine monophosphate (bis-3',5'-c-di-GMP) is a ubiquitous second messenger employed by bacteria. Here, we report that 2',3'-cyclic guanosine monophosphate (2',3'-cGMP) controls the important biological functions, quorum sensing (QS) signaling systems and virulence in Ralstonia solanacearum through the transcriptional regulator RSp0980. This signal specifically binds to RSp0980 with high affinity and thus abolishes the interaction between RSp0980 and the promoters of target genes. In-frame deletion of RSp0334, which contains an evolved GGDEF domain with a LLARLGGDQF motif required to catalyze 2',3'-cGMP to (2',5')(3',5')-cyclic diguanosine monophosphate (2',3'-c-di-GMP), altered the abovementioned important phenotypes through increasing the intracellular 2',3'-cGMP levels. Furthermore, we found that 2',3'-cGMP, its receptor and the evolved GGDEF domain with a LLARLGGDEF motif also exist in the human pathogen Salmonella typhimurium. Together, our work provides insights into the unusual function of the GGDEF domain of RSp0334 and the special regulatory mechanism of 2',3'-cGMP signal in bacteria.

High-resolution temporal profiling of E. coli transcriptional response

Understanding how cells dynamically adapt to their environment is a primary focus of biology research. Temporal information about cellular behavior is often limited by both small numbers of data time-points and the methods used to analyze this data. Here, we apply unsupervised machine learning to a data set containing the activity of 1805 native promoters in E. coli measured every 10 minutes in a high-throughput microfluidic device via fluorescence time-lapse microscopy. Specifically, this data set reveals E. coli transcriptome dynamics when exposed to different heavy metal ions. We use a bioinformatics pipeline based on Independent Component Analysis (ICA) to generate insights and hypotheses from this data. We discovered three primary, time-dependent stages of promoter activation to heavy metal stress (fast, intermediate, and steady). Furthermore, we uncovered a global strategy E. coli uses to reallocate resources from stress-related promoters to growth-related promoters following exposure to heavy metal stress.

A smooth tubercle bacillus from Ethiopia phylogenetically close to the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) includes several human- and animal-adapted pathogens. It is thought to have originated in East Africa from a recombinogenic Mycobacterium canettii-like ancestral pool. Here, we describe the discovery of a clinical tuberculosis strain isolated in Ethiopia that shares archetypal phenotypic and genomic features of M. canettii strains, but represents a phylogenetic branch much closer to the MTBC clade than to the M. canettii strains. Analysis of genomic traces of horizontal gene transfer in this isolate and previously identified M. canettii strains indicates a persistent albeit decreased recombinogenic lifestyle near the emergence of the MTBC. Our findings support that the MTBC emergence from its putative free-living M. canettii-like progenitor is evolutionarily very recent, and suggest the existence of a continuum of further extant derivatives from ancestral stages, close to the root of the MTBC, along the Great Rift Valley.

The metabolic, virulence and antimicrobial resistance profiles of colonising Streptococcus pneumoniae shift after PCV13 introduction in urban Malawi

Streptococcus pneumoniae causes substantial mortality among children under 5-years-old worldwide. Polysaccharide conjugate vaccines (PCVs) are highly effective at reducing vaccine serotype disease, but emergence of non-vaccine serotypes and persistent nasopharyngeal carriage threaten this success. We investigated the hypothesis that following vaccine, adapted pneumococcal genotypes emerge with the potential for vaccine escape. We genome sequenced 2804 penumococcal isolates, collected 4-8 years after introduction of PCV13 in Blantyre, Malawi. We developed a pipeline to cluster the pneumococcal population based on metabolic core genes into "Metabolic genotypes" (MTs). We show that S. pneumoniae population genetics are characterised by emergence of MTs with distinct virulence and antimicrobial resistance (AMR) profiles. Preliminary in vitro and murine experiments revealed that representative isolates from emerging MTs differed in growth, haemolytic, epithelial infection, and murine colonisation characteristics. Our results suggest that in the context of PCV13 introduction, pneumococcal population dynamics had shifted, a phenomenon that could further undermine vaccine control and promote spread of AMR.

Interaction networks of Escherichia coli replication proteins under different bacterial growth conditions

In this work we analyzed protein-protein interactions (PPIs) formed by E. coli replication proteins under three disparate bacterial growth conditions. The chosen conditions corresponded to fast exponential growth, slow exponential growth and growth cessation at the stationary phase. We performed affinity purification coupled with mass spectrometry (AP-MS) of chromosomally expressed proteins (DnaA, DnaB, Hda, SeqA, DiaA, DnaG, HolD, NrdB), tagged with sequential peptide affinity (SPA) tag. Composition of protein complexes was characterized using MaxQuant software. To filter out unspecific interactions, we employed double negative control system and we proposed qualitative and quantitative data analysis strategies that can facilitate hits identification in other AP-MS datasets. Our motivation to undertake this task was still insufficient understanding of molecular mechanisms coupling DNA replication to cellular growth. Previous works suggested that such control mechanisms could involve physical interactions of replication factors with metabolic or cell envelope proteins. However, the dynamic replication protein interaction network (PIN) obtained in this study can be used to characterize links between DNA replication and various cellular processes in other contexts.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

Construction of Streptomyces coelicolor A3(2) mutants that exclusively produce NA4/NA6 intermediates of agarose metabolism through mutation induction

NA4/NA6, an intermediate degradation product of β-agarase, is a high value-added product with anticancer, anti-obesity, and anti-diabetic effects. Therefore, a method that enables the efficient production of NA4/NA6 would be useful from economic and medical perspectives. In this study, we aimed to generate a Streptomyces coelicolor A3(2) mutant M22-2C43 that produces NA4/NA6 as a final product; this method serves as a more efficient alternative to the enzymatic conversion of β-agarase for the generation of these products. The M22-2C43 strain was generated through two rounds of mutagenesis and screening for increased β-agarase activity and effective production of NA4/NA6. We assembled the complete genomes of two mutants, M22 and M22-2C43, which were identified following a two-round screening. Large and small genetic changes were found in these two mutants, including the loss of two plasmids present in wild-type S. coelicolor A3(2) and chromosome circularization of mutant M22-2C43. These findings suggest that mutant M22-2C43 can produce NA4/NA6 as a degradation product due to functional inactivation of the dagB gene through a point mutation (G474A), ultimately preventing further degradation of NA4/NA6 to NA2. To our knowledge, this is the first report of a microbial strain that can effectively produce NA4/NA6 as the main degradation product of β-agarase, opening the door for the use of this species for the large-scale production of this valuable product.

Impact of transient acquired hypermutability on the inter- and intra-species competitiveness of Pseudomonas aeruginosa

Once acquired, hypermutation is unrelenting, and in the long-term, leads to impaired fitness due to its cumulative impact on the genome. This raises the question of why hypermutators arise so frequently in microbial ecosystems. In this work, we explore this problem by examining how the transient acquisition of hypermutability affects inter- and intra-species competitiveness, and the response to environmental insults such as antibiotic challenge. We do this by engineering Pseudomonas aeruginosa to allow the expression of an important mismatch repair gene, mutS, to be experimentally controlled over a wide dynamic range. We show that high levels of mutS expression induce genomic stasis (hypomutation), whereas lower levels of induction lead to progressively higher rates of mutation. Whole-genome sequence analyses confirmed that the mutational spectrum of the inducible hypermutator is similar to the distinctive profile associated with mutS mutants obtained from the airways of people with cystic fibrosis (CF). The acquisition of hypermutability conferred a distinct temporal fitness advantage over the wild-type P. aeruginosa progenitor strain, in both the presence and the absence of an antibiotic selection pressure. However, over a similar time-scale, acquisition of hypermutability had little impact on the population dynamics of P. aeruginosa when grown in the presence of a competing species (Staphylococcus aureus). These data indicate that in the short term, acquired hypermutability primarily confers a competitive intra-species fitness advantage.

Hydrogenotrophic methanogenesis is the key process in the obligately syntrophic consortium of the anaerobic ameba Pelomyxa schiedti

Pelomyxa is a genus of anaerobic amoebae that live in consortia with multiple prokaryotic endosymbionts. Although the symbionts represent a large fraction of the cellular biomass, their metabolic roles have not been investigated. Using single-cell genomics and transcriptomics, we have characterized the prokaryotic community associated with P. schiedti, which is composed of two bacteria, Candidatus Syntrophus pelomyxae (class Deltaproteobacteria) and Candidatus Vesiculincola pelomyxae (class Clostridia), and a methanogen, Candidatus Methanoregula pelomyxae. Fluorescence in situ hybridization and electron microscopy showed that Ca. Vesiculincola pelomyxae is localized inside vesicles, whereas the other endosymbionts occur freely in the cytosol, with Ca. Methanoregula pelomyxae enriched around the nucleus. Genome and transcriptome-based reconstructions of the metabolism suggests that the cellulolytic activity of P. schiedti produces simple sugars that fuel its own metabolism and the metabolism of a Ca. Vesiculincola pelomyxae, while Ca. Syntrophus pelomyxae energy metabolism relies on degradation of butyrate and isovalerate from the environment. Both species of bacteria and the ameba use hydrogenases to transfer the electrons from reduced equivalents to hydrogen, a process that requires a low hydrogen partial pressure. This is achieved by the third endosymbiont, Ca. Methanoregula pelomyxae, which consumes H2 and formate for methanogenesis. While the bacterial symbionts can be successfully eliminated by vancomycin treatment without affecting the viability of the amoebae, treatment with 2-bromoethanesulfonate, a specific inhibitor of methanogenesis, killed the amoebae, indicating the essentiality of the methanogenesis for this consortium.

Inferring bacterial transmission dynamics using deep sequencing genomic surveillance data

Identifying and interrupting transmission chains is important for controlling infectious diseases. One way to identify transmission pairs - two hosts in which infection was transmitted from one to the other - is using the variation of the pathogen within each single host (within-host variation). However, the role of such variation in transmission is understudied due to a lack of experimental and clinical datasets that capture pathogen diversity in both donor and recipient hosts. In this work, we assess the utility of deep-sequenced genomic surveillance (where genomic regions are sequenced hundreds to thousands of times) using a mouse transmission model involving controlled spread of the pathogenic bacterium Citrobacter rodentium from infected to naïve female animals. We observe that within-host single nucleotide variants (iSNVs) are maintained over multiple transmission steps and present a model for inferring the likelihood that a given pair of sequenced samples are linked by transmission. In this work we show that, beyond the presence and absence of within-host variants, differences arising in the relative abundance of iSNVs (allelic frequency) can infer transmission pairs more precisely. Our approach further highlights the critical role bottlenecks play in reserving the within-host diversity during transmission.