Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2006, p. 3252-3258, Vol. 72, No. 5
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.5.3252-3258.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Sung-Keun Rhee,2,
Chi-Yong Ahn,1
So-Ra Ko,1
Gang-Guk Choi,1
Jin-Woo Bae,3
Yong-Ha Park,3 and
Hee-Mock Oh1,3*
Environmental Biotechnology Laboratory, Korea Research Institute of Bioscience and Biotechnology, 52 Eoeun-dong, Yuseong-gu, Daejeon 305-333, Republic of Korea,1 Department of Microbiology and Biotechnology, Chungbuk National University, 12 Gaeshin-dong, Heungduk-gu, Cheongju, Republic of Korea,2 Biological Resource Center, Korea Research Institute of Bioscience and Biotechnology, 52 Eoeun-dong, Yuseong-gu, Daejeon 305-333, Republic of Korea3
Received 10 October 2005/ Accepted 22 February 2006
|
|
|---|
|
|
|---|
Of the functional genes used for the taxonomic study of cyanobacterial strains, including the cpcBA intergenic spacer (IGS) (7, 15, 17, 20, 22), nifH (7), rpoC1 (24), and gyrB (24), the cpcBA IGS is specific to cyanobacteria and has been widely used for the phylogenetic analysis of pure cyanobacterial culture strains. Baker et al. (1, 2) recently employed a PCR amplification method to analyze the cpcBA IGSs from environmental samples, using a primer set previously designed by Neilan et al. (20), and found a limited cyanobacterial diversity. Although the primer set was originally designed with six cpcBA IGSs to study the genetic diversity of several pure culture strains, cpcBA IGS sequence information from various other cyanobacteria has also been deposited in public databases for potential enhanced primer design.
Cyanobacterial blooms capable of producing microcystins are a seasonal problem every summer in the Daechung Reservoir, which is a representative large eutrophic lake in Korea (21). Accordingly, to further elucidate the composition and dynamics of cyanobacteria during bloom, this study investigated the cpcBA IGS diversity in addition to physicochemical and biological factors. To analyze the cpcBA IGS diversity, new degenerate primers were designed based on more than 300 cpcBA IGS sequences that are currently available in public databases. Finally, the cyanobacterial diversity derived from the cpcBA IGS analysis was compared with that determined by 16S rRNA gene PCR-denaturing gradient gel electrophoresis (DGGE).
|
|
|---|
Water quality analysis.
The water temperature, pH, and conductivity were all measured in situ, using a YSI meter (63/100 FT; YSI Inc., Yellow Springs, OH), while the dissolved oxygen (DO) and turbidity were measured with a DO meter (95/100 FT; YSI Inc., Yellow Springs, OH) and a turbidimeter (DRT-15CE; HF Scientific Inc., Fort Meyers, FL), respectively. The Secchi depth was measured using a Secchi disk. The total N (TN) and P (TP) were determined after persulfate oxidation to nitrate (6) and orthophosphate (18), respectively. The resulting nitrate was then determined by a second-derivative method (4), while the orthophosphate was determined using an ascorbic acid method (8). The total dissolved nitrogen (TDN) and total dissolved phosphorus (TDP) were determined after filtering the water sample through a GF/C filter (Whatman Ltd., Maidstone, United Kingdom) and persulfate oxidation. The total particulate nitrogen (TPN) and total particulate phosphorus (TPP) were obtained by subtracting the TDN from the TN and the TDP from the TP, respectively.
The samples used for plankton identification and enumeration were preserved in Lugol's solution and enumerated with a hemocytometer (Fuchs-Rosenthal; Paul Marienfeld GmbH & Co., Lauda-Königshofen, Germany) under phase-contrast microscopy (Microphot-FXA; Nikon Corp., Tokyo, Japan). The chlorophyll a was extracted using a chloroform-methanol mixture (2:1 [vol/vol]) and measured with a fluorometer (Turner 450; Barnstead/Thermolyne, Dubuque, IA) (29). To determine the picocyanobacterial chlorophyll a, the large cyanobacteria were removed using a filter with a 3-µm pore size, and only the fraction that remained on a 0.2-µm filter was used (10).
Primer design for PCR amplification of cpcBA IGS.
To monitor the diversity and community composition of the cyanobacteria based on the cpcBA IGS, a PCR primer set was specifically designed for cyanobacterial cpcBA IGS amplification. Among all known cpcBA IGS sequences (more than 300 sequences), 24 representative cpcBA IGSs from diverse cyanobacteria (Cylindrospermopsis raciborskii, GenBank accession number AF426795; Cylindrospermopsis raciborskii, AF426792; Arthrospira sp. strain Maxima, AJ401168; Microcystis aeruginosa EAWAG175, AJ003181; Fischerella sp., M75599; Nostoc sp. strain PCC 7120, X05239; Anabaena lemmermannii BC Ana 0018, AY886908; Nodularia spumigena nsb105, AF101444; Synechococcus elongatus, D13173; Synechocystis sp. strain PCC 9413, AF068771; Pseudanabaena sp., M99426; Calothrix sp., M36276; Fremyella diplosiphon, X07012; Synechocystis sp. strain PCC 6803, U34930; Chroococcus dispersus, AJ003184; Agmenellum quadruplicatum, K02660; Synechococcus sp. strain WH7803, X59809; Synechococcus sp. strain PCC 7942, AB008546; Synechococcus elongatus PCC 6301, AP008231; Anacystis nidulans, M94218; Oscillatoria sp. strain PCC 6304, AJ401186; Planktothrix rubescens BC-Pla 9307, AJ131820; Planktothrix sp. strain FP1, AF212923; Gloeobacter violaceus PCC 7421, BA000045; Microcystis sp. strain KLL MG-K, AY524850; Microcystis sp. strain KLL MB-J, AY524849; Microcystis sp. strain KLL MB-K, AY524848; Aphanizomenon sp. Norman Lake isolate, AJ243969; and Aphanizomenon flos-aquae, AJ243971) were aligned by using Clustal X (27). In addition, two conserved regions of cpcB and cpcA were chosen to design the primer set (named CPC1F and CPC1R, respectively) used in this study.
The genomic DNAs from these bacteria were extracted as described previously (13), and the PCR was performed with a GeneAmp PCR System 9700 thermal cycler (Applied Biosystems, Foster City, CA) using a 20-µl (total volume) reaction mixture containing 1x PCR buffer (50 mM KCl, 10 mM Tris-HCl, 0.1% Triton X-100, pH 9.0), four deoxyribonucleoside triphosphates at concentrations of 1 mM, 1.5 mM MgCl2, each primer at a concentration of 1 µM, 4 µg of bovine serum albumin (Roche Diagnostics Corp., Indianapolis, IN), and 2.5 U of Taq DNA polymerase (Perkin-Elmer, Norwalk, CT). The amplification conditions were 25 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s, with a final extension step of 72°C for 7 min.
cpcBA IGS and 16S rRNA gene analysis of bloom samples.
For genomic DNA extraction from the bloom samples, about 500 ml of the water samples was filtered using 0.22-µm membrane filter (Millipore System, Bedford, MA) and the dozens of resulting filter papers were stored at 65°C until required. The bulk community DNA was directly extracted from the filter by grinding the frozen samples and treating them with sodium dodecyl sulfate for cell lysis (13). The conditions used for cpcBA IGS amplification from the environmental samples were the same as those used in the PCR primer evaluation described above. The amplified PCR products were directly ligated into a pCR II vector obtained from Invitrogen (San Diego, CA). The ligation and transformation were carried out as described previously (28). A total of 100 clones per sample were screened, and 50 positive clones per sample were randomly selected and sequenced using the primers CPC1F and CPC1R. The diversity coverage by the clone libraries was analyzed using a rarefaction method, as previously described (25).
The cpcBA IGS sequences obtained in this study were aligned with all those available in current databases by using the Clustal X program (27) and then edited using BioEdit (12). The initial phylogenetic trees were based on all available sequences and constructed using the neighbor-joining DNA distance program (23) in MEGA 2 (16), with bootstrap values based on 1,000 replications (11). Based on the initial phylogenetic results, appropriate subsets of cpcBA IGS sequences were selected and subjected to a final phylogenetic analysis.
Bacterial diversity analysis based on 16S rRNA gene PCR-DGGE was performed as described by Ishii and Fukui (14) and Muyzer et al. (19). A 560-bp (bp 341 to 907 in Escherichia coli numbering) fragment of the 16S rRNA gene from four samples used in the cpcBA IGS analysis was amplified and analyzed using DGGE. The bands were excised from the DGGE gel and incubated in distilled water for 24 h at 4°C. The eluent was then reamplified for sequencing with the primers without a GC clamp.
Nucleotide sequence accession numbers.
The GenBank accession numbers for the clones are as follows: P1-5R, AY942919; P1-30, AY942910; P2-20, AY942916; P2-31, AY942914; P2-48, AY942917; P2-52, AY942907; P3-26, AY942908; P3-36, AY942918; P3-51, AY942909; P3-59, AY942915; P4-9, AY942906; P4-21, AY942905; P4-25, AY942912; P4-29, AY942913; P4-30, AY942911; P4-56, AY942920; 16S-1, AY942899; 16S-2, AY942896; 16S-3, AY942902; 16S-4, AY942903; 16S-5, AY942894; 16S-6, AY942895; 16S-7, AY942897; 16S-8, AY942900; 16S-9, AY942904; 16S-10, AY942901; and 16S-11, AY942898.
|
|
|---|
![]() View larger version (32K): [in a new window] |
FIG. 1. Changes in chlorophyll a, TPP, and TPN concentrations (A) and number of total phytoplankton (B) in Daechung Reservoir. The sampling was carried out at 1-week intervals. The samples between the dashed lines were used for further molecular analysis. The error bars indicate standard deviations (n = 3).
|
|
View this table: [in a new window] |
TABLE 1. Biogeochemical characteristics of bloom water
|
Design of primers for PCR amplification of the cpcBA IGS.
Although the reported cpcB sequences are quite diverse (nucleotide similarity of 266 bp at the 3' end of the cpcB gene, 72.3 ± 7.4% [mean ± standard deviation]; range of pairwise similarity, 56 to 98%), several highly conserved regions were useful for the primer design. The priming site of PCbF (20), positioned at nucleotide 250 according to the numbering for Microcystis sp. strain KLL MB-J (accession no. AY524849) cpcB, was selected for the design of the new forward primer CPC1F (5'-GGCKGCYTGYYTRCGYGACATGGA-3'). Since several mismatches were found in the PCbF primer, degeneracy was applied to increase the coverage of CPC1F. In addition, a new priming site near the 5' end of cpcA (nucleotide 43 according to the numbering for Microcystis sp. strain KLL MB-J [accession no. AY524849] cpcA) was explored for the reverse primer PC1R (5'-GCHGATWCYCAAGGNCGYTT-3'). Among the primer sets previously reported (3, 20) and those designed in this study, the CPC1F-CPC1R pair produced the most successful amplification of representative test strains from the cyanobacterial collections of PCC, UTEX, NIES, SAG, and the Korean Collection for Type Cultures (KCTC) (data not shown). As shown in Fig. 2, cpcBA IGSs were amplified from diverse cyanobacteria. As expected, due to size variations in the IGS region, the amplicon size ranged between 400 and 500 bp, indicating that the newly designed primer set was able to cover cpcBA IGSs from diverse cyanobacteria.
![]() View larger version (40K): [in a new window] |
FIG. 2. Photographs of ethidium bromide-stained gels showing amplification products of cpcBA IGSs from reference cyanobacteria. The strains used were as follows: 1, Aphanothece nidulans KCTC AG10041; 2, Aphanothece naegelii KCTC AG10042; 3, Microcystis aeruginosa UTEX 2385; 4, Microcystis sp. strain PCC 7806; 5, Synechococcus sp. strain PCC 7002; 6, Synechocystis sp. strain PCC 6803; 7, Merismopedia tenuissima NIES 230; 8, Arthrospira maxima SAG 49.88; 9, Arthrospira platensis NIES 39; 10, Oscillatoria tenuis UTEX 1566; 11, Oscillatoria tenuis NIES 33; 12, Planktothrix agardhii NIES 204; 13, Anabaena flos-aquae UTEX 2517; 14, Anabaena affinis KCTC AG10008; 15, Anabaena sp. strain KCTC AG10059; 16, Aphanizomenon flos-aquae NIES 81; 17, Nodularia spumigena UTEX 2092; 18, Nostoc sp. strain PCC 7120; 19, Chlorogloeopsis sp. strain PCC 9212. Lanes M, markers (1-kb ladder). The numbers on the left indicate the sizes of the markers after electrophoresis.
|
![]() View larger version (39K): [in a new window] |
FIG. 3. Phylogenetic relationships of cpcB genes cloned from Daechung Reservoir. The reference sequences were from GenBank. The nucleotide sequences (267 bp from the 3' end) from cpcB were used for the phylogenetic analysis. cpcB genes with a similarity of more than 95% are grouped. A representative clone from each group is indicated in parentheses. The bar indicates 0.05 substitution per nucleotide position. The local bootstrap probabilities are indicated at the nodes if they are larger than 50%.
|
Dynamics of the cyanobacterial community during bloom.
To monitor the changes in the cyanobacterial community during cyanobacterial bloom, the relative composition of the cpcBA IGS sequences in each sample was analyzed. Groups D, E, and F were the major components of the cyanobacterial community on 19 August 2003, as shown in Fig. 4. The two most significant changes around the bloom peak on 2 September 2003 were the gradual decrease to extinction of group E and gradual increase and sudden disappearance of group D on 9 September. Meanwhile, group F steadily increased and finally predominated on 9 September. Interestingly, after the bloom peak, the cyanobacterial community profile did not return to that observed before the peak, as group F accounted for up to 85% of all the clones. However, the total cyanobacterial count, including Microcystis, decreased (down to 30% of that at the bloom peak), as shown in Table 1.
![]() View larger version (48K): [in a new window] |
FIG. 4. Distribution dynamics of the cpcBA IGS sequence group during cyanobacterial bloom in Daechung Reservoir. The frequency distribution was calculated based on 50 clones analyzed from each sample. The dates of sampling are indicated.
|
Comparison of cpcBA IGS and 16S rRNA gene analyses for cyanobacterial bloom study.
The cpcBA IGS-based method was compared with a 16S rRNA gene-based method using 16S rRNA gene PCR-DGGE. The best matches with the 16S rRNA gene sequences from the DGGE bands were identified using BLAST (Table 2). As expected, 16S rRNA genes for bacteria and chloroplasts were amplified together with those for cyanobacteria. Dynamic changes in the bacterial community were observed during the bloom period in the DGGE profile (Fig. 5), and based on a comparison of the phylogenies of the 16S rRNA and cpcBA IGS and their succession patterns, tentative correlations were found between the cpcBA IGS sequence groups and cyanobacterial DGGE bands. For example, the sequences of DGGE band 10 and cpcBA IGS group F, which closely matched those of Microcystis (>98% and 100% similarities, respectively), were present throughout the bloom period, while the sequences of DGGE band 8 and cpcBA IGS group E, which closely matched those of Aphanizomenon and Anabaena, occupied a major proportion on 19 August 2003. Group E and band 8 were present as a minor component in the sample taken from 26 August to 9 September 2003. However, no Aphanizomenon- or Anabaena-like microorganisms were observed by microscopy on 19 August 2003 (Table 1). Considering that Aphanizomenon and Anabaena were observed in over 10% of cyanobacteria on 12 August 2003, Aphanizomenon and Anabaena were presumed to have existed on 19 August, even at low densities, but were not observed by microscopic enumeration. These results indicate that molecular techniques are more sensitive than microscopic observation for community study. Another reason could be the different sample volumes for each analysis, because about a 10-times-greater volume was used for molecular methods than for microscopy. Small flock formation by rarely occurring Aphanizomenon and Anabaena and their heterogeneous distribution could have amplified such a discrepancy.
|
View this table: [in a new window] |
TABLE 2. Identities of bands obtained from DGGE analysis of bloom samplesa
|
![]() View larger version (69K): [in a new window] |
FIG. 5. 16S rRNA gene-DGGE profile during cyanobacterial bloom in Daechung Reservoir. The reproducibility of the DGGE profile was confirmed by three replications of the PCR-DGGE of the genomic DNA. The dates of sampling are indicated.
|
Improvement of water quality management depends on the development of proper methods to analyze cyanobacterial diversity and determine the cause of toxin production from blooms. In the present study, a molecular analysis of the cpcBA IGS was shown to be efficient for analyzing cyanobacterial diversity, not least because it can exclude other bacteria. Thus, the present protocol could complement conventional tools, such as morphology-based microscopic analyses, the direct determination of toxins, and 16S rRNA gene analyses. Furthermore, since the IGS length of each group of cyanobacteria is different, the proposed cpcBA IGS amplification protocol could also be directly coupled to length polymorphism analysis techniques, such as terminal restriction fragment length polymorphism, to facilitate rapid monitoring of the cyanobacterial community.
S.-G. Kim and S.-K. Rhee contributed equally to this work. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»