Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2009, p. 5121-5130, Vol. 75, No. 15
0099-2240/09/$08.00+0 doi:10.1128/AEM.00447-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Ecology and Evolutionary Biology, UCB 334, University of Colorado, Boulder, Colorado 80309,1 Cooperative Institute for Research in Environmental Sciences, UCB 216, University of Colorado, Boulder, Colorado 80309,2 National Center for Atmospheric Research, Boulder, Colorado 80301,3 Department of Computer Science, UCB 430, University of Colorado, Boulder, Colorado,4 Desert Research Institute, Storm Peak Laboratory, 2215 Raggio Parkway, Reno, Nevada 89512,5 Department of Chemistry and Biochemistry, UCB 215, University of Colorado, Boulder, Colorado 803096
Received 23 February 2009/ Accepted 2 June 2009
|
|
|---|
|
|
|---|
Microbes are often considered passive inhabitants of the atmosphere, dispersing via airborne dust particles. However, recent studies suggest that many atmospheric microbes may be metabolically active (3, 4, 64), even up to altitudes of 20,000 m (34). Some airborne microbes may alter atmospheric conditions directly by acting as cloud condensation nuclei (7, 25, 56) and/or ice nuclei (IN) (19, 41, 56, 57, 61); this hypothesis is supported by the observation that most ice nuclei in snow samples are inactivated by a 95°C heat treatment (16, 17). However, the overall contribution of airborne microbes to atmospheric processes such as ice nucleation remains unclear.
The best-studied ice-nucleating microbes are gram-negative bacteria that have also been isolated from leaf surfaces, including Pseudomonas syringae, Pseudomonas fluorescens, Erwinia herbicola, Xanthomonas campestri, and Sphingomonas spp. (45). These bacteria have been cultured extensively, and their ice-nucleating activity has been traced to a membrane-bound glycoprotein (40, 42, 70). However, their specific influence on atmospheric processes remains, at this point, largely anecdotal. Less is known about the ice-nucleating activities of fungi, but a few studies have shown that fungi can be effective ice nucleators, capable of initiating ice nucleation at temperatures as high as –2°C (41, 61). At this point, all known ice-nucleating microorganisms are amenable to culture-based studies, but given that the vast majority of microorganisms have yet to be cultured, it is likely that other ice-nucleating microbes remain undiscovered.
The work presented here addresses three overarching questions. (i) Are microbial abundances altered by changes in atmospheric conditions? (ii) How is the diversity and composition of airborne microbial communities influenced by changes in atmospheric conditions? (iii) Can we identify known and novel ice-nucleating microbes in the atmosphere by testing for correlations between taxa abundances and the concentrations of biological ice nuclei? To address these questions, we combined epifluorescence microscopy, tagged pyrosequencing, Sanger sequencing, and an ice nucleation assay with atmospheric measurements to characterize the microbial communities at a high-elevation research site.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Description of samples collected and corresponding meteorological conditions
|
Airborne microbial sample collection.
Air samples were captured on triplicate sterile 0.22-µm cellulose nitrate filters (Fisher Scientific, Pittsburgh, PA) via vacuum filtration with a flow rate of 7.5 liters min–1. The total volume of air that passed through each filter during each sampling period ranged from 2.5 m3 to 5.4 m3. One filter from each triplicate set was used for each of the following assays: total microbial abundance measurements via epifluorescence microscopy, DNA extraction and the associated microbial community composition analyses, and ice nucleation assays. Additionally, two fresh snow samples were collected following overnight storms at SPL that occurred during the sampling period (Table 1). Snow was melted at 4°C in sterile bags, and 1 liter of snowmelt was filtered through each of three 0.22-µm cellulose nitrate filters. The same three assays described above were also performed on the two collected snow samples.
Total microbial abundance.
Aerosol particles were shaken from the filters at approximately 200 rpm in a small petri dish in 8 ml of high-performance liquid chromatography (HPLC)-grade water for 2 h at 4°C. A subset of the filters was examined under a microscope after this shaking process to assure that most visible particles were removed from the filter. While the efficiency of particle removal was not evaluated for every individual sample, all sample filters were treated identically, and therefore any bias associated with this method should have been held constant across the sample set. The microbial particles suspended in HPLC-grade water were stained with 4',6'-diamidino-2-phenylindole (DAPI), a DNA binding dye (Kirkegaard and Perry, Gaithersburg, MD) at a final working concentration of 500 µg ml–1 and counted at 1,000x magnification using a Nikon Eclipse E400 epifluorescence microscope following the protocol described by Hernandez et al. (36). Intact cells were counted on 25-mm-diameter black polycarbonate filters with pore sizes of 0.22 µm (GEI-W&PT, Trevose, PA). Microbial abundance is expressed as cells per cubic meter of air, taking into account the dilution, the flow rate, and the number of hours sampled.
Ice nucleation counts.
Particles were removed from the filters using the same protocol described in the total particle count section above. The drop-freeze assay was performed using a modified version of the protocol described by Vali (72). Drop-freeze plates were assembled by coating an aluminum foil-covered copper plate with a solution of 5% (wt/vol) petrolatum in xylene. After the xylene had evaporated, 10-µl drops were arrayed (96 drops per sample) in an 8 by 12 format across the surface of the plate. Extra columns of 10-µl, HPLC-grade water drops (not containing sample) were placed on each side of the array as a plate control. The temperature of the plate was then cooled with a Peltier device at a constant rate of 1°C every 30 s until all of the water drops froze. The drop-freeze plate preparation was considered to be robust if the HPLC-grade water drops did not freeze above –20°C. All sample drops freezing at –10°C and above were recorded as containing ice nuclei and were used to calculate the abundance of ice nuclei in each sample. The total concentration of ice nuclei in each sample was calculated using equation 13 under the cumulative nucleus spectra section of Vali's 1971 method (72). The cumulative nucleus concentration calculation is sensitive to the total number of drops used as well as the volume of each drop in the assay. For example, if the sample is aliquoted into a large number of small-volume drops, more total ice nuclei will be observed than if the same sample were aliquoted into fewer, larger-volume drops. The calculation described for Vali's 1971 method (72) determines the total concentration of ice nuclei per ml of sample volume. The total number of ice nuclei was then calculated per cubic meter of air or per liter of snowmelt for the air and snow samples, respectively. A subset of samples was reanalyzed using replicate filters in order to demonstrate the repeatability of the assay, and the ice nucleation numbers varied by less than 1% between replicate filters.
Total microbial community composition via Sanger sequencing.
We examined the overall microbial community composition of bacterial and fungal communities in air and snow samples using a method similar to that of Lauber et al. (44) by constructing clone libraries of partial rRNA genes followed by Sanger sequencing. DNA from microbial cells trapped on the filters was extracted using the Ultra-Clean plant DNA isolation kit (MoBio Laboratories, Carlsbad, CA). The filters were cut into strips using sterile scissors and forceps under a laminar flow hood and loaded into the bead tube of the DNA extraction kit and heated to 65°C for 10 min, followed by 2 min of vortexing. The remaining steps of the extraction were carried out according to the manufacturer's instructions. Ribosomal sequences from all three domains (Archaea, Bacteria, and Eukarya) were amplified using the universal primer set 515f and 1391r, which amplifies an 850- to 1,100-bp portion of the small-subunit rRNA gene (6). PCRs were performed using the Platinum PCR Supermix (Invitrogen, Carlsbad, CA). Primers were added to a final concentration of 0.2 µM, and 2 µl of DNA template was used in each 25-µl reaction mixture. PCR amplification was performed with an initial denaturation at 94°C for 3 min followed by 35 cycles of 94°C for 30 s, 50°C for 30 s, and 72°C for 30 s with a final extension at 72°C for 10 min. Each sample was amplified in triplicate, and the amplicons were pooled by sample and run on a 1.2% agarose gel. The bands were purified with the Qiaquick gel extraction kit (Qiagen, Carlsbad, CA). PCR amplicons were cloned into the TOPO TA cloning vector, transformed into TOP10 chemically competent cells (Invitrogen, Carlsbad, CA), and grown on L-agar supplemented with 50 µg ml–1 ampicillin. Following overnight incubation at 37°C, colonies were picked and M13 colony PCR was performed using identical PCR conditions to those described above. PCR products were sent to Agencourt Bioscience (Beverly, MA) for single-pass sequencing. Negative control filters were extracted and PCR amplified and no-template PCRs were run. These reactions produced no visible PCR product on a SYBR green-stained agarose gel.
Sanger sequence analysis.
Sequences were binned into major taxonomic groups (i.e., bacteria, fungi, and plants) using the BLAST algorithm (2) against the NCBI nonredundant database (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Sequences with expect (E) values greater than 1e–100 were not included in the analysis. Of the 672 Sanger sequences having an average length of 750 bp, 12% were of low quality or chimeric and were therefore removed from the analysis. Approximately 7% of the quality sequences were most closely matched to plants and these sequences were not further analyzed. Bacterial and fungal sequences were analyzed separately, as there is currently no single database that can simultaneously classify both of these groups with confidence. The bacterial sequences were aligned using the NAST aligner (minimum length, 600 bp; minimum identity, 75%), checked for chimeras using Bellerophon (21) against the Greengenes database (20), and classified with the Greengenes Classify tool (http://greengenes.lbl.gov/cgi-bin/nph-classify.cgi) using the Hugenholtz taxonomy. More detailed phylogenetic analysis of specific bacterial groups was performed with the maximum likelihood tree-building program RAxML-7.0.4 (68) using the GTR nucleotide substitution model, and a consensus tree was built from the highest-scoring tree of 100 bootstrapped trees. Fungal sequences were aligned using MUSCLE (26) and classified by performing a local BLAST analysis against the comprehensive 16S/18S Silva database, r96_ssu.fasta file (http://www.arb-silva.de/no_cache/download/archive/release_96_exports/). The number of unique fungal phylotypes was defined as those sequences sharing
97% sequence similarity, estimated with the FastGroupII algorithm (74).
Bacterial community composition via tagged pyrosequencing.
We used a newly developed bar-coded pyrosequencing procedure to comprehensively survey the airborne bacterial communities. This protocol is identical to the protocol described by Fierer et al. (32), including both the PCR conditions and primer sequences. Negative controls (both no-template control and template control from unused filters) were included at all steps of the process to check for contamination. PCR amplicons from each sample were pooled at approximately equimolar concentrations into a single tube used for the pyrosequencing reaction. Samples were pyrosequenced at the Environmental Genomics Core Facility at the University of South Carolina on a 454 Life Sciences genome sequence FLX (Roche) machine.
Pyrosequencing phylogenetic analysis.
Sequences were analyzed and processed according to the methods of Hamady et al. (35). Only those sequences that were >200 bp in length with a quality score of >25 with no ambiguous characters were included in the analysis (37). Sequences were assigned to samples by examining the 12-bp bar code (32). Phylotypes were identified using Megablast to identify connected components (nearest neighbor), sets of similar sequences (with the parameters E value of 1e–8, minimum coverage of 99%, and minimum pairwise identity of 97%). A representative sequence was chosen from each phylotype by selecting the most highly connected sequence, i.e., the sequence that had the most hits more significant than the BLAST threshold compared to other sequences in the data set (46). The set of all representative sequences was aligned using NAST (21) (with the parameters of minimum alignment length of 190 and sequence identity of 70%) with a PH lanemask (http://greengenes.lbl.gov/) to screen out hypervariable regions of the sequence. Sequences were checked for chimeras using an in-house BLAST-based chimera checking program, which basically checks to see if the sequence ends hit different phyla. A related neighbor-joining tree was built using Clear-cut (67), employing the Kimura correction. Unweighted UniFrac (50, 51) was run using the resulting tree and the sequences annotated by environment type. The UniFrac algorithm provides an estimate of the overall phylogenetic distance between each pair of communities and therefore avoids some of the pitfalls associated with phylotype-based community analyses (51). Taxonomic identity of the phylotypes was assigned with BLAST against the Greengenes database (20) by using an E value cutoff of 1e–10 and the Hugenholtz taxonomy.
Statistical analysis.
UniFrac was used to quantitatively compare the phylogenetic structures of the individual air and snow samples (51). The UniFrac algorithm was applied to a neighbor-joining tree containing all sequences across all samples together with an environment file used to match sequence identity numbers (IDs) to their corresponding sample ID. Relationships between microbial community phylogenetic distance, as determined using the UniFrac algorithm, and meteorological parameters were investigated using the PRIMER software package (www.primer-e.com) (18). Nonmetric multidimensional scaling was used to visualize the UniFrac phylogenetic distance between the air and snow bacterial communities in each pair of environments. The RELATE function in PRIMER was used to evaluate correlations between community phylogenetic distance and each of the individual weather parameters. The effect of each individual weather parameter (Table 1) on community phylogenetic distance was considered significant if P was <0.05. The variations in total particle abundance, total microbial abundance, and the cumulative ice nuclei between clear or cloudy skies were assessed with Student's t test using Statview v5.01.
Nucleotide sequence accession numbers.
The Sanger sequences have been deposited in the GenBank nr database. The bacterial sequences have accession numbers FJ746985 through FJ747245, and the fungal sequences have accession numbers FJ747246 through FJ747530. The authors are currently working on depositing the pyrosequencing data in the GenBank short-read archive; please contact the corresponding author if the pyrosequencing data for this study are desired.
|
|
|---|
Particle counts in relation to total microbial abundance and total ice nuclei.
The number of total particles and DAPI-stained particles (i.e., cells) remained relatively constant across all nine sampling periods. Total microbial counts ranged from 9.6 x 105 to 6.6 x 106 cells m–3 of air, and the concentrations of total particles with diameters of 530 nm to more than 5 µm, as determined via APS size fractionation, ranged between 8.2 x 107 and 2.3 x 108 particles m–3 (Fig. 1A). There was no significant difference between clear and cloudy air samples with respect to either cell concentrations (P = 0.5) or total particle counts (P = 0.10). However, the total number of ice nuclei was far more variable across the 2-week sampling period than either the cell counts or total particle counts, ranging from 0 IN/m3 to 91 IN/m3 (Fig. 1). When samples were binned into clear (n = 4) or cloudy (n = 5) air samples, the mean IN concentration was significantly greater in the cloudy air samples (P < 0.05) (Fig. 1A). Similarly, the mean number of IN was positively correlated with relative humidity measurements across all of the sampling events (Fig. 1B).
![]() View larger version (15K): [in a new window] |
FIG. 1. Total APS particle counts, total microbial abundance, and total number of ice nuclei observed in clear and cloudy skies and in replicate snow samples. (A) The left y axis shows APS particle counts (gray bars) (P = 0.09 for clear versus cloudy) and DAPI counts (DNA-containing particles) (black bars) (P = 0.47 for clear versus cloudy) for samples collected in clear air, cloudy air, and from fresh snow. The right y axis shows the number of IN observed per volume of air or melted snow (white bars) (P = 0.05 between cloudy and clear). The air samples are in units of particles or IN m–3 of air, while the snow samples are in units of particles or IN per liter of snow melt. There were no statistics run between the air and snow samples, as the units cannot be compared, and there were no total particle measurements, as the APS instrument was only used for atmospheric aerosols. (B) Relationship between the mean relative humidity during each of the air sampling periods and the total number of ice nuclei m–3 of air.
|
![]() View larger version (55K): [in a new window] |
FIG. 2. Coarse level taxonomic description of the nine air samples and two snow samples. The numbers in parentheses represent the total number of sequenced clones from each sample (E < 1e–100) that were nonchimeric with an E value of <1e–100. Refer to Table 1 for a complete description of each sample ID listed on the x axis.
|
![]() View larger version (41K): [in a new window] |
FIG. 3. Composition of the fungal communities as determined using full-length Sanger sequencing. Numbers in parentheses indicate the number of sequences obtained from each sample. Refer to Table 1 for a complete description of each sample listed on the x axis.
|
97% identity). Sequences were also grouped at the 95% and 99% sequence similarity levels, which yielded 571 and 1,229 unique phylotypes, respectively. Rarefaction curves were constructed for the most and least diverse samples at each of these sequence similarity levels (Fig. 4). Rarefaction curves at all three similarity levels did not include an asymptote, suggesting that we did not survey the full extent of airborne bacterial diversity, even with an average of 407 sequences per sample (Fig. 4).
![]() View larger version (24K): [in a new window] |
FIG. 4. Rarefaction curves at the 99%, 97%, and 95% sequence similarity levels. Curves shown are from the most diverse (sample ID 4-4.CL.N; solid lines) and least diverse (sample ID 3-24.CL.D; dashed lines) of the nine airborne bacterial samples.
|
![]() View larger version (35K): [in a new window] |
FIG. 5. Most abundant bacterial groups identified using bar-coded pyrosequencing. Numbers in parentheses indicate the number of sequences obtained from each sample (E < 1e–100) that were nonchimeric with an E value of <1e–100 from a single tagged pyrosequencing run. Refer to Table 1 for a complete description of each sample ID listed on the x axis. Proteobacterial groups are designated by the Greek symbols , β, and .
|
![]() View larger version (5K): [in a new window] |
FIG. 6. Nonmetric multidimensional scaling plot of the UniFrac distance matrix from the tagged pyrosequence bacterial data set. These plots compress the multidimensional UniFrac distance matrix to two dimensions. (A) Air and snow data. Stress value, 0.01. (B) Data for air only; IDs correspond to the IDs listed in Table 1. Stress value, 0.07.
|
|
|
|---|
Recently, primary biological particles have been shown to contribute a large fraction to the total atmospheric aerosol abundance (38). Of the total aerosol concentration at SPL, the number of particles with biological characteristics (fluorescing in the presence of the DNA binding dye DAPI) accounted for approximately 5 to 10% of the particles measured by the APS (Fig. 1A). This microbial aerosol fraction is similar to estimates provided in the recent literature. For example, air samples that were collected in Germany (55), the South Atlantic (53), and at a remote location in Siberia (54) contained microbial fractions (percentage of fluorescing particles in relation to the total particle count) ranging from 17% to 30%. While our study did not include measurements spanning the different seasons, the microbial aerosol fraction was relatively constant over the 2-week sampling period (Fig. 1A), which is consistent with the temporal stability observed by the Matthias-Maser (53, 54, 55) and Jaenicke (38) studies mentioned above.
Bacteria/fungi ratios.
The relative abundance of the major taxonomic groups (plant, bacteria, and fungi) varied little over the 2-week sampling period (Fig. 2). The plant sequences were minor constituents of the nine libraries constructed from the air samples, with fungal and bacterial sequences having approximately equal representation in the libraries (Fig. 2). The 1:1 ratio of bacterial and fungal sequences observed in these air samples is consistent with that reported by Wilson et al. (73), who observed a similar ratio in a single air sample collected in California (73). The bacteria/fungi ratios at SPL differed from those of both a recent ribosomal sequencing-based study (33) and an earlier culture-based survey (24), as both studies observed much higher variation in their bacterial to fungal ratios over time. These differences may be attributed to differences in sampling location and/or the time of sampling. The air collected by Fierer et al. (33) was collected at head height on a college campus, and the di Giorgio et al. (24) work was conducted at two cities in France over a year. Both studies were performed at locations that are likely influenced by a number of anthropogenic disturbances that could potentially cause large shifts in microbial abundance and composition over short time periods, whereas the air at SPL is comparatively undisturbed with regard to surface-level microbial sources. While there are approximately equal numbers of bacterial and fungal sequences in our airborne samples, it is important to recognize that these sequence ratios do not necessarily indicate that bacteria and fungi have equivalent biomass concentrations in the atmosphere, as these taxa differ with respect to their cell sizes and the number of small-subunit rRNA copies per cell.
Fungal community composition.
At the level of taxonomic resolution afforded by the 18S rRNA gene, the fungal sequences appear to be more variable than the bacterial communities (Fig. 3), although some of this variability could be attributed to the small number of sequences used to assess the airborne fungal communities. The majority of the samples collected at SPL were dominated by the Dothideales and Eurotiomycetes classes of fungi (Fig. 3), which were also observed to be the dominant fungi colonizing sediment exposed by a receding glacier (39). Culturable fungi that are commonly thought to be abundant in the atmosphere, such as Alternaria, Cladosporium, and Penicillium spp. (15, 66), were not found in the air samples at SPL. These differences may be attributed to culturing bias, although we did not compare cultivation-dependent and cultivation-independent surveys of fungal diversity in this study.
Bacterial diversity.
Just as the relative abundances of airborne bacteria and fungi were moderately constant over time, the types of bacteria present in the atmosphere sampled at SPL also changed very little over the 2-week sampling period (Fig. 5). Although the 4,864 bacterial sequences generated from the bacteria-specific pyrosequencing run represent one of the more comprehensive surveys of airborne bacterial diversity to date (with an average of 407 sequences per sample), we still did not survey the full extent of bacterial diversity within the individual air samples, as evidenced by the lack of asymptotes in Fig. 4. The bacterial diversity levels observed at SPL are considerably lower than those reported in other molecular-based studies of the atmosphere (13) and in other environments such as in soil (31). The differences in diversity levels between the various airborne landscapes may be attributed to sampling location, as Despres et al. (22, 23) observed lower bacterial diversity levels in the air of the high-alpine environment compared to their urban and rural sampling locales (22, 23).
Bacterial community composition.
The bacterial taxa identified from the nine air samples collected at SPL were similar to those found in other studies where air samples were also dominated by Betaproteobacteria and Gammaproteobacteria (22, 23, 31, 52, 63). In particular, we found that most of the Betaproteobacteria sequences were comamonads (members of the order Burkholderiales) (Fig. 5), which are a metabolically and phenotypically diverse group of bacteria that is relatively common in the atmosphere (13, 33). The dominant Gammaproteobacteria taxon, the Moraxellaceae, was highly represented across the nine air samples (Fig. 5), with Psychrobacter spp. being the most abundant genus-level group within this family. Interestingly, members of the genus Psychrobacter have previously been isolated from various low-temperature environments, such as Antarctic soils (10), Antarctic sea ice (11, 12), and Siberian permafrost (60).
Bacteria in the atmosphere are distinct from those found in fresh snow.
Previous studies have used snow samples to predict microbial abundances and the numbers of ice-nucleating particles of biological and nonbiological origin present in the atmosphere (16, 17). This approach would seem legitimate, as precipitation events are known to scavenge and collect particles that eventually get deposited onto the ground (62). However, the bacterial communities found in the air above SPL appear to be distinct when compared to the bacterial communities found in freshly fallen snow (Fig. 5 and 6). The two snow samples harbored distinct bacterial communities, although this variation was likely due to the exposure of sample 4-3.SN to intense sunlight for approximately 2 hours prior to sampling, which may explain the overrepresentation of cyanobacterial sequences in this sample. Snow sample 4-7.SN was much more representative of the bacterial communities found in the air, although there were still a number of bacterial groups found in this snow sample that were not present in any of the air samples (Fig. 5), indicating that freshly fallen snow is not truly representative of the microbial life residing in the atmosphere.
Shifts in total IN in response to shifting weather conditions.
While the types and abundances of microorganisms in the air samples were fairly constant over time and not strongly influenced by changes in atmospheric conditions, the abundance of ice-nucleating particles was far more variable, increasing with increases in relative humidity (Fig. 1). Other studies have also shown that ice nucleation rates in cirrus and in wave clouds respond positively to increases in humidity (71). The ice nuclei counted in this study were counted at temperatures equal to or higher than –10°C, which has recently been used as a temperature cutoff to define ice nucleators of biological origin (16, 17). The most potent nonbiological ice-nucleating particles generally initiate freezing at temperatures below –10°C (65). Based on these studies, we made the assumption that the IN counted via the drop-freeze assay were likely to have been initiated by bacteria, fungi, and/or plant pollen. Throughout our bacterial community analysis, we identified ribosomal DNA gene sequences with homology to the following known ice-nucleating bacteria: P. syringae, P. fluorescens, E. herbicola, and Sphingomonas spp. (Fig. 7). However, these bacterial taxa were not common in either our Sanger sequence or pyrosequence data sets. For instance, the Pseudomonadaceae family, which contains the most extensively studied ice-nucleating species, represented less than 1% of the snow pyrosequences and only 3% of the air pyrosequence data set. However, sequences matching the Psychrobacter genus were relatively common in our data sets, with similar abundance levels across all nine air samples (Fig. 5). Members of this genus have recently been shown to initiate ice nucleation at temperatures as high as –2°C (60), which rivals the ice-nucleating temperature of the well-studied bacterium P. syringae. The phylogeny presented in Fig. 7 illustrates the phylogenetic relatedness between the 16S rRNA sequences generated from this current study to 16S rRNA sequences of many of the verified ice-nucleating bacterial isolates.
![]() View larger version (31K): [in a new window] |
FIG. 7. Phylogenetic distribution of representative full-length Sanger sequences that have similarities to the 16S rRNA gene of known ice-nucleating bacteria. Sequences were aligned using the Greengenes NAST aligner, and the maximum-likelihood tree was made using RAxML-7.0.4. A Streptomyces sequence was used as the outgroup. Each Sanger sequence is represented by its sample ID (see Table 1 for sample descriptions).
|
This study was supported by the CIRES Innovative Research Program. Support was also provided by Steamboat and Resort Corporation. The Desert Research Institute is an equal opportunity service provider and employer and a permittee of the Medicine-Bow Routt National Forests.
Published ahead of print on 5 June 2009. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»