Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

Genetic characterization of the mitochondrial DNA from Lepeophtheirus salmonis (Crustacea; Copepoda). A new gene organization revealed

The mitochondrial DNA (mtDNA) from the salmon louse, Lepeophtheirus salmonis, is 15445 bp. It includes the genes coding for cytochrome B (Cyt B), ATPase subunit 6 and 8 (A6 and A8), NADH dehydrogenase subunits 1-6 and 4L (ND1, ND2, ND3, ND4, ND4L, ND5 and ND6), cytochrome c oxidase subunits I-III (COI, COII and COIII), two rRNA genes (12S rRNA and 16S rRNA) and 22 tRNAs. Two copies of tRNA-Lys are present in the mtDNA of L. salmonis, while tRNA-Cys was not identified. Both DNA strands contain coding regions in the salmon louse, in contrast to the other copepod characterized Tigriopus japonicus, but only a few genes overlap. In vertebrates, ND4 and ND4L are transcribed as one bicistronic mRNA, and are therefore localized together. The same organization is also found in crustaceans, with the exceptions of T. japonicus, Neocalanus cristatus and L. salmonis that deviate from this pattern. Another exception of the L. salmonis mtDNA is that A6 and A8 do not overlap, but are separated by several genes. The protein-coding genes have a bias towards AT-rich codons. The mitochondrial gene order in L. salmonis differs significantly from the copepods T. japonicus, Eucalanus bungii, N. cristatus and the other 13 crustaceans previously characterized. Furthermore, the mitochondrial rRNA genes are encoded on opposite strands in L. salmonis. This has not been found in any other arthropods, but has been reported in two starfish species. In a phylogenetic analysis, using an alignment of mitochondrial protein sequences, L. salmonis groups together with T. japonicus, being distant relatives to the other crustaceans.     

Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.     

Unusual mitochondrial genome structure of the freshwater goby Odontobutis platycephala: rearrangement of tRNAs and an additional non-coding region

Mitochondrial DNA was isolated from the Korean freshwater gobioid fish Odontobutis platycephala by long-polymerase chain reaction with conserved primers and this mtDNA was sequenced by primer walking using flanking sequences as sequencing primers. The resultant O. platycephala mtDNA sequence was found to be 17 588 bp in size with a mostly conserved structural organization when compared with that of other teleost fish. Rearrangements of tRNAs (tRNA-Ser, tRNA-Leu, tRNA-His) and an additional non-coding region (533 bp) were present between the ND4 and ND5 genes. In the present paper, the basic characteristics of the O. platycephala mitochondrial genome is reported, including its structural organization, base composition of rRNAs, tRNAs and protein-encoding genes, characteristics of mitochondrial tRNAs and the peculiar rearrangement features of some parts of the mtDNA. Phylogenetic analysis performed using the cytochrome b gene sequences of 16 Korean freshwater fishes (15 gobioids) with the Bayesian algorithm showed that O. platycephala forms a clade (1·00 of posterior probability) with other species of Odontobutis . This suggests that the observed rearrangement between the ND4 and ND5 genes in the O. platycephala mitogenome reflects independent events.     

A low rate of replacement substitutions in two major Ovis aries mitochondrial genomes

Two mitochondrial DNA molecules which represent major Ovis aries mtDNA haplogroups were cloned and comparatively sequenced to assess the degree of intraspecific variation. A total of 9623 bp that correspond to 58% of both mitochondrial genomes were determined. The control region, the Cyt b , ND2, ND3, ND4L, COIII and 12 tRNA genes, including the origin of L-strand replication, were completely characterized. Partial sequence information was obtained from the 12S and 16S rRNA and an additional six protein coding and six tRNA genes. The control regions of the two mtDNAs showed a nucleotide divergence of 4·34% while coding regions differed by 0·44%. The number of sheep coding region substitutions was similar to values observed in intraspecific comparisons of mitochondrial DNAs that represent remote points in genealogical trees of mice and humans. However, replacement substitutions were only observed at ∼30% of the rate in mice and ∼20% of the rate in humans. Nucleotide substitutions with a potential for phenotypic effects were found in the 12S and 16S rRNA and in the ND1 and COIII genes.     

Nucleotide Sequence and Gene Organization of the Starfish Asterina Pectinifera Mitochondrial Genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCU)(Ser), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.     

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. [Reviewing Editor: Dr. Yues van de Peer] F. Lamberti: Deceased, 2004     

Phylogenetic relationships of silver crucian carp Carassius auratus gibelio, C. auratus cuvieri, crucian carp C. carassius, and common carp Cyprinus carpio as inferred from mitochondrial DNA variation

PCR-RFLP analysis of the ND3/ND4L/ND4 and 12S/16S rRNA regions and nucleotide sequence variation of the cytochrome b gene were used to study the mtDNA divergence in species of the family Cyprinidae, to examine the phylogenetic relationships of the species, and to identify their taxonomic status. The results indicated that an ancestral form diverged into silver crucian carp and crucian carp after its separation from the common carp lineage. The divergence of continental Carassius auratus gibelio and Japanese C. auratus cuvieri occurred more recently. Two well distinguishable mtDNA phylogroups, suggesting divergent evolution, were observed in continental C. auratus gibelio populations. The divergence was possibly related to the formation of two silver crucian carp groups with different types of reproduction, triploid gynogenetic and diploid gonochoric. At the same time, the results supported the high probability of current genetic exchange between the forms. In view of these findings and high morphological similarity of the two forms, they were not considered to be separate species.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

The complete nucleotide sequence of the mitochondrial genome of Tetraodon nigroviridis.

The fresh water pufferfish Tetraodon nigroviridis is a model organism for studying evolution of genome and gene functions, but its mitochondrial genome (mtDNA) sequence is still not available. We determined the complete nucleotide sequence of its mtDNA using shotgun sequencing. The T. nigroviridis mtDNA was 16,462 bp, and contained 13 protein coding genes, 22 tRNAs, 2 rRNAs and a major non-coding region. The gene order was identical to the common type of vertebrate mtDNA, whereas the G + C content in the sense strand was 46.9%, much higher than most other fish species. One hundred and three SNPs were detected in the control region of the mtDNA of 35 individuals, a majority of SNPs were detected in the 5' end of the control region. A phylogenetic study including 21 fish species was performed on concatenated amino acid sequences of 12 protein coding genes, and revealed that the T. nigroviridis was clustered with Fugu rubripes into a group. The complete mtDNA sequence and SNPs in its control region will be useful in studying fish evolution, in differentiating different Tetraodon species and in analyzing genetic diversity within T. nigroviridis.     

A compositional segmentation of the human mitochondrial genome is related to heterogeneities in the guanine mutation rate

We applied a hidden Markov model segmentation method to the human mitochondrial genome to identify patterns in the sequence, to compare these patterns to the gene structure of mtDNA and to see whether these patterns reveal additional characteristics important for our understanding of genome evolution, structure and function. Our analysis identified three segmentation categories based upon the sequence transition probabilities. Category 2 segments corresponded to the tRNA and rRNA genes, with a greater strand-symmetry in these segments. Category 1 and 3 segments covered the protein- coding genes and almost all of the non-coding D-loop. Compared to category 1, the mtDNA segments assigned to category 3 had much lower guanine abundance. A comparison to two independent databases of mitochondrial mutations and polymorphisms showed that the high substitution rate of guanine in human mtDNA is largest in the category 3 segments. Analysis of synonymous mutations showed the same pattern. This suggests that this heterogeneity in the mutation rate is partly independent of respiratory chain function and is a direct property of the genome sequence itself. This has important implications for our understanding of mtDNA evolution and its use as a ‘molecular clock’ to determine the rate of population and species divergence.     

Nucleotide sequence of nine protein-coding genes and 22 tRNAs in the mitochondrial DNA of the sea starPisaster ochraceus

Summary We have cloned and sequenced over 9 kb of the mitochondrial genome from the sea starPisaster ochraceus. Within a continuous 8.0-kb fragment are located the genes for NADH dehydrogenase subunits 1, 2, 3, and 4L (ND1, ND2, ND3, and ND4L), cytochrome oxidase subunits I, II, and III (COI, COII, and COIII), and adenosine triphosphatase subunits 6 and 8 (ATPase 6 and ATPase 8). This large fragment also contains a cluster of 13 tRNA genes between ND1 and COI as well as the genes for isoleucine tRNA between ND1 and ND2, arginine tRNA between COI and ND4L, lysine tRNA between COII and ATPase 8, and the serine (UCN) tRNA between COIII and ND3. The genes for the other five tRNAs lie outside this fragment. The gene for phenylalanine tRNA is located between cytochrome b and the 12S ribosomal genes. The genes for tRNA^glu and tRNA^thr are 3 to the 12S ribosomal gene. The tRNAs for histidine and serine (AGN) are adjacent to each other and lie between ND4 and ND5. These data confirm the novel gene order in mitochondrial DNA (mtDNA) of sea stars and delineate additional distinctions between the sea star and other mtDNA molecules.     

Phylogenetic relationships of silver crusian carp Carassius auratus gibelio, C. auratus cuvieri, crucian carp Carassius carassius, and common carp Cyprinus carpio as inferred from mitochondrial DNA variation

PCR-RFLP analysis of the ND3/ND4L/ND4 and 12S/16S rRNA regions and nucleotide sequence variation of the cytochrome b gene were used to study the mtDNA divergence in species of the family Cyprinidae, to examine the phylogenetic relationships of the species, and to identify their taxonomic status. The results indicated that an ancestral form diverged into silver crucian carp and crucian carp after its separation from the common carp lineage. The divergence of continental Carassius auratus gibelio and Japanese C. auratus cuvieri occurred more recently. Two well distinguishable mtDNA phylogroups, suggesting divergent evolution, were observed in continental C. auratus gibelio populations. The divergence was possibly related to the formation of two silver crucian carp groups with different types of reproduction, triploid gynogenetic and diploid gonochoric. At the same time, the results supported the high probability of current genetic exchange between the forms. In view of these findings and high morphological similarity of the two forms, they were not considered to be separate species.     

The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: fivefold tandem repetition of a coding region.

We sequenced across all of the gene boundaries in the mitochondrial genome of the cattle tick, Boophilus microplus, to determine the arrangement of its genes. The mtDNA of B. microplus has a coding region, composed of tRNA(Glu) and 60 bp of the 3' end of ND1, that is repeated five times. Boophilus microplus is the first coelomate animal known to have more than two copies of a coding sequence. The mitochondrial genome of B. microplus has other unusual features, including (1) reduced T arms in tRNAs, (2) an AT bias in codon use, (3) two control regions that have evolved in concert, (4) three gene rearrangements, and (5) a stem-loop between tRNA(Gln) and tRNA(Phe). The short T arms and small control regions (CRs) of B. microplus and other ticks suggest strong selection for small genomes. Imprecise termination of replication beyond its origin, which can account for the evolution of tandem repeats of coding regions in other mitochondrial genomes, cannot explain the evolution of the fivefold repeated sequence in the mitochondrial genome of B. microplus. Instead, slipped-strand mispairing or recombination are the most plausible explanations for the evolution of these tandem repeats.