BACKGROUND
the origin of mitochondria by endosymbiosis has emerged as a pivotal event in the evolution of eukaryotes.
all eukaryote groups that have been studied bear a derivative of this endosymbiont, and for most the resulting mitochondrion is central to energy metabolism as well as providing several other anabolic and catabolic functions  <cit> .
a relict, though functionally essential, mitochondrial genome  persists in all but a few anaerobic eukaryotes, and the genes in these genomes firmly identify the original endosymbiont as an α-proteobacterium  <cit> .
the jakobid flagellate reclinomonas americana has the least derived mitochondrial genome characterized to date, with at least  <dig> genes encoded on a single, circular-mapping  <dig> kb chromosome  <cit> .
more typically mitochondrial genomes have been reduced to 40– <dig> genes arranged on either circular- or linear-mapping chromosomes of 15– <dig> kb   <cit> .
in some eukaryotic groups, however, the mtdna has been modified more substantially, resulting in extremes in genome structure.
for example, trypanosomatid mtdna consists of a few dozen large circular molecules and several thousand minicircles that encode guide rnas that participate in extensive u insertion/deletion rna editing  <cit> .
diplonemid mitochondria also contain multiple circular mtdna molecules, each encoding gene fragments that are trans-spliced to generate functional transcripts  <cit> .
another example is the mtdna in the ichthyosporean, amoebidium parasiticum: in this case, mitochondrial genes are fragmented and dispersed over several hundred linear chromosomes, totaling >  <dig> kb  <cit> .
over the diversity of eukaryotes, mitochondrial genomes exhibit other interesting characteristics, including the use of a number of different non-standard genetic codes, many of which involve alterations in start and, more rarely, stop codons  <cit> .
one large group in which particularly interesting mitochondrial genome variation has been found is alveolates.
three major phyla make up alveolates: ciliates, apicomplexans, and dinoflagellates, with apicomplexans and dinoflagellates being sister clades to the exclusion of ciliates  <cit> .
within alveolates, ciliate mtdna is the most conventional, consisting of a linear molecule, 40– <dig> kb in length, that codes for many of the standard mitochondrial proteins found in other organisms  <cit> .
by contrast, the mtdna of the apicomplexan genus plasmodium is the smallest known, consisting of a linear,  <dig> kb tandem repeat  <cit>  with only three protein-coding genes: cytochrome oxidase subunit  <dig> , cytochrome oxidase subunit  <dig>  and cytochrome b .
in addition, ciliate mtdna encodes two ribosomal rnas , but the corresponding apicomplexan genes are fragmented to an unprecedented degree and scattered about the genome  <cit> .
to date, dinoflagellate mtdnas have been the least well studied of alveolate mitochondrial genomes, with existing data pointing to a genome exhibiting several eccentricities.
the first sequences isolated were four copies of cox <dig> from crypthecodinium cohnii, each of which was found to occur in a unique genomic context  <cit> .
southern blots demonstrated multiple different copies of this gene that varied in abundance, suggesting the c. cohnii mitochondrial genome is not as streamlined as in apicomplexans.
subsequently, cob and cox <dig> have been found as well, and multiple, sometimes fragmented copies of these genes have now been reported from diverse dinoflagellates   <cit> .
most unexpected, however, was the demonstration that protein-coding transcripts are heavily edited at the rna level in diverse dinoflagellates  <cit> , unlike the case in either apicomplexans or ciliates.
to gain greater insight into the nature of dinoflagellate mitochondrial genomes, we have generated a large body of mitochondrial genomic and transcriptional data for two distantly related dinoflagellate species, c. cohnii and karlodinium micrum.
these data encompass more than  <dig> mtdna fragments totaling >  <dig> kb, and more than  <dig> mitochondrial transcripts.
this new information highlights several novel features of the organization and expression of the dinoflagellate mitochondrial genome, and concurrent studies in two additional distantly related dinoflagellates, amphidinium carterae  <cit>  and oxyrrhis marina  <cit> , corroborate a number of our findings.
together, these data reinforce the conclusion that the dinoflagellate mitochondrial genome has been substantially reorganized since the divergence of dinoflagellates and apicomplexans from a common ancestor.
RESULTS
genomic sequence reveals a complex mitochondrial genome
crypthecodinium cohnii
previously reported c. cohnii cox <dig> sequences indicated multiple copies of the gene with different flanking sequences  <cit> .
to test if this genomic complexity extends to other c. cohnii mitochondrial genes, we sequenced multiple genomic clones containing cob and/or cox <dig>  a library of ecori restriction fragments constructed from a fraction enriched in mtdna was screened using a c. cohnii cob gene probe, obtained by pcr.
this screen recovered a cob clone linked to a 57-bp cox <dig> fragment, which itself was used to probe for cox3-containing clones.
in total,  <dig> clones were characterized , ranging in size from  <dig>  kb to  <dig>  kb .
end sequencing and restriction mapping identified six unique cob-containing clones, and three unique cox3-containing clones.
four clones were completely sequenced .
the largest clone, pc3# <dig>  , contains a complete or nearly complete cob gene , followed by three other identifiable sequences: a 49-bp stretch identical to a sequence previously found in a cox1-containing clone  <cit> ; a 113-bp cox <dig> segment; and a 99-bp large subunit  rrna sequence corresponding to mitochondrial lsug in apicomplexans  <cit> .
two additional cob clones were sequenced, pcb# <dig>  and pcb# <dig> .
both encode cob, but with different flanking sequences than in pc3# <dig> .
pcb# <dig> contains unique 3' sequence immediately after the cob repeat, whereas pcb# <dig> contains additional common sequence with pc3# <dig>  for ~ <dig> kb before unique sequence occurs .
amongst these clones, we observed two different 5'-flanking sequences and three different 3'-flanking sequences .
this arrangement recapitulates the organization of cox <dig> in c. cohnii mtdna  <cit> , i.e., a central repeat  containing most of the cob orf) flanked by different arrays of unique upstream and downstream sequences.
partial sequencing of the remaining clones revealed an additional unique 5'-flanking sequence  and one additional unique 3'-flanking sequence  in the immediate vicinity of the cob orf .
of the three cob-containing clones described above, only pcb# <dig> encodes a complete cytochrome b  protein .
pc3# <dig>  and pcb# <dig> share an alternative 3' sequence that predicts a cob c-terminal sequence lacking  <dig> amino acid residues compared with the pcb#2-predicted cob as well as the corresponding plasmodium falciparum cob.
this suggests that the pc3# <dig>  and pcb# <dig> cob orfs represent pseudogenes.
variable 3' coding sequences were also seen previously for c. cohnii cox <dig>  with some coding sequences also truncated compared to other dinoflagellate sequences  <cit> .
one cox3-containing clone  was also sequenced, but it was found not to encode an intact cox <dig> gene.
instead, this clone encoded  <dig> bp identical in sequence to the portion of pc3# <dig>  that included the 113-bp cox <dig> segment and the 49-bp cox <dig> sequence .
this clone was also flanked by unique sequences, providing further evidence that mitochondrial genes occur in multiple genomic contexts in c. cohnii.
to further investigate the arrangements and relative numbers of mtdna elements, southern hybridization analysis was performed using region-specific probes.
as shown in figure  <dig>  probes were generated specific to: the cob coding sequence ; two cob 3'-flanking regions ; the cox <dig> sequence ; and the rrna sequence lsug .
these probes were hybridized against a mtdna-enriched fraction hydrolyzed by ecori.
with the 'cob' probe, a strong signal was detected at  <dig>  kb and weaker signals at  <dig> ,  <dig> ,  <dig> , and  <dig>  kb .
this result is consistent with dominant ecori clones being  <dig>  kb, and with multiple genomic contexts for cob.
probing with 3' flanking sequence 'cb1' revealed a similar banding pattern to that generated by the 'cob' probe, indicating that this region is typically contiguous with the cob coding sequence.
probing with 'cb3' presented a very different profile, with  <dig> bands ranging in size from  <dig>  to  <dig>  kb and of varying intensity .
the cb <dig> sequence evidently occurs in numerous ecori fragments, some without cob.
probing with 'cox3' and 'rnl' also revealed multiple bands with varying intensity , again indicating that these mtdna elements are present in several different genomic arrangements.
together these southern data verify the existence of multiple copies of c. cohnii mtdna elements occurring in different contexts, and indicate that up to  <dig> different arrangements occur for some of these elements.
karlodinium micrum
putative mitochondrial genes were identified from a survey of  <dig> k. micrum expressed sequence tag  sequences assembled into  <dig> unique clusters  <cit> .
oligoadenylation of mitochondrial gene transcripts is known from other organisms  <cit> , and this also appears to be the case in dinoflagellates as the poly-dependent k. micrum survey also contained many cdnas for mitochondrial genes.
mitochondrial sequences were identified by homology to genes in other systems, and all such cdnas were fully sequenced.
using this strategy we identified sequences representing the three protein-encoding genes found in c. cohnii: cox <dig> , cob  and cox <dig> .
the average a+t content of these sequences was 69% , consistent with their being encoded in the mitochondrion.
we found no other mitochondrial protein-coding sequences exhibiting the strong a+t biases suggestive of an origin from mtdna .
several short cdna sequences, however, with high similarity to the fragmented apicomplexan mitochondrial rrnas  <cit>   were identified.
these correspond to apicomplexan lsu rrna fragments lsua, rna <dig>  lsue, lsug and rna <dig> , small subunit  rrna fragment rna <dig> , and an rna  that has yet to be assigned to either the lsu or ssu rrna.
while these sequences have a lesser a+t bias  compared with the mitochondrial protein-encoding sequences, the high similarity of these sequences to their apicomplexan counterparts , and known oligoadenylation of these transcripts in apicomplexans  <cit> , strongly implicates these sequences as additional elements of the k. micrum mtdna.
with these  <dig> mtdna tags, we used pcr to generate genomic sequences corresponding to each gene and regions linking them, with the aim of assembling large portions of k. micrum mtdna sequence.
intergenic sequence recovered by this approach was used to provide further priming sites to extend the sampling of k. micrum mtdna.
in addition to amplification of individual genes, a total of  <dig> distinct gene linkage products were generated and fully sequenced .
this analysis yielded a sequence in which mitochondrial genes were linked to one another in many different contexts.
gene fragments were also common, as were mtdnas with three or four distinct fragments or tandem repeats .
in total, cob sequences were found in at least six mutually exclusive linkages, cox <dig> in five, cox <dig> in four, lsue in nine, rna <dig> in six, rna <dig> in five and rna <dig> in one.
additionally, two large cdnas  provided further evidence of multiple copies of mitochondrial genes and gene fragments linked in novel arrangements.
ef <dig>  for example, contains the lsug coding sequence, a second partial lsug unit within a 170-bp repeat, the lsua sequence, the rna <dig> sequence, and an internal fragment of the cox <dig> gene .
these cdnas also indicate that polycistronic transcription occurs in dinoflagellate mitochondria.
intergenic sequences from the pcr clones were examined for additional coding elements by comparison to publicly available databases, specifically searching against k. micrum ests as well as comparing the intergenic regions to one another.
no identifiable genes were found, but one cdna sequence  was represented in one mtdna clone, implicating this sequence as an additional transcriptional unit of the mitochondrial genome .
comparison of intergenic sequences to one another revealed numerous dispersed repeated sequences with either 100% or very high degrees of identity .
overall, data from k. micrum are consistent with those from c. cohnii, both pointing to a complex genome organization evidently underpinned by a high level of recombination within dinoflagellate mitochondria.
inverted repeats in mtdna
previous analysis of c. cohnii cox <dig> identified many short inverted repeats in flanking, non-coding sequences  <cit> .
we have applied a similar analysis to the c. cohnii cob- and cox3-containing sequences, as well as the k. micrum mtdna data, and find a very similar pattern of repeat features, although we also note some differences between the two taxa.
within the c. cohnii sequences, we screened for inverted repeats of different length and distance between them, and found two distinct but prevalent classes of this element type.
the first class is similar to those previously described  <cit> , and consists of very closely spaced, small inverted repeats .
these inverted repeats occur almost exclusively within non-coding sequence, with the only exceptions being at the very extremities of genes .
a second class of inverted repeats consists of longer repeat elements  no more than  <dig> nucleotides apart.
such inverted repeats are also prevalent in c. cohnii mtdna, and are almost exclusively features of the non-coding sequences .
analysis of k. micrum mtdna showed that inverted repeats are also a feature of intergenic sequences; however, in this case only the larger class of inverted repeats was found, with none of the smaller, closely spaced inverted repeats occurring in any of the mtdna sequences .
again these repeats are almost exclusively located within intergenic regions, with genic inverted repeats only occasionally present, within gene extremities.
no equivalent inverted repeats were found in a random sample of  <dig> k. micrum nucleus-encoded gene sequences .
the sequences of repeated elements in both c. cohnii and k. micrum are consistent with secondary structures such as stem loops and hairpins, and in both cases the repeated elements that could form such stem structures are typically g+c rich, in spite of the a+t bias of these organelle genomes.
the inverted repeats described here are also distinct from secondary structural elements of the rrnas  that typically consist of imperfect inverted repeats.
densely packed inverted repeats, primarily in intergenic regions, was also recently described from a. carterae mtdna  <cit> .
in this case, imperfect inverted repeats were predicted to form stems of 50– <dig> nucleotides, with at-rich loops of ~10– <dig> nucleotides.
while inverted repeats therefore appear to be a consistent feature of dinoflagellate mitochondrial genomes, the elaboration of these elements is variable between taxa, with shorter repeats only present in c. cohnii.
mitochondrial gene transcripts lack stop and start codons
extensive substitutional rna editing of transcripts occurs in dinoflagellate mitochondria, so exactly where an open reading frame begins and ends can only be tentatively inferred from genomic dna.
accordingly we used k. micrum cdnas, and publicly available mrna sequences from several other dinoflagellates, to identify the ends of all three protein-coding genes.
absence of stop codons
oligoadenylation of transcripts apparently occurs upstream of any canonical stop codon in all protein-encoding transcripts analyzed, and for only one gene does oligoadenylation create an in-frame canonical stop codon.
this lack of encoded stop codons applies to transcripts for cob, cox <dig> and cox <dig> represented from multiple species.
all  <dig> cob transcripts from k. micrum are oligoadenylated at the same point, which corresponds to the expected c-terminus of cob homologues , but does not include an in-frame stop.
the 3' ends of transcripts from four other dinoflagellates  are oligoadenylated at precisely the same position .
for cox <dig>  the mrna sequences from four taxa  are all oligoadenylated at the same position, where the protein sequence is predicted to terminate ; once again, none of these encode a stop codon.
the k. micrum cox <dig> cdnas present an even more interesting situation.
five of nine cdnas are oligoadenylated approximately  <dig> codons upstream of the predicted c-terminus, and without an in-frame stop codon .
however, another four cdnas are oligoadenylated a further  <dig> nucleotides downstream; these cdnas encode amino acid sequence with high similarity to the c-terminus of cox <dig>  in this case, oligoadenylation follows a u residue creating an in-frame uaa stop codon.
the generation of an in-frame stop codon concomitant with oligoadenylation is also apparent in amphidinium cox <dig> mrna; however, as in k. micrum, other cox <dig> amphidinium transcripts are oligoadenylated prematurely, within a few bases of the premature oligoadenylation site in k. micrum cdnas .
alternative oligoadenylation sites have also been reported for cox <dig> transcripts in the dinoflagellate g. polyedra  <cit> .
a potential alternative stop codon was sought among these transcript data by looking for a codon that occurs exclusively in the 3' region of these coding sequences.
however, no such candidate codon could be identified either within or between the taxa surveyed, nor is there any evidence for use of a non-standard genetic code .
moreover, oligoadenylation consistently occurred at the position where the protein sequence is expected to terminate, leaving little or no apparent untranslated region .
alternative start codons
dependence on a standard atg start codon also is apparently relaxed in dinoflagellate mitochondria.
from multiple dinoflagellate species mrnas for the three protein-coding genes extend beyond conserved n-termini, suggesting these transcripts are likely to be full length, but all lack a plausible n-terminal aug .
existing genomic sequences corroborate the lack of initiating atgs.
transcript data for cox <dig> from three species  and cox <dig> from k. micrum are all apparently full length based on protein alignments and all lack an aug in the terminal region .
the corresponding genomic region upstream of k. micrum cox <dig> does not contain an in-frame atg until  <dig> nucleotides upstream of the conserved sequence, and  <dig> stop codons fall between them, supporting the likely absence of an atg from this gene.
genomic sequences for c. cohnii cox <dig>  however, do contain an in-frame atg ~ <dig> codons upstream of n-terminal sequence conservation seen among dinoflagellates.
while it is possible that this particular atg serves as the initiator codon in this taxon, the lack of any sequence conservation with the corresponding k. micrum sequence within this 13-residue stretch  suggests that this might also represent a chance atg within the 5' utr.
k. micrum cob mrnas do encode an aug close to the site where sequence conservation with other cob proteins begins, but on close inspection there is conserved sequence upstream of this codon .
further, cob from the early-diverging member of the dinoflagellates, oxyrrhis marina, lacks this aug or any other upstream of this region  <cit> .
in mrnas of all other available species  there is strong conservation of the four predicted amino acid residues upstream of this atg , further suggesting that translation likely initiates upstream of it .
the conservative change of this second residue, v to l, among dinoflagellate taxa  supports the inference that this region likely represents protein-coding sequence rather than utr.
some conservation of this sequence with plasmodium cob is also apparent .
none of the four apparently full-length k. micrum cob genomic sequences encodes an additional atg codon between this region of conservation and the next in-frame stop codon , and the same situation is seen in a p. piscicida cob sequence.
the c. cohnii genomic sequences are the only cases to date where potential atg codons do occur in this upstream sequence .
however, two of these occur well upstream of any 5'-sequence conservation among dinoflagellates, and would represent unusually long  and divergent cob proteins in these cases .
trans-splicing of cox3
included among the k. micrum cox <dig> cdnas were four inferred to be full length  based on protein alignments , and five inferred to be prematurely oligoadenylated at nucleotide  <dig>  despite the fact that the longer cdna is likely the functional cox <dig> mrna, a genomic copy corresponding to it could not be amplified from genomic dna using multiple primer combinations .
the longest product obtained from genomic dna corresponded to nucleotides 50– <dig> of the full-length cox <dig> sequence.
six genomic fragments containing cox <dig> sequence were obtained by amplifying between genes, and these suggest that the gene is fragmented in the genome .
notably, three unique cox <dig> genomic sequences are truncated at nucleotide  <dig>  precisely where the short cox <dig> transcripts are oligoadenylated .
immediately downstream is a stop codon, and subsequently no further sequence similarity to cox <dig>  similarly, the only genomic sequences found to encode the 3' end of the long transcript are 5'-truncated at nucleotide  <dig>  with sequence unrelated to cox <dig> upstream of this point .
taken together, these data suggest that the long cox <dig> transcript is the product of trans-splicing, where nucleotides 1– <dig> are joined to nucleotides 718– <dig> arising from two different genomic fragments.
the intervening five nucleotides  are all a residues in the full-length cox <dig> transcript, suggesting that trans-splicing occurs within the oligo tail of the upstream transcript.
mitochondrial rrnas are fragmented in a similar pattern as in apicomplexans
ssu and lsu rrnas are encoded in all characterized mtdnas; however, until recently  <cit>  no mitochondrial rrna sequences had been described from dinoflagellates.
in this study we have identified several discrete, short sequences with strong similarity to components of the highly fragmented rrnas of apicomplexans  <cit>  .
from k. micrum, we obtained cdna sequences representing five lsu rrna fragments , one ssu rrna fragment , and one unassigned rrna fragment , all of which correspond to known transcriptional units of the plasmodium mitochondrial genome.
we also identified an additional lsu rrna fragment, lsuf, as well as lsue and lsug, from an est survey we previously conducted in heterocapsa triquetra  <cit> .
alignment of lsua, lsue, lsuf, lsug and rna <dig> to their plasmodium lsu homologues is shown in figure  <dig>  ssu rrna fragment rna <dig> and unassigned fragment rna <dig> share 66% and 74% sequence identity to plasmodium homologues, respectively.
for each fragment, multiple cdnas were sequenced , and oligoadenylation was found to occur at a consistent site .
although these cdnas are all relatively short, the 5' ends could not be definitively determined from these cdnas because the 5'-lengths were variable.
further, genomic copies  encoded conserved sequence upstream of the 5' ends of cdnas of lsue and lsug .
for c. cohnii, the lsug sequence identified on ecori clone pc3# <dig>  was analyzed by 3' race and the site of oligoadenylation was shown to be identical to that in the corresponding k. micrum and h. triquetra cdnas .
northern analysis of c. cohnii rna showed a single lsug-positive band at ~ <dig> nucleotides  <cit> .
this size corresponds well with the limit of conservation among lsu rrna sequences, as well as the size of the plasmodium lsug.
c. cohnii lsue was also amplified and the ends determined by 5'-cdna sequencing and 3' race .
northern hybridization against mitochondrial rna confirmed the presence of an ~ <dig> nucleotide rna species  <cit> .
the oligoadenylation sites for mitochondrial rrna fragments are identical among dinoflagellates, and either identical or within a few nucleotides of those observed in plasmodium .
the 5' ends of these sequences, whether defined experimentally  or by sequence conservation, are also very similar to those of their plasmodium counterparts.
the only possible exception is k. micrum rna <dig>  where the sole cdna obtained contained substantial upstream  and downstream  sequence compared to the region with similarity to plasmodium rna <dig>  however it is possible that this cdna represents an unprocessed precursor, and accordingly further work is required to substantiate the size of this putative rrna fragment.
secondary structure predictions for dinoflagellate sequences lsua, lsue, lsuf, lsug, rna <dig> and putative rna <dig>  all indicate that the expected folding and intermolecular base pairings occur , and these fragments are likely to contribute to a viable reconstituted lsu rrna, as for plasmodium.
rna editing
protein-coding genes
rna editing has been described for cox <dig>  cob and cox <dig> transcripts from diverse dinoflagellates, including the cob mrna of k. micrum  <cit> .
comparison of k. micrum cdna and corresponding mtdna sequences for the three genes identified here confirms this conclusion for transcripts of cob, and further shows that cox <dig> and cox <dig> transcripts are also edited.
the average density of editing of the cox <dig> transcripts is one substitution per  <dig> nucleotides and this value is consistent with other studies in different species  <cit> .
by contrast, editing in cox <dig> transcripts is over twice as dense, at one substitution per  <dig> nucleotides, making cox <dig> the most heavily edited gene transcript in dinoflagellates.
editing of cob mrna lies in between these extremes, at one substitution per  <dig> nucleotides.
in the case of cox <dig> transcripts, four types of substitutional changes were detected at  <dig> sites.
of these, 48% were a to g substitutions, followed by u to c  and smaller proportions of c to u and g to c edits .
this observation is consistent with cox <dig> mrna editing occurring in other species, where most  of the reported changes are a to g and u to c substitutions  <cit> .
so far, g to c changes have only been observed in mtdna-encoded mrnas of dinoflagellates, whereas a to g changes have only been reported in nucleus-encoded mrnas.
cox <dig> mrna editing types are generally consistent with those observed in cox <dig> and cob mrnas.
five types of substitutional changes were observed at  <dig> sites, of which 42% were a to g changes, followed by c to u and u to c edits , as well as three g to a edits  and a single g to c edit .
for both cox <dig> and cox <dig> mrnas, the majority of substitutions occur at the first or second positions of affected codons , and over 90% of editing events result in a change in predicted amino acid.
in k. micrum cox <dig> mrna , editing also removes a uag codon, which is typically a stop codon but is apparently unassigned in dinoflagellates.
analysis of the  <dig> cdnas corresponding to cox <dig> and cob offers further insight into the process of rna editing in dinoflagellates.
despite overall uniformity of transcript editing, some cdnas exhibit pre-edited states.
k. micrum cox <dig> and cob contain  <dig> and  <dig> editing sites, respectively, with the cdnas analyzed here representing in total  <dig> and  <dig> potential editing events, respectively.
however at nine of these sites in the cox <dig> cdnas, and five in the cob cdnas, the pre-edited nucleotide occurs, indicating  <dig> % and  <dig> % 'non-edits', respectively.
these 'non-edits' were present in only a few cdnas , suggesting that the great majority of cdnas represent mature transcripts.
the pre-edited sites are scattered throughout the transcripts where they are found, occur between other edited sites, and in no obvious order in any sequence.
these pre-edited sites may indicate editing failures, in which case such transcripts could give rise to defective translation products.
alternatively, they may represent editing intermediates.
if the latter is the case these data suggest that editing does not occur in a linear sequence along each transcript.
pre-edited mitochondrial cdnas have also recently been found in a. carterae mtdna  <cit> .
rrna transcripts
comparisons of rrna cdnas to genomic sequences are constrained by the smaller sizes of these sequences , in particular where pcr has been used to amplify genomic sequence a greater portion of this sequence represented primer binding sites and therefore cannot be used in such a comparison.
nevertheless, from the available data, there is no evidence of editing of rna <dig>  rna <dig> or rna <dig>  for lsue, complete genomic sequence  was available from the internal regions of five pcr fragments, with the majority of the sequence available from a further four pcr products using lsue primers.
these sequences were identical to the cdnas except for three consecutive nucleotides that were absent in two of the three lsue cdnas obtained from the est survey.
to test this anomaly, a further five cdnas were independently generated, and these all contained the three nucleotides, and therefore were identical to genomic lsue sequences and to one of the original est sequences.
these results suggest that the three-nucleotide deletions seen in two cdnas represent a rare artifact, likely generated during reverse transcription, and that k. micrum lsue is likely also not edited.
there was, however, evidence of substitutional editing for lsua and lsug.
in both cases genomic copies of these sequences differed from transcripts: in lsug at eight positions and in lsua at six positions .
consistent with the protein-coding genes, these substitutions consist mainly of a to g , c to u  and u to c  substitutions, with one case of c to g. given that dinoflagellate mitochondrial genes occur in multiple copies, recovery of further, independently isolated copies of these genes will be required to substantiate these inferences of rrna editing.
evidence for rrna editing has also recently been reported with the dinoflagellate a. catenella, where two inferred editing events were identified for the 'lsue-like' rrna  <cit> .
discussion
prior to this study our view of the dinoflagellate mitochondrial genome was gleaned from relatively sparse molecular data obtained from several diverse dinoflagellate taxa.
these data nevertheless provided a tantalizing view of a mitochondrial genome displaying several eccentricities.
coding sequences for entire or partial versions of cox <dig>  cob or cox <dig> have been shown to occur in multiple copies and in different genomic contexts in c. cohnii  <cit> , g. polyedra  <cit> , p. piscicida  <cit> , and a. catenella  <cit> .
these data paint a picture of dinoflagellate genomes in sharp contrast to the minimalist  <dig> kb apicomplexan mtdna, which encodes single copies of these genes, tightly packed together  <cit> .
similarly, extensive rna editing has been described in mrnas from diverse dinoflagellates  <cit> , a process that does not occur in apicomplexans.
in this study we have generated a much more comprehensive body of mitochondrial genomic and transcript data for two dinoflagellate species, c. cohnii and k. micrum, and these data are bolstered by a concurrent mitochondrial genomic study of the dinoflagellate a. carterae  <cit> .
together, these results reinforce the view that the dinoflagellate mitochondrial genome has diverged radically in form from that of apicomplexans, despite the persistence of some intriguing similarities.
mitochondrial genome content and form
compared to the complement of  <dig> to  <dig> genes in the mitochondrial genome of ciliates  <cit> , the most basal member of the phylum alveolata, the very low information content of apicomplexan mtdna  clearly shows that there has been considerable mitochondrial gene loss and/or relocation to the nucleus during alveolate evolution.
we infer that much, if not all, of this gene relocation must have occurred prior to the last common ancestor of dinoflagellates and apicomplexans.
in est surveys, we have only identified the same three protein-coding genes ; moreover, we found no other mitochondrial orfs of known function in >  <dig> kb and >  <dig> kb of mtdna sequence from k. micrum and c. cohnii, respectively.
these findings are consistent with the previous demonstration that cox <dig>  an otherwise nearly ubiquitous component of mitochondrial genomes, has been relocated to the nucleus in both apicomplexans and dinoflagellates  <cit> .
the only additional genes we identified are ones representing the mitochondrial ssu and lsu rrnas, which together with cox <dig> and cob are universally present in mtdna.
no trna genes have been found linked to mtdna sequences, similar to apicomplexans, where trnas are apparently imported into mitochondria from the cytoplasm  <cit> .
dinoflagellates and apicomplexans also share the characteristic of highly fragmented ssu and lsu rrnas.
fragmentation of mitochondrial rrna genes has been documented in the mitochondrial genomes of several eukaryotes, including ciliates  <cit> , several green algae  <cit>  and a fungus  <cit> .
the degree of fragmentation in apicomplexan mitochondrial rrna is more extreme than in these other cases, with  <dig> fragments for the ssu and lsu rrnas reported to date, coding regions for which are rearranged and interspersed with other genes in the genome  <cit> .
from within three disparate dinoflagellate taxa we have identified eight rrna fragments similar to fragments in p. falciparum, and three of these rrna species have also recently been reported from two further taxa, a. catenella and o. marina  <cit> .
the dinoflagellate rrna fragments mostly appear to correspond to their plasmodium counterparts in length and sequence termini, suggesting that a stable level and pattern of fragmentation has been inherited from the common ancestor of dinoflagellates and apicomplexans.
given that ciliate mitochondrial rrnas are comparatively intact , the extreme fragmentation in dinoflagellates and apicomplexans must have occurred since their divergence from ciliates.
despite a similar gene content the arrangement of dinoflagellate and apicomplexan mitochondrial genomes is radically different.
where the apicomplexan genome is relatively simple and compact, the dinoflagellate mitochondrial genome is complex, with multiple copies of each gene imbedded within different genomic contexts.
gene fragments and non-coding regions are also repeated, altogether suggesting a great deal of recombination in the genome, which is also consistent with the lack of sequence divergence among the multiple copies of these elements.
shotgun sequence data recently published for the a. carterae mitochondrion corroborate this picture of a recombining complex genome, and further suggest that the majority of the mitochondrial genome  might be non-coding  <cit> .
gene expression in dinoflagellate mitochondria
within the k. micrum est survey, long cdnas that encoded several mitochondrial genes or gene fragments  were noted.
by contrast, most mitochondrial cdnas we recovered encoded a single gene, suggesting the longer transcripts may be rapidly processed into shorter molecules.
polycistronic transcripts up to  <dig>  kb are also known from apicomplexan mtdna, these are rapidly processed to short, single-gene transcripts  <cit> .
interestingly, the polycistronic transcripts from k. micrum are not edited, indicating that rna editing acts on the individual gene transcripts.
the use of alternative initiation codons in dinoflagellate mitochondrial genes is consistent with what is seen in the mitochondria of other alveolates.
in plasmodium species, cox <dig> and cox <dig> lack an in-frame atg, and while cob does contain a atg near the initiation site, it is uncertain whether initiation occurs at this site or upstream of it  <cit>  .
att and ata have been proposed as alternative initiator codons in plasmodium species  <cit>  .
several mitochondrial genes from the ciliate tetrahymena pyriformis also apparently use alternative initiation codons of the form atn or ntg: in the case of cob an atg within eight codons of the predicted n-terminus is apparently ignored, with gtg used in its place  <cit> .
thus, there are precedents for reliance on codons other than atg for translation initiation within alveolates.
potential initiator atn/ntg codons exist in all three karlodinium mitochondrial genes; however, a broader survey of dinoflagellates or analysis of protein sequences will be necessary to identify the most likely candidates.
an absence of stop codons is more unusual.
in t. pyriformis all mitochondrial protein-coding genes terminate with taa  <cit> .
tga encodes tryptophan  and tag is simply not used.
all three plasmodium mitochondrial protein-coding genes also use taa  <cit> .
by contrast many dinoflagellate mitochondrial gene transcripts appear to lack any termination codon.
with only a single known exception , transcripts are oligoadenylated upstream of any of the standard termination codons, and rna editing does not generate an in-frame stop.
further, in none of the transcripts is a sense codon uniquely localized in the 3' region in such a way as to suggest that it serves as an alternative terminator .
the oligoadenylation of k. micrum cox <dig> mrna does produces a uaa codon, as is also the case for cox <dig> transcripts for a. carterae and o. marina  <cit> , that suggests that cox <dig>  unlike cox <dig> or cob, might utilize conventional stop.
such a mechanism for reconstituting a functional uaa is known to occur in some mammalian mitochondrial transcripts  <cit> .
it is unclear how the mitochondrial translation machinery might cope with the absence of termination codons.
release factors that are essential for disassembly of the ribosome usually recognize specific codons, so the absence of these codons could block ribosome disassembly.
there are precedents in other mitochondrial systems for the lack of termination codons: transcripts of two plant mitochondrial genes have been shown to be oligoadenylated upstream of in-frame stops  <cit> .
proteins encoded by both of these genes can be detected, indicating that the corresponding transcripts are successfully translated.
in human mitochondria, a rare mutation has been shown to ablate a stop codon, and yet the corresponding protein is still detectable is these cell lines  <cit> .
eubacteria are known to be able to rescue damaged mrna molecules that have lost their termination codon by use of a specialized rna with properties of both a trna and an mrna  <cit> .
these so-called tmrnas restart protein synthesis by providing a terminal mrna section that encodes a functional stop codon.
it has been speculated that an equivalent system might be used in plant and animal mitochondrial systems where mrnas lack stop codons  <cit> .
indeed, tmrna-like rna species have been identified in the mitochondria of jakobid flagellates such as r. americana; however, these rnas lack the terminal mrna-like segment of a conventional tmrna  <cit> .
moreover, the c-terminal tag provided by a tmrna normally targets the modified protein for degradation rather than for function  <cit> .
whatever the actual mechanism of translation termination in dinoflagellate mitochondria, it appears to present a clear difference with respect to protein synthesis termination in ciliate and apicomplexan mitochondria.
lastly, we have found a likely case of trans-splicing of dinoflagellate mitochondrial transcripts, which adds a further layer of complexity to genome organization and expression in these organelles.
while we cannot conclusively eliminate the possibility of a complete cox <dig> coding sequence in dinoflagellates we have not been able to detect an intact gene.
this negative result is consistent with all other studies to date, which report only partial cox <dig> sequences from five different dinoflagellate taxa  <cit>  .
all available data from genomic fragments and transcripts suggest that the complete cox <dig> transcript is generated by trans-splicing.
such trans-splicing has not been reported for either apicomplexan or ciliate mitochondria.
in ciliates nad <dig> is split into two segments  <cit>  but they are independently transcribed, and there is no evidence of splicing of the corresponding transcripts to create a continuous, complete nad <dig> orf  <cit> .
trans-splicing occurs in plant mitochondria  <cit> , but in these cases the coding breakpoints are flanked by group ii intron elements, which form secondary structures that mediate the splicing events.
we have no evidence of group ii introns in dinoflagellate mtdna, but we do note that the intergenic sequences contain numerous inverted repeats consistent with extensive secondary structure, which might conceivably facilitate splicing events.
the unique nature of the dinoflagellate trans-splicing is also evident from the inclusion of five a residues at the splice boundary that appear to derive from the oligo tail of the upstream fragment.
the removal of any downstream sequence by oligoadenylation prior to splicing argues against the involvement of a cis-acting element such as a group ii intron in the splicing process.
it is conceivable that oligoadenylation of the short 5' cox <dig> transcript could serve as a degradation signal for these short transcripts, as has been observed in human mitochondria  <cit> .
however, lack of a complete cox <dig> coding sequence, coupled with the fact that the site of oligoadenylation corresponds with the break in coding sequence of 5' and 3' cox <dig> portions, suggests that the short cox <dig> transcripts are important intermediates in the generation of the complete cox <dig> transcripts.
rna editing
the rna editing observed in k. micrum cox <dig>  cob and cox <dig> mrnas is consistent with the level and type of editing observed in cox <dig> and cob mrnas in other dinoflagellate species  <cit> , with the exception that cox <dig> is even more heavily edited than either of cox <dig> or cob.
while some editing sites are conserved, others are unique to certain taxa, suggesting that new editing sites are constantly evolving in dinoflagellates.
in this study we also found evidence in k. micrum of editing of rrna fragments lsug and lsua.
rna editing of a. catenella lsue has also recently been reported  <cit> .
at present the data are insufficient to assess the conservation of rrna editing sites among taxa; however, two inferred editing sites in a. catenella lsue are not edited in k. micrum, suggesting that rrna editing sites are constantly evolving as with those in protein-coding genes.
whether rna editing plays some functional role in dinoflagellate mitochondria is unclear.
from analysis of protein-coding genes in several dinoflagellates, lin et al  <cit>  noted that the majority of editing events are to either a c or g, thus generating a net reduction in a+u content from the bias of ~70% for the coding sequences.
we observe this trend also in karlodinium protein-coding sequences.
this re-tailoring of mrna sequences might better accommodate the suite of nucleus-encoded trnas that are likely imported from the cytoplasm, and which typically participate in the decoding of nucleus-encoded mrnas having a more balanced a+u content  <cit> .
ribosomal rna is also sensitive to a+u content, with secondary structure elements such as hairpin loops better stabilized by g-c than by a-u pairs; thus, helical regions tend to be relatively more g+c rich than other rrna domains.
while the available data for rrna editing are limited , it is interesting that the editing types in rrnas have an overall neutral impact on a+u content.
indeed the a+t content of mitochondrial genomic sequence specifying rrnas is already much reduced  compared to that of the protein-coding genes.
this observation might add weight to the notion that editing helps correct  the a+t skew of protein-coding genes.
the mechanism of rna editing in dinoflagellate mitochondria is also unknown; however the possibility of a guide rna -assisted mechanism, similar to that employed in trypanosomatid mitochondria  <cit> , has recently been suggested  <cit> .
nash et al  <cit>  report that gene fragments encoded in mitochondria sometimes encode the 'corrected' nucleotide at an inferred editing site .
thus such fragments could encode templates that direct the editing events of full-length transcripts.
we analyzed the k. micrum data for similar evidence of post-edited nucleotides represented in gene fragments.
from five fragments  that span  <dig> editing sites across the three protein-encoding genes, only one site in one of the fragments corresponds to a 'corrected' nucleotide seen in cdnas at an inferred editing site .
an independent copy of the cob genomic sequence verified that this nucleotide difference is genuine .
hence this might represent an example of an editing template in k. micrum; however, if grnas are responsible for all editing events, a very large number of additional fragments must exist to direct the remainder of the changes.
clearly further work is required to shed light on the mechanism of rna editing in dinoflagellates.
future directions
a key question that remains is whether the observed diversity of dinoflagellate mitochondrial genes, gene fragments, and repetitive elements derives from a single mtdna molecule or from multiple chromosomes.
a similar scenario of mitochondrial genes occurring as multiple copies and fragments is seen in the ichthyosporean a. parasiticum, a unicellular organism closely related to animals  <cit> .
in this protist, several hundred small linear chromosomes constitute the mitochondrial genome, each encoding a smattering of genes and partial genes.
diplonemids, members of the phylum euglenozoa, also contain fragmented genes on separate circular mitochondrial chromosomes  <cit> .
it is unknown whether either of these unusual situations applies to the organization of the dinoflagellate mitochondrial genome; however, in this regard we make two preliminary observations.
one is that long-range pcr was unable to generate longer contiguous sequences linking the many mtdna elements we report in this study.
rather, additional short unique gene linkages were obtained, and it is clear that we have yet to sample the full diversity of gene combinations.
secondly, the presence of individual genes in partial tandem repeats  is consistent with minicircles, as seen in dinoflagellate plastid genomes  <cit> .
if these cases represent true minicircles, we have been unable to amplify a corresponding sequence to close these circles .
it is also possible, of course, that the tandem repeats that we observe are simply a consequence of further recombination events, and the high diversity of gene combinations.
CONCLUSIONS
a greater depth of sampling of dinoflagellate mitochondrial dna and mrna has provided a clearer view of a complex genome and many peculiarities of gene expression.
we find that the dinoflagellate mitochondrial genome shares several features in common with the mtdna of its apicomplexan sister lineage, but also many novel characteristics.
features in common for the two lineages are:  a very high level of gene relocation from the mitochondrion,  extensive rrna gene fragmentation and dispersal, and  use of non-standard initiation codons.
features unique to dinoflagellates are:  gene copy number expansion and reorganization,  loss of stop codons from protein-coding genes,  mrna trans-splicing, and  rna editing of protein-coding and rrna transcripts.
these data demonstrate a remarkable burst of organelle genome evolution in dinoflagellates following divergence from apicomplexa, and also challenge our understanding of the mechanistic details of genome maintenance and expression, most notably translation termination.

 the origin of mitochondria by endosymbiosis has emerged as a pivotal event in the evolution of eukaryotes. all eukaryote groups that have been studied bear a derivative of this endosymbiont, and for most the resulting mitochondrion is central to energy metabolism as well as providing several other anabolic and catabolic functions. a relict, though functionally essential, mitochondrial genome persists in all but a few anaerobic eukaryotes, and the genes in these genomes firmly identify the original endosymbiont as an α-proteobacterium. the jakobid flagellate reclinomonas americana has the least derived mitochondrial genome characterized to date, with at least genes encoded on a single, circular-mapping kb chromosome. more typically mitochondrial genomes have been reduced to 40– genes arranged on either circular- or linear-mapping chromosomes of 15– kb. in some eukaryotic groups, however, the mtdna has been modified more substantially, resulting in extremes in genome structure. for example, trypanosomatid mtdna consists of a few dozen large circular molecules and several thousand minicircles that encode guide rnas that participate in extensive u insertion/deletion rna editing. diplonemid mitochondria also contain multiple circular mtdna molecules, each encoding gene fragments that are trans-spliced to generate functional transcripts. another example is the mtdna in the ichthyosporean, amoebidium parasiticum: in this case, mitochondrial genes are fragmented and dispersed over several hundred linear chromosomes, totaling > kb. over the diversity of eukaryotes, mitochondrial genomes exhibit other interesting characteristics, including the use of a number of different non-standard genetic codes, many of which involve alterations in start and, more rarely, stop codons. one large group in which particularly interesting mitochondrial genome variation has been found is alveolates. three major phyla make up alveolates: ciliates, apicomplexans, and dinoflagellates, with apicomplexans and dinoflagellates being sister clades to the exclusion of ciliates. within alveolates, ciliate mtdna is the most conventional, consisting of a linear molecule, 40– kb in length, that codes for many of the standard mitochondrial proteins found in other organisms. by contrast, the mtdna of the apicomplexan genus plasmodium is the smallest known, consisting of a linear, kb tandem repeat with only three protein-coding genes: cytochrome oxidase subunit, cytochrome oxidase subunit and cytochrome b. in addition, ciliate mtdna encodes two ribosomal rnas, but the corresponding apicomplexan genes are fragmented to an unprecedented degree and scattered about the genome. to date, dinoflagellate mtdnas have been the least well studied of alveolate mitochondrial genomes, with existing data pointing to a genome exhibiting several eccentricities. the first sequences isolated were four copies of cox from crypthecodinium cohnii, each of which was found to occur in a unique genomic context. southern blots demonstrated multiple different copies of this gene that varied in abundance, suggesting the c. cohnii mitochondrial genome is not as streamlined as in apicomplexans. subsequently, cob and cox have been found as well, and multiple, sometimes fragmented copies of these genes have now been reported from diverse dinoflagellates. most unexpected, however, was the demonstration that protein-coding transcripts are heavily edited at the rna level in diverse dinoflagellates, unlike the case in either apicomplexans or ciliates. to gain greater insight into the nature of dinoflagellate mitochondrial genomes, we have generated a large body of mitochondrial genomic and transcriptional data for two distantly related dinoflagellate species, c. cohnii and karlodinium micrum. these data encompass more than mtdna fragments totaling > kb, and more than mitochondrial transcripts. this new information highlights several novel features of the organization and expression of the dinoflagellate mitochondrial genome, and concurrent studies in two additional distantly related dinoflagellates, amphidinium carterae and oxyrrhis marina, corroborate a number of our findings. together, these data reinforce the conclusion that the dinoflagellate mitochondrial genome has been substantially reorganized since the divergence of dinoflagellates and apicomplexans from a common ancestor.  genomic sequence reveals a complex mitochondrial genome crypthecodinium cohnii previously reported c. cohnii cox sequences indicated multiple copies of the gene with different flanking sequences. to test if this genomic complexity extends to other c. cohnii mitochondrial genes, we sequenced multiple genomic clones containing cob and/or cox a library of ecori restriction fragments constructed from a fraction enriched in mtdna was screened using a c. cohnii cob gene probe, obtained by pcr. this screen recovered a cob clone linked to a 57-bp cox fragment, which itself was used to probe for cox3-containing clones. in total, clones were characterized, ranging in size from kb to kb. end sequencing and restriction mapping identified six unique cob-containing clones, and three unique cox3-containing clones. four clones were completely sequenced. the largest clone, pc3#, contains a complete or nearly complete cob gene, followed by three other identifiable sequences: a 49-bp stretch identical to a sequence previously found in a cox1-containing clone; a 113-bp cox segment; and a 99-bp large subunit rrna sequence corresponding to mitochondrial lsug in apicomplexans. two additional cob clones were sequenced, pcb# and pcb#. both encode cob, but with different flanking sequences than in pc3#. pcb# contains unique 3' sequence immediately after the cob repeat, whereas pcb# contains additional common sequence with pc3# for ~ kb before unique sequence occurs. amongst these clones, we observed two different 5'-flanking sequences and three different 3'-flanking sequences. this arrangement recapitulates the organization of cox in c. cohnii mtdna, i.e., a central repeat containing most of the cob orf) flanked by different arrays of unique upstream and downstream sequences. partial sequencing of the remaining clones revealed an additional unique 5'-flanking sequence and one additional unique 3'-flanking sequence in the immediate vicinity of the cob orf. of the three cob-containing clones described above, only pcb# encodes a complete cytochrome b protein. pc3# and pcb# share an alternative 3' sequence that predicts a cob c-terminal sequence lacking amino acid residues compared with the pcb#2-predicted cob as well as the corresponding plasmodium falciparum cob. this suggests that the pc3# and pcb# cob orfs represent pseudogenes. variable 3' coding sequences were also seen previously for c. cohnii cox with some coding sequences also truncated compared to other dinoflagellate sequences. one cox3-containing clone was also sequenced, but it was found not to encode an intact cox gene. instead, this clone encoded bp identical in sequence to the portion of pc3# that included the 113-bp cox segment and the 49-bp cox sequence. this clone was also flanked by unique sequences, providing further evidence that mitochondrial genes occur in multiple genomic contexts in c. cohnii. to further investigate the arrangements and relative numbers of mtdna elements, southern hybridization analysis was performed using region-specific probes. as shown in figure probes were generated specific to: the cob coding sequence; two cob 3'-flanking regions; the cox sequence; and the rrna sequence lsug. these probes were hybridized against a mtdna-enriched fraction hydrolyzed by ecori. with the'cob' probe, a strong signal was detected at kb and weaker signals at,,, and kb. this result is consistent with dominant ecori clones being kb, and with multiple genomic contexts for cob. probing with 3' flanking sequence'cb1' revealed a similar banding pattern to that generated by the'cob' probe, indicating that this region is typically contiguous with the cob coding sequence. probing with'cb3' presented a very different profile, with bands ranging in size from to kb and of varying intensity. the cb sequence evidently occurs in numerous ecori fragments, some without cob. probing with'cox3' and'rnl' also revealed multiple bands with varying intensity, again indicating that these mtdna elements are present in several different genomic arrangements. together these southern data verify the existence of multiple copies of c. cohnii mtdna elements occurring in different contexts, and indicate that up to different arrangements occur for some of these elements. karlodinium micrum putative mitochondrial genes were identified from a survey of k. micrum expressed sequence tag sequences assembled into unique clusters. oligoadenylation of mitochondrial gene transcripts is known from other organisms, and this also appears to be the case in dinoflagellates as the poly-dependent k. micrum survey also contained many cdnas for mitochondrial genes. mitochondrial sequences were identified by homology to genes in other systems, and all such cdnas were fully sequenced. using this strategy we identified sequences representing the three protein-encoding genes found in c. cohnii: cox, cob and cox. the average a+t content of these sequences was 69%, consistent with their being encoded in the mitochondrion. we found no other mitochondrial protein-coding sequences exhibiting the strong a+t biases suggestive of an origin from mtdna. several short cdna sequences, however, with high similarity to the fragmented apicomplexan mitochondrial rrnas were identified. these correspond to apicomplexan lsu rrna fragments lsua, rna lsue, lsug and rna, small subunit rrna fragment rna, and an rna that has yet to be assigned to either the lsu or ssu rrna. while these sequences have a lesser a+t bias compared with the mitochondrial protein-encoding sequences, the high similarity of these sequences to their apicomplexan counterparts, and known oligoadenylation of these transcripts in apicomplexans, strongly implicates these sequences as additional elements of the k. micrum mtdna. with these mtdna tags, we used pcr to generate genomic sequences corresponding to each gene and regions linking them, with the aim of assembling large portions of k. micrum mtdna sequence. intergenic sequence recovered by this approach was used to provide further priming sites to extend the sampling of k. micrum mtdna. in addition to amplification of individual genes, a total of distinct gene linkage products were generated and fully sequenced. this analysis yielded a sequence in which mitochondrial genes were linked to one another in many different contexts. gene fragments were also common, as were mtdnas with three or four distinct fragments or tandem repeats. in total, cob sequences were found in at least six mutually exclusive linkages, cox in five, cox in four, lsue in nine, rna in six, rna in five and rna in one. additionally, two large cdnas provided further evidence of multiple copies of mitochondrial genes and gene fragments linked in novel arrangements. ef for example, contains the lsug coding sequence, a second partial lsug unit within a 170-bp repeat, the lsua sequence, the rna sequence, and an internal fragment of the cox gene. these cdnas also indicate that polycistronic transcription occurs in dinoflagellate mitochondria. intergenic sequences from the pcr clones were examined for additional coding elements by comparison to publicly available databases, specifically searching against k. micrum ests as well as comparing the intergenic regions to one another. no identifiable genes were found, but one cdna sequence was represented in one mtdna clone, implicating this sequence as an additional transcriptional unit of the mitochondrial genome. comparison of intergenic sequences to one another revealed numerous dispersed repeated sequences with either 100% or very high degrees of identity. overall, data from k. micrum are consistent with those from c. cohnii, both pointing to a complex genome organization evidently underpinned by a high level of recombination within dinoflagellate mitochondria. inverted repeats in mtdna previous analysis of c. cohnii cox identified many short inverted repeats in flanking, non-coding sequences. we have applied a similar analysis to the c. cohnii cob- and cox3-containing sequences, as well as the k. micrum mtdna data, and find a very similar pattern of repeat features, although we also note some differences between the two taxa. within the c. cohnii sequences, we screened for inverted repeats of different length and distance between them, and found two distinct but prevalent classes of this element type. the first class is similar to those previously described, and consists of very closely spaced, small inverted repeats. these inverted repeats occur almost exclusively within non-coding sequence, with the only exceptions being at the very extremities of genes. a second class of inverted repeats consists of longer repeat elements no more than nucleotides apart. such inverted repeats are also prevalent in c. cohnii mtdna, and are almost exclusively features of the non-coding sequences. analysis of k. micrum mtdna showed that inverted repeats are also a feature of intergenic sequences; however, in this case only the larger class of inverted repeats was found, with none of the smaller, closely spaced inverted repeats occurring in any of the mtdna sequences. again these repeats are almost exclusively located within intergenic regions, with genic inverted repeats only occasionally present, within gene extremities. no equivalent inverted repeats were found in a random sample of k. micrum nucleus-encoded gene sequences. the sequences of repeated elements in both c. cohnii and k. micrum are consistent with secondary structures such as stem loops and hairpins, and in both cases the repeated elements that could form such stem structures are typically g+c rich, in spite of the a+t bias of these organelle genomes. the inverted repeats described here are also distinct from secondary structural elements of the rrnas that typically consist of imperfect inverted repeats. densely packed inverted repeats, primarily in intergenic regions, was also recently described from a. carterae mtdna. in this case, imperfect inverted repeats were predicted to form stems of 50– nucleotides, with at-rich loops of ~10– nucleotides. while inverted repeats therefore appear to be a consistent feature of dinoflagellate mitochondrial genomes, the elaboration of these elements is variable between taxa, with shorter repeats only present in c. cohnii. mitochondrial gene transcripts lack stop and start codons extensive substitutional rna editing of transcripts occurs in dinoflagellate mitochondria, so exactly where an open reading frame begins and ends can only be tentatively inferred from genomic dna. accordingly we used k. micrum cdnas, and publicly available mrna sequences from several other dinoflagellates, to identify the ends of all three protein-coding genes. absence of stop codons oligoadenylation of transcripts apparently occurs upstream of any canonical stop codon in all protein-encoding transcripts analyzed, and for only one gene does oligoadenylation create an in-frame canonical stop codon. this lack of encoded stop codons applies to transcripts for cob, cox and cox represented from multiple species. all cob transcripts from k. micrum are oligoadenylated at the same point, which corresponds to the expected c-terminus of cob homologues, but does not include an in-frame stop. the 3' ends of transcripts from four other dinoflagellates are oligoadenylated at precisely the same position. for cox the mrna sequences from four taxa are all oligoadenylated at the same position, where the protein sequence is predicted to terminate; once again, none of these encode a stop codon. the k. micrum cox cdnas present an even more interesting situation. five of nine cdnas are oligoadenylated approximately codons upstream of the predicted c-terminus, and without an in-frame stop codon. however, another four cdnas are oligoadenylated a further nucleotides downstream; these cdnas encode amino acid sequence with high similarity to the c-terminus of cox in this case, oligoadenylation follows a u residue creating an in-frame uaa stop codon. the generation of an in-frame stop codon concomitant with oligoadenylation is also apparent in amphidinium cox mrna; however, as in k. micrum, other cox amphidinium transcripts are oligoadenylated prematurely, within a few bases of the premature oligoadenylation site in k. micrum cdnas. alternative oligoadenylation sites have also been reported for cox transcripts in the dinoflagellate g. polyedra. a potential alternative stop codon was sought among these transcript data by looking for a codon that occurs exclusively in the 3' region of these coding sequences. however, no such candidate codon could be identified either within or between the taxa surveyed, nor is there any evidence for use of a non-standard genetic code. moreover, oligoadenylation consistently occurred at the position where the protein sequence is expected to terminate, leaving little or no apparent untranslated region. alternative start codons dependence on a standard atg start codon also is apparently relaxed in dinoflagellate mitochondria. from multiple dinoflagellate species mrnas for the three protein-coding genes extend beyond conserved n-termini, suggesting these transcripts are likely to be full length, but all lack a plausible n-terminal aug. existing genomic sequences corroborate the lack of initiating atgs. transcript data for cox from three species and cox from k. micrum are all apparently full length based on protein alignments and all lack an aug in the terminal region. the corresponding genomic region upstream of k. micrum cox does not contain an in-frame atg until nucleotides upstream of the conserved sequence, and stop codons fall between them, supporting the likely absence of an atg from this gene. genomic sequences for c. cohnii cox however, do contain an in-frame atg ~ codons upstream of n-terminal sequence conservation seen among dinoflagellates. while it is possible that this particular atg serves as the initiator codon in this taxon, the lack of any sequence conservation with the corresponding k. micrum sequence within this 13-residue stretch suggests that this might also represent a chance atg within the 5' utr. k. micrum cob mrnas do encode an aug close to the site where sequence conservation with other cob proteins begins, but on close inspection there is conserved sequence upstream of this codon. further, cob from the early-diverging member of the dinoflagellates, oxyrrhis marina, lacks this aug or any other upstream of this region. in mrnas of all other available species there is strong conservation of the four predicted amino acid residues upstream of this atg, further suggesting that translation likely initiates upstream of it. the conservative change of this second residue, v to l, among dinoflagellate taxa supports the inference that this region likely represents protein-coding sequence rather than utr. some conservation of this sequence with plasmodium cob is also apparent. none of the four apparently full-length k. micrum cob genomic sequences encodes an additional atg codon between this region of conservation and the next in-frame stop codon, and the same situation is seen in a p. piscicida cob sequence. the c. cohnii genomic sequences are the only cases to date where potential atg codons do occur in this upstream sequence. however, two of these occur well upstream of any 5'-sequence conservation among dinoflagellates, and would represent unusually long and divergent cob proteins in these cases. trans-splicing of cox3 included among the k. micrum cox cdnas were four inferred to be full length based on protein alignments, and five inferred to be prematurely oligoadenylated at nucleotide despite the fact that the longer cdna is likely the functional cox mrna, a genomic copy corresponding to it could not be amplified from genomic dna using multiple primer combinations. the longest product obtained from genomic dna corresponded to nucleotides 50– of the full-length cox sequence. six genomic fragments containing cox sequence were obtained by amplifying between genes, and these suggest that the gene is fragmented in the genome. notably, three unique cox genomic sequences are truncated at nucleotide precisely where the short cox transcripts are oligoadenylated. immediately downstream is a stop codon, and subsequently no further sequence similarity to cox similarly, the only genomic sequences found to encode the 3' end of the long transcript are 5'-truncated at nucleotide with sequence unrelated to cox upstream of this point. taken together, these data suggest that the long cox transcript is the product of trans-splicing, where nucleotides 1– are joined to nucleotides 718– arising from two different genomic fragments. the intervening five nucleotides are all a residues in the full-length cox transcript, suggesting that trans-splicing occurs within the oligo tail of the upstream transcript. mitochondrial rrnas are fragmented in a similar pattern as in apicomplexans ssu and lsu rrnas are encoded in all characterized mtdnas; however, until recently no mitochondrial rrna sequences had been described from dinoflagellates. in this study we have identified several discrete, short sequences with strong similarity to components of the highly fragmented rrnas of apicomplexans. from k. micrum, we obtained cdna sequences representing five lsu rrna fragments, one ssu rrna fragment, and one unassigned rrna fragment, all of which correspond to known transcriptional units of the plasmodium mitochondrial genome. we also identified an additional lsu rrna fragment, lsuf, as well as lsue and lsug, from an est survey we previously conducted in heterocapsa triquetra. alignment of lsua, lsue, lsuf, lsug and rna to their plasmodium lsu homologues is shown in figure ssu rrna fragment rna and unassigned fragment rna share 66% and 74% sequence identity to plasmodium homologues, respectively. for each fragment, multiple cdnas were sequenced, and oligoadenylation was found to occur at a consistent site. although these cdnas are all relatively short, the 5' ends could not be definitively determined from these cdnas because the 5'-lengths were variable. further, genomic copies encoded conserved sequence upstream of the 5' ends of cdnas of lsue and lsug. for c. cohnii, the lsug sequence identified on ecori clone pc3# was analyzed by 3' race and the site of oligoadenylation was shown to be identical to that in the corresponding k. micrum and h. triquetra cdnas. northern analysis of c. cohnii rna showed a single lsug-positive band at ~ nucleotides. this size corresponds well with the limit of conservation among lsu rrna sequences, as well as the size of the plasmodium lsug. c. cohnii lsue was also amplified and the ends determined by 5'-cdna sequencing and 3' race. northern hybridization against mitochondrial rna confirmed the presence of an ~ nucleotide rna species. the oligoadenylation sites for mitochondrial rrna fragments are identical among dinoflagellates, and either identical or within a few nucleotides of those observed in plasmodium. the 5' ends of these sequences, whether defined experimentally or by sequence conservation, are also very similar to those of their plasmodium counterparts. the only possible exception is k. micrum rna where the sole cdna obtained contained substantial upstream and downstream sequence compared to the region with similarity to plasmodium rna however it is possible that this cdna represents an unprocessed precursor, and accordingly further work is required to substantiate the size of this putative rrna fragment. secondary structure predictions for dinoflagellate sequences lsua, lsue, lsuf, lsug, rna and putative rna all indicate that the expected folding and intermolecular base pairings occur, and these fragments are likely to contribute to a viable reconstituted lsu rrna, as for plasmodium. rna editing protein-coding genes rna editing has been described for cox cob and cox transcripts from diverse dinoflagellates, including the cob mrna of k. micrum. comparison of k. micrum cdna and corresponding mtdna sequences for the three genes identified here confirms this conclusion for transcripts of cob, and further shows that cox and cox transcripts are also edited. the average density of editing of the cox transcripts is one substitution per nucleotides and this value is consistent with other studies in different species. by contrast, editing in cox transcripts is over twice as dense, at one substitution per nucleotides, making cox the most heavily edited gene transcript in dinoflagellates. editing of cob mrna lies in between these extremes, at one substitution per nucleotides. in the case of cox transcripts, four types of substitutional changes were detected at sites. of these, 48% were a to g substitutions, followed by u to c and smaller proportions of c to u and g to c edits. this observation is consistent with cox mrna editing occurring in other species, where most of the reported changes are a to g and u to c substitutions. so far, g to c changes have only been observed in mtdna-encoded mrnas of dinoflagellates, whereas a to g changes have only been reported in nucleus-encoded mrnas. cox mrna editing types are generally consistent with those observed in cox and cob mrnas. five types of substitutional changes were observed at sites, of which 42% were a to g changes, followed by c to u and u to c edits, as well as three g to a edits and a single g to c edit. for both cox and cox mrnas, the majority of substitutions occur at the first or second positions of affected codons, and over 90% of editing events result in a change in predicted amino acid. in k. micrum cox mrna, editing also removes a uag codon, which is typically a stop codon but is apparently unassigned in dinoflagellates. analysis of the cdnas corresponding to cox and cob offers further insight into the process of rna editing in dinoflagellates. despite overall uniformity of transcript editing, some cdnas exhibit pre-edited states. k. micrum cox and cob contain and editing sites, respectively, with the cdnas analyzed here representing in total and potential editing events, respectively. however at nine of these sites in the cox cdnas, and five in the cob cdnas, the pre-edited nucleotide occurs, indicating % and %'non-edits', respectively. these'non-edits' were present in only a few cdnas, suggesting that the great majority of cdnas represent mature transcripts. the pre-edited sites are scattered throughout the transcripts where they are found, occur between other edited sites, and in no obvious order in any sequence. these pre-edited sites may indicate editing failures, in which case such transcripts could give rise to defective translation products. alternatively, they may represent editing intermediates. if the latter is the case these data suggest that editing does not occur in a linear sequence along each transcript. pre-edited mitochondrial cdnas have also recently been found in a. carterae mtdna. rrna transcripts comparisons of rrna cdnas to genomic sequences are constrained by the smaller sizes of these sequences, in particular where pcr has been used to amplify genomic sequence a greater portion of this sequence represented primer binding sites and therefore cannot be used in such a comparison. nevertheless, from the available data, there is no evidence of editing of rna rna or rna for lsue, complete genomic sequence was available from the internal regions of five pcr fragments, with the majority of the sequence available from a further four pcr products using lsue primers. these sequences were identical to the cdnas except for three consecutive nucleotides that were absent in two of the three lsue cdnas obtained from the est survey. to test this anomaly, a further five cdnas were independently generated, and these all contained the three nucleotides, and therefore were identical to genomic lsue sequences and to one of the original est sequences. these results suggest that the three-nucleotide deletions seen in two cdnas represent a rare artifact, likely generated during reverse transcription, and that k. micrum lsue is likely also not edited. there was, however, evidence of substitutional editing for lsua and lsug. in both cases genomic copies of these sequences differed from transcripts: in lsug at eight positions and in lsua at six positions. consistent with the protein-coding genes, these substitutions consist mainly of a to g, c to u and u to c substitutions, with one case of c to g. given that dinoflagellate mitochondrial genes occur in multiple copies, recovery of further, independently isolated copies of these genes will be required to substantiate these inferences of rrna editing. evidence for rrna editing has also recently been reported with the dinoflagellate a. catenella, where two inferred editing events were identified for the'lsue-like' rrna. discussion prior to this study our view of the dinoflagellate mitochondrial genome was gleaned from relatively sparse molecular data obtained from several diverse dinoflagellate taxa. these data nevertheless provided a tantalizing view of a mitochondrial genome displaying several eccentricities. coding sequences for entire or partial versions of cox cob or cox have been shown to occur in multiple copies and in different genomic contexts in c. cohnii, g. polyedra, p. piscicida, and a. catenella. these data paint a picture of dinoflagellate genomes in sharp contrast to the minimalist kb apicomplexan mtdna, which encodes single copies of these genes, tightly packed together. similarly, extensive rna editing has been described in mrnas from diverse dinoflagellates, a process that does not occur in apicomplexans. in this study we have generated a much more comprehensive body of mitochondrial genomic and transcript data for two dinoflagellate species, c. cohnii and k. micrum, and these data are bolstered by a concurrent mitochondrial genomic study of the dinoflagellate a. carterae. together, these results reinforce the view that the dinoflagellate mitochondrial genome has diverged radically in form from that of apicomplexans, despite the persistence of some intriguing similarities. mitochondrial genome content and form compared to the complement of to genes in the mitochondrial genome of ciliates, the most basal member of the phylum alveolata, the very low information content of apicomplexan mtdna clearly shows that there has been considerable mitochondrial gene loss and/or relocation to the nucleus during alveolate evolution. we infer that much, if not all, of this gene relocation must have occurred prior to the last common ancestor of dinoflagellates and apicomplexans. in est surveys, we have only identified the same three protein-coding genes; moreover, we found no other mitochondrial orfs of known function in > kb and > kb of mtdna sequence from k. micrum and c. cohnii, respectively. these findings are consistent with the previous demonstration that cox an otherwise nearly ubiquitous component of mitochondrial genomes, has been relocated to the nucleus in both apicomplexans and dinoflagellates. the only additional genes we identified are ones representing the mitochondrial ssu and lsu rrnas, which together with cox and cob are universally present in mtdna. no trna genes have been found linked to mtdna sequences, similar to apicomplexans, where trnas are apparently imported into mitochondria from the cytoplasm. dinoflagellates and apicomplexans also share the characteristic of highly fragmented ssu and lsu rrnas. fragmentation of mitochondrial rrna genes has been documented in the mitochondrial genomes of several eukaryotes, including ciliates, several green algae and a fungus. the degree of fragmentation in apicomplexan mitochondrial rrna is more extreme than in these other cases, with fragments for the ssu and lsu rrnas reported to date, coding regions for which are rearranged and interspersed with other genes in the genome. from within three disparate dinoflagellate taxa we have identified eight rrna fragments similar to fragments in p. falciparum, and three of these rrna species have also recently been reported from two further taxa, a. catenella and o. marina. the dinoflagellate rrna fragments mostly appear to correspond to their plasmodium counterparts in length and sequence termini, suggesting that a stable level and pattern of fragmentation has been inherited from the common ancestor of dinoflagellates and apicomplexans. given that ciliate mitochondrial rrnas are comparatively intact, the extreme fragmentation in dinoflagellates and apicomplexans must have occurred since their divergence from ciliates. despite a similar gene content the arrangement of dinoflagellate and apicomplexan mitochondrial genomes is radically different. where the apicomplexan genome is relatively simple and compact, the dinoflagellate mitochondrial genome is complex, with multiple copies of each gene imbedded within different genomic contexts. gene fragments and non-coding regions are also repeated, altogether suggesting a great deal of recombination in the genome, which is also consistent with the lack of sequence divergence among the multiple copies of these elements. shotgun sequence data recently published for the a. carterae mitochondrion corroborate this picture of a recombining complex genome, and further suggest that the majority of the mitochondrial genome might be non-coding. gene expression in dinoflagellate mitochondria within the k. micrum est survey, long cdnas that encoded several mitochondrial genes or gene fragments were noted. by contrast, most mitochondrial cdnas we recovered encoded a single gene, suggesting the longer transcripts may be rapidly processed into shorter molecules. polycistronic transcripts up to kb are also known from apicomplexan mtdna, these are rapidly processed to short, single-gene transcripts. interestingly, the polycistronic transcripts from k. micrum are not edited, indicating that rna editing acts on the individual gene transcripts. the use of alternative initiation codons in dinoflagellate mitochondrial genes is consistent with what is seen in the mitochondria of other alveolates. in plasmodium species, cox and cox lack an in-frame atg, and while cob does contain a atg near the initiation site, it is uncertain whether initiation occurs at this site or upstream of it. att and ata have been proposed as alternative initiator codons in plasmodium species. several mitochondrial genes from the ciliate tetrahymena pyriformis also apparently use alternative initiation codons of the form atn or ntg: in the case of cob an atg within eight codons of the predicted n-terminus is apparently ignored, with gtg used in its place. thus, there are precedents for reliance on codons other than atg for translation initiation within alveolates. potential initiator atn/ntg codons exist in all three karlodinium mitochondrial genes; however, a broader survey of dinoflagellates or analysis of protein sequences will be necessary to identify the most likely candidates. an absence of stop codons is more unusual. in t. pyriformis all mitochondrial protein-coding genes terminate with taa. tga encodes tryptophan and tag is simply not used. all three plasmodium mitochondrial protein-coding genes also use taa. by contrast many dinoflagellate mitochondrial gene transcripts appear to lack any termination codon. with only a single known exception, transcripts are oligoadenylated upstream of any of the standard termination codons, and rna editing does not generate an in-frame stop. further, in none of the transcripts is a sense codon uniquely localized in the 3' region in such a way as to suggest that it serves as an alternative terminator. the oligoadenylation of k. micrum cox mrna does produces a uaa codon, as is also the case for cox transcripts for a. carterae and o. marina, that suggests that cox unlike cox or cob, might utilize conventional stop. such a mechanism for reconstituting a functional uaa is known to occur in some mammalian mitochondrial transcripts. it is unclear how the mitochondrial translation machinery might cope with the absence of termination codons. release factors that are essential for disassembly of the ribosome usually recognize specific codons, so the absence of these codons could block ribosome disassembly. there are precedents in other mitochondrial systems for the lack of termination codons: transcripts of two plant mitochondrial genes have been shown to be oligoadenylated upstream of in-frame stops. proteins encoded by both of these genes can be detected, indicating that the corresponding transcripts are successfully translated. in human mitochondria, a rare mutation has been shown to ablate a stop codon, and yet the corresponding protein is still detectable is these cell lines. eubacteria are known to be able to rescue damaged mrna molecules that have lost their termination codon by use of a specialized rna with properties of both a trna and an mrna. these so-called tmrnas restart protein synthesis by providing a terminal mrna section that encodes a functional stop codon. it has been speculated that an equivalent system might be used in plant and animal mitochondrial systems where mrnas lack stop codons. indeed, tmrna-like rna species have been identified in the mitochondria of jakobid flagellates such as r. americana; however, these rnas lack the terminal mrna-like segment of a conventional tmrna. moreover, the c-terminal tag provided by a tmrna normally targets the modified protein for degradation rather than for function. whatever the actual mechanism of translation termination in dinoflagellate mitochondria, it appears to present a clear difference with respect to protein synthesis termination in ciliate and apicomplexan mitochondria. lastly, we have found a likely case of trans-splicing of dinoflagellate mitochondrial transcripts, which adds a further layer of complexity to genome organization and expression in these organelles. while we cannot conclusively eliminate the possibility of a complete cox coding sequence in dinoflagellates we have not been able to detect an intact gene. this negative result is consistent with all other studies to date, which report only partial cox sequences from five different dinoflagellate taxa. all available data from genomic fragments and transcripts suggest that the complete cox transcript is generated by trans-splicing. such trans-splicing has not been reported for either apicomplexan or ciliate mitochondria. in ciliates nad is split into two segments but they are independently transcribed, and there is no evidence of splicing of the corresponding transcripts to create a continuous, complete nad orf. trans-splicing occurs in plant mitochondria, but in these cases the coding breakpoints are flanked by group ii intron elements, which form secondary structures that mediate the splicing events. we have no evidence of group ii introns in dinoflagellate mtdna, but we do note that the intergenic sequences contain numerous inverted repeats consistent with extensive secondary structure, which might conceivably facilitate splicing events. the unique nature of the dinoflagellate trans-splicing is also evident from the inclusion of five a residues at the splice boundary that appear to derive from the oligo tail of the upstream fragment. the removal of any downstream sequence by oligoadenylation prior to splicing argues against the involvement of a cis-acting element such as a group ii intron in the splicing process. it is conceivable that oligoadenylation of the short 5' cox transcript could serve as a degradation signal for these short transcripts, as has been observed in human mitochondria. however, lack of a complete cox coding sequence, coupled with the fact that the site of oligoadenylation corresponds with the break in coding sequence of 5' and 3' cox portions, suggests that the short cox transcripts are important intermediates in the generation of the complete cox transcripts. rna editing the rna editing observed in k. micrum cox cob and cox mrnas is consistent with the level and type of editing observed in cox and cob mrnas in other dinoflagellate species, with the exception that cox is even more heavily edited than either of cox or cob. while some editing sites are conserved, others are unique to certain taxa, suggesting that new editing sites are constantly evolving in dinoflagellates. in this study we also found evidence in k. micrum of editing of rrna fragments lsug and lsua. rna editing of a. catenella lsue has also recently been reported. at present the data are insufficient to assess the conservation of rrna editing sites among taxa; however, two inferred editing sites in a. catenella lsue are not edited in k. micrum, suggesting that rrna editing sites are constantly evolving as with those in protein-coding genes. whether rna editing plays some functional role in dinoflagellate mitochondria is unclear. from analysis of protein-coding genes in several dinoflagellates, lin et al noted that the majority of editing events are to either a c or g, thus generating a net reduction in a+u content from the bias of ~70% for the coding sequences. we observe this trend also in karlodinium protein-coding sequences. this re-tailoring of mrna sequences might better accommodate the suite of nucleus-encoded trnas that are likely imported from the cytoplasm, and which typically participate in the decoding of nucleus-encoded mrnas having a more balanced a+u content. ribosomal rna is also sensitive to a+u content, with secondary structure elements such as hairpin loops better stabilized by g-c than by a-u pairs; thus, helical regions tend to be relatively more g+c rich than other rrna domains. while the available data for rrna editing are limited, it is interesting that the editing types in rrnas have an overall neutral impact on a+u content. indeed the a+t content of mitochondrial genomic sequence specifying rrnas is already much reduced compared to that of the protein-coding genes. this observation might add weight to the notion that editing helps correct the a+t skew of protein-coding genes. the mechanism of rna editing in dinoflagellate mitochondria is also unknown; however the possibility of a guide rna -assisted mechanism, similar to that employed in trypanosomatid mitochondria, has recently been suggested. nash et al report that gene fragments encoded in mitochondria sometimes encode the'corrected' nucleotide at an inferred editing site. thus such fragments could encode templates that direct the editing events of full-length transcripts. we analyzed the k. micrum data for similar evidence of post-edited nucleotides represented in gene fragments. from five fragments that span editing sites across the three protein-encoding genes, only one site in one of the fragments corresponds to a'corrected' nucleotide seen in cdnas at an inferred editing site. an independent copy of the cob genomic sequence verified that this nucleotide difference is genuine. hence this might represent an example of an editing template in k. micrum; however, if grnas are responsible for all editing events, a very large number of additional fragments must exist to direct the remainder of the changes. clearly further work is required to shed light on the mechanism of rna editing in dinoflagellates. future directions a key question that remains is whether the observed diversity of dinoflagellate mitochondrial genes, gene fragments, and repetitive elements derives from a single mtdna molecule or from multiple chromosomes. a similar scenario of mitochondrial genes occurring as multiple copies and fragments is seen in the ichthyosporean a. parasiticum, a unicellular organism closely related to animals. in this protist, several hundred small linear chromosomes constitute the mitochondrial genome, each encoding a smattering of genes and partial genes. diplonemids, members of the phylum euglenozoa, also contain fragmented genes on separate circular mitochondrial chromosomes. it is unknown whether either of these unusual situations applies to the organization of the dinoflagellate mitochondrial genome; however, in this regard we make two preliminary observations. one is that long-range pcr was unable to generate longer contiguous sequences linking the many mtdna elements we report in this study. rather, additional short unique gene linkages were obtained, and it is clear that we have yet to sample the full diversity of gene combinations. secondly, the presence of individual genes in partial tandem repeats is consistent with minicircles, as seen in dinoflagellate plastid genomes. if these cases represent true minicircles, we have been unable to amplify a corresponding sequence to close these circles. it is also possible, of course, that the tandem repeats that we observe are simply a consequence of further recombination events, and the high diversity of gene combinations.  a greater depth of sampling of dinoflagellate mitochondrial dna and mrna has provided a clearer view of a complex genome and many peculiarities of gene expression. we find that the dinoflagellate mitochondrial genome shares several features in common with the mtdna of its apicomplexan sister lineage, but also many novel characteristics. features in common for the two lineages are: a very high level of gene relocation from the mitochondrion, extensive rrna gene fragmentation and dispersal, and use of non-standard initiation codons. features unique to dinoflagellates are: gene copy number expansion and reorganization, loss of stop codons from protein-coding genes, mrna trans-splicing, and rna editing of protein-coding and rrna transcripts. these data demonstrate a remarkable burst of organelle genome evolution in dinoflagellates following divergence from apicomplexa, and also challenge our understanding of the mechanistic details of genome maintenance and expression, most notably translation termination. 
the mitochondrial genome of apicomplexans is uniquely reduced in gene content and size, encoding only three proteins and two ribosomal rnas within a highly compacted kb dna. the dinoflagellate mitochondrial genome shares the same gene complement and fragmentation of rrna genes with its apicomplexan counterpart. dinoflagellate mitochondrial genomes have been comparatively poorly studied: limited available data suggest some similarities with apicomplexan mitochondrial genomes but an even more radical type of genomic organization. here, we investigate structure, content and expression of dinoflagellate mitochondrial genomes. from two dinoflagellates, crypthecodinium cohnii and karlodinium micrum, we generated over kb of mitochondrial genomic data that indicate a reduced gene content paralleling that of mitochondrial genomes in apicomplexans, i.e., only three protein-encoding genes and at least eight conserved components of the highly fragmented large and small subunit rrnas. polycistronic transcripts were found, standard start codons are absent, and oligoadenylation occurs upstream of stop codons, resulting in the absence of termination codons. unlike in apicomplexans, dinoflagellate mitochondrial genes occur in multiple copies, often as gene fragments, and in numerous genomic contexts. analysis of cdnas suggests several novel aspects of dinoflagellate mitochondrial gene expression. rna substitutional editing, a process previously identified for mrnas in dinoflagellate mitochondria, is also implicated in rrna expression.
 dinoflagellates comprise an ecologically significant and diverse eukaryotic phylum that is sister to the phylum containing apicomplexan endoparasites. the mitochondrial genome of apicomplexans is uniquely reduced in gene content and size, encoding only three proteins and two ribosomal rnas within a highly compacted kb dna. dinoflagellate mitochondrial genomes have been comparatively poorly studied: limited available data suggest some similarities with apicomplexan mitochondrial genomes but an even more radical type of genomic organization. here, we investigate structure, content and expression of dinoflagellate mitochondrial genomes.  from two dinoflagellates, crypthecodinium cohnii and karlodinium micrum, we generated over kb of mitochondrial genomic data that indicate a reduced gene content paralleling that of mitochondrial genomes in apicomplexans, i.e., only three protein-encoding genes and at least eight conserved components of the highly fragmented large and small subunit rrnas. unlike in apicomplexans, dinoflagellate mitochondrial genes occur in multiple copies, often as gene fragments, and in numerous genomic contexts. analysis of cdnas suggests several novel aspects of dinoflagellate mitochondrial gene expression. polycistronic transcripts were found, standard start codons are absent, and oligoadenylation occurs upstream of stop codons, resulting in the absence of termination codons. transcripts of at least one gene, cox are apparently trans-spliced to generate full-length mrnas. rna substitutional editing, a process previously identified for mrnas in dinoflagellate mitochondria, is also implicated in rrna expression.  the dinoflagellate mitochondrial genome shares the same gene complement and fragmentation of rrna genes with its apicomplexan counterpart. however, it also exhibits several unique characteristics. most notable are the expansion of gene copy numbers and their arrangements within the genome, rna editing, loss of stop codons, and use of trans-splicing. 