Sequence variation and structural conservation in the D-loop region and flanking genes of mitochondrial DNA from Japanese pond frogs

The nucleotide sequences of the D-loop region and its flanking genes of the mitochondrial DNA (mtDNA) from Japanese pond frogs were determined by the methods of PCR, cloning, and sequencing. The frogs belonged to two species, one subspecies, and one local race. The gene arrangements adjacent to the D-loop region were analyzed. The frogs shared a unique mitochondrial gene order that was found in Rana catesbeiana; i.e., cyt b--D-loop region--tRNA(Leu(CUN))--tRNA(Thr)--tRNA(Pro)--tRNA(Phe)--12S rRNA. The arrangements of the three tRNA genes of these frogs were different from those of X. laevis, a species which has the same overall structure as in mammals. Highly repetitive sequences with repeat units (16-bp or 17-bp sequence specific for each taxon) were found in the D-loop region. The length of repetitive sequences varied from 0.6 kbp to 1.2 kbp, and caused the extensive size variation in mtDNA. Several short sequence elements such as putative TAS, OH, CSB-1, and CSB-2 were found in the D-loop region of these frogs. The sequences of these short regulatory elements were conserved in R. catesbeiana, X. laevis, and also in human. The comparison of sequence divergences of the D-loop region and its adjacent genes among various taxa revealed that the rates of nucleotide substitutions depend on genes. The nucleotide sequences of the 3'-side segment of the D-loop region were the most variable among taxa, whereas those of the tRNA and 12S rRNA genes were the most conservative.     

The mitochondrial NADH dehydrogenase subunit 6 (ND6) gene in Murres: relevance to phylogenetic and population studies among birds.

The nucleotide sequences of the mitochondrial ND6 and tRNA(Glu) genes and part of the displacement loop region in two closely related seabird species are presented. A chicken type gene organization in which the tRNA(Glu), ND6, and displacement loop are localized next to each other was found in these species and suggests that this is a conserved feature of avian mitochondrial DNA. The nucleotide and amino acid divergences of ND6 at different taxonomic levels are assessed, and its relevance to phylogenetic studies in birds is discussed.     

Mitochondrial DNA in the sea urchin Arbacia lixula: evolutionary inferences from nucleotide sequence analysis.

From the stirodont Arbacia lixula we determined the sequence of 5,127 nucleotides of mitochondrial DNA (mtDNA) encompassing 18 tRNAs, two complete coding genes, parts of three other coding genes, and part of the 12S ribosomal RNA (rRNA). The sequence confirms that the organization of mtDNA is conserved within echinoids. Furthermore, it underlines the following peculiar features of sea urchin mtDNA: the clustering of tRNAs, the short noncoding regulatory sequence, and the separation by the ND1 and ND2 genes of the two rRNA genes. Comparison with the orthologous sequences from the camarodont species Paracentrotus lividus and Strongylocentrotus purpuratus revealed that (1) echinoids have an extra piece on the amino terminus of the ND5 gene that is probably the remnant of an old leucine tRNA gene; (2) third-position codon nucleotide usage has diverged between A. lixula and the camarodont species to a significant extent, implying different directional mutational pressures; and (3) the stirodont-camarodont divergence occurred twice as long ago as did the P. lividus-S. purpuratus divergence.     

Complete nucleotide sequence and gene rearrangement of the mitochondrial genome of the bell-ring frog, Buergeria buergeri (family Rhacophoridae)

In this study we determined the complete nucleotide sequence (19,959 bp) of the mitochondrial DNA of the rhacophorid frog Buergeria buergeri. The gene content, nucleotide composition, and codon usage of B. buergeri conformed to those of typical vertebrate patterns. However, due to an accumulation of lengthy repetitive sequences in the D-loop region, this species possesses the largest mitochondrial genome among all the vertebrates examined so far. Comparison of the gene organizations among amphibian species (Rana, Xenopus, salamanders and caecilians) revealed that the positioning of four tRNA genes and the ND5 gene in the mtDNA of B. buergeri diverged from the common vertebrate gene arrangement shared by Xenopus, salamanders and caecilians. The unique positions of the tRNA genes in B. buergeri are shared by ranid frogs, indicating that the rearrangements of the tRNA genes occurred in a common ancestral lineage of ranids and rhacophorids. On the other hand, the novel position of the ND5 gene seems to have arisen in a lineage leading to rhacophorids (and other closely related taxa) after ranid divergence. Phylogenetic analysis based on nucleotide sequence data of all mitochondrial genes also supported the gene rearrangement pathway.     

Mitochondrial genome of the cockscomb pearl mussel Cristaria plicata (Bivalvia, Unionoida, Unionidae)

The complete mitochondrial genome sequence of the cockscomb pearl mussel Cristaria plicata, which is an endangered species in South Korea, was sequenced. The circle genome (15,708 bp in size) consists of 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes. There were 26 noncoding regions (NCs) found throughout the mitogenome of C. plicata, ranging in size from 2 to 327 bp, and the two largest NC regions, NC286 and NC326, were found between ND5 and tRNA(Gln) (286 bp) and between tRNA(Glu) and ND2 (326 bp), respectively. The 13 mitochondrial protein-coding genes of a female individual of C. plicata collected from Korea (15,708 bp) were compared to those of the Chinese individual (15,712 bp) published before. The result showed that ND3 is the most conserved with 100% nucleotide similarity, and each of the other protein-coding genes has ca. 99%, respectively. The two largest NCs among 26 NCs have totally 98% nucleotide similarity between Korean and Chinese ones.     

Nucleotide sequence of a Xenopus laevis mitochondrial DNA fragment containing the D-loop, flanking tRNA genes and the apocytochrome b gene

Extensive corrections of the nucleotide sequence of the Xenopus laevis mitochondrial (mt) displacement (D) loop and surrounding genes [Wong et al., Nucl. Acids Res. 11 (1983) 4977-4995] are reported, including addition of two stretches of nucleotides and 60 scattered modifications. The additional sequences presented here correspond to the apocytochrome b gene, the tRNAGlu gene and part of URF6. This allows us to propose a conformational model for the X. laevis apocytochrome b protein and also permits comparisons with mammalian mtDNA. The D-loop sequence is poorly conserved except for sequences involved in the regulation of the mt genome (conserved sequence blocks and the DNA polymerase stop sequences). On the other hand, all genes show marked conservation both of their nucleotide sequence and their respective location on the mt genome. Organization of the genetic information described for mammalian mtDNA also holds for the X. laevis mtDNA. This result strongly suggests that all animal vertebrate mtDNAs have followed the same evolutionary pathway.     

The complete mitochondrial genome of the marbled rockfish Sebastiscus marmoratus (Scorpaeniformes, Scorpaenidae): genome characterization and phylogenetic considerations

The complete mitochondrial genome sequence of the marbled rockfish Sebastiscus marmoratus (Scorpaeniformes, Scorpaenidae) was determined and phylogenetic analysis was conducted to elucidate the evolutionary relationship of the marbled rockfish with other Sebastinae species. This mitochondrial genome, consisting of 17301 bp, is highly similar to that of most other vertebrates, containing the same gene order and an identical number of genes or regions, including 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and one putative control region. Most of the genes are encoded on the H-strand, while the ND6 and seven tRNA genes (for Gln, Ala, Asn, Tyr, Ser (UCA), Glu, and Pro) are encoded on the L-strand. The reading frame of two pairs of genes overlapped on the same strand (the ATPase 8 and 6 genes overlapped by ten nucleotides; ND4L and ND4 genes overlapped by seven nucleotides). The possibly nonfunctional light-strand replication origin folded into a typical stem-loop secondary structure and a conserved motif (5'-GCCGG-3') was found at the base of the stem within the tRNA(Cys) gene. An extent termination-associated sequence (ETAS) and conserved sequence blocks (CSB) were identified in the control region, except for CSB-1; unusual long tandem repeats were found at the 3' end of the control region. Phylogenetic analyses supported the view that Sebastinae comprises four genera (Sebates, Hozukius, Helicolenus, and Sebasticus).     

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.     

Evolution of the WANCY region in amniote mitochondrial DNA

In most vertebrate mitochondrial genomes, the site for initiation of light-strand replication, OL, is found within a cluster of five transfer RNA (tRNA) genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), and tRNA(Tyr)). This region and part of the adjacent cytochrome c oxydase subunit I (COI) gene were sequenced for two crocodilian, two turtle, and one snake species and for Sphenodon punctatus; part of the adjacent nicotinamide adenine dinucleotide dehydrogenase subunit 2 (ND2) gene was also sequenced for the crocodilian and turtle species. All had the typical vertebrate gene order. The turtles and the snake have a lengthy noncoding sequence between the tRNA(Asn) and tRNA(Cys) genes that we assumed to be homologous to the mammalian OL. The crocodilians and Sphenodon lack such a sequence, a condition they share with birds. Most proposed phylogenies for the amniotes require that OL at this position was lost at least twice during their diversification or was evolved independently more than once. Within the five tRNA genes, frequencies of substitutions are much higher in loops than in stems. Many loops vary dramatically in size among the species; in the most extreme case, the D-arm of the Sphenodon tRNA(Cys) is a "D-arm replacement" loop of seven nucleotides. Frequency of transitions in stems is relatively uniform across tRNAs, but frequency of transversions varies greatly. Mismatches in stems are infrequent, and their relative frequency in a specific tRNA is unrelated to the frequency of substitution in the corresponding gene. Several features of mammalian mitochondrial tRNAs are conserved in WANCY tRNAs throughout amniotes. The inferred initiation codon for COI is GTG in crocodilians, turtles, and the snake, a condition they share with fishes, certain amphibians, and birds. TTG appears to be the initiation codon for COI in Sphenodon; if correct, this would be a novel initiation codon for vertebrate mitochondrial DNA. Phylogenetic analyses of the inferred amino acid sequences of ND2 and COI support the sister-group relationship of birds and crocodilians and suggest that mammals are an early derived lineage within the amniotes.      

The complete mitochondrial genome of silver croaker Argyrosomus argentatus (Perciforems; Sciaenidae): genome characterization and phylogenetic consideration

The complete mitochondrial genome sequence of the silver croaker, Argyrosomus argentatus, was obtained by using LA-PCR and sequencing. The mitogenome is 16485 bp in length, consists of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and a non-coding control region like those found in other vertebrates, with the gene order similar to that of typical teleosts. Most of the genes of A. argentatus were encoded on the H-strand, while the ND6 and eight tRNA (Gln, Ala, Asn, Cys, Tyr, Ser (UCN), Glu and Pro)) genes were encoded on the L-strand. The reading frames of two pairs of genes overlapped: ATPase8 and 6 and ND4L and ND4 by ten and seven nucleotides, respectively. The origin of L-strand replication in A. argentatus was in a cluster of five tRNA genes (WANCY) and was 46 nucleotides in length. The conserved motif (5'-GCGGG-3') was found at the base of the stem within the tRNA(Cys) gene. Within the control region, we identified all of the conserved motifs except for CSB-F.     

Mitochondrial sequence evolution in spiders: intraspecific variation in tRNAs lacking the TPsiC Arm

Analyses of mitochondrial DNA sequences from three species of Habronattus jumping spiders (Chelicerata: Arachnida: Araneae) reveal unusual inferred tRNA secondary structures and gene arrangements, providing new information on tRNA evolution within chelicerate arthropods. Sequences from the protein-coding genes NADH dehydrogenase subunit 1 (ND1), cytochrome oxidase subunit I (COI), and subunit II (COII) were obtained, along with tRNA, tRNA, and large-subunit ribosomal RNA (16S) sequences; these revealed several peculiar features. First, inferred secondary structures of tRNA and, likely, tRNA, lack the TPsiC arm and the variable arm and therefore do not form standard cloverleaf structures. In place of these arms is a 5-6-nt T arm-variable loop (TV) replacement loop such as that originally described from nematode mitochondrial tRNAs. Intraspecific variation occurs in the acceptor stem sequences in both tRNAs. Second, while the proposed secondary structure of the 3' end of 16S is similar to that reported for insects, the sequence at the 5' end is extremely divergent, and the entire gene is truncated about 300 nt with respect to Drosophila yakuba. Third, initiation codons appear to consist of ATY (ATT and ATC) and TTG for ND1 and COII, respectively. Finally, Habronattus shares the same ND1-tRNA-16S gene arrangement as insects and crustaceans, thus illustrating variation in a tRNA gene arrangement previously proposed as a character distinguishing chelicerates from insects and crustaceans.     

Complete mitochondrial genome of the rabbitfish Siganus fuscescens (Perciformes, Siganidae).

We determined the complete nucleotide sequence of the mitochondrial genome for the rabbitfish Siganus fuscescens (Perciformes, Siganidae). This mitochondrial genome, consisting of 16,491 base pairs (bp), included 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a noncoding control region similar those found in other vertebrates; the gene order was identical to that of typical vertebrates. Most of the genes of S. fuscescens were encoded on the H-strand, while the ND6 and eight tRNA (Gln, Ala, Asn, Cys, Tyr, Ser [UCN], Glu, and Pro) genes were encoded on the L-strand. The reading frames of ATPase 8 and 6 and those of ND4L and ND4 overlapped by ten and seven nucleotides, respectively. All mitochondrial protein-coding genes began with an ATG start codon, except for CO1, which started with GTG. Open reading frames of S. fuscescens ended with TAA (ND1, CO1, ATPase 8, ND4L, ND5 and ND6), and the remainder had incomplete stop codons, either TA (ATPase 6 and CO3) or T (ND2, CO2, ND3, ND4, and Cytb). The origin of L-strand replication in S. fuscescens was located in a cluster of five tRNA genes (WANCY) and was 34 nucleotides in length. A major noncoding region between the tRNA-Pro and tRNA-Phe genes (828 bp) was considered to be the control region (D-loop). Within this sequence, we identified a conserved sequence block characteristic of this region. The rabbitfish was grouped with Siganus canaliculatus in most parsimony analyses, which showed 100% bootstrap support for their divergence. These findings are useful for inferring phylogenetic relationships and identification within the suborder Acanthuroidei.     

Molecular phylogenetic relationships of pond frogs distributed in the Palearctic region inferred from DNA sequences of mitochondrial 12S ribosomal RNA and cytochrome b genes

The evolutionary relationships of pond frogs distributed in the Far East and Europe were investigated by analyses of nucleotide sequences of mitochondrial 12S ribosomal RNA (12S rRNA) and cytochrome b (cyt b) genes. The nucleotide sequences of a 412-bp segment of the 12S rRNA gene and a 534-bp segment of the cyt b gene were determined by the PCR-direct sequencing method using 19 frogs belonging to six species and one subspecies distributed in the Palearctic region. Phylogenetic trees were constructed by the neighbor-joining and maximum-likelihood methods using Rana catesbeiana or Xenopus laevis as an outgroup. The 412-bp segment of the 12S rRNA gene contained 65 variable sites including gap sites, and the 534-bp segment of the cyt b gene contained 160 variable sites. The nucleotide sequence divergences of the 12S rRNA gene were 0.25-4.83% within the Far Eastern frogs, 0.25-6.22% within the European frogs, and 8.74-11.24% between the Far Eastern and the European frogs, whereas those of the cyt b gene were 3.64-14.73% within the Far Eastern frogs, 0.38-14.42% within the European frogs, and 16.53-23.58% between the Far Eastern and the European frogs. Although most nucleotide substitutions were at the third codon position of the cyt b gene and were silent mutations, 4 amino acid replacements occurred within the Far Eastern frogs, 4 within the European frogs, and 11 between the Far Eastern and the European frogs. The phylogenetic trees constructed from the nucleotide sequence divergences showed slightly different topologies for the 12S rRNA and cyt b genes. R. esculenta from Ukraine was closely related to R. lessonae from Luxembourg in both the 12S rRNA and the cyt b gene sequences.     

Unusual genetic codes and a novel gene structure for tRNA(AGYSer) in starfish mitochondrial DNA

The nucleotide sequence of a 3849-bp fragment of starfish mitochondrial genome was determined. The genes for NADH dehydrogenase subunits 3, 4, 5, and COIII, and three kinds of (tRNA(UCNSer), tRNA(His), and tRNA(AGYSer) were identified by comparing with the genes of other animal mitochondria so far elucidated. The gene arrangement of starfish mitochondrial genome was different from those of vertebrate and insect mitochondrial genomes. Comparison of the protein-encoding nucleotide sequences of starfish mitochondria with those of other animal mitochondria suggested a unique genetic code in starfish mitochondrial genome; both AGA and AGG (arginine in the universal code) code for serine, AUA (isoleucine in the universal code but methionine in most mitochondrial systems) for isoleucine, and AAA (lysine) for asparagine. It was also inferred that these AGA and AGG codons are decoded by serine tRNA(AGYSer) originally corresponding to AGC and AGU codons. This situation is similar to the case of Drosophila mitochondrial genome. Variations in the use of AGA and AGG codons were discussed on the basis of the evolution of animals and decoding capacity of various tRNA(AGYSer) species possessing different sizes of the dihydrouridine (D) arm.     

三斑海马线粒体基因组核苷酸全序列结构与分析

通过PCR扩增并测序获得了三斑海马(Hippocampus trimaculatus)线粒体DNA(mt DNA)全序列。三斑海马线粒体基因组全序列长度为16 534 bp(Gen Bank登录号为KJ956892),编码37个基因,包括13个蛋白编码基因、22个t RNA基因和2个r RNA基因。非编码区域包括1个控制区(D-loop)及一个轻链复制起始区域。大部分基因由H-链编码,包括14个t RNA基因、2个r RNA基因、12个蛋白编码基因;只有ND6和8个t RNA基因是在L-链编码。预测的22个t RNA基因的二级结构均为典型的三叶草状。基因间隔一般1~14 bp不等。此外,还存在7处碱基重叠,其中,4处是鱼类和脊椎动物典型的基因重叠位点。总的碱基含量分别为,A 32.7%,C 23.4%,G 14.6%,T 29.3%,A+T含量为62.0%。其线粒体基因组序列的结构与脊椎动物的典型结构近似。邻接法和贝叶斯法构建的三斑海马系统进化树的拓扑结构相似,这与现有的三斑海马的系统演化地位一致。本研究为海马的进化研究以及保护工作提供了基础数据。     

Mitochondrial genome of the eurasian otterLutra lutra (Mammalia, Carnivora, Mustelidae)

We determined the complete mitochondrial genome of the Eurasian otterLutra lutra, which is an endangered species in Korea. The circle genome (16,536 bp in size) consists of 13 protein-coding, 22 tRNA, and 2 rRNA genes, and a control region, as found in other metazoan animals. Out of the 37 genes, 28 are encoded on the H-strand, and the nine (ND6 and 8 tRNA genes) on the L-strand. Three overlaps among the 13 protein-coding genes were found: ATP8-ATP6, ND4L-ND4, and ND5-ND6. A control region (1090 bp) including the origin of H-strand replication (OH), TAS (a conserved motif TACAT-16bp-ATGTA) and CSB (CSB-1, CSB-2. and CSB-3) was observed between tRNA-Pro and tRNA-Phe genes, and O_L, with 36 highly conserved nucleotides between tRNA-Asn (N) and tRNA-Cys (C) within a cluster of five tRNA genes (WANCY), as typically found in vertebrates. The other important characteristics of theL. lutra mitochondrial genome were described in detail. In addition, a maximum likelihood and Bayesian trees of 9 mustelid species and 1 outgroup were reconstructed based on the nucleotide sequences of 11 protein-coding genes excluding ATP8 and ND6. It showed that Lutrinae formed a monophyletic group with Mustelinae that is not monophyletic. Within the subfamily Lutrinae,L. lutra andEnhydra lutris were grouped together and thenLontra canadentis placed as a sister of the clade. The present result is the first complete mitochondrial genome sequence reported from the genusLutra, and is applicable to molecular phylogenetic, phylogeographic, conservation biological studies for mustelid members. In particular, exploration of sequence variations of the control region may be helpful for analyzing inter-and intra-species variations in the genusLutra.     

A novel mitochondrial gene order in the crinoid echinoderm Florometra serratissima

The complete nucleotide sequence of the mitochondrial genome of the crinoid Florometra serratissima has been determined. It is a circular DNA molecule, 16,005 bp in length, containing the genes for 13 proteins, small and large ribosomal RNAs, and 22 transfer RNAs (tRNAs). Three regions of unassigned sequence (UAS) greater than 73 bp have been located. The largest, UAS I, is 432 bp long and exhibits sequence similarity to the putative mitochondrial control regions seen in other animals. UAS II (77 bp) and UAS III (73 bp) are located between the 5' ends of coding sequences and may play roles as bidirectional promoters. Analyses of nucleotide composition revealed that the major peptide-encoding strand is high in T and low in C. This bias is reflected in a specific pattern of codon usage. Molecular phylogenetic analyses based on cytochrome c oxidase (COI, COII, and COIII) amino acid and nucleotide sequences did not resolve all the relationships between echinoderm classes. The overall animal mitochondrial gene content has been maintained in the crinoid, but there is extensive rearrangement with respect to both the echinoid and the asteroid mtDNA gene maps. Florometra serratissima has a novel genome organization in a segment containing most of the tRNA genes, large and small rRNA genes, and the NADH dehydrogenase subunit 1 and 2 genes. Potential pathways and mechanisms for gene rearrangements between mitochondrial gene maps of echinoderm classes and vertebrates are discussed as indicators of early deuterostome phylogeny.