ClonalFrameML: efficient inference of recombination in whole bacterial genomes

Insight Into the Origin and Evolution of the Vibrio parahaemolyticus Pandemic Strain

A strain of Vibrio parahaemolyticus that emerged in 1995 caused the first known pandemic involving this species. This strain comprises clonal autochthonous ocean-dwelling bacteria whose evolution has occurred in the ocean environment. The low sequence diversity in this population enabled the discovery of information on its origin and evolution that has been hidden in bacterial clones that have evolved over a long period. Multilocus sequencing and microarray analysis, together with phylogenetic analysis, of pandemic and pre-pandemic isolates has suggested that the founder clone was an O3:K6 non-pathogenic strain that initially acquired a toxRS/new region and subsequently acquired at least seven novel genomic islands. Sequencing and comparison of whole genomes later confirmed these early observations, and it confirmed that most of the genetic changes occurred via gene conversion involving horizontally transmitted DNA. The highly clonal population rapidly diversified, especially in terms of antigenicity, and 27 serotypes have already been reported. Comparisons of the core genomes derived from the founder clone indicate that there are only a few hundred single-nucleotide variations between isolates. However, when the whole genome is considered (the core plus non-core genome and from any clonal frame), the amount of DNA with a different clonal frame can reach up to 4.2% and the number of single-nucleotide variations can reach several hundred thousand. Altogether, these and previous observations based on multilocus sequence typing, microarray analysis, and whole-genome sequencing indicate the large contribution made by DNA with different clonal genealogy to genome diversification. The evidence also indicates that horizontal gene transfer (HGT) caused the emergence of new pathogens. Furthermore, the extent of HGT seems to depend on the vicissitudes of the life of each bacterium, as exemplified by differences in thousands of base pairs acquired by HGT among almost identical genetic isolates.

Local genes for local bacteria: Evidence of allopatry in the genomes of transatlantic Campylobacter populations

The genetic structure of bacterial populations can be related to geographical locations of isolation. In some species, there is a strong correlation between geographical distance and genetic distance, which can be caused by different evolutionary mechanisms. Patterns of ancient admixture in Helicobacter pylori can be reconstructed in concordance with past human migration, whereas in Mycobacterium tuberculosis it is the lack of recombination that causes allopatric clusters. In Campylobacter, analyses of genomic data and molecular typing have been successful in determining the reservoir host species, but not geographical origin. We investigated biogeographical variation in highly recombining genes to determine the extent of clustering between genomes from geographically distinct Campylobacter populations. Whole‐genome sequences from 294 Campylobacter isolates from North America and the UK were analysed. Isolates from within the same country shared more recently recombined DNA than isolates from different countries. Using 15 UK/American closely matched pairs of isolates that shared ancestors, we identify regions that have frequently and recently recombined to test their correlation with geographical origin. The seven genes that demonstrated the greatest clustering by geography were used in an attribution model to infer geographical origin which was tested using a further 383 UK clinical isolates to detect signatures of recent foreign travel. Patient records indicated that in 46 cases, travel abroad had occurred <2 weeks prior to sampling, and genomic analysis identified that 34 (74%) of these isolates were of a non‐UK origin. Identification of biogeographical markers in Campylobacter genomes will contribute to improved source attribution of clinical Campylobacter infection and inform intervention strategies to reduce campylobacteriosis.

Whole-genome sequencing to investigate a non-clonal melioidosis cluster on a remote Australian island

Melioidosis is a tropical disease caused by the bacterium Burkholderia pseudomallei. Outbreaks are uncommon and can generally be attributed to a single point source and strain. We used whole-genome sequencing to analyse B. pseudomallei isolates collected from an historical 2-year long case cluster that occurred in a remote northern Australian indigenous island community, where infections were previously linked to a contaminated communal water supply. We analysed the genome-wide relatedness of the two most common multilocus sequence types (STs) involved in the outbreak, STs 125 and 126. This analysis showed that although these STs were closely related on a whole-genome level, they demonstrated evidence of multiple recombination events that were unlikely to have occurred over the timeframe of the outbreak. Based on epidemiological and genetic data, we also identified two additional patients not previously associated with this outbreak. Our results confirm the previous hypothesis that a single unchlorinated water source harbouring multiple B. pseudomallei strains was linked to the outbreak, and that increased melioidosis risk in this community was associated with Piper methysticum root (kava) consumption.

SNVPhyl: a single nucleotide variant phylogenomics pipeline for microbial genomic epidemiology

The recent widespread application of whole-genome sequencing (WGS) for microbial disease investigations has spurred the development of new bioinformatics tools, including a notable proliferation of phylogenomics pipelines designed for infectious disease surveillance and outbreak investigation. Transitioning the use of WGS data out of the research laboratory and into the front lines of surveillance and outbreak response requires user-friendly, reproducible and scalable pipelines that have been well validated. Single Nucleotide Variant Phylogenomics (SNVPhyl) is a bioinformatics pipeline for identifying high-quality single-nucleotide variants (SNVs) and constructing a whole-genome phylogeny from a collection of WGS reads and a reference genome. Individual pipeline components are integrated into the Galaxy bioinformatics framework, enabling data analysis in a user-friendly, reproducible and scalable environment. We show that SNVPhyl can detect SNVs with high sensitivity and specificity, and identify and remove regions of high SNV density (indicative of recombination). SNVPhyl is able to correctly distinguish outbreak from non-outbreak isolates across a range of variant-calling settings, sequencing-coverage thresholds or in the presence of contamination. SNVPhyl is available as a Galaxy workflow, Docker and virtual machine images, and a Unix-based command-line application. SNVPhyl is released under the Apache 2.0 license and available at http://snvphyl.readthedocs.io/ or at https://github.com/phac-nml/snvphyl-galaxy.

A novel isolation protocol and probe-based RT-PCR for diagnosis of gastric infections with the zoonotic pathogen Helicobacter suis

BACKGROUND

Helicobacter suis is a very fastidious microorganism associated with gastritis, gastric ulcers, and mucosa-associated lymphoid tissue lymphoma in humans. In vitro isolation of this agent from human patients has so far been unsuccessful.

MATERIALS AND METHODS

A probe-based real-time PCR (RT-PCR) for the rapid detection of H. suis in gastric biopsies was developed. Secondly, a mouse-passage-based protocol was optimized for isolation of low numbers of viable H. suis bacteria. Mice were inoculated with different numbers of viable H. suis (10² -10⁸ ) and kept for 4 weeks to allow multiplication of this pathogen.

RESULTS

The probe-based real-time PCR (RT-PCR) exhibited a high degree of diagnostic specificity and analytical sensitivity, high linear correlations (r² between 0.995 and 0.999), and high amplification efficiencies (>90%) for H. suis. No cross-reactivity was detected with human, porcine, non-human primate, and murine DNA nor with DNA from other bacteria including Helicobacter spp. and Campylobacter spp. H. suis was successfully re-isolated from the stomach of mice inoculated with at least 10⁴ viable H. suis, using a biphasic medium (pH 5), consisting of Brucella agar with Brucella broth on top, both supplemented with vitox supplement, Campylobacter-selective supplement, amphotericin (5 μg/mL), HCl (0.05%), fetal bovine serum (20%), and linezolid (5 μg/mL). Linezolid was necessary to inhibit proliferation of contaminants, including lactobacilli.

CONCLUSION

The methods described above can be implemented for detection or isolation of H. suis from human gastric biopsies.

Population Genomics of Legionella longbeachae and Hidden Complexities of Infection Source Attribution

Legionella longbeachae is the primary cause of legionellosis in Australasia and Southeast Asia and an emerging pathogen in Europe and the United States; however, our understanding of the population diversity of L. longbeachae from patient and environmental sources is limited. We analyzed the genomes of 64 L. longbeachae isolates, of which 29 were from a cluster of legionellosis cases linked to commercial growing media in Scotland in 2013 and 35 were non-outbreak-associated isolates from Scotland and other countries. We identified extensive genetic diversity across the L. longbeachae species, associated with intraspecies and interspecies gene flow, and a wide geographic distribution of closely related genotypes. Of note, we observed a highly diverse pool of L. longbeachae genotypes within compost samples that precluded the genetic establishment of an infection source. These data represent a view of the genomic diversity of L. longbeachae that will inform strategies for investigating future outbreaks.

Contribution to Clostridium Difficile Transmission of Symptomatic Patients With Toxigenic Strains Who Are Fecal Toxin Negative

Background

The role of symptomatic patients who are toxigenic strain positive (TS+) but fecal toxin negative (FT-) in transmission of Clostridium difficile is currently unknown.

Methods

We investigated the contribution of symptomatic TS+/FT- and TS+/FT+ patients in C. difficile transmission in 2 UK regions. From 2-step testing, all glutamate dehydrogenase (GDH)-positive specimens, regardless of fecal toxin result, from Oxford (April 2012 through April 2013) and Leeds (July 2012 through April 2013) microbiology laboratories underwent culture and whole-genome sequencing (WGS), using WGS to identify toxigenic strains. Plausible sources for each TS+/FT+ case, including TS+/FT- and TS+/FT+ patients, were determined using WGS, with and without hospital admission data.

Results

A total of 1447 of 12772 (11%) fecal samples were GDH positive, 866 of 1447 (60%) contained toxigenic C. difficile, and fecal toxin was detected in 511 of 866 (59%), representing 235 Leeds and 191 Oxford TS+/FT+ cases. TS+/FT+ cases were 3 times more likely to be plausibly acquired from a previous TS+/FT+ case than a TS+/FT- patient. Fifty-one of 265 (19%) TS+/FT+ cases diagnosed >3 months into the study were genetically related (≤2 single-nucleotide polymorphisms) to ≥1 previous TS+/FT+ case or TS+/FT- patient: 27 (10%) to only TS+/FT+ cases, 9 (3%) to only TS+/FT- patients, and 15 (6%) to both. Only 10 of 265 (4%) were genetically related to a previous TS+/FT+ or TS+/FT- patient and shared the same ward simultaneously or within 28 days.

Conclusions

Symptomatic TS+/FT- patients were a source of C. difficile transmission, although they accounted for less onward transmission than TS+/FT+ cases. Although transmission from symptomatic patients with either fecal toxin status accounted for a low overall proportion of new cases, both groups should be infection control targets.

The Bacterial Sequential Markov Coalescent

Bacteria can exchange and acquire new genetic material from other organisms directly and via the environment. This process, known as bacterial recombination, has a strong impact on the evolution of bacteria, for example, leading to the spread of antibiotic resistance across clades and species, and to the avoidance of clonal interference. Recombination hinders phylogenetic and transmission inference because it creates patterns of substitutions (homoplasies) inconsistent with the hypothesis of a single evolutionary tree. Bacterial recombination is typically modeled as statistically akin to gene conversion in eukaryotes, i.e., using the coalescent with gene conversion (CGC). However, this model can be very computationally demanding as it needs to account for the correlations of evolutionary histories of even distant loci. So, with the increasing popularity of whole genome sequencing, the need has emerged for a faster approach to model and simulate bacterial genome evolution. We present a new model that approximates the coalescent with gene conversion: the bacterial sequential Markov coalescent (BSMC). Our approach is based on a similar idea to the sequential Markov coalescent (SMC)-an approximation of the coalescent with crossover recombination. However, bacterial recombination poses hurdles to a sequential Markov approximation, as it leads to strong correlations and linkage disequilibrium across very distant sites in the genome. Our BSMC overcomes these difficulties, and shows a considerable reduction in computational demand compared to the exact CGC, and very similar patterns in simulated data. We implemented our BSMC model within new simulation software FastSimBac. In addition to the decreased computational demand compared to previous bacterial genome evolution simulators, FastSimBac provides more general options for evolutionary scenarios, allowing population structure with migration, speciation, population size changes, and recombination hotspots. FastSimBac is available from https://bitbucket.org/nicofmay/fastsimbac, and is distributed as open source under the terms of the GNU General Public License. Lastly, we use the BSMC within an Approximate Bayesian Computation (ABC) inference scheme, and suggest that parameters simulated under the exact CGC can correctly be recovered, further showcasing the accuracy of the BSMC. With this ABC we infer recombination rate, mutation rate, and recombination tract length of Bacillus cereus from a whole genome alignment.

Ecological and Evolutionary Insights into Xanthomonas citri Pathovar Diversity

Citrus canker, caused by Xanthomonas citri pv. citri, is a serious disease of citrus plants worldwide. Earlier phylogenetic studies using housekeeping genes revealed that X. citri pv. citri is related to many other pathovars, which can be collectively referred as Xanthomonas citri pathovars (XCPs). From the present study, we report the genome sequences of 18 XCPs and compared them with four XCPs available in the public domain. In a tree based on phylogenomic marker genes, all the XCPs form a monophyletic cluster, suggesting their origin from a common ancestor. Phylogenomic analysis using the type strain further established that all the XCPs belong to one species. Clonal analysis of the core genome revealed the presence of two major lineages within this monophyletic cluster consisting of some clonal variants. Incidentally, the majority of these XCPs were first noticed in India, corroborating their clonal relationship and their common origin. Comparative analysis revealed an open pan-genome and the role of interstrain genomic flux of these XCPs since their diversification from a common ancestor. Even though there are wide variations in type III gene effectomes, we identified three core effectors which can be valuable in resistance-breeding programs. Overall, genomic examination of ecological relatives allowed us to dissect the tremendous genomic potential of X. citri species to rapidly evolve into specialized strains infecting diverse crop plants.IMPORTANCE Host specialization is one of the characteristic features of highly evolved pathogens such as the Xanthomonas group of phytopathogenic bacteria. Since the hosts involve staple crops and economically important fruits such as citrus, detailed understanding of the diversity and evolution of such strains infecting diverse plants is important for quarantine purposes. In the present study, we carried out genomic investigation of members of a phylogenetically and ecologically defined group of Xanthomonas strains pathogenic to diverse plants, including citrus. This group includes the oldest Xanthomonas pathovars and also recently emerged pathovars in a particular country where they are endemic. Our high-throughput genomic study has provided novel insights into the evolution of a unique lineage consisting of serious pathogens and their ecological relatives, suggesting the nature, scope, and pattern of rapid and recent diversification. Further, from the level of species to that of clonal variants, the study revealed interesting genomic patterns in diversification of a Xanthomonas lineage and perhaps will inspire careful study of the host range of the included pathovars.

Target-Specific Assay for Rapid and Quantitative Detection of Mycobacterium chimaera DNA

Mycobacterium chimaera is an opportunistic environmental mycobacterium belonging to the Mycobacterium avium-M. intracellulare complex. Although most commonly associated with pulmonary disease, there has been growing awareness of invasive M. chimaera infections following cardiac surgery. Investigations suggest worldwide spread of a specific M. chimaera clone, associated with contaminated hospital heater-cooler units used during the surgery. Given the global dissemination of this clone, its potential to cause invasive disease, and the laboriousness of current culture-based diagnostic methods, there is a pressing need to develop rapid and accurate diagnostic assays specific for M. chimaera Here, we assessed 354 mycobacterial genome sequences and confirmed that M. chimaera is a phylogenetically coherent group. In silico comparisons indicated six DNA regions present only in M. chimaera We targeted one of these regions and developed a TaqMan quantitative PCR (qPCR) assay for M. chimaera with a detection limit of 100 CFU/ml in whole blood spiked with bacteria. In vitro screening against DNA extracted from 40 other mycobacterial species and 22 bacterial species from 21 diverse genera confirmed the in silico-predicted specificity for M. chimaera Screening 33 water samples from heater-cooler units with this assay highlighted the increased sensitivity of PCR compared to culture, with 15 of 23 culture-negative samples positive by M. chimaera qPCR. We have thus developed a robust molecular assay that can be readily and rapidly deployed to screen clinical and environmental specimens for M. chimaera.

Origin and Evolution of the Kiwifruit Canker Pandemic

Recurring epidemics of kiwifruit (Actinidia spp.) bleeding canker disease are caused by Pseudomonas syringae pv. actinidiae (Psa). In order to strengthen understanding of population structure, phylogeography, and evolutionary dynamics, we isolated Pseudomonas from cultivated and wild kiwifruit across six provinces in China. Based on the analysis of 80 sequenced Psa genomes, we show that China is the origin of the pandemic lineage but that strain diversity in China is confined to just a single clade. In contrast, Korea and Japan harbor strains from multiple clades. Distinct independent transmission events marked introduction of the pandemic lineage into New Zealand, Chile, Europe, Korea, and Japan. Despite high similarity within the core genome and minimal impact of within-clade recombination, we observed extensive variation even within the single clade from which the global pandemic arose.

Genomic analysis of ST88 community-acquired methicillin resistant Staphylococcus aureus in Ghana

Background

The emergence and evolution of community-acquired methicillin resistant Staphylococcus aureus (CA-MRSA) strains in Africa is poorly understood. However, one particular MRSA lineage called ST88, appears to be rapidly establishing itself as an “African” CA-MRSA clone. In this study, we employed whole genome sequencing to provide more information on the genetic background of ST88 CA-MRSA isolates from Ghana and to describe in detail ST88 CA-MRSA isolates in comparison with other MRSA lineages worldwide.

Methods

We first established a complete ST88 reference genome (AUS0325) using PacBio SMRT sequencing. We then used comparative genomics to assess relatedness among 17 ST88 CA-MRSA isolates recovered from patients attending Buruli ulcer treatment centres in Ghana, three non-African ST88s and 15 other MRSA lineages.

Results

We show that Ghanaian ST88 forms a discrete MRSA lineage (harbouring SCCmec-IV [2B]). Gene content analysis identified five distinct genomic regions enriched among ST88 isolates compared with the other S. aureus lineages. The Ghanaian ST88 isolates had only 658 core genome SNPs and there was no correlation between phylogeny and geography, suggesting the recent spread of this clone. The lineage was also resistant to multiple classes of antibiotics including β-lactams, tetracycline and chloramphenicol.

Discussion

This study reveals that S. aureus ST88-IV is a recently emerging and rapidly spreading CA-MRSA clone in Ghana. The study highlights the capacity of small snapshot genomic studies to provide actionable public health information in resource limited settings. To our knowledge this is the first genomic assessment of the ST88 CA-MRSA clone.

Patterns and processes of Mycobacterium bovis evolution revealed by phylogenomic analyses

Mycobacterium bovis is an important animal pathogen worldwide that parasitizes wild and domesticated vertebrate livestock as well as humans. A comparison of the five M. bovis complete genomes from the United Kingdom, South Korea, Brazil, and the United States revealed four novel large-scale structural variations of at least 2,000 bp. A comparative phylogenomic study including 2,483 core genes of 38 taxa from eight countries showed conflicting phylogenetic signal among sites. By minimizing this effect, we obtained a tree that better agrees with sampling locality. Results supported a relatively basal position of African strains (all isolated from Homo sapiens ), confirming that Africa was an important region for early diversification and that humans were one of the earliest hosts. Selection analyses revealed that functional categories such as “Lipid transport and metabolism,” “Cell cycle control, cell division, chromosome partitioning” and “Cell motility” were significant for the evolution of the group, besides other categories previously described, showing importance of genes associated with virulence and cholesterol metabolism in the evolution of M. bovis . PE/PPE genes, many of which are known to be associated with virulence, were major targets for large-scale polymorphisms, homologous recombination, and positive selection, evincing for the first time a plethora of evolutionary forces possibly contributing to differential adaptability in M. bovis . By assuming different priors, US strains originated and started to diversify around 150–5,210 ya. By further analyzing the largest set of US genomes to date (76 in total), obtained from 14 host species, we detected that hosts were not clustered in clades (except for a few cases), with some faster-evolving strains being detected, suggesting fast and ongoing reinfections across host species, and therefore, the possibility of new bovine tuberculosis outbreaks.

Population structure of multidrug resistant Klebsiella oxytoca within hospitals across the UK and Ireland identifies sharing of virulence and resistance genes with K. pneumoniae

Klebsiella oxytoca, a member of the Enterobacteriaceae, is a gram-negative pathogenic bacterium of environmental origin, which can cause infection in healthcare settings. Outbreaks of multidrug-resistant K. oxytoca infection have been increasingly reported in hospitalized patients. Despite the growing importance of this pathogen, there is limited knowledge about the population structure and epidemiology of antimicrobial resistant K. oxytoca. We investigated the population structure and genomic basis of antimicrobial resistance of 41 multidrug resistant K. oxytoca isolates recovered from bloodstream infections across the UK and Ireland. Our results show that K. oxytoca has a highly diverse population, which is composed of several distinct clades, and we identified one recent expansion of a clone within our dataset. Although the K. oxytoca genomes are clearly distinct from the genomes of a global collection of Klebsiella pneumoniae complex, pre-dominantly composed of K. pneumoniae, we found evidence for sharing of core genes through recombination, as well as the exchange of accessory antimicrobial resistance and virulence factor genes between the species. Our findings also suggest that the different K. oxytoca clades have acquired antimicrobial resistance and virulence factor genes independently. This highlights the clinical and therapeutic importance of genetic flexibility in K. oxytoca and the relevance of this in its role as an opportunistic pathogen.

Inferring Ancestral Recombination Graphs from Bacterial Genomic Data

Homologous recombination is a central feature of bacterial evolution, yet it confounds traditional phylogenetic methods. While a number of methods specific to bacterial evolution have been developed, none of these permit joint inference of a bacterial recombination graph and associated parameters. In this article, we present a new method which addresses this shortcoming. Our method uses a novel Markov chain Monte Carlo algorithm to perform phylogenetic inference under the ClonalOrigin model. We demonstrate the utility of our method by applying it to ribosomal multilocus sequence typing data sequenced from pathogenic and nonpathogenic Escherichia coli serotype O157 and O26 isolates collected in rural New Zealand. The method is implemented as an open source BEAST 2 package, Bacter, which is available via the project web page at http://tgvaughan.github.io/bacter.

Effects of control interventions on Clostridium difficile infection in England: an observational study

BACKGROUND

The control of Clostridium difficile infections is an international clinical challenge. The incidence of C difficile in England declined by roughly 80% after 2006, following the implementation of national control policies; we tested two hypotheses to investigate their role in this decline. First, if C difficile infection declines in England were driven by reductions in use of particular antibiotics, then incidence of C difficile infections caused by resistant isolates should decline faster than that caused by susceptible isolates across multiple genotypes. Second, if C difficile infection declines were driven by improvements in hospital infection control, then transmitted (secondary) cases should decline regardless of susceptibility.

METHODS

Regional (Oxfordshire and Leeds, UK) and national data for the incidence of C difficile infections and antimicrobial prescribing data (1998-2014) were combined with whole genome sequences from 4045 national and international C difficile isolates. Genotype (multilocus sequence type) and fluoroquinolone susceptibility were determined from whole genome sequences. The incidence of C difficile infections caused by fluoroquinolone-resistant and fluoroquinolone-susceptible isolates was estimated with negative-binomial regression, overall and per genotype. Selection and transmission were investigated with phylogenetic analyses.

FINDINGS

National fluoroquinolone and cephalosporin prescribing correlated highly with incidence of C difficile infections (cross-correlations >0·88), by contrast with total antibiotic prescribing (cross-correlations <0·59). Regionally, C difficile decline was driven by elimination of fluoroquinolone-resistant isolates (approximately 67% of Oxfordshire infections in September, 2006, falling to approximately 3% in February, 2013; annual incidence rate ratio 0·52, 95% CI 0·48-0·56 vs fluoroquinolone-susceptible isolates: 1·02, 0·97-1·08). C difficile infections caused by fluoroquinolone-resistant isolates declined in four distinct genotypes (p<0·01). The regions of phylogenies containing fluoroquinolone-resistant isolates were short-branched and geographically structured, consistent with selection and rapid transmission. The importance of fluoroquinolone restriction over infection control was shown by significant declines in inferred secondary (transmitted) cases caused by fluoroquinolone-resistant isolates with or without hospital contact (p<0·0001) versus no change in either group of cases caused by fluoroquinolone-susceptible isolates (p>0·2).

INTERPRETATION

Restricting fluoroquinolone prescribing appears to explain the decline in incidence of C difficile infections, above other measures, in Oxfordshire and Leeds, England. Antimicrobial stewardship should be a central component of C difficile infection control programmes.

FUNDING

UK Clinical Research Collaboration (Medical Research Council, Wellcome Trust, National Institute for Health Research); NIHR Oxford Biomedical Research Centre; NIHR Health Protection Research Unit on Healthcare Associated Infection and Antimicrobial Resistance (Oxford University in partnership with Public Health England [PHE]), and on Modelling Methodology (Imperial College, London in partnership with PHE); and the Health Innovation Challenge Fund.

Evolutionary origins of the emergent ST796 clone of vancomycin resistant Enterococcus faecium

From early 2012, a novel clone of vancomycin resistant Enterococcus faecium (assigned the multi locus sequence type ST796) was simultaneously isolated from geographically separate hospitals in south eastern Australia and New Zealand. Here we describe the complete genome sequence of Ef_aus0233, a representative ST796 E. faecium isolate. We used PacBio single molecule real-time sequencing to establish a high quality, fully assembled genome comprising a circular chromosome of 2,888,087 bp and five plasmids. Comparison of Ef_aus0233 to other E. faecium genomes shows Ef_aus0233 is a member of the epidemic hospital-adapted lineage and has evolved from an ST555-like ancestral progenitor by the accumulation or modification of five mosaic plasmids and five putative prophage, acquisition of two cryptic genomic islands, accrued chromosomal single nucleotide polymorphisms and a 80 kb region of recombination, also gaining Tn1549 and Tn916, transposons conferring resistance to vancomycin and tetracycline respectively. The genomic dissection of this new clone presented here underscores the propensity of the hospital E. faecium lineage to change, presumably in response to the specific conditions of hospital and healthcare environments.

Population genomic insights into variation and evolution of Xanthomonas oryzae pv. oryzae

Xanthomonas oryzae pv. oryzae ( Xoo) is a serious pathogen of rice causing bacterial leaf blight disease. Resistant varieties and breeding programs are being hampered by the emergence of highly virulent strains. Herein we report population based whole genome sequencing and analysis of 100 Xoo strains from India. Phylogenomic analysis revealed the clustering of Xoo strains from India along with other Asian strains, distinct from African and US Xo strains. The Indian Xoo population consists of a major clonal lineage and four minor but highly diverse lineages. Interestingly, the variant alleles, gene clusters and highly pathogenic strains are primarily restricted to minor lineages L-II to L-V and in particularly to lineage L-III. We could also find the association of an expanded CRISPR cassette and a highly variant LPS gene cluster with the dominant lineage. Molecular dating revealed that the major lineage, L-I is youngest and of recent origin compared to remaining minor lineages that seems to have originated much earlier in the past. Further, we were also able to identify core effector genes that may be helpful in efforts towards building durable resistance against this pathogen.

Genome Analysis of Clostridium difficile PCR Ribotype 014 Lineage in Australian Pigs and Humans Reveals a Diverse Genetic Repertoire and Signatures of Long-Range Interspecies Transmission

Clostridium difficile PCR ribotype (RT) 014 is well-established in both human and porcine populations in Australia, raising the possibility that C. difficile infection (CDI) may have a zoonotic or foodborne etiology. Here, whole genome sequencing and high-resolution core genome phylogenetics were performed on a contemporaneous collection of 40 Australian RT014 isolates of human and porcine origin. Phylogenies based on MLST (7 loci, STs 2, 13, and 49) and core orthologous genes (1260 loci) showed clustering of human and porcine strains indicative of very recent shared ancestry. Core genome single nucleotide variant (SNV) analysis found 42% of human strains showed a clonal relationship (separated by ≤2 SNVs in their core genome) with one or more porcine strains, consistent with recent inter-host transmission. Clones were spread over a vast geographic area with 50% of the human cases occurring without recent healthcare exposure. These findings suggest a persistent community reservoir with long-range dissemination, potentially due to agricultural recycling of piggery effluent. We also provide the first pan-genome analysis for this lineage, characterizing its resistome, prophage content, and in silico virulence potential. The RT014 is defined by a large "open" pan-genome (7587 genes) comprising a core genome of 2296 genes (30.3% of the total gene repertoire) and an accessory genome of 5291 genes. Antimicrobial resistance genotypes and phenotypes varied across host populations and ST lineages and were characterized by resistance to tetracycline [tetM, tetA(P), tetB(P) and tetW], clindamycin/erythromycin (ermB), and aminoglycosides (aph3-III-Sat4A-ant6-Ia). Resistance was mediated by clinically important mobile genetic elements, most notably Tn6194 (harboring ermB) and a novel variant of Tn5397 (harboring tetM). Numerous clinically important prophages (Siphoviridae and Myoviridae) were identified as well as an uncommon accessory gene regulator locus (agr3). Conservation in the pathogenicity locus and S-layer correlated with ST affiliation, further extending the concept of clonal C. difficile lineages. This study provides novel insights on the genetic variability and strain relatedness of C. difficile RT014, a lineage of emerging One Health importance. Ongoing molecular and genomic surveillance of strains in humans, animals, food, and the environment is imperative to identify opportunities to reduce the overall CDI burden.

Supra-operonic clusters of functionally related genes (SOCs) are a source of horizontal gene co-transfers

Adaptation of bacteria occurs predominantly via horizontal gene transfer (HGT). While it is widely recognized that horizontal acquisitions frequently encompass multiple genes, it is unclear what the size distribution of successfully transferred DNA segments looks like and what evolutionary forces shape this distribution. Here, we identified 1790 gene family pairs that were consistently co-gained on the same branches across a phylogeny of 53 E. coli strains. We estimated a lower limit of their genomic distances at the time they were transferred to their host genomes; this distribution shows a sharp upper bound at 30 kb. The same gene-pairs can have larger distances (up to 70 kb) in other genomes. These more distant pairs likely represent recent acquisitions via transduction that involve the co-transfer of excised prophage genes, as they are almost always associated with intervening phage-associated genes. The observed distribution of genomic distances of co-transferred genes is much broader than expected from a model based on the co-transfer of genes within operons; instead, this distribution is highly consistent with the size distribution of supra-operonic clusters (SOCs), groups of co-occurring and co-functioning genes that extend beyond operons. Thus, we propose that SOCs form a basic unit of horizontal gene transfer.

Reconstructing the Ancestral Relationships Between Bacterial Pathogen Genomes

Following recent developments in DNA sequencing technology, it is now possible to sequence hundreds of whole genomes from bacterial isolates at relatively low cost. Analyzing this growing wealth of genomic data in terms of ancestral relationships can reveal many interesting aspects of the evolution, ecology, and epidemiology of bacterial pathogens. However, reconstructing the ancestry of a sample of bacteria remains challenging, especially for the majority of species where recombination is frequent. Here, we review and describe the computational techniques currently available to infer ancestral relationships, including phylogenetic methods that either ignore or account for the effect of recombination, as well as model-based and model-free phylogeny-independent approaches.

Insidious Risk of Severe Mycobacterium chimaera Infection in Cardiac Surgery Patients

BACKGROUND

An urgent UK investigation was launched to assess risk of invasive Mycobacterium chimaera infection in cardiothoracic surgery and a possible association with cardiopulmonary bypass heater-cooler units following alerts in Switzerland and The Netherlands.

METHODS

Parallel investigations were pursued: (1) identification of cardiopulmonary bypass-associated M. chimaera infection through national laboratory and hospital admissions data linkage; (2) cohort study to assess patient risk; (3) microbiological and aerobiological investigations of heater-coolers in situ and under controlled laboratory conditions; and (4) whole-genome sequencing of clinical and environmental isolates.

RESULTS

Eighteen probable cases of cardiopulmonary bypass-associated M. chimaera infection were identified; all except one occurred in adults. Patients had undergone valve replacement in 11 hospitals between 2007 and 2015, a median of 19 months prior to onset (range, 3 months to 5 years). Risk to patients increased after 2010 from <0.2 to 1.65 per 10000 person-years in 2013, a 9-fold rise for infections within 2 years of surgery (rate ratio, 9.08 [95% CI, 1.81-87.76]). Endocarditis was the most common presentation (n = 11). To date, 9 patients have died. Investigations identified aerosol release through breaches in heater-cooler tanks. Mycobacterium chimaera and other pathogens were recovered from water and air samples. Phylogenetic analysis found close clustering of strains from probable cases.

CONCLUSIONS

We identified low but escalating risk of severe M. chimaera infection associated with heater-coolers with cases in a quarter of cardiothoracic centers. Our investigations strengthen etiological evidence for the role of heater-coolers in transmission and raise the possibility of an ongoing, international point-source outbreak. Active management of heater-coolers and heightened clinical awareness are imperative given the consequences of infection.

Identification of Klebsiella capsule synthesis loci from whole genome data

Klebsiella pneumoniae is a growing cause of healthcare-associated infections for which multi-drug resistance is a concern. Its polysaccharide capsule is a major virulence determinant and epidemiological marker. However, little is known about capsule epidemiology since serological typing is not widely accessible and many isolates are serologically non-typeable. Molecular typing techniques provide useful insights, but existing methods fail to take full advantage of the information in whole genome sequences. We investigated the diversity of the capsule synthesis loci (K-loci) among 2503 K. pneumoniae genomes. We incorporated analyses of full-length K-locus nucleotide sequences and also clustered protein-encoding sequences to identify, annotate and compare K-locus structures. We propose a standardized nomenclature for K-loci and present a curated reference database. A total of 134 distinct K-loci were identified, including 31 novel types. Comparative analyses indicated 508 unique protein-encoding gene clusters that appear to reassort via homologous recombination. Extensive intra- and inter-locus nucleotide diversity was detected among the wzi and wzc genes, indicating that current molecular typing schemes based on these genes are inadequate. As a solution, we introduce Kaptive, a novel software tool that automates the process of identifying K-loci based on full locus information extracted from whole genome sequences (https://github.com/katholt/Kaptive). This work highlights the extensive diversity of Klebsiella K-loci and the proteins that they encode. The nomenclature, reference database and novel typing method presented here will become essential resources for genomic surveillance and epidemiological investigations of this pathogen.

A genome-wide association study identifies a horizontally transferred bacterial surface adhesin gene associated with antimicrobial resistant strains

Carbapenems are a class of last-resort antibiotics; thus, the increase in bacterial carbapenem-resistance is a serious public health threat. Acinetobacter baumannii is one of the microorganisms that can acquire carbapenem-resistance; it causes severe nosocomial infection, and is notoriously difficult to control in hospitals. Recently, a machine-learning approach was first used to analyze the genome sequences of hundreds of susceptible and resistant A. baumannii strains, including those carrying commonly acquired resistant mechanisms, to build a classifier that can predict strain resistance. A complementary approach is to explore novel genetic elements that could be associated with the antimicrobial resistance of strains, independent of known mechanisms. Therefore, we carefully selected A. baumannii strains, spanning various genotypes, from public genome databases, and conducted the first genome-wide association study (GWAS) of carbapenem resistance. We employed a recently developed method, capable of identifying any kind of genetic variation and accounting for bacterial population structure, and evaluated its effectiveness. Our study identified a surface adhesin gene that had been horizontally transferred to an ancestral branch of A. baumannii, as well as a specific region of that gene that appeared to accumulate multiple individual variations across the different branches of carbapenem-resistant A. baumannii strains.

The Establishment and Diversification of Epidemic-Associated Serogroup W Meningococcus in the African Meningitis Belt, 1994 to 2012

Epidemics of invasive meningococcal disease (IMD) caused by meningococcal serogroup A have been eliminated from the sub-Saharan African so-called "meningitis belt" by the meningococcal A conjugate vaccine (MACV), and yet, other serogroups continue to cause epidemics. Neisseria meningitidis serogroup W remains a major cause of disease in the region, with most isolates belonging to clonal complex 11 (CC11). Here, the genetic variation within and between epidemic-associated strains was assessed by sequencing the genomes of 92 N. meningitidis serogroup W isolates collected between 1994 and 2012 from both sporadic and epidemic IMD cases, 85 being from selected meningitis belt countries. The sequenced isolates belonged to either CC175 (n = 9) or CC11 (n = 83). The CC11 N. meningitidis serogroup W isolates belonged to a single lineage comprising four major phylogenetic subclades. Separate CC11 N. meningitidis serogroup W subclades were associated with the 2002 and 2012 Burkina Faso epidemics. The subclade associated with the 2012 epidemic included isolates found in Burkina Faso and Mali during 2011 and 2012, which descended from a strain very similar to the Hajj (Islamic pilgrimage to Mecca)-related Saudi Arabian outbreak strain from 2000. The phylogeny of isolates from 2012 reflected their geographic origin within Burkina Faso, with isolates from the Malian border region being closely related to the isolates from Mali. Evidence of ongoing evolution, international transmission, and strain replacement stresses the importance of maintaining N. meningitidis surveillance in Africa following the MACV implementation. IMPORTANCE Meningococcal disease (meningitis and bloodstream infections) threatens millions of people across the meningitis belt of sub-Saharan Africa. A vaccine introduced in 2010 protects against Africa's then-most common cause of meningococcal disease, N. meningitidis serogroup A. However, other serogroups continue to cause epidemics in the region-including serogroup W. The rapid identification of strains that have been associated with prior outbreaks can improve the assessment of outbreak risk and enable timely preparation of public health responses, including vaccination. Phylogenetic analysis of newly sequenced serogroup W strains isolated from 1994 to 2012 identified two groups of strains linked to large epidemics in Burkina Faso, one being descended from a strain that caused an outbreak during the Hajj pilgrimage in 2000. We find that applying whole-genome sequencing to meningococcal disease surveillance collections improves the discrimination among strains, even within a single nation-wide epidemic, which can be used to better understand pathogen spread.

Listeria monocytogenes sequence type 1 is predominant in ruminant rhombencephalitis

Listeria (L.) monocytogenes is an opportunistic pathogen causing life-threatening infections in diverse mammalian species including humans and ruminants. As little is known on the link between strains and clinicopathological phenotypes, we studied potential strain-associated virulence and organ tropism in L. monocytogenes isolates from well-defined ruminant cases of clinical infections and the farm environment. The phylogeny of isolates and their virulence-associated genes were analyzed by multilocus sequence typing (MLST) and sequence analysis of virulence-associated genes. Additionally, a panel of representative isolates was subjected to in vitro infection assays. Our data suggest the environmental exposure of ruminants to a broad range of strains and yet the strong association of sequence type (ST) 1 from clonal complex (CC) 1 with rhombencephalitis, suggesting increased neurotropism of ST1 in ruminants, which is possibly related to its hypervirulence. This study emphasizes the importance of considering clonal background of L. monocytogenes isolates in surveillance, epidemiological investigation and disease control.

gmos: Rapid Detection of Genome Mosaicism over Short Evolutionary Distances

Prokaryotic and viral genomes are often altered by recombination and horizontal gene transfer. The existing methods for detecting recombination are primarily aimed at viral genomes or sets of loci, since the expensive computation of underlying statistical models often hinders the comparison of complete prokaryotic genomes. As an alternative, alignment-free solutions are more efficient, but cannot map (align) a query to subject genomes. To address this problem, we have developed gmos (Genome MOsaic Structure), a new program that determines the mosaic structure of query genomes when compared to a set of closely related subject genomes. The program first computes local alignments between query and subject genomes and then reconstructs the query mosaic structure by choosing the best local alignment for each query region. To accomplish the analysis quickly, the program mostly relies on pairwise alignments and constructs multiple sequence alignments over short overlapping subject regions only when necessary. This fine-tuned implementation achieves an efficiency comparable to an alignment-free tool. The program performs well for simulated and real data sets of closely related genomes and can be used for fast recombination detection; for instance, when a new prokaryotic pathogen is discovered. As an example, gmos was used to detect genome mosaicism in a pathogenic Enterococcus faecium strain compared to seven closely related genomes. The analysis took less than two minutes on a single 2.1 GHz processor. The output is available in fasta format and can be visualized using an accessory program, gmosDraw (freely available with gmos).

Multilocus Sequence Typing of Strains of Bacterial Spot of Lettuce Collected in the United States

Studies on genetic diversity and recombination in bacterial pathogens are providing a better understanding of the mechanisms shaping bacterial diversity, which can affect disease control. Xanthomonas campestris pv. vitians, causal agent of bacterial leaf spot of lettuce, is a threat to the worldwide lettuce industry. We examined the genetic variation within a sample of 83 strains from California, Florida, and Ohio using multilocus sequence typing of six housekeeping genes, totaling 2.7 kb. Additionally, polymorphism in two virulence-related genes, hrpB2 and a putative glycosyl hydrolase, were examined. Based on housekeeping genes, we found three genetic groups of strains that were all able to induce the disease. These included strains collected from weeds and irrigation water that had haplotypes identical to strains from diseased lettuce. High linkage disequilibrium across the sequenced loci indicates that the pathogen is predominantly clonal but recombination has contributed to the observed sequence variation. Although there was significant genetic variation in X. campestris pv. vitians within and among sampled states, identical haplotypes were observed across all three states. This finding suggests that seedborne inoculum may contribute to the diversity of X. campestris pv. vitians in the United States. Knowledge of the genetic structure of the pathogen may be used for developing resistant lettuce varieties.

Evolutionary study of Yersinia genomes deciphers emergence of human pathogenic species

On record, there are 17 species in the Yersinia genus, of which three are known to be pathogenic to human. While the chromosomal and pYV (or pCD1) plasmid-borne virulence genes as well as pathogenesis of these three species are well studied, their genomic evolution is poorly understood. Our study aims to predict the key evolutionary events that led to the emergence of pathogenic Yersinia species by analyzing gene gain-and-loss, virulence genes, and “Clustered regularly-interspaced short palindromic repeats”. Our results suggest that the most recent ancestor shared by the human pathogenic Yersinia was most probably an environmental species that had adapted to the human body. This might have led to ecological specialization that diverged Yersinia into ecotypes and distinct lineages based on differential gene gain-and-loss in different niches. Our data also suggest that Y. pseudotuberculosis group might be the donor of the ail virulence gene to Y. enterocolitica. Hence, we postulate that evolution of human pathogenic Yersinia might not be totally in parallel, but instead, there were lateral gene transfer events. Furthermore, the presence of virulence genes seems to be important for the positive selection of virulence plasmid. Our studies provide better insights into the evolutionary biology of these bacteria.

Architecture of a Species: Phylogenomics of Staphylococcus aureus

A deluge of whole-genome sequencing has begun to give insights into the patterns and processes of microbial evolution, but genome sequences have accrued in a haphazard manner, with biased sampling of natural variation that is driven largely by medical and epidemiological priorities. For instance, there is a strong bias for sequencing epidemic lineages of methicillin-resistant Staphylococcus aureus (MRSA) over sensitive isolates (methicillin-sensitive S. aureus: MSSA). As more diverse genomes are sequenced the emerging picture is of a highly subdivided species with a handful of relatively clonal groups (complexes) that, at any given moment, dominate in particular geographical regions. The establishment of hegemony of particular clones appears to be a dynamic process of successive waves of replacement of the previously dominant clone. Here we review the phylogenomic structure of a diverse range of S. aureus, including both MRSA and MSSA. We consider the utility of the concept of the 'core' genome and the impact of recombination and horizontal transfer. We argue that whole-genome surveillance of S. aureus populations could lead to better forecasting of antibiotic resistance and virulence of emerging clones, and a better understanding of the elusive biological factors that determine repeated strain replacement.

Putatively novel serotypes and the potential for reduced vaccine effectiveness: capsular locus diversity revealed among 5405 pneumococcal genomes

The pneumococcus is a leading global pathogen and a key virulence factor possessed by the majority of pneumococci is an antigenic polysaccharide capsule ('serotype'), which is encoded by the capsular (cps) locus. Approximately 100 different serotypes are known, but the extent of sequence diversity within the cps loci of individual serotypes is not well understood. Investigating serotype-specific sequence variation is crucial to the design of sequence-based serotyping methodology, understanding pneumococcal conjugate vaccine (PCV) effectiveness and the design of future PCVs. The availability of large genome datasets makes it possible to assess population-level variation among pneumococcal serotypes and in this study 5405 pneumococcal genomes were used to investigate cps locus diversity among 49 different serotypes. Pneumococci had been recovered between 1916 and 2014 from people of all ages living in 51 countries. Serotypes were deduced bioinformatically, cps locus sequences were extracted and variation was assessed within the cps locus, in the context of pneumococcal genetic lineages. Overall, cps locus sequence diversity varied markedly: low to moderate diversity was revealed among serogroups/types 1, 3, 7, 9, 11 and 22; whereas serogroups/types 6, 19, 23, 14, 15, 18, 33 and 35 displayed high diversity. Putative novel and/or hybrid cps loci were identified among all serogroups/types apart from 1, 3 and 9. This study demonstrated that cps locus sequence diversity varied widely between serogroups/types. Investigation of the biochemical structure of the polysaccharide capsule of major variants, particularly PCV-related serotypes and those that appear to be novel or hybrids, is warranted.

Recurrent evolution of host and vector association in bacteria of the Borrelia burgdorferi sensu lato species complex

Background

The Borrelia burgdorferi sensu lato (s.l.) species complex consists of tick-transmitted bacteria and currently comprises approximately 20 named and proposed genospecies some of which are known to cause Lyme Borreliosis. Species have been defined via genetic distances and ecological niches they occupy. Understanding the evolutionary relationship of species of the complex is fundamental to explaining patterns of speciation. This in turn forms a crucial basis to frame testable hypotheses concerning the underlying processes including host and vector adaptations.

Results

Illumina Technology was used to obtain genome-wide sequence data for 93 strains of 14 named genospecies of the B. burgdorferi species complex and genomic data already published for 18 additional strain (including one new species) was added. Phylogenetic reconstruction based on 114 orthologous single copy genes shows that the genospecies represent clearly distinguishable taxa with recent and still ongoing speciation events apparent in Europe and Asia. The position of Borrelia species in the phylogeny is consistent with host associations constituting a major driver for speciation. Interestingly, the data also demonstrate that vector associations are an additional driver for diversification in this tick-borne species complex. This is particularly obvious in B. bavariensis, a rodent adapted species that has diverged from the bird-associated B. garinii most likely in Asia. It now consists of two populations one of which most probably invaded Europe following adaptation to a new vector (Ixodes ricinus) and currently expands its distribution range.

Conclusions

The results imply that genotypes/species with novel properties regarding host or vector associations have evolved recurrently during the history of the species complex and may emerge at any time. We suggest that the finding of vector associations as a driver for diversification may be a general pattern for tick-borne pathogens. The core genome analysis presented here provides an important source for investigations of the underlying mechanisms of speciation in tick-borne pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3016-4) contains supplementary material, which is available to authorized users.

A Primer on Infectious Disease Bacterial Genomics

The number of large-scale genomics projects is increasing due to the availability of affordable high-throughput sequencing (HTS) technologies. The use of HTS for bacterial infectious disease research is attractive because one whole-genome sequencing (WGS) run can replace multiple assays for bacterial typing, molecular epidemiology investigations, and more in-depth pathogenomic studies. The computational resources and bioinformatics expertise required to accommodate and analyze the large amounts of data pose new challenges for researchers embarking on genomics projects for the first time. Here, we present a comprehensive overview of a bacterial genomics projects from beginning to end, with a particular focus on the planning and computational requirements for HTS data, and provide a general understanding of the analytical concepts to develop a workflow that will meet the objectives and goals of HTS projects.

Emergence of endemic MLST non-typeable vancomycin-resistant Enterococcus faecium

BACKGROUND

Enterococcus faecium is a major nosocomial pathogen causing significant morbidity and mortality worldwide. Assessment of E. faecium using MLST to understand the spread of this organism is an important component of hospital infection control measures. Recent studies, however, suggest that MLST might be inadequate for E. faecium surveillance.

OBJECTIVES

To use WGS to characterize recently identified vancomycin-resistant E. faecium (VREfm) isolates non-typeable by MLST that appear to be causing a multi-jurisdictional outbreak in Australia.

METHODS

Illumina NextSeq and Pacific Biosciences SMRT sequencing platforms were used to determine the genome sequences of 66 non-typeable E. faecium (NTEfm) isolates. Phylogenetic and bioinformatics analyses were subsequently performed using a number of in silico tools.

RESULTS

Sixty-six E. faecium isolates were identified by WGS from multiple health jurisdictions in Australia that could not be typed by MLST due to a missing pstS allele. SMRT sequencing and complete genome assembly revealed a large chromosomal rearrangement in representative strain DMG1500801, which likely facilitated the deletion of the pstS region. Phylogenomic analysis of this population suggests that deletion of pstS within E. faecium has arisen independently on at least three occasions. Importantly, the majority of these isolates displayed a vancomycin-resistant genotype.

CONCLUSIONS

We have identified NTEfm isolates that appear to be causing a multi-jurisdictional outbreak in Australia. Identification of these isolates has important implications for MLST-based typing activities designed to monitor the spread of VREfm and provides further evidence supporting the use of WGS for hospital surveillance of E. faecium.

Distribution of Bexsero® Antigen Sequence Types (BASTs) in invasive meningococcal disease isolates: Implications for immunisation

Serogroup B is the only major disease-associated capsular group of Neisseria meningitidis for which no protein-polysaccharide conjugate vaccine is available. This has led to the development of multi-component protein-based vaccines that target serogroup B invasive meningococcal disease (IMD), including Bexsero®, which was implemented for UK infants in 2015, and Trumenba®. Given the diversity of meningococcal protein antigens, post-implementation surveillance of IMD isolates, including characterisation of vaccine antigens, is essential for assessing the effectiveness of such vaccines. Whole genome sequencing (WGS), as realised in the Meningitis Research Foundation Meningococcus Genome Library (MRF-MGL), provides a rapid, comprehensive, and cost-effective approach to this. To facilitate the surveillance of the antigen targets included in Bexsero® (fHbp, PorA, NHBA and NadA) for protective immunity, a Bexsero® Antigen Sequence Type (BAST) scheme, based on deduced peptide sequence variants, was implemented in the PubMLST.org/neisseria database, which includes the MRF-MGL and other isolate collections. This scheme enabled the characterisation of vaccine antigen variants and here the invasive meningococci isolated in Great Britain and Ireland in the epidemiological years 2010/11 to 2013/14 are analysed. Many unique BASTs (647) were present, but nine of these accounted for 39% (775/1966) of isolates, with some temporal and geographic differences in BAST distribution. BASTs were strongly associated with other characteristics, such as serogroup and clonal complex (cc), and a significant increase in BAST-2 was associated with increased prevalence of serogroup W clonal complex 11 meningococci. Potential coverage was assessed by the examination of the antigen peptide sequences present in the vaccine and epidemiological dataset. There were 22.8-30.8% exact peptide matches to Bexsero® components and predicted coverage of 66.1%, based on genotype-phenotype modelling for 63.7% of serogroup B isolates from 2010/14 in UK and Ireland. While there are many caveats to this estimate, it lies within the range of other published estimates.

Genomics of Natural Populations of Staphylococcus aureus

Staphylococcus aureus is a major human pathogen and an important cause of livestock infections. The first S. aureus genomes to be published, 15 years ago, provided the first view of genome structure and gene content. Since then, thousands of genomes from a wide array of strains from different sources have been sequenced. Comparison of these sequences has resulted in broad insights into population structure, bacterial evolution, clone emergence and expansion, and the molecular basis of niche adaptation. Furthermore, this information is now being applied clinically in outbreak investigations to inform infection control measures and to determine appropriate treatment regimens. In this review, we summarize some of the broad insights into S. aureus biology gained from the analysis of genomes and discuss future directions and opportunities in this dynamic field of research.

Whole-genome sequencing to determine transmission of Neisseria gonorrhoeae: an observational study

BACKGROUND

New approaches are urgently required to address increasing rates of gonorrhoea and the emergence and global spread of antibiotic-resistant Neisseria gonorrhoeae. We used whole-genome sequencing to study transmission and track resistance in N gonorrhoeae isolates.

METHODS

We did whole-genome sequencing of isolates obtained from samples collected from patients attending sexual health services in Brighton, UK, between Jan 1, 2011, and March 9, 2015. We also included isolates from other UK locations, historical isolates from Brighton, and previous data from a US study. Samples from symptomatic patients and asymptomatic sexual health screening underwent nucleic acid amplification testing; positive samples and all samples from symptomatic patients were cultured for N gonorrhoeae, and resulting isolates were whole-genome sequenced. Cefixime susceptibility testing was done in selected isolates by agar incorporation, and we used sequence data to determine multi-antigen sequence types and penA genotypes. We derived a transmission nomogram to determine the plausibility of direct or indirect transmission between any two cases depending on the time between samples: estimated mutation rates, plus diversity noted within patients across anatomical sites and probable transmission pairs, were used to fit a coalescent model to determine the number of single nucleotide polymorphisms expected.

FINDINGS

1407 (98%) of 1437 Brighton isolates between Jan 1, 2011, and March 9, 2015 were successfully sequenced. We identified 1061 infections from 907 patients. 281 (26%) of these infections were indistinguishable (ie, differed by zero single nucleotide polymorphisms) from one or more previous cases, and 786 (74%) had evidence of a sampled direct or indirect Brighton source. We observed multiple related samples across geographical locations. Of 1273 infections in Brighton (including historical data), 225 (18%) were linked to another case elsewhere in the UK, and 115 (9%) to a case in the USA. Four lineages initially identified in Brighton could be linked to 70 USA sequences, including 61 from a lineage carrying the mosaic penA XXXIV allele, which is associated with reduced cefixime susceptibility.

INTERPRETATION

We present a whole-genome-sequencing-based tool for genomic contact tracing of N gonorrhoeae and demonstrate local, national, and international transmission. Whole-genome sequencing can be applied across geographical boundaries to investigate gonorrhoea transmission and to track antimicrobial resistance.

FUNDING

Oxford National Institute for Health Research Health Protection Research Unit and Biomedical Research Centre.

Population genomics of the symbiotic plasmids of sympatric nitrogen-fixing Rhizobium species associated with Phaseolus vulgaris

Cultivated common beans are the primary protein source for millions of people around the world who subsist on low-input agriculture, enabled by the symbiotic N2 -fixation these legumes perform in association with rhizobia. Within a single agricultural plot, multiple Rhizobium species can nodulate bean roots, but it is unclear how genetically isolated these species remain in sympatry. To better understand this issue, we sequenced and compared the genomes of 33 strains isolated from the rhizosphere and root nodules of a particular bean variety grown in the same agricultural plot. We found that the Rhizobium species we observed coexist with low genetic recombination across their core genomes. Accessory plasmids thought to be necessary for the saprophytic lifestyle in soil show similar levels of genetic isolation, but with higher rates of recombination than the chromosomes. However, the symbiotic plasmids are extremely similar, with high rates of recombination and do not appear to have co-evolved with the chromosome or accessory plasmids. Therefore, while Rhizobium species are genetically isolated units within the microbial community, a common symbiotic plasmid allows all Rhizobium species to engage in symbiosis with the same host in a single agricultural plot.

Comparative Genomics Analysis of Streptomyces Species Reveals Their Adaptation to the Marine Environment and Their Diversity at the Genomic Level

Over 200 genomes of streptomycete strains that were isolated from various environments are available from the NCBI. However, little is known about the characteristics that are linked to marine adaptation in marine-derived streptomycetes. The particularity and complexity of the marine environment suggest that marine streptomycetes are genetically diverse. Here, we sequenced nine strains from the Streptomyces genus that were isolated from different longitudes, latitudes, and depths of the South China Sea. Then we compared these strains to 22 NCBI downloaded streptomycete strains. Thirty-one streptomycete strains are clearly grouped into a marine-derived subgroup and multiple source subgroup-based phylogenetic tree. The phylogenetic analyses have revealed the dynamic process underlying streptomycete genome evolution, and lateral gene transfer is an important driving force during the process. Pan-genomics analyses have revealed that streptomycetes have an open pan-genome, which reflects the diversity of these streptomycetes and guarantees the species a quick and economical response to diverse environments. Functional and comparative genomics analyses indicate that the marine-derived streptomycetes subgroup possesses some common characteristics of marine adaptation. Our findings have expanded our knowledge of how ocean isolates of streptomycete strains adapt to marine environments. The availability of streptomycete genomes from the South China Sea will be beneficial for further analysis on marine streptomycetes and will enrich the South China Sea's genetic data sources.

Maternal colonization with Streptococcus agalactiae and associated stillbirth and neonatal disease in coastal Kenya

Streptococcus agalactiae (group B streptococcus, GBS) causes neonatal disease and stillbirth, but its burden in sub-Saharan Africa is uncertain. We assessed maternal recto-vaginal GBS colonization (7,967 women), stillbirth and neonatal disease. Whole-genome sequencing was used to determine serotypes, sequence types and phylogeny. We found low maternal GBS colonization prevalence (934/7,967, 12%), but comparatively high incidence of GBS-associated stillbirth and early onset neonatal disease (EOD) in hospital (0.91 (0.25-2.3)/1,000 births and 0.76 (0.25-1.77)/1,000 live births, respectively). However, using a population denominator, EOD incidence was considerably reduced (0.13 (0.07-0.21)/1,000 live births). Treated cases of EOD had very high case fatality (17/36, 47%), especially within 24 h of birth, making under-ascertainment of community-born cases highly likely, both here and in similar facility-based studies. Maternal GBS colonization was less common in women with low socio-economic status, HIV infection and undernutrition, but when GBS-colonized, they were more probably colonized by the most virulent clone, CC17. CC17 accounted for 267/915 (29%) of maternal colonizing (265/267 (99%) serotype III; 2/267 (0.7%) serotype IV) and 51/73 (70%) of neonatal disease cases (all serotype III). Trivalent (Ia/II/III) and pentavalent (Ia/Ib/II/III/V) vaccines would cover 71/73 (97%) and 72/73 (99%) of disease-causing serotypes, respectively. Serotype IV should be considered for inclusion, with evidence of capsular switching in CC17 strains.

Regulation of genetic flux between bacteria by restriction-modification systems

Restriction-modification (R-M) systems are often regarded as bacteria's innate immune systems, protecting cells from infection by mobile genetic elements (MGEs). Their diversification has been recently associated with the emergence of particularly virulent lineages. However, we have previously found more R-M systems in genomes carrying more MGEs. Furthermore, it has been suggested that R-M systems might favor genetic transfer by producing recombinogenic double-stranded DNA ends. To test whether R-M systems favor or disfavor genetic exchanges, we analyzed their frequency with respect to the inferred events of homologous recombination and horizontal gene transfer within 79 bacterial species. Genetic exchanges were more frequent in bacteria with larger genomes and in those encoding more R-M systems. We created a recognition target motif predictor for Type II R-M systems that identifies genomes encoding systems with similar restriction sites. We found more genetic exchanges between these genomes, independently of their evolutionary distance. Our results reconcile previous studies by showing that R-M systems are more abundant in promiscuous species, wherein they establish preferential paths of genetic exchange within and between lineages with cognate R-M systems. Because the repertoire and/or specificity of R-M systems in bacterial lineages vary quickly, the preferential fluxes of genetic transfer within species are expected to constantly change, producing time-dependent networks of gene transfer.

Inferences from tip-calibrated phylogenies: a review and a practical guide

Molecular dating of phylogenetic trees is a growing discipline using sequence data to co‐estimate the timing of evolutionary events and rates of molecular evolution. All molecular‐dating methods require converting genetic divergence between sequences into absolute time. Historically, this could only be achieved by associating externally derived dates obtained from fossil or biogeographical evidence to internal nodes of the tree. In some cases, notably for fast‐evolving genomes such as viruses and some bacteria, the time span over which samples were collected may cover a significant proportion of the time since they last shared a common ancestor. This situation allows phylogenetic trees to be calibrated by associating sampling dates directly to the sequences representing the tips (terminal nodes) of the tree. The increasing availability of genomic data from ancient DNA extends the applicability of such tip‐based calibration to a variety of taxa including humans, extinct megafauna and various microorganisms which typically have a scarce fossil record. The development of statistical models accounting for heterogeneity in different aspects of the evolutionary process while accommodating very large data sets (e.g. whole genomes) has allowed using tip‐dating methods to reach inferences on divergence times, substitution rates, past demography or the age of specific mutations on a variety of spatiotemporal scales. In this review, we summarize the current state of the art of tip dating, discuss some recent applications, highlight common pitfalls and provide a ‘how to’ guide to thoroughly perform such analyses.

Identifying lineage effects when controlling for population structure improves power in bacterial association studies

Bacteria pose unique challenges for genome-wide association studies because of strong structuring into distinct strains and substantial linkage disequilibrium across the genome(1,2). Although methods developed for human studies can correct for strain structure(3,4), this risks considerable loss-of-power because genetic differences between strains often contribute substantial phenotypic variability(5). Here, we propose a new method that captures lineage-level associations even when locus-specific associations cannot be fine-mapped. We demonstrate its ability to detect genes and genetic variants underlying resistance to 17 antimicrobials in 3,144 isolates from four taxonomically diverse clonal and recombining bacteria: Mycobacterium tuberculosis, Staphylococcus aureus, Escherichia coli and Klebsiella pneumoniae. Strong selection, recombination and penetrance confer high power to recover known antimicrobial resistance mechanisms and reveal a candidate association between the outer membrane porin nmpC and cefazolin resistance in E. coli. Hence, our method pinpoints locus-specific effects where possible and boosts power by detecting lineage-level differences when fine-mapping is intractable.

Detection of homologous recombination in closely related strains

Detection of recombination events in a bacterial genome is both important from the evolutionary point of view, and of practical interest. Indeed, homologous recombination (HR) plays a major role in the exchange of antigenic determinants between strains. There exist statistical methods to detect recently recombined segments in whole-genome sequences that use a high local density of substitutions as a signal of HR events with a source outside considered strains. However, it is difficult to detect the HR events within a set of strains, which represent whole species diversity, due to a low number of substitutions in recombined segments and high level of diversity of strains. Here, we analyzed HR in 20 Escherichia coli (E. coli) strains to define what fraction of segments with a high substitution rate were introduced in a genome by HR. For detection of HR, we used the segmentation, performed by the adaptive weights smoothing (AWS) algorithm. It detects sharp changes in the structure of observed data analyzing only qualitative structural information. We validated the approach on simulated data, applied it to the analysis of E. coli strains, and determined the recombination rates between phylogroups.

Genomic Analysis of Salmonella enterica Serovar Typhimurium Characterizes Strain Diversity for Recent U.S. Salmonellosis Cases and Identifies Mutations Linked to Loss of Fitness under Nitrosative and Oxidative Stress

Salmonella enterica serovar Typhimurium is one of the most common S. enterica serovars associated with U.S. foodborne outbreaks. S. Typhimurium bacteria isolated from humans exhibit wide-ranging virulence phenotypes in inbred mice, leading to speculation that some strains are more virulent in nature. However, it is unclear whether increased virulence in humans is related to organism characteristics or initial treatment failure due to antibiotic resistance. Strain diversity and genetic factors contributing to differential human pathogenicity remain poorly understood. We reconstructed phylogeny, resolved genetic population structure, determined gene content and nucleotide variants, and conducted targeted phenotyping assays for S. Typhimurium strains collected between 1946 and 2012 from humans and animals in the United States and abroad. Strains from recent U.S. salmonellosis cases were associated with five S. Typhimurium lineages distributed within three phylogenetic clades, which are not restricted by geography, year of acquisition, or host. Notably, two U.S. strains and four Mexican strains are more closely related to strains associated with human immunodeficiency virus (HIV)-infected individuals in sub-Saharan Africa than to other North American strains. Phenotyping studies linked variants specific to these strains in hmpA and katE to loss of fitness under nitrosative and oxidative stress, respectively. These results suggest that U.S. salmonellosis is caused by diverse S. Typhimurium strains circulating worldwide. One lineage has mutations in genes affecting fitness related to innate immune system strategies for fighting pathogens and may be adapting to immunocompromised humans by a reduction in virulence capability, possibly due to a lack of selection for its maintenance as a result of the worldwide HIV epidemic.

Nontyphoidal Salmonella bacteria cause an estimated 1.2 million illnesses annually in the United States, 80 million globally, due to ingestion of contaminated food or water. Salmonella Typhimurium is one of the most common serovars associated with foodborne illness, causing self-limiting gastroenteritis and, in approximately 5% of infected patients, systemic infection. Although some S. Typhimurium strains are speculated to be more virulent than others, it is unknown how strain diversity and genetic factors contribute to differential human pathogenicity. Ours is the first study to examine the diversity of S. Typhimurium associated with recent cases of U.S. salmonellosis and to provide some initial correlation between observed genotypes and phenotypes. Definition of specific S. Typhimurium lineages based on such phenotype/genotype correlations may identify strains with greater capability of associating with specific food sources, allowing outbreaks to be more quickly identified. Additionally, defining simple correlates of pathogenesis may have predictive value for patient outcome.

Repeated local emergence of carbapenem-resistant Acinetobacter baumannii in a single hospital ward

We recently reported a dramatic increase in the prevalence of carbapenem-resistant Acinetobacter baumannii infections in the intensive care unit (ICU) of a Vietnamese hospital. This upsurge was associated with a specific oxa23-positive clone that was identified by multilocus VNTR analysis. Here, we used whole-genome sequence analysis to dissect the emergence of carbapenem-resistant A. baumannii causing ventilator-associated pneumonia (VAP) in the ICU during 2009–2012. To provide historical context and distinguish microevolution from strain introduction, we compared these genomes with those of A. baumannii asymptomatic carriage and VAP isolates from this same ICU collected during 2003–2007. We identified diverse lineages co-circulating over many years. Carbapenem resistance was associated with the presence of oxa23, oxa40, oxa58 and ndm1 genes in multiple lineages. The majority of resistant isolates were oxa23-positive global clone GC2; fine-scale phylogenomic analysis revealed five distinct GC2 sublineages within the ICU that had evolved locally via independent chromosomal insertions of oxa23 transposons. The increase in infections caused by carbapenem-resistant A. baumannii was associated with transposon-mediated transmission of a carbapenemase gene, rather than clonal expansion or spread of a carbapenemase-harbouring plasmid. Additionally, we found evidence of homologous recombination creating diversity within the local GC2 population, including several events resulting in replacement of the capsule locus. We identified likely donors of the imported capsule locus sequences amongst the A. baumannii isolated on the same ward, suggesting that diversification was largely facilitated via reassortment and sharing of genetic material within the localized A. baumannii population.

Refined analyses suggest that recombination is a minor source of genomic diversity in Pseudomonas aeruginosa chronic cystic fibrosis infections

Chronic bacterial airway infections in people with cystic fibrosis (CF) are often caused by Pseudomonas aeruginosa, typically showing high phenotypic diversity amongst co-isolates from the same sputum sample. Whilst adaptive evolution during chronic infections has been reported, the genetic mechanisms underlying the observed rapid within-population diversification are not well understood. Two recent conflicting reports described very high and low rates of homologous recombination in two closely related P. aeruginosa populations from the lungs of different chronically infected CF patients. To investigate the underlying cause of these contrasting observations, we combined the short read datasets from both studies and applied a new comparative analysis. We inferred low rates of recombination in both populations. The discrepancy in the findings of the two previous studies can be explained by differences in the application of variant calling techniques. Two novel algorithms were developed that filter false-positive variants. The first algorithm filters variants on the basis of ambiguity within duplications in the reference genome. The second omits probable false-positive variants at regions of non-homology between reference and sample caused by structural rearrangements. As gains and losses of prophage or genomic islands are frequent causes of chromosomal rearrangements within microbial populations, this filter has broad appeal for mitigating false-positive variant calls. Both algorithms are available in a Python package.

Utility of Whole-Genome Sequencing in Characterizing Acinetobacter Epidemiology and Analyzing Hospital Outbreaks

Acinetobacter baumannii frequently causes nosocomial infections and outbreaks. Whole-genome sequencing (WGS) is a promising technique for strain typing and outbreak investigations. We compared the performance of conventional methods with WGS for strain typing clinical Acinetobacter isolates and analyzing a carbapenem-resistant A. baumannii (CRAB) outbreak. We performed two band-based typing techniques (pulsed-field gel electrophoresis and repetitive extragenic palindromic-PCR), multilocus sequence type (MLST) analysis, and WGS on 148 Acinetobacter calcoaceticus-A. baumannii complex bloodstream isolates collected from a single hospital from 2005 to 2012. Phylogenetic trees inferred from core-genome single nucleotide polymorphisms (SNPs) confirmed three Acinetobacter species within this collection. Four major A. baumannii clonal lineages (as defined by MLST) circulated during the study, three of which are globally distributed and one of which is novel. WGS indicated that a threshold of 2,500 core SNPs accurately distinguished A. baumannii isolates from different clonal lineages. The band-based techniques performed poorly in assigning isolates to clonal lineages and exhibited little agreement with sequence-based techniques. After applying WGS to a CRAB outbreak that occurred during the study, we identified a threshold of 2.5 core SNPs that distinguished nonoutbreak from outbreak strains. WGS was more discriminatory than the band-based techniques and was used to construct a more accurate transmission map that resolved many of the plausible transmission routes suggested by epidemiologic links. Our study demonstrates that WGS is superior to conventional techniques for A. baumannii strain typing and outbreak analysis. These findings support the incorporation of WGS into health care infection prevention efforts.

Practical Approaches for Detecting Selection in Microbial Genomes

Microbial genome evolution is shaped by a variety of selective pressures. Understanding how these processes occur can help to address important problems in microbiology by explaining observed differences in phenotypes, including virulence and resistance to antibiotics. Greater access to whole-genome sequencing provides microbiologists with the opportunity to perform large-scale analyses of selection in novel settings, such as within individual hosts. This tutorial aims to guide researchers through the fundamentals underpinning popular methods for measuring selection in pathogens. These methods are transferable to a wide variety of organisms, and the exercises provided are designed for researchers with any level of programming experience.

Within-host evolution of bacterial pathogens

Whole-genome sequencing has opened the way for investigating the dynamics and genomic evolution of bacterial pathogens during the colonization and infection of humans. The application of this technology to the longitudinal study of adaptation in an infected host--in particular, the evolution of drug resistance and host adaptation in patients who are chronically infected with opportunistic pathogens--has revealed remarkable patterns of convergent evolution, suggestive of an inherent repeatability of evolution. In this Review, we describe how these studies have advanced our understanding of the mechanisms and principles of within-host genome evolution, and we consider the consequences of findings such as a potent adaptive potential for pathogenicity. Finally, we discuss the possibility that genomics may be used in the future to predict the clinical progression of bacterial infections and to suggest the best option for treatment.