*Legionella* is a waterborne bacterium that causes Community-Acquired Pneumonia (CAP) and hospital-acquired pneumonia (HAP).[@ref1] After water distribution systems, cooling towers are the leading sources of Legionnaires' disease outbreaks.[@ref2] These structures, which are used to cool warm water by heat transfer, have a temperature range (20 to 50°C) ideal for growth of *Legionella*. The water vapor released by cooling towers may contain *Legionella*, which can cause human infections. The Centers for Disease Control and Prevention (CDC) recommend periodic routine culturing of cooling tower water. The Occupational Safety and Health Administration (OSHA) recommends immediate action (namely, disinfection or pH adjustment) if *Legionella* concentrations exceed 10 Colony Forming Units (CFU)/mL. Despite the various disinfectants used, cooling tower water contamination issues persist worldwide due to inadequate maintenance and monitoring protocols.[@ref3] Multiple countries mandate *Legionella* reporting to prevent the spread of Legionellosis. The regulations and surveillance systems implemented include the International Health Regulations (IHR) (WHO 2007), European Working Group for *Legionella* Infections (EWGLI, <http://www.ewgli.org>), and European Legionnaires' Disease Surveillance Network (ELDSNet) (<http://ecdc.europa.eu>).[@ref4] In developing countries, the incidence of healthcare-associated infections (HAI) is 2- to 20-fold higher than that in developed countries due to the lack of surveillance systems.[@ref5] Determining the optimal *Legionella* monitoring protocol and methods is required to improve our understanding of the behavior and responses to external factors, such as disinfectants, of *Legionella* and its relationship with other microorganisms.[@ref3],[@ref6] *Legionella* spp. surveillance is conducted primarily by conventional culturing methods. However, these have limited sensitivity, a long incubation period, and the possibility of false negative results due to bacteria's being in a viable but not cultivable (VBNC) state.[@ref7] In contrast, nucleic acid-based methods and metagenomics have had a major impact on environmental and clinical microbiology and have in many cases replaced traditional culture dependent methods.[@ref8] Comparison of conventional culture and molecular methods would enhance our ability to rapidly detect *Legionella*. We thus compared standard conventional culture, real-time polymerase chain reaction (RT-PCR), and 16S rRNA next generation sequencing for *Legionella* surveillance and detection. Moreover, the bacterial community composition of cooling tower water was evaluated.

Methods {#sec1-2}
=======

Water sampling {#sec2-1}
--------------

This is an environmental research study where a total of 12 water samples were collected at 4 time points from December 2015 to May 2016, from a cooling tower at King Fahd Hospital of the University (KFHU) (Imam Abdulrahman Bin Faisal University). At each time point, three water samples were taken using sterile 1-L glass bottles containing sodium thiosulfate, and transferred directly to the microbiology laboratory of the Institute for Research and Medical Consultations (IRMC) for microbiological testing and molecular analyses.

Cultivation of Legionella {#sec2-2}
-------------------------

Cooling tower water samples were subject to *Legionella* isolation and enumeration on glycine vancomycin pyrophosphate cycloheximide (GVPC) medium according to International Standard protocol ISO 11731 (ISO, 1998).[@ref9] Initially, neat and 1:10 diluted (1/10) samples were spread onto agar plates. Next, the water samples were concentrated by filtering 500-mL samples through 0.22-µm polycarbonate filters (Millipore). The bacteria were resuspended in 5 mL of sterile water by sonicating for 4 min (Grant, Cambridge, UK), and 0.2 mL of the neat and 1:10-diluted concentrates were spread on GVPC agar. Next, the concentrates were subjected to heat (50°C for 30 min) and acid (0.2 M HCl--KCl for 5 min) treatments and then spread on GPVC agar.

For each water sample, 6 GVPC plates were streaked and incubated at 37°C for 7 days. Colonies were counted after 3 and 7 days. Colonies suspected to be *Legionella* were incubated on BCYE agar with and without L-cysteine at 37°C for 3 days. Isolates that grew only on BCYE with L-cysteine were confirmed as *Legionella*. Results are expressed in CFU/L.

DNA extraction {#sec2-3}
--------------

To extract DNA, 500 mL of each water sample were filtered through 0.22-µm polyvinylidene difluoride (PVDF) filters (Millipore), and the filters were crushed using sterile pestles on dry ice. The debris was immersed in tissue lysis buffer and subjected to DNA extraction following the instructions included in the High Purity PCR Template Preparation Kit (Roche).

Legionella quantification by real-time PCR {#sec2-4}
------------------------------------------

*Legionella* was quantified by real-time PCR amplification of the 16S rRNA gene (386 bp) using a Rotor-Gene RG3000 thermocycler (Corbett) and *Legionella* mericon kit (Qiagen). The PCR reaction mixtures (20 µL) were prepared according to the manufacturer's protocol under the following amplification conditions: 95°C for 5 min, then 40 cycles of 95°C for 15 seconds, 60°C for 15 seconds and 72°C for 10 seconds. Four standards, together with triplicate positive (included in the kit) and negative (DNA replaced with PCR-grade water) controls, were included in each PCR reaction. Rotor-Gene software (version 6) was used to analyze the data. Results were expressed as the mean genome units (GU)/L of the 3 samples obtained per time point.

Sequencing of the 16S rRNA gene {#sec2-5}
-------------------------------

Two cooling tower water samples were analyzed by 16S rRNA next-generation sequencing. Two libraries were prepared according to the 16S metagenomic sequencing library preparation guide. A 450-bp fragment of the V3--V4 region of the 16S rRNA was amplified using Illumina barcoded primers (F515/R806; 5'-GTGCCAGCMGCCGCGGTAA-3' /5'-GGACTACHVGGGTWTCTAAT-3') on the Illumina Miseq platform. Indices adapted to the barcodes were added using the Nextera XT DNA Sample Preparation Index Kit (Illumina). Following amplification, the samples were purified using the AMPure XP Beads Kit (Beckman Coulter). After quantification by real-time PCR using the KAPA Library Quantification Kit (Applied Biosystems), the libraries were subjected to cluster generation and 250-bp paired-end sequencing on Miseq using the Miseq Reagent Nano Kit, version 2 (Illumina).

Sequence analyses {#sec2-6}
-----------------

The fastq files generated were combined and their quality checked using FastQC software (\>20). The sequences were next denoised using MOTHUR v. 1.36.0 software.[@ref10] Low-quality reads (\<150 bp) and chimeras were discarded. A distance matrix was created using a cutoff of 0.2, and sequences were classified into operational taxonomic units (OTU) at a 97% similarity level. All sequences representing less than 0.005% of the total sequence reads were discarded to minimize the sequencing error rate. Each representative sequence was aligned to the SILVA bacteria database using BLASTn. Sequences identities of ≥95% and ≥97% were used as criteria for genus and species, respectively.

Accession numbers {#sec2-7}
-----------------

The sequences were deposited in GenBank under the accession numbers KY077116--KY077141 for W1 and KY348802--KY348828 for W2.

Results {#sec1-3}
=======

Legionella enumeration {#sec2-8}
----------------------

Culturing on GVPC medium resulted in non-detection of *Legionella* spp. in all of the water samples. In contrast, *Legionella* concentrations in the water samples determined by real-time PCR ranged from 5.11 × 10^5^ to 3.09 × 10^6^ GU/L (**[Figure 1](#F1){ref-type="fig"}**). Three water samples showed very close results (5.11x10^5^ GU/l, 5.72x10^5^ GU/l and 8.74x10^5^ GU/l). These concentrations were slightly different from the forth sample (3.09x10^6^ GU/l).

![Legionella concentrations (mean ± 2SD) in cooling tower water by real-time PCR.](SaudiMedJ-39-137-g001){#F1}

Bacterial community composition {#sec2-9}
-------------------------------

The bacterial community composition of 2 water samples from cooling tower drift was determined. Total of 10,238 high-quality sequence was obtained from the W1 library and 11,164 high-quality sequence was obtained from the W2 library. The majority of these sequences were assigned to Proteobacteria (W1, 63%; W2, 60%), Bacteroidetes (W1, 16%; W2, 25%), and Gemmatimonadetes (W1, 11%; W2, 4%) (**[Figure 2](#F2){ref-type="fig"}**).

![Bacterial community composition of water libraries W1 and W2.](SaudiMedJ-39-137-g002){#F2}

The following Proteobacteria subclasses were detected in the W1 and W2 libraries: Alphaproteobacteria (W1, 31%; W2, 34%), Betaproteobacteria (W1, 15%; W2, 13%), Gammaproteobacteria (W1, 10%; W2, 4%), and Deltaproteobacteria (W1, 5%; W2, 7%). Bacteriodetes included in majority Sphingobacteria (W1, 23%; W2, 23%) (**[Figure 3](#F3){ref-type="fig"}**).

![Bacterial phyla detected in water libraries W1 and W2.](SaudiMedJ-39-137-g003){#F3}

The sequences in library W1 was assigned to 14 genera and the sequences in library W2 was assigned to 15 genera (**[Figure 4](#F4){ref-type="fig"}**). Unclassified bacteria comprised 13.5% of library W1 and 14.6% of library W2. *Hydrogenophaga* predominated in both libraries (relative abundance, 12%). The rest of the genera (12) were shared in majority between the 2 water samples where *Sediminibacterium, sphingopyxis, Algoriphagus, Legionella, Lacibacter* and *Vampirovibrio* exceeded 6% of relative abundance.

![Bacterial genera detected in water libraries W1 and W2.](SaudiMedJ-39-137-g004){#F4}

Discussion {#sec1-4}
==========

Cultivable Legionella spp {#sec2-10}
-------------------------

No *Legionella* spp. were detected by culturing. This may be due to: (1) sample preparation (concentration, sonication, and acid/heat treatment, (2) the VBNC state, or (3) the presence of bacteria that inhibit growth of *Legionella* on GVPC medium.[@ref7],[@ref11] According to the OSHA guidelines, a negative culture result is justification for not taking any action. In contrast, detection of 1,000--10,000 CFU/L requires immediate disinfection and review of control measures according to the EWGLI (2011).[@ref12]

Total Legionella spp {#sec2-11}
--------------------

Real-time PCR indicated the presence of 10^5^ to 10^6^ GU/L *Legionella spp*. in all water samples. These concentrations include both viable and non-viable *Legionella*, which hampers assessment of the impact on health (**[Figure 1](#F1){ref-type="fig"}**). These results are in agreement with previous reports of high *Legionella* concentrations in natural water storage and rivers, as well as in constructed hot water systems and cooling towers.[@ref13] Moreover, the results of culture and quantitative PCR tend to show discrepancies, possibly due to protection of *Legionella spp*. within protozoa-like amoebae (detected by PCR) and their inability to grow on culture media.[@ref14],[@ref15] These bacteria use protozoa as a shield against disinfectants and environmental stresses such as low nutrient levels; such organisms are termed amoebae-resisting bacteria (ARB).[@ref3] Few studies of *Legionella* monitoring in water systems have considered amoebae.[@ref14],[@ref16] Future work should assess the involvement of free-living amoebae (FLA) with *Legionella spp*. in cooling tower water.

Bacterial community composition {#sec2-12}
-------------------------------

The bacterial community composition of 2 water samples from the cooling tower drift was evaluated. Total of 10,238 high-quality sequence was obtained from library W1 and 11,164 high-quality sequence was obtained from library W2; *Proteobacteria* predominated in both. This phylum typically constitutes the major portion of the bacterial flora of cooling towers.[@ref17]

Of the detected genera, 78% (13.5 ≤ relative abundance ≥ 3.4%) were common to both libraries. The predominant bacterium was the autotrophic hydrogen-oxidizing *Hydrogenophaga*.[@ref18] This bacterium is naturally present in soil and water and frequently detected in cooling towers, indicating high metabolic activity in cooling tower water.[@ref17] *Sphingopyxis* was also detected at high abundance (W1, 9%; W2, 9%). *Sphingopyxis* is a gram-negative bacterium frequently detected in water systems; it belongs to the family *Sphingomonadaceae*, which exhibits high level antimicrobial resistance and is a pioneer in biofilm formation.[@ref19],[@ref20] *Sphingopyxis* is also a possible pathogen in drinking-water distribution systems.[@ref21]

Soil bacteria including *Sediminibacterium* (W1, 11.2%; W2, 9%), *Algoriphagus* (W1, 9%; W2, 7.9%), and *Lacibacter* (W1, 7.9%; W2, 7.9%) were also detected.[@ref22] The predominant genera were similar among the water samples, indicating high reproducibility, despite detection of unique bacteria in each water sample (relative abundance ≤3.4%) (W1, *Limnobacter* (3.4%) and *Verrucomicrobia* (1.1%); W2, *Methylibium* (2.2%), *Aquicella* (3.4%), and *Porphyrobacter* (1.1%)).

Limitation of the study {#sec2-13}
-----------------------

*Legionella* spp. were detected in both of the water libraries (W1, 9%; W2, 5.6%). The 450-bp 16S rRNA technique used in this study was not capable of identifying *Legionella* to the species level. The genus *Legionella* includes pathogenic species such as *L. pneumophila*, the causative agent of Legionnaires' disease, which can be spread via aerosols through cooling tower drift.[@ref23] Therefore, to identify bacteria to the species level, it is necessary to sequence the full length of the 16S rRNA gene.

In conclusion, in this study, 2 molecular methods, but not a conventional culture method, successfully detected *Legionella* in cooling tower water. Until scientifically validated as primary methods of detection, nucleic acid-based techniques should be applied in conjunction with culturing for prescreening and monitoring of *Legionella* in water samples. Future work should also involve monitoring of free-living amoebae in cooling tower water.

*The authors would like to thank the technical staff at the Institute for Research and Medical Consultations (Mrs. Charmeen Vilanueva Helaga) for technical support*.

**Disclosure.** This study was funded by Imam Abdulraman Bin Faisal University, Dammam, Kingdom of Saudi Arabia (Grant number UD2014-009-MED)
