Analysis of Bottlenecks in Experimental Models of Infection

Within-plant genetic drift to control virus adaptation to host resistance genes

Manipulating evolutionary forces imposed by hosts on pathogens like genetic drift and selection could avoid the emergence of virulent pathogens. For instance, increasing genetic drift could decrease the risk of pathogen adaptation through the random fixation of deleterious mutations or the elimination of favorable ones in the pathogen population. However, no experimental proof of this approach is available for a plant-pathogen system. We studied the impact of pepper (Capsicum annuum) lines carrying the same major resistance gene but contrasted genetic backgrounds on the evolution of Potato virus Y (PVY). The pepper lines were chosen for the contrasted levels of genetic drift (inversely related to Ne, the effective population size) they exert on PVY populations, as well as for their contrasted resistance efficiency (inversely related to the initial replicative fitness, Wi, of PVY in these lines). Experimental evolution was performed by serially passaging 64 PVY populations every month on six contrasted pepper lines during seven months. These PVY populations exhibited highly divergent evolutionary trajectories, ranging from viral extinctions to replicative fitness gains. The sequencing of the PVY VPg cistron, where adaptive mutations are likely to occur, allowed linking these replicative fitness gains to parallel adaptive nonsynonymous mutations. Evolutionary trajectories were well explained by the genetic drift imposed by the host. More specifically, Ne, Wi and their synergistic interaction played a major role in the fate of PVY populations. When Ne was low (i.e. strong genetic drift), the final PVY replicative fitness remained close to the initial replicative fitness, whereas when Ne was high (i.e. low genetic drift), the final PVY replicative fitness was high independently of the replicative fitness of the initially inoculated virus. We show that combining a high resistance efficiency (low Wi) and a strong genetic drift (low Ne) is the best solution to increase resistance durability, that is, to avoid virus adaptation on the long term.

Shedding light on bacteria-host interactions with the aid of TnSeq approaches

Bacteria are highly adaptable and grow in diverse niches, where they often interact with eukaryotic organisms. These interactions with different hosts span the entire spectrum from symbiosis to pathogenicity and thus determine the lifestyle of the bacterium. Knowledge of the genetic determinants involved in animal and plant host colonization by pathogenic and mutualistic bacteria is not only crucial to discover new drug targets for disease management but also for developing novel biostimulant strategies. In the last decades, significant progress in genome-wide high-throughput technologies such as transposon insertion sequencing has led to the identification of pathways that enable efficient host colonization. However, the extent to which similar genes play a role in this process in different bacteria is yet unclear. This review highlights the commonalities and specificities of bacterial determinants important for bacteria-host interaction.

Temperature dependence of the mutation rate towards antibiotic resistance

Objectives

Environmental conditions can influence mutation rates in bacteria. Fever is a common response to infection that alters the growth conditions of infecting bacteria. Here we examine how a temperature change, such as is associated with fever, affects the mutation rate towards antibiotic resistance.

Methods

We used a fluctuation test to assess the mutation rate towards antibiotic resistance in Escherichia coli at two different temperatures: 37°C (normal temperature) and 40°C (fever temperature). We performed this measurement for three different antibiotics with different modes of action: ciprofloxacin, rifampicin and ampicillin.

Results

In all cases, the mutation rate towards antibiotic resistance turned out to be temperature dependent, but in different ways. Fever temperatures led to a reduced mutation rate towards ampicillin resistance and an elevated mutation rate towards ciprofloxacin and rifampicin resistance.

Conclusions

This study shows that the mutation rate towards antibiotic resistance is impacted by a small change in temperature, such as associated with fever. This opens a new avenue to mitigate the emergence of antibiotic resistance by coordinating the choice of an antibiotic with the decision of whether or not to suppress fever when treating a patient. Hence, optimized combinations of antibiotics and fever suppression strategies may be a new weapon in the battle against antibiotic resistance.

Inducible transposon mutagenesis for genome-scale forward genetics

Transposon insertion sequencing (Tn-seq) is a powerful method for genome-scale functional genetics in bacteria. However, its effectiveness is often limited by a lack of mutant diversity, caused by either inefficient transposon delivery or stochastic loss of mutants due to population bottlenecks. Here, we introduce "InducTn-seq", which leverages inducible mutagenesis for temporal control of transposition. InducTn-seq generates millions of transposon mutants from a single colony, enabling the sensitive detection of subtle fitness defects and transforming binary classifications of gene essentiality into a quantitative fitness measurement across both essential and non-essential genes. Using a mouse model of infectious colitis, we show that InducTn-seq bypasses a highly restrictive host bottleneck to generate a diverse transposon mutant population from the few cells that initiate infection, revealing the role of oxygen-related metabolic plasticity in pathogenesis. Overall, InducTn-seq overcomes the limitations of traditional Tn-seq, unlocking new possibilities for genome-scale forward genetic screens in bacteria.

Interferons and tuft cell numbers are bottlenecks for persistent murine norovirus infection

Noroviruses (NoVs) are a leading cause of viral gastroenteritis. Despite global clinical relevance, our understanding of how host factors, such as antiviral cytokines interferons (IFNs), modulate NoV population dynamics is limited. Murine NoV (MNoV) is a tractable in vivo model for the study of host regulation of NoV. A persistent strain of MNoV, CR6, establishes a reservoir in intestinal tuft cells for chronic viral shedding in stool. However, the influence of host innate immunity and permissive cell numbers on viral population dynamics is an open question. We generated a pool of 20 different barcoded viruses (CR6BC) by inserting 6-nucleotide barcodes at the 3' position of the NS4 gene and used this pool as our viral inoculum for in vivo infections of different mouse lines. We found that over the course of persistent CR6 infection, shed virus was predominantly colon-derived, and viral barcode richness decreased over time irrespective of host immune status, suggesting that persistent infection involves a series of reinfection events. In mice lacking the IFN-λ receptor, intestinal barcode richness was enhanced, correlating with increased viral intestinal replication. IL-4 treatment, which increases tuft cell numbers, also increased barcode richness, indicating the abundance of permissive tuft cells to be a bottleneck during CR6 infection. In mice lacking type I IFN signaling (Ifnar1-/-) or all IFN signaling (Stat1-/-), barcode diversity at extraintestinal sites was dramatically increased, implicating different IFNs as critical bottlenecks at specific tissue sites. Of interest, extraintestinal barcodes were overlapping but distinct from intestinal barcodes, indicating that disseminated virus represents a distinct viral population than that replicating in the intestine. Barcoded viruses are a valuable tool to explore the influence of host factors on viral diversity in the context of establishment and maintenance of infection as well as dissemination and have provided important insights into how NoV infection proceeds in immunocompetent and immunocompromised hosts.

Transmission Bottleneck Size Estimation from De Novo Viral Genetic Variation

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, these approaches have the potential to substantially underestimate true transmission bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arise de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these 2 respiratory viruses.

Effects of plant tissue permeability on invasion and population bottlenecks of a phytopathogen

Pathogen genetic diversity varies in response to environmental changes. However, it remains unclear whether plant barriers to invasion could be considered a genetic bottleneck for phytopathogen populations. Here, we implement a barcoding approach to generate a pool of 90 isogenic and individually barcoded Ralstonia solanacearum strains. We used 90 of these strains to inoculate tomato plants with different degrees of physical permeability to invasion (intact roots, wounded roots and xylem inoculation) and quantify the phytopathogen population dynamics during invasion. Our results reveal that the permeability of plant roots impacts the degree of population bottleneck, genetic diversity, and composition of Ralstonia populations. We also find that selection is the main driver structuring pathogen populations when barriers to infection are less permeable, i.e., intact roots, the removal of root physical and immune barriers results in the predominance of stochasticity in population assembly. Taken together, our study suggests that plant root permeability constitutes a bottleneck for phytopathogen invasion and genetic diversity.

CRISPR Interference-Based Functional Small RNA Genomics

Small RNAs (sRNAs) are versatile regulators universally present in species across the prokaryotic kingdom, yet their functional characterization remains a major bottleneck. Gene inactivation through random transposon insertion has proven extremely valuable in discovering hidden gene functions. However, this approach is biased toward long genes and usually results in the underrepresentation of sRNA mutants. In contrast, CRISPR interference (CRISPRi) harnesses guide RNAs to recruit cleavage-deficient Cas nucleases to specific DNA loci. The ensuing steric hindrance inhibits RNA polymerase assembly at-or migration along-predefined genes, allowing for targeted knockdown screens without major length bias. In this chapter, we provide a detailed protocol for CRISPRi-based functional screening of bacterial sRNAs. Using the abundant microbiota species Bacteroides thetaiotaomicron as a model, we describe the design and generation of a guide library targeting the full intergenic sRNA repertoire of this organism and its application to identify sRNA knockdown-associated fitness effects. Our protocol is generic and thus suitable for the systematic assessment of sRNA-associated phenotypes in a wide range of bacterial species and experimental conditions. We expect CRISPRi-based functional genomics to boost sRNA research in understudied bacterial taxa, for instance, members of the gut microbiota.

How colonization bottlenecks, tissue niches, and transmission strategies shape protozoan infections

Protozoan pathogens such as Plasmodium spp., Leishmania spp., Toxoplasma gondii, and Trypanosoma spp. are often associated with high-mortality, acute and chronic diseases of global health concern. For transmission and immune evasion, protozoans have evolved diverse strategies to interact with a range of host tissue environments. These interactions are linked to disease pathology, yet our understanding of the association between parasite colonization and host homeostatic disruption is limited. Recently developed techniques for cellular barcoding have the potential to uncover the biology regulating parasite transmission, dissemination, and the stability of infection. Understanding bottlenecks to infection and the in vivo tissue niches that facilitate chronic infection and spread has the potential to reveal new aspects of parasite biology.

Organ agar serves as physiologically relevant alternative for in vivo bacterial colonization

Animal models for host-microbial interactions have proven valuable, yielding physiologically relevant data that may be otherwise difficult to obtain. Unfortunately, such models are lacking or nonexistent for many microbes. Here, we introduce organ agar, a straightforward method to enable the screening of large mutant libraries while avoiding physiological bottlenecks. We demonstrate that growth defects on organ agar were translatable to bacterial colonization deficiencies in a murine model. Specifically, we present a urinary tract infection agar model to interrogate an ordered library of Proteus mirabilis transposon mutants, with accurate prediction of bacterial genes critical for host colonization. Thus, we demonstrate the ability of ex vivo organ agar to reproduce in vivo deficiencies. Organ agar was also useful for identifying previously unknown links between biosynthetic genes and swarming motility. This work provides a readily adoptable technique that is economical and uses substantially fewer animals. We anticipate this method will be useful for a wide variety of microorganisms, both pathogenic and commensal, in a diverse range of model host species.

Construction and application of the conditionally essential gene knockdown library in Klebsiella pneumoniae to screen potential antimicrobial targets and virulence genes via Mobile-CRISPRi-seq

Klebsiella pneumoniae is a ubiquitous human pathogen, and its clinical treatment faces two major challenges: multidrug resistance and the pathogenesis of hypervirulent K. pneumoniae. The discovery and study of conditionally essential (CE) genes that can function as potential antimicrobial targets has always been a research concern due to their restriction in the development of novel antibiotics. However, the lack of essential functional genomic data has hampered the study of the mechanisms of essential genes related to antimicrobial susceptibility. In this study, we developed a pooled CE genes mobile clustered regularly interspaced short palindromic repeat (CRISPR) interference screening method (Mobile-CRISPRi-seq) for K. pneumoniae to identify genes that play critical roles in antimicrobial fitness in vitro and host immunity in vivo. Targeting 870 predicted CE genes in K. pneumoniae, Mobile-CRISPRi-seq uncovered the depletion of tetrahydrofolate synthesis pathway genes folB and folP under trimethoprim pressure. Our screening also identified genes waaE and fldA related to polymyxin and β-lactam susceptibility by applying a screening strategy based on Mobile-CRISPRi-seq and comparative genomics. Furthermore, using a mouse infection model and Mobile-CRISPRi-seq, multiple virulence genes were identified, and among these genes, pal, yciS, and ribB were demonstrated to contribute to the pathogenesis of K. pneumoniae. This study provides a simple, rapid, and effective platform for screening potential antimicrobial targets and virulence genes in K. pneumoniae, and this broadly applicable system can be expanded for high-throughput functional gene study in multiple pathogenic bacteria, especially in gram-negative bacteria. IMPORTANCE The discovery and investigation of conditionally essential (CE) genes that can function as potential antimicrobial targets has always been a research concern because of the restriction of antimicrobial targets in the development of novel antibiotics. In this study, we developed a pooled CE gene-wide mobile clustered regularly interspaced short palindromic repeat (CRISPR) interference sequencing (Mobile-CRISPRi-seq) strategy in Klebsiella pneumoniae to identify genes that play critical roles in the fitness of antimicrobials in vitro and host immunity in vivo. The data suggest a robust tool to screen for loss-of-function phenotypes in a pooled gene knockdown library in K. pneumoniae, and Mobile-CRISPRi-seq may be expanded to multiple bacteria for screening and identification of genes with crucial roles in the fitness of antimicrobials and hosts.

Promotion of plasmid maintenance by heterogeneous partitioning of microbial communities

Transferable plasmids play a critical role in shaping the functions of microbial communities. Previous studies suggested multiple mechanisms underlying plasmid persistence and abundance. Here, we focus on the interplay between heterogeneous community partitioning and plasmid fates. Natural microbiomes often experience partitioning that creates heterogeneous local communities with reduced population sizes and biodiversity. Little is known about how population partitioning affects the plasmid fate through the modulation of community structure. By modeling and experiments, we show that heterogeneous community partitioning can paradoxically promote the persistence of a plasmid that would otherwise not persist in a global community. Among the local communities created by partitioning, a minority will primarily consist of members able to transfer the plasmid fast enough to support its maintenance by serving as a local plasmid haven. Our results provide insights into plasmid maintenance and suggest a generalizable approach to modulate plasmid persistence for engineering and medical applications.

Innate immune responses yield tissue-specific bottlenecks that scale with pathogen dose

To cause infection, pathogens must overcome bottlenecks imposed by the host immune system. These bottlenecks restrict the inoculum and largely determine whether pathogen exposure results in disease. Infection bottlenecks therefore quantify the effectiveness of immune barriers. Here, using a model of Escherichia coli systemic infection, we identify bottlenecks that tighten or widen with higher inoculum sizes, revealing that the efficacy of innate immune responses can increase or decrease with pathogen dose. We term this concept "dose scaling". During E. coli systemic infection, dose scaling is tissue specific, dependent on the lipopolysaccharide (LPS) receptor TLR4, and can be recapitulated by mimicking high doses with killed bacteria. Scaling therefore depends on sensing of pathogen molecules rather than interactions between the host and live bacteria. We propose that dose scaling quantitatively links innate immunity with infection bottlenecks and is a valuable framework for understanding how the inoculum size governs the outcome of pathogen exposure.

Transmission bottleneck size estimation from de novo viral genetic variation

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, allele frequencies can change dramatically over the course of an individual's infection, such that sites that are polymorphic in the donor at the time of transmission may not be polymorphic in the donor at the time of sampling and allele frequencies at donor-polymorphic sites may change dramatically over the course of a recipient's infection. Because of this, transmission bottleneck sizes estimated using allele frequencies observed at a donor's polymorphic sites may be considerable underestimates of true bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arose de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these two respiratory viruses, using an approach that does not tend to underestimate transmission bottleneck sizes.

Neisseria gonorrhoeae MtrCDE Efflux Pump During In Vivo Experimental Genital Tract Infection in Men and Mice Reveals the Presence of Within-Host Colonization Bottleneck

The MtrCDE efflux pump of Neisseria gonorrhoeae exports a wide range of antimicrobial compounds that the gonococcus encounters at mucosal surfaces during colonization and infection. Here, we evaluate the role of this efflux pump system in strain FA1090 in human male urethral infection with a Controlled Human Infection Model. Using the strategy of competitive multi-strain infection with wild-type FA1090 and an isogenic mutant strain that does not contain a functional MtrCDE pump, we found that the presence of the efflux pump during human experimental infection did not confer a competitive advantage. This finding is in contrast to previous findings in female mice, which demonstrated that gonococci of strain FA19 lacking a functional MtrCDE pump had a significantly reduced fitness compared to the wild type strain in the lower genital tract of female mice. We conducted competitive infections in female mice with FA19 and FA1090 strains, including mutants that do not assemble a functional Mtr efflux pump, demonstrating the fitness advantage provided byt the MtrCDE efflux pump during infection of mice is strain dependent. Our data indicate that new gonorrhea treatment strategies targeting the MtrCDE efflux pump functions may not be universally efficacious in naturally occurring infections. Owing to the equal fitness of FA1090 strains in men, our experiments unexpectedly demonstrated the likely presence of an early colonization bottleneck of N. gonorrhoeae in the human male urethra.

TRIAL REGISTRATION

Clinicaltrials.gov NCT03840811 .

Host genetics highlights IFN-γ-dependent Toxoplasma genes encoding secreted and non-secreted virulence factors in in vivo CRISPR screens

Secreted virulence factors of Toxoplasma to survive in immune-competent hosts have been extensively explored by classical genetics and in vivo CRISPR screen methods, whereas their requirements in immune-deficient hosts are incompletely understood. Those of non-secreted virulence factors are further enigmatic. Here we develop an in vivo CRISPR screen system to enrich not only secreted but also non-secreted virulence factors in virulent Toxoplasma-infected C57BL/6 mice. Notably, combined usage of immune-deficient Ifngr1-/- mice highlights genes encoding various non-secreted proteins as well as well-known effectors such as ROP5, ROP18, GRA12, and GRA45 as interferon-γ (IFN-γ)-dependent virulence genes. The screen results suggest a role of GRA72 for normal GRA17/GRA23 localization and the IFN-γ-dependent role of UFMylation-related genes. Collectively, our study demonstrates that host genetics can complement in vivo CRISPR screens to highlight genes encoding IFN-γ-dependent secreted and non-secreted virulence factors in Toxoplasma.

Innate immune responses yield tissue-specific bottlenecks that scale with pathogen dose

To cause infection, pathogens must overcome bottlenecks imposed by the host immune system. These bottlenecks restrict the inoculum and largely determine whether pathogen exposure results in disease. Infection bottlenecks therefore quantify the effectiveness of immune barriers. Here, using a model of Escherichia coli systemic infection, we identify bottlenecks that tighten or widen with higher inoculum sizes, revealing that the efficacy of innate immune responses can increase or decrease with pathogen dose. We term this concept "dose scaling". During E. coli systemic infection, dose scaling is tissue specific, dependent on the LPS receptor TLR4, and can be recapitulated by mimicking high doses with killed bacteria. Scaling is therefore due to sensing of pathogen molecules rather than interactions between the host and live bacteria. We propose that dose scaling quantitatively links innate immunity with infection bottlenecks and is a valuable framework for understanding how the inoculum size governs the outcome of pathogen exposure.

Organ agar serves as physiologically relevant alternative for in vivo colonization

Animal models for host-microbial interactions have proven valuable, yielding physiologically relevant data that may be otherwise difficult to obtain. Unfortunately, such models are lacking or nonexistent for many microbes. Here, we introduce organ agar, a straightforward method to enable the screening of large mutant libraries while avoiding physiological bottlenecks. We demonstrate that growth defects on organ agar were translatable to colonization deficiencies in a murine model. Specifically, we present a urinary tract infection agar model to interrogate an ordered library of Proteus mirabilis transposon mutants, with accurate prediction of bacterial genes critical for host colonization. Thus, we demonstrate the ability of ex vivo organ agar to reproduce in vivo deficiencies. This work provides a readily adoptable technique that is economical and uses substantially fewer animals. We anticipate this method will be useful for a wide variety of microorganisms, both pathogenic and commensal, in a diverse range of model host species.

Novel method to recover Salmonella enterica cells for Tn-Seq approaches from lettuce leaves and agricultural environments using combination of sonication, filtration, and dialysis membrane

Salmonella enterica in agricultural environments has become an important concern, due to its potential transmission to humans and the associated public health risks. To identify genes contributing to Salmonella adaptation to such environments, transposon sequencing has been used in recent years. However, isolating Salmonella from atypical hosts, such as plant leaves, can pose technical challenges due to low bacterial content and the difficulty to separate an adequate number of bacteria from host tissues. In this study, we describe a modified methodology using a combination of sonication and filtration to recover S. enterica cells from lettuce leaves. We successfully recovered over a total of 3.5 × 106Salmonella cells in each biological replicate from two six-week old lettuce leaves, 7 days after infiltration with a Salmonella suspension of 5 × 107 colony forming units (CFU)/mL. Moreover, we have developed a dialysis membrane system as an alternative method for recovering bacteria from culture medium, mimicking a natural environment. Inoculating 107 CFU/mL of Salmonella into the media based on plant (lettuce and tomato) leaf and diluvial sand soil, a final concentration of 109.5 and 108.5 CFU/mL was obtained, respectively. One millilitre of the bacterial suspension after 24 h incubation at 28 °C using 60 rpm agitation was pelleted, corresponding to 109.5 and 108.5 cells from leaf- or soil-based media. The recovered bacterial population, from both lettuce leaves and environment-mimicking media, can adequately cover a presumptive library density of 106 mutants. In conclusion, this protocol provides an effective method to recover a Salmonella transposon sequencing library from in planta and in vitro systems. We expect this novel technique to foster the study of Salmonella in atypical hosts and environments, as well as other comparable scenarios.

Quantitative dose-response analysis untangles host bottlenecks to enteric infection

Host bottlenecks prevent many infections before the onset of disease by eliminating invading pathogens. By monitoring the diversity of a barcoded population of the diarrhea causing bacterium Citrobacter rodentium during colonization of its natural host, mice, we determine the number of cells that found the infection by establishing a replicative niche. In female mice the size of the pathogen's founding population scales with dose and is controlled by a severe yet slow-acting bottleneck. Reducing stomach acid or changing host genotype modestly relaxes the bottleneck without breaking the fractional relationship between dose and founders. In contrast, disrupting the microbiota causes the founding population to no longer scale with the size of the inoculum and allows the pathogen to infect at almost any dose, indicating that the microbiota creates the dominant bottleneck. Further, in the absence of competition with the microbiota, the diversity of the pathogen population slowly contracts as the population is overtaken by bacteria having lost the critical virulence island, the locus of enterocyte effacement (LEE). Collectively, our findings reveal that the mechanisms of protection by colonization bottlenecks are reflected in and can be generally defined by the impact of dose on the pathogen's founding population.

Gastric acid and escape to systemic circulation represent major bottlenecks to host infection by Citrobacter rodentium

The gastrointestinal (GI) environment plays a critical role in shaping enteric infections. Host environmental factors create bottlenecks, restrictive events that reduce the genetic diversity of invading bacterial populations. However, the identity and impact of bottleneck events on bacterial infection are largely unknown. We used Citrobacter rodentium infection of mice, a model of human pathogenic Escherichia coli infections, to examine bacterial population dynamics and quantify bottlenecks to host colonization. Using Sequence Tag-based Analysis of Microbial Populations (STAMP) we characterized the founding population size (Nb') and relatedness of C. rodentium populations at relevant tissue sites during early- and peak-infection. We demonstrate that the GI environment severely restricts the colonizing population, with an average Nb' of only 12-43 lineages (of 2,000+ inoculated) identified regardless of time or biogeographic location. Passage through gastric acid and escape to the systemic circulation were identified as major bottlenecks during C. rodentium colonization. Manipulating such events by increasing gastric pH dramatically increased intestinal Nb'. Importantly, removal of the stomach acid barrier had downstream consequences on host systemic colonization, morbidity, and mortality. These findings highlight the capability of the host GI environment to limit early pathogen colonization, controlling the population of initial founders with consequences for downstream infection outcomes.

Transposon mutagenesis in oral streptococcus

Oral streptococci are gram-positive facultative anaerobic bacteria that are normal inhabitants of the human oral cavity and play an important role in maintaining oral microecological balance and pathogenesis. Transposon mutagenesis is an effective genetic manipulation strategy for studying the function of genomic features. In order to study cariogenic related genes and crucial biological element genes of oral Streptococcus, transposon mutagenesis was widely used to identify functional genes. With the advent of next-generation sequencing (NGS) technology and the development of transposon random mutation library construction methods, transposon insertion sequencing (TIS) came into being. Benefiting from high-throughput advances in NGS, TIS was able to evaluate the fitness contribution and essentiality of genetic features in the bacterial genome. The application of transposon mutagenesis, including TIS, to oral streptococci provided a massive amount of valuable detailed linkage data between genetic fitness and genetic backgrounds, further clarify the processes of colonization, virulence, and persistence and provides a more reliable basis for investigating relationships with host ecology and disease status. This review focuses on transposon mutagenesis, including TIS, and its applicability in oral streptococci.

Construction of an Ordered Transposon Library for Uropathogenic Proteus mirabilis HI4320

Ordered transposon libraries are a valuable resource for many bacterial species, especially those with difficult methods for generating targeted genetic mutations. Here, we present the construction of an ordered transposon library for the bacterial urinary tract pathogen Proteus mirabilis strain HI4320. This library will facilitate future studies into P. mirabilis biology. For large experimental screens, it may be used to overcome bottleneck constraints and avoid biased outcomes resulting from gene length. For smaller studies, the library allows sidestepping the laborious construction of single targeted mutants. This library, containing 18,432 wells, was condensed into a smaller library containing 1,728 mutants. Each selected mutant had a single transposon insertion in an open reading frame, covering 45% of predicted genes encoded by P. mirabilis HI4320. This coverage was lower than expected and was due both to library wells with no mapped insertions and a surprisingly high proportion of mixed clones and multiple transposon insertion events. We offer recommendations for improving future library construction and suggestions for how to use this P. mirabilis library resource. IMPORTANCE Ordered libraries facilitate large genetic screens by guaranteeing high genomic coverage with a minimal number of mutants, and they can save time and effort by reducing the need to construct targeted mutations. This resource is now available for P. mirabilis, a common and complicating agent of catheter-associated urinary tract infection. We also present obstacles encountered during library construction with the goal to aid others who would like to construct ordered transposon libraries in other species.

Guided by the principles of microbiome engineering: Accomplishments and perspectives for environmental use

Although the accomplishments of microbiome engineering highlight its significance for the targeted manipulation of microbial communities, knowledge and technical gaps still limit the applications of microbiome engineering in biotechnology, especially for environmental use. Addressing the environmental challenges of refractory pollutants and fluctuating environmental conditions requires an adequate understanding of the theoretical achievements and practical applications of microbiome engineering. Here, we review recent cutting-edge studies on microbiome engineering strategies and their classical applications in bioremediation. Moreover, a framework is summarized for combining both top-down and bottom-up approaches in microbiome engineering toward improved applications. A strategy to engineer microbiomes for environmental use, which avoids the build-up of toxic intermediates that pose a risk to human health, is suggested. We anticipate that the highlighted framework and strategy will be beneficial for engineering microbiomes to address difficult environmental challenges such as degrading multiple refractory pollutants and sustain the performance of engineered microbiomes in situ with indigenous microorganisms under fluctuating conditions.

Characterizing viral within-host diversity in fast and non-equilibrium demo-genetic dynamics

High-throughput sequencing has opened the route for a deep assessment of within-host genetic diversity that can be used, e.g., to characterize microbial communities and to infer transmission links in infectious disease outbreaks. The performance of such characterizations and inferences cannot be analytically assessed in general and are often grounded on computer-intensive evaluations. Then, being able to simulate within-host genetic diversity across time under various demo-genetic assumptions is paramount to assess the performance of the approaches of interest. In this context, we built an original model that can be simulated to investigate the temporal evolution of genotypes and their frequencies under various demo-genetic assumptions. The model describes the growth and the mutation of genotypes at the nucleotide resolution conditional on an overall within-host viral kinetics, and can be tuned to generate fast non-equilibrium demo-genetic dynamics. We ran simulations of this model and computed classic diversity indices to characterize the temporal variation of within-host genetic diversity (from high-throughput amplicon sequences) of virus populations under three demographic kinetic models of viral infection. Our results highlight how demographic (viral load) and genetic (mutation, selection, or drift) factors drive variations in within-host diversity during the course of an infection. In particular, we observed a non-monotonic relationship between pathogen population size and genetic diversity, and a reduction of the impact of mutation on diversity when a non-specific host immune response is activated. The large variation in the diversity patterns generated in our simulations suggests that the underlying model provides a flexible basis to produce very diverse demo-genetic scenarios and test, for instance, methods for the inference of transmission links during outbreaks.

Mutation in the Two-Component System Regulator YycH Leads to Daptomycin Tolerance in Methicillin-Resistant Staphylococcus aureus upon Evolution with a Population Bottleneck

Adaptive laboratory evolution (ALE) is a useful tool to study the evolution of antibiotic tolerance in bacterial populations under diverse environmental conditions. The role of population bottlenecks in the evolution of tolerance has been investigated in Escherichia coli, but not in a more clinically relevant pathogen, methicillin-resistant Staphylococcus aureus (MRSA). In this study, we used ALE to evolve MRSA under repetitive daptomycin treatment and incorporated population bottlenecks following antibiotic exposure. We observed that the populations finally attained a tolerance mutation in the yycH gene after 2 weeks of evolution with population bottlenecks, and additional mutations in yycI and several other genes further increased the tolerance level. The tolerant populations also became resistant to another glycopeptide antibiotic, vancomycin. Through proteomics, we showed that yycH and yycI mutations led to the loss of function of the proteins and downregulated the WalKR two-component system and the downstream players, including the autolysin Atl and amidase Sle1, which are important for cell wall metabolism. Overall, our study offers new insights into the evolution of daptomycin tolerance under population bottlenecking conditions, which are commonly faced by pathogens during infection; the study also identified new mutations conferring daptomycin tolerance and revealed the proteome alterations in the evolved tolerant populations. IMPORTANCE Although population bottlenecks are known to influence the evolutionary dynamics of microbial populations, how such bottlenecks affect the evolution of tolerance to antibiotics in a clinically relevant methicillin-resistant S. aureus (MRSA) pathogen are still unclear. Here, we performed in vitro evolution of MRSA under cyclic daptomycin treatment and applied population bottlenecks following the treatment. We showed that under these experimental conditions, MRSA populations finally attained mutations in yycH, yycI, and several other genes that led to daptomycin tolerance. The discovered yycH and yycI mutations caused early termination of the genes and loss of function of the proteins, and they subsequently downregulated the expression of proteins controlled by the WalKR two-component system, such as Atl and Sle1. In addition, we compared our proteomics data with multiple studies on distinct daptomycin-tolerant MRSA mutants to identify proteins with a consistent expression pattern that could serve as biological markers for daptomycin tolerance in MRSA.

The ADP-Heptose Biosynthesis Enzyme GmhB is a Conserved Gram-Negative Bacteremia Fitness Factor

Klebsiella pneumoniae is a leading cause of Gram-negative bacteremia, which is a major source of morbidity and mortality worldwide. Gram-negative bacteremia requires three major steps: primary site infection, dissemination to the blood, and bloodstream survival. Because K. pneumoniae is a leading cause of health care-associated pneumonia, the lung is a common primary infection site leading to secondary bacteremia. K. pneumoniae factors essential for lung fitness have been characterized, but those required for subsequent bloodstream infection are unclear. To identify K. pneumoniae genes associated with dissemination and bloodstream survival, we combined previously and newly analyzed insertion site sequencing (InSeq) data from a murine model of bacteremic pneumonia. This analysis revealed the gene gmhB as important for either dissemination from the lung or bloodstream survival. In Escherichia coli, GmhB is a partially redundant enzyme in the synthesis of ADP-heptose for the lipopolysaccharide (LPS) core. To characterize its function in K. pneumoniae, an isogenic knockout strain (ΔgmhB) and complemented mutant were generated. During pneumonia, GmhB did not contribute to lung fitness and did not alter normal immune responses. However, GmhB enhanced bloodstream survival in a manner independent of serum susceptibility, specifically conveying resistance to spleen-mediated killing. In a tail-vein injection of murine bacteremia, GmhB was also required by K. pneumoniae, E. coli, and Citrobacter freundii for optimal fitness in the spleen and liver. Together, this study identifies GmhB as a conserved Gram-negative bacteremia fitness factor that acts through LPS-mediated mechanisms to enhance fitness in blood-filtering organs.

Genome-wide analysis of Brucella melitensis genes required throughout intranasal infection in mice

Brucellae are facultative intracellular Gram-negative coccobacilli that chronically infect various mammals and cause brucellosis. Human brucellosis is among the most common bacterial zoonoses and the vast majority of cases are attributed to B. melitensis. Using transposon sequencing (Tn-seq) analysis, we showed that among 3369 predicted genes of the B. melitensis genome, 861 are required for optimal growth in rich medium and 186 additional genes appeared necessary for survival of B. melitensis in RAW 264.7 macrophages in vitro. As the mucosal immune system represents the first defense against Brucella infection, we investigated the early phase of pulmonary infection in mice. In situ analysis at the single cell level indicates a succession of killing and growth phases, followed by heterogenous proliferation of B. melitensis in alveolar macrophages during the first 48 hours of infection. Tn-seq analysis identified 94 additional genes that are required for survival in the lung at 48 hours post infection. Among them, 42 genes are common to RAW 264.7 macrophages and the lung conditions, including the T4SS and purine synthesis genes. But 52 genes are not identified in RAW 264.7 macrophages, including genes implicated in lipopolysaccharide (LPS) biosynthesis, methionine transport, tryptophan synthesis as well as fatty acid and carbohydrate metabolism. Interestingly, genes implicated in LPS synthesis and β oxidation of fatty acids are no longer required in Interleukin (IL)-17RA^-/- mice and asthmatic mice, respectively. This demonstrates that the immune status determines which genes are required for optimal survival and growth of B. melitensis in vivo.

Brucellosis is one of the most widespread bacterial zoonoses worldwide. Using transposon sequencing (Tn-seq) analysis, we showed that among 3369 predicted genes of the Brucella melitensis genome, 861 are required for optimal growth in rich medium and 186 additional genes appeared necessary for survival of B. melitensis in RAW 264.7 macrophages in vitro. We also investigated the early phase of pulmonary infection in mice and identified 94 additional genes that are required for survival in the lung at 48 hours post infection. Among them, 42 genes are common to RAW 264.7 macrophages and the lung conditions, including the T4SS and purine synthesis genes. But 52 genes are not identified in RAW 264.7 macrophages, including genes implicated in lipopolysaccharide (LPS) biosynthesis, methionine transport, tryptophan synthesis as well as fatty acid and carbohydrate metabolism. Interestingly, genes implicated in LPS synthesis and β oxidation of fatty acids are no longer required in Interleukin (IL)-17RA^-/- mice and asthmatic mice, respectively. Our work demonstrates that both the immune status and the nature of the infected cell type determines which genes are required for optimal survival and growth of B. melitensis in vivo.

Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections

During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.

The bacterium Staphylococcus aureus lives harmlessly on our skin and noses. However, occasionally, it gets into our blood and internal organs, such as our bones and joints, where it causes severe, long-lasting infections that are difficult to treat.

Over time, S. aureus acquire characteristics that help them to adapt to different locations, such as transitioning from the nose to the blood, and avoid being killed by antibiotics. Previous studies have identified changes, or ‘mutations’, in genes that are likely to play an important role in this evolutionary process. One of these genes, called accessory gene regulator (or agr for short), has been shown to control the mechanisms S. aureus use to infect cells and disseminate in the body. However, it is unclear if there are changes in other genes that also help S. aureus adapt to life inside the human body.

To help resolve this mystery, Giulieri et al. collected 2,500 samples of S. aureus from almost 400 people. This included bacteria harmlessly living on the skin or in the nose, as well as strains that caused an infection. Gene sequencing revealed a small number of genes, referred to as ‘adaptive genes’, that often acquire mutations during infection. Of these, agr was the most commonly altered. However, mutations in less well-known genes were also identified: some of these genes are related to resistance to antibiotics, while others are involved in chemical processes that help the bacteria to process nutrients.

Most mutations were caused by random errors being introduced in to the bacteria’s genetic code which stopped genes from working. However, in some cases, genes were turned off by small fragments of DNA moving around and inserting themselves into different parts of the genome.

This study highlights a group of genes that help S. aureus to thrive inside the body and cause severe and prolonged infections. If these results can be confirmed, it may help to guide which antibiotics are used to treat different infections. Furthermore, understanding which genes are important for infection could lead to new strategies for eliminating this dangerous bacterium.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

Genetically distant bacteriophages select for unique genomic changes in Enterococcus faecalis

The human microbiota harbors diverse bacterial and bacteriophage (phage) communities. Bacteria evolve to overcome phage infection, thereby driving phage evolution to counter bacterial resistance. Understanding how phages select for genetic alterations in medically relevant bacteria is important as phages become established biologics for the treatment of multidrug‐resistant (MDR) bacterial infections. Before phages can be widely used as standalone or combination antibacterial therapies, we must obtain a deep understanding of the molecular mechanisms of phage infection and how host bacteria alter their genomes to become resistant. We performed coevolution experiments using a single Enterococcus faecalis strain and two distantly related phages to determine how phage pressure impacts the evolution of the E. faecalis genome. Whole‐genome sequencing of E. faecalis following continuous exposure to these two phages revealed mutations previously demonstrated to be essential for phage infection. We also identified mutations in genes previously unreported to be associated with phage infection in E. faecalis. Intriguingly, there was only one shared mutation in the E. faecalis genome that was selected by both phages tested, demonstrating that infection by two genetically distinct phages selects for diverse variants. This knowledge serves as the basis for the continued study of E. faecalis genome evolution during phage infection and can be used to inform the design of future therapeutics, such as phage cocktails, intended to target MDR E. faecalis.

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

Bacterial Lymphatic Metastasis in Infection and Immunity

Lymphatic vessels permeate tissues around the body, returning fluid from interstitial spaces back to the blood after passage through the lymph nodes, which are important sites for adaptive responses to all types of pathogens. Involvement of the lymphatics in the pathogenesis of bacterial infections is not well studied. Despite offering an obvious conduit for pathogen spread, the lymphatic system has long been regarded to bar the onward progression of most bacteria. There is little direct data on live virulent bacteria, instead understanding is largely inferred from studies investigating immune responses to viruses or antigens in lymph nodes. Recently, we have demonstrated that extracellular bacterial lymphatic metastasis of virulent strains of Streptococcus pyogenes drives systemic infection. Accordingly, it is timely to reconsider the role of lymph nodes as absolute barriers to bacterial dissemination in the lymphatics. Here, we summarise the routes and mechanisms by which an increasing variety of bacteria are acknowledged to transit through the lymphatic system, including those that do not necessarily require internalisation by host cells. We discuss the anatomy of the lymphatics and other factors that influence bacterial dissemination, as well as the consequences of underappreciated bacterial lymphatic metastasis on disease and immunity.

The Infectious Dose Shapes Vibrio cholerae Within-Host Dynamics

During infection, the rates of pathogen replication, death, and migration affect disease progression, dissemination, transmission, and resistance evolution. Here, we follow the population dynamics of Vibrio cholerae in a mouse model by labeling individual bacteria with one of >500 unique, fitness-neutral genomic tags. Using the changes in tag frequencies and CFU numbers, we inform a mathematical model that describes the within-host spatiotemporal bacterial dynamics. This allows us to disentangle growth, death, forward, and retrograde migration rates continuously during infection. Our model has robust predictive power across various experimental setups. The population dynamics of V. cholerae shows substantial spatiotemporal heterogeneity in replication, death, and migration. Importantly, we find that the niche available to V. cholerae in the host increases with inoculum size, suggesting cooperative effects during infection. Therefore, it is not enough to consider just the likelihood of exposure (50% infectious dose) but rather the magnitude of exposure to predict outbreaks. IMPORTANCE Determining the rates of bacterial migration, replication, and death during infection is important for understanding how infections progress. Separately measuring these rates is often difficult in systems where multiple processes happen simultaneously. Here, we use next-generation sequencing to measure V. cholerae migration, replication, death, and niche size along the mouse gastrointestinal tract. We show that the small intestine of the mouse is a heterogeneous environment, and the population dynamic characteristics change substantially between adjacent gut sections. Our approach also allows us to characterize the effect of inoculum size on these processes. We find that the niche size in mice increases with the infectious dose, hinting at cooperative effects in larger inocula. The dose-response relationship between inoculum size and final pathogen burden is important for the infected individual and is thought to influence the progression of V. cholerae epidemics.

Pathogen clonal expansion underlies multiorgan dissemination and organ-specific outcomes during murine systemic infection

The dissemination of pathogens through blood and their establishment within organs lead to severe clinical outcomes. However, the within-host dynamics that underlie pathogen spread to and clearance from systemic organs remain largely uncharacterized. In animal models of infection, the observed pathogen population results from the combined contributions of bacterial replication, persistence, death, and dissemination, each of which can vary across organs. Quantifying the contribution of each these processes is required to interpret and understand experimental phenotypes. Here, we leveraged STAMPR, a new barcoding framework, to investigate the population dynamics of extraintestinal pathogenic Escherichia coli, a common cause of bacteremia, during murine systemic infection. We show that while bacteria are largely cleared by most organs, organ-specific clearance failures are pervasive and result from dramatic expansions of clones representing less than 0.0001% of the inoculum. Clonal expansion underlies the variability in bacterial burden between animals, and stochastic dissemination of clones profoundly alters the pathogen population structure within organs. Despite variable pathogen expansion events, host bottlenecks are consistent yet highly sensitive to infection variables, including inoculum size and macrophage depletion. We adapted our barcoding methodology to facilitate multiplexed validation of bacterial fitness determinants identified with transposon mutagenesis and confirmed the importance of bacterial hexose metabolism and cell envelope homeostasis pathways for organ-specific pathogen survival. Collectively, our findings provide a comprehensive map of the population biology that underlies bacterial systemic infection and a framework for barcode-based high-resolution mapping of infection dynamics.

Refined Quantification of Infection Bottlenecks and Pathogen Dissemination with STAMPR

Pathogen population dynamics during infection are critical determinants of infection susceptibility and define patterns of dissemination. However, deciphering these dynamics, particularly founding population sizes in host organs and patterns of dissemination between organs, is difficult because measuring bacterial burden alone is insufficient to observe these patterns. Introduction of allelic diversity into otherwise identical bacteria using DNA barcodes enables sequencing-based measurements of these parameters, in a method known as STAMP (Sequence Tag-based Analysis of Microbial Populations). However, bacteria often undergo unequal expansion within host organs, resulting in marked differences in the frequencies of barcodes in input and output libraries. Here, we show that these differences confound STAMP-based analyses of founding population sizes and dissemination patterns. We present STAMPR, a successor to STAMP, which accounts for such population expansions. Using data from systemic infection of barcoded extraintestinal pathogenic E. coli, we show that this new framework, along with the metrics it yields, enhances the fidelity of measurements of bottlenecks and dissemination patterns. STAMPR was also validated on an independent barcoded Pseudomonas aeruginosa data set, uncovering new patterns of dissemination within the data. This framework (available at https://github.com/hullahalli/stampr_rtisan), when coupled with barcoded data sets, enables a more complete assessment of within-host bacterial population dynamics. IMPORTANCE Barcoded bacteria are often employed to monitor pathogen population dynamics during infection. The accuracy of these measurements is diminished by unequal bacterial expansion rates. Here, we develop computational tools to circumvent this limitation and establish additional metrics that collectively enhance the fidelity of measuring within-host pathogen founding population sizes and dissemination patterns. These new tools will benefit future studies of the dynamics of pathogens and symbionts within their respective hosts and may have additional barcode-based applications beyond host-microbe interactions.

Population Bottlenecks Strongly Affect the Evolutionary Dynamics of Antibiotic Persistence

Bacterial persistence is a potential cause of antibiotic therapy failure. Antibiotic-tolerant persisters originate from phenotypic differentiation within a susceptible population, occurring with a frequency that can be altered by mutations. Recent studies have proven that persistence is a highly evolvable trait and, consequently, an important evolutionary strategy of bacterial populations to adapt to high-dose antibiotic therapy. Yet, the factors that govern the evolutionary dynamics of persistence are currently poorly understood. Theoretical studies predict far-reaching effects of bottlenecking on the evolutionary adaption of bacterial populations, but these effects have never been investigated in the context of persistence. Bottlenecking events are frequently encountered by infecting pathogens during host-to-host transmission and antibiotic treatment. In this study, we used a combination of experimental evolution and barcoded knockout libraries to examine how population bottlenecking affects the evolutionary dynamics of persistence. In accordance with existing hypotheses, small bottlenecks were found to restrict the adaptive potential of populations and result in more heterogeneous evolutionary outcomes. Evolutionary trajectories followed in small-bottlenecking regimes additionally suggest that the fitness landscape associated with persistence has a rugged topography, with distinct trajectories toward increased persistence that are accessible to evolving populations. Furthermore, sequencing data of evolved populations and knockout libraries after selection reveal various genes that are potentially involved in persistence, including previously known as well as novel targets. Together, our results do not only provide experimental evidence for evolutionary theories, but also contribute to a better understanding of the environmental and genetic factors that guide bacterial adaptation to antibiotic treatment.

Quantification of Brucella abortus population structure in a natural host

Cattle are natural hosts of the intracellular pathogen Brucella abortus, which inflicts a significant burden on the health and reproduction of these important livestock. The primary routes of infection in field settings have been described, but it is not known how the bovine host shapes the structure of B. abortus populations during infection. We utilized a library of uniquely barcoded B. abortus strains to temporally and spatially quantify population structure during colonization of cattle through a natural route of infection. Introducing 10⁸ bacteria from this barcoded library to the conjunctival mucosa resulted in expected levels of local lymph node colonization at a 1-wk time point. We leveraged variance in strain abundance in the library to demonstrate that only 1 in 10,000 brucellae introduced at the site of infection reached a parotid lymph node. Thus, cattle restrict the overwhelming majority of B. abortus introduced via the ocular conjunctiva at this dose. Individual strains were spatially restricted within the host tissue, and the total B. abortus census was dominated by a small number of distinct strains in each lymph node. These results define a bottleneck that B. abortus must traverse to colonize local lymph nodes from the conjunctival mucosa. The data further support a model in which a small number of spatially isolated granulomas founded by unique strains are present at 1 wk postinfection. These experiments demonstrate the power of barcoded transposon tools to quantify infection bottlenecks and to define pathogen population structure in host tissues.

barCoder: a tool to generate unique, orthogonal genetic tags for qPCR detection

BACKGROUND

Tracking dispersal of microbial populations in the environment requires specific detection methods that discriminate between the target strain and all potential natural and artificial interferents, including previously utilized tester strains. Recent work has shown that genomic insertion of short identification tags, called "barcodes" here, allows detection of chromosomally tagged strains by real-time PCR. Manual design of these barcodes is feasible for small sets, but expansion of the technique to larger pools of distinct and well-functioning assays would be significantly aided by software-guided design.

RESULTS

Here we introduce barCoder, a bioinformatics tool that facilitates the process of creating sets of uniquely identifiable barcoded strains. barCoder utilizes the genomic sequence of the target strain and a set of user-specified PCR parameters to generate a list of suggested barcode "modules" that consist of binding sites for primers and probes, and appropriate spacer sequences. Each module is designed to yield optimal PCR amplification and unique identification. Optimal amplification includes metrics such as ideal melting temperature and G+C content, appropriate spacing, and minimal stem-loop formation; unique identification includes low BLAST hits against the target organism, previously generated barcode modules, and databases (such as NCBI). We tested the ability of our algorithm to suggest appropriate barcodes by generating 12 modules for Bacillus thuringiensis serovar kurstaki-a simulant for the potential biowarfare agent Bacillus anthracis-and three each for other potential target organisms with variable G+C content. Real-time PCR detection assays directed at barcodes were specific and yielded minimal cross-reactivity with a panel of near-neighbor and potential contaminant materials.

CONCLUSIONS

The barCoder algorithm facilitates the generation of synthetically barcoded biological simulants by (a) eliminating the task of creating modules by hand, (b) minimizing optimization of PCR assays, and (c) reducing effort wasted on non-unique barcode modules.

Clonal Interference and Mutation Bias in Small Bacterial Populations in Droplets

Experimental evolution studies have provided key insights into the fundamental mechanisms of evolution. One striking observation is that parallel and convergent evolution during laboratory evolution can be surprisingly common. However, these experiments are typically performed with well-mixed cultures and large effective population sizes, while pathogenic microbes typically experience strong bottlenecks during infection or drug treatment. Yet, our knowledge about adaptation in very small populations, where selection strength and mutation supplies are limited, is scant. In this study, wild-type and mutator strains of the bacterium Escherichia coli were evolved for about 100 generations towards increased resistance to the β-lactam antibiotic cefotaxime in millifluidic droplets of 0.5 µL and effective population size of approximately 27,000 cells. The small effective population size limited the adaptive potential of wild-type populations, where adaptation was limited to inactivating mutations, which caused the increased production of outer-membrane vesicles, leading to modest fitness increases. In contrast, mutator clones with an average of ~30-fold higher mutation rate adapted much faster by acquiring both inactivating mutations of an outer-membrane porin and particularly inactivating and gain-of-function mutations, causing the upregulation or activation of a common efflux pump, respectively. Our results demonstrate how in very small populations, clonal interference and mutation bias together affect the choice of adaptive trajectories by mediating the balance between high-rate and large-benefit mutations.

RESTAMP - Rate estimates by sequence-tag analysis of microbial populations

Microbial division rates determine the speed of mutation accumulation and thus the emergence of antimicrobial resistance. Microbial death rates are affected by antibiotic action and the immune system. Therefore, measuring these rates has advanced our understanding of host-pathogen interactions and antibiotic action. Several methods based on marker-loss or few inheritable neutral markers exist that allow estimating microbial division and death rates, each of which has advantages and limitations. Technical bottlenecks, i.e., experimental sampling events, during the experiment can distort the rate estimates and are typically unaccounted for or require additional calibration experiments. In this work, we introduce RESTAMP (Rate Estimates by Sequence Tag Analysis of Microbial Populations) as a method for determining bacterial division and death rates. This method uses hundreds of fitness neutral sequence barcodes to measure the rates and account for experimental bottlenecks at the same time. We experimentally validate RESTAMP and compare it to established plasmid loss methods. We find that RESTAMP has a number of advantages over plasmid loss or previous marker based techniques. (i) It enables to correct the distortion of rate estimates by technical bottlenecks. (ii) Rate estimates are independent of the sequence tag distribution in the starting culture allowing the use of an arbitrary number of tags. (iii) It introduces a bottleneck sensitivity measure that can be used to maximize the accuracy of the experiment. RESTAMP allows studying microbial population dynamics with great resolution over a wide dynamic range and can thus advance our understanding of host-pathogen interactions or the mechanisms of antibiotic action.

Bottleneck size and selection level reproducibly impact evolution of antibiotic resistance

During antibiotic treatment, the evolution of bacterial pathogens is fundamentally affected by bottlenecks and varying selection levels imposed by the drugs. Bottlenecks-that is, reductions in bacterial population size-lead to an increased influence of random effects (genetic drift) during bacterial evolution, and varying antibiotic concentrations during treatment may favour distinct resistance variants. Both aspects influence the process of bacterial evolution during antibiotic therapy and thereby treatment outcome. Surprisingly, the joint influence of these interconnected factors on the evolution of antibiotic resistance remains largely unexplored. Here we combine evolution experiments with genomic and genetic analyses to demonstrate that bottleneck size and antibiotic-induced selection reproducibly impact the evolutionary path to resistance in pathogenic Pseudomonas aeruginosa, one of the most problematic opportunistic human pathogens. Resistance is favoured-expectedly-under high antibiotic selection and weak bottlenecks, but-unexpectedly-also under low antibiotic selection and severe bottlenecks. The latter is likely to result from a reduced probability of losing favourable variants through drift under weak selection. Moreover, the absence of high resistance under low selection and weak bottlenecks is caused by the spread of low-resistance variants with high competitive fitness under these conditions. We conclude that bottlenecks, in combination with drug-induced selection, are currently neglected key determinants of pathogen evolution and outcome of antibiotic treatment.

An experimental design tool to optimize inference precision in data-driven mathematical models of bacterial infections in vivo

The management of bacterial diseases calls for a detailed knowledge about the dynamic changes in host-bacteria interactions. Biological insights are gained by integrating experimental data with mechanistic mathematical models to infer experimentally unobservable quantities. This inter-disciplinary field would benefit from experiments with maximal information content yielding high-precision inference. Here, we present a computationally efficient tool for optimizing experimental design in terms of parameter inference in studies using isogenic-tagged strains. We study the effect of three experimental design factors: number of biological replicates, sampling timepoint selection and number of copies per tagged strain. We conduct a simulation study to establish the relationship between our optimality criterion and the size of parameter estimate confidence intervals, and showcase its application in a range of biological scenarios reflecting different dynamics patterns observed in experimental infections. We show that in low-variance systems with low killing and replication rates, predicting high-precision experimental designs is consistently achieved; higher replicate sizes and strategic timepoint selection yield more precise estimates. Finally, we address the question of resource allocation under constraints; given a fixed number of host animals and a constraint on total inoculum size per host, infections with fewer strains at higher copies per strain lead to higher-precision inference.

Exploration of Bacterial Bottlenecks and Streptococcus pneumoniae Pathogenesis by CRISPRi-Seq

Streptococcus pneumoniae is an opportunistic human pathogen that causes invasive diseases, including pneumonia, with greater health risks upon influenza A virus (IAV) co-infection. To facilitate pathogenesis studies in vivo, we developed an inducible CRISPR interference system that enables genome-wide fitness testing in one sequencing step (CRISPRi-seq). We applied CRISPRi-seq to assess bottlenecks and identify pneumococcal genes important in a murine pneumonia model. A critical bottleneck occurs at 48 h with few bacteria causing systemic infection. This bottleneck is not present during IAV superinfection, facilitating identification of pneumococcal pathogenesis-related genes. Top in vivo essential genes included purA, encoding adenylsuccinate synthetase, and the cps operon required for capsule production. Surprisingly, CRISPRi-seq indicated no fitness-related role for pneumolysin during superinfection. Interestingly, although metK (encoding S-adenosylmethionine synthetase) was essential in vitro, it was dispensable in vivo. This highlights advantages of CRISPRi-seq over transposon-based genetic screens, as all genes, including essential genes, can be tested for pathogenesis potential.

Elucidating host-microbe interactions in vivo by studying population dynamics using neutral genetic tags

Host–microbe interactions are highly dynamic in space and time, in particular in the case of infections. Pathogen population sizes, microbial phenotypes and the nature of the host responses often change dramatically over time. These features pose particular challenges when deciphering the underlying mechanisms of these interactions experimentally, as traditional microbiological and immunological methods mostly provide snapshots of population sizes or sparse time series. Recent approaches – combining experiments using neutral genetic tags with stochastic population dynamic models – allow more precise quantification of biologically relevant parameters that govern the interaction between microbe and host cell populations. This is accomplished by exploiting the patterns of change of tag composition in the microbe or host cell population under study. These models can be used to predict the effects of immunodeficiencies or therapies (e.g. antibiotic treatment) on populations and thereby generate hypotheses and refine experimental designs. In this review, we present tools to study population dynamics in vivo using genetic tags, explain examples for their implementation and briefly discuss future applications.

Adapt or Perish: Evolutionary Rescue in a Gradually Deteriorating Environment

We investigate the evolutionary rescue of a microbial population in a gradually deteriorating environment, through a combination of analytical calculations and stochastic simulations. We consider a population destined for extinction in the absence of mutants, which can survive only if mutants sufficiently adapted to the new environment arise and fix. We show that mutants that appear later during the environment deterioration have a higher probability to fix. The rescue probability of the population increases with a sigmoidal shape when the product of the carrying capacity and of the mutation probability increases. Furthermore, we find that rescue becomes more likely for smaller population sizes and/or mutation probabilities if the environment degradation is slower, which illustrates the key impact of the rapidity of environment degradation on the fate of a population. We also show that our main conclusions are robust across various types of adaptive mutants, including specialist and generalist ones, as well as mutants modeling antimicrobial resistance evolution. We further express the average time of appearance of the mutants that do rescue the population and the average extinction time of those that do not. Our methods can be applied to other situations with continuously variable fitnesses and population sizes, and our analytical predictions are valid in the weak-to-moderate mutation regime.

Adaptation of Arginine Synthesis among Uropathogenic Branches of the Escherichia coli Phylogeny Reveals Adjustment to the Urinary Tract Habitat

Uropathogenic Escherichia coli (UPEC) is the most common cause of human urinary tract infection (UTI). Population bottlenecks during early stages of UTI make high-throughput screens impractical for understanding clinically important later stages of UTI, such as persistence and recurrence. As UPEC is hypothesized to be adapted to these later pathogenic stages, we previously identified 29 genes evolving under positive selection in UPEC. Here, we found that 8 of these genes, including argI (which is involved in arginine biosynthesis), are important for persistence in a mouse model of UTI. Deletion of argI and other arginine synthesis genes resulted in (i) arginine auxotrophy and (ii) defects in persistent UTI. Replacement of a B2 clade argI with a non-B2 clade argI complemented arginine auxotrophy, but the resulting strain remained attenuated in its ability to cause persistent bacteriuria. Thus, argI may have a second function during UTI that is not related to simple arginine synthesis. This study demonstrates how variation in metabolic genes can impact virulence and provides insight into the mechanisms and evolution of bacterial virulence.

Urinary tract infections (UTIs) are predominantly caused by uropathogenic Escherichia coli (UPEC). UPEC pathogenesis requires passage through a severe population bottleneck involving intracellular bacterial communities (IBCs) that are clonal expansions of a single invading UPEC bacterium in a urothelial superficial facet cell. IBCs occur only during acute pathogenesis. The bacteria in IBCs form the founder population that develops into persistent extracellular infections. Only a small fraction of UPEC organisms proceed through the IBC cycle, regardless of the inoculum size. This dramatic reduction in population size precludes the utility of genomic mutagenesis technologies for identifying genes important for persistence. To circumvent this bottleneck, we previously identified 29 positively selected genes (PSGs) within UPEC and hypothesized that they contribute to virulence. Here, we show that 8 of these 29 PSGs are required for fitness during persistent bacteriuria. Conversely, 7/8 of these PSG mutants showed essentially no phenotype in acute UTI. Deletion of the PSG argI leads to arginine auxotrophy. Relative to the other arg genes, argI in the B2 clade (which comprises most UPEC strains) of E. coli has diverged from argI in other E. coli clades. Replacement of argI in a UPEC strain with a non-UPEC argI allele complemented the arginine auxotrophy but not the persistent bacteriuria defect, showing that the UPEC argI allele contributes to persistent infection. These results highlight the complex roles of metabolic pathways during infection and demonstrate that evolutionary approaches can identify infection-specific gene functions downstream of population bottlenecks, shedding light on virulence and the genetic evolution of pathogenesis.

Nutrient levels and trade-offs control diversity in a serial dilution ecosystem

Microbial communities feature an immense diversity of species and this diversity is linked to outcomes ranging from ecosystem stability to medical prognoses. Yet the mechanisms underlying microbial diversity are under debate. While simple resource-competition models don't allow for coexistence of a large number of species, it was recently shown that metabolic trade-offs can allow unlimited diversity. Does this diversity persist with more realistic, intermittent nutrient supply? Here, we demonstrate theoretically that in serial dilution culture, metabolic trade-offs allow for high diversity. When a small amount of nutrient is supplied to each batch, the serial dilution dynamics mimic a chemostat-like steady state. If more nutrient is supplied, community diversity shifts due to an 'early-bird' effect. The interplay of this effect with different environmental factors and diversity-supporting mechanisms leads to a variety of relationships between nutrient supply and diversity, suggesting that real ecosystems may not obey a universal nutrient-diversity relationship.

Within-host spatiotemporal dynamic of systemic salmonellosis: Ways to track infection, reaction to vaccination and antimicrobial treatment

During the last two decades our understanding of the complex in vivo host-pathogen interactions has increased due to technical improvements and new research tools. The rapid advancement of molecular biology, flow cytometry and microscopy techniques, combined with mathematical modelling, have empowered in-depth studies of systemic bacterial infections across scales from single molecules, to cells, to organs and systems to reach the whole organism level. By tracking subpopulations of bacteria in vivo using molecular or fluorescent tags, it has been possible to reconstruct the spread of infection within and between organs, allowing unprecedented quantification of the effects of antimicrobial treatment and vaccination. This review illustrates recent advances in the study of heterogeneous traits of the infection process and illustrate approaches to investigate the reciprocal interactions between antimicrobial treatments, bacterial growth/death as well as inter- and intra-organ spread. We also discuss how vaccines impact the in vivo behaviour of bacteria and how these findings can guide vaccine design and rational antimicrobial treatment.

A decade of advances in transposon-insertion sequencing

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

In this Review, several experts discuss progress in the decade since the development of transposon-based approaches for bacterial genetic screens. They describe how advances in both experimental technologies and analytical strategies are resulting in insights into diverse biological processes.