Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2007, p. 4522-4531, Vol. 73, No. 14
0099-2240/07/$08.00+0 doi:10.1128/AEM.02387-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Civil and Environmental Engineering, Arizona State University, Tempe, Arizona 85287,1 Department of Chemical Engineering, Environmental Engineering Program, Yale University, New Haven, Connecticut 065202
Received 10 October 2006/ Accepted 13 May 2007
|
|
|---|
|
|
|---|
Direct confirmation of biosolid wind aerosolization may not be obtainable by common wastewater bioaerosol analysis, the focus of which has been coliform detection (14, 25). The large reduction of coliforms during sewage sludge stabilization, coupled with the inability of many gram-negative microorganisms to survive environmental and collection-associated stresses (8, 43), makes coliform bacteria a poor choice for tracking biosolid-derived aerosols. Indeed, the release of biosolid bioaerosols has only been documented at the source of overt aerosol-producing activities such as biosolid spreading, loading, or disk incorporation into soils (4, 31, 38). In these cases, the large aerosol mass concentrations of biosolids (ca. 500 to 1,000 µg/m3) were documented more readily by culturing spores of thermotolerant Clostridia and bacteriophages than by coliform detection (11). The potentially lower aerosol mass loading produced during wind aerosolization events suggests the need for a sensitive, culture-independent method to track these aerosols. Culture-independent analysis circumvents sampling stress and bioaerosol die-off and allows for the inclusion of highly enriched, uncultured microorganisms as indicators. Additionally, DNA-based sequence comparison provides more specific and definitive information from which source tracking can be conducted.
The central hypothesis of the following research was that DNA-based analysis can be used to specifically and sensitively confirm the unintended aerosolization of land-applied biosolids during high-wind events. Three biosolid source tracking methods based on the 16S rRNA gene sequences of Clostridia, Chloroflexi, and Euryarchaeota ("Clostridia, Chloroflexi, and Euryarchaeota methods" hereafter) were tested and then applied to downwind and upwind control aerosol samples collected during high- and low-wind events at class B biosolid land application sites.
|
|
|---|
Three aerosol measurement stands were located at the downwind edge of the land application area and were spaced 20 m apart. Each sample stand contained two liquid impinger biosamplers (SKC, Inc., Eighty Four, PA) and one aerosol filter sampler. All samplers were suspended 1.5 m above ground and pointed directly into the wind. For upwind control, two sampling stands, each containing two liquid impingers and one filter sampler, were located approximately 150 m upwind of the biosolid application area. One real-time monitor for particulate matter smaller than 10 µm in aerodynamic diameter (PM10) (DustTrak aerosol monitor; TSI, Inc., St. Paul, MN) was located on a downwind sample stand, while another was located with the upwind samplers. Meteorological conditions were monitored every minute with a Weather monitor II (Davis Instrument Corporation, Hayward, CA); these conditions included temperature, relative humidity, wind speed, and wind direction.
Aerosol, biosolid, and soil sample collection.
Liquid impinger biosamplers were filled with 20 ml of sterile phosphate-buffered saline (PBS; 30 mM phosphate buffer, 125 mM NaCl [pH 7.2]) and collected aerosols at 12.5 liters/min. Impingers were refilled with sterile PBS during experimentation (120-min duration) to maintain volume. Impingers were also shielded from sunlight with aluminum foil. Impinger contents were used in both culture analysis and PCR-based analyses. Microbial aerosol sampling for PCR-based analyses was also carried out using filter samplers that contained 1.0-µm-pore-size polycarbonate track-etched (PCTE) filters (GE Osmonics, Minnetonka, MN) mounted on open face holders and operated at a flow rate of 33 liters/min. For long-term sample storage at 80°C, impinger contents were transferred into sterile glass tubes or filters were transferred into a sterile 5-cm-diameter petri dish and sealed with Teflon tape. Flow rates for all sampling equipment were calibrated prior to each field experiment.
For bulk biosolid and soil analysis, approximately 50 g (wet weight) of material was collected at the time and location of aerosol sampling. Each sample was a composite, containing 10 g to 20 g of bulk material from at least three different field locations. Soil samples were taken in fields adjacent to where biosolids had been applied. For long-term storage, bulk biosolids and soil samples were archived at 80°C. Moisture content of all bulk samples was determined gravimetrically by oven drying 10 g of samples at 105°C for 24 h. Soil particle size distribution was determined via sieve analysis for the largest particles and hydrometer analysis for particles smaller than 75 µm, in accordance with previously described methods (2).
MST methods.
Microbial source tracking (MST) targets that were unique to and highly enriched in biosolids were first identified by constructing phylogenetic clone libraries for bulk biosolid and aerosol samples taken directly downwind of a biosolid spreading operation, and ambient aerosol samples taken upwind from biosolid spreading operations. (Further details on these aerosol samples are available in reference 28.) Based on the results of these libraries and a review of the literature concerning anaerobic digester microbial ecology, three MST methods were chosen. These methods include (i) the culture of sulfite-reducing Clostridia followed by PCR amplification of colonies to obtain 16S rRNA gene sequences; (ii) DNA extraction, PCR amplification, and sequencing of the 16S rRNA gene for a bacterium of the class Chloroflexi (hereafter "Chloroflexi sp.") commonly found in bulk biosolids; and, (iii) DNA extraction, nested PCR, and terminal restriction fragment length polymorphism (tRFLP) analysis of 16S rRNA gene sequences of Euryarchaeota specific to biosolids. These methods were first applied to a variety of bulk biosolid and bulk soil samples to determine their average rate of correct classification and then applied to aerosol samples taken at land application sites during high-wind events.
Culture analysis.
The enumeration of sulfite-reducing Clostridia was performed using a modified membrane filtration technique (35), where cells filtered onto 0.22-µm Durapore membranes were anaerobically incubated at 37°C for 48 h on an antibiotic-supplemented egg yolk-free tryptose-sulfite-cycloserine (TSC) agar. For bulk biosolid analysis, microorganisms were extracted from 10 g (wet weight) of biosolids by being mixed with 100 ml of 0.25x Ringer's solution (38 mM NaCl, 1.4 mM KCl, 1.1 mM CaCl2, 0.6 mM NaHCO3) according to previously described methods (26). For aerosol samples, up to 35 ml of aerosol impinger collection fluid (pooled liquid impinger contents from the same sampler stand) was filtered to provide a detection level of approximately 1 CFU/m3.
DNA isolation.
DNA was extracted from frozen biosolid, soil, and aerosol samples. For aerosol filter samples, PCTE filters were eluted in the 5-ml petri dish used for storage by adding 3 ml sterile PBS with 0.05% Tween and a Teflon bar (4 cm long, 3-mm radius). This mixture was agitated on an orbital shaker at 36 rpm for 12 h. The eluant was concentrated by filtering through a 1-cm2 section of a 0.22-µm-pore-size PCTE filter. For aerosol liquid impinger samples, 5 to 10 ml of sample was also concentrated by filtering through a 1-cm2 section of a PCTE filter. After aerosol sample concentration, the 1-cm2 PCTE filters were placed in a 2-ml bead-beating tube and cells were extracted from the filters by vortexing at top speed for 2 min in 200 µl of extraction buffer I (150 mM Na2EDTA, 225 mM NaCl, 0.05% Tween 20 [pH 8.5]). DNA was then isolated by bead beating to lyse cells, phenol-chloroform extraction to separate DNA, and column filtration for DNA cleanup and concentration. This method followed a previously described method (30), with the exception that no lysozyme was added and the bead-beating rate was set at 2,500 rpm. The method described above was also used to isolate DNA from bulk soil and biosolid samples. In these cases, 200 mg of material was placed in bead-beating tubes with 200 µl of extraction buffer I.
PCR amplification.
16S rRNA genes were PCR amplified from genomic DNA extracted from bulk biosolid, bulk soil, and aerosol samples. All primers used in this study are listed in Table 1. PCR for the bulk biosolid and aerosol universal clone libraries was carried out with both bacterial (27F and 907R) and universal (515F and 1392R) primer sets. In addition to the two primer sets, both bulk biosolid and aerosol libraries were constructed from two sets of biosolid and aerosol samples collected on different days. PCR amplification was carried out in a 50-µl-volume reaction mixture containing 5 µl of 10x buffer (50 mM KCl, 10 mM Tris-HCl), 10 µl of 25 mM MgCl2, 1 µl of deoxynucleoside triphosphate (dNTP) mix (0.2 mM each dATP, dCTP, dGTP, and dTTP), 5 µl of 10-mg/ml bovine serum albumin (BSA), 4 µl of 5 µM each of the bacterial primers 27F and 907R or the universal primers 515F and 1392R, 2 µl of extracted DNA template, and 5 U of Taq DNA polymerase (Eppendorf, Westbury, NY). The following thermocycler (Eppendorf Mastercycler; Perkin-Elmer, Inc., Boston, MA) temperature profile was used (for both bacterial and universal primers): hot start for 12 min at 94°C followed by 20 cycles of 92°C for 30 s, 65°C minus 1°C/cycle for 30 s, and 72°C for 90 s. This initial amplification was followed by 20 additional cycles at 92°C for 30 s, 45°C for 90 s, and 72°C for 90 s and a final extension at 72°C for 20 min. PCR for colonies of sulfite-reducing Clostridia was performed using universal primers 515F and 1392R in accordance with reagents and the temperature profile described above. DNA template was supplied by picking single black colonies of Clostridia from plates and suspending cells in the 50-µl PCR mix. For the Chloroflexi sp., 16S rRNA gene amplification was performed with forward primer CH608 and reverse primer CH1124. PCR was carried out with 5 µl of extracted aerosol, bulk soil, or bulk biosolid DNA template in a 50-µl PCR mixture containing 5 µl of 10x PCR buffer (50 mM KCl, 10 mM Tris-HCl), 3 µl of 1.5 mM MgCl2, 1 µl dNTP (0.2 mM each dATP, dCTP, dGTP, and dTTP), 5 µl of 10-mg/ml BSA, 4 µl of 5 µM each primer, and 2.5 U Taq DNA polymerase (Eppendorf). The following temperature program was used: hot start of 5 min at 80°C, after which the Taq DNA polymerase was added; followed by 5 min at 94°C; followed by 34 cycles of 60 s at 94°C, 90 s at 56°C, and 90 s at 72°C. A final extension at 72°C for 6 min concluded the amplification. Finally, amplification of 16S rRNA genes Euryarchaeota were performed using nested PCR. The primary primer set was primer A21F specific for the domain Archaea and the universal primer 1392R. The secondary primer set included MB301F (6-carboxyfluorescein labeled at the 5' end) which circumscribes the genera Methanobacterium, Methanobrevibacter, and Methanosphaera (41) and primer AR912R specific for the domain Archaea. The primary amplification PCR mixture contained 5 µl of DNA template in 50 µl PCR mix containing 10 µl of 5x PCR buffer (Phusion HF buffer), 1 µl dNTP (0.2 mM each dATP, dCTP, dGTP, and dTTP), 1.5 µl of dimethyl sulfoxide, 2.5 µl of 5 µM each primer, and 1 U Phusion DNA polymerase (Finnzymes, Espoo, Finland). PCR for the primary amplification was performed using the following temperature program: 30 s at 98°C, followed by 34 cycles of 10 s at 98°C, 30 s at 59°C, and 30 s at 72°C. A final extension at 72°C for 7 min completed the program. The secondary PCR was carried out with 2 µl of the primary product in 50 µl PCR mix. PCR was performed as described for the primary PCR, with the exception of an annealing temperature of 66°C. In all PCRs, amplicon lengths shown in Table 1 were confirmed by gel electrophoresis. Sterile filtered water-negative controls were included in all PCRs, and no amplification was observed in these controls.
|
View this table: [in a new window] |
TABLE 1. PCR primers used in this study
|
All retrieved sequences were compared with available 16S rRNA gene sequences by using the Basic Local Alignment Search Tool (BLAST) and the National Center of Biotechnology Information (NCBI) database (1). Sequences were aligned using CLUSTALX (40), and construction of phylogenetic trees was performed by the neighbor-joining method using MEGA2 software (21).
tRFLP analysis.
A 2.5-µl master mix containing 0.3 µl restriction enzyme solution (5,000 U/ml SfcI; New England Bolas, Ipswich, MA), 2 µl of 10x NE buffer (20 mM Tris-acetate, 50 mM potassium acetate, 10 mM magnesium acetate, 1 mM dithiothreitol), and 0.2 µl BSA (100 µg/ml) was added to 17.5 µl of sterile water containing 0.3 µg PCR of product. Samples were incubated at 25°C for 16 h and then inactivated by heating to 65°C for 20 min. Digested PCR product was purified using the QIAquick nucleotide removal kit (QIAGEN, Inc., Valencia, CA) according to the manufacturer's protocol and eluted in 30 µl of 10 mM Tris buffer. A 3-µl volume containing 5 parts loading buffer (50 mg/ml blue dextran, 25 mM EDTA, GeneScan-1000 [ROX] red; Applied Biosystems, Foster City, CA) and 1 part standard (GeneScan-1000 [ROX] red; Applied Biosystems, Foster City, CA) was mixed with 2 µl of the cleaned digested product and loaded into a gel sequencer (no. 377; Applied Biosystems, Foster City, CA) to separate terminal fragments. Fragment data were viewed using Genescan and Genotyper software (Applied Biosystems).
Nucleotide sequence accession numbers.
Sequences obtained from the bulk biosolid, bulk soil, and aerosol samples have been deposited in the NCBI GenBank database under accession no. EF029242 to EF029507 (biosolid and aerosol clone libraries) and EF029508 to EF029712 (MST method development and application).
|
|
|---|
![]() View larger version (20K): [in a new window] |
FIG. 1. Phylogenetic distribution of universal and bacterial clones derived from (a) bulk biosolids, (b) aerosols downwind of biosolid spreading operations, and (c) ambient control aerosols.
|
![]() View larger version (31K): [in a new window] |
FIG. 2. 16S rRNA gene-based phylogenetic tree for selected isolates of Chloroflexi. The tree includes clones sequenced from biosolids (B) and soil (S), downwind aerosols (D), and selected type strains isolated from wastewater environments. The tree is based on the neighbor-joining method; the bar indicates 10% estimated sequence divergence. Closed circles represent >50% bootstrap values, and open circles represent <50% bootstrap values. The tree is rooted using E. coli (AM184252) as the out-group (not shown).
|
![]() View larger version (11K): [in a new window] |
FIG. 3. 16S rRNA gene-based phylogenetic tree of clones of Euryarchaeota sequenced from biosolids (B) and soil (S) and selected type strains. Branch I contains 22 sequences from biosolids. The tree is based on the neighbor-joining method. Closed circles represent >50% bootstrap values, and open circles represent <50% bootstrap values. The tree is rooted using Geoglobus ahangari (AF220165) as the out-group.
|
MST verification.
Each MST method was tested against a suite of field-collected bulk biosolid and soil samples to confirm both presence in biosolids and absence in surrounding soils. Soils were used for comparison, as they are the major source of non-biosolid aerosols in the agricultural areas considered.
Sulfite-reducing Clostridia were successfully cultured from all biosolid samples and in 43% of soil samples, although concentrations between the two bulk matrices were markedly different. The average concentrations (mean ± standard deviation) were 3.72 x 106 ± 6.65 x 106 CFU/dry g in biosolids and 1.54 ± 2.27 CFU/dry g for soil samples. The ecologies of culturable sulfite-reducing Clostridia in biosolids and soils were also dissimilar, as shown in Fig. 4. Sequences >99% similar to Clostridium bifermentans (92% of the 78 biosolid sequences) dominated biosolid populations, while sequences similar to Clostridium subterminale (56%) dominated in bulk soil samples. C. bifermentans was identified in only 1 of 15 soil samples (Table 2). Clostridium perfringens strains were commonly isolated from both biosolids and soils. Based on these results, the 16S rRNA gene sequence for C. bifermentans was chosen as a marker for biosolid source tracking. Similarly, the Chloroflexi sp. and Euryarchaeota MST methods were applied to the same group of bulk soil and biosolid samples. Using specific primers, sequences of the Chloroflexi sp. were found in all 14 biosolid samples tested and in 2 of the 14 bulk soil samples analyzed; these sequences were >99% similar to the sequences of Chloroflexi sp. retrieved from the universal biosolid and source aerosol clone libraries. Finally, all bulk biosolid samples analyzed by the tRFLP method produced 251-bp fragments, while only 1 of the 10 soil samples considered showed these fragments (Fig. 5a and b). The average rates of correct classification (ARCC) (defined here as 100 the percentage of positive responses in soils) were 93.3%, 85.7%, and 90% for the Clostridia, Chloroflexi sp., and Euryarchaeota methods, respectively. No single soil sample was positive for all three MST methods, yielding an ARCC of 100% when all MST methods were used to test a sample. Table 2 summarizes the results for MST method testing of the bulk biosolid and soil samples.
![]() View larger version (12K): [in a new window] |
FIG. 4. 16S rRNA gene-based phylogenetic tree sequences from cultured sulfite-reducing Clostridia isolated from biosolids (B) and soil (S) and selected type strains. Branch I contains C. bifermentans sequences from biosolids and from soils, and branch II contains C. subterminale sequences from soil. The tree is based on the neighbor-joining method. Closed circles represent >50% bootstrap values, and open circles represent <50% bootstrap values. The tree is rooted using Geoglobus ahangari (AF220165) as the out-group (not shown).
|
|
View this table: [in a new window] |
TABLE 2. Results of MST testing in bulk soil and bulk biosolids
|
![]() View larger version (9K): [in a new window] |
FIG. 5. Characteristic tRFLP electrophoretic output for the Euryarchaeota MST method. The 251-bp fragment length and fluorescence intensity of amplified 16S rRNA genes are shown for biosolids (a) and were not present in sequences amplified from soils (b).
|
|
View this table: [in a new window] |
TABLE 3. Conditions during field experiments
|
![]() View larger version (16K): [in a new window] |
FIG. 6. Concentrations of sulfite-reducing (red.) Clostridia in aerosols during high-wind events for (a) no-disking and (b) disking scenarios. During the no-disking scenario, two of the three dates for low wind showed no detection and all seven upwind controls showed no detection. For disking, all three upwind controls showed no detection.
|
|
View this table: [in a new window] |
TABLE 4. Results of MST application to aerosol samples
|
![]() View larger version (23K): [in a new window] |
FIG. 7. Time-resolved wind speed and PM10 concentrations during wind aerosolization experiments. Measurements were taken downwind of a field in which biosolids had recently been disk incorporated.
|
|
|
|---|
No MST method produced correct classification rates of 100%. Several factors may have contributed. First, to be representative, soil samples were taken in fields adjacent to and approximately 100 m upwind of the biosolid land application location. Deposition of aerosols produced during nearby biosolid land application or other dairy or livestock activities also could not be unambiguously accounted for. While available records from biosolid land application contractors did not indicate previous biosolid application at these control sites, the possibility of prior, undocumented biosolid or animal compost application could not be definitively excluded. Although these uncertainties do not allow for a strictly controlled field experiment, they are inherent to both the agricultural setting and the biosolid land application process. For these reasons, a multiple-MST approach may be required to produce very high rates (greater than 95%) of correct classification.
Application of MST methods to aerosol samples.
The application of MST methods to aerosols demonstrated that biosolid-derived aerosols could be emitted from land-applied fields during high-wind events. For high-wind events in both the no-disking and disking scenarios, the majority of downwind aerosol samples contained 16S rRNA gene sequences originating from biosolids (13 of 23 downwind). When MST methods were applied to aerosol samples taken at upwind control locations, only 1 of 29 samples was positive. PM10 concentrations (only in disked fields) and culturable analysis of sulfite-reducing Clostridia in aerosol samples also support the conclusion of wind aerosolization (Fig. 6 and 7). The lower ARCC for field-measured aerosols compared to the ARCC for bulk soils may be due to the influence of detection levels: aerosols provide approximately 1 mg of material for analysis, while sample material in bulk soils is unlimited.
Aerosol and bulk biosolid concentrations of sulfite-reducing Clostridia can also be used to estimate the biosolid content of aerosols. Using the average concentration of sulfite-reducing Clostridia in bulk biosolids sampled (3.72 x 106 ± 6.65 x 106 CFU/dry g), the 1- to 20-CFU/m3 range of aerosolized sulfite-reducing Clostridia corresponds to a 0.1- to 2-µg/m3 source aerosol concentration of pure biosolids and suggests limits of detection below this value. This aerosol concentration corresponds to less than 20% of total PM10 in the no-disking cases and less than 2% of the total PM10 in the disking scenario. These biosolid-derived PM10 values are markedly lower than the 500- to 1,000-µg/m3 source concentrations reported during biosolid spreading and disking of dewatered biosolids (29, 31). However, these low values should not be viewed as insignificant since the source emission area could be spread over several square kilometers and produce an equal or greater biosolid mass emission rate into the environment (mass biosolids emitted per time) than those produced by spreading and disking.
The low aerosol concentrations of released biosolids required sensitive MST methods. The direct PCR methods (Chloroflexi sp. and Euryarchaeota) appeared to be more sensitive than the Clostridia method as higher percentages of biosolid-positive results were produced in downwind aerosol samples using the Chloroflexi sp. and Euryarchaeota methods than the Clostridia method. Based on inhibition and inefficiencies in sample processing, realistic PCR detection levels of 100 to 200 cells have previously been established in environmental aerosol samples (32). In contrast, only one spore of sulfite-reducing Clostridia is necessary to form a detectable colony. These differences are negated when considering the enrichment in anaerobically digested biosolids of Chloroflexi sp. and Euryarchaeota (>1% total cells) compared to sulfite-reducing Clostridia (<0.01% total cells) (7, 23, 34). The two culture-independent methods also have greater potential for increased sensitivity. Higher volume filtration may be used for sampling without regard to collection inactivation. In addition, progress in filter extraction efficiency, PCR inhibition removal, and quantitative PCR will continue to improve the sensitivity and quantitativeness of culture-independent aerosol detection.
Finally, downwind concentrations of sulfite-reducing Clostridia and MST results were similar for the no-disking and disking experiments. However, PM10 aerosolization was markedly higher during disking experiments than with the no-disking scenario. These differences underscore the importance of understanding how the land application process affects ambient soil aerosolization. The threshold wind velocity required to aerosolize soil particles is depressed in soils with sandy textures, low soil moisture, and high levels of disturbance (18). Disking clearly disturbed soil structure, and thus larger concentrations of PM10 were released in the disking scenario compared to the no-disking scenario. A consequence of the similar aerosol concentrations of Clostridia for different PM10 levels is that the percentage of biosolid PM10 to total PM10 was approximately 10 times greater in the no-disking case. During disking, the incorporation of biosolids into the subsurface should render most of them unavailable for aerosolization and may partially explain the differences in these ratios. Before a comprehensive set of recommendations on how to reduce wind aerosolization can be made to practitioners, specific mechanistic wind aerosolization studies that include the MST methods developed here will be required to understand how climate, soil, biosolid characteristics, and land application influence wind aerosolization of biosolids. In addition, the three MST methods developed should be tested in new environments using different biosolids to ensure their applicability under all relevant conditions.
Three MST methods have been developed to track the aerosolization of anaerobically digested class B land-applied biosolids. Application of each MST method to field aerosol samples in conjunction with field PM10 measurements provides corroborating evidence that biosolids can become aerosolized during high-wind events. While the results presented here demonstrate the utility of the three biosolid source tracking methods, these methods and results may not be applicable to every land application location, given the diversity of wastewater sources, biosolid processing methods, and land application environments. However, the biosolid and field conditions of this study represent common land application practice. Anaerobic digestion of sewage sludge is the most common method for sludge stabilization (39), and land application is most prevalent in the western and midwestern United States, where arid soil conditions prevail (3). The application of culture-independent MST techniques to aerosols is a new field that has the broader potential for understanding the airborne exposure route and tracking the local, regional, and global transport of aerosols.
We thank Scott Bingham from Arizona State University for assistance with tRFLP and the City of Phoenix for assistance with land application site access.
Published ahead of print on 18 May 2007. ![]()
|
|
|---|
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»