Recombination drives genome evolution in outbreak-related Legionella pneumophila isolates

Legionella pneumophila is a strictly environmental pathogen and the etiological agent of legionellosis. It is known that non-vertical processes have a major role in the short-term evolution of pathogens, but little is known about the relevance of these and other processes in environmental bacteria. We report the whole-genome sequencing of 69 L. pneumophila strains linked to recurrent outbreaks in a single location (Alcoy, Spain) over 11 years. We found some examples where the genome sequences of isolates of the same sequence type and outbreak did not cluster together and were more closely related to sequences from different outbreaks. Our analyses identify 16 recombination events responsible for almost 98% of the SNPs detected in the core genome and an apparent acceleration in the evolutionary rate. These results have profound implications for the understanding of microbial populations and for public health interventions in Legionella outbreak investigations.

[1]  Teresa M. Coque,et al.  Antibiotic resistance shaping multi-level population biology of bacteria , 2013, Front. Microbiol..

[2]  N. Fry,et al.  Addition of neuA, the Gene Encoding N-Acylneuraminate Cytidylyl Transferase, Increases the Discriminatory Ability of the Consensus Sequence-Based Scheme for Typing Legionella pneumophila Serogroup 1 Strains , 2007, Journal of Clinical Microbiology.

[3]  A. Wilm,et al.  LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets , 2012, Nucleic acids research.

[4]  Julian Parkhill,et al.  Whole-genome sequencing for analysis of an outbreak of meticillin-resistant Staphylococcus aureus: a descriptive study , 2013, The Lancet. Infectious Diseases.

[5]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[6]  Jukka Corander,et al.  Dense genomic sampling identifies highways of pneumococcal recombination , 2014, Nature Genetics.

[7]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[8]  P. Glaser,et al.  Shaping a bacterial genome by large chromosomal replacements, the evolutionary history of Streptococcus agalactiae , 2008, Proceedings of the National Academy of Sciences.

[9]  Mark Achtman,et al.  Evolutionary History of Salmonella Typhi , 2006, Science.

[10]  M. Pirmohamed,et al.  Emergence and global spread of epidemic healthcare-associated Clostridium difficile , 2012, Nature Genetics.

[11]  M. Suchard,et al.  Bayesian Phylogenetics with BEAUti and the BEAST 1.7 , 2012, Molecular biology and evolution.

[12]  M. Glickman,et al.  Mycobacterium tuberculosis nuoG Is a Virulence Gene That Inhibits Apoptosis of Infected Host Cells , 2007, PLoS pathogens.

[13]  Carl T. Bergstrom,et al.  Bacteria are different: observations, interpretations, speculations, and opinions about the mechanisms of adaptive evolution in prokaryotes. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[14]  Ruifu Yang,et al.  Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis , 2012, Proceedings of the National Academy of Sciences.

[15]  S. Fontana,et al.  Cluster of travel-associated Legionnaires disease in Lazise, Italy, July to August 2011. , 2011, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[16]  V. Laza-Stanca,et al.  An outbreak of Legionnaires' disease associated with a display spa pool in retail premises, Stoke-on-Trent, United Kingdom, July 2012. , 2012, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[17]  Steven J. M. Jones,et al.  Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. , 2011, The New England journal of medicine.

[18]  C. Lück,et al.  Extension of the Legionella pneumophila sequence-based typing scheme to include strains carrying a variant of the N-acylneuraminate cytidylyltransferase gene. , 2014, Clinical microbiology and infection : the official publication of the European Society of Clinical Microbiology and Infectious Diseases.

[19]  J. Etienne,et al.  Consensus Sequence-Based Scheme for Epidemiological Typing of Clinical and Environmental Isolates of Legionella pneumophila , 2005, Journal of Clinical Microbiology.

[20]  C. Fraser,et al.  Fuzzy species among recombinogenic bacteria , 2005, BMC Biology.

[21]  E. Feil,et al.  A Dominant Clone of Leptospira interrogans Associated with an Outbreak of Human Leptospirosis in Thailand , 2007, PLoS neglected tropical diseases.

[22]  F. Blattner,et al.  Mauve: multiple alignment of conserved genomic sequence with rearrangements. , 2004, Genome research.

[23]  Mihai Pop,et al.  Bacillus anthracis comparative genome analysis in support of the Amerithrax investigation , 2011, Proceedings of the National Academy of Sciences.

[24]  P. Alifano,et al.  Hypermutation in pathogenic bacteria: frequent phase variation in meningococci is a phenotypic trait of a specialized mutator biotype. , 1999, Molecular cell.

[25]  J. Ernst,et al.  Mycobacterium tuberculosis inhibits neutrophil apoptosis, leading to delayed activation of naive CD4 T cells. , 2012, Cell host & microbe.

[26]  E. Birney,et al.  Velvet: algorithms for de novo short read assembly using de Bruijn graphs. , 2008, Genome research.

[27]  X. Didelot,et al.  A comparison of homologous recombination rates in bacteria and archaea , 2009, The ISME Journal.

[28]  James H. Bullard,et al.  The origin of the Haitian cholera outbreak strain. , 2011, The New England journal of medicine.

[29]  J. Townsend,et al.  Assessing the Probability of Detection of Horizontal Gene Transfer Events in Bacterial Populations , 2012, Front. Microbio..

[30]  Alexandros Stamatakis,et al.  RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models , 2006, Bioinform..

[31]  Andries J. van Tonder,et al.  Evolutionary dynamics of Clostridium difficile over short and long time scales , 2010, Proceedings of the National Academy of Sciences.

[32]  F. González-Candelas,et al.  Quantifying nonvertical inheritance in the evolution of Legionella pneumophila. , 2011, Molecular biology and evolution.

[33]  F. González-Candelas,et al.  Legionellosis Outbreak Associated with Asphalt Paving Machine, Spain, 2009 , 2010, Emerging infectious diseases.

[34]  G. Smith,et al.  A pilot study of rapid whole-genome sequencing for the investigation of a Legionella outbreak , 2013, BMJ Open.

[35]  A. Underwood,et al.  Comparison of the Legionella pneumophila population structure as determined by sequence-based typing and whole genome sequencing , 2013, BMC Microbiology.

[36]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[37]  G F Mallison,et al.  Legionnaires' disease: description of an epidemic of pneumonia. , 1977, The New England journal of medicine.

[38]  F. González-Candelas,et al.  Population structure and recombination in environmental isolates of Legionella pneumophila. , 2007, Environmental Microbiology.

[39]  Jessica L Miller,et al.  The Type I NADH Dehydrogenase of Mycobacterium tuberculosis Counters Phagosomal NOX2 Activity to Inhibit TNF-α-Mediated Host Cell Apoptosis , 2010, PLoS pathogens.

[40]  Barry S. Fields,et al.  Legionella and Legionnaires' Disease: 25 Years of Investigation , 2002, Clinical Microbiology Reviews.

[41]  S. Thorn,et al.  Public health response to an outbreak of Legionnaires' disease in Edinburgh, United Kingdom, June 2012. , 2012, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[42]  Julian Parkhill,et al.  Rapid whole-genome sequencing for investigation of a neonatal MRSA outbreak. , 2012, The New England journal of medicine.

[43]  J. Parkhill,et al.  Whole-genome sequencing to identify transmission of Mycobacterium abscessus between patients with cystic fibrosis: a retrospective cohort study , 2013, The Lancet.

[44]  Jukka Corander,et al.  Hyper-Recombination, Diversity, and Antibiotic Resistance in Pneumococcus , 2009, Science.

[45]  D. W. Kim,et al.  Shigella sonnei genome sequencing and phylogenetic analysis indicate recent global dissemination from Europe , 2012, Nature Genetics.

[46]  S. Mondot,et al.  Population Structure of the Fish-Pathogenic Bacterium Flavobacterium psychrophilum , 2008, Applied and Environmental Microbiology.

[47]  C. Buchrieser,et al.  Genome dynamics in Legionella: the basis of versatility and adaptation to intracellular replication. , 2013, Cold Spring Harbor perspectives in medicine.

[48]  Albert J. Vilella,et al.  VariScan: Analysis of evolutionary patterns from large-scale DNA sequence polymorphism data , 2005, Bioinform..

[49]  N. Perna,et al.  progressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement , 2010, PloS one.

[50]  W. Doolittle,et al.  On the origin of prokaryotic species. , 2009, Genome research.

[51]  Toni Gabaldón,et al.  trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses , 2009, Bioinform..

[52]  J. Gogarten,et al.  Horizontal gene transfer : genomes in flux , 2009 .

[53]  Korbinian Strimmer,et al.  APE: Analyses of Phylogenetics and Evolution in R language , 2004, Bioinform..

[54]  G. Sundin,et al.  General and inducible hypermutation facilitate parallel adaptation in Pseudomonas aeruginosa despite divergent mutation spectra , 2012, Proceedings of the National Academy of Sciences.

[55]  J. Rothberg,et al.  Prospective Genomic Characterization of the German Enterohemorrhagic Escherichia coli O104:H4 Outbreak by Rapid Next Generation Sequencing Technology , 2011, PloS one.

[56]  P. Yiallouros,et al.  First outbreak of nosocomial Legionella infection in term neonates caused by a cold mist ultrasonic humidifier. , 2013, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[57]  N. Fry,et al.  Sequence-Based Typing of Legionella pneumophila Serogroup 1 Offers the Potential for True Portability in Legionellosis Outbreak Investigation , 2003, Journal of Clinical Microbiology.

[58]  Andrés Moya,et al.  Legionella pneumophila pangenome reveals strain-specific virulence factors , 2010, BMC Genomics.

[59]  M. Cetron,et al.  A community outbreak of Legionnaires' disease linked to hospital cooling towers: an epidemiological method to calculate dose of exposure. , 1999, International journal of epidemiology.

[60]  M. Evans,et al.  A community outbreak of Legionnaires' disease in South Wales, August-September 2010. , 2010, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[61]  J. Burton,et al.  Rapid Pneumococcal Evolution in Response to Clinical Interventions , 2011, Science.

[62]  Ning Ma,et al.  BLAST+: architecture and applications , 2009, BMC Bioinformatics.

[63]  Omar E. Cornejo,et al.  The Population and Evolutionary Dynamics of Homologous Gene Recombination in Bacteria , 2009, PLoS genetics.

[64]  Martin Vingron,et al.  TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing , 2002, Bioinform..

[65]  M. Lipsitch,et al.  Population genomics of post-vaccine changes in pneumococcal epidemiology , 2013, Nature Genetics.

[66]  D. Orozco,et al.  Clinical Study of an Outbreak of Legionnaire's Disease in Alcoy, Southeastern Spain , 2002, European Journal of Clinical Microbiology and Infectious Diseases.