The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): introns, a paucity of tRNA genes, and a near-standard genetic code

The complete mitochondrial genome sequence of the hydrothermal vent galatheid crab Shinkaia crosnieri (Crustacea: Decapoda: Anomura): a novel arrangement and incomplete tRNA suite

Background

Metazoan mitochondrial genomes usually consist of the same 37 genes. Such genes contain useful information for phylogenetic analyses and evolution modelling. Although complete mitochondrial genomes have been determined for over 1,000 animals to date, hydrothermal vent species have, thus far, remained excluded due to the scarcity of collected specimens.

Results

The mitochondrial genome of the hydrothermal vent galatheid crab Shinkaia crosnieri is 15,182 bp in length, and is composed of 13 protein-coding genes, two ribosomal RNA genes and only 18 transfer RNA genes. The total AT content of the genome, as is typical for decapods, is 72.9%. We identified a non-coding control region of 327 bp according to its location and AT-richness. This is the smallest control region discovered in crustaceans so far. A mechanism of cytoplasmic tRNA import was addressed to compensate for the four missing tRNAs. The S. crosnieri mitogenome exhibits a novel arrangement of mitochondrial genes. We investigated the mitochondrial gene orders and found that at least six rearrangements from the ancestral pancrustacean (crustacean + hexapod) pattern have happened successively. The codon usage, nucleotide composition and bias show no substantial difference with other decapods. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of the 13 protein-coding genes prove consistent with the previous classification based upon their morphology.

Conclusion

The present study will supply considerable data of use for both genomic and evolutionary research on hydrothermal vent ecosystems. The mitochondrial genetic characteristics of decapods are sustained in this case of S. crosnieri despite the absence of several tRNAs and a number of dramatic rearrangements. Our results may provide evidence for the immigrating hypothesis about how vent species originate.

Group II introns break new boundaries: presence in a bilaterian's genome

Group II introns are ribozymes, removing themselves from their primary transcripts, as well as mobile genetic elements, transposing via an RNA intermediate, and are thought to be the ancestors of spliceosomal introns. Although common in bacteria and most eukaryotic organelles, they have never been reported in any bilaterian animal genome, organellar or nuclear. Here we report the first group II intron found in the mitochondrial genome of a bilaterian worm. This location is especially surprising, since animal mitochondrial genomes are generally distinct from those of plants, fungi, and protists by being small and compact, and so are viewed as being highly streamlined, perhaps as a result of strong selective pressures for fast replication while establishing germ plasm during early development. This intron is found in the mtDNA of an annelid worm, (an undescribed species of Nephtys), where the complete sequence revealed a 1819 bp group II intron inside the cox1 gene. We infer that this intron is the result of a recent horizontal gene transfer event from a viral or bacterial vector into the mitochondrial genome of Nephtys sp. Our findings hold implications for understanding mechanisms, constraints, and selective pressures that account for patterns of animal mitochondrial genome evolution

The mitochondrial genome of Hydra oligactis (Cnidaria, Hydrozoa) sheds new light on animal mtDNA evolution and cnidarian phylogeny

The 16,314-nuceotide sequence of the linear mitochondrial DNA (mtDNA) molecule of Hydra oligactis (Cnidaria, Hydrozoa)--the first from the class Hydrozoa--has been determined. This sequence contains genes for 13 energy pathway proteins, small and large subunit rRNAs, and methionine and tryptophan tRNAs, as is typical for cnidarians. All genes have the same transcriptional orientation and their arrangement in the genome is similar to that of the jellyfish Aurelia aurita. In addition, a partial copy of cox1 is present at one end of the molecule in a transcriptional orientation opposite to the rest of the genes, forming a part of inverted terminal repeat characteristic of linear mtDNA and linear mitochondrial plasmids. The sequence close to at least one end of the molecule contains several homonucleotide runs as well as small inverted repeats that are able to form strong secondary structures and may be involved in mtDNA maintenance and expression. Phylogenetic analysis of mitochondrial genes of H. oligactis and other cnidarians supports the Medusozoa hypothesis but also suggests that Anthozoa may be paraphyletic, with octocorallians more closely related to the Medusozoa than to the Hexacorallia. The latter inference implies that Anthozoa is paraphyletic and that the polyp (rather than a medusa) is the ancestral body type in Cnidaria.

The complete mitochondrial genomes of needle corals, Seriatopora spp. (Scleractinia: Pocilloporidae): an idiosyncratic atp8, duplicated trnW gene, and hypervariable regions used to determine species phylogenies and recently diverged populations

Complete DNA sequences were determined for the mitochondrial (mt) genomes of the needle corals, Seriatopora caliendrum (17,011bp) and S. hystrix (17,060bp). Gene arrangement of the Seriatopora mt genomes is similar to the 14 currently published scleractinian mitogenomes with three unusual features, including an idiosyncratic atp8, a duplicated trnW (tRNA(TRP)), and a putative control region located between atp6 and nad4. Atp8, located between duplicate trnW genes, showed relatively low amino acid similarity (25.6-34.6%) with those of published scleractinian corals. A reverse-transcription polymerase chain reaction confirmed the transcription of this novel atp8 gene in Seriatopora. A duplicated trnW was detected in the region close to the cox1 gene and shares the highly conserved primary and secondary structure of its original counterpart. The intergenic spacer between atp6 and nad4, which contains several distinct repeated elements, is being designated as the putative control region in the Seriatopora mt genomes. Evaluation of the molecular evolution of several protein-coding genes and intergenic spacers showed 3- to 4-fold higher divergence rates among populations or between species than those published for scleractinian mt genomes. This study not only successfully revealed the phylogenies of S. hystrix and S. caliendrum from the West Pacific Ocean by mtDNA, but also highlighted the potential utilities of mt hypervariable regions in phylogenetic construction below the species level for Seriatopora.

The mitochondrial genome of Pocillopora (Cnidaria: Scleractinia) contains two variable regions: the putative D-loop and a novel ORF of unknown function

The complete mitochondrial genomes of two individuals attributed to different morphospecies of the scleractinian coral genus Pocillopora have been sequenced. Both genomes, respectively 17,415 and 17,422 nt long, share the presence of a previously undescribed ORF encoding a putative protein made up of 302 amino acids and of unknown function. Surprisingly, this ORF turns out to be the second most variable region of the mitochondrial genome (1% nucleotide sequence difference between the two individuals) after the putative control region (1.5% sequence difference). Except for the presence of this ORF and for the location of the putative control region, the mitochondrial genome of Pocillopora is organized in a fashion similar to the other scleractinian coral genomes published to date. For the first time in a cnidarian, a putative second origin of replication is described based on its secondary structure similar to the stem-loop structure of O(L), the origin of L-strand replication in vertebrates.

Novel group I introns encoding a putative homing endonuclease in the mitochondrial cox1 gene of Scleractinian corals

Analyses of mitochondrial sequences revealed the existence of a group I intron in the cytochrome oxidase subunit 1 (cox1) gene in 13 of 41 genera (20 out of 73 species) of corals conventionally assigned to the suborder Faviina. With one exception, phylogenies of the coral cox1 gene and its intron were concordant, suggesting at most two insertions and many subsequent losses. The coral introns were inferred to encode a putative homing endonuclease with a LAGLI-DADG motif as reported for the cox1 group I intron in the sea anemone Metridium senile. However, the coral and sea anemone cox1 group I introns differed in several aspects, such as the intron insertion site and sequence length. The coral cox1 introns most closely resemble the mitochondrial cox1 group I introns of a sponge species, which also has the same insertion site. The coral introns are also more similar to the introns of several fungal species than to that of the sea anemone (although the insertion site differs in the fungi). This suggests either a horizontal transfer between a sponge and a coral or independent transfers from a similar fungal donor (perhaps one with an identical insertion site that has not yet been discovered). The common occurrence of this intron in corals strengthens the evidence for an elevated abundance of group I introns in the mitochondria of anthozoans.

The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria:Anthozoa:Antipatharia) supports classification of antipatharians within the subclass Hexacorallia

Black corals comprise a globally distributed shallow- and deep-water taxon whose phylogenetic position within the Anthozoa has been debated. We sequenced the complete mitochondrial genome of the antipatharian Chrysopathes formosa to further evaluate its phylogenetic relationships. The circular mitochondrial genome (18,398 bp) consists of 13 energy pathway protein-coding genes and two ribosomal RNAs, but only two transfer RNA genes (trnM and trnW), as well as a group I intron within the nad5 gene that contains the only copies of nad1 and nad3. No novel genes were found in the antipatharian mitochondrial genome. Gene order and genome content are most similar to those of the sea anemone Metridium senile (subclass Hexacorallia), with differences being the relative location of three contiguous genes (cox2-nad4-nad6) and absence (from the antipatharian) of a group I intron within the cox1 gene. Phylogenetic analyses of multiple protein-coding genes support classifying the Antipatharia within the subclass Hexacorallia and not the subclass Ceriantipatharia; however, the sister-taxon relationships of black corals within Hexacorallia remain inconclusive.

Mitochondrial Genome of Savalia savaglia (Cnidaria, Hexacorallia) and Early Metazoan Phylogeny

Mitochondrial genomes have recently become widely used in animal phylogeny, mainly to infer the relationships between vertebrates and other bilaterians. However, only 11 of 723 complete mitochondrial genomes available in the public databases are of early metazoans, including cnidarians (Anthozoa, mainly Scleractinia) and sponges. Although some cnidarians (Medusozoa) are known to possess atypical linear mitochondrial DNA, the anthozoan mitochondrial genome is circular and its organization is similar to that of other metazoans. Because the phylogenetic relationships among Anthozoa as well as their relation to other early metazoans still need to be clarified, we tested whether sequencing the complete mitochondrial genome of Savalia savaglia, an anthozoan belonging to the order Zoantharia (=Zoanthidea), could be useful to infer such relationships. Compared to other anthozoans, S. savaglia's genome is unusually long (20,766 bp) due to the presence of several noncoding intergenic regions (3691 bp). The genome contains all 13 protein coding genes commonly found in metazoans, but like other Anthozoa it lacks most of the tRNAs. Phylogenetic analyses of S. savaglia mitochondrial sequences show Zoantharia branching closely to other Hexacorallia, either as a sister group to Actiniaria or as a sister group to Actiniaria and Scleractinia. The close relationships suggested between Zoantharia and Actiniaria are reinforced by strong similarities in their gene order and the presence of similar introns in the COI and ND5 genes. Our study suggests that mitochondrial genomes can be a source of potentially valuable information on the phylogeny of Hexacorallia and may provide new insights into the evolution of early metazoans.

Morphological and molecular characterisation of Abyssoanthus nankaiensis,, a new family, new genus and new species of deep-sea zoanthid (Anthozoa:Hexacorallia:Zoantharia) from a north-west Pacific methane cold seep

The order Zoantharia (= Zoanthiniaria, Zoanthidea) remains one of the most taxonomically neglected anthozoan groups, primarily owing to the difficulty of examining the internal morphology of sand-encrusted zoanthids, upon which classification has been largely based until very recently. Additionally, relatively few zoanthid species (all belonging to the genus Epizoanthus) have been properly described from the deep-sea. Numerous individuals of an unidentified sediment-encrusting zoanthid-like species were observed and sampled during a Shinkai 6500 deep-sea submersible dive at a methane cold seep (depth 3259 m) off Muroto at the Nankai Trough, Japan. Unlike previously described deep-sea zoanthids, Abyssoanthus nankaiensis, gen. et sp. nov. (Abyssoanthidae fam. nov.) was found to be non-colonial, free-living (non-commensal), and, uniquely, on mudstone in the vicinity of a methane cold seep. Morphologically, A. nankaiensis, gen. et sp. nov. is characterised by its relatively uniform polyp diameter from oral to aboral end with 19–22 mesenteries. Additionally, DNA (mitochondrial 16S rDNA and mitochondrial cytochrome oxidase c subunit I DNA) sequences obtained from these samples also unambiguously place this specimen in a previously undescribed and new family within the order Zoantharia. This is the first reported zoanthid species from a methane cold seep or other so-called ‘extreme’ environment, and the first molecular characterisation of any such deep-sea zoanthid.

Mitochondrial genome of the moon jelly Aurelia aurita (Cnidaria, Scyphozoa): A linear DNA molecule encoding a putative DNA-dependent DNA polymerase

The 16,937-nuceotide sequence of the linear mitochondrial DNA (mt-DNA) molecule of the moon jelly Aurelia aurita (Cnidaria, Scyphozoa) - the first mtDNA sequence from the class Scypozoa and the first sequence of a linear mtDNA from Metazoa - has been determined. This sequence contains genes for 13 energy pathway proteins, small and large subunit rRNAs, and methionine and tryptophan tRNAs. In addition, two open reading frames of 324 and 969 base pairs in length have been found. The deduced amino-acid sequence of one of them, ORF969, displays extensive sequence similarity with the polymerase [but not the exonuclease] domain of family B DNA polymerases, and this ORF has been tentatively identified as dnab. This is the first report of dnab in animal mtDNA. The genes in A. aurita mtDNA are arranged in two clusters with opposite transcriptional polarities; transcription proceeding toward the ends of the molecule. The determined sequences at the ends of the molecule are nearly identical but inverted and lack any obvious potential secondary structures or telomere-like repeat elements. The acquisition of mitochondrial genomic data for the second class of Cnidaria allows us to reconstruct characteristic features of mitochondrial evolution in this animal phylum.

Putative cross-kingdom horizontal gene transfer in sponge (Porifera) mitochondria

BACKGROUND

The mitochondrial genome of Metazoa is usually a compact molecule without introns. Exceptions to this rule have been reported only in corals and sea anemones (Cnidaria), in which group I introns have been discovered in the cox1 and nad5 genes. Here we show several lines of evidence demonstrating that introns can also be found in the mitochondria of sponges (Porifera).

RESULTS

A 2,349 bp fragment of the mitochondrial cox1 gene was sequenced from the sponge Tetilla sp. (Spirophorida). This fragment suggests the presence of a 1143 bp intron. Similar to all the cnidarian mitochondrial introns, the putative intron has group I intron characteristics. The intron is present in the cox1 gene and encodes a putative homing endonuclease. In order to establish the distribution of this intron in sponges, the cox1 gene was sequenced from several representatives of the demosponge diversity. The intron was found only in the sponge order Spirophorida. A phylogenetic analysis of the COI protein sequence and of the intron open reading frame suggests that the intron may have been transmitted horizontally from a fungus donor.

CONCLUSION

Little is known about sponge-associated fungi, although in the last few years the latter have been frequently isolated from sponges. We suggest that the horizontal gene transfer of a mitochondrial intron was facilitated by a symbiotic relationship between fungus and sponge. Ecological relationships are known to have implications at the genomic level. Here, an ecological relationship between sponge and fungus is suggested based on the genomic analysis.

GenDecoder: genetic code prediction for metazoan mitochondria

Although the majority of the organisms use the same genetic code to translate DNA, several variants have been described in a wide range of organisms, both in nuclear and organellar systems, many of them corresponding to metazoan mitochondria. These variants are usually found by comparative sequence analyses, either conducted manually or with the computer. Basically, when a particular codon in a query-species is linked to positions for which a specific amino acid is consistently found in other species, then that particular codon is expected to translate as that specific amino acid. Importantly, and despite the simplicity of this approach, there are no available tools to help predicting the genetic code of an organism. We present here GenDecoder, a web server for the characterization and prediction of mitochondrial genetic codes in animals. The analysis of automatic predictions for 681 metazoans aimed us to study some properties of the comparative method, in particular, the relationship among sequence conservation, taxonomic sampling and reliability of assignments. Overall, the method is highly precise (99%), although highly divergent organisms such as platyhelminths are more problematic. The GenDecoder web server is freely available from http://darwin.uvigo.es/software/gendecoder.html.

The complete sequence of the mitochondrial genome of Nautilus macromphalus (Mollusca: Cephalopoda)

BACKGROUND

Mitochondria contain small genomes that are physically separate from those of nuclei. Their comparison serves as a model system for understanding the processes of genome evolution. Although complete mitochondrial genome sequences have been reported for more than 600 animals, the taxonomic sampling is highly biased toward vertebrates and arthropods, leaving much of the diversity yet uncharacterized.

RESULTS

The mitochondrial genome of the bellybutton nautilus, Nautilus macromphalus, a cephalopod mollusk, is 16,258 nts in length and 59.5% A+T, both values that are typical of animal mitochondrial genomes. It contains the 37 genes that are almost universally found in animal mtDNAs, with 15 on one DNA strand and 22 on the other. The arrangement of these genes can be derived from that of the distantly related Katharina tunicata (Mollusca: Polyplacophora) by a switch in position of two large blocks of genes and transpositions of four tRNA genes. There is strong skew in the distribution of nucleotides between the two strands, and analysis of this yields insight into modes of transcription and replication. There is an unusual number of non-coding regions and their function, if any, is not known; however, several of these demark abrupt shifts in nucleotide skew, and there are several identical sequence elements at these junctions, suggesting that they may play roles in transcription and/or replication. One of the non-coding regions contains multiple repeats of a tRNA-like sequence. Some of the tRNA genes appear to overlap on the same strand, but this could be resolved if the polycistron were cleaved at the beginning of the downstream gene, followed by polyadenylation of the product of the upstream gene to form a fully paired structure.

CONCLUSION

Nautilus macromphalus mtDNA contains an expected gene content that has experienced few rearrangements since the evolutionary split between cephalopods and polyplacophorans. It contains an unusual number of non-coding regions, especially considering that these otherwise often are generated by the same processes that produce gene rearrangements. The skew in nucleotide composition between the two strands is strong and associated with the direction of transcription in various parts of the genomes, but a comparison with K. tunicata implies that mutational bias during replication also plays a role. This appears to be yet another case where polyadenylation of mitochondrial tRNAs restores what would otherwise be an incomplete structure.

Mitochondrial genome of Trichoplax adhaerens supports placozoa as the basal lower metazoan phylum

Mitochondrial genomes of multicellular animals are typically 15- to 24-kb circular molecules that encode a nearly identical set of 12-14 proteins for oxidative phosphorylation and 24-25 structural RNAs (16S rRNA, 12S rRNA, and tRNAs). These genomes lack significant intragenic spacers and are generally without introns. Here, we report the complete mitochondrial genome sequence of the placozoan Trichoplax adhaerens, a metazoan with the simplest known body plan of any animal, possessing no organs, no basal membrane, and only four different somatic cell types. Our analysis shows that the Trichoplax mitochondrion contains the largest known metazoan mtDNA genome at 43,079 bp, more than twice the size of the typical metazoan mtDNA. The mitochondrion's size is due to numerous intragenic spacers, several introns and ORFs of unknown function, and protein-coding regions that are generally larger than those found in other animals. Not only does the Trichoplax mtDNA have characteristics of the mitochondrial genomes of known metazoan outgroups, such as chytrid fungi and choanoflagellates, but, more importantly, it shares derived features unique to the Metazoa. Phylogenetic analyses of mitochondrial proteins provide strong support for the placement of the phylum Placozoa at the root of the Metazoa.

Naked corals: skeleton loss in Scleractinia

Stony corals, which form the framework for modern reefs, are classified as Scleractinia (Cnidaria, Anthozoa, and Hexacorallia) in reference to their external aragonitic skeletons. However, persistent notions, collectively known as the "naked coral" hypothesis, hold that the scleractinian skeleton does not define a natural group. Three main lines of evidence have suggested that some stony corals are more closely related to one or more of the soft-bodied hexacorallian groups than they are to other scleractinians: (i) morphological similarities; (ii) lack of phylogenetic resolution in molecular analyses of scleractinians; and (iii) discrepancy between the commencement of a diverse scleractinian fossil record at 240 million years ago (Ma) and a molecule-based origination of at least 300 Ma. No molecular evidence has been able to clearly reveal relationships at the base of a well supported clade composed of scleractinian lineages and the nonskeletonized Corallimorpharia. We present complete mitochondrial genome data that provide strong evidence that one clade of scleractinians is more closely related to Corallimorpharia than it is to a another clade of scleractinians. Thus, the scleractinian skeleton, which we estimate to have originated between 240 and 288 Ma, was likely lost in the ancestry of Corallimorpharia. We estimate that Corallimorpharia originated between 110 and 132 Ma during the late- to mid-Cretaceous, coinciding with high levels of oceanic CO(2), which would have impacted aragonite solubility. Corallimorpharians escaped extinction from aragonite skeletal dissolution, but some modern stony corals may not have such fortunate fates under the pressure of increased anthropogenic CO(2) in the ocean.

Complete mitochondrial DNA sequence of a Conoidean gastropod, Lophiotoma (Xenuroturris) cerithiformis: gene order and gastropod phylogeny

We have determined the first complete nucleotide sequence of the mitochondrial genome of a venomous mollusc, the Conoidean gastropod, Lophiotoma (Xenuroturris) cerithiformis. It is 15,380 nucleotide pairs (ntp) and encodes 13 proteins, two ribosomal RNAs and 22 tRNAs of the mitochondrion's own protein synthesizing system. The protein mRNAs, ribosomal RNAs and 13 of the tRNAs are transcribed from the same strand, the remaining tRNAs from the other strand. The longest segment of unassigned sequence is 139 ntp and includes a 82 ntp segment that is a perfect inverted repeat sequence of 37 ntp separated by 8 nt. The gene arrangement of L. cerithiformis mtDNA shows remarkable similarity to the gene arrangements of mtDNAs of the vetigastropod Haliotis rubra, the polyplacophoran Katharina tunicata and the cephalopod Octopus vulgaris, but differs dramatically from the gene arrangements found in the mtDNAs of pulmonate and opisthobranch gastropods, as well as mtDNAs of bivalves and scaphopods. A single sixteen gene inversion that distinguishes L. cerithiformis mtDNA from mtDNAs of H. rubra, K. tunicata and O. vulgaris is shared by mtDNA of a littorinomorph gastropod Littorina saxitalis, suggesting a close relationship of conoidean and littorinomorph gastropods.

No variation and low synonymous substitution rates in coral mtDNA despite high nuclear variation

Background

The mitochondrial DNA (mtDNA) of most animals evolves more rapidly than nuclear DNA, and often shows higher levels of intraspecific polymorphism and population subdivision. The mtDNA of anthozoans (corals, sea fans, and their kin), by contrast, appears to evolve slowly. Slow mtDNA evolution has been reported for several anthozoans, however this slow pace has been difficult to put in phylogenetic context without parallel surveys of nuclear variation or calibrated rates of synonymous substitution that could permit quantitative rate comparisons across taxa. Here, I survey variation in the coding region of a mitochondrial gene from a coral species (Balanophyllia elegans) known to possess high levels of nuclear gene variation, and estimate synonymous rates of mtDNA substitution by comparison to another coral (Tubastrea coccinea).

Results

The mtDNA surveyed (630 bp of cytochrome oxidase subunit I) was invariant among individuals sampled from 18 populations spanning 3000 km of the range of B. elegans, despite high levels of variation and population subdivision for allozymes over these same populations. The synonymous substitution rate between B. elegans and T. coccinea (0.05%/site/10⁶ years) is similar to that in most plants, but 50–100 times lower than rates typical for most animals. In addition, while substitutions to mtDNA in most animals exhibit a strong bias toward transitions, mtDNA from these corals does not.

Conclusion

Slow rates of mitochondrial nucleotide substitution result in low levels of intraspecific mtDNA variation in corals, even when nuclear loci vary. Slow mtDNA evolution appears to be the basal condition among eukaryotes. mtDNA substitution rates switch from slow to fast abruptly and unidirectionally. This switch may stem from the loss of just one or a few mitochondrion-specific DNA repair or replication genes.

Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements

Phylogenetic relationships among the metazoan phyla are the subject of an ongoing controversy. Analysis of mitochondrial gene arrangements is a powerful tool to investigate these relationships; however, its previous application outside of individual animal phyla has been hampered by the lack of informative out-group data. To address this shortcoming, we determined complete mitochondrial DNA sequences for the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. With sponges as an outgroup, we investigated phylogenetic relationships of nine bilaterian phyla using both breakpoint analysis of global mitochondrial gene arrangements and maximum parsimony analysis of mitochondrial gene adjacencies. Our results provide strong support for a group that includes protostome (but not deuterostome) coelomate, pseudocoelomate, and acoelomate animals, thus clearly rejecting the Coelomata hypothesis. Two other groups of bilaterian animals, Lophotrochozoa and Ambulacraria, are also supported by our analyses. However, due to the remarkable stability of mitochondrial gene arrangements in Deuterostomia and the Ecdysozoa, conclusions on their evolutionary history cannot be drawn.

The Mitochondrial Genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): Considerable Economization in the Length and Structural Features of Encoded Genes

The complete sequence of the mitochondrial genome of the plant parasitic nematode Xiphinema americanum sensu stricto has been determined. At 12626bp it is the smallest metazoan mitochondrial genome reported to date. Genes are transcribed from both strands. Genes coding for 12 proteins, 2 rRNAs and 17 putative tRNAs (with the tRNA-C, I, N, S1, S2 missing) are predicted from the sequence. The arrangement of genes within the X. americanum mitochondrial genome is unique and includes gene overlaps. Comparisons with the mtDNA of other nematodes show that the small size of the X. americanum mtDNA is due to a combination of factors. The two mitochondrial rRNA genes are considerably smaller than those of other nematodes, with most of the protein encoding and tRNA genes also slightly smaller. In addition, five tRNAs genes are absent, lengthy noncoding regions are not present in the mtDNA, and several gene overlaps are present.

Novel Repetitive Structures, Deviant Protein-Encoding Sequences andUnidentified ORFs in the Mitochondrial Genome of the BrachiopodLingula anatina

Complete sequence determination of the brachiopod Lingula anatina mtDNA (28,818 bp) revealed an organization that is remarkably atypical for an animal mt-genome. In addition to the usual set of 37 animal mitochondrial genes, which make up only 57% (16,555 bp) of the entire sequence, the genome contains lengthy unassigned sequences. All the genes are encoded in the same DNA strand, generally in a compact way, whereas the overall gene order is highly divergent in comparison with known animal mtDNA. Individual genes are generally longer and deviate considerably in sequence from their homologues in other animals. The genome contains two major repeat regions, in which 11 units of unassigned sequences and six genes (atp8, trnM, trnQ, trnV, and part of cox2 and nad2) are found in repetition, in the form of nested direct repeats of unparalleled complexity. One of the repeat regions contains unassigned repeat units dispersed among several unique sequences, novel repetitive structure for animal mtDNAs. Each of those unique sequences contains an open reading frame for a polypeptide between 80 and 357 amino acids long, potentially encoding a functional molecule, but none of them has been identified with known proteins. In both repeat regions, tRNA genes or tRNA gene-like sequences flank major repeated units, supporting the view that those structures play a role in the mitochondrial gene rearrangements. Although the intricate repeated organization of this genome can be explained by recurrent tandem duplications and subsequent deletions mediated by replication errors, other mechanisms, such as nonhomologous recombinations, appear to explain certain structures more easily.

Sequencing and comparing whole mitochondrial genomes of animals

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, which can be especially powerful. We describe here the protocols commonly used for physically isolating mitochondrial DNA (mtDNA), for amplifying these by polymerase chain reaction (PCR) or rolling circle amplification (RCA), for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences with determining and comparing complete mitochondrial DNA sequences.

Mitochondrial Genetic Codes Evolve to Match Amino Acid Requirements of Proteins

Mitochondria often use genetic codes different from the standard genetic code. Now that many mitochondrial genomes have been sequenced, these variant codes provide the first opportunity to examine empirically the processes that produce new genetic codes. The key question is: Are codon reassignments the sole result of mutation and genetic drift? Or are they the result of natural selection? Here we present an analysis of 24 phylogenetically independent codon reassignments in mitochondria. Although the mutation-drift hypothesis can explain reassignments from stop to an amino acid, we found that it cannot explain reassignments from one amino acid to another. In particular--and contrary to the predictions of the mutation-drift hypothesis--the codon involved in such a reassignment was not rare in the ancestral genome. Instead, such reassignments appear to take place while the codon is in use at an appreciable frequency. Moreover, the comparison of inferred amino acid usage in the ancestral genome with the neutral expectation shows that the amino acid gaining the codon was selectively favored over the amino acid losing the codon. These results are consistent with a simple model of weak selection on the amino acid composition of proteins in which codon reassignments are selected because they compensate for multiple slightly deleterious mutations throughout the mitochondrial genome. We propose that the selection pressure is for reduced protein synthesis cost: most reassignments give amino acids that are less expensive to synthesize. Taken together, our results strongly suggest that mitochondrial genetic codes evolve to match the amino acid requirements of proteins.

Polymorphism in Nucleotide Sequence of Mitochondrial Intergenic Region in Scleractinian Coral (Galaxea fascicularis)

A region of 826 bp that is unlikely to code for a protein, ribosomal RNA, or transfer RNA was identified between the cytochrome b and NADH-ubiquinone oxidoreductase chain 2 loci in the mitochondrial DNA of the scleractinian reef coral Galaxea fascicularis. Nucleotide sequences were determined in a part (625 bp) of this intergenic region in 95 individuals collected at 9 sites in the Ryukyu Archipelago in southwestern Japan. A total of 8 haplotypes were found, and a deletion of 290 bp was found in 3 of them. Significant differences were found in frequencies of the haplotypes at 3 sampling sites. The presence or absence of the deletion was highly correlated with the hard or soft morphotype. The deletion was found in the majority of hard-type colonies, but in a small fraction of soft-type individuals.

Forces maintaining organellar genomes: is any as strong as genetic code disparity or hydrophobicity?

It remains controversial why mitochondria and chloroplasts retain the genes encoding a small subset of their constituent proteins, despite the transfer of so many other genes to the nucleus. Two candidate obstacles to gene transfer, suggested long ago, are that the genetic code of some mitochondrial genomes differs from the standard nuclear code, such that a transferred gene would encode an incorrect amino acid sequence, and that the proteins most frequently encoded in mitochondria are generally very hydrophobic, which may impede their import after synthesis in the cytosol. More recently it has been suggested that both these interpretations suffer from serious "false positives" and "false negatives": genes that they predict should be readily transferred but which have never (or seldom) been, and genes whose transfer has occurred often or early, even though this is predicted to be very difficult. Here I consider the full known range of ostensibly problematic such genes, with particular reference to the sequences of events that could have led to their present location. I show that this detailed analysis of these cases reveals that they are in fact wholly consistent with the hypothesis that code disparity and hydrophobicity are much more powerful barriers to functional gene transfer than any other. The popularity of the contrary view has led to the search for other barriers that might retain genes in organelles even more powerfully than code disparity or hydrophobicity; one proposal, concerning the role of proteins in redox processes, has received widespread support. I conclude that this abandonment of the original explanations for the retention of organellar genomes has been premature. Several other, relatively minor, obstacles to gene transfer certainly exist, contributing to the retention of relatively many organellar genes in most lineages compared to animal mtDNA, but there is no evidence for obstacles as severe as code disparity or hydrophobicity. One corollary of this conclusion is that there is currently no reason to suppose that engineering nuclear versions of the remaining mammalian mitochondrial genes, a feat that may have widespread biomedical relevance, should require anything other than sequence alterations obviating code disparity and causing modest reductions in hydrophobicity without loss of enzymatic function.

Mitochondria of Protists

Over the past several decades, our knowledge of the origin and evolution of mitochondria has been greatly advanced by determination of complete mitochondrial genome sequences. Among the most informative mitochondrial genomes have been those of protists (primarily unicellular eukaryotes), some of which harbor the most gene-rich and most eubacteria-like mitochondrial DNAs (mtDNAs) known. Comparison of mtDNA sequence data has provided insights into the radically diverse trends in mitochondrial genome evolution exhibited by different phylogenetically coherent groupings of eukaryotes, and has allowed us to pinpoint specific protist relatives of the multicellular eukaryotic lineages (animals, plants, and fungi). This comparative genomics approach has also revealed unique and fascinating aspects of mitochondrial gene expression, highlighting the mitochondrion as an evolutionary playground par excellence.

Variation in coding (NADH dehydrogenase subunits 2, 3, and 6) and noncoding intergenic spacer regions of the mitochondrial genome in Octocorallia (Cnidaria: Anthozoa)

Low rates of evolution in cnidarian mitochondrial genes such as COI and 16S rDNA have hindered molecular systematic studies in this important invertebrate group. We sequenced fragments of 3 mitochondrial protein-coding genes (NADH dehydrogenase subunits ND2, ND3 and ND6) as well as the COI-COII intergenic spacer, the longest noncoding region found in the octocoral mitochondrial genome, to determine if any of these regions contain levels of variation sufficient for reconstruction of phylogenetic relationships among genera of the anthozoan subclass Octocorallia. Within and between the soft coral families Alcyoniidae and Xeniidae, sequence divergence in the genes ND2 (539 bp), ND3 (102 bp), and ND6 (444 bp) ranged from 0.5% to 12%, with the greatest pairwise distances between the 2 families. The COI-COII intergenic spacer varied in length from 106 to 122 bp, and pairwise sequence divergence values ranged from 0% to 20.4%. Phylogenetic trees constructed using each region separately were poorly resolved. Better phylogenetic resolution was obtained in a combined analysis using all 3 protein-coding regions (1085 bp total). Although relationships among some pairs of species and genera were well supported in the combined analysis, the base of the alcyoniid family tree remained an unresolved polytomy. We conclude that variation in the NADH subunit coding regions is adequate to resolve phylogenetic relationships among families and some genera of Octocorallia, but insufficient for most species - or population-level studies. Although the COI-COII intergenic spacer exhibits greater variability than the protein-coding regions and may contain useful species-specific markers, its short length limits its phylogenetic utility.

Mitochondrial genome data support the basal position of Acoelomorpha and the polyphyly of the Platyhelminthes

We determined 9.7, 5.2, and 6.8 kb, respectively, of the mitochondrial genomes of the acoel Paratomella rubra, the nemertodermatid Nemertoderma westbladi, and the free-living rhabditophoran platyhelminth Microstomum lineare. The identified gene arrangements are unique among metazoans, including each other, sharing no more than one or two single gene boundaries with a few distantly related taxa. Phylogenetic analysis of the amino acid sequences inferred from the sequenced genes confirms that the acoelomorph flatworms (acoels+nemertodermatids) do not belong to the Platyhelminthes, but are, instead, the most basal extant bilaterian group. Therefore, the Platyhelminthes, as traditionally constituted, is a polyphyletic phylum.

Identification of chaetognaths as protostomes is supported by the analysis of their mitochondrial genome

Determining the phylogenetic position of enigmatic phyla such as Chaetognatha is a longstanding challenge for biologists. Chaetognaths (or arrow worms) are small, bilaterally symmetrical metazoans. In the past decades, their relationships within the metazoans have been strongly debated because of embryological and morphological features shared with the two main branches of Bilateria: the deuterostomes and protostomes. Despite recent attempts based on molecular data, the Chaetognatha affinities have not yet been convincingly defined. To answer this fundamental question, we determined the complete mitochondrial DNA genome of Spadella cephaloptera. We report three unique features: it is the smallest metazoan mitochondrial genome known and lacks both atp8 and atp6 and all tRNA genes. Furthermore phylogenetic reconstructions show that Chaetognatha belongs to protostomes. This implies that some embryological characters observed in chaetognaths, such as a gut with a mouth not arising from blastopore (deuterostomy) and a mesoderm derived from archenteron (enterocoely), could be ancestral features (plesiomorphies) of bilaterians.

The mitochondrial genome of Paraspadella gotoi is highly reduced and reveals that chaetognaths are a sister group to protostomes

We report the complete mtDNA sequence from a member of the phylum Chaetognatha (arrow worms). The Paraspadella gotoi mtDNA is highly unusual, missing 23 of the genes commonly found in animal mtDNAs, including atp6, which has otherwise been found universally to be present. Its 14 genes are unusually arranged into two groups, one on each strand. One group is punctuated by numerous noncoding intergenic nucleotides although the other group is tightly packed, having no noncoding nucleotides, leading to speculation that there are two transcription units with differing modes of expression. The phylogenetic position of the Chaetognatha within the Metazoa has long been uncertain, with conflicting or equivocal results from various morphological analyses and rRNA sequence comparisons. Comparisons here of amino acid sequences from mitochondrially encoded proteins give a single most parsimonious tree that supports a position of Chaetognatha as sister to the protostomes studied here. From this analysis, one can more clearly interpret the patterns of evolution of various developmental features, especially regarding the embryological fate of the blastopore.

Phylogenetic analyses among octocorals (Cnidaria): mitochondrial and nuclear DNA sequences (lsu-rRNA, 16S and ssu-rRNA, 18S) support two convergent clades of branching gorgonians

Gorgonian octocorals lack corroborated hypotheses of phylogeny. This study reconstructs genealogical relationships among some octocoral species based on published DNA sequences from the large ribosomal subunit of the mitochondrial RNA (lsu-rRNA, 16S: 524bp and 21 species) and the small subunit of the nuclear RNA (ssu-rRNA, 18S: 1815bp and 13 spp) using information from insertions-deletions (INDELS) and the predicted secondary structure of the lsu-rRNA (16S). There were seven short (3-10bp) INDELS in the 18S with consistent phylogenetic information. The INDELS in the 16S corresponded to informative signature sequences homologous to the G13 helix found in Escherichia coli. We found two main groups of gorgonian octocorals using a maximum parsimony analysis of the two genes. One group corresponds to deep-water taxa including species from the suborders Calcaxonia and Scleraxonia characterized by an enlargement of the G13 helix. The second group has species from Alcyoniina, Holaxonia and again Scleraxonia characterized by insertions in the 18S. Gorgonian corals, branching colonies with a gorgonin-containing flexible multilayered axis (Holaxonia and Calcaxonia), do not form a monophyletic group. These corroborated results from maternally inherited (16S) and biparentally inherited (18S) genes support a hypothesis of independent evolution of branching in the two octocoral clades.

Unique mitochondrial genome architecture in unicellular relatives of animals

Animal mtDNAs are typically small (approximately 16 kbp), circular-mapping molecules that encode 37 or fewer tightly packed genes. Here we investigate whether similarly compact mitochondrial genomes are also present in the closest unicellular relatives of animals, i.e., choanoflagellate and ichthyosporean protists. We find that the gene content and architecture of the mitochondrial genomes of the choanoflagellate Monosiga brevicollis, the ichthyosporean Amoebidium parasiticum, and Metazoa are radically different from one another. The circular-mapping choanoflagellate mtDNA with its long intergenic regions is four times as large and contains two times as many protein genes as do animal mtDNAs, whereas the ichthyosporean mitochondrial genome totals >200 kbp and consists of several hundred linear chromosomes that share elaborate terminal-specific sequence patterns. The highly peculiar organization of the ichthyosporean mtDNA raises questions about the mechanism of mitochondrial genome replication and chromosome segregation during cell division in this organism. Considering that the closest unicellular relatives of animals possess large, spacious, gene-rich mtDNAs, we posit that the distinct compaction characteristic of metazoan mitochondrial genomes occurred simultaneously with the emergence of a multicellular body plan in the animal lineage.

Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria)

Mitochondrial genes have been used extensively in population genetic and phylogeographical analyses, in part due to a high rate of nucleotide substitution in animal mitochondrial DNA (mtDNA). Nucleotide sequences of anthozoan mitochondrial genes, however, are virtually invariant among conspecifics, even at third codon positions of protein-coding sequences. Hence, mtDNA markers are of limited use for population-level studies in these organisms. Mitochondrial gene sequence divergence among anthozoan species is also low relative to that exhibited in other animals, although higher level relationships can be resolved with these markers. Substitution rates in anthozoan nuclear genes are much higher than in mitochondrial genes, whereas nuclear genes in other metazoans usually evolve more slowly than, or similar to, mitochondrial genes. Although several mechanisms accounting for a slow rate of sequence evolution have been proposed, there is not yet a definitive explanation for this observation. Slow evolution and unique characteristics may be common in primitive metazoans, suggesting that patterns of mtDNA evolution in these organisms differ from that in other animal systems.

The closest unicellular relatives of animals

Molecular phylogenies support a common ancestry between animals (Metazoa) and Fungi, but the evolutionary descent of the Metazoa from single-celled eukaryotes (protists) and the nature and taxonomic affiliation of these ancestral protists remain elusive. We addressed this question by sequencing complete mitochondrial genomes from taxonomically diverse protists to generate a large body of molecular data for phylogenetic analyses. Trees inferred from multiple concatenated mitochondrial protein sequences demonstrate that animals are specifically affiliated with two morphologically dissimilar unicellular protist taxa: Monosiga brevicollis (Choanoflagellata), a flagellate, and Amoebidium parasiticum (Ichthyosporea), a fungus-like organism. Statistical evaluation of competing evolutionary hypotheses confirms beyond a doubt that Choanoflagellata and multicellular animals share a close sister group relationship, originally proposed more than a century ago on morphological grounds. For the first time, our trees convincingly resolve the currently controversial phylogenetic position of the Ichthyosporea, which the trees place basal to Choanoflagellata and Metazoa but after the divergence of Fungi. Considering these results, we propose the new taxonomic group Holozoa, comprising Ichthyosporea, Choanoflagellata, and Metazoa. Our findings provide insight into the nature of the animal ancestor and have broad implications for our understanding of the evolutionary transition from unicellular protists to multicellular animals.

The Mitochondrial Genome of the Sipunculid Phascolopsis gouldii Supports Its Association with Annelida Rather than Mollusca

We determined the sequence of about half (7,470 nt) of the mitochondrial genome of the sipunculid Phascolopsis gouldii, the first representative of this phylum to be so studied. All of the 19 identified genes are transcribed from the same DNA strand. The arrangement of these genes is remarkably similar to that of the oligochaete annelid Lumbricus terrestris. Comparison of both the inferred amino acid sequences and the gene arrangements of a variety of diverse metazoan taxa reveals that the phylum Sipuncula is more closely related to Annelida than to Mollusca. This requires reinterpretation of the homology of several embryological features and of patterns of animal body plan evolution.

Differential localization of nuclear-encoded tRNAs between the cytosol and mitochondrion in Leishmania tarentolae

All mitochondrial tRNAs of the kinetoplastid protozoan Leishmania tarentolae are encoded in the nucleus and are imported from the cytosol into the mitochondrion. We previously reported the partitioning of five tRNAs and found that all were shared between the two compartments to different extents. To increase our knowledge of the tRNAs of this organism, and to attempt to understand the signals involved in their subcellular localization, a method to RT-PCR amplify new tRNAs was developed. Various tRNAs were 3' polyadenylated and reverse transcribed with a sequence-tagged primer. The cDNA was tagged by ligation to an anchor oligonucleotide, and the resulting double-tagged cDNA was amplified by PCR. Four new tRNAs were obtained, bringing to 20 the total number of L. tarentolae tRNAs identified to date. The subcellular localization of 17 tRNAs was quantitatively analyzed by two-dimensional gel electrophoresis and northern hybridization. In general, the previously suggested operational classification of tRNAs into three groups (mainly cytosolic, mainly mitochondrial, and shared between the two compartments) is still valid, but the relative abundance of each tRNA in the cytosol or mitochondrion varied greatly as did the level of expression.

The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa

We sequenced the complete mitochondrial DNA (mtDNA) of the articulate brachiopod Terebratalia transversa. The circular genome is 14,291 bp in size, relatively small compared with other published metazoan mtDNAs. The 37 genes commonly found in animal mtDNA are present; the size decrease is due to the truncation of several tRNA, rRNA, and protein genes, to some nucleotide overlaps, and to a paucity of noncoding nucleotides. Although the gene arrangement differs radically from those reported for other metazoans, some gene junctions are shared with two other articulate brachiopods, Laqueus rubellus and Terebratulina retusa. All genes in the T. transversa mtDNA, unlike those in most metazoan mtDNAs reported, are encoded by the same strand. The A+T content (59.1%) is low for a metazoan mtDNA, and there is a high propensity for homopolymer runs and a strong base-compositional strand bias. The coding strand is quite G+T-rich, a skew that is shared by the confamilial (laqueid) species L. rubellus but is the opposite of that found in T. retusa, a cancellothyridid. These compositional skews are strongly reflected in the codon usage patterns and the amino acid compositions of the mitochondrial proteins, with markedly different usages being observed between T. retusa and the two laqueids. This observation, plus the similarity of the laqueid noncoding regions to the reverse complement of the noncoding region of the cancellothyridid, suggests that an inversion that resulted in a reversal in the direction of first-strand replication has occurred in one of the two lineages. In addition to the presence of one noncoding region in T. transversa that is comparable with those in the other brachiopod mtDNAs, there are two others with the potential to form secondary structures; one or both of these may be involved in the process of transcript cleavage.

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% A+T). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.

Animal phylogeny and the ancestry of bilaterians: inferences from morphology and 18S rDNA gene sequences

Insight into the origin and early evolution of the animal phyla requires an understanding of how animal groups are related to one another. Thus, we set out to explore animal phylogeny by analyzing with maximum parsimony 138 morphological characters from 40 metazoan groups, and 304 18S rDNA sequences, both separately and together. Both types of data agree that arthropods are not closely related to annelids: the former group with nematodes and other molting animals (Ecdysozoa), and the latter group with molluscs and other taxa with spiral cleavage. Furthermore, neither brachiopods nor chaetognaths group with deuterostomes; brachiopods are allied with the molluscs and annelids (Lophotrochozoa), whereas chaetognaths are allied with the ecdysozoans. The major discordance between the two types of data concerns the rooting of the bilaterians, and the bilaterian sister-taxon. Morphology suggests that the root is between deuterostomes and protostomes, with ctenophores the bilaterian sister-group, whereas 18S rDNA suggests that the root is within the Lophotrochozoa with acoel flatworms and gnathostomulids as basal bilaterians, and with cnidarians the bilaterian sister-group. We suggest that this basal position of acoels and gnathostomulids is artifactal because for 1,000 replicate phylogenetic analyses with one random sequence as outgroup, the majority root with an acoel flatworm or gnathostomulid as the basal ingroup lineage. When these problematic taxa are eliminated from the matrix, the combined analysis suggests that the root lies between the deuterostomes and protostomes, and Ctenophora is the bilaterian sister-group. We suggest that because chaetognaths and lophophorates, taxa traditionally allied with deuterostomes, occupy basal positions within their respective protostomian clades, deuterostomy most likely represents a suite of characters plesiomorphic for bilaterians.

Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAS and has a gene arrangement relatable to those of coelomate metazoans

The complete mitochondrial DNA (mtDNA) of the nematode Trichinella spiralis has been amplified in four overlapping fragments and 16,656 bp of its sequence has been determined. This sequence contains the 37 genes typical of metazoan mtDNAs, including a putative atp8, which is absent from all other nematode mtDNAs examined. The genes are transcribed from both mtDNA strands and have an arrangement relatable to those of coelomate metazoans, but not to those of secernentean nematodes. All protein genes appear to initiate with ATN codons, typical for metazoans. Neither TTG nor GTT start codons, inferred for several genes of other nematodes, were found. The 22 T. spiralis tRNA genes fall into three categories: (i) those with the potential to form conventional "cloverleaf" secondary structures, (ii) those with TPsiC arm + variable arm replacement loops, and (iii) those with DHU-arm replacement loops. Mt-tRNA(R) has a 5'-UCG-3' anticodon, as in most other metazoans, instead of the very unusual 5'-ACG-3' present in the secernentean nematodes. The sequence also contains a large repeat region that is polymorphic in size at the population and/or individual level.

The mitochondrial genome of the brachiopod Laqueus rubellus

The complete nucleotide sequence of the 14,017-bp mitochondrial (mt) genome of the articulate brachiopod Laqueus rubellus is presented. Being one of the smallest of known mt genomes, it has an extremely compact gene organization. While the same 13 polypeptides, two rRNAs, and 22 tRNAs are encoded as in most other animal mtDNAs, lengthy noncoding regions are absent, with the longest apparent intergenic sequence being 54 bp in length. Gene-end sequence overlaps are prevalent, and several stop codons are abbreviated. The genes are generally shorter, and three of the protein-coding genes are the shortest among known homologues. All of the tRNA genes indicate size reduction in either or both of the putative TPsiC and DHU arms compared with standard tRNAs. Possession of a TV (TPsiC arm-variable loop) replacement loop is inferred for tRNA(R) and tRNA(L-tag). The DHU arm appears to be unpaired not only in tRNA(S-tct) and tRNA(S-tga), but also in tRNA(C), tRNA(I), and tRNA(T), a novel condition. All the genes are encoded in the same DNA strand, which has a base composition rich in thymine and guanine. The genome has an overall gene arrangement drastically different from that of any other organisms so far reported, but contains several short segments, composed of 2-3 genes, which are found in other mt genomes. Combined cooccurrence of such gene assortments indicates that the Laqueus mt genome is similar to the annelid Lumbricus, the mollusc Katharina, and the octocoral Sarcophyton mt genomes, each with statistical significance. Widely accepted schemes of metazoan phylogeny suggest that the similarity with the octocoral could have arisen through a process of convergent evolution, while it appears likely that the similarities with the annelid and the mollusc reflect phylogenetic relationships.

Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis

The use of some multiple-sequence alignments in phylogenetic analysis, particularly those that are not very well conserved, requires the elimination of poorly aligned positions and divergent regions, since they may not be homologous or may have been saturated by multiple substitutions. A computerized method that eliminates such positions and at the same time tries to minimize the loss of informative sites is presented here. The method is based on the selection of blocks of positions that fulfill a simple set of requirements with respect to the number of contiguous conserved positions, lack of gaps, and high conservation of flanking positions, making the final alignment more suitable for phylogenetic analysis. To illustrate the efficiency of this method, alignments of 10 mitochondrial proteins from several completely sequenced mitochondrial genomes belonging to diverse eukaryotes were used as examples. The percentages of removed positions were higher in the most divergent alignments. After removing divergent segments, the amino acid composition of the different sequences was more uniform, and pairwise distances became much smaller. Phylogenetic trees show that topologies can be different after removing conserved blocks, particularly when there are several poorly resolved nodes. Strong support was found for the grouping of animals and fungi but not for the position of more basal eukaryotes. The use of a computerized method such as the one presented here reduces to a certain extent the necessity of manually editing multiple alignments, makes the automation of phylogenetic analysis of large data sets feasible, and facilitates the reproduction of the final alignment by other researchers.

Mitochondrial genome evolution and the origin of eukaryotes

Recent results from ancestral (minimally derived) protists testify to the tremendous diversity of the mitochondrial genome in various eukaryotic lineages, but also reinforce the view that mitochondria, descendants of an endosymbiotic alpha-Proteobacterium, arose only once in evolution. The serial endosymbiosis theory, currently the most popular hypothesis to explain the origin of mitochondria, postulates the capture of an alpha-proteobacterial endosymbiont by a nucleus-containing eukaryotic host resembling extant amitochondriate protists. New sequence data have challenged this scenario, instead raising the possibility that the origin of the mitochondrion was coincident with, and contributed substantially to, the origin of the nuclear genome of the eukaryotic cell. Defining more precisely the alpha-proteobacterial ancestry of the mitochondrial genome, and the contribution of the endosymbiotic event to the nuclear genome, will be essential for a full understanding of the origin and evolution of the eukaryotic cell as a whole.

Evolution of organellar genomes

Accumulating molecular data, particularly complete organellar genome sequences, continue to advance our understanding of the evolution of mitochondrial and chloroplast DNAs. Although the notion of a single primary origin for each organelle has been reinforced, new models have been proposed that tie the acquisition of mitochondria more closely to the origin of the eukaryotic cell per se than is implied by classic endosymbiont theory. The form and content of the ancestral proto-mitochondrial and proto-chloroplast genomes are becoming clearer but unusual patterns of organellar genome structure and organization continue to be discovered. The 'single-gene circle' arrangement recently reported for dinoflagellate chloroplast genomes is a notable example of a highly derived organellar genome.

Complete DNA sequence of the mitochondrial genome of the ascidian Halocynthia roretzi (Chordata, Urochordata)

The complete nucleotide sequence of the 14,771-bp-long mitochondrial (mt) DNA of a urochordate (Chordata)-the ascidian Halocynthia roretzi-was determined. All the Halocynthia mt-genes were found to be located on a single strand, which is rich in T and G rather than in A and C. Like nematode and Mytilus edulis mtDNAs, that of Halocynthia encodes no ATP synthetase subunit 8 gene. However, it does encode an additional tRNA gene for glycine (anticodon TCT) that enables Halocynthia mitochondria to use AGA and AGG codons for glycine. The mtDNA carries an unusual tRNA(Met) gene with a TAT anticodon instead of the usual tRNA(Met)(CAT) gene. As in other metazoan mtDNAs, there is not any long noncoding region. The gene order of Halocynthia mtDNA is completely different from that of vertebrate mtDNAs except for tRNA(His)-tRNA(Ser)(GCU), suggesting that evolutionary change in the mt-gene structure is much accelerated in the urochordate line compared with that in vertebrates. The amino acid sequences of Halocynthia mt-proteins deduced from their gene sequences are quite different from those in other metazoans, indicating that the substitution rate in Halocynthia mt-protein genes is also accelerated.

Intracellular gene transfer in action: dual transcription and multiple silencings of nuclear and mitochondrial cox2 genes in legumes

The respiratory gene cox2, normally present in the mitochondrion, was previously shown to have been functionally transferred to the nucleus during flowering plant evolution, possibly during the diversification of legumes. To search for novel intermediate stages in the process of intracellular gene transfer and to assess the evolutionary timing and frequency of cox2 transfer, activation, and inactivation, we examined nuclear and mitochondrial (mt) cox2 presence and expression in over 25 legume genera and mt cox2 presence in 392 genera. Transfer and activation of cox2 appear to have occurred during recent legume evolution, more recently than previously inferred. Many intermediate stages of the gene transfer process are represented by cox2 genes in the studied legumes. Nine legumes contain intact copies of both nuclear and mt cox2, although transcripts could not be detected for some of these genes. Both cox2 genes are transcribed in seven legumes that are phylogenetically interspersed with species displaying only nuclear or mt cox2 expression. Inactivation of cox2 in each genome has taken place multiple times and in a variety of ways, including loss of detectable transcripts or transcript editing and partial to complete gene loss. Phylogenetic evidence shows about the same number (3-5) of separate inactivations of nuclear and mt cox2, suggesting that there is no selective advantage for a mt vs. nuclear location of cox2 in plants. The current distribution of cox2 presence and expression between the nucleus and mitochondrion in the studied legumes is probably the result of chance mutations silencing either cox2 gene.

Evolutionary genomics in Metazoa: the mitochondrial DNA as a model system

One of the most important aspects of mitochondrial (mt) genome evolution in Metazoa is constancy of size and gene content of mtDNA, whose plasticity is maintained through a great variety of gene rearrangements probably mediated by tRNA genes. The trend of mtDNA to maintain the same genetic structure within a phylum (e.g., Chordata) is generally accepted, although more recent reports show that a considerable number of transpositions are observed also between closely related organisms. Base composition of mtDNA is extremely variable. Genome GC content is often low and, when it increases, the two complementary bases distribute asymmetrically, creating, particularly in vertebrates, a negative GC-skew. In mammals, we have found coding strand base composition and average degree of gene conservation to be related to the asymmetric replication mechanism of mtDNA. A quantitative measurement of mtDNA evolutionary rate has revealed that each of the various components has a different evolutionary rate. Non-synonymous rates are gene specific and fall in a range comparable to that of nuclear genes, whereas synonymous rates are about 22-fold higher in mt than in nuclear genes. tRNA genes are among the most conserved but, when compared to their nuclear counterparts, they evolve 100 times faster. Finally, we describe some molecular phylogenetic reconstructions which have produced unexpected outcomes, and might change our vision of the classification of living organisms.

Complete sequence of the mitochondrial DNA of the red alga Porphyra purpurea. Cyanobacterial introns and shared ancestry of red and green algae

The mitochondrial DNA (mtDNA) of Porphyra purpurea, a circular-mapping genome of 36,753 bp, has been completely sequenced. A total of 57 densely packed genes has been identified, including the basic set typically found in animals and fungi, as well as seven genes characteristic of protist and plant mtDNAs and specifying ribosomal proteins and subunits of succinate:ubiquinone oxidoreductase. The mitochondrial large subunit rRNA gene contains two group II introns that are extraordinarily similar to those found in the cyanobacterium Calothrix sp, suggesting a recent lateral intron transfer between a bacterial and a mitochondrial genome. Notable features of P. purpurea mtDNA include the presence of two 291-bp inverted repeats that likely mediate homologous recombination, resulting in genome rearrangement, and of numerous sequence polymorphisms in the coding and intergenic regions. Comparative analysis of red algal mitochondrial genomes from five different, evolutionarily distant orders reveals that rhodophyte mtDNAs are unusually uniform in size and gene order. Finally, phylogenetic analyses provide strong evidence that red algae share a common ancestry with green algae and plants.