Accurate strain identification is essential for anyone working with bacteria. For many species, multilocus sequence typing (MLST) is considered the “gold standard” of typing, but it is traditionally performed in an expensive and time-consuming manner. As the costs of whole-genome sequencing (WGS) continue to decline, it becomes increasingly available to scientists and routine diagnostic laboratories. Currently, the cost is below that of traditional MLST. The new challenges will be how to extract the relevant information from the large amount of data so as to allow for comparison over time and between laboratories. Ideally, this information should also allow for comparison to historical data. We developed a Web-based method for MLST of 66 bacterial species based on WGS data. As input, the method uses short sequence reads from four sequencing platforms or preassembled genomes. Updates from the MLST databases are downloaded monthly, and the best-matching MLST alleles of the specified MLST scheme are found using a BLAST-based ranking method. The sequence type is then determined by the combination of alleles identified. The method was tested on preassembled genomes from 336 isolates covering 56 MLST schemes, on short sequence reads from 387 isolates covering 10 schemes, and on a small test set of short sequence reads from 29 isolates for which the sequence type had been determined by traditional methods. The method presented here enables investigators to determine the sequence types of their isolates on the basis of WGS data. This method is publicly available at www.cbs.dtu.dk/services/MLST .

Laboratory evolution experiments have led to important findings relating organism adaptation and genomic evolution. However, continuous monitoring of long-term evolution has been lacking for natural systems, limiting our understanding of these processes in situ. Here we characterize the evolutionary dynamics of a lineage of a clinically important opportunistic bacterial pathogen,Pseudomonas aeruginosa, as it adapts to the airways of several individual cystic fibrosis patients over 200,000 bacterial generations, and provide estimates of mutation rates of bacteria in a natural environment. In contrast to predictions based on in vitro evolution experiments, we document limited diversification of the evolving lineage despite a highly structured and complex host environment. Notably, the lineage went through an initial period of rapid adaptation caused by a small number of mutations with pleiotropic effects, followed by a period of genetic drift with limited phenotypic change and a genomic signature of negative selection, suggesting that the evolving lineage has reached a major adaptive peak in the fitness landscape. This contrasts with previous findings of continued positive selection from long-term in vitro evolution experiments. The evolved phenotype of the infecting bacteria further suggests that the opportunistic pathogen has transitioned to become a primary pathogen for cystic fibrosis patients.

Genome sequencing of bacterial pathogens has advanced our understanding of their evolution, epidemiology, and response to antibiotic therapy. However, we still have only a limited knowledge of the molecular changes in in vivo evolving bacterial populations in relation to long-term, chronic infections. For example, it remains unclear what genes are mutated to facilitate the establishment of long-term existence in the human host environment, and in which way acquisition of a hypermutator phenotype with enhanced rates of spontaneous mutations influences the evolutionary trajectory of the pathogen. Here we perform a retrospective study of the DK2 clone type of P. aeruginosa isolated from Danish patients suffering from cystic fibrosis (CF), and analyze the genomes of 55 bacterial isolates collected from 21 infected individuals over 38 years. Our phylogenetic analysis of 8,530 mutations in the DK2 genomes shows that the ancestral DK2 clone type spread among CF patients through several independent transmission events. Subsequent to transmission, sub-lineages evolved independently for years in separate hosts, creating a unique possibility to study parallel evolution and identification of genes targeted by mutations to optimize pathogen fitness (pathoadaptive mutations). These genes were related to antibiotic resistance, the cell envelope, or regulatory functions, and we find that the prevalence of pathoadaptive mutations correlates with evolutionary success of co-evolving sub-lineages. The long-term co-existence of both normal and hypermutator populations enabled comparative investigations of the mutation dynamics in homopolymeric sequences in which hypermutators are particularly prone to mutations. We find a positive exponential correlation between the length of the homopolymer and its likelihood to acquire mutations and identify two homopolymer-containing genes preferentially mutated in hypermutators. This homopolymer facilitated differential mutagenesis provides a novel genome-wide perspective on the different evolutionary trajectories of hypermutators, which may help explain their emergence in CF infections. Pseudomonas aeruginosa is the dominating pathogen of chronic airway infections in patients with cystic fibrosis (CF). Although bacterial long-term persistence in CF hosts involves mutation and selection of genetic variants with increased fitness in the CF lung environment, our understanding of the within-host evolutionary processes is limited. Here, we performed a retrospective study of the P. aeruginosa DK2 clone type, which is a transmissible clone isolated from chronically infected Danish CF patients over a period of 38 years. Whole-genome analysis of DK2 isolates enabled a fine-grained reconstruction of the recent evolutionary history of the DK2 lineage and an identification of bacterial genes targeted by mutations to optimize pathogen fitness. The identification of such pathoadaptive genes gives new insight into how the pathogen evolves under the selective pressures of the host immune system and drug therapies. Furthermore, isolates with increased rates of mutation (hypermutator phenotype) emerged in the DK lineage. While this phenotype may accelerate evolution, we also show that hypermutators display differential mutagenesis of certain genes which enable them to follow alternative evolutionary pathways. Overall, our study identifies genes important for bacterial persistence and provides insight into the different mutational mechanisms that govern the adaptive genetic changes.

Microbial population polymorphisms are commonly observed in natural environments, including long-term infected hosts. However, the underlying processes promoting and stabilizing diversity are difficult to unravel and are not well understood. Here, we use chronic infection of cystic fibrosis airways by the opportunistic pathogenPseudomonas aeruginosaas a system for investigating bacterial diversification processes during the course of infection. We analyze clonal bacterial isolates sampled during a 32-year period and map temporal and spatial variations in population diversity to different infection sites within the infected host. We show that the ancestral infecting strain diverged into distinct sublineages, each with their own functional and genomic signatures and rates of adaptation, immediately after initial colonization. The sublineages coexisted in the host for decades, suggesting rapid evolution of stable population polymorphisms. Critically, the observed generation and maintenance of population diversity was the result of partitioning of the sublineages into physically separated niches in the CF airway. The results reveal a complex within-host population structure not previously realized and provide evidence that the heterogeneity of the highly structured and complex host environment promotes the evolution and long-term stability of pathogen population diversity during infection.IMPORTANCEWithin-host pathogen evolution and diversification during the course of chronic infections is of importance in relation to therapeutic intervention strategies, yet our understanding of these processes is limited. Here, we investigate intraclonal population diversity inP. aeruginosaduring chronic airway infections in cystic fibrosis patients. We show the evolution of a diverse population structure immediately after initial colonization, with divergence into multiple distinct sublineages that coexisted for decades and occupied distinct niches. Our results suggest that the spatial heterogeneity in CF airways plays a major role in relation to the generation and maintenance of population diversity and emphasize that a single isolate in sputum may not represent the entire pathogen population in the infected individual. A more complete understanding of the evolution of distinct clonal variants and their distribution in different niches could have positive implications for efficient therapy.