BACKGROUND
dunaliella salina teodoresco  <cit>  is one of the best-studied unicellular green algae  <cit> .
this is not only because d. salina is halotolerant, thriving in extreme saline environments  <cit> , but also because it can produce large quantities of β-carotene  in lipid globules located within the chloroplast  <cit> .
these traits make d. salina a model organism for investigating the evolution of salt adaptation  <cit>  and an attractive "cell factory" for the commercial production of β-carotene  <cit> .
although a great deal is known about the physiology and biochemistry of d. salina  <cit> , very little is known about the genomes it carries, especially those within its organelles.
until now, nothing was known about the size, conformation, or gene complement of either the mitochondrial or plastid genomes of d. salina  even though the sequences of these genomes are essential to the development of new d. salina technologies, such as a viable plastid transformation system  <cit> .
research on green-algal organelle genomes has led to significant advancements in genetic engineering.
the first stable transformation of a plastid genome was achieved in  <dig> using the unicellular green alga chlamydomonas reinhardtii  <cit>  and, soon after, the first example of recombinant protein expression in a plastid was also achieved using c. reinhardtii  <cit> .
since then, many techniques for plastid engineering have been first developed for green algae and then adapted for use in land plants  <cit> .
given the relatively close evolutionary proximity of c. reinhardtii and d. salina  <cit> , it is reasonable to assume that many of the technologies for c. reinhardtii plastid transformation might be transferable to d. salina.
various groups have attempted to transform d. salina  <cit> ; however, a lack of plastid-genome sequence data has prevented successful plastid transformation.
therefore, the first step in developing an efficient and reliable plastid transformation system for d. salina is to sequence its organelle genomes.
d. salina is an attractive alga for organelle genome research and plastome engineering for a variety of reasons: i) various strains and geographical isolates of d. salina are readily available from algal culture collections around the world; ii) d. salina is relatively easy to grow and maintain -- it is one of the few microalgae that are being cultivated currently on a large scale; iii) d. salina lacks a rigid cell wall, facilitating organelle dna extraction; iv) d. salina is unicellular, with only a single plastid, making it easier, as compared with multicellular species, to develop homoplasmic lines of plastid transformants; and v) being a close relative of the model green algae c. reinhardtii and volvox carteri, means d. salina is an ideal species for comparative plant studies, especially comparative genomics, because the united states department of energy joint genome institute  is sequencing, or has sequenced, the c. reinhardtii, v. carteri, and d. salina nuclear genomes.
in  <dig>  the doe jgi began sequencing the d. salina strain ccap  19/ <dig> nuclear genome, which is approximately  <dig> megabases  in length .
d. salina was selected for genome sequencing because of its potential as a feedstock producer for biofuels production  <cit>  and its model status for studying saline adaptation.
all of the d. salina whole genome shotgun sequencing  trace files that the doe jgi produced are publicly available at the genbank trace archive  <cit> ; soon a complete assembly of the d. salina nuclear dna  will be made public.
the fact that d. salina is the third chlamydomonadalean alga for which there is a genome sequencing project, reaffirms that the chlamydomonadales are emerging as the one of foremost lineages for comparative genomics.
the taxonomic position of the dunaliella genus is still under debate  <cit> ; however, it is often placed within the chlamydomonadales .
for the purpose of this study, our definition of the chlamydomonadales follows that of lewis and mccourt  <cit> , which includes the dunaliella genus, and is equivalent to both the basal-bodies-clockwise group   <cit>  and the volvocales sensu nakada et al.
<cit> .
notably, some strains of d. salina were incorrectly identified in the past, which resulted in the deposition of inaccurately labeled dna sequence data in public databases.
moreover, there exists a debate regarding the delineation of the species d. salina teod.
and dunaliella bardawil avron et ben-amotz  <cit> .
according to borowitzka and borowitzka  <cit>  and borowitzka and siva  <cit> , the species d. bardawil is a nomen nudum; however, the name d. bardawil is still in use.
given the above issues, one should exercise great caution when using dna sequence data in public databases that are said to have originated from d. salina.
as of november  <dig>   <dig>  complete and almost complete organelle dna sequences are available for eight chlamydomonadalean algae, amounting to two plastid and eight mitochondrial genome sequences  <cit> ; moreover, comprehensive genetic maps and limited sequence data are available for the plastid genomes of an additional four taxa  <cit> .
the general features of these organelle genome data, as well as the species from which they are derived, are summarized in table  <dig>  many of the available chlamydomonadalean organelle genome sequences are atypical in one way or another, having extreme sizes , unusual conformations , and/or severely biased nucleotide composition .
furthermore, there can be extensive size, conformational and/or compositional differences among the organelle genomes of closely related chlamydomonadalean species  <cit> .
this wide assortment of genome architectures makes the chlamydomonadales an ideal lineage for studying genome evolution  <cit> .
note: "--", data not available.
clades are defined by nakada et al.
<cit> .
a intronic open reading frames were not considered as coding dna.
b these data vary because of the presence/absence of optional introns.
c these data are based on almost-complete genome sequences.
d mtdna consists of two fragments; data are based on the concatenation of these fragments.
e these data are based on gel-electrophoresis and southern-blot analyses.
f the circular conformation of the v. carteri mtdna is based on genome-assembly data  <cit>  and needs to be confirmed by gel electrophoresis experiments.
g this strain is formally known as polytomella sag 63- <dig>
most of our knowledge of chlamydomonadalean organelle genomes comes from species within the reinhardtinia clade .
indeed, six of the eight chlamydomonadalean algae for which significant organelle dna sequence data are available come from this clade , including c. reinhardtii and v. carteri.
given that the reinhardtinia clade contains only a small fraction of the species diversity within the chlamydomonadales, it would be intriguing to explore the organelle genomes of algae from other chlamydomonadalean clades.
it would be particularly interesting to see if chlamydomonadalean algae from outside the reinhardtinia clade have large, bloated plastid genomes.
both the c. reinhardtii and v. carteri plastid genomes, the only complete  plastid dna sequences that are available from the chlamydomonadales, are among the largest and most noncoding-dna dense plastid genomes observed to date, with sizes of  <dig> and ~ <dig> kilobases , respectively  <cit> .
the forces driving these genomes towards distention are unknown, but they may be connected to the combined effects of a low mutation rate and a low effective population size  <cit> .
here we present the complete mitochondrial dna  and plastid dna  sequences of d. salina strain ccap 19/ <dig> -- a member of the dunaliellinia clade  <cit> .
the salient features of these genomes are described and compared with other chlamydomonadalean organelle genomes, particularly those of c. reinhardtii and v. carteri.
the evolutionary and biotechnological implications of these sequences are discussed.
the overreaching goals of this study are to use the d. salina organelle dna data to test contemporary theories on genome evolution and to lay the foundation for a d. salina plastid transformation system.
RESULTS
overview of the d. salina organelle genomes
the organelle genome sequences of d. salina were assembled using publicly available trace files that the doe jgi d. salina nuclear genome sequencing project produced .
genetic maps of the d. salina organelle genomes are shown in figures  <dig>  and  <dig> ; for comparison, these two figures also include the corresponding genetic maps from c. reinhardtii and volvox carteri.
table  <dig> outlines the general features of the d. salina organelle genomes, including their length, coding and noncoding dna contents and nucleotide compositions, and compares these statistics to those from other chlamydomonadalean organelle dnas.
a venn diagram highlighting the differences in gene content among available mtdna sequences from the chlamydomonadales is presented in figure  <dig>  a schematic compilation comparing the amounts of noncoding dna in the d. salina organelle genomes with those from other completely sequenced  organelle genomes is shown in figure  <dig>  and analyses of the repetitive elements within the d. salina organelle dna are summarized in figure  <dig> and supplementary figures s <dig> and s <dig> .
size, conformation, and nucleotide composition
the d. salina mitochondrial and plastid genomes are  <dig>  and  <dig> kb, respectively, and assemble as circular molecules -- an observation that adds further support to the hypothesis that linear mitochondrial genomes in chlamydomonadalean algae are restricted to species within the reinhardtinia clade .
the mitochondrial genome of d. salina is small relative to those of non-chlamydomonadalean green algae, which are, on average,  <dig>  kb; however, it is still  <dig> - <dig>  kb larger than all other available chlamydomonadalean mtdnas, except for that of v. carteri, which is ~ <dig> kb .
the size of the d. salina plastid genome is more pronounced than its mitochondrial counterpart, being the largest ptdna sequenced thus far.
its nearest rivals are the  <dig>  kb ptdna of the chlorophycean green alga stigeoclonium helveticum and the  <dig>  kb plastid genome of the geranium pelargonium × hortorum  <cit> .
large plastid genomes are a common theme among chlamydomonadalean algae: the c. reinhardtii plastid genome, at  <dig>  kb  <cit> , is the fourth largest completely sequenced ptdna, partial sequence data indicate that the v. carteri ptdna is ~ <dig> kb in length , and gel electrophoresis results place the plastid genomes of chlamydomonas gelatinosa  and chlamydomonas moewusii  at ~ <dig> kb and ~ <dig> kb in length, respectively  <cit> .
the impressive size of the d. salina plastid genome  is a reflection of a prodigious noncoding dna content rather than an unusually large gene repertoire.
within the d. salina ptdna, a pair of inverted repeats, each with a length of  <dig>  kb, divides the genome into a large  and a small single-copy region , referred to as the lsc and ssc regions.
the d. salina inverted repeats are  <dig>  kb shorter than their c. reinhardtii counterparts.
this size discrepancy occurs because the c. reinhardtii inverted repeats contain psba, a gene that is located in the ssc region of the d. salina ptdna.
unlike in d. salina, the ssc and lsc regions of the c. reinhardtii ptdna are virtually indistinguishable with sizes of ~ <dig> kb.
the precise lengths of the lsc, ssc, and inverted repeat regions for the v. carteri plastid genome are unknown; however, preliminary size estimates place them at > <dig> kb  <cit> .
southern blot analyses and partial sequence data indicate that the inverted repeats of the c. moewusii ptdna may be upwards of  <dig> kb  <cit> .
the gc content of the d. salina organelle dnas is  <dig> %  and  <dig> % , which is unremarkable in relation to other archaeplastidial organelle genomes .
however, they are still the most gc-poor  organelle dnas observed within the chlamydomonadales, which is significant because the chlamydomonadales are one of the few lineages out of all eukaryotes known to contain species with gc-rich mitochondrial genomes   <cit> .
the gc content among the different regions of the d. salina mitochondrial and plastid genomes is relatively constant: 33% and 34% for coding dna; 34% and 32% for introns and intronic open reading frames ; and 37% and 31% for intergenic regions.
for the different codon-site positions of the mtdna and ptdna protein-coding regions, the gc content is approximately 38% and 42% ; 38% and 52% ; and 19% and 13% .
cumulative gc-skew analyses,  of the d. salina mtdna show a strong positive correlation with the transcriptional orientation , reflecting the slightly higher gc content of the coding and intronic regions relative to the noncoding mtdna.
the same analysis of the ptdna gives a more disordered plot, but one typical of ptdna, because of the many shifts in transcriptional polarity throughout the genome.
coding content
like other chlamydomonadalean species, d. salina has a severely diminished mtdna gene content of only  <dig> genes, which represent seven proteins, two rrnas, and three trnas.
outside of the chlamydomonadales, the only species known to have more reduced mtdna gene contents are found in the phyla apicomplexa and dinoflagellata and arguably some species within the supergroup excavata  <cit> .
various studies have tried to explain why chlamydomonadalean algae have such reduced mtdna gene contents  <cit> , but at present no straightforward answer to this question exists.
the d. salina mtdna gene inventory mirrors those of chlorogonium elongatum and chlamydomonas eugametos, but it shows some differences to those of reinhardtinia-clade algae.
these differences, which can be visualized on the venn diagram in figure  <dig>  involve changes in trna-coding content and in the number of rrna-coding fragments found on the genome.
for example, the d. salina mtdna encodes three trnas, whereas polytomella mtdna contains only trnm, and the mitochondrial rrns and rrnl genes of d. salina are divided into three  and six  coding modules, whereas in available reinhardtinia-clade mtdnas the rrns and rrnl genes are fragmented into at least four and eight coding modules, respectively.
given the similarities among the d. salina, c. elongatum, and c. eugametos mitochondrial genomes, these findings add further appreciation for the stability of mtdna gene content among chlamydomonadalean species outside of the reinhardtinia clade and underscore the instability of mtdna gene content among reinhardtinia-clade taxa .
the d. salina plastid genome is much more gene rich than its mitochondrial counterpart, with  <dig> genes -- five of which are duplicates found in the inverted repeats.
when ignoring these duplicates, there are a total of  <dig> protein-,  <dig> rrna-, and  <dig> trna-coding genes.
this gene content is reduced from those of green-plant species outside the chlamydomonadales, which on average have  <dig> ptdna-encoded genes, representing  <dig> proteins,  <dig> rrnas, and  <dig> trnas.
it appears that chlamydomonadalean algae, at some point during their evolution, went through a major reduction in ptdna  coding content relative to most other photosynthetic eukaryotes.
the d. salina ptdna gene repertoire is identical to those of c. reinhardtii and v. carteri with the following minor exceptions: i) the d. salina plastid genome encodes three copies of trni -- one more than the c. reinhardtii and v. carteri ptdnas; ii) d. salina, like c. reinhardtii, has two ptdna copies of trne, whereas v. carteri has only one; iii) for d. salina and v. carteri, the rps <dig> gene is represented by a single open reading frame, whereas for c. reinhardtii, rps <dig> is fragmented into two adjacent open reading frames ; and iv) the d. salina ptdna does not contain the roaa-like gene , which is present in the v. carteri ptdna.
preliminary investigations of the plastids from moewusinia-clade algae  <cit>  indicate that their ptdna gene contents are similar to those of d. salina, c. reinhardtii, and v. carteri.
altogether, these findings suggest that ptdna gene content is uniform throughout the chlamydomonadales, save for some minor differences in the number of trna-coding genes.
gene order
all  <dig> genes on the d. salina mitochondrial genome are encoded on the same strand , a characteristic shared by the three other circular-mapping chlamydomonadalean mtdnas sequenced thus far .
sequence data from chlamydomonadalean algae whose mtdnas map as linear molecules, such as c. reinhardtii and polytomella spp., reveal genomes that have two unequally sized gene clusters  with opposing transcriptional polarities, which proceed outwards toward the ends of the genome .
based on these and the above observations, it is reasonable to assume that the ancestral chlamydomonadalean mtdna mapped as a circular molecule with ~ <dig> genes, all of which were encoded on the same strand, and that the events giving rise to linear mtdna were connected with, or resulted in, a shift in transcriptional orientation of approximately one-third of the genes.
the mtdna gene order of d. salina is unique and differs from those of other chlamydomonadalean algae.
very few conserved mtdna gene clusters are shared among d. salina and other chlamydomonadalean species , but this is not surprising considering that mtdna gene arrangements can vary significantly even among closely related species in this group  <cit> .
previous reports suggest that homologous or illegitimate recombination between mtdna repeats is causing mitochondrial genome rearrangements in chlamydomonadalean algae  <cit> .
the d. salina mitochondrial genome does contain minor amounts of repetitive dna ; these repeats may be catalysts for genome rearrangements.
in c. reinhardtii, mitochondrial genes are organized into operons, which are first transcribed into polycistronic primary transcripts and then subsequently processed into mature monocistronic units via endo- and exonucleolytic cleavage  <cit> .
a scan of the d. salina mitochondrial genetic map reveals clusters of tightly packed genes separated by large stretches of noncoding dna .
these gene clusters may reflect the layout of operons in the genome, a theory supported by the fact that they are punctuated by regions of noncoding dna that can be folded into secondary structures.
in contrast to the mtdna, genes in the d. salina plastid genome are found on both strands and occur in small groups of two to four genes , which are distributed among the lsc and ssc regions and the inverted repeats .
the former two regions contain approximately  <dig> and  <dig> genes, respectively, whereas the inverted repeats contain only five genes -- fewer than any chlamydomonadalean inverted repeat explored heretofore.
the following additional genes are observed in the inverted repeats of other chlamydomonadalean algae: psba , rbcl , and atpb .
regions of gene synteny between the ptdna of d. salina and those of c. reinhardtii and v. carteri  are highlighted in gray on figure  <dig>  the allocation of d. salina genes into small clusters is consistent with what is known for the c. reinhardtii plastid genome, where genes appear to be transcribed into monocistronic and dicistronic transcripts  <cit>  rather than the larger polycistronic transcripts that are observed for the mtdna.
thus, regions of gene colinearity between d. salina and c. reinhardtii  may represent conserved transcriptional units.
introns and intergenic regions
one of the more salient features of the d. salina organelle genomes is their noncoding dna content: 58% of the mtdna and  <dig> % of the ptdna consist of either intergenic or intronic dna.
these values approach those of the v. carteri mitochondrial  and plastid  genomes, which are currently the most inflated organelle dna sequences from the chlorophyta  .
in fact, next to v. carteri, the d. salina ptdna has a greater noncoding dna composition than any other plastid genome sequenced to date, exceeding that of the legume trifolium subterraneum  and c. reinhardtii    <cit> .
however, one would expect the unsequenced plastid genomes of c. gelatinosa and c. moewusii, based on their estimated sizes  <cit> , to have more noncoding dna than d. salina but less than that of v. carteri .
interestingly, both the mitochondrial and plastid genomes of d. salina have equally large noncoding dna densities .
this observation goes against what is seen in c. reinhardtii where the mtdna and ptdna have opposing architectures , but it is consistent with the v. carteri organelle genomes, which are both distended with noncoding dna .
the noncoding dna in the d. salina organelle genomes can be subdivided into two categories: intergenic regions, which make up  <dig>  kb  of the mtdna and  <dig>  kb  of the ptdna, and introns and intronic orfs, which together represent  <dig>  kb  and  <dig>  kb  of mitochondrial and plastid genomes, respectively.
for the d. salina ptdna, it is sometimes difficult to distinguish between intergenic dna and intronic dna because intron-like sequences  are found in many of the intergenic regions .
altogether,  <dig> putative group-i introns were found in the mtdna  and  <dig> putative introns were discerned in the ptdna:  <dig> within genes  and  <dig> within intergenic regions .
see figures  <dig> and  <dig> as well as supplementary table  <dig>  for a comparison of the organelle genome intron content of d. salina with those of c. reinhardtii and v. carteri.
note, because of the inverted repeats,  <dig> of the  <dig> introns in the ptdna are duplicates .
seventeen of the gene-located ptdna introns contain orfs , whereas no orfs were found within the ptdna intergenic introns.
the remnants of eight intronic orfs  were found in the intergenic ptdna regions; these pseudo orfs, which are often located adjacent to intergenic introns , appear to be nonfunctional because they contain frameshifts in their coding regions.
all eight of the pseudo orfs show sequence similarity to genes that are typically found in either group-i or group-ii introns , such as genes coding for integrase-, maturase-, reverse-transcriptase- and endonuclease-like proteins.
intergenic introns have been identified in other genomes  <cit> , but until now they had never been observed in chlamydomonadalean organelle dna, or, to the best of our knowledge, any other green-algal organelle genomes.
most of the intergenic introns are highly derived and could only be identified using domain v, which is the most conserved secondary structure element of group-ii introns  <cit> .
further experiments will need to be performed to confirm that the intergenic introns are functional  rather than inert sequences.
if they are functional, then it would imply that many of the intergenic regions of the d. salina plastid genome are transcribed.
there is also the possibility that the individual intergenic introns represent the fragments of larger introns that assemble after translation .
however, secondary structure modeling of the intergenic introns gave no obvious indications that this was the case.
the intron/gene ratios for the d. salina mitochondrial and plastid genomes are  <dig>  and  <dig> , respectively.
these values are much larger than those of other chlamydomonadalean organelle genomes, which range from  <dig> to  <dig>  for available mitochondrial genomes, and are less than  <dig>  for the two available ptdna sequences.
notably, the d. salina ptdna intron/gene ratio exceeds the average value for land-plant mitochondrial genomes   <cit> , which are considered to be among the most intron-dense organelle dnas.
the d. salina organelle dnas contain significantly more introns than their c. reinhardtii and v. carteri counterparts .
in stark contrast to the  <dig> introns found in the d. salina mtdna, c. reinhardtii and v. carteri have five  and four  mitochondrial introns, respectively.
moreover, for c. reinhardtii all five introns are "optional" and, as of yet, the maximum number found in a single strain is three  <cit> .
a similar trend is observed for the ptdna: when counting duplicate genes only once, the plastid genome of c. reinhardtii has six introns  and that of v. carteri has eight ; both these values are significantly less than the  <dig> unique putative introns found in the d. salina plastid genome.
interestingly, for all of the genes that contain introns in the c. reinhardtii and v. carteri organelle dnas, their homologues in d. salina also contain introns, with the exception of cema, which is intronless in d. salina but contains a group-ii intron in v. carteri.
the number of introns per gene and the intron insertion sites can differ among these three algae.
other notable features of the d. salina noncoding organelle dna include three pseudogenes  in the plastid genome.
these pseudogenes, whose functional copies are also present in the ptdna, were classified as such, not because they contain frameshifts in their coding sequence or because they appear highly degenerate relative to their functional counterparts, but because they are missing the first half or first two-thirds of their coding sequences.
furthermore, atpbψ and chllψ are located immediately downstream of their functional copies, and in both cases a group-i intron is sandwiched between the functional gene and the pseudogene .
repeats
unlike the mtdna, which is relatively devoid of repeats, the d. salina plastid genome abounds with repetitive elements.
the difference in repeat content between the mitochondrial and plastid genomes can be visualized by comparing their respective dotplot similarity matrices, which are shown in supplementary figures s <dig>  and s <dig>  .
looking at the ptdna dotplot, it is apparent that the ptdna repeats are found in intergenic regions, introns, and in some of the longer protein-coding genes, such as ftsh, rpoc <dig>  and ycf <dig>  nucleotide blast analyses of the d. salina plastid genome indicate that there are upwards of  <dig>  repeats in the ptdna, forming approximately  <dig> repeat subclasses.
with some exceptions, these repeats range from 30- <dig> nt in length and are 70-90% at.
some shorter  gc-rich repeat elements  were also identified.
the high degree of sequence similarity among the different ptdna repeats is attributable to homopolymer runs rather than a recurring sequence motif.
the c. reinhardtii and v. carteri plastid genomes are also rich in repetitive dna and southern blot analyses suggest that the c. gelatinosa ptdna is as well  <cit> .
presumably, c. pitschmannii plastid genome, the smallest ptdna observed from the chlorophyceae , contains fewer repeats than other chlamydomonadalean ptdnas  <cit> , implying a connection between repeat content and plastid genome size.
in a general sense, the d. salina ptdna repeats are analogous to the short dispersed repeats described for the c. reinhardtii ptdna  <cit> , but they lack the consistent motif of the palindromic elements of the v. carteri plastid genome  <cit> .
the d. salina mtdna dotplot reveals a mostly blank matrix, with the exception of some small diagonal lines, which correspond to palindromic repeats .
clusters of palindromic repeats are found in  <dig> different noncoding regions of the mitochondrial genome.
the genomic breadth of these clusters can be seen in figure  <dig>  where purple asterisks pinpoint the precise location of these repeat clusters.
a consensus sequence of one of the more frequently occurring clusters is depicted in figure  <dig>  the mean length of the palindromic repeat clusters is  <dig> nt and, on average, they have a gc content of 50%.
each cluster contains approximately three palindromic elements , and the individual palindromes within each cluster range from 23- <dig> nt in length, and can be either at- or gc-rich .
ten of the eighteen intergenic regions that are found in the mtdna contain either one or two palindromic clusters , an arrangement that suggests that the palindromes may play a role in gene processing.
palindromic repeats have been described in other chlamydomonadalean mitochondrial genomes, including those of c. reinhardtii, v. carteri, and polytomella capuana  <cit> , and in many cases chlamydomonadalean mtdna palindromic repeats have been implicated in rna processing  <cit> .
a notable observation is that the mitochondrial palindromic repeats of both d. salina and v. carteri often contain the sequence 5'-ttta-3'  in the loops of their hairpin structures .
organelle dna from other d. salina strains
prior to this study, organelle dna sequence data for d. salina was limited to  <dig> genbank entries amounting to  <dig>  kb of ptdna, divided over six protein-coding loci .
most of these  <dig> entries do not list the strains of d. salina that were used for sequencing and only a few are associated with published articles.
alignments of the ccap 19/ <dig> ptdna data generated here with the d. salina ptdna data available at genbank reveal a significant amount of sequence divergence: pairwise nucleotide diversity values varied from  <dig> % to  <dig> %, which suggests that none of the  <dig> d. salina entries are derived from strain ccap 19/ <dig> and that at least some of the sequences at genbank labeled d. salina are not the same species as ccap 19/ <dig>  in comparison, the average silent-site ptdna diversity among seven different north american isolates of c. reinhardtii was estimated to be  <dig> %  <cit> , and a recent study on seven v. carteri geographical isolates found πsilent of the ptdna to be ~ <dig> %  <cit> .
the high degree of sequence diversity among d. salina genbank entries appears to be a reflection of the known issue of misidentification of dunaliella salina isolates.
dunaliella researchers should be mindful of the high levels of sequence divergence between ccap 19/ <dig> and the available d. salina genbank entries if using these data for sequence-based studies, such as designing pcr primers.
paving the way towards plastid transformation
the development of a reliable, high-efficiency genetic transformation system for d. salina is an important objective for the dunaliella research community, especially when considering the many industrial applications that this technology could provide.
genetic transformation of the d. salina nuclear genome has been successful  <cit> ; however, a lack of ptdna sequence data has prevented successful attempts at transforming the d. salina plastid genome -- although an unsuccessful attempt was made to transform the plastid of dunaliella tertiolecta  <cit> .
for dunaliella species, there are significant advantages to plastid transformation over nuclear transformation , some of which relate to the fact that ptdna is polyploid and experiences high levels of gene expression.
now that the d. salina ptdna sequence is available, scientists will be able to use these data to develop plastid transformation vectors targeting specific regions of the d. salina ptdna.
plastid transformation occurs via homologous recombination between an engineered vector and a selected region of the plastid genome.
moreover, promoter and 3' utr regions of genes  <cit>  can now be used to design vectors for d. salina with improved expression levels.
in principle, transgenes can be integrated into any site of the plastid genome, but transcriptionally-active intergenic regions are ideal  <cit> .
one of the most frequently used integration sites is the trni-trna intergenic spacer, which is found in the inverted repeat of most plastid genomes, including that of d. salina.
given its popularity, the trni-trna intergenic region would be an ideal site for early attempts at transforming the d. salina ptdna.
however, the discovery of intergenic introns in the d. salina ptdna is an indication that many of the intergenic regions are transcriptionally active, which, if true, should allow for a diverse range of transformation targeting sites.
CONCLUSIONS
the d. salina organelle genomes are large, intron-dense molecules comprised predominantly of noncoding nucleotides.
repetitive elements punctuate the noncoding regions of the mitochondrial and plastid genomes, but are much more prevalent in the ptdna.
overall, the discovery of putative intergenic introns in the d. salina ptdna adds a new layer of complexity to the diverse repertoire of organelle genome architectures found in the chlamydomonadales.
already a model organism for synthesizing β-carotene and studying salt adaptation, d. salina is now an ideal species for investigating organelle genome expansion.
the high level of repetitive elements found in the plastid genome of d. salina may mirror the high level of repetitive sequences that is expected to be present in the nucdna .
publication of the plastid genome of d. salina is expected to result in major advances of plastid engineering with generation and use of transgenic d. salina strains for a number of new applications in the fields of biofuels as well as in vaccine antigene and biopharmaceutical production.
the complete sequence of the d. salina nuclear genome will be made available soon, placing d. salina in a group of a select few photosynthetic eukaryotes for which complete genome sequences from all three genetic compartments are available.

 dunaliella salina teodoresco is one of the best-studied unicellular green algae. this is not only because d. salina is halotolerant, thriving in extreme saline environments, but also because it can produce large quantities of β-carotene in lipid globules located within the chloroplast. these traits make d. salina a model organism for investigating the evolution of salt adaptation and an attractive "cell factory" for the commercial production of β-carotene. although a great deal is known about the physiology and biochemistry of d. salina, very little is known about the genomes it carries, especially those within its organelles. until now, nothing was known about the size, conformation, or gene complement of either the mitochondrial or plastid genomes of d. salina even though the sequences of these genomes are essential to the development of new d. salina technologies, such as a viable plastid transformation system. research on green-algal organelle genomes has led to significant advancements in genetic engineering. the first stable transformation of a plastid genome was achieved in using the unicellular green alga chlamydomonas reinhardtii and, soon after, the first example of recombinant protein expression in a plastid was also achieved using c. reinhardtii. since then, many techniques for plastid engineering have been first developed for green algae and then adapted for use in land plants. given the relatively close evolutionary proximity of c. reinhardtii and d. salina, it is reasonable to assume that many of the technologies for c. reinhardtii plastid transformation might be transferable to d. salina. various groups have attempted to transform d. salina; however, a lack of plastid-genome sequence data has prevented successful plastid transformation. therefore, the first step in developing an efficient and reliable plastid transformation system for d. salina is to sequence its organelle genomes. d. salina is an attractive alga for organelle genome research and plastome engineering for a variety of reasons: i) various strains and geographical isolates of d. salina are readily available from algal culture collections around the world; ii) d. salina is relatively easy to grow and maintain -- it is one of the few microalgae that are being cultivated currently on a large scale; iii) d. salina lacks a rigid cell wall, facilitating organelle dna extraction; iv) d. salina is unicellular, with only a single plastid, making it easier, as compared with multicellular species, to develop homoplasmic lines of plastid transformants; and v) being a close relative of the model green algae c. reinhardtii and volvox carteri, means d. salina is an ideal species for comparative plant studies, especially comparative genomics, because the united states department of energy joint genome institute is sequencing, or has sequenced, the c. reinhardtii, v. carteri, and d. salina nuclear genomes. in the doe jgi began sequencing the d. salina strain ccap 19/ nuclear genome, which is approximately megabases in length. d. salina was selected for genome sequencing because of its potential as a feedstock producer for biofuels production and its model status for studying saline adaptation. all of the d. salina whole genome shotgun sequencing trace files that the doe jgi produced are publicly available at the genbank trace archive; soon a complete assembly of the d. salina nuclear dna will be made public. the fact that d. salina is the third chlamydomonadalean alga for which there is a genome sequencing project, reaffirms that the chlamydomonadales are emerging as the one of foremost lineages for comparative genomics. the taxonomic position of the dunaliella genus is still under debate; however, it is often placed within the chlamydomonadales. for the purpose of this study, our definition of the chlamydomonadales follows that of lewis and mccourt, which includes the dunaliella genus, and is equivalent to both the basal-bodies-clockwise group and the volvocales sensu nakada et al.. notably, some strains of d. salina were incorrectly identified in the past, which resulted in the deposition of inaccurately labeled dna sequence data in public databases. moreover, there exists a debate regarding the delineation of the species d. salina teod. and dunaliella bardawil avron et ben-amotz. according to borowitzka and borowitzka and borowitzka and siva, the species d. bardawil is a nomen nudum; however, the name d. bardawil is still in use. given the above issues, one should exercise great caution when using dna sequence data in public databases that are said to have originated from d. salina. as of november complete and almost complete organelle dna sequences are available for eight chlamydomonadalean algae, amounting to two plastid and eight mitochondrial genome sequences; moreover, comprehensive genetic maps and limited sequence data are available for the plastid genomes of an additional four taxa. the general features of these organelle genome data, as well as the species from which they are derived, are summarized in table many of the available chlamydomonadalean organelle genome sequences are atypical in one way or another, having extreme sizes, unusual conformations, and/or severely biased nucleotide composition. furthermore, there can be extensive size, conformational and/or compositional differences among the organelle genomes of closely related chlamydomonadalean species. this wide assortment of genome architectures makes the chlamydomonadales an ideal lineage for studying genome evolution. note: "--", data not available. clades are defined by nakada et al.. a intronic open reading frames were not considered as coding dna. b these data vary because of the presence/absence of optional introns. c these data are based on almost-complete genome sequences. d mtdna consists of two fragments; data are based on the concatenation of these fragments. e these data are based on gel-electrophoresis and southern-blot analyses. f the circular conformation of the v. carteri mtdna is based on genome-assembly data and needs to be confirmed by gel electrophoresis experiments. g this strain is formally known as polytomella sag 63- most of our knowledge of chlamydomonadalean organelle genomes comes from species within the reinhardtinia clade. indeed, six of the eight chlamydomonadalean algae for which significant organelle dna sequence data are available come from this clade, including c. reinhardtii and v. carteri. given that the reinhardtinia clade contains only a small fraction of the species diversity within the chlamydomonadales, it would be intriguing to explore the organelle genomes of algae from other chlamydomonadalean clades. it would be particularly interesting to see if chlamydomonadalean algae from outside the reinhardtinia clade have large, bloated plastid genomes. both the c. reinhardtii and v. carteri plastid genomes, the only complete plastid dna sequences that are available from the chlamydomonadales, are among the largest and most noncoding-dna dense plastid genomes observed to date, with sizes of and ~ kilobases, respectively. the forces driving these genomes towards distention are unknown, but they may be connected to the combined effects of a low mutation rate and a low effective population size. here we present the complete mitochondrial dna and plastid dna sequences of d. salina strain ccap 19/ -- a member of the dunaliellinia clade. the salient features of these genomes are described and compared with other chlamydomonadalean organelle genomes, particularly those of c. reinhardtii and v. carteri. the evolutionary and biotechnological implications of these sequences are discussed. the overreaching goals of this study are to use the d. salina organelle dna data to test contemporary theories on genome evolution and to lay the foundation for a d. salina plastid transformation system.  overview of the d. salina organelle genomes the organelle genome sequences of d. salina were assembled using publicly available trace files that the doe jgi d. salina nuclear genome sequencing project produced. genetic maps of the d. salina organelle genomes are shown in figures and; for comparison, these two figures also include the corresponding genetic maps from c. reinhardtii and volvox carteri. table outlines the general features of the d. salina organelle genomes, including their length, coding and noncoding dna contents and nucleotide compositions, and compares these statistics to those from other chlamydomonadalean organelle dnas. a venn diagram highlighting the differences in gene content among available mtdna sequences from the chlamydomonadales is presented in figure a schematic compilation comparing the amounts of noncoding dna in the d. salina organelle genomes with those from other completely sequenced organelle genomes is shown in figure and analyses of the repetitive elements within the d. salina organelle dna are summarized in figure and supplementary figures s and s. size, conformation, and nucleotide composition the d. salina mitochondrial and plastid genomes are and kb, respectively, and assemble as circular molecules -- an observation that adds further support to the hypothesis that linear mitochondrial genomes in chlamydomonadalean algae are restricted to species within the reinhardtinia clade. the mitochondrial genome of d. salina is small relative to those of non-chlamydomonadalean green algae, which are, on average, kb; however, it is still - kb larger than all other available chlamydomonadalean mtdnas, except for that of v. carteri, which is ~ kb. the size of the d. salina plastid genome is more pronounced than its mitochondrial counterpart, being the largest ptdna sequenced thus far. its nearest rivals are the kb ptdna of the chlorophycean green alga stigeoclonium helveticum and the kb plastid genome of the geranium pelargonium × hortorum. large plastid genomes are a common theme among chlamydomonadalean algae: the c. reinhardtii plastid genome, at kb, is the fourth largest completely sequenced ptdna, partial sequence data indicate that the v. carteri ptdna is ~ kb in length, and gel electrophoresis results place the plastid genomes of chlamydomonas gelatinosa and chlamydomonas moewusii at ~ kb and ~ kb in length, respectively. the impressive size of the d. salina plastid genome is a reflection of a prodigious noncoding dna content rather than an unusually large gene repertoire. within the d. salina ptdna, a pair of inverted repeats, each with a length of kb, divides the genome into a large and a small single-copy region, referred to as the lsc and ssc regions. the d. salina inverted repeats are kb shorter than their c. reinhardtii counterparts. this size discrepancy occurs because the c. reinhardtii inverted repeats contain psba, a gene that is located in the ssc region of the d. salina ptdna. unlike in d. salina, the ssc and lsc regions of the c. reinhardtii ptdna are virtually indistinguishable with sizes of ~ kb. the precise lengths of the lsc, ssc, and inverted repeat regions for the v. carteri plastid genome are unknown; however, preliminary size estimates place them at > kb. southern blot analyses and partial sequence data indicate that the inverted repeats of the c. moewusii ptdna may be upwards of kb. the gc content of the d. salina organelle dnas is % and %, which is unremarkable in relation to other archaeplastidial organelle genomes. however, they are still the most gc-poor organelle dnas observed within the chlamydomonadales, which is significant because the chlamydomonadales are one of the few lineages out of all eukaryotes known to contain species with gc-rich mitochondrial genomes. the gc content among the different regions of the d. salina mitochondrial and plastid genomes is relatively constant: 33% and 34% for coding dna; 34% and 32% for introns and intronic open reading frames; and 37% and 31% for intergenic regions. for the different codon-site positions of the mtdna and ptdna protein-coding regions, the gc content is approximately 38% and 42%; 38% and 52%; and 19% and 13%. cumulative gc-skew analyses, of the d. salina mtdna show a strong positive correlation with the transcriptional orientation, reflecting the slightly higher gc content of the coding and intronic regions relative to the noncoding mtdna. the same analysis of the ptdna gives a more disordered plot, but one typical of ptdna, because of the many shifts in transcriptional polarity throughout the genome. coding content like other chlamydomonadalean species, d. salina has a severely diminished mtdna gene content of only genes, which represent seven proteins, two rrnas, and three trnas. outside of the chlamydomonadales, the only species known to have more reduced mtdna gene contents are found in the phyla apicomplexa and dinoflagellata and arguably some species within the supergroup excavata. various studies have tried to explain why chlamydomonadalean algae have such reduced mtdna gene contents, but at present no straightforward answer to this question exists. the d. salina mtdna gene inventory mirrors those of chlorogonium elongatum and chlamydomonas eugametos, but it shows some differences to those of reinhardtinia-clade algae. these differences, which can be visualized on the venn diagram in figure involve changes in trna-coding content and in the number of rrna-coding fragments found on the genome. for example, the d. salina mtdna encodes three trnas, whereas polytomella mtdna contains only trnm, and the mitochondrial rrns and rrnl genes of d. salina are divided into three and six coding modules, whereas in available reinhardtinia-clade mtdnas the rrns and rrnl genes are fragmented into at least four and eight coding modules, respectively. given the similarities among the d. salina, c. elongatum, and c. eugametos mitochondrial genomes, these findings add further appreciation for the stability of mtdna gene content among chlamydomonadalean species outside of the reinhardtinia clade and underscore the instability of mtdna gene content among reinhardtinia-clade taxa. the d. salina plastid genome is much more gene rich than its mitochondrial counterpart, with genes -- five of which are duplicates found in the inverted repeats. when ignoring these duplicates, there are a total of protein-, rrna-, and trna-coding genes. this gene content is reduced from those of green-plant species outside the chlamydomonadales, which on average have ptdna-encoded genes, representing proteins, rrnas, and trnas. it appears that chlamydomonadalean algae, at some point during their evolution, went through a major reduction in ptdna coding content relative to most other photosynthetic eukaryotes. the d. salina ptdna gene repertoire is identical to those of c. reinhardtii and v. carteri with the following minor exceptions: i) the d. salina plastid genome encodes three copies of trni -- one more than the c. reinhardtii and v. carteri ptdnas; ii) d. salina, like c. reinhardtii, has two ptdna copies of trne, whereas v. carteri has only one; iii) for d. salina and v. carteri, the rps gene is represented by a single open reading frame, whereas for c. reinhardtii, rps is fragmented into two adjacent open reading frames; and iv) the d. salina ptdna does not contain the roaa-like gene, which is present in the v. carteri ptdna. preliminary investigations of the plastids from moewusinia-clade algae indicate that their ptdna gene contents are similar to those of d. salina, c. reinhardtii, and v. carteri. altogether, these findings suggest that ptdna gene content is uniform throughout the chlamydomonadales, save for some minor differences in the number of trna-coding genes. gene order all genes on the d. salina mitochondrial genome are encoded on the same strand, a characteristic shared by the three other circular-mapping chlamydomonadalean mtdnas sequenced thus far. sequence data from chlamydomonadalean algae whose mtdnas map as linear molecules, such as c. reinhardtii and polytomella spp., reveal genomes that have two unequally sized gene clusters with opposing transcriptional polarities, which proceed outwards toward the ends of the genome. based on these and the above observations, it is reasonable to assume that the ancestral chlamydomonadalean mtdna mapped as a circular molecule with ~ genes, all of which were encoded on the same strand, and that the events giving rise to linear mtdna were connected with, or resulted in, a shift in transcriptional orientation of approximately one-third of the genes. the mtdna gene order of d. salina is unique and differs from those of other chlamydomonadalean algae. very few conserved mtdna gene clusters are shared among d. salina and other chlamydomonadalean species, but this is not surprising considering that mtdna gene arrangements can vary significantly even among closely related species in this group. previous reports suggest that homologous or illegitimate recombination between mtdna repeats is causing mitochondrial genome rearrangements in chlamydomonadalean algae. the d. salina mitochondrial genome does contain minor amounts of repetitive dna; these repeats may be catalysts for genome rearrangements. in c. reinhardtii, mitochondrial genes are organized into operons, which are first transcribed into polycistronic primary transcripts and then subsequently processed into mature monocistronic units via endo- and exonucleolytic cleavage. a scan of the d. salina mitochondrial genetic map reveals clusters of tightly packed genes separated by large stretches of noncoding dna. these gene clusters may reflect the layout of operons in the genome, a theory supported by the fact that they are punctuated by regions of noncoding dna that can be folded into secondary structures. in contrast to the mtdna, genes in the d. salina plastid genome are found on both strands and occur in small groups of two to four genes, which are distributed among the lsc and ssc regions and the inverted repeats. the former two regions contain approximately and genes, respectively, whereas the inverted repeats contain only five genes -- fewer than any chlamydomonadalean inverted repeat explored heretofore. the following additional genes are observed in the inverted repeats of other chlamydomonadalean algae: psba, rbcl, and atpb. regions of gene synteny between the ptdna of d. salina and those of c. reinhardtii and v. carteri are highlighted in gray on figure the allocation of d. salina genes into small clusters is consistent with what is known for the c. reinhardtii plastid genome, where genes appear to be transcribed into monocistronic and dicistronic transcripts rather than the larger polycistronic transcripts that are observed for the mtdna. thus, regions of gene colinearity between d. salina and c. reinhardtii may represent conserved transcriptional units. introns and intergenic regions one of the more salient features of the d. salina organelle genomes is their noncoding dna content: 58% of the mtdna and % of the ptdna consist of either intergenic or intronic dna. these values approach those of the v. carteri mitochondrial and plastid genomes, which are currently the most inflated organelle dna sequences from the chlorophyta. in fact, next to v. carteri, the d. salina ptdna has a greater noncoding dna composition than any other plastid genome sequenced to date, exceeding that of the legume trifolium subterraneum and c. reinhardtii. however, one would expect the unsequenced plastid genomes of c. gelatinosa and c. moewusii, based on their estimated sizes, to have more noncoding dna than d. salina but less than that of v. carteri. interestingly, both the mitochondrial and plastid genomes of d. salina have equally large noncoding dna densities. this observation goes against what is seen in c. reinhardtii where the mtdna and ptdna have opposing architectures, but it is consistent with the v. carteri organelle genomes, which are both distended with noncoding dna. the noncoding dna in the d. salina organelle genomes can be subdivided into two categories: intergenic regions, which make up kb of the mtdna and kb of the ptdna, and introns and intronic orfs, which together represent kb and kb of mitochondrial and plastid genomes, respectively. for the d. salina ptdna, it is sometimes difficult to distinguish between intergenic dna and intronic dna because intron-like sequences are found in many of the intergenic regions. altogether, putative group-i introns were found in the mtdna and putative introns were discerned in the ptdna: within genes and within intergenic regions. see figures and as well as supplementary table for a comparison of the organelle genome intron content of d. salina with those of c. reinhardtii and v. carteri. note, because of the inverted repeats, of the introns in the ptdna are duplicates. seventeen of the gene-located ptdna introns contain orfs, whereas no orfs were found within the ptdna intergenic introns. the remnants of eight intronic orfs were found in the intergenic ptdna regions; these pseudo orfs, which are often located adjacent to intergenic introns, appear to be nonfunctional because they contain frameshifts in their coding regions. all eight of the pseudo orfs show sequence similarity to genes that are typically found in either group-i or group-ii introns, such as genes coding for integrase-, maturase-, reverse-transcriptase- and endonuclease-like proteins. intergenic introns have been identified in other genomes, but until now they had never been observed in chlamydomonadalean organelle dna, or, to the best of our knowledge, any other green-algal organelle genomes. most of the intergenic introns are highly derived and could only be identified using domain v, which is the most conserved secondary structure element of group-ii introns. further experiments will need to be performed to confirm that the intergenic introns are functional rather than inert sequences. if they are functional, then it would imply that many of the intergenic regions of the d. salina plastid genome are transcribed. there is also the possibility that the individual intergenic introns represent the fragments of larger introns that assemble after translation. however, secondary structure modeling of the intergenic introns gave no obvious indications that this was the case. the intron/gene ratios for the d. salina mitochondrial and plastid genomes are and, respectively. these values are much larger than those of other chlamydomonadalean organelle genomes, which range from to for available mitochondrial genomes, and are less than for the two available ptdna sequences. notably, the d. salina ptdna intron/gene ratio exceeds the average value for land-plant mitochondrial genomes, which are considered to be among the most intron-dense organelle dnas. the d. salina organelle dnas contain significantly more introns than their c. reinhardtii and v. carteri counterparts. in stark contrast to the introns found in the d. salina mtdna, c. reinhardtii and v. carteri have five and four mitochondrial introns, respectively. moreover, for c. reinhardtii all five introns are "optional" and, as of yet, the maximum number found in a single strain is three. a similar trend is observed for the ptdna: when counting duplicate genes only once, the plastid genome of c. reinhardtii has six introns and that of v. carteri has eight; both these values are significantly less than the unique putative introns found in the d. salina plastid genome. interestingly, for all of the genes that contain introns in the c. reinhardtii and v. carteri organelle dnas, their homologues in d. salina also contain introns, with the exception of cema, which is intronless in d. salina but contains a group-ii intron in v. carteri. the number of introns per gene and the intron insertion sites can differ among these three algae. other notable features of the d. salina noncoding organelle dna include three pseudogenes in the plastid genome. these pseudogenes, whose functional copies are also present in the ptdna, were classified as such, not because they contain frameshifts in their coding sequence or because they appear highly degenerate relative to their functional counterparts, but because they are missing the first half or first two-thirds of their coding sequences. furthermore, atpbψ and chllψ are located immediately downstream of their functional copies, and in both cases a group-i intron is sandwiched between the functional gene and the pseudogene. repeats unlike the mtdna, which is relatively devoid of repeats, the d. salina plastid genome abounds with repetitive elements. the difference in repeat content between the mitochondrial and plastid genomes can be visualized by comparing their respective dotplot similarity matrices, which are shown in supplementary figures s and s. looking at the ptdna dotplot, it is apparent that the ptdna repeats are found in intergenic regions, introns, and in some of the longer protein-coding genes, such as ftsh, rpoc and ycf nucleotide blast analyses of the d. salina plastid genome indicate that there are upwards of repeats in the ptdna, forming approximately repeat subclasses. with some exceptions, these repeats range from 30- nt in length and are 70-90% at. some shorter gc-rich repeat elements were also identified. the high degree of sequence similarity among the different ptdna repeats is attributable to homopolymer runs rather than a recurring sequence motif. the c. reinhardtii and v. carteri plastid genomes are also rich in repetitive dna and southern blot analyses suggest that the c. gelatinosa ptdna is as well. presumably, c. pitschmannii plastid genome, the smallest ptdna observed from the chlorophyceae, contains fewer repeats than other chlamydomonadalean ptdnas, implying a connection between repeat content and plastid genome size. in a general sense, the d. salina ptdna repeats are analogous to the short dispersed repeats described for the c. reinhardtii ptdna, but they lack the consistent motif of the palindromic elements of the v. carteri plastid genome. the d. salina mtdna dotplot reveals a mostly blank matrix, with the exception of some small diagonal lines, which correspond to palindromic repeats. clusters of palindromic repeats are found in different noncoding regions of the mitochondrial genome. the genomic breadth of these clusters can be seen in figure where purple asterisks pinpoint the precise location of these repeat clusters. a consensus sequence of one of the more frequently occurring clusters is depicted in figure the mean length of the palindromic repeat clusters is nt and, on average, they have a gc content of 50%. each cluster contains approximately three palindromic elements, and the individual palindromes within each cluster range from 23- nt in length, and can be either at- or gc-rich. ten of the eighteen intergenic regions that are found in the mtdna contain either one or two palindromic clusters, an arrangement that suggests that the palindromes may play a role in gene processing. palindromic repeats have been described in other chlamydomonadalean mitochondrial genomes, including those of c. reinhardtii, v. carteri, and polytomella capuana, and in many cases chlamydomonadalean mtdna palindromic repeats have been implicated in rna processing. a notable observation is that the mitochondrial palindromic repeats of both d. salina and v. carteri often contain the sequence 5'-ttta-3' in the loops of their hairpin structures. organelle dna from other d. salina strains prior to this study, organelle dna sequence data for d. salina was limited to genbank entries amounting to kb of ptdna, divided over six protein-coding loci. most of these entries do not list the strains of d. salina that were used for sequencing and only a few are associated with published articles. alignments of the ccap 19/ ptdna data generated here with the d. salina ptdna data available at genbank reveal a significant amount of sequence divergence: pairwise nucleotide diversity values varied from % to %, which suggests that none of the d. salina entries are derived from strain ccap 19/ and that at least some of the sequences at genbank labeled d. salina are not the same species as ccap 19/ in comparison, the average silent-site ptdna diversity among seven different north american isolates of c. reinhardtii was estimated to be %, and a recent study on seven v. carteri geographical isolates found πsilent of the ptdna to be ~ %. the high degree of sequence diversity among d. salina genbank entries appears to be a reflection of the known issue of misidentification of dunaliella salina isolates. dunaliella researchers should be mindful of the high levels of sequence divergence between ccap 19/ and the available d. salina genbank entries if using these data for sequence-based studies, such as designing pcr primers. paving the way towards plastid transformation the development of a reliable, high-efficiency genetic transformation system for d. salina is an important objective for the dunaliella research community, especially when considering the many industrial applications that this technology could provide. genetic transformation of the d. salina nuclear genome has been successful; however, a lack of ptdna sequence data has prevented successful attempts at transforming the d. salina plastid genome -- although an unsuccessful attempt was made to transform the plastid of dunaliella tertiolecta. for dunaliella species, there are significant advantages to plastid transformation over nuclear transformation, some of which relate to the fact that ptdna is polyploid and experiences high levels of gene expression. now that the d. salina ptdna sequence is available, scientists will be able to use these data to develop plastid transformation vectors targeting specific regions of the d. salina ptdna. plastid transformation occurs via homologous recombination between an engineered vector and a selected region of the plastid genome. moreover, promoter and 3' utr regions of genes can now be used to design vectors for d. salina with improved expression levels. in principle, transgenes can be integrated into any site of the plastid genome, but transcriptionally-active intergenic regions are ideal. one of the most frequently used integration sites is the trni-trna intergenic spacer, which is found in the inverted repeat of most plastid genomes, including that of d. salina. given its popularity, the trni-trna intergenic region would be an ideal site for early attempts at transforming the d. salina ptdna. however, the discovery of intergenic introns in the d. salina ptdna is an indication that many of the intergenic regions are transcriptionally active, which, if true, should allow for a diverse range of transformation targeting sites.  the d. salina organelle genomes are large, intron-dense molecules comprised predominantly of noncoding nucleotides. repetitive elements punctuate the noncoding regions of the mitochondrial and plastid genomes, but are much more prevalent in the ptdna. overall, the discovery of putative intergenic introns in the d. salina ptdna adds a new layer of complexity to the diverse repertoire of organelle genome architectures found in the chlamydomonadales. already a model organism for synthesizing β-carotene and studying salt adaptation, d. salina is now an ideal species for investigating organelle genome expansion. the high level of repetitive elements found in the plastid genome of d. salina may mirror the high level of repetitive sequences that is expected to be present in the nucdna. publication of the plastid genome of d. salina is expected to result in major advances of plastid engineering with generation and use of transgenic d. salina strains for a number of new applications in the fields of biofuels as well as in vaccine antigene and biopharmaceutical production. the complete sequence of the d. salina nuclear genome will be made available soon, placing d. salina in a group of a select few photosynthetic eukaryotes for which complete genome sequences from all three genetic compartments are available. 
in fact, the d. salina plastid genome, at kb, is the largest complete plastid dna sequence currently deposited in genbank, and both the mitochondrial and plastid genomes have unprecedentedly high intron densities for organelle dna: ~ and ~ introns per gene, respectively. moreover, what appear to be the relics of genes, introns, and intronic open reading frames are found scattered throughout the intergenic ptdna regions -- a trait without parallel in other characterized organelle genomes and one that gives insight into the mechanisms and modes of expansion of the d. salina ptdna. this is because d. salina can produce massive amounts of β-carotene, which can be collected for commercial purposes, and because of its potential as a feedstock for biofuels production. although the biochemistry and physiology of d. salina have been studied in great detail, virtually nothing is known about the genomes it carries, especially those within its mitochondrion and plastid. this study presents the complete mitochondrial and plastid genome sequences of d. salina and compares them with those of the model green algae chlamydomonas reinhardtii and volvox carteri. the d. salina organelle genomes are large, circular-mapping molecules with ~60% noncoding dna, placing them among the most inflated organelle dnas sampled from the chlorophyta. they also suggest that the events giving rise to the expanded ptdna architecture of d. salina and other chlamydomonadales may have occurred early in the evolution of this lineage. although interesting from a genome evolution standpoint, the d. salina organelle dna sequences will aid in the development of a viable plastid transformation system for this model alga, and they will complement the forthcoming d. salina nuclear genome sequence, placing d. salina in a group of a select few photosynthetic eukaryotes for which complete genome sequences from all three genetic compartments are available. these findings confirm the notion that chlamydomonadalean algae have some of the most extreme organelle genomes of all eukaryotes.
 dunaliella salina teodoresco, a unicellular, halophilic green alga belonging to the chlorophyceae, is among the most industrially important microalgae. this is because d. salina can produce massive amounts of β-carotene, which can be collected for commercial purposes, and because of its potential as a feedstock for biofuels production. although the biochemistry and physiology of d. salina have been studied in great detail, virtually nothing is known about the genomes it carries, especially those within its mitochondrion and plastid. this study presents the complete mitochondrial and plastid genome sequences of d. salina and compares them with those of the model green algae chlamydomonas reinhardtii and volvox carteri.  the d. salina organelle genomes are large, circular-mapping molecules with ~60% noncoding dna, placing them among the most inflated organelle dnas sampled from the chlorophyta. in fact, the d. salina plastid genome, at kb, is the largest complete plastid dna sequence currently deposited in genbank, and both the mitochondrial and plastid genomes have unprecedentedly high intron densities for organelle dna: ~ and ~ introns per gene, respectively. moreover, what appear to be the relics of genes, introns, and intronic open reading frames are found scattered throughout the intergenic ptdna regions -- a trait without parallel in other characterized organelle genomes and one that gives insight into the mechanisms and modes of expansion of the d. salina ptdna.  these findings confirm the notion that chlamydomonadalean algae have some of the most extreme organelle genomes of all eukaryotes. they also suggest that the events giving rise to the expanded ptdna architecture of d. salina and other chlamydomonadales may have occurred early in the evolution of this lineage. although interesting from a genome evolution standpoint, the d. salina organelle dna sequences will aid in the development of a viable plastid transformation system for this model alga, and they will complement the forthcoming d. salina nuclear genome sequence, placing d. salina in a group of a select few photosynthetic eukaryotes for which complete genome sequences from all three genetic compartments are available. 