Genomic analysis of Chlamydia psittaci from a multistate zoonotic outbreak in two chicken processing plants.

Four Chlamydia psittaci isolates were recovered from clinical specimens from ill workers during a multistate outbreak at two chicken processing plants. Whole genome sequencing analyses revealed high similarity to C. psittaci genotype D. The isolates differed from each other by only two single nucleotide polymorphisms, indicating a common source.

Corrigendum: Phylogeny of Vibrio vulnificus From the Analysis of the Core-Genome: Implications for Intra-Species Taxonomy.

[This corrects the article DOI: 10.3389/fmicb.2017.02613.].

Comparative genome analysis reveals a complex population structure of Legionella pneumophila subspecies.

The majority of Legionnaires' disease (LD) cases are caused by Legionella pneumophila, a genetically heterogeneous species composed of at least 17 serogroups. Previously, it was demonstrated that L. pneumophila consists of three subspecies: pneumophila, fraseri and pascullei. During an LD outbreak investigation in 2012, we detected that representatives of both subspecies fraseri and pascullei colonized the same water system and that the outbreak-causing strain was a new member of the least represented subspecies pascullei. We used partial sequence based typing consensus patterns to mine an international database for additional representatives of fraseri and pascullei subspecies. As a result, we identified 46 sequence types (STs) belonging to subspecies fraseri and two STs belonging to subspecies pascullei. Moreover, a recent retrospective whole genome sequencing analysis of isolates from New York State LD clusters revealed the presence of a fourth L. pneumophila subspecies that we have termed raphaeli. This subspecies consists of 15 STs. Comparative analysis was conducted using the genomes of multiple members of all four L. pneumophila subspecies. Whereas each subspecies forms a distinct phylogenetic clade within the L. pneumophila species, they share more average nucleotide identity with each other than with other Legionella species. Unique genes for each subspecies were identified and could be used for rapid subspecies detection. Improved taxonomic classification of L. pneumophila strains may help identify environmental niches and virulence attributes associated with these genetically distinct subspecies.

Development of a multiplex TaqMan real-time PCR assay for the detection of Chlamydia psittaci and Chlamydia pneumoniae in human clinical specimens.

Diagnosis of Chlamydia psittaci and Chlamydia pneumoniae infections has traditionally relied on serological assays. We developed a multiplex real-time PCR assay for detection of C. psittaci, C. pneumoniae and an internal control. Results of this assay demonstrated 100% concordance compared to results of previously tested human clinical specimens.

Correction: Comprehensive bioinformatics analysis of Mycoplasma pneumoniae genomes to investigate underlying population structure and type-specific determinants.

[This corrects the article DOI: 10.1371/journal.pone.0174701.].

Legionnaires' Disease Outbreak Caused by Endemic Strain of Legionella pneumophila, New York, New York, USA, 2015.

During the summer of 2015, New York, New York, USA, had one of the largest and deadliest outbreaks of Legionnaires' disease in the history of the United States. A total of 138 cases and 16 deaths were linked to a single cooling tower in the South Bronx. Analysis of environmental samples and clinical isolates showed that sporadic cases of legionellosis before, during, and after the outbreak could be traced to a slowly evolving, single-ancestor strain. Detection of an ostensibly virulent Legionella strain endemic to the Bronx community suggests potential risk for future cases of legionellosis in the area. The genetic homogeneity of the Legionella population in this area might complicate investigations and interpretations of future outbreaks of Legionnaires' disease.

Characterization of Legionella Species from Watersheds in British Columbia, Canada.

Legionella spp. present in some human-made water systems can cause Legionnaires' disease in susceptible individuals. Although legionellae have been isolated from the natural environment, variations in the organism's abundance over time and its relationship to aquatic microbiota are poorly understood. Here, we investigated the presence and diversity of legionellae through 16S rRNA gene amplicon and metagenomic sequencing of DNA from isolates collected from seven sites in three watersheds with varied land uses over a period of 1 year. Legionella spp. were found in all watersheds and sampling sites, comprising up to 2.1% of the bacterial community composition. The relative abundance of Legionella tended to be higher in pristine sites than in sites affected by agricultural activity. The relative abundance levels of Amoebozoa, some of which are natural hosts of legionellae, were similarly higher in pristine sites. Compared to other bacterial genera detected, Legionella had both the highest richness and highest alpha diversity. Our findings indicate that a highly diverse population of legionellae may be found in a variety of natural aquatic sources. Further characterization of these diverse natural populations of Legionella will help inform prevention and control efforts aimed at reducing the risk of Legionella colonization of built environments, which could ultimately decrease the risk of human disease. IMPORTANCE Many species of Legionella can cause Legionnaires' disease, a significant cause of bacterial pneumonia. Legionella in human-made water systems such as cooling towers and building plumbing systems are the primary sources of Legionnaires' disease outbreaks. In this temporal study of natural aquatic environments, Legionella relative abundance was shown to vary in watersheds associated with different land uses. Analysis of the Legionella sequences detected at these sites revealed highly diverse populations that included potentially novel Legionella species. These findings have important implications for understanding the ecology of Legionella and control measures for this pathogen that are aimed at reducing human disease.

Mycoplasma hominis Infections Transmitted Through Amniotic Tissue Product.

Background: Mycoplasma hominis is a commensal genitourinary tract organism that can cause infections outside the genitourinary tract. We investigated a cluster of M. hominis surgical site infections in patients who underwent spine surgery, all associated with amniotic tissue linked to a common donor. Methods: Laboratory tests of tissue product from the donor, including culture, quantitative real-time polymerase chain reaction (qPCR), and whole-genome sequencing were performed. Use of this amniotic tissue product was reviewed. A multistate investigation to identify additional cases and locate any unused products was conducted. Results: Twenty-seven tissue product vials from a donor were distributed to facilities in 7 states; at least 20 vials from this donor were used in 14 patients. Of these, 4 of 14 (29%) developed surgical site infections, including 2 M. hominis infections. Mycoplasma hominis was detected by culture and qPCR in 2 unused vials from the donor. Sequencing indicated >99% similarity between patient and unopened vial isolates. For 5 of 27 (19%) vials, the final disposition could not be confirmed. Conclusions: Mycoplasma hominis was transmitted through amniotic tissue from a single donor to 2 recipients. Current routine donor screening and product testing does not detect all potential pathogens. Clinicians should be aware that M. hominis can cause surgical site infections, and may not be detected by routine clinical cultures. The lack of a standardized system to track tissue products in healthcare facilities limits the ability of public health agencies to respond to outbreaks and investigate other adverse events associated with these products.

Comprehensive bioinformatics analysis of Mycoplasma pneumoniae genomes to investigate underlying population structure and type-specific determinants.

Mycoplasma pneumoniae is a significant cause of respiratory illness worldwide. Despite a minimal and highly conserved genome, genetic diversity within the species may impact disease. We performed whole genome sequencing (WGS) analysis of 107 M. pneumoniae isolates, including 67 newly sequenced using the Pacific BioSciences RS II and/or Illumina MiSeq sequencing platforms. Comparative genomic analysis of 107 genomes revealed >3,000 single nucleotide polymorphisms (SNPs) in total, including 520 type-specific SNPs. Population structure analysis supported the existence of six distinct subgroups, three within each type. We developed a predictive model to classify an isolate based on whole genome SNPs called against the reference genome into the identified subtypes, obviating the need for genome assembly. This study is the most comprehensive WGS analysis for M. pneumoniae to date, underscoring the power of combining complementary sequencing technologies to overcome difficult-to-sequence regions and highlighting potential differential genomic signatures in M. pneumoniae.

Complete Genome Sequences of Legionella pneumophila subsp. fraseri Strains Detroit-1 and Dallas 1E.

We report here the complete genome sequences of two of the earliest known strains of Legionella pneumophila subsp. fraseri Detroit-1 is serogroup 1 and was isolated from a lung biopsy specimen in 1977. Dallas 1E is serogroup 5 and was isolated in 1978 from a cooling tower.

Complete Genome Sequence of Mycoplasma pneumoniae Type 2 Reference Strain FH Using Single-Molecule Real-Time Sequencing Technology.

Mycoplasma pneumoniae type 2 strain FH was previously sequenced with Illumina (FH-Illumina) and 454 (FH-454) technologies according to Xiao et al. (2015) and Krishnakumar et al. (2010). Comparative analyses revealed differences in genomic content between these sequences, including a 6-kb region absent from the FH-454 submission. Here, we present a complete genome sequence of FH sequenced with the Pacific Biosciences RSII platform.

Phylogeny of Vibrio vulnificus from the Analysis of the Core-Genome: Implications for Intra-Species Taxonomy.

Vibrio vulnificus (Vv) is a multi-host pathogenic species currently subdivided into three biotypes (Bts). The three Bts are human-pathogens, but only Bt2 is also a fish-pathogen, an ability that is conferred by a transferable virulence-plasmid (pVvbt2). Here we present a phylogenomic analysis from the core genome of 80 Vv strains belonging to the three Bts recovered from a wide range of geographical and ecological sources. We have identified five well-supported phylogenetic groups or lineages (L). L1 comprises a mixture of clinical and environmental Bt1 strains, most of them involved in human clinical cases related to raw seafood ingestion. L2 is formed by a mixture of Bt1 and Bt2 strains from various sources, including diseased fish, and is related to the aquaculture industry. L3 is also linked to the aquaculture industry and includes Bt3 strains exclusively, mostly related to wound infections or secondary septicemia after farmed-fish handling. Lastly, L4 and L5 include a few strains of Bt1 associated with specific geographical areas. The phylogenetic trees for ChrI and II are not congruent to one another, which suggests that inter- and/or intra-chromosomal rearrangements have been produced along Vv evolution. Further, the phylogenetic trees for each chromosome and the virulence plasmid were also not congruent, which also suggests that pVvbt2 has been acquired independently by different clones, probably in fish farms. From all these clones, the one with zoonotic capabilities (Bt2-Serovar E) has successfully spread worldwide. Based on these results, we propose a new updated classification of the species based on phylogenetic lineages rather than on Bts, as well as the inclusion of all Bt2 strains in a pathovar with the particular ability to cause fish vibriosis, for which we suggest the name "piscis."

Development of a multiplex taqMan real-time PCR assay for typing of Mycoplasma pneumoniae based on type-specific indels identified through whole genome sequencing.

We developed a multiplex real-time PCR assay for simultaneously detecting M. pneumoniae and typing into historically-defined P1 types. Typing was achieved based on the presence of short type-specific indels identified through whole genome sequencing. This assay was 100% specific compared to existing methods and may be useful during epidemiologic investigations.

Dynamics of genome change among Legionella species.

Legionella species inhabit freshwater and soil ecosystems where they parasitize protozoa. L. pneumonphila (LP) serogroup-1 (Lp1) is the major cause of Legionnaires' Disease (LD), a life-threatening pulmonary infection that can spread systemically. The increased global frequency of LD caused by Lp and non-Lp species underscores the need to expand our knowledge of evolutionary forces underlying disease pathogenesis. Whole genome analyses of 43 strains, including all known Lp serogroups 1-17 and 17 emergent LD-causing Legionella species (of which 33 were sequenced in this study) in addition to 10 publicly available genomes, resolved the strains into four phylogenetic clades along host virulence demarcations. Clade-specific genes were distinct for genetic exchange and signal-transduction, indicating adaptation to specific cellular and/or environmental niches. CRISPR spacer comparisons hinted at larger pools of accessory DNA sequences in Lp than predicted by the pan-genome analyses. While recombination within Lp was frequent and has been reported previously, population structure analysis identified surprisingly few DNA admixture events between species. In summary, diverse Legionella LD-causing species share a conserved core-genome, are genetically isolated from each other, and selectively acquire genes with potential for enhanced virulence.

Genomic Analysis Reveals Novel Diversity among the 1976 Philadelphia Legionnaires' Disease Outbreak Isolates and Additional ST36 Strains.

Legionella pneumophila was first recognized as a cause of severe and potentially fatal pneumonia during a large-scale outbreak of Legionnaires' disease (LD) at a Pennsylvania veterans' convention in Philadelphia, 1976. The ensuing investigation and recovery of four clinical isolates launched the fields of Legionella epidemiology and scientific research. Only one of the original isolates, "Philadelphia-1", has been widely distributed or extensively studied. Here we describe the whole-genome sequencing (WGS), complete assembly, and comparative analysis of all Philadelphia LD strains recovered from that investigation, along with L. pneumophila isolates sharing the Philadelphia sequence type (ST36). Analyses revealed that the 1976 outbreak was due to multiple serogroup 1 strains within the same genetic lineage, differentiated by an actively mobilized, self-replicating episome that is shared with L. pneumophila str. Paris, and two large, horizontally-transferred genomic loci, among other polymorphisms. We also found a completely unassociated ST36 strain that displayed remarkable genetic similarity to the historical Philadelphia isolates. This similar strain implies the presence of a potential clonal population, and suggests important implications may exist for considering epidemiological context when interpreting phylogenetic relationships among outbreak-associated isolates. Additional extensive archival research identified the Philadelphia isolate associated with a non-Legionnaire case of "Broad Street pneumonia", and provided new historical and genetic insights into the 1976 epidemic. This retrospective analysis has underscored the utility of fully-assembled WGS data for Legionella outbreak investigations, highlighting the increased resolution that comes from long-read sequencing and a sequence type-matched genomic data set.

Complete Genome Sequences of the Historical Legionella pneumophila Strains OLDA and Pontiac.

Here, we report the complete genome sequences of Legionella pneumophila serogroup 1 strains OLDA and Pontiac, which predate the 1976 Philadelphia Legionnaires' disease outbreak. Strain OLDA was isolated in 1947 from an apparent sporadic case, and strain Pontiac caused an explosive outbreak at a Michigan health department in 1968.

Complete Genome Sequences of Three Outbreak-Associated Legionella pneumophila Isolates.

We report here the complete genome sequences of three Legionella pneumophila isolates that are associated with a Legionnaires' disease outbreak in New York in 2012. Two clinical isolates (D7630 and D7632) and one environmental isolate (D7631) were recovered from this outbreak. A single isolate-specific virulence gene was found in D7632. These isolates were included in a large study evaluating the genomic resolution of various bioinformatics approaches for L. pneumophila serogroup 1 isolates.

Three Genome Sequences of Legionella pneumophila subsp. pascullei Associated with Colonization of a Health Care Facility.

Here, we report the complete genome sequences of three Legionella pneumophila subsp. pascullei strains (including both serogroup 1 and 5 strains) that were found in the same health care facility in 1982 and 2012.

Genomic Resolution of Outbreak-Associated Legionella pneumophila Serogroup 1 Isolates from New York State.

UNLABELLED: A total of 30 Legionella pneumophila serogroup 1 isolates representing 10 separate legionellosis laboratory investigations ("outbreaks") that occurred in New York State between 2004 and 2012 were selected for evaluation of whole-genome sequencing (WGS) approaches for molecular subtyping of this organism. Clinical and environmental isolates were available for each outbreak and were initially examined by pulsed-field gel electrophoresis (PFGE). Sequence-based typing alleles were extracted from WGS data yielding complete sequence types (ST) for isolates representing 8 out of the 10 outbreaks evaluated in this study. Isolates from separate outbreaks sharing the same ST also contained the fewest differences in core genome single nucleotide polymorphisms (SNPs) and the greatest proportion of identical allele sequences in a whole-genome multilocus sequence typing (wgMLST) scheme. Both core SNP and wgMLST analyses distinguished isolates from separate outbreaks, including those from two outbreaks sharing indistinguishable PFGE profiles. Isolates from a hospital-associated outbreak spanning multiple years shared indistinguishable PFGE profiles but displayed differences in their genome sequences, suggesting the presence of multiple environmental sources. Finally, the rtx gene demonstrated differences in the repeat region sequence among ST1 isolates from different outbreaks, suggesting that variation in this gene may be useful for targeted molecular subtyping approaches for L. pneumophila This study demonstrates the utility of various genome sequence analysis approaches for L. pneumophila for environmental source attribution studies while furthering the understanding of Legionella ecology. IMPORTANCE: We demonstrate that whole-genome sequencing helps to improve resolution of Legionella pneumophila isolated during laboratory investigations of legionellosis compared to traditional subtyping methods. These data can be important in confirming the environmental sources of legionellosis outbreaks. Moreover, we evaluated various methods to analyze genome sequence data to help resolve outbreak-related isolates.