Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 23, Number 5—May 2017
Research

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

Rodrigo Bacigalupe, Diane Lindsay, Giles Edwards, and J. Ross FitzgeraldComments to Author 
Author affiliations: The Roslin Institute, University of Edinburgh, Midlothian, Scotland, UK (R. Bacigalupe, J.R. Fitzgerald); Glasgow Royal Infirmary, Glasgow, Scotland, UK (D. Lindsay, G. Edwards)

Cite This Article

Abstract

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.

Legionellosis presents as 2 clinically distinct forms: an influenza-like illness called Pontiac fever and a severe pneumonia known as Legionnaires’ disease (1). In Europe and the United States, most legionellosis cases are caused by Legionella pneumophila serogroup 1 (1,2); <5% of cases are caused by nonpneumophila Legionella spp. (3,4). In Australasia, New Zealand, and some countries in Asia, infections caused by L. longbeachae occur at comparable levels to infections caused by L. pneumophila (57). Unlike L. pneumophila infections, which are typically linked to artificial water systems, L. longbeachae infections are associated with exposure to soil, compost, and potting mixes (8).

The number of legionellosis cases caused by L. longbeachae is increasing worldwide (7), with a notable rise reported across Europe (911). Within the United Kingdom, most L. longbeachae infections have been identified in Scotland, where 6 cases were diagnosed during 2008–2012 (12) and another 6 were diagnosed in the summer of 2013 and represented a singular increased incidence or cluster with all patients requiring intensive care hospitalization (11). Epidemiologic investigation revealed that most patients from the 2013 cluster were avid gardeners, and L. longbeachae was isolated from respiratory secretions and from samples of the growing media they had used for gardening before becoming ill (11,12). However, an investigation into the provenance of the growing media did not reveal a single commercial or manufacturing source that would suggest a common origin for the L. longbeachae associated with the outbreak (11).

Molecular typing methods used to discriminate between L. longbeachae and other Legionella spp. and between the 2 L. longbeachae serogroups have limited efficacy, and although considerable evidence supports growing media as a source for L. longbeachae infections (13,14), there is still a lack of genetic evidence for an epidemiologic link. Furthermore, a population genomic study involving large numbers of L. pneumophila isolates has been conducted (15,16), but the same has not been done for L. longbeachae, so the diversity of environmental and pathogenic genotypes and the relationship between them remains unknown for L. longbeachae. To examine the etiology of the 2013 cluster of legionellosis cases in Scotland in the context of L. longbeachae species diversity, we analyzed the genomes of 70 Legionella spp. isolates from 4 countries over 18 years.

Materials and Methods

Bacterial Isolates

We sequenced 65 isolates that had previously been identified as L. longbeachae. These isolates were obtained during 1996–2014 from several patients, growing media samples (including compost and soil), and a hot water supply. Of these isolates, 55 were from Scotland (29 from the 2013 cluster of infections and 26 from other clinical and environmental samples) and 10 were from patients and environmental compost samples in New Zealand (Technical Appendix Table).

In our analysis, we also included all publicly available genome sequences for L. longbeachae: L. longbeachae NSW150 (serogroup 1) and L. longbeachae C-4E7 (serogroup 2) isolated from patients in Australia; and L. longbeachae D-4968 (serogroup 1), L. longbeachae ATCC39642 (serogroup 1), and L. longbeachae 98072 (serogroup 2) isolated from patients in the United States (1719). We sequenced multiple isolates (n = 2 to 5) for each of 3 patients and their linked growing media samples from the 2013 outbreak in Scotland and for 2 additional compost samples. The species of all isolates had been determined by serotyping or macrophage infectivity potentiator (mip) gene sequencing (20,21).

Bacterial Culture, Genomic DNA Isolation, and Whole-Genome Sequencing

We cultured Legionella spp. isolates in a microaerophilic and humid environment at 37°C on BCYE (buffered charcoal yeast extract) agar plates for 48 h. We then picked individual colonies from the plates and grew them in ACES-buffered yeast extract broth containing Legionella BCYE Growth Supplement (Oxoid Ltd., Basingstoke, UK) with shaking at 37°C for 24–48 h. We extracted genomic DNA from fresh cultures by using the QIAGEN DNeasy Blood and Tissue Kit (QIAGEN Benelux B.V., Venlo, the Netherlands).

We prepared sequencing libraries by using the Nextera XT kit for MiSeq or HiSeq (all from Illumina, San Diego, CA, USA) sequencing at Edinburgh Genomics, University of Edinburgh (Edinburgh, Scotland, UK). For each isolate, one 2 × 250–bp or two 2 × 200–bp paired-end sequencing runs were carried out using the MiSeq and HiSeq technologies, respectively. Raw reads were quality checked using FastQC v0.10.1 (22), and primers were trimmed by using Cutadapt (23). We used wgsim software (24) to simulate sequence reads for publicly available, complete whole-genome sequences.

Bioinformatic Analysis and Data Deposition

A detailed description of the bioinformatic analyses is available in the online Technical Appendix. The sequence data for the 65 genomes of Legionella spp. sequenced in this study were deposited in the SRA database (accession no. PRJEB14754).

Results

Limitations of Current Typing Approaches for Legionella spp. Identification

Figure 1

Thumbnail of 16S rRNA gene–based phylogenetic tree of the sequenced genomes and all the cultured and type Legionella spp. strains available in the ribosomal database project (http://rdp.cme.msu.edu/), as accessed in May 2015. Scale bar indicates the mean number of nucleotide substitutions per site.

Figure 1. 6S rRNA gene–based phylogenetic tree of the sequenced genomes and all the cultured and type Legionella spp. strains available in the ribosomal database project (http://rdp.cme.msu.edu/), as accessed in May 2015....

We sequenced 65 isolates obtained from several patients and environmental samples over 18 years in different countries and previously identified as L. longbeachae. To confirm the species identity of the Legionella isolates, we constructed a phylogenetic tree that included all Legionella type strains for which cultures are available, based on the 16S rRNA gene sequence (25). We also built phylogenetic trees based on the whole-genome content and core-genome diversity. For each approach, 64 of the 70 isolates examined co-segregated within the L. longbeachae–specific clade, 4 isolates clustered with Legionella anisa, and 2 belonged to a separate clade that was distinct from all known Legionella spp. (Figure 1; Technical Appendix Figures 1, 2). The species identities were further supported by determination of the average nucleotide identity values (Technical Appendix Figure 3), a widely used method for bacterial species delineation based on genomic relatedness (26). Of note, L. anisa is the most common nonpneumophila Legionella spp. in Europe (2729). In addition, L. longbeachae isolates 13.8642 (from a compost sample from Scotland) and 13.8295 (from a patient in New Zealand) belong to a putative novel Legionella spp. Overall, the data indicate that current serotyping methods and mip gene sequencing are limited in their capacity to identify L. longbeachae to the species level.

To investigate the genetic relatedness of L. longbeachae strains associated with the 2013 outbreak to temporally and geographically distinct isolates, we constructed a core genome–based neighbor-joining tree of the 64 confirmed L. longbeachae isolates obtained from 4 countries over 18 years (Technical Appendix Figure 4). This phylogenetic tree presents a comet-like pattern, with 2 distinct clades separated by 9,911 single-nucleotide polymorphisms, representing the major serogroups (serogroups 1 and 2) previously identified for L. longbeachae (20), each containing isolates from patient and environmental samples from different years. In contrast with findings from a previous analysis of 2 isolates of L. longbeachae serogroup 1 (20), we observed a higher diversity among the 56 isolates within serogroup 1 (Technical Appendix Figures 1, 4); this finding is not unexpected, given the difference in the number of genomes examined. Nevertheless, compared with isolates from the same serogroup in other Legionella spp., such as L. pneumophila serogroup 1 (2% polymorphism) (20), L. longbeachae serogroup 1 exhibits a lower diversity (<0.1% polymorphism). Although serogroup 1 and 2 clades contained isolates from Scotland, Australasia, and the United States, 96% of the isolates from Scotland (including all of the 2013 outbreak isolates) belonged to serogroup 1, suggesting that serogroup 1 may be more clinically relevant in Scotland than in some other countries where L. longbeachae is a more established cause of legionellosis. However, analysis of more isolates from different countries would be required to investigate this observation further.

Effect of Recombination on L. longbeachae Serogroup 1 Population Structure

Figure 2

Thumbnail of Core genome–based maximum-likelihood phylogeny of Legionella longbeachae serogroup 1 isolates corrected for recombination; source, country, year of isolation, relatedness and plasmid carriage are indicated. Related isolates are shown in the same color; those from the 2013 outbreak are indicated by gray. Isolates from the same patient are clustered together but do not co-segregate with cognate compost samples. Scale bar indicates the mean number of nucleotide substitutions per site.

Figure 2. Core genome–based maximum-likelihood phylogeny of Legionella longbeachae serogroup 1 isolates corrected for recombination; source, country, year of isolation, relatedness and plasmid carriage are indicated. Related isolates are shown in the same...

It is established that recombination has played a key role in shaping the evolutionary history of L. pneumophila, but its effect on L. longbeachae population structure is unknown (22,30). This knowledge is critical because for highly recombinant bacteria, recombination networks may represent evolutionary relationships more explicitly than traditional phylogenetic trees. Therefore, we constructed a recombination network of all serogroup 1 isolates by using the neighbor-net algorithm of SplitsTree4 (31). The resultant network displayed a reticulate topology with an extensive reticulated background from which clusters of isolates emerge, supporting an evolutionary history involving recombination (p< 0.01 by ϕ test) (32), followed by clonal expansion and subsequent additional recombination events among some lineages (Technical Appendix Figure 5). Using BratNextGen (33), we identified a total of 94 predicted recombination events affecting more than half of the core genome (1.74 Mb of 3.36 Mb) and representing recent and ancient recombination events of different sizes (range 1,350 bp–350 Kbp) distributed across the phylogeny (Technical Appendix Figure 6). Given the reported limitation in sensitivity of BratNextGen for the identification of all recombination events (34), we also used ClonalFrameML (35), an algorithm that uses maximum likelihood inference to simultaneously detect recombination in bacterial genomes and account for it in phylogenetic reconstruction. The estimated average length of the recombined fragments was 8,047 bp, and the ratio of recombination to mutation was 1.42, indicating a greater role for recombination over mutation in the diversification of L. longbeachae. This estimate is in accordance with early estimates for L. pneumophila based on multiple gene sequence data (36), but it is low compared with recent estimates based on whole-genome sequence data [recombination to mutation ratios of 16.8 (30) or 47.93 (37)]. Differences in the clonal diversity of Legionella spp. sequence datasets used to determine recombination rates could affect the estimates. Reconstruction of the phylogeny after removal of all predicted recombinant sequences resulted in a tree with largely similar clusters of isolates but with reduced branch lengths and variation in the position of nodes deep in the phylogeny (Figure 2).

Accessory Genome Analysis Indicates Extensive Interspecies and Intraspecies Gene Flow

Figure 3

Thumbnail of Legionella longbeachae plasmid analysis: contigs networks reconstructions for 6 representative L. longbeachae types of plasmid content. The networks of the contigs representing the main chromosome and plasmids comprising the genome obtained by using PLACNET (38), a program enabling reconstruction of plasmids from whole-genome sequence datasets. The sizes of the contig nodes (in gray) are proportional to their lengths; continuous lines correspond to scaffold links. Dashed lines repre

Figure 3. Legionella longbeachae plasmid analysis: contigs networks reconstructions for 6 representative L. longbeachae types of plasmid content. The networks of the contigs representing the main chromosome and plasmids comprising the genome obtained...

The extent to which horizontal gene transfer occurs among L. longbeachae isolates and between L. longbeachae and other Legionella spp. is unknown. In our study, the pangenome of L. longbeachae represented by the 56 serogroup 1 isolates was 6,890 genes, including a core genome of 2,574 genes; the average gene content was 3,558 genes per strain. The accessory genome, which included only strain-dependent genes varied from 809 to 1,155 genes, depending on the strain. A parsimony clustering analysis based on the presence or absence of all genes classified the isolates in a manner distinct from that in a core genome–based maximum-likelihood tree, suggesting extensive horizontal gene transfer among L. longbeachae isolates (Technical Appendix Figures 1, 2). BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) analysis of all assembled contigs was used to filter for plasmid-related homologous sequences, revealed 2 major plasmids: pLLO, described previously in L. longbeachae NSW150 (20), and pLELO, originally identified in L. pneumophila subsp. pneumophila (22). Of the 55 serogroup 1 isolates, 36 contained sequences for the pLLO and pLELO plasmids. Of note, the distribution of these plasmids among the L. longbeachae isolates correlated with the gene content–based clustering, whereas the distribution of plasmids in the core genome–based tree was independent of the phylogeny (Figure 2). In addition, 11 isolates appeared to contain plasmids with sequences homologous to those for pLLO and pLELO, which is indicative of recombinant forms of the plasmid. Further examination of plasmid diversity using a modified version of PLACNET (38), a program enabling reconstruction of plasmids from whole-genome sequence datasets, confirmed that some plasmids consisted of a mosaic of recombinant fragments homologous to pLELO, pLLO, or other unknown plasmids (Figure 3). Overall, these data indicate the high prevalence of specific plasmids among L. longbeachae isolates and reveal extensive recombination and horizontal gene transfer among different Legionella spp (39). The high prevalence of plasmids in L. longbeachae is notable, considering these elements may be less common in L. pneumophila (30).

Figure 4

Thumbnail of Variation in gene content between environmental and patient Legionella longbeachae samples. A) Increase in pangenome size with every addition of a L. longbeachae genome. Environmental isolates pangenomes (green) are larger and continue increasing after the addition of 5 genomes, consistent with an open pangenome, but the within-patient pangenome plateaus quickly, consistent with a more closed pangenome. B) Average gene content of environmental isolates is significantly higher than t

Figure 4. Variation in gene content between environmental and patient Legionella longbeachae samples. A) Increase in pangenome size with every addition of a L. longbeachae genome. Environmental isolates pangenomes (green) are larger and...

To examine the possibility that clinical and environmental isolates of L. longbeachae contained genomic differences reflecting their distinct origins, we compared their accessory genome content. For isolates obtained from a single patient sample, the accessory genome was highly conserved compared with those for environmental isolates from a single compost sample or closely related environmental isolates from distinct compost samples (Figure 4, panel A). In addition, considering the average gene content of all sequenced isolates (28 clinical and 27 environmental), the gene content for L. longbeachae from growing media samples (3,586 genes) was significantly higher than that for isolates from patients (3,533 genes; 2-sample t-test, t = 2.5213; d.f. = 53; p = 0.01474) (Figure 4, panel B). The data imply that gene loss occurs during human infection or that L. longbeachae strains with reduced gene content have enhanced human infectivity. However, we did not identify a specific enriched gene or functional category in clinical or environmental samples (data not shown).

Source Attribution Confounded by Complex Serogroup 1 Populations within Environmental Samples

Having accounted for the influence of recombination on the phylogeny of L. longbeachae, we investigated the diversity of isolates associated with 5 patients and their linked compost samples obtained during 2008–2014, including 3 patients from the 2013 outbreak in Scotland. Of note, isolates from the 2013 outbreak were distributed across several subclades of the tree, indicating that the infections were caused by different strains (Figure 2). However, all isolates from a single patient clustered together, consistent with a monoclonal etiology of each infection. Of note, for all 5 patients, clinical isolates were not closely allied to the environmental isolates obtained from linked compost samples, and therefore a genetic link between patient and compost samples could not be established. Most subclades included isolates of diverse geographic origin, consistent with a wide distribution for L. longbeachae strains; however, 3 L. longbeachae isolates originating from Australasia (strains 13.8294, 13.8293, and NSW150) belonged to their own region-specific cluster (Figure 2).

We hypothesized that the lack of genetic relatedness between L. longbeachae isolates from patients and linked compost samples could be explained by a highly diverse population of L. longbeachae in growing media samples compounded by a sampling strategy consisting of a single sequenced isolate. All 5 compost samples for which we had >1 isolate contained isolates distributed across multiple clades in the phylogenetic tree. In particular, 5 isolates from the same growing media sample linked to a patient infected in Edinburgh in 2014 were distributed across 4 distinct clades, demonstrating that within a single environmental sample, considerable species diversity may be represented (Figure 2). Taken together, these data suggest that for future outbreak investigations, extensive sampling of environmental samples may be required to identify genotypes responsible for episodes of legionellosis infection, if indeed they are present.

Discussion

Our findings reveal the population genomic structure for L. longbeachae, an emerging pathogen in Europe and the United States, and includes a genome-scale investigation into an outbreak of L. longbeachae legionellosis. We provide evidence for extensive recombination and lateral gene transfer among L. longbeachae, including the presence of widely distributed mosaic plasmids that have likely recombined with plasmids from other Legionella spp., suggesting an ecologic overlap or shared habitat. Our analysis highlights the need to account for recombination events when determining the genetic relatedness of L. longbeachae isolates.

Our application of whole-genome sequencing for diagnostic purposes revealed the misidentification, using current serotyping methods, of several L. anisa isolates as L. longbeachae and led to the identification of a putative novel Legionella sp. linked to legionellosis. These findings highlight the limitations of current typing methods for differentiation of Legionella spp. and accurate identification of legionellosis etiology.

We used whole-genome sequencing to attempt to establish a genetic link between legionellosis infections and associated compost samples. Our inability to establish a link probably reflects the traditional strategy of single isolate sampling, which when applied to a highly diverse pool of L. longbeachae genotypes fails to detect the infecting genotype. We suggest that the approach to investigating the source of future legionellosis cases linked to growing media will require a radical revision of sampling protocols to maximize the chances of isolating the infecting strain, if present. Taken together, our findings provide a view of the population structure of L. longbeachae and highlight the complexities of tracing the origin of legionellosis associated with growing media. Overall, our findings demonstrate the resolution afforded by whole-genome sequencing for understanding the biology underpinning legionellosis and provide information that should be considered for future epidemiologic investigations.

Mr. Bacigalupe is a PhD candidate at the Roslin Institute, University of Edinburgh. His primary research focuses on the evolution, adaptation, and outbreak dynamics of bacterial pathogens.

Top

Acknowledgments

We are grateful to Carmen Buchrieser for providing the original sequence reads for L. longbeachae strains ATCC39642, 98072, and C-4E7. We thank David Harte for supplying the New Zealand strains.

Funding was provided by the Chief Scientist’s Office Scotland (grant ETM/421 to J.R.F.) and by the Biotechnology and Biological Sciences Research Council (ISP3 grant BB/J004227/1 to J.R.F.).

Top

References

  1. Fields  BS, Benson  RF, Besser  RE. Legionella and Legionnaires’ disease: 25 years of investigation. Clin Microbiol Rev. 2002;15:50626. DOIPubMed
  2. European Centre for Disease Prevention and Control. Surveillance report. Legionnaires’ disease in Europe, 2010. 2012 [cited 2016 Jul 9]. http://ecdc.europa.eu/en/publications/publications/sur-legionnaires-disease-surveillance-2010.pdf
  3. Joseph  CA, Ricketts  KD; European Working Group for Legionella Infections. Legionnaires disease in Europe 2007-2008. Euro Surveill. 2010;15:19493.PubMed
  4. Marston  BJ, Lipman  HB, Breiman  RF. Surveillance for Legionnaires’ disease. Risk factors for morbidity and mortality. Arch Intern Med. 1994;154:241722. DOIPubMed
  5. Li  JS, O’Brien  ED, Guest  C. A review of national legionellosis surveillance in Australia, 1991 to 2000. Commun Dis Intell Q Rep. 2002;26:4618.PubMed
  6. Cramp  GJ, Harte  D, Douglas  NM, Graham  F, Schousboe  M, Sykes  K. An outbreak of Pontiac fever due to Legionella longbeachae serogroup 2 found in potting mix in a horticultural nursery in New Zealand. Epidemiol Infect. 2010;138:1520. DOIPubMed
  7. Whiley  H, Bentham  R. Legionella longbeachae and legionellosis. Emerg Infect Dis. 2011;17:57983. DOIPubMed
  8. Yu  VL, Plouffe  JF, Pastoris  MC, Stout  JE, Schousboe  M, Widmer  A, et al. Distribution of Legionella species and serogroups isolated by culture in patients with sporadic community-acquired legionellosis: an international collaborative survey. J Infect Dis. 2002;186:1278. DOIPubMed
  9. García  C, Ugalde  E, Campo  AB, Miñambres  E, Kovács  N. Fatal case of community-acquired pneumonia caused by Legionella longbeachae in a patient with systemic lupus erythematosus. Eur J Clin Microbiol Infect Dis. 2004;23:1168. DOIPubMed
  10. den Boer  JW, Yzerman  EPF, Jansen  R, Bruin  JP, Verhoef  LPB, Neve  G, et al. Legionnaires’ disease and gardening. Clin Microbiol Infect. 2007;13:8891. DOIPubMed
  11. Potts  A, Donaghy  M, Marley  M, Othieno  R, Stevenson  J, Hyland  J, et al. Cluster of Legionnaires disease cases caused by Legionella longbeachae serogroup 1, Scotland, August to September 2013. Euro Surveill. 2013;18:20656. DOIPubMed
  12. Lindsay  DSJ, Brown  AW, Brown  DJ, Pravinkumar  SJ, Anderson  E, Edwards  GF. Legionella longbeachae serogroup 1 infections linked to potting compost. J Med Microbiol. 2012;61:21822. DOIPubMed
  13. Steele  TW, Lanser  J, Sangster  N. Isolation of Legionella longbeachae serogroup 1 from potting mixes. Appl Environ Microbiol. 1990;56:4953.PubMed
  14. Koide  M, Arakaki  N, Saito  A. Distribution of Legionella longbeachae and other legionellae in Japanese potting soils. J Infect Chemother. 2001;7:2247. DOIPubMed
  15. Reuter  S, Harrison  TG, Köser  CU, Ellington  MJ, Smith  GP, Parkhill  J, et al. A pilot study of rapid whole-genome sequencing for the investigation of a Legionella outbreak. BMJ Open. 2013;3:e002175. DOIPubMed
  16. Rao  C, Benhabib  H, Ensminger  AW. Phylogenetic reconstruction of the Legionella pneumophila Philadelphia-1 laboratory strains through comparative genomics. PLoS One. 2013;8:e64129. DOIPubMed
  17. Cazalet  C, Gomez-Valero  L, Rusniok  C, Lomma  M, Dervins-Ravault  D, Newton  HJ, et al. Analysis of the Legionella longbeachae genome and transcriptome uncovers unique strategies to cause Legionnaires’ disease. PLoS Genet. 2010;6:e1000851. DOIPubMed
  18. Kozak  NA, Buss  M, Lucas  CE, Frace  M, Govil  D, Travis  T, et al. Virulence factors encoded by Legionella longbeachae identified on the basis of the genome sequence analysis of clinical isolate D-4968. J Bacteriol. 2010;192:103044. DOIPubMed
  19. Gomez-Valero  L, Rusniok  C, Jarraud  S, Vacherie  B, Rouy  Z, Barbe  V, et al. Extensive recombination events and horizontal gene transfer shaped the Legionella pneumophila genomes. BMC Genomics. 2011;12:536. DOIPubMed
  20. Ratcliff  RM, Lanser  JA, Manning  PA, Heuzenroeder  MW. Sequence-based classification scheme for the genus Legionella targeting the mip gene. J Clin Microbiol. 1998;36:15607.PubMed
  21. Fallon  RJ, Abraham  WH. Experience with heat-killed antigens of L. longbeachae serogroups 1 and 2, and L. jordanis in the indirect fluorescence antibody test. Zentralbl Bakteriol Mikrobiol Hyg A. 1983;255:814.PubMed
  22. Babraham Bioinformatics. FastQC. 2010 [cited 2016 Jul 9]. http://www.bioinformatics.babraham.ac.uk/projects/fastqc
  23. Martin  M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal. 2011;17:10–12.
  24. GitHub, Inc. wgsim Read Simulator [cited 2016 Jul 9]. https://github.com/lh3/wgsim
  25. Cole  JR, Wang  Q, Fish  JA, Chai  B, McGarrell  DM, Sun  Y, et al. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014;42(D1):D63342. DOIPubMed
  26. Kim  M, Oh  HS, Park  SC, Chun  J. Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int J Syst Evol Microbiol. 2014;64:34651. DOIPubMed
  27. Health Protection Scotland. Surveillance report: legionellosis in Scotland 2013–2014. 2015 Sep 1 [cited 2015 Aug 15]. http://www.hps.scot.nhs.uk/resp/wrdetail.aspx?id=65135&wrtype=6
  28. van der Mee-Marquet  N, Domelier  AS, Arnault  L, Bloc  D, Laudat  P, Hartemann  P, et al. Legionella anisa, a possible indicator of water contamination by Legionella pneumophila. J Clin Microbiol. 2006;44:569. DOIPubMed
  29. Svarrer  CW, Uldum  SA. The occurrence of Legionella species other than Legionella pneumophila in clinical and environmental samples in Denmark identified by mip gene sequencing and matrix-assisted laser desorption ionization time-of-flight mass spectrometry. Clin Microbiol Infect. 2012;18:10049. DOIPubMed
  30. Underwood  AP, Jones  G, Mentasti  M, Fry  NK, Harrison  TG. Comparison of the Legionella pneumophila population structure as determined by sequence-based typing and whole genome sequencing. BMC Microbiol. 2013;13:302. DOIPubMed
  31. Huson  DH, Bryant  D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:25467. DOIPubMed
  32. Bruen  TC, Philippe  H, Bryant  D. A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006;172:266581. DOIPubMed
  33. Marttinen  P, Hanage  WP, Croucher  NJ, Connor  TR, Harris  SR, Bentley  SD, et al. Detection of recombination events in bacterial genomes from large population samples. Nucleic Acids Res. 2012;40:e6. DOIPubMed
  34. de Been  M, van Schaik  W, Cheng  L, Corander  J, Willems  RJ. Recent recombination events in the core genome are associated with adaptive evolution in Enterococcus faecium. Genome Biol Evol. 2013;5:152435. DOIPubMed
  35. Didelot  X, Wilson  DJ. ClonalFrameML: efficient inference of recombination in whole bacterial genomes. PLOS Comput Biol. 2015;11:e1004041. DOIPubMed
  36. Coscollá  M, Comas  I, González-Candelas  F. Quantifying nonvertical inheritance in the evolution of Legionella pneumophila. Mol Biol Evol. 2011;28:9851001. DOIPubMed
  37. Sánchez-Busó  L, Comas  I, Jorques  G, González-Candelas  F. Recombination drives genome evolution in outbreak-related Legionella pneumophila isolates. Nat Genet. 2014;46:120511. DOIPubMed
  38. Lanza  VF, de Toro  M, Garcillán-Barcia  MP, Mora  A, Blanco  J, Coque  TM, et al. Plasmid flux in Escherichia coli ST131 sublineages, analyzed by plasmid constellation network (PLACNET), a new method for plasmid reconstruction from whole genome sequences. PLoS Genet. 2014;10:e1004766. DOIPubMed
  39. Cazalet  C, Rusniok  C, Brüggemann  H, Zidane  N, Magnier  A, Ma  L, et al. Evidence in the Legionella pneumophila genome for exploitation of host cell functions and high genome plasticity. Nat Genet. 2004;36:116573. DOIPubMed

Top

Figures

Top

Cite This Article

DOI: 10.3201/eid2305.161165

Table of Contents – Volume 23, Number 5—May 2017

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

J. Ross Fitzgerald, The Roslin Institute, University of Edinburgh, Midlothian, EH259RG Scotland, UK

Send To

character(s) remaining.

Comment submitted successfully, thank you for your feedback.

Top

Page created: April 14, 2017
Page updated: April 14, 2017
Page reviewed: April 14, 2017
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external