Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention20202420332200008-168910.3201/eid1603.081689ResearchResearchNew Endemic Legionella pneumophila Serogroup I Clones, Ontario, CanadaL. pneumophila Serogroup I Clones, CanadaTijetNathalieTangPatrickRomilowychMyaDuncanCarlaNgVictoriaFismanDavid N.JamiesonFrancesLowDonald E.GuyardCyrilOntario Agency for Health Protection and Promotion, Toronto, Ontario, Canada (N. Tijet, P. Tang, M. Romilowych, C. Duncan, D.N. Fisman, F. Jamieson, D.E. Low, C. Guyard)Research Institute of the Hospital for Sick Children, Toronto (V. Ng, D.N. Fisman)Australian National University, Canberra, Australian Capital Territory, Australia (V. Ng)University of Toronto, Toronto (F. Jamieson, D.E. Low, C. Guyard)Mount Sinai Hospital, Toronto (F. Jamieson, D.E. Low, C. Guyard)Address for correspondence: Cyril Guyard, Ontario Agency for Health Protection and Promotion, 81 Resources Rd, Toronto, Ontario, M6H 2S2, Canada; email: cyril.guyard@ontario.ca32010163447454

Identifying geographic distribution can improve surveillance and clinical testing procedures.

The water-borne pathogen Legionella pneumophila serogroup 1 (Lp1) is the most commonly reported etiologic agent of legionellosis. To examine the genetic diversity, the long-term epidemiology, and the molecular evolution of Lp1 clinical isolates, we conducted sequence-based typing on a collection of clinical isolates representing 3 decades of culture-confirmed legionellosis in Ontario, Canada. Analysis showed that the population of Lp1 in Ontario is highly diverse and combines lineages identified worldwide with local strains. Identical types were identified in sporadic and outbreak-associated strains. In the past 15 years, the incidence of some lineages distributed worldwide has tended to decrease, and local endemic clones and lineages have emerged. Comparative geographic distribution analysis suggests that some lineages are specific to eastern North America. These findings have general clinical implications for the study of Lp1 molecular evolution and for the identification of Lp1 circulating strains in North America.

Keywords: Legionella pneumophila serogroup 1molecular epidemiologysequence-based typingphylogenyclonal analysisemerging clonesbacteriaCanadaresearch

Legionella species are implicated in 2 clinical syndromes: Legionnaires’ disease (LD) and Pontiac fever, which are collectively known as legionellosis. Pontiac fever is a self-limited, influenza-like illness, whereas Legionnaires’ disease is a common cause of serious bacterial pneumonia (1,2).

Among the 52 species and 70 serogroups of Legionella species (3), L. pneumophila is the major cause of sporadic and outbreak legionellosis (91.5%), and serogroup 1 is the predominant serotype (84.2%) (4). In industrialized countries, L. pneumophila is the second most common pathogen detected in cases of community-acquired pneumonia that requires patient admission to intensive care units (5,6).

During an outbreak or after the detection of sporadic cases, appropriate identification and typing methods are essential for epidemiologic investigations. Adequate typing methods are also crucial to determine the degree of relatedness of bacteria and to enable the reconstruction of microevolutionary events (7). On the basis of analysis of 7 loci, a standard sequence-based method for the typing of L. pneumophila serogroup 1 (Lp1) was developed by the European Working Group for Legionella Infections (EWGLI) (8,9).

In previous population-based studies, a Legionella sequence-based typing (SBT) scheme was used to analyze clinical strains either from Europe or with limited time-span coverage (1012). In the present study, we applied SBT to examine the genetic diversity, the long-term epidemiology, and the molecular evolution of Lp1 clinical isolates using a population-based collection that encompassed isolates from 30 years of culture-confirmed legionellosis cases in Ontario.

MethodsSource of Isolates

Legionellosis is a notifiable disease in Ontario (population 13 million persons). Since 1978, the diagnosis of Legionella infections has been centralized at the Ontario Public Health Laboratory (OPHL). This laboratory serves as the Legionella reference laboratory and performs all testing for outbreak investigations and most testing of clinical specimens. Therefore, isolates analyzed in this study are representative of the strains isolated in Ontario in the past 3 decades. Information available in the Ontario database includes dates of onset of illness, patient’s age and sex, and city and hospital or healthcare facility from which specimens were submitted (13). No specimens were submitted for Legionella isolation from 1978 through 1979. From 1980 through 1985, a mean ± SD of 424.1 ± 281.3 specimens was submitted for isolation every year. From 1986 through 2007, a mean ± SD of 1,783.5 ± 258.4 specimens was submitted for isolation every year. The mean ± SD number of Lp1 culture-confirmed cases/year during the study period was 7.4 ± 3.5. The proportion of culture-confirmed case-patients with L. pneumophila infection remained stable during the period of analysis, and 66% of the isolates were Lp1 (13).

Lung tissues, bronchial-alveolar lavage specimens, or sputum specimens were homogenized by using a tissue grinder, streaked on buffered charcoal yeast extract agar plates and incubated at 37°C (3–7 days). Species and serogroups were confirmed by direct immunofluorescent antibody assay and slide-agglutination (14,15). Isolates (n = 217) were stored at –80°C in trypticase soy broth supplemented with 5% horse blood. Twenty-three isolates obtained in 1996 and 1997 could not be used and were omitted from our analysis. Outbreaks were defined as >2 cases that were submitted from the same hospital or healthcare facility or with links to a common source with onset during the same 30-day period.

Sequence-based Typing

SBT using loci flaA, pilE, asd, mip, mompS, proA, and neuA was conducted according to the EWGLI scheme (8,9). Automated contig-assembly and base-calling of DNA sequence traces were performed by using the EWGLI sequence quality tool (16). The sequences obtained from this work are available in the EWGLI-SBT database (www.ewgli.org)

Phylogenetic and Allelic Diversity Analyses

Multiple sequence alignments of concatenated DNA sequences and phylogram construction were carried out with ClustalW2 (www.ebi.ac.uk/Tools/clustalw2/index.html) by using the neighbor-joining method with 1,000 bootstrap replicates (17). Clonal analyses were performed by using eBURST3 (http://eBURST.mlst.net) with a group definition set to 6 identical alleles and sequence type (ST) allelic profiles were clustered with the unweighted pair group method with arithmetic mean (UPGMA) algorithm by using splits Tree4 (18). The standardized index of association (IAS) and the mean genetic diversity were calculated with LIAN 3.5 (19).

Ontario Health Regions and Rates Calculations

The 36 public health units of Ontario were aggregated into 7 health regions (OHRs) with populations ranging from ≈0.5 to 2 million persons: Toronto, South West, Central South, Central West, Central East, East, and North. Average rates were calculated by dividing disease counts by the Statistics Canada population estimates (20,21). OHR population estimates were not available before 1995 and were estimated from 1990 through 1995 and 2006 through 2007 by linear extrapolation.

Statistical Analyses and Mapping

Statistical analyses were performed by using STATA (StataCorp, College Station, TX, USA). Thematic maps were created by using ArcGIS (ESRI, Redlands, CA, USA).

ResultsLp1 Sequence-based Typing

The 194 isolates, collected from 1978 through 2007, were resolved into 62 STs (Table A1). Seven STs were represented by at least 10 isolates, 13 STs consisted of groups containing 2–4 isolates, and 42 STs were represented by 1 isolate. In comparison to the EWGLI dataset, 42 STs have only been reported in North America, and 41 STs are unique to the province of Ontario. The ST with the largest number of isolates was ST1 (n = 31). STs previously reported in the EWGLI database and responsible for >9 cases in Ontario include ST36 (n = 10), ST37 (n = 21), ST42 (n = 10), and ST62 (n = 16). Two STs, specific to the province of Ontario, were detected in >9 legionellosis cases: ST211 (n = 15) and ST222 (n = 13). STs resolved from outbreak isolates were confirmed to be epidemiologically concordant since related isolates were assigned identical STs. ST211 strains were obtained from patients in 1 outbreak (n = 2) in 1993 and from 13 patients with sporadic cases. Seven of the 13 ST222 isolates were recovered from a legionellosis outbreak at a long-term care facility in Ontario in 2005 (10). ST226 was also differentiated from strains responsible for a suspected outbreak (n = 2) and is specific to Ontario.

Across the 7 loci, 99 alleles were identified. Three new alleles were found, 2 of which (asd 32 and proA 33) were identified in a ST357 strain isolated in 2002 from a patient with a sporadic case. The third new allele (mompS 52) was identified in a ST358 strain isolated in the South West OHR. At the individual loci level, the total number of alleles ranged from 10 at flaA to 21 at mompS. Because the population of Lp1 clinical isolates found in Ontario appeared to be distinct from the isolates reported in the EWGLI database, we performed linkage analyses and looked at genetic diversity of our dataset compared to that of the EWGLI dataset. Linkage analyses showed that the IAS for the complete dataset was 0.4913. This value is comparable with the IAS of 0.494 in previous studies and suggests linkage disequilibrium in the dataset obtained from the population of clinical Lp1 isolates in Ontario (22). With a value of 0.8041 ± 0.0155, the mean ± SD genetic diversity of our dataset was similar to the mean genetic diversity of the EWGLI dataset.

Phylogenetic Analysis

The population structure of Lp1 clinical isolates from Ontario was analyzed by using 62 concatenated sequences of 7 loci and compared with results of cluster analysis deduced from SBT allelic profiles. Three major clusters were visually identified from the phylogenetic tree (Figure 1). All clusters contain isolates from outbreak and sporadic cases. None of the identified clusters or subgroups were exclusively formed with strains identified with STs specifically reported in Ontario. This suggests that Ontario strains are phylogenetically related to strains found in the EWGLI dataset. Phylogenetic cluster I (n = 39) included the epidemic strain ST1 and 7 STs of sporadic cases. With 114 isolates and 46 STs, cluster II was the largest and most diverse group from the dataset. In this cluster, Ontario outbreak strains ST37 and ST211 were subgrouped with ST36 (Philadelphia strain). Cluster III comprised 8 STs and 41 isolates. With the exception of ST222 (n = 13), none of the STs grouped in this cluster were reported to be outbreak strains.

Phylogenetic analysis of flaA, pilE, asd, mip, mompS, proA, and neuA concatenated sequences from the 62 Legionella pneumophila serogroup 1 sequence types (STs) identified in Ontario. The tree was constructed with ClustalW2 (www.ebi.ac.uk/Tools/clustalw2/index.html) and the neighbor-joining method with 1,000 bootstrap replicates. Scale bar indicates genetic distances between sequences. STs in boldface were detected in outbreaks.

Next, the UPGMA algorithm was used to construct a dendrogram based on a matrix of pairwise allelic differences between the 62 STs of our dataset (Figure 2). The topology of the UPGMA dendrogram was partially congruent with the neighbor-joining tree based on allelic sequences. The UPGMA dendrogram contains 3 major clusters of related STs arbitrarily named A, B, and C (Figure 2). Cluster B contains all isolates of cluster II except ST210 and ST199, which grouped with cluster C. In contrast, STs found in clusters I and III were separated into clusters A and C. ST1 and ST52 clustered in a separate branch at the base of the dendrogram, (Figure 2), which suggests that they could be phylogenetically distant from other STs. However, this divergence was not observed with the neighbor-joining method. Based on this finding, for the rest of the analysis, we considered cluster II as a well-defined phylogenetic group and clusters I/III were analyzed as a single group.

Dendrogram created by the unweighted pair group method with arithmetic mean method based on the 62 allelic profiles of 194 Legionella pneumophila serogroup 1 isolates. Clonal groups (CGs) identified by eBURST (http://eBURST.mlst.net) are indicated with solid lines, and STs included in CGs are in boldface. Ontario STs included in clonal complexes (CC) identified by comparative eBURST analysis with the European Working Group for Legionella Infections database are indicated with dashed lines. The 3 major clusters are indicated on the right of the figure with bold lines. The number of strains isolated in Ontario is indicated below CG. Scale bar indicates linkage distances.

Identification of Clonal Lineages

The eBURST clonal analysis of our strains showed that the province of Ontario presents a semiclonal population with 27 single isolates and 11 clonal groups (CGs) (Figure 2). With 54 isolates and 10 STs, CG1 was the clonal group with the largest number of isolates and STs. This clonal group (27.8% of Ontario isolates) contained STs that were reported elsewhere (ST36, ST37, and ST104) but also STs that were unique to Ontario (ST193, ST195, ST196, ST197, ST211, and ST229). The founder of CG1 was predicted to be ST36 (bootstrap confidence [BC] = 68%), and the predominant single locus variant of this group was ST37 (n = 21). Members of CG1 were recovered from both sporadic and outbreak cases. CG2 (n = 7) only contained isolates with STs that are unique to Ontario and the ancestor of this group was predicted to be ST209 (BC = 28%). All isolates of CG2 were obtained from sporadic cases. Other clonal groups unique to Ontario included CG4 (n = 17), CG6 (n = 4) and CG7 (n = 2). Each of the 5 remaining clonal groups contained only 2 STs with combinations of STs specific to Ontario or previously reported in the EWGLI database.

We next did an eBURST comparative analysis of the SBT dataset from Ontario with the EWGLI database. STs detected in Ontario were only clustered in 17 of the 44 clonal complexes (CC) identified with the EWGLI database. This suggests that >60% of the clonal groups of the EWGLI dataset are absent from Ontario. Comparative eBURST analysis showed that CG1 is part of CC-B, which is the most diversified clonal complex in the international database (59 STs) (Figure 3). In addition to STs included in CG1, Ontario isolates ST40 and ST202 grouped in CC-B. Philadelphia strain ST36 was predicted to be the founder (BC = 91%) of this clonal complex. The high number of STs clustered in CC-B (12 STs) suggests that strains belonging to this group are evolving in Ontario (Figure 2).

Representation of Legionella pneumophila serogroup 1 clonal complexes (CC) A, B, C, D obtained by comparative eBURST (http://eBURST.mlst.net) analysis between the Ontario collection and the European Working Group for Legionella Infections database. Each circle represents a single sequence type (ST). Size of the circle is proportional to the number of isolates. Dark circles represent predicted founder of each CC. Labels in boldface indicate STs found in both datasets, regular black characters indicate STs absent from the Ontario collection, light gray characters indicate STs exclusively found in Ontario. Solid lines represent single-locus variants.

Isolates from CG11 grouped within a CC comprising the highest number of EWGLI isolates (n = 490). The predicted founder of this clonal complex is ST1 and despite its high number of isolates, it comprises only 35 STs. Similarly, Ontario CG11 (n = 33) contained only 2 STs which suggests that strains belonging to these clonal groups have limited genetic variability.

Isolates from CG4 (n = 17) grouped with 2 North American clinical strains ST276 and ST289 in CC-A (Figure 3). Another clonal complex of interest was CC-C accounting for 16 STs (Figure 3). This complex included isolates grouped in CG2 and CG8 as well as 7 STs that were not identified in Ontario. Single isolates ST45, ST109, and ST230 from the Ontario dataset also grouped in CC-C. Finally, CC-D of the comparative eBURST analysis was identical to CG6 (Figure 2).

Temporal Trends of Lp1 Culture-confirmed Legionellosis in Ontario

During 1978–2007, differences could be observed in the distribution of specific STs, clonal complexes and phylogroups. During 1981–1994, ST1 strains were regularly isolated (n = 29) with case numbers ranging from 0 to 5 (peaking in 1983) (Figure 4, panel A). During this time period, ST1 caused 3 outbreaks and a significant increase in ST1 occurrence was observed (incidence rate ratio [IRR] 15.37, 95% confidence interval [CI] 3.67–64.43, p<0.001). In contrast, after 1994, prevalence of ST1 isolates decreased markedly (IRR 0.07, CI 0.02–0.27, p<0.001). On average, from 1978 through 2007, we observed a significant decrease of 9% per year of ST1 strains in Ontario (p<0.001). From 1995 through 2007, only 2 legionellosis ST1 culture-confirmed cases were reported.

Prevalence of Legionella pneumophila serogroup 1 sequence type 1 (ST1) (A) and ST211 (B) endemic strains in Ontario. Black bar sections indicate proportion of strains from isolated cases and white bar sections indicate proportion of isolates from outbreaks.

In contrast, other STs have emerged in the past 20 years. ST47 was detected 3 times during 2003–2006. This strain was not isolated before 2003 in Ontario and it corresponds to the ST of the Lorraine strain. This emerging strain is highly prevalent in France where it was reported as the cause of 2 major outbreaks (11). Two other emerging strains that are unique to Ontario are ST211 and ST222 (Figure 4, panel B). ST211 was first isolated in 1989 accounting for 12.5% of clinical isolates. It was regularly isolated from 1990 through 2006, and sporadic cases peaked in 1999 (23.1%). ST222 was first reported in 1999, and the prevalence of this strain has significantly increased (IRR 1.30 per year, CI 1.12–1.53, p<0.001). Excluding outbreak isolates, ST222 accounted for 11.1% to 15.4% of clinical isolates in 1999, 2000, 2006, and 2007.

At the clonal complex level, some groups of strains have recently emerged in the province of Ontario. Ontario strains of CC-A were not detected in Ontario before 1992 and oscillated from 10% to 25% from 1998 through 2007 (Figure 5, panel A). This observation is consistent with the emergence of ST222, which is a major contributor of this clonal complex. Chronological evolution of Ontario isolates belonging to CC-C is also noteworthy because it appears to be an emerging clonal complex (Figure 5, panel B). From 1985 through 1988, the maximum incidence of isolates from this clonal complex was 15.4%. CC-C isolates were not reported in 1989 and 1995, but the incidence of CC-C isolates gradually increased from 1998 until 2007 (when it peaked at 37.5%).

Incidence of Ontario Legionella pneumophila serogroup 1 isolates from clonal complexes (CC) A and C. CC-A (A) and CC-C (B) were identified by eBURST (http://eBURST.mlst.net) comparative analysis using the Ontario and the European Working Group for Legionella Infections international databases. Black bar sections indicate proportion of strains isolated during sporadic cases. White bar sections indicate proportion of outbreak isolates.

Geographic Distribution of Lp1 Clinical Strains Isolated in Ontario

The geographic distribution of some individual STs, clonal complexes and phylogenetic clusters was not homogenous. With geographic ratios between 37.5% and 100%, the 7 STs that caused >10 clinical cases were all prevalent in the Toronto OHR. ST1 was widely distributed: South West, Central South, Central East, East, and Toronto OHRs. ST37 was found in all OHRs except the East OHR. ST62 was also homogeneously distributed in all OHRs with the exception of the North OHR. In contrast, distributions of ST36 and ST42 were not homogenous because they were only identified in the Central South, Central West, East, and Toronto OHRs. Distribution of ST211 was limited to the Toronto OHR. From 1978 through 2006, ST211 was identified 8 times in the same hospital. Excluding suspected linked cases, ST211 was exclusively detected 7 times in the Toronto OHR. Despite their recent emergence in 1999, ST222 strains were reported in multiple OHRs, including South West, Central West, Toronto, and East.

Comparative eBURST analysis between Ontario and EWGLI datasets grouped Ontario single isolates ST208 (n = 2) with ST257. ST257 was identified in a clinical case of legionellosis from New Hampshire, USA. This small clonal complex may be geographically restricted to eastern North America. Similarly, geographic distribution of CC-A (CG4, ST276, and ST289) was restricted to eastern North America. ST276 and ST289 were only reported in the states of New York and Connecticut (23). With the exception of the North OHR, strains belonging to CC-B were reported in all OHRs, although a high prevalence of CC-B isolates were identified in the Toronto OHR (63.2%). Without significant geographic prevalence, CC-C isolates were identified in all OHRs. CC-D comprised STs only reported in Ontario (Figure 3).

Geographic distributions of Ontario major phylogenetic groups were analyzed by mapping average rates of culture-confirmed Lp1 cases according to OHRs from 1990 through 2007 (Figure 6). Rates of clusters based on sequence-based types and phylogroups appeared to be partially dependent on geographic location. As expected from our distribution analyses of individual STs and clonal complexes, the Toronto OHR was more likely to have legionellosis cases caused by cluster I/III (0.062/100,000 persons/year) and cluster II (0.08/100,000 persons/year) than all other OHRs. With rates ranging from 0.02/100,000 person years for the North OHR to 0.08/100,000 person-years for the Toronto OHR, isolates from cluster II were unevenly reported in all OHRs. In contrast, legionellosis caused by cluster I/III were not identified in the North OHR and the rate for the Central OHR was only 0.01/100,000 person-years, which is 3.3× less than the rate reported for cluster II.

Geographic distribution of phylogenetic clusters II and I/III from 1990 through 2007. Rates are cases of infection with Legionella pneumophilia serogroup 1 clones per 100,000 persons per year. The province of Ontario was divided into 7 health regions (OHRs) with populations ranging from ≈0.5 to 2 million persons: Toronto, South West (SW), Central South (CS), Central West (CW), Central East (CE), East, and North.

Discussion

This report represents the first large-scale population-based SBT analysis of Lp1 clinical isolates within North America over a 30-year period. Sixty-two STs were identified among the isolates of the Ontario collection, which reflects a high degree of genetic diversity of Lp1 clinical isolates. Forty-one STs were unique to Ontario. Thus, the population of clinical Lp1 of this province consists of a combination of widely distributed and local isolates.

Although most sporadic cases were caused by isolates with a unique ST, some STs, like ST1, were found to be responsible for sporadic cases and outbreaks cases. This finding is in agreement with the recent identification of the Paris strain in sporadic and outbreak cases (24). Moreover, we have identified additional clones, specifically ST37, ST211, ST222, and ST226, which were also detected in sporadic and outbreaks cases. Although some small clonal groups such as CG2 exclusively comprised sporadic cases, our comparative analysis found no specific correlation between clonal complexes or phylogenetic clusters and ability to cause sporadic or outbreak cases.

A notable finding of our study is the high proportion of ST1 isolates identified, which supports the hypothesis that some specific Lp1 clones have gained widespread dissemination. This hypothesis was proposed after analyses of protein polymorphism and pulsed-field gel electrophoresis showed similar patterns from isolates distributed worldwide (25,26). Recently, the Paris strain was suggested to be one of these worldwide distributed strains because it has been identified in many European countries in patients and environmental samples (24,27). In France, the Paris strain was identified in 12.2% of culture-confirmed cases from 1998 through 2002, and it was shown to be the most prevalent endemic strain (28). Our study suggests that the notion of worldwide distributed strains could be broadened to include other sequence types such as ST36 or ST37. These 2 STs have been reported in European countries for clinical and environmental isolates and comprised 16% of STs from culture-confirmed cases in Ontario.

In previous studies, the clinical predominance and large distribution of ST1 suggested that it is a stable clone, well adapted to environmental survival or to host infection (28,29). Surprisingly, although ST1 was identified in 16.5% of the culture-confirmed cases of legionellosis over the past 30 years, our study also shows that the incidence of ST1 strains has decreased dramatically during the past 12 years. Because clonal analysis suggests that ST1 presents a limited genetic variability in our geographic area, we can hypothesize that its ability to colonize the environment or to be isolated by culture or its virulence might have been impaired in the recent years. In contrast, endemic clonal groups and clinical strains like ST211 have emerged in our geographic area in the past 15 years. A surveillance study recently reported a new endemic Lorraine strain (ST47) emerged in France (11). ST47 was only recovered 3 times over 30 years in Ontario, but CC-C, comprising ST47, is an emerging clonal complex in Ontario. Our analysis suggests that, globally, ST1 strains are being replaced by other emerging strains or clonal complexes.

Geographic distribution analysis of culture-confirmed population rates suggests that strains from cluster II are largely distributed in Ontario, whereas clusters I/III were mostly reported in the OHRs in close proximity to Lake Ontario. This finding could reflect differences in ecologic niches (either combined with degree of adaption of organisms to cause human disease or not). Some endemic emerging STs and clonal groups are exclusively detected in Ontario, in eastern North America, or in both. In the United States, the census regions with highest incidence rates for legionellosis are East North central and Middle Atlantic, at the proximity of the Great Lakes (30). Legionella species are abundant in surface waters and the Great Lakes ecosystem might represent an ideal ecologic niche for these bacteria. This hypothesis is in agreement with the identification of clonal complexes comprising isolates exclusively originating from eastern North America. This finding contrasts with findings of a recent population structure analysis of L. pneumophila that used allelic profiles from the EWGLI database that could not identify eBURST groups containing profiles originating from a single geographic area (22).

In conclusion, we showed that the population of clinical Lp1 in the province of Ontario is a combination of worldwide distributed and local strains. Our population of isolates might represent more severe cases as human respiratory samples are more frequently taken from patients requiring hospitalization, but the decreased prevalence of some clones and the emergence of local group of isolates suggest that the population of Lp1 has evolved or adapted to its environment during the past 30 years. Further research is required to explain the changing incidence of these STs and to investigate the fitness of emerging strains or clonal groups. Outcomes of this research will be helpful to improve surveillance programs for legionellosis as well as to ensure adequacy of clinical testing procedures with circulating strains.

Suggested citation for this article: Tijet N, Tang P, Romilowych M, Duncan C, Ng V, David N, et al. New endemic Legionella pneumophila serogroup 1 clones, Ontario, Canada. Emerg Infect Dis [serial on the Internet]. 2010 Mar [date cited]. http://dx.doi.org/10.3201/eid1603.081689

Acknowledgments

We thank the European Working Group for Legionella Infections for assistance and for authorizing access to EWGLI database, and David J Farrell for reviewing the manuscript.

This study was supported by funding from the innovation program of Ontario Ministry of Government Services (no. 07007) and the Ontario Ministry of Health and Long term Care.

Sequence types and allelic profiles of clinical <italic>Legionella pneumophila</italic> serogoup 1 isolates from Ontario, 1978–2007
Sequence type
Allelic profile
Isolation dates
No. isolates
% Outbreak isolates
Unique to Ontario?
11,4,3,1,1,1,11981–20043119.35No
222,3,6,10,2,1,619871No
232,3,9,10,2,1,619931No
363,4,1,1,14,9,11986–200510No
373,4,1,1,14,9,111985–2007219.52No
403,6,1,14,14,9,111991–19922No
424,7,11,3,11,12,91989–200710No
444,8,11,10,10,12,219841No
455,1,22,26,6,10,1220041No
475,10,22,15,6,2,62003–20063No
521,10,3,1,1,1,119911No
597,6,17,3,13,11,111982–19924No
628,10,3,15,18,1,61986–200616No
8412,9,2,5,3,17,1519801No
9412,8,11,5,20,12,219991No
1043,10,1,1,14,9,119951No
1095,1,22,15,6,10,620071No
11512,29,2,.5,3,17,1520061No
1212,10,3,10,9,4,620051No
1434,17,11,23,5,12,1919891No
15411,14,16,16,15,13,220001No
1838,10,9,15,18,1,619801Yes
1922,10,5,5,35,5,619811Yes
1933,10,1,28,14,9,619811Yes
19412,8,11,2,11,12,219811Yes
1953,10,3,1,14,9,11984–19922Yes
1963,10,1,28,14,9,1119851Yes
1973,4,1,1,29,9,119851Yes
19811,14,3,1,15,13,919861Yes
1996,10,17,28,4,14,1119861Yes
20012,15,11,10,11,12,219861Yes
2023,13,1,28,1,9,319871Yes
2035,8,22,10,6,10,219881Yes
2041,10,3,3,9,4,919881Yes
2055,10,22,15,6,1,61985–19862Yes
20812,8,11,5,8,12,221986–19902Yes
2095,1,22,10,6,10,191986–20074Yes
2106,10,15,10,9,14,101988–20032Yes
2113,10,1,1,14,9,111989–20061513.33Yes
2122,10,5,5,35,5,91991–20073Yes
2132,19,5,10,18,1,21992–20053Yes
2142,10,9,10,10,1,619921Yes
2152,10,14,10,21,14,319931Yes
21812,9,2,21,3,17,1519931Yes
2195,10,22,10,6,5,619941Yes
2203,6,1,3,9,11,919951Yes
22111,14,16,1,15,13,919951Yes
2222,19,5,10,18,1,101999–20071353.85Yes
2235,1,22,10,6,10,219981yes
2244,8,11,16,42,12,219991Yes
2252,10,3,25,9,4,1320011Yes
22612,8,11,26,20,12,22001–2002366.66Yes
2272,19,3,10,18,1,220021Yes
2284,8,22,16,29,12,21995–20022Yes
2295,4,1,1,14,9,120021Yes
2305,1,22,30,6,10,1020021Yes
2318,10,3,15,21,12,2020071Yes
2325,1,22,23,6,10,1920071Yes
2343,4,1,1,14,9,919991Yes
33312,31,2,5,3,17,151988–19954Yes
35711,14,32,3,7,33,920021Yes
35812,3,9,13,52,1,619981Yes

Dr Tijet is a scientist at Ontario Agency for Health Protection and Promotion, Toronto, Ontario, Canada. Her recent research interest is the molecular epidemiology of infectious diseases and development of molecular diagnostic tools for detecting and characterizing pathogens.

ReferencesFraser DW, Tsai TR, Orenstein W, Parkin WE, Beecham HJ, Sharrar RG, Legionnaires’ disease: description of an epidemic of pneumonia. N Engl J Med. 1977;297:118997335244Glick TH, Gregg MB, Berman B, Mallison G, Rhodes WW Jr, Kassanoff I Pontiac fever. An epidemic of unknown etiology in a health department: I. Clinical and epidemiologic aspects. Am J Epidemiol. 1978;107:14960623097Euzeby JP List of prokaryotic names with standing in nomenclature—genus Legionella [cited 2008 Oct 29]. http://www.bacterio.cict.fr/l/legionella.htmlYu VL, Plouffe JF, Pastoris MC, Stout JE, Schousboe M, Widmer A, 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 10.1086/34108712089674Vergis EN, Akbas E, Yu VL Legionella as a cause of severe pneumonia. Semin Respir Crit Care Med. 2000;21:29530416088740Smith P, Moore M, Alexander N, Hicks L, O’Loughlin R Surveillance for travel-associated Legionnaires disease—United States, 2005–2006. MMWR Morb Mortal Wkly Rep. 2007;56:1261318059257Feil EJ, Enright MC Analyses of clonality and the evolution of bacterial pathogens. Curr Opin Microbiol. 2004;7:30813 10.1016/j.mib.2004.04.00215196500Ratzow S, Gaia V, Helbig JH, Fry NK, Luck PC 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. J Clin Microbiol. 2007;45:19658 10.1128/JCM.00261-0717409215Gaia V, Fry NK, Afshar B, Luck PC, Meugnier H, Etienne J, Consensus sequence-based scheme for epidemiological typing of clinical and environmental isolates of Legionella pneumophila. J Clin Microbiol. 2005;43:204752 10.1128/JCM.43.5.2047-2052.200515872220Gilmour MW, Bernard K, Tracz DM, Olson AB, Corbett CR, Burdz T, Molecular typing of a Legionella pneumophila outbreak in Ontario, Canada. J Med Microbiol. 2007;56:33641 10.1099/jmm.0.46738-017314363Ginevra C, Forey F, Campese C, Reyrolle M, Che D, Etienne J, Lorraine strain of Legionella pneumophila serogroup 1, France. Emerg Infect Dis. 2008;14:6735 10.3201/eid1404.07096118394295Borchardt J, Helbig JH, Luck PC Occurrence and distribution of sequence types among Legionella pneumophila strains isolated from patients in Germany: common features and differences to other regions of the world. Eur J Clin Microbiol Infect Dis. 2008;27:2936 10.1007/s10096-007-0392-317909867Ng V, Tang P, Jamieson F, Guyard C, Low DE, Fisman DN Laboratory-based evaluation of legionellosis epidemiology in Ontario, Canada, 1978 to 2006. BMC Infect Dis. 2009;9:68 10.1186/1471-2334-9-6819460152Thacker WL, Wilkinson HW, Benson RF Comparison of slide agglutination test and direct immunofluorescence assay for identification of Legionella isolates. J Clin Microbiol. 1983;18:111386358251Cherry WB, Pittman B, Harris PP, Hebert GA, Thomason BM, Thacker L, Detection of Legionnaires disease bacteria by direct immunofluorescent staining. J Clin Microbiol. 1978;8:32938359594Ewing B, Hillier L, Wendl MC, Green P Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8:175859521921Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:29478 10.1093/bioinformatics/btm40417846036Huson DH, Bryant D Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:25467 10.1093/molbev/msj03016221896Hudson RR Analytical results concerning linkage disequilibrium in models with genetic transformation and conjugation. J Evol Biol. 1994;7:53548 10.1046/j.1420-9101.1994.7050535.xStatsCanada Population and dwelling counts, for Canada, Provinces and Territories, and census divisions, 2006 and 2001 censuses [cited 2008 Nov 5]. http://www12.statcan.ca/english/census06/data/popdwell/Table.cfm?T=702&PR=35&SR=1&S=3&O=D; 2007StatsCanada Estimates of population (2001 Census and administrative data), by age group and sex, Canada, Provinces, Territories, health regions (June 2005 boundaries) and peer groups, annual (number) [2008 Nov 11]. http://cansim2.statcan.ca/cgi-win/CNSMCGI.EXE?Lang=E&ArrayId=01095315&Array_Pick=1&Detail=1&ResultTemplate=CII/CII___&RootDir=CII/&TblDetail=1&C2SUB=HEALTH 2007Edwards MT, Fry NK, Harrison TG Clonal population structure of Legionella pneumophila inferred from allelic profiling. Microbiology. 2008;154:85264 10.1099/mic.0.2007/012336-018310031Kozak NA, Benson RF, Brown E, Alexander NT, Taylor TH Jr, Shelton BG, Distribution of lag-1 alleles and sequence-based types among Legionella pneumophila serogroup 1 clinical and environmental isolates in the United States. J Clin Microbiol. 2009;47:252535 10.1128/JCM.02410-0819553574Cazalet C, Jarraud S, Ghavi-Helm Y, Kunst F, Glaser P, Etienne J, Multigenome analysis identifies a worldwide distributed epidemic Legionella pneumophila clone that emerged within a highly diverse species. Genome Res. 2008;18:43141 10.1101/gr.722980818256241Selander RK, McKinney RM, Whittam TS, Bibb WF, Brenner DJ, Nolte FS, Genetic structure of populations of Legionella pneumophila. J Bacteriol. 1985;163:1021374030689Drenning SD, Stout JE, Joly JR, Yu VL Unexpected similarity of pulsed-field gel electrophoresis patterns of unrelated clinical isolates of Legionella pneumophila, serogroup 1. J Infect Dis. 2001;183:62832 10.1086/31852511170989Lawrence C, Reyrolle M, Dubrou S, Forey F, Decludt B, Goulvestre C, Single clonal origin of a high proportion of Legionella pneumophila serogroup 1 isolates from patients and the environment in the area of Paris, France, over a 10-year period. J Clin Microbiol. 1999;37:2652510405416Aurell H, Etienne J, Forey F, Reyrolle M, Girardo P, Farge P, Legionella pneumophila serogroup 1 strain Paris: endemic distribution throughout France. J Clin Microbiol. 2003;41:33202 10.1128/JCM.41.7.3320-3322.200312843082Amemura-Maekawa J, Kura F, Chang B, Watanabe H Legionella pneumophila serogroup 1 isolates from cooling towers in Japan form a distinct genetic cluster. Microbiol Immunol. 2005;49:10273316365527Neil K, Berkelman R Increasing incidence of legionellosis in the United States, 1990–2005: changing epidemiologic trends. Clin Infect Dis. 2008;47:5919 10.1086/59055718665818