Insights from Genomics into Bacterial Pathogen Populations

Download PDF České info

Bacterial pathogens impose a heavy burden of disease on human populations worldwide. The gravest threats are posed by highly virulent respiratory pathogens, enteric pathogens, and HIV-associated infections. Tuberculosis alone is responsible for the deaths of 1.5 million people annually. Treatment options for bacterial pathogens are being steadily eroded by the evolution and spread of drug resistance. However, population-level whole genome sequencing offers new hope in the fight against pathogenic bacteria. By providing insights into bacterial evolution and disease etiology, these approaches pave the way for novel interventions and therapeutic targets. Sequencing populations of bacteria across the whole genome provides unprecedented resolution to investigate (i) within-host evolution, (ii) transmission history, and (iii) population structure. Moreover, advances in rapid benchtop sequencing herald a new era of real-time genomics in which sequencing and analysis can be deployed within hours in response to rapidly changing public health emergencies. The purpose of this review is to highlight the transformative effect of population genomics on bacteriology, and to consider the prospects for answering abiding questions such as why bacteria cause disease.

Published in the journal: Insights from Genomics into Bacterial Pathogen Populations. PLoS Pathog 8(9): e32767. doi:10.1371/journal.ppat.1002874
Category: Review
doi: https://doi.org/10.1371/journal.ppat.1002874

Summary

Introduction

Bacteria are the most abundant group of organisms, and a major source of human disease and mortality. Bacterial cells account for most of the earth's biomass [1], and the 100 trillion microbial residents of the human body outnumber human cells 10 to 1 [2]. Bacteria that cause pneumonia, diarrhea, and tuberculosis are leading causes of death worldwide [3], [4]. In countries with a low overall burden of infectious disease such as the United States, bacteria are nevertheless responsible for more than 60% of the deaths attributable to communicable disease, with hospital-associated infections, HIV-associated infections, and tuberculosis most prominent (Table 1).

**Tab. 1. Major bacterial causes of death: World and United States.**

Since the introduction of the earliest antibiotics, bacteria have evolved resistance [5]. Treatment options continue to be eroded by the spread of antibiotic resistance [6], not only in countries with advanced health care infrastructure, but globally [7], [8]. However, advances in DNA sequencing capacity offer hope in the fight against pathogenic bacteria because the ability to sequence populations of bacterial genomes is illuminating our understanding of bacterial evolution and virulence. Ultimately these insights will underpin translational research into improved medical practice, drug and vaccine targets, and public health policy.

High-throughput whole genome sequencing (Figure 1) represents a genuine step change for the study of bacterial populations because current approaches are based on the analysis of gene fragments amounting to just a few thousandths the total length of the genome [9], [10]. Population genomics offers unprecedented sensitivity for the detection of rare genetic variants, vastly improved resolution for population studies, and direct sequencing of functionally relevant loci. As a result, it is driving new understanding of within-host evolution, transmission, and population structure. Moreover, the advent of rapid benchtop sequencing is changing the way that microbiology is conducted, signaling a new era of real-time genomics and disseminated collaborative analysis in response to rapidly changing public health emergencies.

**Fig. 1. An example workflow for high-throughput whole genome sequencing in bacteria.**

Within-Host Evolution

Successful colonization of a host is essential to the lifecycle of the pathogen, and the dynamics of the host-pathogen interaction determine the outcome of the interaction, including the severity of disease. DNA/RNA sequencing has greatly advanced the understanding of viral dynamics during infection [11], [12], including the ability to predict disease progression [13], [14], but progress in bacteria has lagged behind, owing to much larger genomes and sparser genetic variation [15]. However, whole genome sequencing in populations of bacteria colonizing individual hosts is shedding new light on the host-pathogen interaction, and the dynamics of bacterial evolution within the host.

At the whole genome scale, genetic variation has been discovered in singly infected hosts colonized by species as disparate as Mycobacterium tuberculosis [16], Salmonella enterica [17], and Staphylococcus aureus [18]. The absolute number of variable sites detected in singly infected hosts is small, frequently fewer than 10 single nucleotide polymorphisms (SNPs), although this varies by species and depends on the number of genomes sequenced and the time elapsed between sampling. Other forms of genetic variation observable at the species level [19] have also been detected, including short insertions and deletions (indels), and variation in the presence or absence of mobile elements such as prophages [16]–[18].

The real-time mutation rate is a key factor in determining the potential for bacterial pathogens to adapt to the host immune system or drug intervention. Traditional estimates of bacterial substitution rates over geological timescales predict fewer than 0.01 mutations per megabase (Mb) per year [20], [21]. Yet laboratory estimates and limited sequencing of longitudinal samples suggest rates 100 or 1,000 times faster [22]–[26]. These estimates have been put to the test by whole genome sequencing, yielding within-host mutation rates ranging from 0.1/Mb/year in Mycobacterium tuberculosis [16] through 2.7/Mb/year in Staphylococcus aureus [27] to 19/Mb/year in Helicobacter pylori [28]. This supports the conclusion that short-term substitution rates in bacteria are several orders of magnitude faster than long-term rates [23], [26], a finding that may be explained by the delayed action of purifying selection [29], [30]. In other words, over longer evolutionary periods the substitution rate depends on selection as well as mutation. It also demonstrates the potential for bacteria to adapt within the host. For example, the discovery that the genome-wide mutation rate in latent tuberculosis infection is similar to that in active disease may explain reports that found treating even latent infections with the antibiotic isoniazid was a risk factor for the emergence of isoniazid resistance [16].

Many bacterial pathogens are common constituents of the body's natural flora [31]. Evolution during colonization may trigger a transition from healthy carriage to invasive disease. For example, 27% of adults carry Staphylococcus aureus asymptomatically in the nose [32], a known risk factor for disease [33]. In a study of one long-term carrier who developed a bloodstream infection, the genomes of invasive bacteria were found to possess an excess of mutations that truncated proteins, including a transcriptional regulator implicated in pathogenicity [27], [34]. Although further work would be needed to establish causality, this demonstrates the potential for loss-of-function mutations to induce radical functional change during colonization.

Unusual patterns of mutation in the genome during colonization may signal adaptive change and reveal mechanisms of virulence or immune evasion. A study of a 16-year outbreak of chronic Burkholderia dolosa infection among cystic fibrosis patients revealed evidence for parallel adaptive evolution [35]. Seventeen genes accrued three or more mutations across the 14 patients, of which a significant excess altered the encoded protein. Some of these mutations affected important phenotypes, including oxygen-dependent gene regulation—which may be pertinent to lung infection—antibiotic resistance, and outer membrane synthesis. Mutations not previously implicated in pathogenesis present novel therapeutic targets.

In hosts colonized multiple times by distinct genotypes, whole genome sequencing affords an opportunity to investigate recombination in vivo. Horizontal gene transfer, also known as recombination, is a fundamental process that generates diversity and facilitates the spread of advantageous genes [36], [37]. A longitudinal study of the highly promiscuous gut pathogen Helicobacter pylori identified recombination events from the clustering of SNP differences introduced by the import of DNA from one strain to another [28]. Surprisingly, multiple fragments of around 400 bases appeared to be simultaneously imported in a span up to 20 kilobases long. This pattern of integration was implied by the results of a similar study in Streptococcus pneumoniae [38], demonstrating the power of whole genome sequencing to illuminate molecular mechanisms.

Detection of Transmission Events

Whole genome sequencing offers unprecedented resolution to distinguish degrees of relatedness among bacterial isolates, and this is a powerful tool for microbial forensics [39]. Genome sequencing complements existing epidemiological tools by providing a means to reconstruct recent chains of transmission, identify sequential acquisition of strains by persistent carriers, and identify cryptic outbreaks that might otherwise go unnoticed.

The superiority of genomics over traditional approaches to molecular epidemiology was demonstrated in a study of Staphylococcus aureus ST-239 [40], a widely disseminated multi-drug resistant clonal lineage dominant in much of Asia. Traditional typing methods offer little discriminatory power for subtyping ST-239, but 5,842 SNPs were discovered by whole genome sequencing, revealing detailed geographical structure within the lineage. Against this backdrop of geographical differentiation, examples of recent intercontinental transmission were evident from the clustering of two isolates from England and Denmark among the Thai group. Moreover, a cluster of five isolates sampled over 11 weeks from adjacent blocks of a Thai hospital differed by just 14 SNPs, providing evidence of recent hospital transmission.

Population genomics offers complementary tools to routine outbreak investigation. Following the discovery of two index cases, an outbreak of Mycobacterium tuberculosis was uncovered in British Columbia using contact tracing and social network questionnaires [41]. The transitivity of the social network and inability to distinguish isolates by traditional genotyping prevented identification of the source. Whole genome sequencing revealed two distinct lineages, ruling out transmission between social contacts infected with discordant types. Further epidemiological investigation intimated a complex scenario in which an increase in crack cocaine usage triggered simultaneous outbreaks that were sustained by key members of a high-risk social network.

In some cases, the direction of transmission may be discernible from patterns of relatedness and associated epidemiological information. In their study of chronic Burkholderia dolosa infection among cystic fibrosis patients [35], the authors used the chronological accumulation of mutations to discriminate donors from recipients in the transmission network. By the same method, they were even able to infer repeated transmission from the airways to the bloodstream within patients (Figure 2). In a persistent Escherichia coli infection of members of a household over three years [42], whole genome sequencing revealed at least six transmission events between family members including at least two zoonotic transmissions to the family dog.

**Fig. 2. Whole genome sequencing reveals within-host evolution and recent transmission between patients.**

Multiple transmission events resulting in serial acquisition by a single host can be distinguished from persistent or relapsing infection using whole genome sequencing. This is useful in infections such as invasive nontyphoidal Salmonella, a common cause of severe and recurring bloodstream infections among HIV-infected adults in Africa. A study of invasive nontyphoidal Salmonella in 14 Malawian patients discriminated recrudescent (i.e., relapsing) infection from multiple infection on the basis of relatedness inferred from genome-wide SNP differences [17]. Recrudescence accounted for 78% of recurring infections, although recrudescence and multiple infection in the same patient was not uncommon.

Historical Patterns of Transmission

In addition to revealing fine-grained genetic structure that is informative about recent transmission, genomics offers unrivalled precision for reconstructing historical patterns of spread. With comprehensive sampling, we can identify the geographical and temporal origin of pandemics and the dominant transmission routes responsible for global dissemination. For example, a study of Yersinia pestis used genome sequencing to assist in the discovery of 933 SNPs subsequently typed in 286 global isolates [43]. Based on the diversity and juxtaposition of isolates close to the root of the tree, the authors concluded that the origin of plague more than 2,600 years ago occurred in or near China.

Understanding the circumstances under which epidemics emerge and take hold may help to manage contemporary threats and prevent future outbreaks. The history of the seventh and current pandemic of Vibrio cholerae was pieced together using population genomics [44]. An analysis of global isolates revealed three partially overlapping waves of pandemic cholera originating in the Bay of Bengal during the 1950s and leading to a succession of geographically restricted epidemics. Each wave was characterized by a particular armory of genetic elements including distinct forms of the cholera toxin and, from the second wave onwards, SXT/R391 integrative and conjugative elements that confer antibiotic resistance.

Sequencing ancient bacteria is a particularly powerful tool for investigating historical transmission. To reconstruct the history of leprosy, SNPs discovered by whole genome sequencing were typed in over 400 isolates of Mycobacterium leprae, including bacteria isolated from skeletal remains recovered from leprosy graveyards in and around Europe [45]. The paleomicrobiological samples resembled modern European isolates, supporting a model in which leprosy arose in East Africa before dispersing east and west by traders along the Silk Road [46]. The provenance of Yersinia pestis was investigated by sequencing bacteria isolated from teeth disinterred from the East Smithfield burial ground for Black Death victims in London [47]. The reconstructed genome closely resembled the most recent common ancestor of modern plague in humans, suggesting that the Black Death was the main historical event antecedent to contemporary plague worldwide (Figure 3).

**Fig. 3. Patterns of historical transmission reconstructed by whole genome sequencing.**

Zoonosis is a major source of emerging infectious disease, with wildlife the most frequent origin [48]. In the United States, leprosy is rare and most infected individuals have a history of foreign residence. Yet a third of cases in Texas and Louisiana had no such explanation [49]. Genome sequencing and SNP typing revealed that a distinctive strain of leprosy was present in 33 wild armadillos and 26 of 29 unexplained human cases, strongly suggesting zoonotic transmission of Mycobacterium leprae from wild armadillos.

Population Structure, Carriage, and Disease

Many bacterial pathogens cause diseases of varying severity, and some cause no disease at all most of the time (e.g., [32], [33]), constituting a normal part of the body flora [2], [31]. Such observations raise the question of why bacteria cause disease and have led to the notion of the accidental pathogen [50], [51]. Virulence may be accounted for by differences between strains and by the expression of genes encoding toxins, adhesins, and drug resistance, often carried by mobile elements (e.g., [52]–[56]). With whole genome sequencing, populations of virulent and avirulent bacteria can be compared to help explain disease from mechanistic and phylogenetic standpoints. Koch's postulates might be revised [57] to cover the discovery of associations between disease and individual genes or alleles, a process that is likely to accelerate rapidly in the 21^st century. A fuller understanding of bacterial population structure may also help predict the effects of interventions such as vaccination.

Bacteria of the same species isolated from patients with different clinical presentations can be compared directly by whole genome sequencing. Streptococcus pyogenes, also known as Group A Streptococcus, causes benign pharyngitis and invasive disease including scarlet fever. A comparison of around 300 isolates from Ontario indicated that invasive bacteria do not form genetically distinct populations. Rather, closely related bacteria may be invasive or pharyngitis-associated, demonstrating that the ability to cause invasive disease is not restricted to specific strains [58], [59]. Evidence for adaptation in genes involved in capsule synthesis and virulence regulation supported a model in which mutation in vivo plays an important role in progression to invasive disease.

Reconstructing the relationships between strains with different clinical manifestations can reveal the evolutionary origins of disease. A study of disparate Clostridium difficile isolates including representative members of the hypervirulent lineages denoted ribotypes 017, 027, and 078 found they were descended from multiple ancestors in the species tree, consistent with a scenario in which virulence evolved several times during evolution [60]. Whole genome sequencing afforded improved clarity for reconstructing the genealogy of Chlamydia trachomatis and appeared to show that strains causing trachoma, an eye infection, were evolutionarily descended from an ancestor causing urogenital disease [61].

The iatrogenic effect of public health intervention on bacterial pathogen populations, such as the introduction of a novel vaccine or the withdrawal of an antibiotic from agricultural use, is of major importance. Vaccine escape in Streptococcus pneumoniae has been a concern since the introduction of the heptavalent conjugate polysaccharide vaccine PCV7 in 2000 as it protects against many, but not all, serotypes. Two genomic studies found evidence for capsular switching, in which hybrid strains normally covered by the vaccine but expressing non-vaccine serotypes arise through recombination [62], [63]. One such strain has quickly established in the United States, spreading westwards from New England [63].

Real-Time Pathogen Genomics

The relentless demand for higher throughput, lower cost DNA sequencing has spurred dramatic advances in the capacity and rapidity of whole genome sequencing. Benchtop sequencers permit real-time applications of genomics by sequencing small batches of bacteria in a matter of hours. The outbreak of cholera on the island of Haiti in October 2010 provided an early example of the potential for real-time genomics [64]. In the wake of the devastating earthquakes of January 2010 that killed 230,000 people, two million residents were displaced from their homes. Cases were first reported on October 19. By July 2011, 419,511 cases and 5,968 deaths had been reported [65]. Initial investigation found many patients had drunk untreated river water. By November 1, 2010 culturing and PFGE confirmed the pathogen as Vibrio cholerae of probable South Asia origin [66]. First-round whole genome sequencing of Haitian isolates began on November 10 and completed within 2 days [64]. Genomic analysis showed they were essentially identical, but distinct from other cholera circulating in Latin America, instead resembling widely circulating Asia strains [64], a finding that was consistent with possible introduction by United Nations Peacekeeping troops dispatched from Nepal following the earthquakes [65].

Real-time genomics may prove particularly valuable in outbreaks involving newly emerged strains. However, processing genomic data in real time poses considerable analytical challenges. In the May 2011 outbreak of Escherichia coli in Germany, a novel crowd-sourcing experiment was trialed that foretells of the potential of real-time genomics to radically alter the way outbreaks are investigated [67]. The large outbreak was unusual in several aspects: high incidence in adults, greatly increased incidence of hemolytic-uremic syndrome, a preponderance of female patients, and a rare, Shiga toxin-producing serotype not previously linked to outbreaks [67], [68]. A first draft of the genome of an isolate sampled on May 17 was completed within 3 days, then released into the public domain, eliciting curiosity-driven analysis by scientists on four continents [67]. Within a week two dozen reports had been filed on a dedicated open-source wiki. Analysis of this and other strains concluded that the outbreak was caused by the acquisition of a Shiga toxin-encoding prophage and a plasmid bearing an extended-spectrum beta-lactamase gene by an ancestral enteroaggregative strain [67], [68]. The striking virulence of the hybrid may be connected to the atypical presence of three SPATE genes, which are implicated in mucosal damage and intestinal colonization [68].

Two studies from hospitals in the United Kingdom have demonstrated the practical advantages of real-time whole genome sequencing as part of routine outbreak investigation and surveillance. Focusing on the most serious health-care-associated pathogens, C. difficile [69] and S. aureus [69], [70], bacterial samples were isolated from suspected outbreaks in four hospitals. The genomes were sequenced and analyzed within 5 working days of culture, confirming the suspected outbreaks of MRSA (methicillin-resistant S. aureus), but demonstrating that the epidemiologically linked cases of C. difficile infection were in fact genetically distinct. Characterization of the repertoire of resistance and toxin genes provided further information relevant to patient management.

Summary and Perspectives

In the future, population genomics will be central to an improved understanding of the epidemiology, etiology, and evolution of bacterial infectious diseases. However, there are obstacles yet to overcome. Pilot studies have demonstrated the potential genomics has for epidemiological investigation [40], [41], [64], [67]–[70], but creative solutions to the problem of integrating complex epidemiological and genomic data are now required. Currently, genome sequencing relies on culture to yield sufficient bacterial DNA and new technologies are needed to overcome this dependency. If the cost of DNA library preparation can be substantially reduced, genomics will come within reach of public health authorities as a tool for routine surveillance. These and other future challenges are discussed in Box 1.

Box 1. Future Challenges for Pathogen Whole Genome Sequencing

High-throughput whole genome sequencing has been demonstrated to be a practical tool for epidemiological and evolutionary investigation of bacterial pathogens, yet the current technology has certain limitations. The challenge for future advances in sequencing technology is to overcome these problems.

Culture. Reliably sequencing the genomes of individual bacteria requires culture to obtain sufficient quantities of concentrated DNA. This takes time and effort, restricts the approach to culturable organisms, and may introduce artifacts such as in vitro mutation and laboratory cross-contamination. Direct sequencing without culture (e.g., [81]) may in the future relinquish this dependency on culture, but metagenomics approaches present additional challenges for bioinformatics and sequence analysis (see, e.g., [82]).

Library preparation. Exponential increases in the capacity of high-throughput sequencers show no sign of abating. In principle, this should allow the cost of bacterial whole genome sequencing to continue to fall. However, the price per genome also depends on the cost of DNA library preparation, comprising both consumables and labor. Advances in automation and throughput will be required to prevent library preparation becoming a bottleneck, and to reduce the cost sufficiently that bacterial genome sequencing becomes affordable for routine surveillance.

Bioinformatics. The development of bespoke bioinformatics pipelines for bacterial whole genome sequencing represents a considerable investment and a complex set of choices from among the many computational methods on offer. Some degree of normalization is required to ease the burden on users of whole genome sequencing, for example hospital microbiology laboratories, and to promote standardized and replicable workflows.

Platforms. High-throughput sequencing technologies yield large quantities of short read sequences but with substantially elevated error rates, compared to conventional capillary sequencing. The details of sequence length and error profile differ in important ways between platforms. In consequence, different results may be obtained when the same sample is sequenced on different platforms. Improved understanding of the error profiles of different architectures combined with efforts towards quantifying uncertainty in the DNA sequences generated will help minimize discrepancies of this kind.

Genome assembly. De novo assembly is used to join together the short reads of DNA generated by the sequencing machines into longer genome fragments, known as contigs. The ultimate goal is to join all the fragments into a single contig representing the whole bacterial chromosome, known as a closed genome. However, variation in the number of reads sequenced from each part of the genome (the depth of coverage), and the existence of repetitive regions, conspire to prevent this. With longer reads, it should be possible to overcome these problems.

Public databases. To accelerate the pace of discovery and assist collaboration between laboratories, well-organized publicly available databases are required from which bacterial genomes are readily downloaded in convenient formats. Raw data are currently available in short read archives (e.g., http://www.ncbi.nlm.nih.gov/sra and http://www.ebi.ac.uk/ena), but with standardization of bioinformatics processing it should become possible to provide pre-processed data which would dramatically reduce the workload for database users.

Improved understanding of disease etiology helps to direct research into therapies. Genomics is a promising tool for investigating the differences between invasive and non-invasive bacteria at the population and within-host levels [27], [59]–[61]. Tools from human genetics may help in this endeavor. Even so, investigations into bacterial population structure are required to assess the feasibility of genome-wide association studies [71]. Understanding the architecture of traits such as virulence would benefit from the development of high-throughput phenotyping assays. RNA sequencing is one such candidate [72], but differences in gene expression in culture and in vivo are a potential impediment to progress.

Population genomics also promises to improve our understanding of bacterial pathogen evolution. The resolution of whole genome sequencing allows precise calibration of evolutionary rates from longitudinal samples within populations and individual hosts [16]–[18], [28], [35], [40], [47]. This permits the origin of new species to be dated, but the discrepancy between short- and long-term rates requires further explanation [23], [30]. Investigating within-host dynamics will help identify the evolutionary mechanisms involved in disease progression [27], [35]. Sequencing populations of pathogens will reveal the prevalence of virulence factors and drug resistance, and the role of mobile elements in their spread [44]. Ultimately, however, we must pinpoint the evolutionary advantages that bacteria gain from inflicting illnesses if we are to fully understand the causes of bacterial disease.

Zdroje

1. WhitmanWB, ColemanDC, WiebeWJ (1998) Prokaryotes: the unseen majority. Proc Natl Acad Sci U S A 95: 6578–6583.

2. QinJ, LiR, RaesJ, ArumugamM, BurgdorfKS, et al. (2010) A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464: 59–65.

3. World Health Organization (2008) The global burden of disease: 2004 update. Available: http://www.who.int/healthinfo/global_burden_disease. Accessed 10 August 2012.

4. World Health Organization (2012) Global invasive bacterial vaccine preventable diseases (IB-VPD) information and surveillance bulletin. Volume 5. Available: http://www.who.int/nuvi/surveillance/resources/en/index.html. Accessed 10 August 2012.

5. DaviesJ, DaviesD (2010) Origins and evolution of antibiotic resistance. Microbiol Mol Biol Rev 74: 417–433.

6. LivermoreDM (2012) Fourteen years in resistance. Int J Antimicrob Agents 39: 283–294.

7. NordmannP, PoirelL, WalshTR, LivermoreDM (2011) The emerging NDM carbapenemases. Trends Microbiol 19: 588–595.

8. SchmiederR (2012) Insights into antibiotic resistance through metagenomics approaches. Future Microbiol 7: 73–89.

9. MaidenMC, BygravesJA, FeilE, MorelliG, RussellJE, et al. (1998) Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95: 3140–3145.

10. UrwinR, MaidenMC (2003) Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol 11: 479–487.

11. GrenfellBT, PybusOG, GogJR, WoodJL, DalyJM, et al. (2004) Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303: 327–332.

12. PybusOG, RambautA (2009) Evolutionary analysis of the dynamics of viral infectious disease. Nat Rev Genet 10: 540–550.

13. ConnorRI, SheridanKE, CeradiniD, ChoeS, LandauNR (1997) Change in coreceptor use correlates with disease progression in HIV-1–infected individuals. J Exp Med 185: 621–628.

14. LemeyP, Kosakovsky PondSL, DrummondAJ, PybusOG, ShapiroB, et al. (2007) Synonymous substitution rates predict HIV disease progression as a result of underlying replication dynamics. PLoS Comput Biol 3: e29 doi:10.1371/journal.pcbi.0030029.

15. DrakeJW, CharlesworthB, CharlesworthD, CrowJF (1998) Rates of spontaneous mutation. Genetics 148: 1667–1686.

16. FordCB, LinPL, ChaseMR, ShahRR, IartchoukO, et al. (2011) Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection. Nat Genet 43: 482–486.

17. OkoroCK, KingsleyRA, QuailMA, KankwatiraAM, FeaseyNA, et al. (2012) High-resolution single nucleotide polymorphism analysis distinguishes recrudescence and reinfection in recurrent invasive nontyphoidal Salmonella typhimurium disease. Clin Infect Dis 54: 955–963.

18. McAdamPR, HolmesA, TempletonKE, FitzgeraldJR (2011) Adaptive evolution of Staphylococcus aureus during chronic endobronchial infection of a cystic fibrosis patient. PLoS ONE 6: e24301 doi:10.1371/journal.pone.0024301.

19. TouzainF, DenamurE, MédigueC, BarbeV, El KarouiM, et al. (2010) Small variable segments constitute a major type of diversity of bacterial genomes at the species level. Genome Biol 11: R45.

20. OchmanH, WilsonAC (1987) Evolution in bacteria: evidence for a universal substitution rate in cellular genomes. J Mol Evol 26: 74–86.

21. MoranNA, MunsonMA, BaumannP, IshikawaH (1993) A molecular clock in endosymbiotic bacteria is calibrated using insect hosts. Proc R Soc Lond B 253: 167–171.

22. LenskiRE, WinkworthCL, RileyMA (2003) Rates of DNA sequence evolution in experimental populations of Escherichia coli during 20,000 generations. J Mol Evol 56: 498–508.

23. OchmanH (2003) Neutral mutations and neutral substitutions in bacterial genomes. Mol Biol Evol 20: 2091–2096.

24. FalushD, KraftC, TaylorNS, CorreaP, FoxJG, et al. (2001) Recombination and mutation during longterm gastric colonization by Helicobacter pylori: estimates of clock rates, recombination size, and minimal age. Proc Natl Acad Sci U S A 98: 15056–15061.

25. Pérez-LosadaM, CrandallKA, ZenilmanJ, ViscidiRP (2007) Temporal trends in gonococcal population genetics in a high prevalence urban community. Infect Genet Evol 7: 271–278.

26. WilsonDJ, GabrielE, LeatherbarrowAJH, CheesbroughJ, GeeS, et al. (2009) Rapid evolution and the importance of recombination to the gastroenteric pathogen Campylobacter jejuni. Mol Biol Evol 26: 385–397.

27. YoungBC, GolubchikT, BattyEM, FungR, Larner-SvennsonH, et al. (2012) Evolutionary dynamics of Staphylococcus aureus during progression from carriage to disease. Proc Natl Acad Sci U S A 109: 4550–4555.

28. KennemannL, DidelotX, AebischerT, KuhnS, DrescherB, et al. (2011) Helicobacter pylori evolution during human infection. Proc Natl Acad Sci U S A 108: 5033–5038.

29. McDonaldJH, KreitmanM (1991) Adaptive protein evolution at the Adh locus in Drosophila. Nature 351: 652–654.

30. RochaEPC, Maynard SmithJ, HurstLD, HoldenMTG, CooperJE, et al. (2006) Comparisons of dN/dS are time dependent for closely related bacterial genomes. J Theor Biol 239: 226–235.

31. TurnbaughPJ, LeyRE, HamadyM, Fraser-LiggettCM, KnightR, et al. (2007) The human microbiome project. Nature 449: 804–810.

32. WertheimHF, MellesDC, VosMC, van LeeuwenW, van BelkumA, et al. (2005) The role of nasal carriage in Staphylococcus aureus infections. Lancet Infect Dis 5: 751–762.

33. von EiffC, BeckerK, MachkaK, StammerH, PetersG (2001) Nasal carriage as a source of Staphylococcus aureus bacteremia. N Engl J Med 344: 11–16.

34. YangJ, TauschekM, Robins-BrowneRM (2011) Control of bacterial virulence by AraC-like regulators that respond to chemical signals. Trends Microbiol 19: 128–135.

35. LiebermanTD, MichelJB, AingaranM, Potter-BynoeG, RouxD, et al. (2011) Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes. Nat Genet 43: 1275–1280.

36. SmithJM, SmithNH, O'RourkeM, SprattBG (1993) How clonal are bacteria? Proc Natl Acad Sci U S A 90: 4384–4388.

37. WeidenbeckJ, CohanFM (2011) Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches. FEMS Microbiol Rev 5: 957–976.

38. HillerNL, AhmedA, PowellE, MartinDP, EutseyR, et al. (2010) Generation of genic diversity among Streptococcus pneumoniae strains via horizontal gene transfer during a chronic polyclonal pediatric infection. PLoS Pathog 6: e1001108 doi:10.1371/journal.ppat.1001108.

39. BudowleB, SchutzerSE, BreezeRG, KeimPS, MorseSA (2011) editors (2011) Microbial forensics, second edition. Elsevier/Academic Press

40. HarrisSR, FeilEJ, HoldenMTG, QuailMA, NickersonEK, et al. (2010) Evolution of MRSA during hospital transmission and intercontinental spread. Science 327: 469–474.

41. GardyJL, JohnstonJC, Ho SuiSJ, CookVJ, ShahL, et al. (2011) Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N Engl J Med 364: 730–739.

42. ReevesPR, LiuB, ZhouZ, LiD, GuoD, et al. (2011) Rates of mutation and host transmission for an Escherichia coli clone over 3 years. PLoS ONE 6: e26907 doi:10.1371/journal.pone.0026907.

43. MorelliG, SongY, MazzoniCJ, EppingerM, RoumagnacP, et al. (2010) Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity. Nat Genet 42: 1140–1143.

44. MutrejaA, KimDW, ThomsonNR, ConnorTR, LeeJH, et al. (2011) Evidence for several waves of global transmission in the seventh cholera pandemic. Nature 477: 462–465.

45. MonotM, HonoréN, GarnierT, ZidaneN, SherafiD, et al. (2009) Comparative genomic and phylogeographic analysis of Mycobacterium leprae. Nat Genet 41: 1282–1291.

46. MonotM, HonoréN, GarnierT, AraozR, CoppéeJY, et al. (2005) On the origin of leprosy. Science 308: 1040–1042.

47. BosKI, SchuenemannVJ, GoldingGB, BurbanoHA, WaglechnerN, et al. (2011) A draft genome of Yersinia pestis from victims of the Black Death. Nature 478: 506–510.

48. JonesKE, PatelNG, LevyMA, StoreygardA, BalkD, et al. (2008) Global trends in emerging infectious diseases. Nature 451: 990–993.

49. TrumanRW, SinghP, SharmaR, BussoP, RougemontJ, et al. (2011) Probable zoonotic leprosy in the Southern United States. N Engl J Med 364: 1626–1633.

50. MoxonER, JansenVAA (2005) Phage variation: understanding the behaviour of an accidental pathogen. Trends Microbiol 13: 563–565.

51. NandiT, OngC, SinghAP, BoddeyJ, AtkinsT, et al. (2010) A genomic survey of positive selection in Burkholderia pseudomallei provides insights into the evolution of accidental virulence. PLoS Pathog 6: e1000845 doi:10.1371/journal.ppat.1000845.

52. PeacockSJ, MooreCE, JusticeA, KantzanouM, StoryL, et al. (2002) Virulent combinations of adhesin and toxin genes in natural populations of Staphylococcus aureus. Infect Immun 70: 4987–4996.

53. MellesDC, GorkinkRF, BoelensHA, SnijdersSV, PeetersJK, et al. (2004) Natural population dynamics and expansion of pathogenic clones of Staphylococcus aureus. J Clin Invest 114: 1732–1740.

54. DiepBA, CarletonHA, ChangRF, SensabaughGF, Perdreau-RemingtonF (2006) Roles of 34 virulence genes in the evolution of hospital- and community-associated strains of methicillin-resistant Staphylococcus aureus. J Infect Dis 193: 1495–1503.

55. LindsayJA, MooreCE, DayNP, PeacockSJ, WitneyAA, et al. (2006) Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. J Bacteriol 188: 669–676.

56. MalachowaN, DeLeoFR (2010) Mobile genetic elements of Staphylococcus aureus. Cell Mol Life Sci 67: 3057–3071.

57. FredricksDN, RelmanDA (1996) Sequence-based identification of microbial pathogens: a reconsideration of Koch's postulates. Clin Microb Rev 9: 18–33.

58. BeresSB, CarrollRK, SheaPR, SitkiewiczI, Martinez-GutierrezJC, et al. (2010) Molecular complexity of successive bacterial epidemics deconvoluted by comparative pathogenomics. Proc Natl Acad Sci U S A 107: 4371–4376.

59. SheaPR, BeresSB, FloresAR, EwbankAL, Gonzalez-LugoJH, et al. (2011) Distinct signatures of diversifying selection revealed by genome analysis of respiratory tract and invasive bacterial populations. Proc Natl Acad Sci U S A 108: 5039–5044.

60. HeM, SebaihiaM, LawleyTD, StablerRA, DawsonLF, et al. (2010) Evolutionary dynamics of Clostridium difficile over short and long time scales. Proc Natl Acad Sci U S A 107: 7527–7532.

61. HarrisSR, ClarkeIN, Seth-SmithHMB, SolomonAW, CutcliffeLT, et al. (2012) Whole-genome analysis of diverse Chlamydia trachomatis strains identifies phylogenetic relationships masked by current clinical typing. Nat Genet 44: 413–420.

62. CroucherNJ, HarrisSR, FraserC, QuailMA, BurtonJ, et al. (2011) Rapid pneumococcal evolution in response to clinical interventions. Science 331: 430–434.

63. GolubchikT, BrueggemannAB, StreetT, Gertz JrRE, SpencerCCA, et al. (2012) Pneumococcal genome sequencing tracks a vaccine escape variant formed through a multi-fragment recombination event. Nat Genet 44: 352–356.

64. ChinCS, SorensonJ, HarrisJB, RobinsWP, CharlesRC, et al. (2011) The origin of the Haitian cholera outbreak strain. N Engl J Med 364: 33–42.

65. TapperoJW, TauxeRV (2011) Lessons learned during public health response to cholera epidemic in Haiti and the Dominican Republic. Emerg Infect Dis 17: 2087–2093.

66. Centers for Disease Control and Prevention. Press release: laboratory test results of cholera outbreak strain in Haiti announced (1st November 2010).

67. RohdeH, QinJ, CuiY, DongfangL, LomanNJ, et al. (2011) Open-source genomic analysis of shiga-toxin-producing E. coli O104:H4. N Engl J Med 365: 718–724.

68. RaskoDA, WebsterDR, SahlJW, BashirA, BoisenN, et al. (2011) Origins of the E. coli strain causing an outbreak of hemolytic-uremic syndrome in Germany. N Engl J Med 365: 709–717.

69. EyreDW, GolubchikT, GordonNC, BowdenR, PiazzaP, et al. (2012) A pilot study of rapid benchtop sequencing of Staphylococcus aureus and Clostridium difficile for outbreak detection and surveillance. BMJ Open 2: e001124.

70. KöserCU, HoldenMTG, EllingtonMJ, CartwrightEJP, BrownNM, et al. (2012) Rapid whole-genome sequencing for investigation of a neonatal MRSA outbreak. N Engl J Med 366: 2267–2275.

71. SpencerCCA, SuZ, DonnellyP, MarchiniJ (2009) Designing genome-wide association studies: sample size, power, imputation and the choice of genotyping chip. PLoS Genet 5: e1000477 doi:10.1371/journal.pgen.1000477.

72. PickrellJK, MarioniJC, PaiAA, DegnerJF, EngelhardtBE, et al. (2010) Understanding mechanisms underlying expression variation with RNA sequencing. Nature 464: 768–772.

73. LiH, DurbinR (2009) Fast and accurate short read alignment with Burrows-Wheeler Transform. Bioinformatics 25: 1754–1760.

74. LunterG, GoodsonM (2011) Stampy: a statistical algorithm for sensitive and fast mapping of Illumina sequence reads. Genome Res 21: 936–939.

75. ZerbinoDR, BirneyE (2008) Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 18: 821–829.

76. IqbalZ, CaccamoM, TurnerI, FlicekP, McVeanG (2012) De novo assembly and genotyping of variants using colored de Bruijn graphs. Nat Genet 44: 226–232.

77. ThompsonJD, GibsonTJ, PlewniakF, JeanmouginF, HigginsDG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res 25: 4876–4882.

78. KurtzS, PhillipyAL, DelcherM, SmootM, ShumwayM, et al. (2004) Versatile and open source software for comparing large genomes. Genome Biol 5: R12.

79. DarlingAE, MauB, PernaNT (2010) ProgressiveMauve: multiple genome alignment with gene gain, loss, and rearrangement. PLoS ONE 5: e11147.

80. Benedictow OJ (2004) The Black Death 1346–1353: the complete history. Boydell Press.

81. DenefVJ, BanfieldJF (2012) In situ evolutionary rate measurements show ecological success of recently emerged bacterial hybrids. Science 336: 462–466.

82. ChenK, PachterL (2005) Bioinformatics for whole-genome shotgun sequencing of microbial communities. PLoS Comput Biol 1: e24 doi:10.1371/journal.pcbi.0010024.

83. World Health Organization (2008) Projections of mortality and burden of disease 2004–2030, baseline scenario. Available: http://www.who.int/healthinfo/global_burden_disease/projections/en/index.html. Accessed 12 April 2012.

84. World Health Organization (2011) World Health Organization Mortality Database, ICD-10. 24 November 2011 update. Available: http://www.who.int/whosis/mort/download/en/index.html. Accessed 12 April 2012.

85. GetahunH, GunnebergC, GranichR, NunnP (2010) HIV infection-associated tuberculosis: the epidemiology and the response. Clin Infect Dis 50: Suppl 3 S201–S207.