Evidence for polyphyletic origin of the members of the orders of Oscillatoriales and Pleurocapsales as determined by 16S rDNA analysis

Evidence for polyphyletic origin of the members of the orders of Oscillatoriales and Pleurocapsales as determined by 16S rDNA analysis

FEMS Microbiology Letters 201 (2001) 79^82 www.fems-microbiology.org Evidence for polyphyletic origin of the members of the orders of Oscillatoriale...

298KB Sizes 0 Downloads 7 Views

FEMS Microbiology Letters 201 (2001) 79^82

www.fems-microbiology.org

Evidence for polyphyletic origin of the members of the orders of Oscillatoriales and Pleurocapsales as determined by 16S rDNA analysis Tatsuya Ishida a

a;

*, Makoto M. Watanabe b , Junta Sugiyama

a;c

, Akira Yokota

a

Institute of Molecular and Cellular Biosciences, The University of Tokyo, Bunkyo-ku, Tokyo 113-0032, Japan b National Institute for Environmental Studies, 16-2 Onogawa, Tsukuba, Ibaraki 305-0053, Japan c NCIMB Japan Co., Ltd., Kaminakazato O¤ce, 9-2 Sakae-Cho, Kita-Ku, Tokyo 114-0005, Japan Received 16 January 2001; received in revised form 18 May 2001; accepted 19 May 2001 First published online 15 June 2001

Abstract Phylogenetic analyses were carried out on six pleurocapsalean strains and 12 oscillatorialean strains by sequence determination of the 16S rDNA (16S rRNA gene). Although heterocyst-forming strains of the orders Nostocales and Stigonematales were shown to be monophyletic, unicellular strains of the orders Chroococcales and filamentous Oscillatoriales were shown to be polyphyletic, as reported earlier. Moreover, unicellular and baeocyte-forming strains of the order Pleurocapsales, which were thought to be monophyletic, were newly found to be polyphyletic. The results strongly indicate that even morphology was not necessarily reflected in the phylogenetic relationships at the order level. A need exists to reconstruct the taxonomy of cyanobacteria at the order level. ß 2001 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved. Keywords : Cyanobacterium ; Pleurocapsales ; 16S rDNA ; Phylogeny

1. Introduction Cyanobacteria are autotrophic bacteria that perform oxygenic photosynthesis. The taxonomy of cyanobacteria has until now been based mainly on their morphology, and cyanobacteria are now separated into ¢ve orders [1]: Chroococcales (I), Pleurocapsales (II), Oscillatoriales (III), Nostocales (IV), and Stigonematales (V). Phylogenetic analyses of cyanobacteria based on 16S rDNA were carried out by several groups [2^6], which indicated that Chroococcales (I) and Oscillatoriales (III) were polyphyletic. Other orders, Nostocales (IV) and Stigonematales (V), which form heterocysts, were shown to be monophyletic. Of the other orders, monophylicity was indicated for Pleurocapsales (II) which forms internal spores (baeocytes); however, only a few strains were phylogenetically analyzed. Furthermore, the sequences determined in pleu* Corresponding author. Present address: International Patent Organism Depository, National Institute of Advanced Industrial Science and Technology, Higashi 1-1-1, Tukuba, Ibaraki 305-8566, Japan. Tel.: +81 (298) 61-6075; E-mail : [email protected]

rocapsalean strains (II) were incomplete. As in the other orders of cyanobacteria, the taxonomy of Pleurocapsales (II) was based on its morphology [7]. The morphology of cyanobacteria is known to change easily and is not conserved [8]. To grasp a general image of Pleurocapsales (II) and to get more information on Oscillatoriales (III), in this study we investigated the phylogenetic position of the six pleurocapsalean strains (II) and 12 oscillatorialean strains (III) based on 16S rDNA sequences, which have been proved to be very useful in the study of bacterial phylogeny [9]. 2. Materials and methods 2.1. Cyanobacterial strains The strains examined were as follows : Chroococcidiopsis sp. PCC 6712 and C. thermalis PCC 7203; Stanieria sp. PCC 7301 and S. cyanosphaera PCC 7437T ; Pleurocapsa sp. PCC 7319 and PCC 7327 ; Geitlerinema sp. PCC 7105; Leptolyngbya sp. PCC 7104 and PCC 7375; Lyngbya aes-

0378-1097 / 01 / $20.00 ß 2001 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved. PII: S 0 3 7 8 - 1 0 9 7 ( 0 1 ) 0 0 2 4 6 - 4

FEMSLE 10000 29-6-01

80

T. Ishida et al. / FEMS Microbiology Letters 201 (2001) 79^82

tuarii PCC 7419T ; Oscillatoria acuminata PCC 6304 and O. sancta PCC 7515T ; Pseudanabaena sp. PCC 6802, PCC 6903, PCC 7367, PCC 7403, and PCC 7408 ; Symploca sp. PCC 8002. All strains were obtained from the Pasteur Culture Collection (PCC), and all were axenic. Pleurocapsa sp. PCC 7319 and Pseudanabaena sp. PCC 7367 were grown in medium ASN-III [10], Geitlerinema sp. PCC 7105, Leptolyngbya sp. PCC 7375, and Stanieria sp. PCC 7301 in medium ASN-III.B12 [10], Symploca sp. PCC 8002 in medium ASO -III [10], Chroococcidiopsis sp. PCC 6712, C. thermalis PCC 7203, Leptolyngbya sp. PCC 7104, O. acuminata PCC 6304, O. sancta PCC 7515T , Pleurocapsa sp. PCC 7327, Pseudanabaena sp. PCC 6802, PCC 6903, PCC 7403, PCC 7408, and S. cyanosphaera PCC 7437T in medium BG-11 [10], and L. aestuarii PCC 7419T in medium MN [10]. These strains were cultured by incubation in test tubes without aeration at 20³C with £uorescent lamp illumination of 500 lux. The names of strains used fundamentally followed the catalog of PCC [10]. 2.2. PCR ampli¢cation and sequencing The methods for DNA extraction, ampli¢cation, and sequencing of 16S rDNA were the same as those described previously [3], except for PCR oligonucleotide primers and conditions. The almost complete 16S rDNAs from the genomic DNA of the respective strains were ampli¢ed by PCR, using oligonucleotide primers by the combination of 1 (5P-AGAGTTTGATCCTGGCTCAG-3P) and 16 (5PACGCCGACCTAGTGGAGGAA-3P) [11]. PCR was performed at ¢rst for 3 min at 94³C and then 30 cycles with the following features: 1 min at 94³C, 1 min at 55³C, and 2 min at 72³C, followed by a ¢nal elongation step for 10 min at 72³C. 2.3. Alignment and phylogenetic analyses The new 16S rDNA sequences were multiple-aligned using CLUSTAL W, version 1.74 [12] with a selection of cyanobacterial reference sequences obtained from the DNA Data Bank of Japan (DDBJ). The alignment was corrected manually and converted to a distance matrix. A phylogenetic tree was constructed from the distance matrix data by applying the algorithm of the neighbor-joining (NJ) method [13] to Knuc values [14] with multiple substitutions corrected and positions with gaps excluded. To evaluate the robustness of branches in the inferred tree, the bootstrap resampling method of Felsenstein [15] with 1000 replicates was used. Maximum likelihood (ML) analysis was also carried out using the MOLPHY program, version 2.3b3 [16] based on the same alignment data after positions with gaps were excluded. The ML distance matrix was calculated using NucML, and the initial NJ tree was reconstructed by NJdist in the MOLPHY package. The ML tree was ¢nally obtained using NucML with the local rearrangement method from the NJ tree with

the HKY model, and local bootstrap probabilities were estimated by a resampling of the estimated log-likelihood method [17]. 3. Results and discussion The new 16S rDNA sequences were deposited in DDBJ ; accession numbers for each strain are given in Fig. 1. Sequences from 58 cyanobacterial strains and two Grampositive bacterial strains (Bacillus subtilis subsp. subtilis DSM 10 and Clostridium butyricum ATCC 19398) as the outgroups were aligned. Positions with gaps and undetermined and ambiguous sequences were removed. A total of 1214 sites were used for the phylogenetic analysis. Constructed phylogenetic trees, the ML tree (Fig. 1) and the NJ tree (not shown), revealed that all the nostocalean and stigonematalean strains, which are di¡erentiated by the formation of heterocysts in their trichomes, were monophyletic. The unicellular strains, Chroococcales (I) and Pleurocapsales (II), were indicated, as the ¢lamentous oscillatorialean strains (III) were, to be polyphyletic and did not form heterocysts. In previous reports [2^5], except for the pleurocapsalean strains (II), the same results were obtained, and they indicated that the pleurocapsalean strains (II) were monophyletic; however, only a few pleurocapsalean strains (II) were analyzed in that report. In this study, we included more pleurocapsalean strains (II) in the phylogenetic analyses. As shown in Fig. 1, these strains are clearly separated into at least three clusters. The major one consists of six pleurocapsalean strains (II) and one prochloralean strain (VI). The others consist of only one pleurocapsalean strain (II) and strains from other orders ; one was an oscillatorialean strain (III) and the other was a chroococcalean strain (I). In this study, some new sequences were added to cyanobacterial phylogenetic analysis, and it became clear that the order Pleurocapsales (II) is polyphyletic. That is, pleurocapsalean strains (II) of which the phylogenetic positions have been analyzed until now were thought to be placed in a narrow range. Namely, only a part of the cyanobacteria were examined in the phylogenetic analyses, and these analyses were thought to be insu¤cient under the present conditions. At present, there are only Nostocales (IV) and Stigonematales (V) of which the morphological characteristics as the guideline to the order were consistent with the phylogeny (Fig. 1). Cyanobacteria seemed to have diverged at the same time in an early period in the cyanobacterial phylogenetic tree [2^6], and it is di¤cult to elucidate the process of their evolution. This indicates that diversi¢cation of cyanobacteria once occurred during a certain early period. As a result, chroococcalean strains (I) and oscillatorialean strains (III) might be seen in the same cluster on some phylogenetic position in the present phylogenetic trees. It is interesting, however, that the chroococcalean strain (I) Synechococcus sp. PCC 7002

FEMSLE 10000 29-6-01

T. Ishida et al. / FEMS Microbiology Letters 201 (2001) 79^82

81

Fig. 1. ML tree of 58 cyanobacterial strains and two Gram-positive bacterial strains (Bacillus subtilis subsp. subtilis DSM 10 and Clostridium butyricum ATCC 19398) as the outgroups. An alignment of 1214 nucleotides after excluding the positions with gaps was used. Scale bar = 1 base substitution per 10 nucleotide positions. Local bootstrap probabilities are indicated at nodes if larger than 80. The strains of which the 16S rRNA genes were determined in this study are indicated in bold. Accession numbers in the DDBJ, EMBL, and GenBank databases are in parentheses. Order VI: Prochlorales; T: type strain.

FEMSLE 10000 29-6-01

82

T. Ishida et al. / FEMS Microbiology Letters 201 (2001) 79^82

and the oscillatorialean strain (III) Oscillatoria rosea M-220 probably diverged recently. Thus unicellular and ¢lamentous strains are close, and it may help in the determination of morphological change from unicellular to ¢lamentous strain or ¢lamentous to unicellular strain. In this study, we examined not only pleurocapsalean strains (II), but also oscillatorialean strains (III). The 16S rDNA sequences of 12 oscillatorialean strains (III) were determined and their phylogenetic positions analyzed. As shown in Fig. 1, oscillatorialean strains (III) are separated into eight clusters supported by high statistical reliability. Three groups consisted of a few oscillatorialean strains (III), and therefore the sequence of 16S rDNA of more strains should be determined. Although almost oscillatorialean strains (III) did not form a consistent cluster even at the genus level, only the Pseudanabaena species characterized by the presence of a cross-wall [18] formed a single cluster. Thus, all the morphological characteristics at the genus level were not always signi¢cant. Namely, a taxonomic system of cyanobacteria re£ecting phylogenetic relationships should be reconstructed. Acknowledgements This work was supported by Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists. References [1] Castenholz, R.W. and Waterbury, J.B. (1989) Taxa of the cyanobacteria. In: Bergey's Manual of Systematic Bacteriology (Staley, J.T., Bryant, M.P., Pfenning, N. and Holt, J.G., Eds.), Vol. 3, pp. 1727^ 1728. Williams and Wilkins, Baltimore, MD. [2] Nelissen, B., Wilmotte, A., Neefs, J.-M. and De Wachter, R. (1994) Phylogenetic relationships among ¢lamentous helical cyanobacteria investigated on the basis of 16S ribosomal RNA gene sequence analysis. Syst. Appl. Microbiol. 17, 206^210. [3] Ishida, T., Yokota, A. and Sugiyama, J. (1997) Phylogenetic relationships of ¢lamentous cyanobacterial taxa inferred from 16S rRNA sequence divergence. J. Gen. Appl. Microbiol. 43, 237^241.

[4] Honda, D., Yokota, A. and Sugiyama, J. (1999) Detection of seven major evolutionary lineages in cyanobacteria based on the 16S rRNA gene sequence analysis with new sequences of ¢ve marine Synechococcus strains. J. Mol. Evol. 48, 723^739. [5] Turner, S., Pryer, K.M., Miao, V.P.W. and Palmer, J.D. (1999) Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis. J. Eukaryot. Microbiol. 46, 327^338. [6] Garcia-Pichel, F., Lopez-Cortes, A. and Nu«bel, U. (2001) Phylogenetic and morphological diversity of cyanobacteria in soil desert crusts from the Colorado plateau. Appl. Environ. Microbiol. 67, 1902^1910. [7] Waterbury, J.B. (1989) Subsection II. Order Pleurocapsales. In: Bergey's Manual of Systematic Bacteriology (Staley, J.T., Bryant, M.P., Pfenning, N. and Holt, J.G., Eds.), Vol. 3, pp. 1746^1770. Williams and Wilkins, Baltimore, MD. [8] Pearson, J.E. and Kingsbury, J.M. (1966) Culturally induced variation in four morphologically diverse blue-green algae. Am. J. Bot. 53, 192^200. [9] Woese, C.R. (1987) Bacterial evolution. Microbiol. Rev. 51, 221^271. [10] Rippka, R. and Herdman, M. (1992) Pasteur Culture Collection of Cyanobacterial Strains in Axenic Culture. Catalogue and Taxonomic Handbook, Vol. I. Catalogue of Strains. Institut Pasteur, Paris. [11] Wilmotte, A., van der Auwera, G. and De Wachter, R. (1993) Structure of 16S ribosomal RNA of the thermophilic cyanobacterium Chlorogloepsis HTF (`Mastigocladus laminosus HTF') strains PCC 7518, and phylogenetic analysis. FEBS Lett. 317, 96^100. [12] Thompson, J.D., Higgins, D.G. and Gibson, T.J. (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position speci¢c gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673^4680. [13] Saitou, N. and Nei, M. (1987) The neighbor-joining method: A new method for reconstructing ohyligenetic trees. Mol. Biol. Evol. 4, 406^ 425. [14] Kimura, M. (1980) A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16, 111^120. [15] Felsenstein, J. (1985) Con¢dence limits on phylogenies: An approach using the bootstrap. Evolution 39, 783^791. [16] Adachi, J. and Hasegawa, M. (1996) MOLPHY: Programs for Molecular Phylogenetics, Version 2.3. Institute of Statistical Mathematics, Tokyo. [17] Hasegawa, M. and Kishino, H. (1994) Accuracies of the simple methods for estimating the bootstrap probability of a maximum-likelihood tree. Mol. Biol. Evol. 11, 142^145. [18] Castenholz, R.W. (1989) Subsection III. Order Oscillatoriales. In: Bergey's Manual of Systematic Bacteriology (Staley, J.T., Bryant, M.P., Pfenning, N. and Holt, J.G., Eds.), Vol. 3, pp. 1771^1780. Williams and Wilkins, Baltimore, MD.

FEMSLE 10000 29-6-01