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.

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.

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.

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.

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."

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.

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.

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 Resolution of Outbreak-Associated Legionella pneumophila Serogroup 1 Isolates from New York State

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.

Chlamydia psittaci comparative genomics reveals intraspecies variations in the putative outer membrane and type III secretion system genes

Chlamydia psittaci is an obligate intracellular bacterium that can cause significant disease among a broad range of hosts. In humans, this organism may cause psittacosis, a respiratory disease that can spread to involve multiple organs, and in rare untreated cases may be fatal. There are ten known genotypes based on sequencing the major outer-membrane protein gene, ompA, of C. psittaci. Each genotype has overlapping host preferences and virulence characteristics. Recent studies have compared C. psittaci among other members of the Chlamydiaceae family and showed that this species frequently switches hosts and has undergone multiple genomic rearrangements. In this study, we sequenced five genomes of C. psittaci strains representing four genotypes, A, B, D and E. Due to the known association of the type III secretion system (T3SS) and polymorphic outer-membrane proteins (Pmps) with host tropism and virulence potential, we performed a comparative analysis of these elements among these five strains along with a representative genome from each of the remaining six genotypes previously sequenced. We found significant genetic variation in the Pmps and tbl3SS genes that may partially explain differences noted in C. psittaci host infection and disease.

Impact of analytic provenance in genome analysis

Background

Many computational methods are available for assembly and annotation of newly sequenced microbial genomes. However, when new genomes are reported in the literature, there is frequently very little critical analysis of choices made during the sequence assembly and gene annotation stages. These choices have a direct impact on the biologically relevant products of a genomic analysis - for instance identification of common and differentiating regions among genomes in a comparison, or identification of enriched gene functional categories in a specific strain. Here, we examine the outcomes of different assembly and analysis steps in typical workflows in a comparison among strains of Vibrio vulnificus.

Results

Using six recently sequenced strains of V. vulnificus, we demonstrate the "alternate realities" of comparative genomics, and how they depend on the choice of a robust assembly method and accurate ab initio annotation. We apply several popular assemblers for paired-end Illumina data, and three well-regarded ab initio genefinders. We demonstrate significant differences in detected gene overlap among comparative genomics workflows that depend on these two steps. The divergence between workflows, even those using widely adopted methods, is obvious both at the single genome level and when a comparison is performed. In a typical example where multiple workflows are applied to the strain V. vulnificus CECT 4606, a workflow that uses the Velvet assembler and Glimmer gene finder identifies 3275 gene features, while a workflow that uses the Velvet assembler and the RAST annotation system identifies 5011 gene features. Only 3171 genes are identical between both workflows. When we examine 9 assembly/ annotation workflow scenarios as input to a three-way genome comparison, differentiating genes and even differentially represented functional categories change significantly from scenario to scenario.

Conclusions

Inconsistencies in genomic analysis can arise depending on the choices that are made during the assembly and annotation stages. These inconsistencies can have a significant impact on the interpretation of an individual genome's content. The impact is multiplied when comparison of content and function among multiple genomes is the goal. Tracking the analysis history of the data - its analytic provenance - is critical for reproducible analysis of genome data.

Pyrosequencing-based comparative genome analysis of Vibrio vulnificus environmental isolates

Between 1996 and 2006, the US Centers for Disease Control reported that the only category of food-borne infections increasing in frequency were those caused by members of the genus Vibrio. The Gram-negative bacterium Vibrio vulnificus is a ubiquitous inhabitant of estuarine waters, and is the number one cause of seafood-related deaths in the US. Many V. vulnificus isolates have been studied, and it has been shown that two genetically distinct subtypes, distinguished by 16S rDNA and other gene polymorphisms, are associated predominantly with either environmental or clinical isolation. While local genetic differences between the subtypes have been probed, only the genomes of clinical isolates have so far been completely sequenced. In order to better understand V. vulnificus as an agent of disease and to identify the molecular components of its virulence mechanisms, we have completed whole genome shotgun sequencing of three diverse environmental genotypes using a pyrosequencing approach. V. vulnificus strain JY1305 was sequenced to a depth of 33×, and strains E64MW and JY1701 were sequenced to lesser depth, covering approximately 99.9% of each genome. We have performed a comparative analysis of these sequences against the previously published sequences of three V. vulnificus clinical isolates. We find that the genome of V. vulnificus is dynamic, with 1.27% of genes in the C-genotype genomes not found in the E- genotype genomes. We identified key genes that differentiate between the genomes of the clinical and environmental genotypes. 167 genes were found to be specifically associated with environmental genotypes and 278 genes with clinical genotypes. Genes specific to the clinical strains include components of sialic acid catabolism, mannitol fermentation, and a component of a Type IV secretory pathway VirB4, as well as several other genes with potential significance for human virulence. Genes specific to environmental strains included several that may have implications for the balance between self-preservation under stress and nutritional competence.