Bacterial Population Genetics Research Articles

Ectoparasite and bacterial population genetics and community structure indicate extent of bat movement across an island chain.

Few studies have examined the genetic population structure of vector-borne microparasites in wildlife, making it unclear how much these systems can reveal about the movement of their associated hosts. This study examined the complex host–vector–microbe interactions in a system of bats, wingless ectoparasitic bat flies (Nycteribiidae), vector-borne microparasitic bacteria (Bartonella) and bacterial endosymbionts of flies (Enterobacterales) across an island chain in the Gulf of Guinea, West Africa. Limited population structure was found in bat flies and Enterobacterales symbionts compared to that of their hosts. Significant isolation by distance was observed in the dissimilarity of Bartonella communities detected in flies from sampled populations of Eidolon helvum bats. These patterns indicate that, while genetic dispersal of bats between islands is limited, some non-reproductive movements may lead to the dispersal of ectoparasites and associated microbes. This study deepens our knowledge of the phylogeography of African fruit bats, their ectoparasites and associated bacteria. The results presented could inform models of pathogen transmission in these bat populations and increase our theoretical understanding of community ecology in host–microbe systems.

An Open-Source Program (Haplo-ST) for Whole-Genome Sequence Typing Shows Extensive Diversity among Listeria monocytogenes Isolates in Outdoor Environments and Poultry Processing Plants.

A reliable and standardized classification of Listeria monocytogenes is important for accurate strain identification during outbreak investigations. Current whole-genome sequencing (WGS)-based approaches for strain characterization are either difficult to standardize, rendering them less suitable for data exchange, or are not freely available. Thus, we developed a portable and open-source tool, Haplo-ST, to improve standardization and provide maximum discriminatory potential to WGS data tied to a multilocus sequence typing (MLST) framework. Haplo-ST performs whole-genome MLST (wgMLST) for L. monocytogenes while allowing for data exchangeability worldwide. This tool takes in (i) raw WGS reads as input, (ii) cleans the raw data according to user-specified parameters, (iii) assembles genes across loci by mapping to genes from reference strains, and (iv) assigns allelic profiles to assembled genes and provides a wgMLST subtyping for each isolate. Data exchangeability relies on the tool assigning allelic profiles based on a centralized nomenclature defined by the widely used BIGSdb-Lm database. Tests of Haplo-ST's performance with simulated reads from L. monocytogenes reference strains demonstrated high sensitivity (97.5%), and coverage depths of ≥20× were found to be sufficient for wgMLST profiling. We then used Haplo-ST to characterize and differentiate between two groups of L. monocytogenes isolates derived from the natural environment and poultry processing plants. Phylogenetic reconstruction identified lineages within each group, and no lineage specificity was observed with isolate phenotypes (transient versus persistent) or origins. Genetic differentiation analyses between isolate groups identified 21 significantly differentiated loci, potentially enriched for adaptation and persistence of L. monocytogenes within poultry processing plants.IMPORTANCE We have developed an open-source tool (https://github.com/swarnalilouha/Haplo-ST) that provides allele-based subtyping of L. monocytogenes isolates at the whole-genome level. Along with allelic profiles, this tool also generates allele sequences and identifies paralogs, which is useful for phylogenetic tree reconstruction and deciphering relationships between closely related isolates. More broadly, Haplo-ST is flexible and can be adapted to characterize the genome of any haploid organism simply by installing an organism-specific gene database. Haplo-ST also allows for scalable subtyping of isolates; fewer reference genes can be used for low-resolution typing, whereas higher resolution can be achieved by increasing the number of genes used in the analysis. Our tool enabled clustering of L. monocytogenes isolates into lineages and detection of potential loci for adaptation and persistence in food processing environments. Findings from these analyses highlight the effectiveness of Haplo-ST in subtyping and evaluating relationships among isolates in studies of bacterial population genetics.

Consistent Metagenome-Derived Metrics Verify and Delineate Bacterial Species Boundaries.

Longstanding questions relate to the existence of naturally distinct bacterial species and genetic approaches to distinguish them. Bacterial genomes in public databases form distinct groups, but these databases are subject to isolation and deposition biases. To avoid these biases, we compared 5,203 bacterial genomes from 1,457 environmental metagenomic samples to test for distinct clouds of diversity and evaluated metrics that could be used to define the species boundary. Bacterial genomes from the human gut, soil, and the ocean all exhibited gaps in whole-genome average nucleotide identities (ANI) near the previously suggested species threshold of 95% ANI. While genome-wide ratios of nonsynonymous and synonymous nucleotide differences (dN/dS) decrease until ANI values approach ∼98%, two methods for estimating homologous recombination approached zero at ∼95% ANI, supporting breakdown of recombination due to sequence divergence as a species-forming force. We evaluated 107 genome-based metrics for their ability to distinguish species when full genomes are not recovered. Full-length 16S rRNA genes were least useful, in part because they were underrecovered from metagenomes. However, many ribosomal proteins displayed both high metagenomic recoverability and species discrimination power. Taken together, our results verify the existence of sequence-discrete microbial species in metagenome-derived genomes and highlight the usefulness of ribosomal genes for gene-level species discrimination.IMPORTANCE There is controversy about whether bacterial diversity is clustered into distinct species groups or exists as a continuum. To address this issue, we analyzed bacterial genome databases and reports from several previous large-scale environment studies and identified clear discrete groups of species-level bacterial diversity in all cases. Genetic analysis further revealed that quasi-sexual reproduction via horizontal gene transfer is likely a key evolutionary force that maintains bacterial species integrity. We next benchmarked over 100 metrics to distinguish these bacterial species from each other and identified several genes encoding ribosomal proteins with high species discrimination power. Overall, the results from this study provide best practices for bacterial species delineation based on genome content and insight into the nature of bacterial species population genetics.

Defining and combating antibiotic resistance from One Health and Global Health perspectives

Several interconnected human, animal and environmental habitats can contribute to the emergence, evolution and spread of antibiotic resistance, and the health of these contiguous habitats (the focus of the One Health approach) may represent a risk to human health. Additionally, the expansion of resistant clones and antibiotic resistance determinants among human-associated, animal-associated and environmental microbiomes have the potential to alter bacterial population genetics at local and global levels, thereby modifying the structure, and eventually the productivity, of microbiomes where antibiotic-resistant bacteria can expand. Conversely, any change in these habitats (including pollution by antibiotics or by antibiotic-resistant organisms) may influence the structures of their associated bacterial populations, which might affect the spread of antibiotic resistance to, and among, the above-mentioned microbiomes. Besides local transmission among connected habitats-the focus of studies under the One Health concept-the transmission of resistant microorganisms might occur on a broader (even worldwide) scale, requiring coordinated Global Health actions. This Review provides updated information on the elements involved in the evolution and spread of antibiotic resistance at local and global levels, and proposes studies to be performed and strategies to be followed that may help reduce the burden of antibiotic resistance as well as its impact on human and planetary health.

MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies.

Large-scale bacterial population genetics studies are now routine due to cost-effective Illumina short-read sequencing. However, analysing plasmid content remains difficult due to incomplete assembly of plasmids. Bacterial isolates can contain any number of plasmids and assembly remains complicated due to the presence of repetitive elements. Numerous tools have been developed to analyse plasmids but the performance and functionality of the tools are variable. The MOB-suite was developed as a set of modular tools for reconstruction and typing of plasmids from draft assembly data to facilitate characterization of plasmids. Using a set of closed genomes with publicly available Illumina data, the MOB-suite identified contigs of plasmid origin with both high sensitivity and specificity (95 and 88 %, respectively). In comparison, plasmidfinder demonstrated high specificity (99 %) but limited sensitivity (50 %). Using the same dataset of 377 known plasmids, MOB-recon accurately reconstructed 207 plasmids so that they were assigned to a single grouping without other plasmid or chromosomal sequences, whereas plasmidSPAdes was only able to accurately reconstruct 102 plasmids. In general, plasmidSPAdes has a tendency to merge different plasmids together, with 208 plasmids undergoing merge events. The MOB-suite reduces the number of errors but produces more hybrid plasmids, with 84 plasmids undergoing both splits and merges. The MOB-suite also provides replicon typing similar to plasmidfinder but with the inclusion of relaxase typing and prediction of conjugation potential. The MOB-suite is written in Python 3 and is available from https://github.com/phac-nml/mob-suite.

Large-Scale Simulations of Bacterial Populations Over Complex Networks.

The understanding of bacterial population genetics and evolution is crucial in epidemic outbreak studies and pathogen surveillance. However, all epidemiological studies are limited to their sampling capacities which, by being usually biased or limited due to economic constraints, can hamper the real knowledge of the bacterial population structure of a given species. To this end, mathematical models and large-scale simulations can provide a quantitative analytical framework that can be used to assess how or if limited sampling can infer the true population structure. In this article, we address the large-scale simulation of genetic evolution of bacterial populations, using Wright-Fisher model, in the presence of complex host contact networks. We present an efficient approach for large-scale simulations over complex host contact networks, using MapReduce on top of Apache Spark and GraphX API. We evaluate the relation between cluster computing power and simulations speedup and include insights on how bacterial population diversity can be affected by mutation and recombination rates, and network topology.

Neutral Theory, Microbial Practice: Challenges in Bacterial Population Genetics

I detail four major open problems in microbial population genetics with direct implications to the study of molecular evolution: the lack of neutral polymorphism, the modeling of promiscuous genetic exchanges, the genetics of ill-defined populations, and the difficulty of untangling selection and demography in the light of these issues. Together with the historical focus on the study of single nucleotide polymorphism and widespread non-random sampling, these problems limit our understanding of the genetic variation in bacterial populations and their adaptive effects. I argue that we need novel theoretical approaches accounting for pervasive selection and strong genetic linkage to better understand microbial evolution.

New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations.

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is considered one of the most harmful plant diseases in the world. Special attention should be paid to R. pseudosolanacearum phylotype I due to its large host range, its worldwide distribution, and its high evolutionary potential. So far, the molecular epidemiology and population genetics of this bacterium are poorly understood. Until now, the genetic structure of the RSSC has been analyzed on the worldwide and regional scales. Emerging questions regarding evolutionary forces in RSSC adaptation to hosts now require genetic markers that are able to monitor RSSC field populations. In this study, we aimed to evaluate the multilocus variable-number tandem-repeat analysis (MLVA) approach for its ability to discriminate genetically close phylotype I strains and for population genetics studies. We developed a new MLVA scheme (MLVA-7) allowing us to genotype 580 R. pseudosolanacearum phylotype I strains extracted from susceptible and resistant hosts and from different habitats (stem, soil, and rhizosphere). Based on specificity, polymorphism, and the amplification success rate, we selected seven fast-evolving variable-number tandem-repeat (VNTR) markers. The newly developed MLVA-7 scheme showed higher discriminatory power than the previously published MLVA-13 scheme when applied to collections sampled from the same location on different dates and to collections from different locations on very small scales. Our study provides a valuable tool for fine-scale monitoring and microevolution-related study of R. pseudosolanacearum phylotype I populations.IMPORTANCE Understanding the evolutionary dynamics of adaptation of plant pathogens to new hosts or ecological niches has become a key point for the development of innovative disease management strategies, including durable resistance. Whereas the molecular mechanisms underlying virulence or pathogenicity changes have been studied thoroughly, the population genetics of plant pathogen adaptation remains an open, unexplored field, especially for plant-pathogenic bacteria. MLVA has become increasingly popular for epidemiosurveillance and molecular epidemiology studies of plant pathogens. However, this method has been used mostly for genotyping and identification on a regional or global scale. In this study, we developed a new MLVA scheme, targeting phylotype I of the soilborne Ralstonia solanacearum species complex (RSSC), specifically to address the bacterial population genetics on the field scale. Such a MLVA scheme, based on fast-evolving loci, may be a tool of choice for field experimental evolution and spatial genetics studies.

A Novel Approach to Helicobacter pylori Pan-Genome Analysis for Identification of Genomic Islands.

Genomes of a given bacterial species can show great variation in gene content and thus systematic analysis of the entire gene repertoire, termed the pan-genome, is important for understanding bacterial intra-species diversity, population genetics, and evolution. Here, we analyzed the pan-genome from 30 completely sequenced strains of the human gastric pathogen Helicobacter pylori belonging to various phylogeographic groups, focusing on 991 accessory (not fully conserved) orthologous groups (OGs). We developed a method to evaluate the mobility of genes within a genome, using the gene order in the syntenically conserved regions as a reference, and classified the 991 accessory OGs into five classes: Core, Stable, Intermediate, Mobile, and Unique. Phylogenetic networks based on the gene content of Core and Stable classes are highly congruent with that created from the concatenated alignment of fully conserved core genes, in contrast to those of Intermediate and Mobile classes, which show quite different topologies. By clustering the accessory OGs on the basis of phylogenetic pattern similarity and chromosomal proximity, we identified 60 co-occurring gene clusters (CGCs). In addition to known genomic islands, including cag pathogenicity island, bacteriophages, and integrating conjugative elements, we identified some novel ones. One island encodes TerY-phosphorylation triad, which includes the eukaryote-type protein kinase/phosphatase gene pair, and components of type VII secretion system. Another one contains a reverse-transcriptase homolog, which may be involved in the defense against phage infection through altruistic suicide. Many of the CGCs contained restriction-modification (RM) genes. Different RM systems sometimes occupied the same (orthologous) locus in the strains. We anticipate that our method will facilitate pan-genome studies in general and help identify novel genomic islands in various bacterial species.

Strong seed-bank effects in bacterial evolution

Bacterial genomes are mosaics with fragments showing distinct phylogenetic origins or even being unrelated to any other genetic information (ORFan genes). Thus the analysis of bacterial population genetics is in large part a collection of explanations for anomalies in relation to classical population genetic models such as the Wright-Fisher model and the Kingman coalescent that do not adequately describe bacterial population genetics, genomics or evolution. The concept of “species” as an evolutionary coherent biological group that is genetically isolated and shares genetic information through recombination among its members cannot be applied to any bacterial group. Recently, a simple probabilistic model considering the role of strong seed-bank effects in population genetics has been proposed by Blath et al. This model suggests the existence of a genetic pool with high diversity that is not subject to classical selection and extinction. We reason that certain bacterial population genetics anomalies could be explained by the prevalence of strong seed-bank effects among bacteria. To address this possibility we analyzed the genome of the bacterium Azotobacter vinelandii and show that genes that code for functions that are essential for the bacterium biology do not have a relation of ancestry with closely related bacteria, or are ORFan genes. The existence of essential genes that are not inherited from the most recent ancestor cannot be explained by classical population genetics models and is irreconcilable with the current view of genes acquired by horizontal transfer as being accessory or adaptive.

TypOn: the microbial typing ontology.

Bacterial identification and characterization at subspecies level is commonly known as Microbial Typing. Currently, these methodologies are fundamental tools in Clinical Microbiology and bacterial population genetics studies to track outbreaks and to study the dissemination and evolution of virulence or pathogenicity factors and antimicrobial resistance. Due to advances in DNA sequencing technology, these methods have evolved to become focused on sequence-based methodologies. The need to have a common understanding of the concepts described and the ability to share results within the community at a global level are increasingly important requisites for the continued development of portable and accurate sequence-based typing methods, especially with the recent introduction of Next Generation Sequencing (NGS) technologies. In this paper, we present an ontology designed for the sequence-based microbial typing field, capable of describing any of the sequence-based typing methodologies currently in use and being developed, including novel NGS based methods. This is a fundamental step to accurately describe, analyze, curate, and manage information for microbial typing based on sequence based typing methods.

Virulence characteristics and multilocus sequence type of Vibrio parahaemolyticus isolated from clinic

Objective To study on virulence characteristics and multilocus sequence type of Vibrio parahaemolyticus isolated from clinic in Beijing Tongren hospital. Methods Total 152 strains of Vibrio parahaemolyticus isolates were collected from diarrheal patients of outpatients in Beijing Tongren Hospital,Capital Medical University from 2009 to 2011.PCR was used to detect hemolysin gene thermo stable direct themolysin gene (tdh), TDH-related hemolysin gene (trh), typeⅢ secretion system 2 (T3SS2α, T3SS2β) and systematic functional gene (toxRS/new,orf8) for pandemic O3∶K6 clone and its derivatives.The genetic features of these strains were determined by multilocus sequence typing (MLST). Results 96% (146/152) VP harbored tdh gene,2% (3/152) VP harbored trh gene and 100% (152/152) VP harbored T3SS2 gene. In this study there were 107 pandemic strains (both tdh and toxRS/new positive and trh negative),38 pathogenic strains (tdh positive and/or trh positive) and 6 nonpathogenic strains (both tdh and trh negative).All nonpathogenic strains harbored systematic functional gene (toxRS/new,orf8).Only one pathogenic strains harbored orf8 gene.One clone harbored all virulence gene.In this study there were 16 sequence types,and ST3 is the pandemic sequence type,including 113 strains,and four new sequence types were found. Conclusions In this study more than 90% Vibrio haemolyticus harbored tdh gene and ST3 was the pandemic sequence type in Beijing.One can get bacterial pathogenic charateristic and population genetics information by virulence gene testing and MLST.(Chin J Lab Med,2013,36:548-552) Key words: Vibrio parahaemolyticus; Virulence; Multilocus sequence typing; Epidemiology; molecular

Optimal Strategy for Rapid Proteome Re-Arrangements in Bacterial Populations

When a bacterial population is subjected to environmental change, its proteome must be re-arranged for resumption of rapid growth. Upon recurrent environmental changes, rapid proteome responses are required for high average growth rate and fitness value, suggesting high selection pressure for rapid proteome re-arrangements in response to environmental changes. What, then, is the best gene expression strategy for rapid proteome change among bacteria?Here, we demonstrate that the optimal solution to the adaptation problem is to direct all gene expression to synthesis of the proteins in the currently rate limiting set of pathways. The result is obtained by reformulating proteome dynamics in terms of flow-coupled networks in which the rate limiting protein components determine the growth rate. This approach provides a universal frame work for the description of the whole proteome and its adaptation to changing external conditions. It also makes it possible to identify optimal solutions to the adaptation problem, which can be experimentally validated. Under optimal conditions the intrinsically non-linear dynamics of these networks becomes simple and linear with respect to the components of the proteome, while the complex dynamics of the concentrations of the signal molecules requires non-linear mathematics: that what is controlled obeys very simple dynamics while the control systems remain non-linear and complex. This allows for simple, analytical solutions to previously intractable dynamic problems. We are therefore optimistic that the present work, focused on the fitness value of bacterial populations, will contribute to quantitative integration of bacterial physiology and population genetics.

Bacterial Population Genetics in Infectious Disease

Bacterial Population Genetics in Infectious Disease

High-Resolution Melting System to Perform Multilocus Sequence Typing of Campylobacter jejuni

Multi-locus sequence typing (MLST) has emerged as the state-of-the-art method for resolving bacterial population genetics but it is expensive and time consuming. We evaluated the potential of high resolution melting (HRM) to identify known MLST alleles of Campylobacter jejuni at reduced cost and time. Each MLST locus was amplified in two or three sub fragments, which were analyzed by HRM. The approach was investigated using 47 C. jejuni isolates, previously characterized by classical MLST, representing isolates from diverse environmental, animal and clinical sources and including the six most prevalent sequence types (ST) and the most frequent alleles. HRM was then applied to a validation set of 84 additional C. jejuni isolates from chickens; 92% of the alleles were resolved in 35 hours of laboratory time and the cost of reagents per isolate was $20 compared with $100 for sequence-based typing. HRM has the potential to complement sequence-based methods for resolving SNPs and to facilitate a wide range of genotyping studies.

Malate dehydrogenase: A useful phylogenetic marker for the genus Aeromonas

The reconstruction of correct genealogies among biological entities, the estimation of the divergence time between organisms or the study of the different events that occur along evolutionary lineages are not always based on suitable genes. For reliable results, it is necessary to look at full-length sequences of genes under stabilizing selection (neutral or purifying) and behaving as good molecular clocks. In bacteria it has been proved that the malate dehydrogenase gene ( mdh) can be used to determine the inter- and intraspecies divergence, and hence this gene constitutes a potential marker for phylogeny and bacterial population genetics. We have sequenced the full-length mdh gene in 36 type and reference strains of Aeromonas. The species grouping obtained in the phylogenetic tree derived from mdh sequences was in agreement with that currently accepted for the genus Aeromonas. The maximum likelihood models applied to our sequences indicated that the mdh gene is highly conserved among the Aeromonas species and the main evolutionary force acting on it is purifying selection. Only two sites under potential diversifying selection were identified (T 108 and S 193). In order to determine if these two residues could have an influence on the MDH structure, we mapped them in a three-dimensional model constructed from the sequence of A. hydrophila using the human mitochondrial MDH as a template. The presence of purifying selection together with the linear relationship between substitutions and gene divergence makes the mdh an excellent candidate gene for a phylogeny of Aeromonas and probably for other bacterial groups.

On the evolution of epistasis III: The haploid case with mutation

Whether interaction between genes is better represented by synergistic or antagonistic epistasis has been a focus of experimental research in bacterial population genetics. Our previous research on evolution of modifiers of epistasis in diploid systems has indicated that the strength of positive or negative epistasis should increase provided linkage disequilibrium is maintained. Here we study a modifier of epistasis in fitness between two loci in a haploid system. Epistasis is modified in the neighborhood of a mutation–selection balance. We show that when linkage in the three-locus system is tight, an increase in the frequency of a modifier allele that induces either more negative or more positive epistasis is possible. Epistasis here can be measured on either an additive or multiplicative scale.

Molecular Identification, Systematics, and Population Structure of Prokaryotes

Prokaryotic systematics is today perhaps one of the most fascinating and fast advancing fields in microbiology. It has become a widely respected scientific discipline, not least because of the application of new ideas and molecular techniques to bacterial identification and classification. Systematics describes the diversity of prokaryotes and classifies or arranges them into groups based on their common similarities. Although a major reason for establishing systematics has been to provide proper identification of bacterial strains, this discipline has developed into a global, multidisciplinary research field providing a phylogenetic framework for various areas within microbiology, including bacterial taxonomy, microbial ecology, clinical microbiology, population genetics, and environmental microbiology.

Optimizing Antibiotic Therapy

I am in full agreement with Dr. Hibbard and his recommendation that antimicrobials should be used in combination to reduce the likelihood of resistance selection (Microbe, May 2006, p. 209). His explanation of bacterial population genetics clearly highlights the issue. Indeed this concept was originally proposed by Paul Ehrlich as early as 1913 (P. Ehrlich, Lancet 1913, p. 455–453) and has been the backbone of antituberculosis therapy for almost 50 years. However the rationale behind the single drug use, often at inappropriate doses for long durations, is not based on science, but most probably on cost or adverse event savings with the logic that in many infections the antibacterials are used to “support” the host defense in combating infections. In many cases this works well to resolve infections but also serves to select the first-step resistant mutants as Dr. Hibbard describes. Karl Drlica's group have espoused a new concept to tackle this issue, the mutant prevention concept wherein they test bacterial populations to determine how much drug it would take to prevent resistance developing by selection of first–step mutants (Y. Dong, X. Zhao, B. N. Kreiswirth, and K. Drlica, Antimicrob. Agents Chemother. 44:2581–2584, 2000). Today, I suspect cost savings and a wary eye towards adverse events are preventing a more complete embrace of the two-or-more-drugs approach. Moreover, if the most pharmacodynamically effective members of a drug class are used (W. M. Scheld, Emerg. Infect. Dis. 9:1–9, 2003) and shown not to be antagonistic to each other then shorter courses maybe possible and overall costs maybe reduced as fewer relapses and follow-up therapy will be required as resistance is kept at bay. As we have few new antibacterials in development we should recall the advice of our predecessors and use this valuable resource in a more thoughtful way.

Neutral microepidemic evolution of bacterial pathogens

Understanding bacterial population genetics is vital for interpreting the response of bacterial populations to selection pressures such as antibiotic treatment or vaccines targeted at only a subset of strains. The evolution of transmissible bacteria occurs by mutation and localized recombination and is influenced by epidemiological as well as molecular processes. We demonstrate that the observed population genetic structure of three important human pathogens, Streptococcus pneumoniae, Neisseria meningitidis, and Staphylococcus aureus, can be explained by using a simple evolutionary model that is based on neutral mutational drift, modulated by recombination, and which incorporates the impact of epidemic transmission in local populations. The predictions of this neutral "microepidemic" model are found to closely fit observed genetic relatedness distributions of bacteria sampled from their natural population, and it provides estimates of the relative rate of recombination that agree well with empirical estimates. The analysis suggests the emergence of neutral bacterial population structure from overlapping microepidemics within clustered host populations and provides insight into the nature and size distribution of these clusters. These findings challenge the assumption that strains of bacterial pathogens differ markedly in relative fitness.