Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2003, p. 6634-6643, Vol. 69, No. 11
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.11.6634-6643.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Microbial Ecology, Centre for Limnology, Netherlands Institute for Ecology, 3631 AC Nieuwersluis,1 Aquatic Microbiology, Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, 1018 VW Amsterdam, The Netherlands2
Received 5 May 2003/ Accepted 12 August 2003
|
|
|---|
|
|
|---|
For ecological studies of cyanobacteria, it is essential to be able to differentiate between closely related organisms. Some important features, such as nitrogen fixation and toxicity, are not uniformly dispersed across taxa or morphological groups and may vary between species or even between strains of the same species (9, 41). For example, toxic cyanobacterial blooms often require characterization below the genus level. Both toxin-producing and -nonproducing strains are included in most morphologically distinguishable species of the genera Microcystis (18, 30) and Anabaena (5). Shifts in the ratio of toxic to nontoxic genotypes of Microcystis may (partly) explain the often poor correlations of cell counts with toxin content in many field studies (6).
To study the functions and interactions of cyanobacteria, it is necessary to reveal the composition and dynamics of their populations. Also, it is important to relate isolated strains to their counterparts in nature for extrapolation of findings from physiological experiments carried out on cultures to natural conditions. Linking laboratory cultures and field populations is often problematic due to the selectivity inherent to cultivation and to morphological changes occurring after cultivation (27, 31). Denaturing gradient gel electrophoresis (DGGE), a technique for sequence-dependent separation of PCR products (25), can be used to assess the genotypic diversity in environmental samples and to judge the purity and uniqueness of isolated strains. Isolated cultures can be assigned to field populations based on the comparison of their DGGE profiles, and sequence information from profile bands can be used to characterize the organisms that are present. A section of DNA is suitable for DGGE analysis if it can be specifically amplified from the target organisms, has sufficient sequence heterogeneity for the desired resolution, and, preferably, is part of a gene for which a considerable number of sequences have been deposited in sequence databases. DGGE of hetR, a gene assigned to heterocyst differentiation, has been used to study diversity in isolated strains of the cyanobacterial genera Trichodesmium and Nostoc (28, 32). Also, nifH, a gene encoding dinitrogenase reductase in many microorganisms, including cyanobacteria (42), has been used for DGGE analysis of the very diverse functional group of diazotrophic (N2-fixing) organisms (22). An important drawback for the application of both protein-encoding genes, especially for the analysis of communities of cyanobacteria, is that they are present in only a limited number of cyanobacterial genera. Nübel et al. (27) developed primers for specific amplification of a 16S rRNA gene segment from cyanobacteria and plastids, allowing DGGE analysis of cyanobacterial populations (1, 11, 27). However, the taxonomic resolution offered by 16S rRNA genes is insufficient for discrimination of closely related organisms. As a result, research has increasingly focused on the rRNA 16S to 23S internal transcribed spacer (rRNA-ITS). The greater degree of sequence heterogeneity, as well as a considerable number of published rRNA-ITS sequences, makes rRNA-ITS very suitable for high-resolution analysis of cyanobacteria. Restriction enzyme digestion of rRNA-ITS has been used to resolve closely related cyanobacterial strains (21, 23, 26, 32), and direct sequencing has been used to study subgeneric phylogenetic relationships in genera such as Microcystis (30), Trichodesmium (28), and picocyanobacteria (34). Recently, DGGE has been used for analysis of Synechococcus (4) and Aphanizomenon (20) rRNA-ITS sequences.
In the present work, we combined the relatively high sequence variation in rRNA-ITS with the potential of DGGE to separate even small differences in sequence to introduce methods for high-resolution analysis of cultures and populations of cyanobacteria. We developed and tested primers and protocols for selective amplification and DGGE analysis of cyanobacterial genotypes in the presence of DNA from contaminating microorganisms. Their value for analysis of cyanobacterial strains of various genera and of field samples was examined, and the resolution that could be achieved by different rRNA DGGE methods was compared by using closely related Microcystis strains.
|
|
|---|
Strain isolation and culture conditions.
The strains of cyanobacteria used in this study are given in Table 1. The cultures isolated in this study (Microcystis and Gloeotrichia spp.) were obtained by picking single colonies from water samples by using a dissecting microscope and sterile glass Pasteur pipettes with narrow openings, followed by repeated washings in a small volume of O2 medium (8). Alternatively, single Microcystis colonies were grown on plates containing O2 medium solidified with agarose (0.3%, wt/vol). The colonies grown on plates were rendered unialgal through repeated plating. None of the cyanobacterial cultures that were used in this study were axenic. All cyanobacteria able to grow in freshwater medium were maintained in our laboratory in O2 medium (5 ml in 30-ml tubes or 30 ml in 100-ml Erlenmeyer flasks) at 20°C at an irradiance of approximately 15 µmol · m-2 · s-1 for a cycle with a light-to-dark ratio of 16:8 h. A Microcoleus isolate was maintained as decribed by Jonkers et al. (16). The bacteria used as negative controls to test the specificity of the primers were Serratia marcescens (DSM 1636), Acinetobacter calcoaceticus strain BD4 (DSM 586), Erwinia carotovora subsp. carotovora (DSM 30168), Escherichia coli (DSM 423), Rhodococcus erythropolis (DSM 43188), Lactobacillus reuteri (DSM 20016), Lactococcus lactis subsp. lactis (NIZO-81), Bacillus polymyxa (DSM 36), Bacillus subtilis (DSM 10), and Pseudomonas stutzeri (DSM 5190). Cultures were obtained from the German Collection of Microorganisms and Cell Cultures (cultures encoded DSM) and from the Dutch Institute for Dairy Research (encoded NIZO).
|
View this table: [in a new window] |
TABLE 1. Cultures of cyanobacteria used in this study
|
PCR amplification.
Primer sequences and references are given in Table 2. The forward primer for amplification of part of the 16S rRNA gene was slightly modified from that used by Nübel et al. (27) to further improve their theoretical specificity, because the efficiency of amplification and therefore the selectivity of the primer is typically determined by the nature of the 3'-end nucleotides (19, 35). PCR amplification was performed in an MBS 0.5 S thermocycler (ThermoHybaid, Ashford, United Kingdom) in a 25-µl reaction mixture containing approximately 50 ng of DNA, 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 0.01% (wt/vol) of gelatin, 200 µM (each) deoxynucleotide, 1.5 mM MgCl2, 2.5 U of Taq DNA polymerase (Boehringer Mannheim, Mannheim, Germany), and 0.5 µM (each) primer. The temperature cycling conditions for the amplification of part of the 16S rRNA were modified slightly from those of Nübel et al.(27). After preincubation at 94°C for 5 min, a total of 30 cycles were performed at 94°C for 1 min, 60°C for 1 min, and 72°C for 1 min. The temperature cycling was concluded with a final step of 5 min at 72°C. The optimized temperature cycling conditions for the amplification of rRNA-ITS (ITSa, ITSb, and ITSc; Results and Table 2) were as follows. After preincubation at 94°C for 5 min, a total of 30 cycles were performed at 94°C for 1 min, annealing temperature (Ta) for 1 min, and 72°C for 1 min. In the first 20 cycles, Ta decreased by 1°C after every second cycle, from 62°C in the first cycle to 52°C in the twentieth. This touch down procedure was followed to reduce nonspecific annealing of the primers. In the last 10 cycles, Ta was 52°C. The temperature cycling was concluded with a final step of 30 min at 72°C.
|
View this table: [in a new window] |
TABLE 2. Primers used in this study
|
Sequencing of DNA from DGGE bands.
A small piece of gel from the middle of the target band was excised from the DGGE gel and incubated in 50 µl of sterile Milli-Q-purified water for 24 h at 4°C. The eluent was reamplified by using the original primer set and run on DGGE to confirm its identity. For sequence analysis, the eluent was reamplified with reverse primers that had M13 priming sites added to the 5' side of the original primers. The PCR products were purified by using the Concert Rapid PCR Purification System (GibcoBRL Life Technologies, Paisley, United Kingdom), and these products were used as templates for sequencing reactions with the Thermo Sequenase Primer cycle sequencing kit (Amersham Pharmacia Biotech). Sequencing reaction products were analyzed on an ALF Express II sequencer (Amersham Pharmacia Biotech) with CY5 fluorescence-labeled M13 sequence primers or a labeled G-C clamp. The sequences were processed by using the program Sequencher, version 4.0.5 (Gene Codes Corp., Ann Arbor, Mich.), and similarity with sequences deposited in the GenBank, EMBL, and DDBJ databases was checked by using the program BLAST (2) (http://www.ncbi.nlm.nih.gov/BLAST/).
Template mixtures experiment.
ITSa was amplified from the selected cultures, and for each strain the amount of PCR product was estimated from an agarose gel. The PCR products were mixed and diluted in order to obtain template DNA from selected cultures or from mixtures thereof, all in identical final concentrations. This template DNA was used for ITSa amplification, and the PCR products were analyzed on a DGGE gel.
Nucleotide sequence accession numbers.
The sequences were deposited at EMBL and were assigned accession numbers AJ579895 to AJ579906.
|
|
|---|
Reverse primer sites for the amplification of cyanobacterial rRNA-ITS sequences were gathered from previous research (Table 2). Three reverse primers were used, two (373R and ITS3R) targeted at highly conserved sequence motifs in the ITS and one (ULR) at the 5' end of the 23S gene. The combination of the forward primer (CSIF) and the reverse primers yields three primer sets, ITSa (primers CSIF plus 373R), ITSb (CSIF plus ITS3R), and ITSc (CSIF plus ULR) (Table 2). The rationale for using more than one primer combination for amplification of the rRNA-ITS was the trade-off that exists between the broad applicability of the shorter fragments and the different melting behavior and higher amount of sequence information associated with longer fragments (see below). Alignment of cyanobacterial rRNA-ITS sequences deposited in the GenBank, EMBL, and DDBJ databases showed that virtually all cyanobacteria had at least one operon containing a perfect match with the reverse primer 373R, which is targeted to the highly conserved tRNAIle gene. Only some Prochlorococcus strains contained one mismatch. Most deposited cyanobacterial sequences also had a perfect match with the ITS3R primer sequence. Only a number of picocyanobacteria did not contain the target sequence. A BLAST search performed on the reverse primer sequences resulted in the highest similarity scores for cyanobacteria. However, it cannot be concluded that the reverse primers are selective for cyanobacteria, since the database of rRNA ITS and 23S sequences is relatively small (compared to that for 16S sequences) and may be biased for cyanobacteria.
Theoretical melting curves calculated for rRNA-ITS sequences of several major genera (using the program MELT94) showed that the G-C clamp positioned at the forward (5') primer yielded a favorable melting curve, with a lower melting domain at the 3' end.
In vitro specificity and applicability of the primers for rRNA-ITS amplification.
The forward primer presented in Table 2 had a complete or almost complete match with cyanobacterial sequences, but it also had some sequence similarity with several noncyanobacterial 16S rRNA sequences in the database. To confirm the selectivity of the primers for amplification of cyanobacterial rRNA-ITS, we tried to amplify DNA from bacterial strains with some degree of sequence homology to the forward primer (four to five mismatches) and from a few randomly chosen strains. No PCR products were generated from these bacteria either at the optimized annealing temperature (as described below) or at 10° below this temperature.
The phylogenetic coverage within the cyanobacterial phylum of the primers was investigated by using DNA from strains of a range of different cyanobacterial genera as templates for PCR amplification (Table 3). For each primer combination, the highest temperature at which PCR products were generated for all genera was selected as the optimal PCR amplification temperature. As a consequence, for several genera amplification is possible at higher annealing temperatures than those used in our optimized protocol. The number of PCR products and their sizes and relative intensities varied among the different genera (Table 3). For all strains that were tested, ITSa amplification yielded one band of 275 to 350 bp. ITSb amplification resulted in PCR products between 350 and 800 bp for most tested strains, with the exception of Synechocystis, Leptolyngbya, and Lyngbya. One PCR product was formed from the genera Microcystis, Pseudanabaena, Planktothrix, Trichodesmium, Prochlorothrix, and Synechococcus, and two PCR products were formed from Aphanizomenon, Anabaena, Anabaenopsis, Gloeotrichia, Nodularia, Nostoc, and Cylindrospermopsis. ITSc primers successfully amplified DNA from all the strains that were tested, yielding for most genera PCR products ranging between 450 and 900 bp. Prochlorothrix hollandica yielded one very long product of 1,050 bp. Amplification with ITSb or ITSc primers yielded the same number of bands for most strains; only ITSc amplification of Planktothrix strains produced two bands instead of one.
|
View this table: [in a new window] |
TABLE 3. Number and approximate sizes of the PCR products of strains from various cyanobacterial genera amplified with the primers described in this studya
|
![]() View larger version (91K): [in a new window] |
FIG. 1. Discrimination of closely related cyanobacteria. DNA from 20 Microcystis strains was amplified with four different primer combinations (A through D) and separated on DGGE gels. The primer sequences and combinations are given in Table 2. The different primer combinations amplified a segment of the 16S gene (A) or different portions of the rRNA-ITS, i.e., ITSa (B), ITSb (C), or ITSc (D). The Microcystis strains are SAG17.85 (lane 1), V80 (lane 2), V72 (lane 3), V73 (lane 4), CYA43 (lane 5), V67 (lane 6), K29 (lane 7), K50 (lane 8), CYA140 (lane 9), CYA228 (lane 10), PCC7820 (lane 11), PCC7806 (lane 12), V91 (lane 13), Z6 (lane 14), Z11 (lane 15), V89 (lane 16), V28 (lane 17), S2 (lane 18), V88 (lane 19), and V40 (lane 20). Information about the Microcystis strains is given in Table 1. The DGGE gels had a 30 to 40% denaturant concentration gradient, and only the part of the gel containing bands is shown.
|
The presence of several homologous DNA templates may affect PCR amplification and result in heteroduplex (see above) and chimera formation. To investigate the occurrence of these phenomena, we mixed ITSa amplification products from seven Microcystis strains that could be separated by using DGGE and used these mixtures as templates for ITSa amplification. Figure 2 shows the DGGE banding patterns generated from template mixtures of two, three, five, or seven Microcystis strains. In spite of the presence of several related sequences, each strain in the mixtures gave rise to one single band of similar intensity and no additional bands were detected.
![]() View larger version (61K): [in a new window] |
FIG. 2. DGGE profiles of ITSa amplified from DNA template mixtures. Lanes 1 through 4 show ITSa PCR products which were amplified from equal concentrations of template mixtures of Microcystis strains PCC7820 plus K29 (lane 1), PCC7820 plus K29 plus SAG17.85 (lane 2), Z6 plus PCC7820 plus CYA140 plus K29 plus SAG17.85 (lane 3), and S2 plus Z6 plus PCC7820 plus CYA140 plus K29 plus CYA43 plus SAG17.85 (lane 4).
|
![]() View larger version (100K): [in a new window] |
FIG. 3. DGGE profiles of various cyanobacterial genera amplified with ITSa and ITSc primers. The following strains were analyzed: Synechococcus sp. (lane 1), Synechocystis sp. (lane 2), Leptolyngbya sp. (lane 3), Lyngbya sp. (lane 4), Pseudanabaena catenata (lane 5), Trichodesmium erythraeum (lane 6), Anabaena variablis (lane 7), Aphanizomenon flos-aquae (lane 8), Aphanizomenon gracile (lane 9), Planktothrix aghardii (lane 10), Planktothrix prolifica (lane 11), Cylindrospermopsis raciborskii (lane 12), Anabaenopsis arnoldi (lane 13), Gloeotrichia sp. (lane 14), Nostoc sp. (lane 15), Nodularia spumigena (lane 16), and Prochlorothrix hollandica (lane 17). The gels had a denaturant concentration gradient of 25 to 40%, and only the part of the gel containing bands is shown. The lower gel was assembled from two gels.
|
![]() View larger version (83K): [in a new window] |
FIG. 4. DGGE profiles generated from ITSa amplification products from freshwater lake samples. The samples were taken from Lake Zeegerplas on July 20 (lane 1) and August 21 (lane 2) when cyanobacteria were blooming and from Lake Kinselmeer on May 5 (lane 3) during a cyanobacterial bloom and on January 16 in the absence of a visible presence of cyanobacteria (lane 4). Several dominant bands were excised, reamplified, and sequenced.
|
|
|
|---|
Primer 322, described by Wilmotte et al. (40), and primer 16CITS, described by Neilan et al. (26), are targeted to sequences just upstream from the target sequence of forward primer CSIF and thus could produce virtually identical amplification products. Primer 322, however, is targeted to highly conserved bacterial sequences and is therefore suitable only for studies of axenic cultures of cyanobacteria (14, 21). Primer 16CITS is specific for cyanobacteria. However, from our alignment of 16S sequences it appears that primer CSIF has more sequence differences with noncyanobacterial sequences, especially at the 3' site. For primer 16CITS, we found 71 noncyanobacterial sequences with two mismatches (as opposed to none for primer CSIF) and 266 noncyanobacterial sequences with three mismatches (as opposed to two for CSIF). Primer CSIF can therefore be considered more specific for cyanobacteria. Nevertheless, we emphasize that this was based on sequences in the present database and that no experimental comparisons of the primers were made.
Coverage of a broad range of genera for primer CSIF was confirmed by the probe match software tool by using published cyanobacterial sequences and was supported by in vitro tests (Table 3 and Fig. 3). We were able to generate PCR products from Nostoc strain CYA124, even though some deposited sequences from this genus had three mismatches with the primer sequence. Investigations must still be done to determine whether PCR products can be formed from all strains from the genera Nostoc and Spirulina (two or three mismatches in deposited sequences), from Gloeobacter and Calothrix (three mismatches), and from genera without deposited sequences. The generation of amplification products from strains of all tested cyanobacterial genera with primer 373R as the reverse primer and from most genera with primer ITS3R is supported by the universal occurrence of the highly conserved gene for tRNAIle and of certain structural elements in cyanobacteria (7, 13, 34).
The difference in resolution with 16S and ITS DGGE in Microcystis is in agreement with the reported average sequence diversity of less than 1% in 16S (7, 29) compared to up to 7% in rRNA-ITS (26, 30). The high resolution and the single bands that are generated make rRNA-ITS DGGE a valuable method for analysis of this genus. We made the assumption that the multiple bands containing Microcystis sequences in the rRNA-ITS DGGE profiles of some cultures reflect the coisolation of two different Microcystis strains instead of the occurrence of two different rRNA operons. This assumption can be confirmed unequivocally only by repeated plating to make certain that the analyzed cultures are unialgal. This elaborate procedure was not carried out, since we considered the possibility of coisolation of two morphologically indistinguishable organisms far more likely than the occurrence of many different operons in a limited number of deviant strains. Also, all strains that were most likely unialgal (because they originated from culture collections or had been isolated by using plating techniques) gave rise to one band on DGGE. Following this line of reasoning, the explanation for the two sequences in culture collection strain SAG17.85 (Fig. 1C and D) must be that they acquired unnoticed contamination during years of cultivation in the lab. Support for the occurrence of only one operon (or several identical operons) in the genus Microcystis comes from the unambiguous sequences that Otsuka et al. (30) obtained from sequencing the rRNA-ITS of 47 Microcystis strains.
The occurrence of multiple operons in one organism resulting in rRNA-ITS amplicons of different sizes (Table 3) has been found before in the genera Anabaena, Aphanizomenon, Anabaenopsis, Nostoc, Cylindrospermopsis, and Nodularia (13, 14, 26). Sequencing of the rRNA-ITS from Nodularia, Nostoc, Calothrix, and Scytonema strains (3, 7, 13) revealed that not all rRNA operons of cyanobacteria contain tRNAIle genes. However, since other operons in these organisms do contain tRNAIle genes, they are detectable with ITSa amplification (as shown for Nodularia and Nostoc in Fig. 3A). The presence of tRNA genes on only one copy of the rRNA may explain for a number of genera the production of only one band after ITSa amplification (Fig. 3A and Table 3), despite the existence of different operons revealed by ITSb or ITSc amplification (Fig. 3B and Table 3). Alternatively, generation of a single ITSa DGGE band could be explained by identical sequences at the 5' end of different rRNA-ITS operons.
For the generation of DGGE profiles from longer stretches of rRNA-ITS, primer set ITSc is preferable in most cases. Compared to ITSb analysis, the resolution was similar (compare Fig. 1C and D), yet more genera could be amplified and there is more sequence information contained in the amplicons (60 to 100 bp extra; Table 3). Nevertheless, for some genera, primer set ITSb was more suitable for DGGE profiling. For instance, ITSb amplification products of Nodularia yielded sharp DGGE bands (data not shown), whereas ITSc DGGE resulted in diffuse bands unsuitable for analysis (Fig. 3B). Also, DGGE profiles resulting from ITSb amplification of Planktothrix aghardii were relatively simple (one band) compared to the profile resulting from ITSc amplification (three bands). Due to the occurrence of multiple rRNA-ITS operons in one organism, identification of all dominant bands produced with ITSb and ITSc primers is necessary to come to a reliable estimate of the diversity of cyanobacteria in complex DGGE profiles.
It must be emphasized that DGGE diversity profiles do not necessarily reflect the true diversity in the field. DNA extraction efficiency and the number of rRNA operons may vary between genera, and the PCR step has several inherent pitfalls, which complicate the interpretation of DGGE profiles of communities. Artifacts introduced by the PCR amplification reaction of homologous sequences include chimera (17) and heteroduplex formation (15), which may lead to overestimation of the number of organisms, and template annealing (36) and preferential amplification of some DNA templates (33), which may result in a shift in the apparent ratio of the different organisms after PCR amplification. Heteroduplex bands were formed in some cultures suspected of containing more than one Microcystis isolate (Fig. 1) and in some organisms containing multiple operons (Fig. 3). In contrast, deliberate mixing of template DNA from several Microcystis strains showed no evidence of chimera or heteroduplex formation, since no additional bands were detected in the DGGE profiles (Fig. 2). Equal band intensities pointed to the absence of preferential amplification. However, this effect is more likely between genera, since primer binding may differ due to imperfect matching.
Conclusion.
By using the primers and protocols described in this paper, cyanobacterial communities and isolates can be studied at high resolution. The selectivity of the primers makes it possible to focus on cyanobacteria in the presence of other organisms. The rationale for using different rRNA-ITS primer sets has become clear from our results. Each set has its advantages and disadvantages, and the choice of which one to use depends on the study that is to be executed. For studying cyanobacteria, ITSa primers have the broadest applicability and produce the most straightforward DGGE profiles. However, bands from some organisms may end up at identical positions in the gel even though sequence differences do exist (either because sequence differences occur in the nonamplified 3' part of the rRNA-ITS or because different sequences dictate identical melting behavior). The information obtained from DGGE of ITSa may be insufficient to describe an ecosystem in detail or to make sure that a culture contains one cyanobacterial isolate. The longer stretches of rRNA-ITS that are amplified by ITSb and ITSc primers make detection of additional sequence differences possible and thus provide supplementary resolution (Fig. 1). Also, these amplicons contain more sequence information, which is especially interesting when field profiles are analyzed. The disadvantages of ITSb primers are their inability to amplify DNA of some genera (Table 3) and the complex DGGE profiles they produce for some genera. Most sequence information can be retrieved from ITSc primer amplicons. Moreover, for most tested genera, they yielded sharp bands on DGGE gels.
Josje Snoek is gratefully acknowledged for isolation and maintenance of cyanobacterial cultures. Miguel Dionisio-Pires, Olav M. Skulberg, Sven Jansson, and Henk M. Jonkers are gratefully acknowledged for supplying cyanobacterial cultures.
|
|
|---|
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»