Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2008, p. 3030-3037, Vol. 74, No. 10
0099-2240/08/$08.00+0 doi:10.1128/AEM.02760-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

INRA, UR50, Laboratoire de Biotechnologie de l'Environnement, Avenue des Etangs, Narbonne F-11100, France,1 Bouisson Bertrand Laboratoires, Parc Euromédecine, 778 Rue de la Croix Verte, 34 196 Montpellier, France2
Received 7 December 2007/ Accepted 24 March 2008
|
|
|---|
|
|
|---|
In this study, the ecology of Legionella spp. in relation to the bacterial community was investigated in a cooling tower plant. The dynamics of the bacteria's structural diversity during the proliferation of L. pneumophila were analyzed. The dynamics of Legionella spp. were also investigated to determine the changes that arose in the Legionella population during the proliferation of the pathogenic species L. pneumophila. In fact, little is known regarding the Legionella species present in cooling tower systems or their dynamics within the same cooling tower (28). The growth requirements of Legionella, its ability to enter a viable-but-nonculturable state, the association of Legionella with protozoa, and the occurrence of Legionella in biofilms all tend to complicate its detection (7), and cultures on selective media have sometimes failed to isolate Legionella species from environmental samples (13). Cultivation-independent methods have thus been used to characterize the Legionella population in aquatic environments (12, 13, 30, 35). In the present study, PCR methods were used to quantify Legionella spp. and L. pneumophila and to analyze Legionella sp. diversity. The dynamics of Legionella spp., in relationship to those of the dominant bacteria, were followed by single-strand conformation polymorphism (SSCP) in a cooling tower plant over a 9-month period which included several weeks when L. pneumophila proliferated.
|
|
|---|
Chemical treatment.
Isothiazolone was used as a biocide. Volumes of 5 to 65 liters were applied weekly from 26 May 2005 to 28 September 2005 in order to reach concentrations of 100 to 250 ppm. On three dates (28 July, 11 August, and 2 September), the biocide and biodispersant (detergent) were added simultaneously. On these occasions, isothiazolone volumes were 65, 62, and 60 liters, respectively, while the biodispersant volume added was 20 liters on each date.
Sampling.
Water was sampled at different dates at a point located on the supply pipe common to all eight cooling towers, prior to its release at the top of the towers. One liter of water was collected, filtered through 0.45-µm-pore-size filters (Millipore, Billerica, MA), and frozen while awaiting DNA extraction.
ATP measurements.
For ATP measurements, the WaterGiene ATP test was used (Charm Sciences Inc., Lawrence, KS).
DNA extraction.
DNA was extracted directly from filters with the Aquadien extraction kit (Bio-Rad, Hercules, CA). For each sample, all analyses were carried out with the same DNA extract.
Real-time quantitative PCR.
Quantitative PCRs were performed with IQ-Check Legionella and IQ-Check L. pneumophila kits from Bio-Rad by following the stipulations of the manufacturer and using the iCycler IQ apparatus (Bio-Rad). Standard DNA curves were generated by amplification of serial 10-fold dilutions of genomic DNA of L. pneumophila ATCC 33152 in sterilized water and designed by an iCycler IQ apparatus. The L. pneumophila genome mass used was 4.3 fg of DNA. A standard DNA curve was generated for each assay. The cycle threshold (CT) corresponding to the number of cycles at which the reaction becomes exponential, was compared to the standard curve in order to calculate the number of genomic units (GU) in the DNA extract of the samples. The PCR equation was CT = –1/[log (1 + E)] x (log X0 – log X), where E is the efficiency of the PCR, X is the concentration, and X0 is the initial concentration. The standard DNA curve was validated when the slope was between –3.9 and –3.1, corresponding to 80 and 110% PCR efficiency, respectively. The intercept varied from 37 to 40. The Bio-Rad kit contained an internal control present in each amplification mixture to check the presence of inhibitory factors. In the event of PCR inhibition, the sample was further diluted and reanalyzed. Two negative controls were performed for each assay, a negative control for PCR (obtained by replacing the DNA with water) and a negative control for DNA extraction. Both controls had to be negative to validate the assay. PCR results were converted to genomic units per liter. Average values were calculated from duplicates. The detection threshold for Legionella spp. was 30 GU/well, which corresponded to 960 GU liter–1. In the event of PCR inhibition, the dilution of the sample brought the detection threshold to 1,920 GU liter–1 (3.3 log GU liter–1). The detection threshold for L. pneumophila was 640 GU liter–1 (2.8 log GU liter–1). The protocol of quantification by PCR of Legionella spp. and L. pneumophila followed the XP T90-471 standard elaborated by the T90E WG AFNOR (Association Française de Normalization) (6).
PCR-SSCP.
To analyze the overall structure of the bacterial community as a whole, the V3 region of the 16S rRNA gene was amplified with primers W49 (5'-ACGGTCCAGACTCCTACGGG-3', Escherichia coli position F331) and 5'-fluorescein phosphoramidite-labeled W104 (5'-TTACCGCGGCTGCTGGCAC-3', E. coli position R533) (14). PCR amplifications were performed with a Mastercycler thermal cycler (Eppendorf, Hamburg, Germany). The reaction mixtures contained 1x polymerase buffer, 0.2 mM deoxynucleoside triphosphates (dNTPs), 130 ng of each primer, 0.5 U of Pfu Turbo DNA polymerase (Stratagene, La Jolla, CA), 1 µl of genomic DNA, and water added to obtain a final volume of 50 µl. The PCR conditions were an initial denaturation step of 2 min at 94°C; 25 cycles of a three-stage program of 30 s at 94°C, 30 s at 61°C, and 30 s at 72°C; and a final elongation for 10 min at 72°C. The reactions were stopped by cooling the mixture to 4°C. The analysis of Legionella sp. diversity was performed after a nested PCR. First, primers specific to the 16S rRNA gene of Legionella spp. were used to amplify a 653-bp fragment including the V3 region, i.e., LEG-225 (5'-AAGATTAGCCTGCGTCCGAT-3') and LEG-858 (5'-GTCAACTTATCGCGTTTGCT-3') (27). The reaction mixtures contained 1x polymerase buffer, 0.2 mM dNTPs, 200 ng of each primer, 1 U of redTaq DNA polymerase (Sigma-Aldrich, St. Louis, MO), 1 µl of genomic DNA, and water added to obtain a final volume of 50 µl. The PCR conditions were an initial denaturation step of 90 s at 95°C; 30 cycles of 30 s at 95°C, 60 s at 64°C, and 60 s at 72°C; and a final elongation for 5 min at 72°C. Then, 1 µl of the PCR product was used to amplify the V3 16S rRNA gene bacterial region with primers W49 and W104 as described above. Amplification product sizes were confirmed by electrophoresis on a 2% (wt/vol) agarose gel.
SSCP analysis permits the separation of DNA fragments of the same size but with different compositions. One microliter of further diluted PCR products was added to 18 µl of formamide and 1 µl of internal size standard Rox 400 HD (Applied Biosystems, Foster City, CA) diluted 10 times. The sample was then denatured for 5 min at 95°C and placed directly on ice for 10 min. SSCP was performed with the ABI 3130 genetic analyzer (Applied Biosystems) equipped with four 50-cm capillary tubes filled with 5.6% conformation analysis polymer (Applied Biosystems) in corresponding buffer and 10% glycerol. The injection of DNA in capillaries required 5 kV for 3 s. Electrophoresis was carried out at 15 kV and 32°C for 30 min per sample. Raw SSCP data were exported into the easily handled csv format with the Chromagna shareware (developed by Mark J. Miller at the U.S. National Institutes of Health), and statistical analyses were performed with SAFUM (37) and the Matlab 6.5 software (MathWorks).
Identification of Legionella species and eukaryotic species by cloning and sequencing.
Two libraries of Legionella 16S rRNA genes were built, one with DNA from the 22 June sample and the other with DNA from the sample collected on 13 July. For the eukaryotic library, the DNA extract from water collected on 13 July was used. PCR mixtures contained 1x polymerase buffer, 0.2 mM dNTPs, 200 ng of each primer, 1 U of redTaq DNA polymerase, 1 µl of genomic DNA, and water added to obtain a final volume of 50 µl. A 653-bp fragment of Legionella 16S rRNA genes was amplified by PCR with primers LEG-225 and LEG-858 as described above, but with 35 cycles instead of 30. For amplification of a fragment of the eukaryotic 18S ribosomal DNA, primers W16 (5'-CTTAATTTGACTCAACACGG-3') (20) and W176 (5'-GGGCATCACAGACCTGTT-3') were used. The PCR conditions were an initial denaturation step of 2 min at 94°C; 35 cycles of 60 s at 94°C, 60 s at 51°C, and 60 s at 72°C; and a final elongation for 10 min at 72°C. All PCR products were purified with a QIAquick PCR purification kit (Qiagen, Hilden, Germany) in accordance with the manufacturer's instructions. The purified PCR products were cloned and transformed with the pCR4-TOPO plasmid and TOP10 E. coli competent cells, as indicated by the manufacturer (TOPO TA cloning kit; Invitrogen, Carlsbad, CA). Recombinant cells were selected by kanamycin resistance and ccd gene killer inactivation before cultivation at 37°C for 24 h in LB2x medium (tryptone at 20 g liter–1, yeast extract at 10 g liter–1, NaCl at 10 g liter–1). Sequences were obtained from clone culture (Millegen, Toulouse, France). The primer sequences were removed, and the presence of chimerical sequences was checked for with the CHECK-CHIMERA tool available at the Ribosomal Data Project (26) and the Pintail program (5). Sequences were compared with GenBank databases (www.ncbi.nlm.nih.gov/BLAST) by using the BLASTN program (4). They were imported and aligned into the January 2004 ARB database (25). The aligned sequences were added to the ARB tree by using the parsimony tool. A tree gathering sequences from environmental clones and from described Legionella species (613 bp) was then built by neighbor joining (29). The tree was rooted with the sequence from Coxiella burnetii (D89798). Bootstrap analyses (1,000 replicates) were used to assess the robustness of inferred monophyletic groups.
Nucleotide sequence accession numbers.
The 16S rRNA gene sequences of Legionella spp. determined in this study were deposited in the GenBank database under accession numbers EU309480 to EU309490.
|
|
|---|
Abundance of Legionella spp. and L. pneumophila and effect of chemical treatments.
Concentrations of Legionella spp. and L. pneumophila were determined by real-time quantitative PCR from April 2005 to January 2006. Results are gathered in Fig. 1. The threshold of detection by quantitative PCR was 3.3 log GU liter–1 for Legionella spp. and 2.8 log GU liter–1 for L. pneumophila. In April 2005, the concentration of Legionella spp. was high (6.4 log GU liter–1) and the plant was shut down and cleaned. At the subsequent startup of the cooling tower in May, the Legionella sp. concentration was below 3.3 log GU liter–1 and L. pneumophila was not detected. However, after only 5 days, the concentration of Legionella spp. increased to 5.1 log GU liter–1 and then remained between 4.4 and 5.6 log GU liter–1 for several weeks. L. pneumophila was detected in May and June, but its concentration remained under the quantification threshold until the end of June (phase 1). It started to rise in July, and the highest concentrations of L. pneumophila were obtained in July and August, at up to 4.7 log GU liter–1 (phase 2). In September, concentrations of L. pneumophila remained low. The plant was shut down in October, and from the end of October to January (phase 3), L. pneumophila was quantified at a high concentration only once, when the plant was restarted.
![]() View larger version (12K): [in a new window] |
FIG. 1. Concentrations of Legionella spp. ( ) and L. pneumophila ( ) in a cooling tower facility as determined by quantitative PCR. When detection of L. pneumophila occurred but at concentrations below the quantification threshold, results are represented by double arrows. When the concentration was below the detection threshold, results are represented by single arrows. The detection thresholds of Legionella spp. and L. pneumophila are indicated, respectively, by solid and dotted lines (Legionella detection threshold = 3.3 log GU liter–1, L. pneumophila detection threshold = 2.8 log GU liter–1). Stars indicate combined biodispersant and biocide treatments.
|
Dynamics of the microbial community. (i) Bacterial population.
The structure of the bacterial community was examined from May to September by using SSCP fingerprints. SSCP fingerprints were compared by principal-component analysis (PCA). Close SSCP fingerprints gather together on the PCA plots. The comparison is based on the totality of the SSCP signal. The first two dimensions of the PCA plots represent the two components which best highlight the differences between the SSCP fingerprints. Figure 2 presents two PCA plots, one for components 1 and 2 and one for components 1 and 3.
![]() View larger version (22K): [in a new window] |
FIG. 2. PCA of SSCP fingerprints obtained for Legionella spp. (a) and bacteria (b) during phase 1 ( ), phase 2 (), and phase 3 ( ). Samples in boxes were collected after chemical biodispersant and biocide treatment.
|
![]() View larger version (8K): [in a new window] |
FIG. 3. SSCP fingerprints showing different bacterial community structures (resilient species [a], high-diversity SSCP fingerprint [b], and after chemical treatment on 28 July [c]) and Legionella sp. diversity during phases 1 (d) and 2 (e).
|
Overall, the community structure of Legionella spp. varied depending on the time of year, whereas bacterial community structure was composed largely by the same group of a few dominant species and was reestablished within a few days after any disruption. Legionella spp. changed throughout the time period studied, with no reoccurrence of the initial population, whereas the same bacterial flora was observed during different periods of the year. Therefore, there could be no correlation between the dynamics of Legionella spp. and those of the dominant bacteria.
Phylogenetic positioning of Legionella species.
The Legionella spp. present before the proliferation of L. pneumophila (sample collected on 22 June during phase 1) and during the proliferation (sample collected on 13 July during phase 2) were identified after cloning and sequencing of 16S rRNA genes. The 16 sequences obtained from the first sample were coded SEC, and the 34 sequences obtained during phase 2 were coded SED. All 16S rRNA gene sequences had the greatest similarity to 16S rRNA gene sequences of Legionella spp., confirming the specificity of the primers used. The percentages of similarity with the closest match in the GenBank database, as given by the BLASTN program, ranged from 94 to 99%; low similarity percentages were obtained for both libraries. A phylogenetic tree was built with 16S rRNA gene sequences from described Legionella species and from the closest environmental clones (Fig. 4). The diversity of Legionella spp. was high before the proliferation of L. pneumophila, with phylotypes differently positioned on the phylogenetic tree and distantly related. The dominant species were close to L. fallonii (renamed after Legionella-like amoebal pathogen [LLAP] 10), L. fairfieldensis, L. lytica, and other LLAPs. Several sequences (SEC10, SEC15, and SED03) were only distantly related to described species and were associated with uncultured Legionella spp. from drinking water, groundwater, sludge, or sediment. In the sample collected in July, L. fallonii and L. pneumophila were the main species recovered, with, respectively, 47 and 33% of the sequences. Positioning on the phylogenetic tree was not directly related to positioning on the SSCP fingerprints, since SSCP analysis was based on a smaller fragment of the 16S rRNA gene (around 200 bp). However, peak A gathered L. lytica and other LLAPs, as well as L. pneumophila, and peak B could be identified as L. fallonii. The other peaks corresponded mainly to environmental clones. Peak B dominated in July and August and was present from April to September.
![]() View larger version (25K): [in a new window] |
FIG. 4. Phylogenetic tree showing positioning of Legionella spp. present in a cooling tower facility prior to L. pneumophila proliferation (phase 1, SEC clones) and during its proliferation (phase 2, SED clones). The percentage of each phylotype in the corresponding library is shown in parentheses. Cluster 1 includes L. moravica (Z49729), L. quateirensis (Z49732), L. worsleiensis (Z49739), and L. shakespearei (Z49736). Cluster 2 includes L. anisa (X73394), L. parisiensis (Z49731), L. dumoffii (Z32637), L. gormanii (Z32639), L. cherrii (X73404), L. wadsworthii (Z49738), L. steigerwaltii (Z49737), L. tucsonensis (Z32644), and L. bozemanae (Z49719). Cluster 3 includes L. cincinnatiensis (Z49721), L. santicrucis (Z49735), L. longbeachae (AY444741), and L. sainthelensi (X73399). Cluster 4 includes L. jamestowniensis (X73409), L. jordanis (X73396), L. brunensis (X73403), L. birminghamensis (Z49717), and L. quinlivanii (Z49733). Cluster 5 includes L. israelensis (Z32640), L. nautarum (Z49730), L. oakridgensis (X73397), L. impletisoli (AB233209), L. yabuuchiae (AB233212), L. busanensis (AF424887), L. gresilensis (AF122883), and L. beliardensis (AF122884).
|
|
|
|---|
During the spring and summer, Legionella sp. concentrations varied between 4.5 and 5.5 log GU liter–1. The maximal concentrations of L. pneumophila were obtained in July and August and were between 4.5 and 5 log GU liter–1. These concentrations are within the range of values previously reported for cooling towers (36). The high level of L. pneumophila obtained when the plant was restarted in October may be due to fouling of the network and considerable biofilm development during the shutdown period. This result fits in with the fact that cooling towers are implicated in outbreaks of legionellosis, particularly at startup or during construction (7, 10). It has been suggested that there may be a relationship between high Legionella counts in cooling towers and the occurrence of outbreaks of legionellosis (31). In this study, for similar Legionella sp. concentrations, some samples were highly contaminated with L. pneumophila, while in others the pathogenic bacterium was not detected. These results reinforce the existing assumption that the assessment of health risks from cooling towers cannot be reliably based upon single and infrequent Legionella tests (9, 24).
Combining biocides and biodispersants did reduce the concentration of Legionella spp. and of L. pneumophila. However, this effect was transient. Chemical treatment used to control Legionella in human-made water systems does not lead to total eradication of the bacterium, and recolonization occurs as soon as the treatment is interrupted (34). Legionellae are, in fact, protected inside amoebae and in biofilms and can proliferate again, recontaminating the water once the biocides lose their effect (21, 34). Amoebae present during the proliferation of L. pneumophila were identified after amplification of the 18S ribosomal DNA with primers universal for the domain Eukarya. Sequences close to those of Plantyamoeba and Vannella were obtained. These flagellate amoebae are frequently found in freshwater or seawater and are affiliated with Legionella (32). Acanthamoeba was also present. This amoeba is commonly isolated from Legionella-contaminated plumbing systems (32) and is an important host of L. pneumophila in water (7). The survival of L. pneumophila within Acanthamoeba cells during biocide treatment has been reported (21), and the presence of this amoeba could explain why the effect of biocides on Legionella abundance was only transient.
The treatment reduced the concentration of Legionella spp. to such a low level that it was not possible to analyze their diversity. However, when the population grew back to its initial level, the Legionella spp. were the same as prior to treatment. This shows clearly that even if the treatment induced a change in the structure of the Legionella population, its effect was only transient. Concerning dominant bacteria, a major disruption of the community structure due to biocide and biodispersant treatment was observed (samples collected on 29 July and 1 August, i.e., 1 and 4 days after treatment). The differences observed between the other SSCP bacterial fingerprints may relate to the origin of the cells. The similar fingerprints showing a low diversity may correspond to planktonic bacteria present in all samples, and the increase in bacterial diversity may originate from pieces of biofilm detached from the surfaces. This is congruent with the fact that when the biofilm was greater in size, the diversity increased (data not shown). It would be interesting to prove this in a subsequent experiment by comparing the bacterial diversity in water filtered through 10-µm filters with that obtained with 0.45-µm filters. Overall, it appears that the modification of the bacterial flora was probably due, in some cases, to circulating pieces of biofilm, resulting in a greater diversity in the SSCP fingerprints, while in other cases it resulted from chemical treatments, which led to very disrupted fingerprints (Fig. 3).
The suitability of primers LEG-225 and LEG-858 for detecting Legionella spp. in water has been reported previously (12, 13, 27, 35). The present study has shown that these primers can also be used to analyze Legionella diversity by SSCP. Concerning the phylogenetic positioning of Legionella spp., some relationships in the phylogenetic tree were supported moderately by bootstrap values, as previously observed by Carvalho et al. (13).
The high diversity of Legionella spp. in cooling towers has been clearly demonstrated (Fig. 4). This is an important result because the diversity of Legionella spp. in the environment, and especially in cooling tower systems, is poorly documented (28). Several sequences were close to those retrieved from treated surface water supplies and treated groundwater supplies in The Netherlands (clones Tag and S in Fig. 4) (35). Furthermore, LLAPs were particularly well represented in the cooling tower network studied. L. fallonii, which was constantly present in the system and dominated during the summer, was only described fairly recently (1) and has as its type strain LLAP-10T. A recent study of Dutch tap water installations also reported a large proportion of LLAPs (16). Those protozoonotic bacilli, which initially were isolated in coculture with protozoa, were named LLAPs because of their ability to infect and multiply intracellularly within amoebae in the same way that legionellae do (1, 2). Overall, a large proportion of sequences determined in the present study belong to as-yet-uncultured legionellae, a result which highlights the discrepancy between described Legionella species and the Legionella species present in the environment, including human-made systems.
L. pneumophila was one of the two dominant Legionella species during the summer (13 July). However, it was not detected by sequencing at the end of June (22 June). This rapid proliferation of L. pneumophila, in less than 3 weeks, may correspond to an increase in the "background" concentration of L. pneumophila in the environment during the summer, leading to the contamination of the system through water or air. A summer and autumn peak in incidence has been described as a characteristic epidemiological feature of Legionnaires' disease in Europe (11). Furthermore, a seasonal pattern of infection has been reported by Fliermans et al. (18), who injected guinea pigs with sample water collected monthly from a thermally altered lake, the highest frequency of infection by L. pneumophila occurring during the summer months. Such a seasonal variation may well be due to the facts that Legionella spp. multiply faster in the warmer waters of summer and that the greater use of cooling towers in summer provides opportunities for dissemination (11). During the investigation of Legionella colonization in 31 cooling towers in South Australia, it was found that between 60 and 75% of the cooling towers studied were colonized by Legionella during the summer months but only 20 to 30% were colonized during the winter (8). In the present study, the increase in L. pneumophila concentrations, concomitant with stability in the Legionella concentration, indicates that there must have been a decrease in the concentrations of other Legionella species during the same period. Furthermore, peak B dominated the SSCP fingerprints in phase 2, whereas it was nondominant in phase 1. This means that not only did L. pneumophila become one of the dominant species during phase 2, but the relative proportions of other Legionella species also changed significantly. In particular, peaks C and D disappeared from the SSCP fingerprints and the apparent diversity decreased. Factors which might have led to these important changes in the structure of the Legionella population still remain to be identified, i.e., the antagonistic relationships between certain Legionella species and L. pneumophila, differences in temperature growth ranges, changes in amoebal hosts, etc. The results obtained thus open new perspectives for future research. It would also be interesting to analyze Legionella diversity during the proliferation of L. pneumophila in other cooling towers in order to ascertain whether the dynamics observed in the present study also occur in other systems.
The objective of this study was to determine whether the different dynamics of bacteria and of Legionella spp. were correlated in cooling towers. A quantitative correlation was found in some cases (a decrease in ATP levels simultaneous with a decrease in the Legionella sp. concentration after chemical treatment), but it was not qualitative. In particular, no change in the structure of the bacterial community was observed during the proliferation of L. pneumophila. Bacterial diversity was modified by the chemical treatment or the development of biofilms but did not depend on the time period. In contrast, variations in the Legionella population structure were observed over time, with significant changes occurring at the end of June, following a stable phase from mid-May to mid-June. In a previous study of microbial communities in a wide range of aquatic samples containing L. pneumophila, no relationship was found between the occurrence of L. pneumophila and the associated microbiota (15). It was observed, however, that the occurrence of L. pneumophila was possible within a certain range of species richness and diversity. The authors concluded that the relationship between the occurrence of L. pneumophila and the bacteriological characteristics of water is complex and that it may therefore be interesting to concentrate on specific groups of microorganisms. In this study, we used a different approach, focusing our investigation on one cooling tower facility over several months, but we also found no relationship between the dynamics of the Legionella population during the proliferation of L. pneumophila and the dynamics of the dominant bacteria in the system. It may be that the factors that regulate the occurrence of Legionella are different from those that govern overall bacterial populations in cooling tower systems (23).
This study has provided new input on the diversity of Legionella spp. in cooling towers, showing that structural changes in the Legionella population are linked to the time of year, displaying stable periods but also a rapid major shift during the proliferation of L. pneumophila. This study has also demonstrated the application of molecular techniques which enabled us to conduct a comprehensive study of the dynamics of Legionella in relation to the bacterial community as a whole, gathering data on abundance and population structures which cannot easily be obtained by culture-dependent approaches. Controlling the Legionella risk in cooling towers requires exact knowledge of which constituents of the microbial community play a role in the proliferation of the pathogenic species, as well as the importance of these microorganisms compared to abiotic parameters such as the temperature, the mineral composition of the water, or the structure of the network. In this perspective, this study has demonstrated the value of studying the dynamics of the microflora rather than using one-off analyses.
Published ahead of print on 4 April 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»