Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2005, p. 663-671, Vol. 71, No. 2
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.2.663-671.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Departamento de Zoologia and Centro de Neurociências e Biologia Celular,1 Departamento de Bioquímica and Centro de Neurociências e Biologia Celular, Universidade de Coimbra, Coimbra, Portugal2
Received 14 June 2004/ Accepted 2 September 2004
|
|
|---|
|
|
|---|
Legionella spp. are ubiquitous in aqueous environments, where they parasitize and multiply in protozoa. Several amoebae and some ciliated protozoa are known to be potential environmental hosts for legionellae (1). This host-parasite interaction has been extensively studied and is considered to be the fundamental link to the ecology and the pathogenesis of these bacteria. It is now clear that the association of legionellae with protozoa is a major factor regulating the presence of the bacteria in the environment. Indeed, protozoa provide a necessary environment for the growth of legionellae, and they enhance the resistance of these organisms to adverse environmental conditions (2).
Forty-nine species of Legionella with >70 serogroups have been classified to date (26), including five species designated Legionella-like amoebic pathogens that have been isolated only from protozoa (3, 26). Despite the large species diversity,
90% of disease cases are caused by Legionella pneumophila, especially strains belonging to serogroup 1, while the remainder are generally caused by strains of other serogroups of L. pneumophila, L. micdadei, and L. longbeachae (52). However, L. longbeachae constitutes 31% of the isolates involved in community-acquired disease in Australia and New Zealand (52).
The infection caused by the inhalation of legionellae depends on the ability of these organisms to enter and to multiply within alveolar macrophages, causing the destruction of these phagocytes and damage to the pulmonary tissues (21). Some virulence traits, such as high cytopathogenic activity, have been recognized in several L. pneumophila strains that are more virulent than others, but a set of distinctive characteristics correlated with increased pathogenicity could not be clearly defined (4). The Dot/Icm complex, however, is critical for the virulence capability of Legionella. The dot/icm genes, initially identified in L. pneumophila, are also present in all other Legionella spp. studied (4). This complex is thought to constitute a type IV-like secretion system that is capable of injecting effector molecules into the host cell, which allows Legionella to evade the endocytic pathways, probably by delaying the phagosome-lysosome fusion in macrophages, among other roles in the pathogenicity to macrophages (11).
It is not completely clear why specific legionellae are the major causes of disease; some of these organisms could have greater inherent virulence, or they could simply be more abundant in environments that promote their dispersal to humans. Understandably, the vast majority of the studies of the distribution of legionellae have been performed in man-made aquatic environments, since they are generally implicated in the dispersal of legionellae to humans. Legionella species, primarily L. pneumophila, have been frequently sought in air-conditioning cooling towers, potable-water distribution systems, and associated fixtures, where they are found in high numbers and can be dispersed to humans (5, 12, 33, 45).
In these environments, water temperature and other physical and chemical parameters influence colonization (18, 50). Moreover, the association among legionellae, amoebae, and other protozoa generally found in biofilms is known to be very important in promoting their growth (6, 13, 47, 48). The colonization of surface water and other natural environments, such as hydrothermal areas and thermally altered lakes and streams, by Legionella spp. has been less investigated (8, 17, 32, 44). The presence and persistence of Legionella spp. in groundwater is poorly documented, and the role of this environment as a potential natural reservoir of legionellae has not yet been investigated (16). However, the increased use of this kind of water by large human populations concerns public health authorities, since a large number of pathogens are known to contaminate groundwater (31).
We determined the presence of Legionella spp. and identified them for several years in water samples collected from different boreholes in two distinct areas venting hydrothermal water. We also evaluated the persistence of Legionella strains in these groundwater environments by typing the isolates using random amplified polymorphic DNA (RAPD)-PCR. Finally, sequence analysis of a housekeeping gene (rpoB) and a virulence-related gene (dotA) was used to define the genetic organization of the L. pneumophila populations isolated and their relationship with type and reference strains of L. pneumophila.
|
|
|---|
![]() View larger version (26K): [in a new window] |
FIG. 1. Geological formations in the vicinity of the sampling sites in central Portugal. (A) Sampling sites in area 1. (B) Sampling sites in area 2.
|
The water temperature in wells 2A and 2B was 48°C, while in boreholes 2C and 2D it was 45.5°C, with a pH of 7.9.
Sampling and isolation of Legionella spp.
A total of 66 groundwater samples from area 1 and 45 samples from area 2 were collected at irregular intervals over a 7-year period. To evaluate the influence of the plumbing system and the hydraulic pumps present in one borehole on colonization of the groundwater by Legionella, we removed the inner stainless steel pipe that brings water to the surface, the hydraulic pump, and the electrical cable of well 1A. On that occasion, the components were swabbed to collect biofilm samples. Afterwards, the water in the borehole, without the entire plumbing system, was chlorinated to 50 ppm, which was maintained for 48 h by continuously pumping chlorine through a sterile silicone tube that reached the bottom of the well. The pH of the water in the borehole was adjusted to 7.0 by pumping in an HCl solution to increase the disinfecting effect of the chlorine. Water samples were collected before, during, and after this procedure. We also collected water samples at different depths with a disinfected silicone tube connected to a peristaltic pump.
All water samples were recovered in sterile 2-liter containers, transported at ambient temperature, concentrated, and plated within 24 h. The water was filtered through 45-mm-diameter membrane filters (Gelman Supor 200; 0.2-mm pore size). The filters were placed in small sterile plastic bags containing 10 ml of the original water. The bags containing the filters were rubbed manually for a few minutes to remove organisms from the filters. One portion of the concentrated sample (0.1 ml) was spread directly on the surface of buffered charcoal yeast extract (BCYE) medium containing glycine, vancomycin, polymyxin B, and cycloheximide, designated GVPC; other portions of the concentrated sample were subjected to acid and heat treatments before being spread on GVPC (14). The plates were incubated in a normal atmosphere at 37°C for up to 9 days.
Presumptive Legionella sp. colonies were counted and subcultured on BCYE and BCYE without cysteine. All isolates that grew only on BCYE alone were considered putative Legionella strains, assigned separate designations, and maintained at 80°C in 5% (wt/vol) yeast extract with 15% (vol/vol) glycerol.
Identification of Legionella spp.
Identification of Legionella isolates was performed by analysis of the fatty acid methyl ester (FAME) profiles. Cells were harvested, and FAMEs were obtained by saponification, methylation, and extraction, as described previously by Kuykendall et al. (25); separated; identified; and quantified with the MIS library generation software (Microbial ID Inc., Newark, Del.). The identification was made by comparison of the FAME profiles with a previously constructed database that allows the identification of the vast majority of Legionella spp. (15).
Random amplified polymorphic DNA typing.
Several isolates were selected for RAPD-PCR analysis, taking into consideration the species, the site, and the date of sampling. The RAPD-PCR profiles were also compared with those obtained from isolates recovered 12 years before from a spring in area 2 by Marrão et al. (32). We also determined the RAPD profiles of type strains of each species found in the groundwater samples.
Crude cell lysates were used as DNA templates for RAPD typing as described elsewhere (51). Amplification reactions were performed with a total volume of 50 µl and standardized as follows: 1.5 U of Taq polymerase, 1.5 mM MgCl2 (Pharmacia Biotech), 0.2 mM (each) deoxynucleoside triphosphate, 0.6 µM primer OPA3 (5'-AGTCAGCCAC-3'), and 2.0 µl of crude cell lysates. Samples were subjected to 45 cycles of amplification (Perkin-Elmer model 240) as follows: 1 min at 94°C, 1 min at 34°C, and 2 min at 72°C, followed by a final extension step of 7 min at 72°C (40). The fragments were analyzed by electrophoresis in a 2% agarose gel in Tris-acetate-EDTA buffer.
Genetic structure of L. pneumophila groundwater populations.
L. pneumophila isolates from areas 1 and 2 with different RAPD profiles were selected for sequencing of variable regions of the rpoB and dotA genes to infer the genetic structures of these populations. The same regions of L. pneumophila subsp. pascullei strains U8W (ATCC 33737T), U7W (ATCC 33736), and MICU B (ATCC 33735) were also sequenced. The extraction of genomic DNA was carried out as described by Rainey et al. (36). The amplification of the genes and sequencing of the purified PCR products were carried out as described by Ko et al. (23). The purified reaction mixtures were electrophoresed on a model 310 Genetic Analyzer (Applied Biosystems, Foster City, Calif.).
The phylogenetic analysis was performed after manual checking of the quality of the sequences using the Bioedit editor (19) and alignment against sequences of 21 type and reference strains of L. pneumophila obtained from the public databases (Table 1), using the multiple-alignment CLUSTAL X software package (43). The method of Jukes and Cantor (22) was used to calculate evolutionary distances; phylogenetic dendrograms were constructed using the neighbor-joining method (39), and tree topologies were evaluated by performing bootstrap analysis of 1,000 data sets using the MEGA2 package (24). The amino acid sequences were deduced with the MEGA2 package from the 300- and 360-bp DNA sequences for the partial rpoB and dotA gene sequences, respectively.
|
View this table: [in a new window] |
TABLE 1. Accession numbers of rpoB and dotA genes of L. pneumophila type and reference strains used
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Enumeration, identification, and RAPD types of Legionella spp. isolated from groundwater samples collected from borehole 1A in area 1
|
|
View this table: [in a new window] |
TABLE 3. Enumeration, identification, and RAPD types of Legionella spp. isolated during borehole 1A dismantling and disinfection
|
|
View this table: [in a new window] |
TABLE 4. Enumeration, identification, and RAPD types of Legionella spp. isolated from groundwater samples collected in area 2
|
RAPD typing of Legionella isolates.
All 72 L. oakridgensis isolates typed from area 1 had only one RAPD type, designated LO1 (Fig. 2). Only one RAPD type, designated LS1, was obtained from 29 L. sainthelensi isolates typed from this area (Fig. 2). L. pneumophila isolates were less frequently detected in area 1 but revealed higher RAPD pattern diversity. Indeed, from the seven L. pneumophila isolates typed, we obtained three distinct RAPD profiles, namely, LP1, LP2, and LP3 (Fig. 2).
![]() View larger version (62K): [in a new window] |
FIG. 2. Agarose gel electrophoresis of RAPD types obtained from representative isolates of Legionella spp. detected in groundwater and biofilm samples from areas 1 and 2 and respective type strains. M, Marker XIV (Roche Molecular Biochemical). Values to the left are in base pairs.
|
Representative L. pneumophila isolates (32 isolates) recovered 12 years before from area 2 possessed RAPD types LP4 and LP5, like the more recent isolates, but RAPD types LP8 and LP9 were obtained only from L. pneumophila strains isolated during the earlier study. RAPD types LP6 and LP7 were present in the population only during the previous 7 years (Fig. 3). One of the L. londiniensis profiles (RAPD LL1) also had a long persistence, while RAPD LL3 was only encountered 12 years before, after which it was no longer detected (Fig. 3). As expected, the RAPD profiles exhibited by the type strains of these species were distinct from those obtained from our isolates (Fig. 2).
![]() View larger version (83K): [in a new window] |
FIG. 3. Comparison of RAPD types obtained from strains isolated in 1991 and the more persistent RAPD types from 2003 in area 2. M, Marker XIV (Roche Molecular Biochemical). Values to the left are in base pairs.
|
![]() View larger version (33K): [in a new window] |
FIG. 4. Phylogenetic dendrograms based on analyses of partial sequences of dotA (A) and rpoB (B) of the groundwater isolates and of type and references strains of L. pneumophila. The trees were created using the neighbor-joining method. Isolates and the respective RAPD types characterized in this study are indicated in boldface. Branches supported by values of >50% in the bootstrap analysis (1,000 replications) are indicated. Scale bar, 5 inferred nucleotide substitutions per 100 nucleotides.
|
The deduced amino acid sequences from the partial rpoB gene sequences of all isolates and reference strains were the same, despite the nucleotide differences detected. On the other hand, a remarkable difference in the amino acid sequences was deduced from the dotA nucleotide sequences (results not shown). The degree of amino acid diversity was different within each of the three principal dotA gene clusters. The mean diversity index calculated using the Mega2 equation was: 0.026 ± 0.010 for strains included in cluster dotA-A, 0.025 ± 0.009 for strains included in cluster dotA-B, and 0.045 ± 0.013 for strains included in cluster dotA-C.
|
|
|---|
The ability of Legionella spp. to persistently colonize aquifers was fully demonstrated by the consistent recovery of these organisms, over several years, in both areas studied. In area 1, nevertheless, we never detected legionellae in the samples collected from borehole 1B. The inability to detect legionellae in this borehole, which is only 62 m from 1A, where legionellae were always detected, could indicate the existence of localized environments within the aquifer, generating different conditions that would promote the development of Legionella or at least that could sustain these bacteria in a cultivable state. It should be noted that boreholes 1A and 1B have been drilled to different depths (63.5 m versus 307 m, respectively) and that one is only stainless steel jacketed to a depth of 38 m. These differences could influence colonization by and persistence of the organisms or their host protozoans. Our results could also indicate that the circulation of water through the aquifer is limited, leading to localized colonization of the aquifer by the microorganisms. Moreover, legionellae could not be recovered during the first year of sampling from the most recently drilled wells in area 2 (2C and 2D), indicating that colonization by the organisms is slow.
We enumerated legionellae from the water and from the plumbing structures by dismantling borehole 1A, followed by disinfecting the water column, in an attempt to determine the origin of the colonization. The results clearly indicate that legionellae were present in the groundwater, since the organisms were recovered from the water samples collected at different depths without the pumping system only when the chlorine became undetectable, indicating that groundwater had entered the borehole and reached the surface by natural pressure from the aquifer.
The numbers of Legionella spp. detected in groundwater were lower than those generally found in water samples collected from the man-made environments (5, 33, 49), perhaps because of smaller amounts of biofilms found in groundwater compared with the sometimes massive development of biofilms in surface waters and man-made aquatic environments (13, 32).
L. pneumophila represented the major legionella population persistently recovered from area 2, while L. oakridgensis was the most abundant species isolated from area 1. Remarkably, there are very few environments where L. pneumophila is not the predominant Legionella isolate. Indeed, potting soils in Australia are primarily colonized by L. longbeachae, and it has been suggested that this would be its natural habitat (42). The absence of L. pneumophila was also observed in another Portuguese hydrothermal area, where we isolated only L. sainthelensi and another unidentified Legionella sp. (15, 44).
RAPD analyses demonstrate the long persistence of some clones in groundwater. Distinct clones of L. pneumophila are known to be persistent, and it has been hypothesized that some of these clones are better suited to persist in man-made environments (10, 27, 35). We recognized, for example, two L. pneumophila clones and one L. londiniensis clone in the groundwater for 12 years in area 2, as well as the persistence of one clone each of L. oakridgensis and L. sainthelensi over 7 years in area 1. These results support the idea that the interaction between legionellae, the presence of particular protozoa, and physicochemical parameters may lead to the colonization and persistence of specific clones or species in aquatic environments.
The diversity of RAPD profiles was higher in area 2 than in area 1, which may be related to undetermined environmental parameters and/or to characteristics of the L. pneumophila clones that constitute the major population. Indeed, the genetic structure of these L. pneumophila populations, based on the dotA and rpoB sequence analyses, indicates that the majority of the isolates are closely related to a strain designated 1169-MN-H (7). On the other hand, Ko et al. (23), who used the same genes to analyze the structures of populations from clinical cases and man-made environments, found that the vast majority of the organisms clustered with other reference strains of the three subspecies of L. pneumophila.
Recently, the dotA gene, together with icmX, icmW, and icmV, were found to be the most variable genes in the entire dot/icm complex (34). Interestingly, the deduced amino acid sequences of the dotA gene from the isolates included in the dotA-A cluster (Fig. 4A) have a degree of diversity slightly higher than that calculated for the cluster (dotA-B) that included the majority of the reference strains belonging to L. pneumophila subsp. pneumophila. This result led us to speculate that the diversity of at least the dotA allele is higher in strains isolated from groundwater environments than in those recovered from man-made environments. These results appear to confirm the hypothesis that L. pneumophila isolates from clinical cases and man-made environments belong to a restricted subset of all clones in the species as a whole, as anticipated by Selander et al. (41).
The topologies of inferred rpoB and dotA trees are not similar, as previously observed (23), since the relationship between the L. pneumophila subspecies differs in the two trees. While the inferred rpoB tree agrees with phylogenetic and taxonomic analyses based on other methods (9, 20, 46), the inferred dotA tree does not, since the three subspecies cannot be separated from each other by the sequence analysis of this gene. The main discrepancies between the two trees may result from dotA gene exchange involving L. pneumophila (sensu lacto) strains. Gene exchange can occur by two distinct routes: homologous recombination between closely related individuals or lateral gene transfer between organisms belonging to more distant lineages (28). We may envision lateral gene transfer of dotA within all known reference strains of L. pneumophila subsp. fraseri and L. pneumophila subsp. pascullei and some strains of L. pneumophila subsp. pneumophila. Moreover, the fact that different clones (determined by RAPD analysis) have the same dotA or rpoB sequence genes may indicate the occurrence of homologous recombination events. These events may mediate the dispersal of advantageous alleles that may rise to high frequencies among genetically related individuals by means of sporadic selection events and may be an explanation for differences in the persistences of certain strains in a given environment.
The groundwaters investigated in this study are used, without any disinfection procedure, as water supplies in spas for therapeutic purposes, so that persistence of Legionella spp. is certainly a cause for concern for the spa officials and public health authorities. During the period of our investigations, no case of Legionella-related disease was detected in either area. The absence of disease may be related to the low number of legionellae detected in these groundwaters and the very stringent cleaning procedures adopted within both facilities that prevent proliferation in the water distribution systems. Lastly, the uniqueness of the majority of the dotA alleles detected in these L. pneumophila isolates could also lead us to hypothesize that the strains more frequently present in these particular environments are less pathogenic than others more common in man-made environments. This appealing hypothesis requires further investigation.
We thank Isabel Pinhal for helping to isolate legionellae and Fernanda Nobre for FAME analysis.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»