DNA fingerprinting differentiation between β-carotene hyperproducer strains of Dunaliella from around the world
© Olmos et al; licensee BioMed Central Ltd. 2009
Received: 25 September 2008
Accepted: 30 June 2009
Published: 30 June 2009
Dunaliella salina is the most important species of the genus for β-carotene production. Several investigations have demonstrated that D. salina produces more than 10% dry weight of pigment and that the species grows in salt saturated lagoons. High plasticity in the green stage and the almost indistinguishable differences in the red phase make identification and differentiation of species and ecotypes very difficult and time consuming.
In this work, we applied our intron-sizing method to compare the 18S rDNA fingerprint between D. salina (CCAP 19/18), D. salina/bardawil (UTEX LB2538) and β-carotene hyperproducing strains of Dunaliella isolated from salt saturated lagoons in Baja, Mexico. All hyperproducer strains reached β-carotene levels of about 10 pg/cell. Optical microscopy did not allow to differentiate between these Dunaliella strains; however, 18S rDNA fingerprinting methodology allowed us to differentiate D. salina from D. salina/bardawil.
In Baja Mexico we found D. salina and D. salina/bardawil species by using intron-sizing-method. The National Center for Biotechnology Information (NCBI) Dunaliella 18S rDNA gene sequences were analyzed with our methodology and extraordinary correlation was found with experimental results.
Dunaliella was originally described by Teodoresco 1905. Since then, taxonomic studies among Dunaliella have identified several new species. However, even today differentiation among halophilic and carotenogenic Dunaliella species in both green and red stages is difficult and time consuming. In addition, some strains and species in culture collections are misnamed and have given rise to unnecessary strains and species names . This confusion of strains and species names makes comparison of results by different authors difficult . D. salina and D. salina/bardawil are the only reported β-carotene hyperproducer species of the genus (at least 10%), that grows in salt saturated lagoons . However, a controversy still exists about identification of D. salina/bardawil as different species of D. salina . Molecular identification provides a useful tool to distinguish between inter and intra-specific morphologically similar species [4, 5] and mixed populations [6, 7]. Intra-species identification of community members without cultivation, avoid some selective biases associated with pure culture methods . Specie-specific oligonucleotides could be useful to identify either species from culture collections or from natural environments [9, 10]. Slight phylogenetic and taxonomic differences in Dunaliella species can conceal profound differences in their potential for production of metabolites such as carotenoids . To study these organisms one must discriminate them from their presumed non carotenogenic congeners and from among species with similar morphological features .
Results and discussion
Taking into account that D. salina strains isolated from different continents presented the same fingerprinting profile and β-carotene levels, we can be assured that the intron-sizing method developed by our group is trustworthy, rapid and sensitive, for specifically identifying these hyperproducer strains of D. salina. We can conclude that D. salina strains with the following fingerprinting profile of MA1–MA2 = ~2100 bp and DSs-MA2 = ~700 bp, belong to the same species, are distributed worldwide and their β-carotene hyperproduction capacity is well conserved. Thus, in order to facilitate identification, we suggest a sub-classification calling these strains as "D. salina var Teod".
Recently, the National Center for Biotechnology Information (NCBI) has initiated a D. salina genome sequencing project (ID: 18607). In this sense, it is important to know specific identity of target species, particularly if several groups from several countries are involved in this ambitious and time consuming task. Confidence knowing specific identity of target species can avoid wasting time and money.
Additionally, demands for natural products are growing in importance and β-carotene is not an exception [1, 15]. Rapid and specific identification and certification of commercial strains or species can improve production yields and help in the competitiveness by avoiding contamination with undesired species. Cross-contamination even with other hyperproducer strains still can affect the productivity due different environmental and nutritional conditions required for each strain. As an example, in May of 2008 we purchased LB2538 strain from UTEX again, but surprisingly, a fingerprinting profile of "D. salina var Teod" was obtained (MA1–MA2 = ~2100, DSs-MA2 = ~700) instead of the fingerprinting profile of D. salina/bardawil (MA1–MA2 = ~2500, DBs-MA2 = ~1000) which we were expecting (data not shown).
Moreover, contamination with non hyperproducer strains is more dangerous because these strains grow faster than hyperproducers, due to the fact that there are no metabolic efforts required for β-carotene production. This was the case for 19/18 samples where hyperproducer and non hyperproducer species of Dunaliella were identified by our fingerprinting methodology (Fig. 2).
D. salina 19/30 supposedly was the strain isolated by Dr. Ben-Amotz (D. salina/bardawil), however PCR amplification with MA1–MA2 presented neither the ~2500 bp product associated with the original D. salina/bardawil species , nor the ~2100 bp amplification associated with D. salina var Teod species. Nevertheless, the 19/30 strain amplified a PCR product of ~1700 bp using MA1–MA2 conserved primers (Fig. 3) and the species itself never turns red, indicating β-carotene is not hyperproduced by this strain. In addition, Gene Bank had reported a 19/30 strain with ~1700 bp 18S rDNA gene size (EF473749.1) and another 19/30 strain was reported with ~2100 bp 18S rDNA gene size (DQ447648.1). The first report comes from Italy and the second comes from France. This is another good example of cross-contamination that researchers and even certified cultured collection banks can undergo, due to almost indistinguishable taxonomic characteristics between Dunaliella species.
Auxiliary methods like the one we are proposing could be useful for characterizing and certifying international culture collections as well. Taxonomic characterization is important to identify and differentiate species; however, differentiation between subspecies, ecotypes or similar species with taxonomy tools is a difficult and time consuming task and sometimes impossible. For this reason, the 18S rDNA intron-sizing method proposed in this work provides a novel and powerful DNA-fingerprinting methodology to accomplish a specific, rapid and sensitive identification of carotenogenic Dunaliella species.
On the other hand, D. salina LB1644 strain obtained from UTEX did turn to red; however it produced only 1% dry weight of β-carotene. LB1644 strain amplified a ~1700 bp PCR product with MA1–MA2 conserved primers as reported by our group  and recently by gene bank (DQ009765.1), but did not amplify any band with DSs-MA2 or DBs-MA2. All results; 18S rDNA gene size, fingerprinting profile and pigment production capacity, corroborated that LB1644 strain it is not "D. salina var Teod" n or "D. salina/bardawil". In this case, more taxonomic and molecular research must be done to identify this Dunaliella specie. Absence of introns within 18S rDNA gene is an inconvenience for Dunaliella species fingerprinting determination and identification. Fortunately, all hyperproducer species known carry I or II introns within the 18S rDNA gene. Gómez and González , have been working with ITS sequences and some interesting results have been presented that could help in primer design. In addition, it seems to be the rule for hyperproducer species of Dunaliella to have introns within 18S rDNA, thus "D. salina var Teod" has I, D. salina/bardawil has II and D. parva has II .
However, we did find a non carotenogenic strain in hypersaline environments that has two introns within the 18S rDNA (data not shown). According to literature this strain could be D. viridis. Gene Bank sequence analysis showed a D. viridis strain (DQ009776) with an 18S rDNA gene size of ~2500 bp approximately, which corresponds to the size of the strain that we found. Species-specific primer design must be done to molecularly characterize this strain and, wheatear it is possible accomplish a fingerprint profile of it.
Dunaliellagene bank sequence analysis
Dunaliella strain sequences reported in the National Center for Biotechnology Information (NCBI) and utilized in this work.
Dunaliella parva LB 1983
Dead sea, Israel
Dunaliella bardawil LB 2538
Bardawil lagoon/San Quintin
Dunaliella sp. ABRII G2/1
Dunaliella viridis CONC002
Dunaliella bardawil LB 2538
Dunaliella peircei LB 2192
Lake Marina, Cal
Dunaliella salina LB 200
Dirty, salt lake
Dunaliella sp. SPMO 201-5
Salt flat, Oklahoma
Dunaliella salina CCAP 19/30
Di Giuseppe, G.
Dunaliella salina SAG 19-3
Di Giuseppe, G.
Dunaliella salina CCAP 19/3
Dirty salt lake
Di Giuseppe, G.
Dunaliella salina CCAP 19/18
Dunaliella sp. ABRII G2/1
Dunaliella sp. BSF2
Salt flat, Oklahoma
Dunaliella sp. SPMO 200-8
Salt flat, Oklahoma
Dunaliella sp. SPMO 300-5
Salt flat, Oklahoma
Dunaliella sp. SPMO 600-1
Salt flat, Oklahoma
Dunaliella sp. SPMO 601-1
Salt flat, Oklahoma
It would be very interesting to know its β-carotene production capacities, because this strain has its own fingerprinting profile, different from the ones presented by "D. salina var Teod" and D. salina/bardawil (data not shown).
Table 1, also shows two other groups of Dunaliella strains with an 18S rDNA gene size of ~2100 bp, indicating one intron within the 18S rDNA gene. D. salina M84320 strain from group I, reported by Wilcox and coworkers in 1992, was used to design species-specific oligonucleotide (DSs), due its β-carotene hyperproduction capacity . The 18S rDNA fingerprinting profile from M84320 was already found in D. salina strains from around the world, and to make its identification easier we call these strains "D. salina var Teod" (data mentioned before). Furthermore, D. salina var Teod strains hybridized 100% with DSs primer [5'-GCAGGAGAGCTAATAGGA-3'] and hyperproduce β-carotene, data that corroborates and supports our fingerprinting results (data mentioned before).
In addition, DQ009777 strain from group I was reported as D. bardawil, nevertheless this data disagrees with our results for D. salina/bardawil where an 18S rDNA gene size of ~2500 bp was found, instead of ~2100 bp as DQ009777 had. In addition, this strain did not amplify with the species-specific oligonucleotide (DBs) for D. salina/bardawil, but amplified well with DSs specific primers to D. salina var Teod. However, DQ009777 strain reported as D. bardawil by Buchheim is in agreement with the result that we obtained from D. bardawil (LB2538) strain obtained from UTEX in May of 2008. Nevertheless, results obtained from β-carotene hyperproducer strains isolated from hypersaline environments in Baja Mexico, demonstrate that D. salina/bardawil with ~2500 bp 18S rDNA gene size exist, amplify ~1000 bp PCR product with DBs-MA2 and have an exclusive fingerprinting profile (Fig. 3). These results are in agreement with fingerprinting profile found to D. salina/bardawil (LB2538) obtained ten years ago from UTEX. This comparison shows that D. salina/bardawill from now it is not the same than 10 years ago, in the culture collection bank.
On the other hand, even when D. salina (DQ324002) group II strains have ~2100 bp 18S rDNA gene size as group I had, they only hybridize 66% with DSs primer [5'-GCtAcGtGAGAGCTAAgTAGGaAg-3'], change are shown in subscript letters. Group II, therefore, belongs to a different subgroup of "D. salina var Teod" strains. Unfortunately we do not know if group II strains hyperproduce β-carotene. In a personal communication Dr. Buchheim mentioned that they also do not know about β-carotene production capacity of these strains. It will be important to obtain group II strains and test their β-carotene production capacity. At this moment we know that group II Dunaliella strains has a fingerprint of ~2100 bp with MA1–MA2 and ~750 bp with DSsII-MA2. DSsII (5-GAGAGCTAGCAGAGGGTAG-3) oligonucleotide was designed from a region 100% conserved among group II strains, using BLAST program from NCBI (data not shown).
Finally, Gene Bank BLAST analysis showed a group of Dunaliella strains with ~1700 bp 18S rDNA gene size (data not shown). However, these strains were not included for group comparison because they do not have introns within the 18S rDNA gene and do not produce β-carotene, with the exception of the LB1644 Dunaliella strain. It is important to remember that our fingerprint intron-sizing method depends exclusively on introns within the 18S rDNA gene.
PCR identification with 18S rDNAspecific primers
To corroborate the presence of the most biotechnologically important Dunaliella species in isolated samples, a standard PCR protocol was carried out with species-specific D. salina/bardawil (DBs-MA2) and "D. salina var Teod" (DSs-MA2) primers. Because introns of the 18S rDNA gene represent the highest sequence divergence in Dunaliella [11, 12], specific-primers were designed from them. As we expected, using DBs-MA2 and DSs-MA2 pairs of primers, two different DNA products were found in PCR reactions. Fragments around 1000 bp long were found in samples coming from San Quintín. These results corroborates D. salina/bardawil presence in red environmental samples coming from this place (Fig. 4) On the other hand, fragments around 700 bp long were found in samples coming from La Salina, corroborating the presence of "D. salina var Teod" in the sample mentioned above (Fig. 4). Results obtained with conserved and specific primers match and corroborate presences of D. salina/bardawil and "D. salina var Teod" in the environmental samples.
Sequencing 18S rDNAgenes
The PCR product of around 2500 bp coming from the San Quintín sample using MA1–MA2 conserved primers was sequenced and aligned. BLAST results indicated that D. salina/bardawil reported (AF150905), matched 99% with our sequence. In addition, a ~2100 bp long PCR product amplified from the sample coming from La Salina was sequenced and aligned. BLAST results corroborated "D. salina var Teod" reported (M84320) matched 99% with our sequence.
Cell counts gave around 1–3 × 104 cells/ml in samples analyzed using a Neubauer chamber. Salt saturation was evident due to crystal precipitation. Red pigmentation was observed with high intensity, indicating great amounts of β-carotene (data not shown). Temperatures differed according to locations, being cooler in the North part of the Peninsula (20°C), than in the South (22°C). β-carotene quantification of the molecularly identified Dunaliella strains was developed according to conditions mentioned in materials and methods. β-carotene production reached almost same levels in all Dunaliella strains analyzed (~10 pg/cell).
Table 1, shows three groups of Dunaliella strains. The first group belongs to Dunaliella with two introns and approximately 2500 bp 18S rDNA gene size. Four strains had species-specific sequence differentiation between them and as a consequence a species-specific fingerprinting profile for each one (data shown only for D. salina/bardawil). Table 1 also shows two groups of Dunaliella strains with approximately 2100 bp 18S rDNA gene size. According to our results strains with that size contain only one intron within the 18S rDNA gene, and are classified as D. salina strains. Group I D. salina strains (M84320 group) were called "D. salina var Teod", because all strains share 100% sequence conservation with DSs specific oligonucleotide, growth on hypersaline environments and hyperproduce β-carotene. In addition, these strains had been found around the world and all share the same fingerprinting profile (see results and discussion section).
A second subgroup of D. salina strains called group II (DQ324002 group), was identified from Gene Bank sequence analysis (data not shown). More work must be done to known their β-carotene production capacity. Additionally, even when strains of group II had the same 18S rDNA gene size, as group I did, they only maintain 66% sequence conservation with "D. salina var Teod" specific oligonucleotide (DSs).
On the other hand, D. salina/bardawil obtained ten years ago from UTEX it is not the same strain found presently in the culture collection bank. However, a D. salina/bardawil strain with the original fingerprinting profile was found in San Quintín, B.C. Mexico ten years ago and also during this work.
It is important to point out that although culture collection banks may have some cross-contamination, as do collections of some researchers around the world (results and discussion), it is not our intention to point out cross contamination or incorrect affiliation of Dunaliella species. Nevertheless, we want to provide an auxiliary method for rapid, precise and specific identification of Dunaliella species.
Since our first report in 2002 on differentiation of carotenogenic Dunaliella species by the introns sizing method , we have refined, corroborated and deeply advanced the molecular differentiation of Dunaliella species by the fingerprinting profile method. Actually, it is the only method that has been tested on both collection bank strains as well as field samples, that allows proper identification of the physiologically plastic strains of Dunaliella at the level of species.
Dunaliella salina CCAP 19/18 sample was obtained from Sams Research Services LTD, and it was originally isolated from Hypersaline Brines, Hutt Lagoon, Western Australia.
Dunaliella salina/bardawil was obtained from UTEX (LB2538) ten years ago. This strain was isolated from Bardawil lagoon. A strain with the same 18S rDNA fingerprinting profile and β-carotene production capacity was isolated from Baja California, Mexico .
D. salina 19/30 was obtained from sams algae collection. This strain was isolated from Bardawil lagoon. D. salina (LB1644) was obtained from UTEX.
Sample cells were used for morphological identification as Dunaliella strains and for molecular characterization as Dunaliella species. Cells were counted and harvested at 10,000 r.p.m. for 5 minutes. β-carotene production capacity of strains was evaluated by High Performance Liquid Chromatography (data not shown).
Two locations in Baja California peninsula of Mexico were selected for water-sample collection (Fig. 1). Environmental samples were obtained using 50 ml plastic tubes from red hypersaline waters between May and September 2006. Red microalgae strains present in those samples were morphologically identified as Dunaliella strains and then were molecularly characterized. Sample cells were counted and harvested at 10,000 r.p.m. for 5 minutes. β-carotene production capacity of strains was evaluated by HPLC.
DNA purification and PCR amplifications
Isolation of chromosomal DNA from Dunaliella strains was carried out by the method utilized for Escherichia coli . MA1 [5'-CGGGATCCGTAGTCATATGCTTGTCTC-3'] and MA2 [5-GGAATTCCTTCTGCAGGTTCACC-3'] conserved and DSs (D. salina) [5'-GCAGGAGAGCTAATAGGA-3'] and DBs (D. salina/bardawil) [5'-GGGAGTCTTTTTCCACCT-3'] specific oligonucleotides, were designed from 18S rDNA genes and were previously reported by Olmos and coworkers [5, 11]. In those works exons were used to design MA1 and MA2 conserved primers. In addition, DSs and DBs specific primers were designed from intron I, where most variable sequence is contained. PCR reactions were carried out in a total volume of 100 μl containing 50 ng of chromosomal DNA in TE (Tris-EDTA) buffer, pH 8  and 200 ng MA1 and MA2 conserved primers. The amplification was carried out using 25 cycles in a GeneE thermocycler, with a Tm of 52°C to all reactions. One cycle consisted of 1 minute at 95°C, 1 minute at 52°C and 2 minute at 72°C. PCRs with specific primers were combined with MA2 conserved primer (DSs-MA2, DBs-MA2) and reactions were carried out under same conditions.
Sequencing and alignment
MA1–MA2 PCR products from Dunaliella species were utilized to carry out a sequencing reactions, using previously reported Seq1 (5'-GGTTGATCCTGCCAGTAG-3'), Seq2 (5'-CCGGGCATTTTTGTCTGG-3'), Seq3 (5'-CTGCCAGCACCTTATGAG-3'), Seq4 (5'-GGGAGGATTGACAGATTG-3') and Seq5 (5'-GGAAGGAGAAGTCGTAAC-3') primers. DNA sequences were imported to BLAST for strains identification and to MegAlign program from DNAStar, to search for phylogenetic relationship correlations between them. Dunaliella Gene Bank sequences comparison analysis was carried out with BLAST program using 18S rDNA genes with ~2500, ~2100 and ~1700 bp.
β-carotene purification and HPLC quantification
One ml samples were centrifuged at 10,000 r.p.m. for 5 minutes, supernatant was discarded, and pellets were homogenized and disrupted in one ml of methanol. Samples were incubated on ice for one hour in the dark, centrifuged at 10,000 rpm 5 minutes, filtered and kept at -20°C until HPLC analysis.
Filtered samples and commercial standards were analyzed with an 1100 Hewlett Packard HPLC, using a 5 μm reverse phase 300SB-C18 Zorbax column 15 cm long, and 100% isocratic methanol elution with 0.5 ml/minute flux, throughout 25 minutes.
The authors wish to thanks Judy Swan for English language assistance and Abelardo Campos for pictures edition. This work was supported by CICESE project 682112.
- Oren A: A hundred years of Dunaliella research: 1905–2005. Saline systems. 2005, 1: 2-10.1186/1746-1448-1-2.View ArticleGoogle Scholar
- Borowitzka MA, Siva CJ: The taxonomy of the genus Dunaliella (Chlorophyta, Dunaliellales) with emphasis on the marine and halophilic species. J Appl Phycol. 2007, 19: 567-590. 10.1007/s10811-007-9171-x.View ArticleGoogle Scholar
- Ben-Amotz A, Asís A, Avron M: The biotechnology of cultivating Dunaliella for production of β-carotene rich algae. Bioresource Technol. 1991, 38: 233-235. 10.1016/0960-8524(91)90160-L.View ArticleGoogle Scholar
- Gómez PI, González MA: Genetic variation among sever strains of Dunaliella salina (Chlorophyta) with industrial potencial, based on RAPD bandung patterns and on nuclear ITS rDNA sequences. Aquaculture. 2004, 233: 149-162. 10.1016/j.aquaculture.2003.11.005.View ArticleGoogle Scholar
- Olmos SJ, Paniagua MJ, Contreras FR: Molecular identification of Dunaliella sp. utilizing the 18S rDNA gene. Lett Appl Microbiol. 2000, 30: 80-84. 10.1046/j.1472-765x.2000.00672.x.View ArticleGoogle Scholar
- García TA, Olmos SJ: Quantification by fluorescent in situ hybridization of bacteria associated with Litopenaeus vannamei larvae in Mexican shrimp hatchery. Aquaculture. 2007, 262: 211-218. 10.1016/j.aquaculture.2006.10.039.View ArticleGoogle Scholar
- Hernández ZG, Olmos SJ: Identification of bacterial diversity in the oyster Crassostrea gigas by fluorescent in situ hybridization and polymerase chain reaction. J Appl Microbiol. 2006, 100: 664-672. 10.1111/j.1365-2672.2005.02800.x.View ArticleGoogle Scholar
- DeLong EF, Wickham GS, Pace NR: Phylogenetic stains: ribosomal RNA-based probes for the identification of single cells. Science. 1989, 243: 1360-1363. 10.1126/science.2466341.View ArticleGoogle Scholar
- Ki SJ, Han MS: Sequence-based diagnostics and phylogenetic approach of uncultured freshwater dinoglagellate Peridinium (Dinophyceae) species, based on single-cell sequencing of rDNA. J Appl Phycol. 2005, 17: 147-153. 10.1007/s10811-005-7211-y.View ArticleGoogle Scholar
- Ki JS, Han MS: Rapid molecular identification of the harmful freshwater dinoglagellate Peridinium in varios life stages using genus-specific single-cell PCR. J Appl Phycol. 2007, 19: 467-470. 10.1007/s10811-007-9157-8.View ArticleGoogle Scholar
- Olmos SJ, Paniagua MJ, Contreras FR, Trujillo L: Molecular identification of β-carotene hyperproducing strains of Dunaliella from saline environments using species-specific oligonucleotides. Biotechnol Letters. 2002, 24: 365-369. 10.1023/A:1014516920887.View ArticleGoogle Scholar
- Wilcox LW, Lewis A, Fuerst PA, Floyd GL: Group I introns within the nuclear-encoded small-subunit rRNA gene of three green algae. Mol Biol Evol. 1992, 9: 1103-1118.Google Scholar
- Paniagua MJ, Capa RW, Olmos SJ, Gutierrez ML: The carotenogenesis pathway via the isoprenoid-β-carotene interference approach in a new strain of Dunaliella salina isolated from Baja California Mexico. Marine Drugs. 2009, 7: 45-56. 10.3390/md7010045.View ArticleGoogle Scholar
- Raja R, Hema IS, Balasubramanyam D, Rengasamy R: PCR-identification of Dunaliella salina (Volvocales, Chlorophyta) and its growth characteristics. Microbiol Res. 2007, 162 (2): 168-176. 10.1016/j.micres.2006.03.006.View ArticleGoogle Scholar
- Olaizola M: Comercial development of microalgal biotechnology: from the test tube to the marketplace. Biomol Eng. 2003, 20 (4-6): 459-466. 10.1016/S1389-0344(03)00076-5.View ArticleGoogle Scholar
- Sambrook KT, Frisch EF, Maniatis T: Molecular cloning: a laboratory manual. 1989, Cold Spring Harbor, Cold Spring Harbor PressGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.