Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2008, p. 2822-2833, Vol. 74, No. 9
0099-2240/08/$08.00+0 doi:10.1128/AEM.02741-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Center for Ecology and Evolutionary Biology, 5289 University of Oregon, Eugene, Oregon 97403-5289,1 Department of Chemistry and Biochemistry, California Polytechnic State University, San Luis Obispo, California 93407,2 Thermal Biology Institute and Department of Land Resources and Environmental Sciences, Montana State University, Bozeman, Montana 597173
Received 5 December 2007/ Accepted 2 March 2008
|
|
|---|
0.2 to 4.0) and moderately high temperatures of 40 to 56°C. These unicellular algae occur in far-flung volcanic areas throughout the earth. Three genera (Cyanidium, Galdieria, and Cyanidioschyzon) are recognized. The phylogenetic diversity of culture isolates of the Cyanidiales from habitats throughout Yellowstone National Park (YNP), three areas in Japan, and seven regions in New Zealand was examined by using the chloroplast RuBisCO large subunit gene (rbcL) and the 18S rRNA gene. Based on the nucleotide sequences of both genes, the YNP isolates fall into two groups, one with high identity to Galdieria sulphuraria (type II) and another that is by far the most common and extensively distributed Yellowstone type (type IA). The latter is a spherical, walled cell that reproduces by internal divisions, with a subsequent release of smaller daughter cells. This type, nevertheless, shows a 99 to 100% identity to Cyanidioschyzon merolae (type IB), which lacks a wall, divides by "fission"-like cytokinesis into two daughter cells, and has less than 5% of the cell volume of type IA. The evolutionary and taxonomic ramifications of this disparity are discussed. Although the 18S rRNA and rbcL genes did not reveal diversity among the numerous isolates of type IA, chloroplast short sequence repeats did show some variation by location within YNP. In contrast, Japanese and New Zealand strains showed considerable diversity when we examined only the sequences of 18S and rbcL genes. Most exhibited identities closer to Galdieria maxima than to other strains, but these identities were commonly as low as 91 to 93%. Some of these Japanese and New Zealand strains probably represent undescribed species that diverged after long-term geographic isolation. |
|
|---|
4 and the number of species reaching levels below pH 5 is small (37). The order Cyanidiales consists of three recognized genera, Cyanidium, Galdieria, and Cyanidioschyzon (6, 13, 16), referred to colloquially as "cyanidia" in this work. The genera Cyanidium and Cyanidioschyzon are thought to include a single species each, Cyanidium caldarium and Cyanidioschyzon merolae, respectively, whereas the genus Galdieria has been classified into four species, Galdieria sulphuraria, Galdieria maxima, Galdieria partita, and Galdieria daedala, based on cell morphology. The last three species have recently been considered to be strains of G. sulphuraria, an opinion based on the similarities of sequences of the gene encoding the large subunit of ribulose bisphosphate carboxylase (rbcL) (7). The three genera have been distinguished morphologically, based on reproductive patterns, chloroplast and mitochondrion shapes, the presence or absence of a vacuole, cell size, and the presence or absence of a cell wall. Strains of the genus Cyanidioschyzon are smaller (usually 1 or 2 µm in breadth), with an oval, club-like shape. They reproduce by binary cytokinesis and lack rigid cell walls. Galdieria and Cyanidium have more spherical shapes, possess rigid cell walls, and reproduce through the formation of 4 to 32 small daughter cells within a mother cell. Chloroplasts vary from a multilobed shape in G. sulphuraria to a spherical shape in Cyanidium. Variation is also seen in mitochondrial shape. Some types may also be distinguished by physiological characteristics. G. sulphuraria has the unique ability to grow heterotrophically on at least 50 different carbon sources, including several sugars, sugar alcohols, amino acids, and trichloroacetic acid cycle intermediates (2, 15, 26). This character has been viewed as a plausible explanation for its ability to grow in darkness or semidarkness in soil, gravel, and as endolithic populations, presumably by using organic compounds released from other microorganisms in the community (10, 15, 26). Studies have also shown that G. maxima also can grow heterotrophically, although poorly (15). Up to this point, only G. sulphuraria lacks the ability to use nitrate as the sole nitrogen source (12).
Thermoacidic environments are scattered disparately throughout the earth, but cluster mainly in volcanic areas that have constant geothermal activity and low pH levels, due primarily to the presence of sulfuric acid as a result of the biological and abiological oxidation of elemental sulfur or hydrogen sulfide. Since the Cyanidiales tested thus far cannot tolerate desiccation, the geographical isolation of thermoacidic sites may have led to allopatric speciation events within this group over time (3, 4, 13). One geothermal region with many acidic habitats is Yellowstone National Park (YNP), an approximately 9,000-km2 area with numerous hydrothermal locations that vary in pH, temperature, exposure to high solar radiation, water availability, and concentrations of soluble metals and metalloids (25). This wide spectrum of physicochemical factors could exert selection pressures, potentially providing opportunities for distinct ecotypes to arise even in island-like features within YNP. Indeed, evidence of YNP population "endemism" or habitat and temperature specialization has previously been found for the thermophilic cyanobacterium Mastigocladus within YNP (22). The presence of the archaean Sulfolobus (25), as well as temperature specialization within YNP for the thermophilic cyanobacterium and other observations, establishes precedence for our efforts summarized herein, where both genetic and morphological criteria were used to examine the diversity of a large number of axenic cyanidial isolates obtained from many disparate sites within YNP and from larger spatial scales that include Japan and New Zealand. The data are consistent with the suggestion that geographical isolation has led to speciation within the Cyanidiales on a global scale and to the evolution of site-specific ecotypes within YNP.
|
|
|---|
![]() View larger version (36K): [in a new window] |
FIG. 1. Map of YNP, with numbers indicating from which sampling site the majority of cyanidial cultures were isolated. These sites are also indicated by number in Tables 1, 5, and 6. Site numbers and sample sites are as follows: 1, Lemonade Creek; 2, Obsidian Creek; 3, Nymph Creek; 4, Norris Geyser Basin; 5, Sylvan Springs; 6, Artists Paint Pots; 7, River Group; 8, Imperial Spring and Geyser and Fairy Falls Trailhead; 9, Rabbit Creek source; 10, Heart Lake Trail; 11, Mud Volcano; 12, Highland Creek; 13, Sour Creek, Clear Lake area. The distance "as the crow flies" is approximately 65 km from Mammoth to Grant Village.
|
|
View this table: [in a new window] |
TABLE 1. Global positioning system coordinates for all sampling sites in Yellowstone National Park, Japan, and New Zealand
|
|
View this table: [in a new window] |
TABLE 2. Identification, location, location code, environment, temperature, pH, and GenBank accession number of each strain used in this study
|
PCR amplification.
All primers for PCR studies are listed in Table 3. The 18S rRNA gene was amplified from 37 culture isolates from YNP, New Zealand, and Japan. The PCR mixture contained 50 ng of template DNA, 200 µM deoxynucleoside triphosphates (dNTPs), 1x PCR buffer, 2.5 mM MgCl2, 25 pmol primers CdmF and CdmR, and 5 U of Taq polymerase (Promega, Madison, WI). Cycling conditions were 95°C for 5 min, followed by 30 cycles of 94°C for 1 min, 60°C for 1 min, and 74°C for 1 min, and a final extension at 74°C for 10 min. PCR amplification of the rbcL gene was performed on 73 isolates. The reaction mixture contained 50 ng DNA, 200 µM dNTPs, 1x PCR buffer, 2.5 mM MgCl2, 25 pmol rbcL-1F and rbcL-1R, and 2.5 U Taq polymerase (Promega, Madison, WI). PCR products for both rbcL and 18S ribosomal DNA (rDNA) were sequenced bidirectionally by Macrogen, Inc. (Seoul, South Korea).
|
View this table: [in a new window] |
TABLE 3. Primers used in this study
|
Phylogenetic analyses.
All sequences for the 18S rDNA (487 nucleotides), rbcL (516 nucleotides), and cpSSR clones were edited using Sequencher (Gene Codes Corporation, Ann Arbor, MI) and aligned using ClustalX (www-igbmc.u-strasbg.fr/BioInfo/ClustalX/Top.html). Phylogenetic analysis for the rbcL, 18S rDNA, and cpSSR data sets was carried out using maximum parsimony methods with PAUP (http://paup.csit.fsu.edu/about.html). For maximum parsimony methods, 1,000 bootstrap replicates were included in a heuristic search, starting with a random tree and the tree bisection-reconnection branch-swapping algorithm. Percent variation calculations were made by comparing all isolates to the nearest BLAST relative. For Bayesian analysis (MrBayes, version 3.0b4), Metropolis-coupled Markov chain Monte Carlo analysis was performed with a random starting tree and four chains, of which three were hot and one was cold. Trees were sampled every 100 generations for 500,000 total generations by using a general time-reversible model and a gamma distribution.
Accession numbers for the 18S rDNA genes of previously sequenced rhodophytes used here for phylogenetic comparisons are AB090849 for G. sulphuraria L1B, AB090837 for G. sulphuraria A3D, AF441363 for G. sulphuraria ISG, AB090835 for G. sulphuraria FF1B, AF441359 for G. sulphuraria RT-27, AF441371 for G. sulphuraria MSh, AF441360 for G. sulphuraria AZ, AB091229 for G. sulphuraria 107.79, AF441361 for G. sulphuraria 136, AF441375 for G. sulphuraria ELS, AB090829 for G. sulphuraria FF3A, AB090851 for G. sulphuraria 86C, AF441373 for G. sulphuraria J, AF441365 for G. sulphuraria 21.92, AB090840 for G. sulphuraria 1C, AB090832 for G. maxima IPPAS P507, AB090844 for G. partita J5D, AF441372 for G. partita K, AB090839 for G. daedala IPPAS P508, AF441374 for C. caldarium 182, AF441366 for Cyanidium sp., AB158485 for C. merolae 10D, AF441376 for C. merolae 199, AB090833 for C. merolae 61D, AF168623 for Porphyridium aerugineum, and D88387 for Bangia atropurpurea.
Accession numbers for the rbcL gene sequences used for phylogenetic comparisons are AY541305 for G. sulphuraria DBV015NAFG, AY119768 for G. sulphuraria DBV009VTNE, AY541307 for G. sulphuraria DBV021MEVU, AY119769 for G. sulphuraria DBV063AGCS, AY119767 for G. sulphuraria SAG108.79, AY541302 for G. daedala IPPAS P508, AB018008 for G. partita, AY391370 for G. maxima IPPAS P507, AY119766 for C. merolae DBV001NAPS, AY541296 for C. merolae DBV202NAMN, AY541299 for C. merolae DBV182APAS, AY541298 for C. merolae DBV182JAVA, D63675 for C. merolae 10D, AY541297 for C. caldarium DBV019SIPE, AY541301 for Cyanidium sp. Monte Rotaro 20, AY541300 for Cyanidium sp. Monte Rotaro 19, AY119771 for Bangia fuscopurpurea, and AY119772 for Bangiopsis subsimplex. All accession numbers for the isolates used in the present study are listed in Table 2.
|
|
|---|
![]() ![]() View larger version (57K): [in a new window] |
FIG. 2. Maximum parsimony phylogenetic analysis of the cyanidial cultures described in this study. (A) Tree based on rbcL gene sequences using PAUP. The double asterisks represent those cultures with morphologies similar to that of C. merolae. Colors distinguish the various types. The reference strains (e.g., Cyanidium sp. Monte Rotaro 19 and Monte Rotaro 20 and most of the G. sulphuraria clusters [DVB 015, 021, 009, and 063]) may be found in reference 6. (B) Tree based on 18S rDNA gene sequences using PAUP. For both trees, branch lengths are proportional to the number of changes on a given branch, and bootstraps are given for each node above 50%. The double asterisks represent those cultures with morphologies similar to that of C. merolae. Colors of the various cultures correspond to the types indicated in panel A.
|
|
View this table: [in a new window] |
TABLE 4. Nearest BLAST relative for each type of isolate from YNP, Japan, and New Zealand for both rbcL and 18S rDNAa
|
An analysis of the rbcL sequences also showed two major clades among the isolates from Japan, with their closest identities to G. maxima IPPAS P507 (Fig. 2A), which was originally isolated from Kunashir Island of the Kurile Group of Russia near the northeastern coast of Japan. Clade type IIIA cultures were isolated exclusively from Owakudani, a volcanic area in Hakone National Park, Japan, and showed 99% rbcL sequence identity to the G. maxima IPPAS P507 strain, with little intragroup variation (Table 4). The second clade (type IIIB), with only 96% identity to G. maxima IPPAS P507, was made up of isolates from the Kusatsu and Nakabusa hot spring sites in Japan (Table 4). We also found a single isolate from Kusatsu (CCMEE 5680) sharing only 93% identity to G. maxima IPPAS P507 (Fig. 2A). In contrast to the differences found among the rbcL sequences and their geographically distinct locations, the 18S rDNA sequences varied little among all of the Japanese isolates; in a comparison of 18S rRNA sequences, all 24 Japanese strains were nearly 99% identical to the G. maxima strain IPPAS P507 (Fig. 2B).
The New Zealand strains fell into two main groups based on rbcL sequences. Members of the first group, type IV (Fig. 2A and Table 4), were 96% identical to G. sulphuraria SAG 108.79, whereas the closest database match for the type V isolates was G. maxima IPPAS P507 (Fig. 2A and Table 4), albeit at only 93% identity. Intragroup variation ranged from 0.8% to 1.0% (Table 4). There was also a single isolate, referred to as type VI, that was only 91% identical to G. maxima IPPAS P507, its closest match (Fig. 2A). By sequencing the 18S rDNA for a few selected New Zealand isolates, we established that they were all 99% identical to G. maxima strains IPPAS P507 and 61D and did not show the variation revealed by rbcL sequences (Fig. 2B). The 18S sequence is more highly conserved than that of the rbcL gene.
Morphology.
The isolates varied in cell morphology. In some cases, the morphological characteristics corresponded to separate clades identified in the rbcL-based phylogenetic analysis. For example, the G. maxima-like Japanese type IIIA isolates and the sister clades comprising the New Zealand type V and Japanese type IIIB isolates (Fig. 2A) were fairly similar in shape and size at maturity (Fig. 3D, E, and H). Likewise, the G. sulphuraria-like YNP type II and the G. sulphuraria-like New Zealand type IV isolates were also similar in both characteristics (Fig. 3C and G). In a similar fashion, different rbcL-based clades also reflected clear differences in cell size. For example, compare the G. sulphuraria-like New Zealand type IV with the New Zealand type V (Fig. 2A and 3G and H). As another example, the Japanese type IIIC cell was somewhat oval (Fig. 3F), distinguishing it from all other organisms isolated in this study, and this concurs with the fact that its rbcL identity match with any strain in the database was poor (Table 4). It should be remembered, however, that Galdieria-like cells grow to a somewhat indefinite size (depending on culture conditions) before giving rise (by an internal cell division) to 4 to 32 small daughter cells that then grow to a larger size. Therefore, the selected morphotypes (as shown in Fig. 3) cannot be totally relied on to show stable characteristics of each type.
![]() View larger version (77K): [in a new window] |
FIG. 3. The photomicrographs are selected examples of cultures of the various types included in Fig. 1A. Only type 1B illustrates the morphology of a typical C. merolae cell. The spherical cells in different images illustrate the typical morphology of Galdieria. However, the images of Japan types IIIB and IIIC and New Zealand (NZ) type V indicate cells that are somewhat oblong. The type of division for these three strains remains undetermined. CCMEE culture numbers for the types indicated in the panel are as follows: 5508 (A), 5610 (B), 5573 (C), 5665 (D), 5664 (E), 5680 (F), 5706 (G), 5704 (H), and 5713 (I).
|
cpSSR typing.
Recently, we found that cpSSRs (or microsatellites) could be used to identify genetic differences among six cultured Galdieria-like type IA morphotypes that shared identical 18S rRNA gene sequences with C. merolae (17). Therefore, in the current study, we employed an approach similar to the above method to more closely examine these strains for evidence of genotypic variation. Variation in cpSSR repeat number translated into different amplicon lengths that could easily be distinguished by agarose gel electrophoresis, resulting in a classification of type IA isolates into five major variants based on the number of repeats found in one cpSSR, previously referred to as cpSSR9 (17) (Table 5). Variant 1 of cpSSR9 contained 2.4 repeats of a 136-bp repeat unit sequence, variant 2 contained 3.4 repeats, variant 3 contained 4.4 repeats, and variant 4 had 5.4 repeats. The DNA from some type IA isolates did not amplify with the SSR9 primer set (but would amplify 18S rDNA and rbcL, as described above), and so for the purposes of comparison here, the isolates are referred to variant 5 (Table 5). When summarized as a function of distribution within YNP, cpSSR9 variant 2 isolates were found at all sampled sites except at the Imperial Spring steam vent (Fig. 1 and Table 6). SSR9 variants 3 and 5 were both isolated from nearly half of the sites, although they occurred together in only one, that being in Artists Paint Pots (Fig. 1 and Table 6).
|
View this table: [in a new window] |
TABLE 5. Selected strains from YNP types IA and IB categorized into SSR variants and their locations in YNP
|
|
View this table: [in a new window] |
TABLE 6. Distribution of SSR variants from sites throughout the park
|
![]() View larger version (9K): [in a new window] |
FIG. 4. Unrooted maximum parsimony trees of SSR9 variant 2 (A) and variant 3 (B) of selected YNP type IA and IIIA culture isolates. Bootstrap values are given for all nodes exceeding 50%. Branch lengths are proportional to the number of changes on a given branch.
|
|
|
|---|
Genotypic comparisons within YNP.
Comparisons of 18S rDNA and rbcL sequences from the YNP isolates were almost boringly predictable, falling into only two types. The type II YNP isolates were closely related to G. sulphuraria based on both rbcL and 18S rDNA analyses (Fig. 2A and B). Indeed, our isolates grew heterotrophically in the dark on glucose and were not able to use nitrate as the sole nitrogen source (data not shown), traits that are consistent with the G. sulphuraria strains isolated near Naples, Italy, and elsewhere (12). A second and much more intriguing group of YNP isolates was that for the organisms we refer to as type I, which comprised 95% of our more than 120 culture isolates from YNP and were made up of two different subtypes (IA and IB). Both of these classes have more than 99% rbcL and 18S rDNA sequence identity to C. merolae. However, only four of these isolates (type IB) have morphologies consistent with that currently recognized for C. merolae (20). All of our type I strains (type IA) have morphologies and cell division characteristics consistent with those of Galdieria (Fig. 3, compare panels A and B). The isolation of this type of organism in the current study serves to verify the previous reports of similar identity uncertainty by Ferris et al. (10) and Lehr et al. (17). More importantly, the extensive isolation of type IA morphophylotypes from throughout YNP in the current study also serves to demonstrate that these previous reports do not represent fortuitous observations. Rather, this Galdieria-like organism appears to represent the predominant member of the Cyanidiales within YNP. There is precedence for rDNA sequences commonly used for species determination to be the same in two distinctly different phenotypes. In the case of two photosynthetic dinoflagellates (18), the sequenced DNA fragments (internal transcribed spacers 1 and 2, 5.8S, the large subunit [D1/D2], and the small subunit) were identical in both organisms, but numerous stable morphological, physiological, and ecological characteristics of these two organisms were quite unlike.
As the rbcL and 18S rDNA phylotypes were found to be inadequate for describing genotypic diversity among the type IA isolates from YNP, we turned to cpSSRs in an attempt to identify evidence of some genetic diversity among these organisms. SSRs have been useful for resolving subspecies- or ecotype-level relationships in other eukaryotic organisms (29) as well as strain-level diversity among bacteria (8, 34). As in a previous study using cpSSRs for the purpose of seasonal monitoring of type IA cyanidial populations in a single YNP spring (17), we also showed the utility of SSRs for differentiating a small number of type IA culture isolates. In the current study, we were able to further separate YNP type IA isolates based on differences in cpSSR9 frequencies and sequences. An assessment of isolation sites of these variants revealed that there was an almost complete lack of cohabitation between variants 3 and 5 (Table 6), although there was nothing obvious in the within-YNP geographical distribution patterns for these particular cpSSR9 variants that would suggest migration constraints within YNP (compare Table 6 and Fig. 4 with Fig. 1). While at this juncture there is no information that directly links cpSSR sequence variation with important adaptive traits, it is possible that such clustering and the apparent separation of variants 3 and 5 may correspond to variations in the physicochemical characteristics that may occur among acid thermal environments in YNP (1, 4). For example, some habitats in the Norris Geyser Basin contain elevated levels of mercury (28) and we have measured exceptionally high total mercury concentrations in the soils at the extreme site B (250 µg·g–1 soil [analysis provided by Ed Nater, University of Minnesota]). In assessing the otherwise limited information available for the study sites, there was no apparent relationship between cpSSR9 variant types and temperature measured at sampling times (data not shown). Site distribution separation for variants 3 and 5 did not appear to be related to pH, although the variants were both more frequently encountered at sites where pH values were >2.0. There was also some evidence of pH effect for some of the other SSR9 variants. SSR9 variant 2 was found at all sites in YNP and at all pHs, although it was a more frequently encountered isolate (relative to the other cpSSR9 variants) as site pH values increased from 2.0 to 4.0.
Numerous cultures from samples collected at most sites and all cpSSR9 types were capable of growing well on the isolation medium. Nevertheless, a cultivation bias might fail to detect other SSR9 variants and it is possible that not all SSR variants were cultured from each sampling site. Thus, the apparent distribution patterns summarized in Table 6 may merely be coincidental and without additional genetic and physiologic information it is not yet possible to functionally associate SSR variants with environmental conditions. In the future, we hope to use the cpSSR9 locus, other appropriate loci, more detailed physiological characterization, and detailed chemical analyses of the sampling sites to further investigate possible "speciation" events that may have given rise to ecotypes within type IA in YNP.
A recent study of acid hot springs in Italy showed a high level of rbcL variation in the clone libraries and cultures from the Pisciarelli area within the Phlegrean Field complex and other areas of Italy (6). Two different clades of G. sulphuraria were found, one corresponding to isolates from endolithic sites and one from aquatic environments. In addition to G. sulphuraria, clone libraries and culturing revealed C. caldarium type organisms as well as C. merolae types. Some of the Cyanidium phylotypes were nonthermophilic forms from the Monte Rotaro area of Italy (Fig. 2A). The rbcL data from the Italian study support the biogeographic and habitat niche separation of different cyanidial clades. The Italian cyanidial lineages, however, are distinct from the Japan, New Zealand, and Yellowstone rbcL clades.
A recent report describing photosynthetic organisms in a single acidic siliceous crust in the Norris Geyser Basin (YNP) demonstrated that the majority of the 16S rDNA PCR clones exhibited their closest identities to the C. caldarium chloroplast phylotype (98 or 99%) (36). However, the clone libraries did not detect the presence of Galdieria sp. or C. merolae in this one endolithic sample. In contrast, we did not detect any C. caldarium phylotypes through extensive culturing attempts from collections in YNP, Japan, and New Zealand, nor has C. caldarium been detected in 18S rDNA molecular surveys of Dragon Spring in Norris Geyser Basin (17) or in Nymph Creek (10), both in YNP. The paper by Walker et al. (36) claims to have identified a cyanidial sequence closest to C. caldarium from an endolithic sample in the Norris Geyser Basin, but we believe this claim is in error (unpublished data). It is possible that no C. caldarium phylotype exists in YNP. In earlier studies by Brock (4), all cyanidia were referred to as C. caldarium. However, it is now clear from those results, which included a demonstration of dark heterotrophic growth, that G. sulphuraria was the organism being characterized.
Evolutionary hypotheses and questions.
Seckbach (32, 33) has hypothesized an evolutionary order from C. merolae (the simplest morphologically) to C. caldarium to G. sulphuraria by comparing morphological and physiological characteristics. However, a recent study comparing the complete genome of C. merolae strain 10D with expressed sequence tags and high-throughput genomic sequence reads of G. sulphuraria (covering 70% of the genome of G. sulphuraria strain 074W) suggests an alternative hypothesis (2). It appears that approximately 30% of the G. sulphuraria genome has no relation to that of C. merolae. Both organisms have similar varieties of sugar kinases, but only G. sulphuraria has a high number of monosaccharide transporters. The lack of transport proteins helps explain the obligate photoautotrophy of C. merolae, and the presence of such proteins in G. sulphuraria explains its ability to grow on numerous carbon sources. The same genomic comparison (2) also established that C. merolae lacks genes coding for enzymes involved in the synthesis, modification, and degradation of the cell wall and the synthesis of red algal floridosides (for osmoregulation) and also displayed less cellular compartmentalization. Others have suggested that the cell wall architecture of G. sulphuraria may allow it to tolerate drier environments, while C. merolae may be found mainly in aquatic niches that are osmotically stable and, thus, do not require a rigid wall (2, 6). These cell structural differences suggest that the cyanidia may have evolved from a common ancestor possessing a cell wall and heterotrophic capability as in G. sulphuraria (2). In such an evolutionary scenario, G. sulphuraria would have retained the ability to grow heterotrophically in moist, darker sites (e.g., moist acid soils and gravels), where it would more likely encounter organic carbon sources. The ancestor of C. merolae may have colonized osmotically stable, inorganic habitats and lost the ability to synthesize a cell wall or use organic compounds (2).
Potentially, the YNP type IA isolates represent an interesting piece of this evolutionary puzzle. With their cell walls and spherical shapes, they resemble a Galdieria-type organism, yet genetically, they are nearly identical to C. merolae, at least when viewed using the rbcL and 18S loci. Our culture isolates of both type IA and type IB are from acid streams as well as from moist soil environments. In preliminary screens of a sampling of these organisms, none grew heterotrophically on glucose, but otherwise they varied considerably with respect to growth rate, ability to use nitrate as a sole nitrogen source, and their cardinal temperatures (data not shown). And while such differences are certainly noteworthy with respect to the characterization of these organisms, the extent of their genetic differences is as yet unknown and, thus, it would be difficult to reclassify these various strains into separate ecotypes or other taxonomic entities at this juncture.
The finding of this abundant type IA organism in YNP (a possible intermediate between a typical spherical G. sulphuraria-type cell and a typical Cyanidioschyzon) may help in understanding the evolutionary events that may have allowed these two types of organisms to diverge into different habitats. Natural genetic exchange is so far unknown in the Cyanidiales. If this lack is, in fact, true, the evolutionary history of these organisms may be easier to unravel. Minoda et al. (23) have successfully accomplished limited gene transformation in C. merolae through electroporation, presumably facilitated by the lack of a cell wall in this member of the Cyanidiales. According to Yoon et al. (39), the cyanidial group probably evolved as the earliest branch of the red algae in the Precambrian era, earlier than 1.3 billion years before the present. A few members of the Cyanidiales live in nonthermal acidic habitats (14, 31). It is therefore possible that primordial cyanidia evolved into thermophilic types later, through periodic selection and adaptive divergence, thus escaping from the competition within a less extreme acidic environment. The same evolutionary path may have been followed even earlier, i.e., an escape from extensive competition in a neutral environment to a nonthermal acidic environment where, even today, relatively few species of algae occur.
Global distribution of cyanidia.
On a global scale, phylogenetic analyses based on the rbcL gene identified two different Japanese clades most closely related to G. maxima, which was originally isolated from the Kurile Islands off the northeast coast of Hokkaido (Fig. 2A). However, when we compared these organisms based upon their 18S rDNA sequence data, the clade structures partially collapsed (Fig. 2B). A similar discrepancy between rbcL and 18S phylogenies was observed with the isolates from New Zealand. Organisms more closely clustered with G. sulphuraria, based on rbcL comparisons (e.g., New Zealand type IV [5722 WI and 5706 CR]), were more closely related to G. maxima when viewed from the 18S rDNA perspective (Fig. 2, compare panels A and B). The discrepancy between the chloroplast rbcL and nuclear 18S rDNA sequence data may indicate a difference in the time scale of chloroplast evolutionary changes (or at least of the rbcL gene) within the cyanidia of Japan and New Zealand that are independent of the nuclear DNA as characterized by 18S rDNA, which shows little divergence. There is precedence for rates of chloroplast DNA evolution differing from those of nuclear DNA (11, 35). This precedence also hints at lateral transfer of the rbcL gene. The ancestor of the New Zealand type IV strains could have obtained its plastid rbcL gene from a strain somewhat related to G. sulphuraria, while the more conserved, nuclear 18S rDNA sequence indicates phylogenetic identity closest to G. maxima (Table 4). Although no natural mechanism of gene transfer is known in the Cyanidiales, this does not preclude the possibility of its existence (9). The low rbcL sequence identity for both the Japanese type IIIB (96%) and type IIIC (94%) to G. maxima strain IPPAS P507 suggests the possibility of many more ecotypes or species in existence that are somewhat distantly or closely related to G. maxima than were previously known (Table 4). Some of these new taxa may be distinguished by differences in plastid characteristics rather than by nuclear or cytoplasmic dissimilarities. We anticipate further characterization of these isolates by using physiological differences and the sequences of more variable genes and proteins.
The geographically specific rbcL-based and even 18S-based clades of G. sulphuraria and G. maxima lineages suggest the possibility that unique cyanidial types evolved during long-term isolation over great geographical distances between western North America, Japan, and New Zealand (19). Because of the distances and the great ages of the volcanic hot springs in all three regions and because of the fidelity of specific 18S and rbcL genotypes to only one region, allopatric divergence and speciation must be regarded as the most likely basis of the genetic differences among these geographical provinces. This might not be the case if it were not for the fact that the ranges in chemistry of the acidic habitats sampled in all three areas are very similar. It is known from the literature and from our tests (unpublished data) that cyanidia do not tolerate desiccation, and therefore dispersal of viable founder populations over great distances must be a very rare event (38). The most probable vector may be migratory birds that carry moist acidic mud on their feet or pass viable cells through their acidic digestive systems. The latter case has been demonstrated for unicellular green algae, diatoms, cyanobacteria, and others, even in species not inhabiting acid waters and without the aid of resting stages (30).
Even without multilocus analyses, the divergences of the rbcL and 18S genes suggest speciation events that may have taken place over thousands or millions of years. A caveat of the conclusions raised here is that more strains, species, or ecotypes probably exist in Japan and New Zealand and collecting and isolation attempts have not reached saturation. However, we are probably close to saturation in Yellowstone, since cyanidia are readily cultivated and isolates from almost all acidic areas in the park are in the culture collection and have been sequenced using the rbcL and 18S loci.
Our study of the cyanidia of acid hot springs in YNP, Japan, and New Zealand shows that we are only beginning to understand the diversity of this interesting group of primitive red algae in an environment quite atypical of those of all other members of this phylum.
We thank John Varley and Christie Hendrix at the Yellowstone Center for Resources for general assistance and for obtaining research permits. We also thank Dana J. Skorupa for valuable assistance as well as the three unidentified reviewers of the original manuscript for their constructive input.
Published ahead of print on 14 March 2008. ![]()
|
|
|---|
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»