Legionella spp. are the cause of a severe bacterial pneumonia known as Legionnaires' disease (LD). In some cases, current genetic subtyping methods cannot resolve LD outbreaks caused by common, potentially endemic L. pneumophila (Lp) sequence types (ST), which complicates laboratory investigations and environmental source attribution. In the United States (US), ST1 is the most prevalent clinical and environmental Lp sequence type. In order to characterize the ST1 population, we sequenced 289 outbreak and non-outbreak associated clinical and environmental ST1 and ST1-variant Lp strains from the US and, together with international isolate sequences, explored their genetic and geographic diversity. The ST1 population was highly conserved at the nucleotide level; 98% of core nucleotide positions were invariant and environmental isolates unassociated with human disease (n = 99) contained ~65% more nucleotide diversity compared to clinical-sporadic (n = 139) or outbreak-associated (n = 28) ST1 subgroups. The accessory pangenome of environmental isolates was also ~30-60% larger than other subgroups and was enriched for transposition and conjugative transfer-associated elements. Up to ~10% of US ST1 genetic variation could be explained by geographic origin, but considerable genetic conservation existed among strains isolated from geographically distant states and from different decades. These findings provide new insight into the ST1 population structure and establish a foundation for interpreting genetic relationships among ST1 strains; these data may also inform future analyses for improved outbreak investigations.

Cooling towers (CTs) are a leading source of outbreaks of Legionnaires' disease (LD), a severe form of pneumonia caused by inhalation of aerosols containing Legionella bacteria. Accordingly, proper maintenance of CTs is vital for the prevention of LD. The aim of this study was to determine the distribution of Legionella in a subset of regionally diverse US CTs and characterize the associated microbial communities. Between July and September of 2016, we obtained aliquots from water samples collected for routine Legionella testing from 196 CTs located in eight of the nine continental US climate regions. After screening for Legionella by PCR, positive samples were cultured and the resulting Legionella isolates were further characterized. Overall, 84% (164) were PCR-positive, including samples from every region studied. Of the PCR-positive samples, Legionella spp were isolated from 47% (78), L. pneumophila was isolated from 32% (53), and L. pneumophila serogroup 1 (Lp1) was isolated from 24% (40). Overall, 144 unique Legionella isolates were identified; 53% (76) of these were Legionella pneumophila. Of the 76 L. pneumophila isolates, 51% (39) were Lp1. Legionella were isolated from CTs in seven of the eight US regions examined. 16S rRNA amplicon sequencing was used to compare the bacterial communities of CT waters with and without detectable Legionella as well as the microbiomes of waters from different climate regions. Interestingly, the microbial communities were homogenous across climate regions. When a subset of seven CTs sampled in April and July were compared, there was no association with changes in corresponding CT microbiomes over time in the samples that became culture-positive for Legionella. Legionella species and Lp1 were detected frequently among the samples examined in this first large-scale study of Legionella in US CTs. Our findings highlight that, under the right conditions, there is the potential for CT-related LD outbreaks to occur throughout the US.

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.